Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,281,408 G→T 10.0% I20I (ATC→ATA)  PP_3751 ← cyanate transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,281,4080GT10.0% 143.1 / 5.9 50I20I (ATC→ATA) PP_3751cyanate transporter
Reads supporting (aligned to +/- strand):  ref base G (18/27);  new base T (3/2);  total (21/29)
Fisher's exact test for biased strand distribution p-value = 6.38e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.72e-01

GCTCGCTGACCCAGCGTCGCAGCCACGGTTGCAACAGGGGCACCAGGCCCATGCACAGCACCGGTAGCGCCGTCAGCAAGGCCGCCTGCTGGTAACCCAGGCCGGTGCTGGCGCGCATGGGTTCGAGTAACGGGCCGATACTGGTCAGGATCGGCCGCAGGTTCAACGCCAGCAGGATCACCAGCAACAGGTTCAGTACAGCGCTCATTGTGAGTCGGTTCTCTTGTAACAGGAACTGGCCATTATCCCGGCTCGGCTGGCCTGGCAGAAATGAAGAGATTGGATGCAGGTCAGTGGC  >  NC_002947/4281260‑4281557
                                                                                                                                                    |                                                                                                                                                     
gCTCGCTGACCCAGCGTCGCAGCCACGGTTGCAACAGGGGCACCAGGCCCATGCACAGCACCGGTAGCGCCGTCAGCAAGGCCGCCTGCTGGTAACCCAGGCCGGTGCTGGCGCGCATGGGTTCGAGTAACGGGCCGATACTGGTCAGGa                                                                                                                                                      <  1:608335/150‑1 (MQ=255)
     cTGACCCAGCGTCGCAGCCACGGTTGCAACAGGGGCACCAGGCCCATGCACAGCACCGGTAGCGCCGTCAGCAAGGCCGCCTGCTGGTAACCCAGGCCGGTGCTGGCGCGCATGGGTTCGAGTAACGGGCCGATCCTGACCAGTATCGGc                                                                                                                                                 <  2:445346/150‑1 (MQ=255)
       gACCCAGCGTCGCAGCCACGGTTGCAACAGGGGCACCAGGCCCATGCACAGCACCGGTAGCGCCGTCAGCAAGGCCGCCTGCTGGTAACCCAGGCCGGTGCTGGCGCGCATGGGTTCGAGTAACGGGCCGATACTGGTCAGGATCGGCCg                                                                                                                                               <  2:693616/150‑1 (MQ=255)
          ccAGCGTCGCAGCCACGGTTGCAACAGGGGCACCAGGCCCATGCACAGCACCGGTAGCGCCGTCAGCAAGGCCGCCTGCTGGTAACCCAGGCCGGTGCTGGCGCGCATGGGTTCGAGTAACGGGCCGATACTGGTCAGGATCGGCCGCAg                                                                                                                                            >  1:242065/1‑150 (MQ=255)
           cAGCGTCGCAGCCACGGTTGCAACAGGGGCACCAGGCCCATGCACAGCACCGGTAGCGCCGTCAGCAAGGCCGCCTGCTGGTAACCCAGGCCGGTGCTGGCGCGCATGGGTTCGAGTAACGGGCCGATACTGGTCAGGATCGGCCGCAgg                                                                                                                                           >  1:747116/1‑150 (MQ=255)
                tCGCAGCCACGGTTGCAACAGGGGCACCAGGCCCATGCACAGCACCGGTAGCGCCGTCAGCAAGGCCGCCTGCTGGTAACCCAGGCCGGTGCTGGCGCGCATGGGTTCGAGTAACGGGCCGATACTGGTCAGGATCGGCCGCAGGTTCaa                                                                                                                                      <  1:179780/150‑1 (MQ=255)
                       cACGGTTGCAACAGGGGCACCAGGCCCATGCACAGCACCGGTAGCGCCGTCAGCAAGGCCGCCTGCTGGTAACCCAGGCCGGTGCTGGCGCGCATGGGTTCGAGTAACGGGCCGATACTGGTCAGGATCGGCCGCAGGTTCAACGCcagc                                                                                                                               <  1:136126/150‑1 (MQ=255)
                        aCGGTTGCAACAGGGGCACCAGGCCCATGCACAGCACCGGTAGCGCCGTCAGCAAGGCCGCCTGCTGGTAACCCAGGCCGGTGCTGGCGCGCATGGGTTCGAGTAACGGGCCGATACTGGTCAGGATCGGCCGCAGGTTCAACGCcagca                                                                                                                              >  2:1117507/1‑150 (MQ=255)
                        aCGGTTGCAACAGGGGCACCAGGCCCATGCACAGCACCGGTAGCGCCGTCAGCAAGGCCGCCTGCTGGTAACCCAGGCCGGTGCTGGCGCGCATGGGTTCGAGTAACGGGCCGATACTGGTCAGGATCGGCCGCAGGTTCAACGCcagca                                                                                                                              <  1:916406/150‑1 (MQ=255)
                               cAACAGGGGCACCAGGCCCATGCACAGCACCGGTAGCGCCGTCAGCAAGGCCGCCTGCTGGTAACCCAGGCCGGTGCTGGCGCGCATGGGTTCGAGTAACGGGCCGATACTGGTCAGGATCGGCCGCAGGTTCAACGCCAGCAGGATCAc                                                                                                                       <  1:903591/150‑1 (MQ=255)
                                 aCAGGGGCACCAGGCCCATGCACAGCACCGGTAGCGCCGTCAGCAAGGCCGCCTGCTGGTAACCCAGGCCGGTGCTGGCGCGCATGGGTTCGAGTAACGGGCCGATACTGGTCAGGATCGGCCGCAGGTTCAACGCCAGCAGGATCACCa                                                                                                                     >  2:832436/1‑150 (MQ=255)
                                   aGGGGCACCAGGCCCATGCACAGCACCGGTAGCGCCGTCAGCAAGGCCGCCTGCTGGTAACCCAGGCCGGTGCTGGCGCGCATGGGTTCGAGTAACGGGCCGATCCTGACCAGTATCGGccc                                                                                                                                               <  1:139020/122‑2 (MQ=255)
                                   aGGGGCACCAGGCCCATGCACAGCACCGGTAGCGCCGTCAGCAAGGCCGCCTGCTGGTAACCCAGGCCGGTGCTGGCGCGCATGGGTTCGAGTAACGGGCCGATCCTGACCAGTATCGGccc                                                                                                                                               >  2:139020/1‑121 (MQ=255)
                                        cACCAGGCCCATGCACAGCACCGGTAGCGCCGTCAGCAAGGCCGCCTGCTGGTAACCCAGGCCGGTGCTGGCGCGCATGGGTTCGAGTAACGGGCCGATACTGGTCAGGATCGGCCGCAGGTTCAACGCCAGCAGGATCACCAGCAACAg                                                                                                              <  1:302401/150‑1 (MQ=255)
                                         aCCAGGCCCATGCACAGCACCGGTAGCGCCGTCAGCAAGGCCGCCTGCTGGTAACCCAGGCCGGTGCTGGCGCGCATGGGTTCGAGTAACGGGCCGATACTGGTCAGGATCGGCCGCAGGTTCAACGCCAGCAGGATCACCAGCAACAgg                                                                                                             >  2:188993/1‑150 (MQ=255)
                                              gCCCATGCACAGCACCGGTAGCGCCGTCAGCAAGGCCGCCTGCTGGTAACCCAGGCCGGTGCTGGCGCGCATGGGTTCGAGTAACGGGCCGATACTGGTCAGGATCGGCCGCAGGTTCAACGCCAGCAGGATCACCAGCAACAGGTTCAg                                                                                                        <  1:1117507/150‑1 (MQ=255)
                                                ccATGCACAGCACCGGTAGCGCCGTCAGCAAGGCCGCCTGCTGGTAACCCAGGCCGGTGCTGGCGCGCATGGGTTCGAGTAACGGGCCGATACTGGTCAGGATCGGCCGCAGGTTCAACGCcagcag                                                                                                                             <  1:97324/127‑1 (MQ=255)
                                                ccATGCACAGCACCGGTAGCGCCGTCAGCAAGGCCGCCTGCTGGTAACCCAGGCCGGTGCTGGCGCGCATGGGTTCGAGTAACGGGCCGATACTGGTCAGGATCGGCCGCAGGTTCAACGCcagcag                                                                                                                             >  2:97324/1‑127 (MQ=255)
                                                      acaGCACCGGTAGCGCCGTCAGCAAGGCCGCCTGCTGGTAACCCAGGCCGGTGCTGGCGCGCATGGGTTCGAGTAAAAGCACGAAATTGGTCAGGAGCGGCCGCAGGTTAAACCCCCGCAGGATAACCAGCAACAGGTTCAGTACAgcgc                                                                                                >  1:21317/1‑150 (MQ=255)
                                                             cGGTAGCGCCGTCAGCAAGGCCGCCTGCTGGTAACCCAGGCCGGTGCTGGCGCGCATGGGTTCGAGTAACGGGCCGATACTGGTCAGGATCGGCCGCAGGTTCAACGCCAGCAGGATCACCa                                                                                                                     <  2:884052/122‑1 (MQ=255)
                                                             cGGTAGCGCCGTCAGCAAGGCCGCCTGCTGGTAACCCAGGCCGGTGCTGGCGCGCATGGGTTCGAGTAACGGGCCGATACTGGTCAGGATCGGCCGCAGGTTCAACGCCAGCAGGATCACCa                                                                                                                     >  1:884052/1‑122 (MQ=255)
                                                             cGGTAGCGCCGTCAGCAAGGCCGCCTGCTGGTAACCCAGGCCGGTGCTGGCGCGCATGGGTTCGAGTAACGGGCCGATACTGGTCAGGATCGGCCGCAGGTTCAACGCCAGCAGGATCACCAGCAACAGGTTCAGTACAGCGCTCATtgt                                                                                         <  1:211634/150‑1 (MQ=255)
                                                               gTAGCGCCGTCAGCAAGGCCGCCTGCTGGTAACCCAGGCCGGTGCTGGCGCGCATGGGTTCGAGTAACGGGCCGATACTGGTCAGGATCGGCCGCAGGTTCAACGCca                                                                                                                                 <  2:413424/108‑1 (MQ=255)
                                                               gTAGCGCCGTCAGCAAGGCCGCCTGCTGGTAACCCAGGCCGGTGCTGGCGCGCATGGGTTCGAGTAACGGGCCGATACTGGTCAGGATCGGCCGCAGGTTCAACGCca                                                                                                                                 >  1:413424/1‑108 (MQ=255)
                                                                atcagcCGTCAGCAAGGCCGCCTGCTGGTAACCCAGGCCGGTGCTGGCGCGCATGGGTTCGAGTAACGGGCCGATACTGGTCAGGATCGGCCGCAGGTTCAACGCCAGCAGGATCACCAGAAACAGGTTCAGTACAGCGCTCATTGTGAg                                                                                      >  1:347368/5‑150 (MQ=255)
                                                                   cgcCGTCAGCAAGGCCGCCTGCTGGTAACCCAGGCCGGTGCTGGCGCGCATGGGTTCGAGTAACGGGCCGATACTGGTCAGGATCGGCCGCAGGTTCAACGCCAGCAGGATCACCAGCAACAGGTTCAGTACAGCGCTCATTGTGAGTCg                                                                                   <  1:188993/150‑1 (MQ=255)
                                                                         cAGCAAGGCCGCCTGCTGGTAACCCAGGCCGGTGCTGGCGCGCATGGGTTCGAGTAACAGGCCGATACTGGTCAGGATCGGCCGCAGGTTCAACGCCAGCAGGATCACCAGCAACAGGTTCAGTACAGCGCTCATTGTGAGTCGGTtctc                                                                             <  1:1042858/150‑1 (MQ=255)
                                                                             aaGGCCGCCTGCTGGTAACCCAGGCCGGTGCTGGCGCGCATGGGTTCGAGTAACGGGCCGATACTGGTCAGGATCGGCCGCAGGTTCAACGCCAGCAGGATCACCAGCAACAGGTTCAGTACAGCGCTCATTGTGAGTCGGTTCTCTTGt                                                                         >  1:906090/1‑150 (MQ=255)
                                                                              aGGCCGCCTGCTGGTAACCCAGGCCGGTGCTGGCGCGCATGGGTTCGAGTAACGGGCCGATACTGGTCAGGATCGGCCGCAGGTTCAACGCCAGCAGGATCACCAGCAACAGGTTCAGTACAGCGCTCATTGTGAGTCGGTTCTCTTGTa                                                                        <  2:630176/150‑1 (MQ=255)
                                                                                 ccgccTGCTGGTAACCCAGGCCGGTGCTGGCGCGCATGGGTTCGAGTAACGGGCCGATACTGGTCAGGATCGGCCGCAGGTTCAACGCCAGCAGGATCACCAGCAACAGGTTCAGTACAGCGCTCATTGTGAGTCGGTTCTCTTGTAACa                                                                     >  2:957415/1‑150 (MQ=255)
                                                                                     ctgctgGTAACCCAGGCCGGTGCTGGCGCGCATGGGTTCGAGTAACGGGCCGATACTGGTCAGGATCGGCCGCAGGTTCAACGCCAGCAGGATCACCAGCAACAGGTTCAGTACAGCGCTCATTGTGAGTCGGTTCTCTTGTAACAGGaa                                                                 <  2:747116/150‑1 (MQ=255)
                                                                                      tgctgGTAACCCAGGCCGGTGCTGGCGCGCATGGGTTCGAGTAACGGGCCGATACTGGTCAGGATCGGCCGCAGGTTCAACGCCAGCAGGATCACCAGCAACAGGTTCAGTACAGCGCTCATTGTGAGTCGGTTCTCTTGTAACAGGAAc                                                                <  2:1099293/150‑1 (MQ=255)
                                                                                      tgctgGTAACCCAGGCCGGTGCTGGCGCGCATGGGTTCGAGTAACGGGCCGATACTGGTCAGGATCGGCCGCAGGTTCAACGCCAGCAGGATCACCAGCAACAGGTTCAGTACAGCGCTCATTGTGAGTCGGTTCTCTTGTAACAGGAAc                                                                <  2:21317/150‑1 (MQ=255)
                                                                                          ggTAACCCAGGCCGGTGCTGGCGCGCATGGGTTCGAGTAACGGGCCGATACTGGTCAGGATCGGCCGCAGGTTCAACGCCAGCAGGATCACCAGCAACAGGTTCAGTACAGCGCTCATTGTGAGTCGGTTCTCTTGTAACAGGAACTGGc                                                            <  2:488759/150‑1 (MQ=255)
                                                                                             aaCCCAGGCCGGTGCTGGCGCGCATGGGTTCGAGTAACGGGCCGATACTGGTCAGGATCGGCCGCAGGTTCAACGCCAGCAGGATCACCAGCAACAGGTTCAGTACAGCGCTCATTGTGAGTCGGTTCTCTTGTAACa                                                                     >  1:1018193/1‑138 (MQ=255)
                                                                                             aaCCCAGGCCGGTGCTGGCGCGCATGGGTTCGAGTAACGGGCCGATACTGGTCAGGATCGGCCGCAGGTTCAACGCCAGCAGGATCACCAGCAACAGGTTCAGTACAGCGCTCATTGTGAGTCGGTTCTCTTGTAACa                                                                     <  2:1018193/138‑1 (MQ=255)
                                                                                                  aGGCCGGTGCTGGCGCGCATGGGTTCGAGTAACGGGCCGATACTGGTCAGGATCGGCCGCAGGTTCAACGCCAGCAGGATCACCAGCAACAGGTTCAGTACAGCGCTCATTGTGAGTCGGTTCTCTTGTAACAGGAACTGGCCATTATcc                                                    <  2:906090/150‑1 (MQ=255)
                                                                                                    gCCGGTGCTGGCGCGCATGGGTTCGAGTAACGGGCCGATACTGGTCAGGATCGGCCGCAGGTTCAACGCCAGCAGGATCACCAGCAACAGGTTCAGTACAGCGCTCATTGTGAGTCGGTTCTCTTGTAACAGGAACTGGCCATTATCCcg                                                  <  1:832436/150‑1 (MQ=255)
                                                                                                     ccGGTGCTGGGTCGCATGGGTTCGAGTAACGGGCCCAAACTGCTAAAGCTACACACCAGACTCCACGCTTTCCGTATCAGCAGCTATTGAGTTAGTACAGCGCTTATTGTGTGTCGGTTCTCTTGTAACAGGAACTGGCCATTATcccag                                                 <  2:439226/150‑3 (MQ=255)
                                                                                                         tGCTGGCGCGCATGGGTTCGAGTAACGGGCCGATACTGGTCAGGATCGGCCGCAGGTTCAACGCCAGCAGGATCACCAGCAACAGGTTCAGTACAGCGCTCATTGTGAGTCGGt                                                                                 <  2:482320/114‑1 (MQ=255)
                                                                                                         tGCTGGCGCGCATGGGTTCGAGTAACGGGCCGATACTGGTCAGGATCGGCCGCAGGTTCAACGCCAGCAGGATCACCAGCAACAGGTTCAGTACAGCGCTCATTGTGAGTCGGt                                                                                 >  1:482320/1‑114 (MQ=255)
                                                                                                             ggCGCGCATGGGTTCGAGTAACGGGCCGATACTGGTCAGGATCGGCCGCAGGTTCAACGCCAGCAGGATCACCAGCAACAGGTTCAGTACAGCGCTCATTGTGAGTCGGTTCTCTTGTAACAGGAACTGGCCATTATCCCGGCTCGGctg                                         >  1:925402/1‑150 (MQ=255)
                                                                                                                      gggTTCGAGTAACGGGCCGATACTGGTCAGGATCGGCCGCAGGTTCAACGCCAGCAGGATCACCAGCAACAGGTTCAGTACAGCGCTCATTGTGAGTCGGTTCTCTTGTAACAGGAACTGGCCATTATCCCGGCTCGGCTGGCCTGGCAg                                <  1:476438/150‑1 (MQ=255)
                                                                                                                        gTTCGAGTAACGGGCCGATACTGGTCAGGATCGGCCGCAGGTTCAACGCCAGCAGGATCACCAGCAACAGGTTCAGTACAGCGCTCATTGTGAGTCGGTTCTCTTGTAACAGGAACTGGCCATTATCCCGGCTCGGCTGGCCTGGCAGaa                              >  2:679768/1‑150 (MQ=255)
                                                                                                                            gAGTAACGGGCCGATACTGGTCAGGATCGGCCGCAGGTTCAACGCCAGCAGGATCACCAGCAACAGGTTCAGTACAGCGCTCATTGTGAGTCGGTTCTCTTGTAACAGGAACTGCCCATTATCCCGGCTCGGCTGGCCTGGCAGAAATGa                          >  1:578089/1‑150 (MQ=255)
                                                                                                                                   cggCCGATCCTGACCAGTATCGGCCGCAGGTTCAACGCCAGCAGGATCACCAGCAACAGGTTCAGTACAGCGCTAATAGAAGGTCGGTTCTCTTGTAACAGGAACTGGCCATTATCCCGGCACGGATGGCCTGGCAGAAATGAAGAGAtt                   >  1:677463/2‑150 (MQ=255)
                                                                                                                                     gCCGATCCTGACCAGTATCGGCCGCAGGTTCAACGCCAGCAGGATCACCAGCAACAGGTTCAGTACAGCGCTCATtgt                                                                                         >  2:17523/1‑78 (MQ=38)
                                                                                                                                             tGGTCAGGATCGGCCGCAGGTTCAACGCCAGCAGGATCACCAGCAACAGGTTCAGTACAGCGCTCATTGTGAGTCGGTTCTCTTGTAACAGGAACTGGc                                                            <  1:855742/99‑1 (MQ=255)
                                                                                                                                             tGGTCAGGATCGGCCGCAGGTTCAACGCCAGCAGGATCACCAGCAACAGGTTCAGTACAGCGCTCATTGTGAGTCGGTTCTCTTGTAACAGGAACTGGc                                                            >  2:855742/1‑99 (MQ=255)
                                                                                                                                                  aGGATCGGCCGCAGGTTCAACGCCAGCAGGATCACCAGCAACAGGTTCAGTACAGCGCTCATTGTGAGTCGGTTCTCTTGTAACAGGAACTGGCCATTATCCCGGCTCGGCTGGCCTGGCAGAAATGAAGAGATTGGATGCAGGTCAGTg    <  2:242065/150‑1 (MQ=255)
                                                                                                                                                    gATCGGCCGCAGGTTCAACGCCAGCAACAACACCAGCAACAGGTTCAGTACAGCGCTCATTGTGAGTCGGTTCTCTTGTAACAGGAACTGGCCATTATCCCGGCTCGGCTGGCCTGGCAGAAATGAAGAGATTGGATGCAGGTCAGTGGc  >  1:823004/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
GCTCGCTGACCCAGCGTCGCAGCCACGGTTGCAACAGGGGCACCAGGCCCATGCACAGCACCGGTAGCGCCGTCAGCAAGGCCGCCTGCTGGTAACCCAGGCCGGTGCTGGCGCGCATGGGTTCGAGTAACGGGCCGATACTGGTCAGGATCGGCCGCAGGTTCAACGCCAGCAGGATCACCAGCAACAGGTTCAGTACAGCGCTCATTGTGAGTCGGTTCTCTTGTAACAGGAACTGGCCATTATCCCGGCTCGGCTGGCCTGGCAGAAATGAAGAGATTGGATGCAGGTCAGTGGC  >  NC_002947/4281260‑4281557

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑