Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,911,732 G→C 7.4% intergenic (+81/+47) nirD → / ← PP_1707 nitrite reductase/D‑isomer specific 2‑hydroxyacid dehydrogenase family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,911,7320GC 7.4% 544.2 / 5.0 176intergenic (+81/+47)nirD/PP_1707nitrite reductase/D‑isomer specific 2‑hydroxyacid dehydrogenase family protein
Reads supporting (aligned to +/- strand):  ref base G (71/89);  new base C (14/2);  total (85/91)
Fisher's exact test for biased strand distribution p-value = 1.11e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 178 total aligned reads displayed.

GCCGGGCAGAGGTTGCGGGTGTGGCCGGTGCGGTTGAAGGGGGAGGCGGTGGAGCTGGCGGTGGCCTGATACCCCTTTGCTGCCTGTGCCGGCCTCTTCGCGGGCAAGCCCGCTCCTACGGAATCTGCACCGCCTTGTAATACGGCGCCGACCCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCAGCACAAGCCGACTAGATAACCAACCCCTGCGGCGTGCGCCCGGCAAACACATCCACCAGGCTGGCCACCACCATCTCGCCCATGCGCGTGCGTGTCTGTACAGTGGCACTGGCG  >  NC_002947/1911583‑1911881
                                                                                                                                                      |                                                                                                                                                     
     gCAGAGGTTGCGGGTGTGGCCGGTGCGGTTGAAGGGGGAGGCGGTGGAGCTGGCGGTGGCCTGATACCCCTTTGCTGCCTGTGCCGGCCTCTTCGCGGGCAAGCCCCCTCCGACCGAATCAGCACCCCCGTGTAAAACGGCGCCCAcccg                                                                                                                                                  >  1:5126409/1‑149 (MQ=255)
       tgggcgagatggtggtggccagcctggtggatgtgtttgccgggcgcacgccgcaggggttggttatctagtcggcttgtgctGGCCTCTTCGCGGGCAAGCCCGCTCCTACGGAATCTGCACCGCCTTGTAATACGGCGCCGACCCtgt                                                                                                                                                <  1:4282665‑M1/67‑1 (MQ=255)
       agagGTTGCGGGTGTGGCCGGTGCGGTTGAAGGGGGAGGCGGTGGAGCTGGCGGTGGCCTGATACCCCTTTGCTGCCTGTGCCGGCCTCTTCGCGGGCCAGCCCCCTCCGACGGAAACTGCACCCCCAGGTAAAACGGCGCCCAcccagg                                                                                                                                                >  2:1635676/1‑147 (MQ=255)
        gagGTTGCGGGTGTGGCCGGTGCGGTTGAAGGGGGAGGCGGTGGAGCTGGCGGTGGCCTGATACCCCTTTGCTGCCTGTGCCGGCCTCTTCGCGGGCAAGCCCGCTCCTACGGAATCTGCACCGCCTTGTAATACGGCGCCGACCCtgtg                                                                                                                                               <  1:3954388/150‑1 (MQ=255)
           gTTGCGGGTGTGGCCGGTGCGGTTGAAGGGGGAGGCGGTGGAGCTGGCGGTGGCCTGATACCCCTTTGCTGCCTGTGCCGGCCTCTTCGCGGGCAAGCCCGCTCCTACGGAATCTGCACCGCCTTGTAATACGGCGCCGACCCTGTGGGa                                                                                                                                            <  1:2139108/150‑1 (MQ=255)
            ttGCGGGTGTGGCCGGTGCGGTTGAAGGGGGAGGCGGTGGAGCTGGCGGTGGCCTGATACCCCTTTGCTGCCTGTGCCGGCCTCTTCGCGGGGAAGCCCGCTCCTAAGGAAAAGGCACCGCCAAGTAAAAAGGCGCCGCCCCAGGGGGAg                                                                                                                                           >  2:2611863/1‑150 (MQ=255)
            ttGCGGGTGTGGCCGGTGCGGTTGAAGGGGGAGGCGGTGGAGCTGGCGGTGGCCTGATACCCCTTTGCTGCCTGTGCCGGCCTCTTCGCGGGCAAACCCCCGCCAAACGGAACAGGACCGCCTTGTAATAAGGCGCCCACCCGGGGGGAg                                                                                                                                           >  2:38028/1‑150 (MQ=255)
            ttGCGGGTGTGGCCGGTGCGGTTGAAGGGGGAGGCGGTGGAGCTGGCGGTGGCCTGATACCCCTTTGCTGCCTGTGCCGGCCACTTCGCGGGCAAGCCCCCACCGAACGAAAAAGCACCGCCCTGTAATACCGCGCCCACCCCGTGGGAg                                                                                                                                           >  2:38030/1‑150 (MQ=255)
            gagatggtggtggccagcctggtggatgtgtttgccgggcgcacgccgcaggggttggttatctagtcggcttgtgctGGCCTCTTCGCGGGCAAGCCCGCTCCTACGGAATCTGCACCGCCTTGTAATACGGCGCCGACCCTGTGGGAg                                                                                                                                           >  1:5143496‑M1/79‑150 (MQ=255)
              gatggtggtggccagcctggtggatgtgtttgccgggcgcacgccgcaggggttggttatctagtcggcttgtgctGGCCTCTTCGCGGGCAAGCCCGCTCCTACGGAATCTGCACCGCCTTGTAATACGGCGCCGACCCTGTGGGAGCg                                                                                                                                         >  2:2245221‑M1/77‑150 (MQ=255)
              gCGGGTGTGGCCGGTGCGGTTGAAGGGGGAGGCGGTGGAGCTGGCGGTGGCCTGATACCCCTTTGCTGCCTGTGCCGGCCTCTTCGCGGGCAAGCCCGCTCCTACGGAATCTGCACCGCCTTGTAATACGGCGCCGACCCTGTGGGAGCg                                                                                                                                         <  2:5146431/150‑1 (MQ=255)
              gCGGGTGTGGCCGGTGCGGTTGAAGGGGGAGGCGGTGGAGCTGGCGGTGGCCTGATACCCCTTTGCTGCCTGTGCCGGCCACTTCGCGGGCAAGCCCGCTCCTACGGAATCTGCACCGCCTTGTAATACGGCGCCGACCCTGTGGGAGCg                                                                                                                                         <  2:4268093/150‑1 (MQ=255)
               cGGGTGTGGCCGGTGCGGTTGAAGGGGGAGGCGGTGGAGCTGGCGGTGGCCTGATACCCCTTTGCTGCCTGTGCCGGCCTCTTCGCGGGCAAGCCCGCTCCTACGGAATCTGCACCGCCTTGTAATACGGCGCCGACCCTGTGGGAGCgg                                                                                                                                        <  1:1872344/150‑1 (MQ=255)
               cGGGTGTGGCCGGTGCGGTTGAAGGGGGAGGCGGTGGAGCTGGCGGTGGCCTGATACCCCTTTGCTGCCTGTGCCGGCCTCTTCGCGGGCAAGCCCGCTCCTACGGAATCTGCACCGCCTTGTAATACGGCGCCGACCCTGTGGGAGCgg                                                                                                                                        <  1:3421319/150‑1 (MQ=255)
               cGGGTGTGGCCGGTGCGGTTGAAGGGGGAGGCGGTGGAGCTGGCGGTGGCCTGATACCCCTTTGCGGCCTGTGCCGGCCTCTTCGCGGGCAAGCCCGCTCCTACGGAATCTGCACCGCCTTGTAATACGGCGCCGACCCTGTGGGAGCgg                                                                                                                                        <  2:4899253/150‑1 (MQ=255)
               atggtggtggccagcctggtggatgtgtttgccgggcgcacgccgcaggggttggttatctagtcggcttgtgctGGCCTCTTCGCGGGCAAGCCCGCTCCTACGGAATCTGCACCGCCTTGTAATACGGCGCCGACCCTGTGGGAGCgg                                                                                                                                        >  2:2404394‑M1/76‑150 (MQ=255)
               atggtggtggccagcctggtggatgtgtttgccgggcgcacgccgcaggggttggttatctagtcggcttgtgctGGCCTCTTCGCGGGCAAGCCCGCTCCTACGGAATCTGCACCGCCTTGTAATACGGCGCCGACCCTGTGGGAGCgg                                                                                                                                        >  2:2404393‑M1/76‑150 (MQ=255)
                tggtggtggccagcctggtggacgtgtttgccgggcgcacgccgcaggggttggttatctagtcggcttgtgctGGCCTCTTCGCGGGCAAGCCCGCTCCTACGGAATCTGCACCGCCTTGTAATACGGCGCCGACCCTGTGGGAGCggg                                                                                                                                       >  2:1986646‑M1/75‑150 (MQ=255)
                gggTGTGGCCGGTGCGGTTGAAGGGGGAGGCGGTGGAGCTGGCGGTGGCCTGATACCCCTTTGCTGCCTGTGCCGGCCTCTTCGCGGGCAAGCCCGCTCCTACGGAATCTGCACCGCCTTGTAATACGGCGCCGACCCTGTGGGAGCggg                                                                                                                                       <  1:2183097/150‑1 (MQ=255)
                 ggtggtggccagcctggtggatgtgtttgccgggcgcacgccgcaggggttggttatctagtcggcttgtgctGGCCTCTTCGCGGGCAAGCCCGCTCCTACGGAATCTGCACCGCCTTGTAATACGGCGCCGACCCTGTGGGAGCGGGc                                                                                                                                      >  1:3861238‑M1/74‑150 (MQ=255)
                 ggtggtggccagcctggtggatgtgtttgccgggcgcacgccgcaggggttggttatctagtcggcctgtgccGGCCTCTTCGCGGGCAAGCCCGCTCCTACGGAATCTGCACCGCCTTGTAATACGGCGCCGACCCTGTGGGAGCGGGc                                                                                                                                      >  2:2755105‑M1/74‑150 (MQ=255)
                 ggTGTGGCCGGTGCGGTTGAAGGGGGAGGCGGTGGAGCTGGCGGTGGCCTGATACCCCTTTGCTGCCTGTGCCGGCCTCTTCGCGGGCAAGCCCGCTCCTACGGAATCTGCACCGCCTTGTAATACGGCGCCGACCCTGTGGGAGCGGGc                                                                                                                                      <  2:4446183/150‑1 (MQ=255)
                    gtgGCCGGTGCGGTTGAAGGGGGAGGCGGTGGAGCTGGCGGTGGCCTGATACCCCTTTGCTGCCTGTGCCGGCCTCTTCGCGGGCAAGCCCGCTCCTACGGAATCTGCACCGCCTTGTAATACGGCGCCGACCCTGTGGGAGCGGGCTTg                                                                                                                                   <  2:698164/150‑1 (MQ=255)
                    ggtggccagcctggtggatgtgtttgccgggcgcacgccgcaggggttggttatctagtcggcctgtgccGGCCTCTTCGCGGGCAAGCCCGCTCCTACGGAATCTGCACCGCCTTGTAATACGGCGCCGACCCTGTGGGAGCGGGCTTg                                                                                                                                   >  2:4007004‑M1/71‑150 (MQ=255)
                       gCCGGTGCGGTTGAAGGGGGAGGCGGTGGAGCTGGCGGTGGCCTGATACCCCTTTGCTGCCTGTGCCGGCCTCTTCGCGGGCCAGCCCCCACCTACGGAATAAGCACCGCCGTGTAATAAGGCGCCGACCCGGTGGGAGCGGGCTTGccc                                                                                                                                >  1:4143782/1‑150 (MQ=255)
                       gCCGGTGCGGTTGAAGGGGGAGGCGGTGGAGCTGGCGGTGGCCTGATACCCCTTTGCTGCCTGTGCCGGCCTCTTCGCGGGCAAGCCCGCTCCTACGGAATCTGCACCGCCTTGTAATACGGCGCCGACCCTGTGGGAGCGGGCTTGccc                                                                                                                                <  1:4424102/150‑1 (MQ=255)
                       gCCGGTGCGGTTGAAGGGGGAGGCGGTGGAGCTGGCGGTGGCCTGATACCCCTTTGCTGCCTGTGCCGGCCTCTTCGCGGGCAAGCCCGCGCCTACGGAATCTGCACCGCCTTGTAATACGGCGCCGACCCAGTGGGAGCGGGCTTGccc                                                                                                                                >  2:3624803/1‑150 (MQ=255)
                          ggTGCGGTTGAAGGGGGAGGCGGTGGAGCTGGCGGTGGCCTGATACCCCTTTGCTGCCTGTGCCGGCCTCTTCGCGGGCAAGCCCGCACCAACCGAAAAAGCACCGCCGTGTAATACGGCCCCCCCCCAGTGGGCGCGGGCGTGcccccc                                                                                                                             >  1:5076414/1‑147 (MQ=255)
                            tGCGGTTGTTGGGGGAGGGGGTGGAGGTGGCGGTGGCCTGATACCCCTTTGCTGCCTGTGCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACGGAATCTGCACCGCCTTGTAATACGGCGCCGACCCTGTGGGAGCGGGCTTGCCCGCGaa                                                                                                                           <  1:1276197/150‑1 (MQ=39)
                            tGCGGTTGAAGGGGGAGGCGGTGTAGCTGGCGGTGGCCTGATACCCCTTTGCTGCCTGTGCCGGCCTCTTCGCGGGCAAGCCCGCTCCTACGGAATCTGCACCGCCTTGTAATACGGCGCCGACCCTGTGGGAGCGGGCTTGCCCGCGaa                                                                                                                           <  2:771494/150‑1 (MQ=255)
                            tGCGGTTGAAGGGGGAGGCGGTGGAGCTGGCGGTGGCCTGATACCCCTTTGCTGCCTGTGCCGGCCTCTTCGCGGGCAAGCCCGCTCCTACGGAATCTGCACCGCCTTGTAATACGGCGCCGACCCTGTGGGAGCGGGCTTGCCCGCGaa                                                                                                                           <  1:2806065/150‑1 (MQ=255)
                            tGCGGTTGAAGGGGGAGGCGGTGGAGCTGGCGGTGGCCTGATACCCCTTTGCTGCCTGTGCCGGCCTCTTCGCGGGCAAGCCCGCTCCTACGGAATCTGCACCGCCTTGTAATACGGCGCCGACCCTGTGGGAGCGGGCAAGCCcgctcc                                                                                                                           <  2:5094079/150‑4 (MQ=255)
                              cGGTTGAAGGGGGAGGCGGTGGAGCTGGCGGTGGCCTGATACCCCTTTGCTGCCTGTGCCGGCCTCTTCGCGGGCAAGCCCGCTCCTACGTAATCTGCACCGCCTTGTAATACGGCGCCGACCCTGTGGGAGCGGGCTTGCCCGCGAaga                                                                                                                         <  2:1452072/150‑1 (MQ=255)
                                gTTGAAGGGGGAGGCGGTGGAGCTGGCGGTGGCCTGATACCCCTTTGCTGCCTGTGCCGGCCTCTTCGCGGGCCAGCCCGCGCCTACGGAATCTGCGCCGCCTTGTAATACGGCGCCCACCCGGTGGGAGCGGGGGTGCCCGCGAAGAgg                                                                                                                       >  2:4347288/1‑150 (MQ=255)
                                gTTGAAGGGGGAGGCGGTGGAGCTGGCGGTGGCCTGATACCCCTTTGCTGCCTGTGCCGGCCTCTTCGCGGGCAATCCCGCAACGAAGGCATATTCACCGCCATGTAATAACGCGCCGACCCAGTGGGAGCGGGCGTTGCCGCGAAGAgg                                                                                                                       >  2:1141391/1‑150 (MQ=255)
                                  tGAAGGGGGTGGGGGTGGCGGTGGGGGTGGCCTGACACCCCTTTGCTGCCTGTGCGGGCCACTTCGCGGGCAAGCCCGCTCCTACGGAATCTGCACCGCCTTGTAATACGGCGCCGACCCTGTGGGAGGGGGCTTGCCCGCGAAGAGGcc                                                                                                                     <  2:2363680/150‑1 (MQ=38)
                                  tGAAGGGGGAGGCGGTGGAGCTGGCGGTGGCCTGATACCCCTTTGCTGCCTGTGCCGGCCTCTTCGCGGGCAAGCCCGCTCCTACGGAATCTGCACCGCCTTGTAATACGGCGCCGACCCTGTGGGAGCGGGCTTGCCCGCGAAGAGGcc                                                                                                                     <  2:2696311/150‑1 (MQ=255)
                                      gggggCGGGGGTGGAGCGGGGGGGGGCCTGATACCCCTTTGCGGCCTGGGCCGGCCTCTTCGCGGGCACGCCCGCTCCTCCGGAATCTGCACCGCCTTGTAATCCGGCGCCGACCCTGTGGGAGCGGGCTGGCCCGCGAAGAGGCCAGca                                                                                                                 <  2:1685014/150‑1 (MQ=37)
                                       ggggAGGCGGTGGAGCTGGCGGTGGCCTGATACCCCTTTGCTGCCTGTGCCGGCCTCTTCGCGGGCAAGCCCGCTCCTACGGTATGTGCAGCGCCTTGTAATAAGGCGCCGACCCTGTGGGTGCGGGCTTGGCCGCGAAGAGGCCAGcac                                                                                                                >  2:2167305/1‑150 (MQ=39)
                                         gtttgccgggcgcacgccgcaggggttggttatctagtcggcttgtgccGGCCTCTTCGCGGGCAAGCCCGCTCCTACGGAATCTGCACCGCCTTGTAATACGGCGCCGACCCTGTGGGAGCGGGc                                                                                                                                      <  2:2598955‑M1/77‑1 (MQ=255)
                                         gtttgccgggcgcacgccgcaggggttggttatctagtcggcttgtgccGGCCTCTTCGCGGGCAAGCCCGCTCCTACGGAATCTGCACCGCCTTGTAATACGGCGCCGACCCTGTGGGAGCGGGc                                                                                                                                      >  1:2599615‑M1/50‑126 (MQ=255)
                                         ggAGGCGGTGGAGCTGGCGGTGGCCTGATACCCCTTTGCTGCCTGTGCCGGCCTCTTCGCGGGCAAGCCCGCTCCTACGGAATCTGCACCGCCTTGTAATACGGCGCCGACCCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCAGCACaa                                                                                                              <  2:3214155/150‑1 (MQ=255)
                                          gAGGCGGTGGAGCTGGCGGTGGCCTGATACCCCTTTGCTGCCTGTGCCGGCCTCTTCGCGGGCAAGCCCGCTCCTACGGAATCTGCACCGCCGTGTAATACGGCGCCGACCCGGGGGGAGCGGGCTTGCCCGCGAAGGGGCCAGCACAAg                                                                                                             >  2:391872/1‑150 (MQ=255)
                                            ggCGGTGGAGCTGGCGGTGGCCTGATACCCCTTTGCTGCCTGTGCCGGCCTCTTCGCGGGCAAGCCCGCTCCTACGGAATCTGCACCGCCTTGTAATACGGCGCCGACCCGGTGGGAGCGGGCTTGCCCGCGAAGAGGCCCGCACAAGcc                                                                                                           >  2:1752958/1‑150 (MQ=255)
                                              ccgggcgcacgccgccggggttggttatctagtcggcctgtgccGGCCTCTTCGCGGGCAAGCCCGCTCCTACGGAATCTGCACCGCCTTGTAATACGGCGCCGACCCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCAGCACAAGCCGa                                                                                                         >  2:303246‑M1/45‑150 (MQ=255)
                                               ggTGGAGCTGGCGGTGGCCTGATACCCCTTTGCTGCCTGTGCCGGCCTCTTCGCGGGCAAGCCCGCTCCTACGGAATCTGCACCGCCTTGTAATACGGCGCCGACCCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCAGCACAAGCCGAc                                                                                                        <  2:4032106/150‑1 (MQ=255)
                                                 tGGGGCTGGCGGTGGCCTGATACCCCTTTGCTGCCTGTGCCGGCCTCTTCGCGGGCAAGCCCCCGCCTACGGCATATGCACCGCCGTGTAAAAACGCGCCCACCCCGTGGGAGCGGGCGTGCCCCCGAAGAGGCCCGCACAAGCCCAcaa                                                                                                      >  2:5175034/1‑148 (MQ=255)
                                                 ggcgcacgccgcaggggttggttatctagtcggcctgtgccGGCCTCTTCGCGGGCAAGCCCGCTCCTACGGAATCTGCACCGCCTTGTAATACGGCGCCGACCCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCAGCACAAGCCGACTa                                                                                                      >  2:4604013‑M1/42‑150 (MQ=255)
                                                       tGGCGGTGGCCTGATACCCCTTTGCTGCCTGTGCCGGCCTCTTCGCGGGCAAGCCCGCTCCTACGGAATCTGCACCGCCTTGTAATACGGCGCCGACCCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCAGCACAAGCCGACTAGATaac                                                                                                <  2:3748358/150‑1 (MQ=255)
                                                                  tGATACCCCTTTGCTGCCTGTGCCGGCCTCGTCGCGGGCCATCCCCCTCCTACGGGATCTTCACCGCCGTGTAAAACGGCGCCCGCCCGGTGGGAGCGGGCGTTCCCGCGCAGAGGCCAGCACAAACCGACTAGCTAACCAACCCCAgcg                                                                                     >  2:2932877/1‑150 (MQ=255)
                                                                        cccTTTGCTGCCTGTGCCGGCCTCTTCGCGGGCAAGCCCCCTCCTACGGAATCAGCACCGCCGTGTAATACGGCGCCGACCCTGTGGGAGCGGGCTTGCCCGCCAAGAGGCCAGCACAAGCCGACAAGATAAACAACCCCAGCGGCGTgc                                                                               >  2:5084828/1‑150 (MQ=38)
                                                                          cTTTGCTGCCTGTGCCGGCCTCTTCGCGGGCAAGCCCGCTCCTACGGAATCTGCACCGCCTTGTAATACGGCGCCGACCCt                                                                                                                                                  <  1:705870/81‑1 (MQ=255)
                                                                           tttGCTGCCTGTGCCGGCCTCGTCGCGGGCCAGCCCCCGCCGACCGCAACTGCCCCCCCCGGGAAAAACGCGCCCCCCCCGGGGGAGCCGGCGTGCCCCCGCACAGGCCCCCCCAACCCCCCAAGAGAACCCACCCCCGCGGCGGGCGcc                                                                            >  2:1053250/1‑150 (MQ=255)
                                                                            gggctgcCTGTGCCGGCCTCTTCGCGGGCTAGCCCGCTCCTCCGGAATGTGCACCGCCTTGTAATACGGCGCCGACCCGGTGGGAGCGGGCTTGCCCGCGAAGAGGCCAGCACAAGCCGACTAGATAACCAACCCCTGCGGCGTGCGccc                                                                           <  2:4207741/148‑1 (MQ=255)
                                                                             tgctgcCTGTGCCGGCCTCTTCGCGGGCAAGCCCGCTCCTACGGAATCTGCACCGCCTTGTAATACGGCGCCGACCCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCAGCACAAGCCGACTAGATAACCAACCCCTGCGGCGTGCGCCCg                                                                          <  2:4403169/150‑1 (MQ=255)
                                                                             tgctgcCTGTGCCGGCCGCTGGGCGGCCAAGCCCGCTCCTACGGAATCTGCACCGCCTTGTAATACGGCGCCGACCCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCAGCACAAGCCGACTAGATAACCAACCCCTGCGGCGTGCGCCCg                                                                          <  2:5415649/150‑1 (MQ=255)
                                                                             tgctgcCTGGGCGGGCCCCTCCGCGGGCAACCCCGCTCCTCCGGAATCGGCACCGCCTTGTAAGCGGGCGCCGACCCTGTGGGCGCGGGCTGCCCCGCGAAGAGGCCAGCACAAGCCGACTAGATAACCAACCCCTGCGGCGTGCGCCCg                                                                          <  1:867669/150‑1 (MQ=255)
                                                                                  ccTGGGCCGGCCTCTTCGGGGGCAAGCCCGCACCTCCGGAATCTGCACCGCCTTGTAATCCGGCGCCGACCCTGTGGGCGCGGGCTGGCCCGCGAAGAGGCCAGCACAAGCCGACTAGATAACCAACCCCTGCGGCGTGCGCCCGGCaaa                                                                     <  1:1761245/150‑1 (MQ=255)
                                                                                   cTGTGCCGGCCTCTTCGCGGGCCAGCCCGCTCCTCCGGAATCTGCACCGCCTTGTAAGACGGCGCCGACCCTGGGGGCGCGGGCTTGCCCGCGAAGAGGCCAGCACAAGCCGACTAGATAACCAACCCCTGCGGCGTGCGCCCGGCAAac                                                                    <  2:2850262/150‑1 (MQ=39)
                                                                                   cTGTGCCGGCCTCTTCGCGGGCAAGCCCGCTCCTACGGAATCTGCACCGCCTTGTAATACGGCGCCGACCCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCAGCACAAGCCGACTAGATAACCAACCCCTGCGGCGTGCGCCCGGCAAac                                                                    <  2:829401/150‑1 (MQ=255)
                                                                                         gGGCCTCTTCGCGGGCAAGCCCGCTCCTACGGTATATGCACCGCCTTGTATTAGGGCGCCGACCGTGGGGGTGCGGGCTGGCCCGCGAAGAGGCCAGCACAAGCCGACTAGATAACCAACCCCTGCGGCGTGCGCCCGGCAATCACATcc                                                              >  2:1276384‑M1/2‑150 (MQ=255)
                                                                                             CTCTTCGCGGGCAAGCCCGCTCCTACGGAATCTGCACCGCCTTGTAATACGGCGCCGACCCTGTGGGAGCGGGCTTGCCCGCGCAGAGGCCAGCACAAGCCGACTAGATAACCAACCCCTGCGGCGTGCGCCCGGCAAACACATccacca                                                          <  2:5501805‑M1/150‑1 (MQ=255)
                                                                                               CTTCGCGGGCACGCCCGCTCCTACGTAATCTGCACCGCCTTTTATTACGGCGCCGACCCTGTGGGAGGGGGCTTGCCCGCGAAGAGGCCAGCACAAGCCGACTAGATAACCAACCCCTGCGGCGTACGCCCGGCAAACACATCCACCAgg                                                        >  1:2762066‑M1/1‑150 (MQ=255)
                                                                                                     GGGCAAGCCCGCTCCTACGGAATCTGCACCGCCTTGTAATACGGCGCCGACCCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCAGCACAAGCCGACTAGATAACCAACCCCTGCGGCGTGCGCCCGGCAAACACATCCACCAGGCTGGcc                                                  <  2:1520411‑M1/150‑1 (MQ=255)
                                                                                                     GGGCAAGCCCGCTCCTACGGAATCTGCACCGCCTTGTAATACGGCGCCGACCCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCAGCACAAGCCGACTAGATAACCAACCCCTGCGGCGTGCGCCCGGCAAACACATCCACCAGGCTGGcc                                                  <  2:4356581‑M1/150‑1 (MQ=255)
                                                                                                      GGCAAGCCCGCTCCTCCGGAATCTGCACCGCCTTGTAATCCGGCGCCGACCCTGTGGGCGCGGGCTTGCCCGCGAAGAGGCCAGCACAAGCCGACTAGATAACCAACCCCTGCGGCGTGCGCCCGGCAAACACATCCACCAGGCTGGcca                                                 >  2:3040426‑M1/1‑150 (MQ=255)
                                                                                                      GGCAAGCCCGCTCCTACGGAATCTGCACCGCCTTGTAATACGGCGCCGACCCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCAGCACAAGCCGACTAGATAACCAACCCCTGCGGCGTGCGCCCGGCAAACACATCCACCAGGCTGGcca                                                 <  2:3217883‑M1/150‑1 (MQ=255)
                                                                                                        CAAGCCCGCTCCTACGGAATCTGCACCGCCTTGTAATACGGCGCCGACCCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCAGCACAAGCCGACTAGATAACCAACCCCTGCGGCGTGCGCCCGGCAAACACATCCACCAGGCTGGccacc                                               <  1:2505064‑M1/150‑1 (MQ=255)
                                                                                                        CAAGCCCGCTCCTACGGAATCTGCACCGCCTTGTAATACGGCGCCGACCCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCAGCACAAGCCGACTAGATAACCAACCCCTGCGGCGTGCGCCCGGCAAACACATCCACCAGGCTGGccacc                                               <  1:2505320‑M1/150‑1 (MQ=255)
                                                                                                         AAGCCCGCTCCTACGGAATCTGCACCGCCTTGTAATACGGCGCCGACCCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCAGCACAAGCCGACTAGATAACCAACCCCTGCGGCGTGCGCCCGGCAAACACATCCACCAGGCTGGccacca                                              <  2:5433155‑M1/150‑1 (MQ=255)
                                                                                                            CCCGCTCCTACGGAATCTGCACCGCCTTGTAATACGGCGCCGACCCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCAGCACAAGCCGACTAGATAACCAACCCCTGCGGCGTGCGCCCGGCAAACACATCCACCAGGCTGGccaccacca                                           <  2:3322552‑M1/150‑1 (MQ=255)
                                                                                                            CCCGCTCCTACGGAAACTGCACCGCCTTGAAAACCGGCGCCGACCCCGTGGGAGCGGGCTGGCCCGCGAAGAGGCCAGCACAAGCCGACTAGATAACCAACCCCTGCGGCGTGCGCCCGGCAAACACATCCACCAGGCTGGccaccacca                                           >  1:1144861‑M1/1‑150 (MQ=255)
                                                                                                             CCGCTCCTACGGAATCTGCACCGCCTTGTAATACGGCGCCGACCCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCAGCACAAGCCGACTAGATAACCAACCCCTGCGGCGTGCGCCCGGCAAACACATCCACCAGGCTGGCCACCACCAt                                          <  2:2382474‑M1/150‑1 (MQ=255)
                                                                                                             CCGCTCCTACGGAATCTGCACCGCCTTGTAATACGGCGCCGACCCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCAGCACAAGCCGACTAGATAACCAACCCCTGCGGCGTGCGCCCGGCAAACACATCCACCAGGCTGGCCACCACCAt                                          <  2:3889762‑M1/150‑1 (MQ=255)
                                                                                                              CGCTCCTACGGAATCTGCACCGCCTTGTAATACGGCGCCGACCCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCAGCACAAGCCGACTAGATAACCAACCCCTGCGGCGTGCGCCCGGCAAACACATCCACCAGGCTGGCCACCACCAtc                                         >  1:3217883‑M1/1‑150 (MQ=255)
                                                                                                                   CTACGGAATCTGCACCGCCTTGTCATACGGCGCCGACCCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCAGCACAAGCCGACTAGATAACCAACCCCTGCGGCGTGCGCCCGGCAAACACATCCACCAGGCTGGCCACCACCATCTCGcc                                    <  1:4507816‑M1/150‑1 (MQ=255)
                                                                                                                   CTACGGAATCTGCACCGCCTTGTAATACGGCGCCGACCCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCAGCACAAGCCGACTAGATAACCAACCCCTGCGGCGTGCGCCCGGCAAACACATCCACCAGGCTGGCCACCACCATCTCGcc                                    <  1:1721888‑M1/150‑1 (MQ=255)
                                                                                                                      aCGGAATCTGCACCGCCTTGTAATACGGCGCCGACCCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCAGCACAAGCCGACTAGATAACCAACCCCTGCGGCGTGCGCCCGGCAAACACATCCACCAGGCTGGCCACCACCATCTCGCCCa                                  >  2:2436780/1‑150 (MQ=255)
                                                                                                                       cGGAATCTGCACCGCCTTGTAATACGGCGCCGACCCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCAGCACAAGCCGACTAGATAACCAACCCCTGCGGCGTGCGCCCGGCAAACACATCCACCAGGCTGGCCACCACCATCTCGCCCAt                                 >  1:902757/1‑150 (MQ=255)
                                                                                                                         gAATCTGCACCGCCTTGTAATACGGCGCCGACCCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCAGCACAAGCCGACTAGATAACCAACCCCTGCGGCGTGCGCCCGGCAAACACATCCACCAGGCTGGCCACCACCATCTCGCCCATgc                               >  1:5528592/1‑150 (MQ=255)
                                                                                                                         gAATCTGCACCGCCTTGTAATACGGCGCCGACCCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCAGCACAAGCCGACTAGATAACCAACCCCTGCGGCGTGCGCCCGGCAAACACATCCACCAGGCTGGCCACCACCATCTCGCCCATgc                               >  1:2707259/1‑150 (MQ=255)
                                                                                                                              tGCACCGCCTTGTAATACGGCGCCGACCCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCAGCACAAGCCGACTAGATAACCAACCCCTGCGGCGTGCGCCCGGCAAACACATCCACCAGGCTGGCCACCACCATCTCGCCCATGCgcgtg                          >  1:481744/1‑150 (MQ=255)
                                                                                                                               gccccGCCTTTTAATCCGGCGCCCCCCCGGGGGGCGGGGGCGTGCCCGCGAAGAGGCCAGCACAAGCCGACTAGATAACCAACCCCTGCGGCGTGCGCCCGGCAAACACATCCACCAGGCTGGCCACCACCATCTCGCCCATGCgcgtgc                         <  1:2500911/147‑1 (MQ=255)
                                                                                                                               gCACCGCCTTGTAATACGGCGCCGACCCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCAGCACAAGCCGACTAGATAACCAACCCCTgcg                                                                                     >  1:320255/1‑90 (MQ=255)
                                                                                                                                 aCCGCCTTGTAATACGGCGCCGACCCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCAGCACAAGCCGACTAGATAACCAACCCCTGCGGCGTGCGCCCGGCAAACACATCCACCAGGCTGGCCACCACCATCTCGCCCATGCgcgtgcgt                       <  1:5599550/150‑1 (MQ=255)
                                                                                                                                  ccGCCTTGTAATACGGCGCCGACCCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCAGCACAAGCCGACTAGATAACCAACCCCTGCGGCGTGCGCCCGGCAAACACATCCACCAGGCTGGCCACCACCATCTCGCCCATGCGcgtgcgtg                      >  2:4130034/1‑150 (MQ=255)
                                                                                                                                   cGCCTTGTAATACGGCGCCGACCCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCAGCACAAGCCGACTAGATAACCAACCCCTGCGGCGTGCGCCCGGCAAACACATCCACCAGGCTGGCCACCACCATCTCGCCCATGCGCGTGCgtgt                     <  2:4507805/150‑1 (MQ=255)
                                                                                                                                   cGCCTTGTAATACGGCGCCGACCCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCAGCACAAGCCGACTAGATAACCAACCCCTGCGGCGTGCGCCCGGCAAACACATCCACCAGGCTGGCCACCACCATCTCGCCCATGCGCGTGCgtgt                     >  2:4688514/1‑150 (MQ=255)
                                                                                                                                   cGCCTTGTAATACGGCGCCGACCCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCAGCACAAGCCGACTAGATAACCAACCCCTGCGGCGTGCGCCCGGCAAACACATCCACCAGGCTGGCCACCACCATCTCGCCCATGCGCGTGCgtgt                     >  1:486250/1‑150 (MQ=255)
                                                                                                                                   cGCCTTGTAATACGGCGCCGACCCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCAGCACAAGCCGACTAGATAACCAACCCCTGCGGCGTGCGCCCGGCAAACACATCCACCAGGCTGGCCACCACCATCTCGCCCATGCGCGTGCgtgt                     >  2:211056/1‑150 (MQ=255)
                                                                                                                                     ccTTTTAATCGGGCGCCGCCCCTGGGGGCGGGGGCTGGCCCGCGAAGAGGCCAGCACAAGCCGACTAGATAACCAACCCCTGCGGCGTGCGCCCGGCAAACACATCCACCAGGCTGGCCACCACCATCTCGCCCATGCGCGTGCGTGTCt                   <  1:1635676/150‑1 (MQ=255)
                                                                                                                                       ttGTAATCCGGCGCCGCCCCGGGGGGAGCGGGCTGGCCCGCGAAGAGGCCAGCACAAGCCGACTAGATAACCAACCCCTGCGGCGTGCGCCCGGCAAACACATCCACCAGGCTGGCCACCACCATCTCGCCCATGCGCGTGCGTGTCTGt                 <  2:1233486/150‑1 (MQ=255)
                                                                                                                                       ttGTAATACGGCGCCGACCCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCAGCACAAGCCGACTAGATAACCAACCCCTGCGGCGTGCGCCCGGCAAACACATCCACCAGGCTGGCCACCACCATCTCGCCCATGCGCGTGCGTGTCTGt                 >  2:4961366/1‑150 (MQ=255)
                                                                                                                                         gTAATACGGCGCCGACCCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCAGCACAAGCCGACTAGATAACCAACCCCTGCGGCGTGCGCCCGGCAAACACATCCACCAGGCTGGCCACCACCATCTCGCCCATGCGCGTGCGTGTCTGTAc               >  1:343875/1‑150 (MQ=255)
                                                                                                                                            aTACGGCGCCGACCCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCAGCACAAGCCGACTAGATAACCAACCCCTGCGGCGTGCGCCCGGCAAACACATCCACCAGGCTGGCCACCACCATCTCGCCCATGCGCGTGCGTGTCTGTACAGt            <  1:3949055/150‑1 (MQ=255)
                                                                                                                                              aCGGCGCCGACCCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCAGCACAAGCCGACTAGATAACCAACCCCTGCGGCGTGCGCCCGGCAAACACATCCACCAGGCTGGCCACCACCATCTCGCCCATGCGCGTGCGTGTCTGTACAGTgg          >  1:1220994/1‑150 (MQ=255)
                                                                                                                                                ggCGCCGACCCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCAGAACAAGCCCACTAGATAACCAACCCCTGCGGCGTGCGCCCGGCAAACACATCAACCAGGCTGGCCACCACCATCTCGCCCATGCGCGTGCGAGTCTCTACAGTGGcc        >  2:3525982/1‑149 (MQ=39)
                                                                                                                                                 gcgcCGACCCTGTGGGTGGGGGGTTGCCCGCGAAGTGGCCAGCACAAGCCGACTAGATAACCAACCCCTGCGGCGTGCGCCCGGCAAACTCATCCACCAGGCTGGCCTCCTGGTTCTCGCCCATGCGCGTGCGTGTCTGTACAGTGGCAc       <  2:1717910/150‑1 (MQ=255)
                                                                                                                                                     cGACCCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCAGCACAAGCCGACTAGATAACCAACCCCTGCGGCGTGCGCCCGGCAAACACATCCACCAGGCTGGCCACCACCATCTCGCCCATGCGCGTGCGTGTCTGTACAGTGGCACTGgc   >  2:3588150/1‑150 (MQ=255)
                                                                                                                                                      gACCCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCAGCACAAGCCGACTAGATAACCAACCCCTGCGGCGTGCGCCCGGCAAACACATCCACCAGGCTGGCCACCACCATCTCGCCCATGCGCGTGCGTGTCTGTACAGTGGCACTGgcg  >  1:467440/1‑150 (MQ=255)
                                                                                                                                                      |                                                                                                                                                     
GCCGGGCAGAGGTTGCGGGTGTGGCCGGTGCGGTTGAAGGGGGAGGCGGTGGAGCTGGCGGTGGCCTGATACCCCTTTGCTGCCTGTGCCGGCCTCTTCGCGGGCAAGCCCGCTCCTACGGAATCTGCACCGCCTTGTAATACGGCGCCGACCCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCAGCACAAGCCGACTAGATAACCAACCCCTGCGGCGTGCGCCCGGCAAACACATCCACCAGGCTGGCCACCACCATCTCGCCCATGCGCGTGCGTGTCTGTACAGTGGCACTGGCG  >  NC_002947/1911583‑1911881

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