Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,911,733 A→C 5.6% intergenic (+82/+46) 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,7330AC 5.6% 533.9 / 4.6 177intergenic (+82/+46)nirD/PP_1707nitrite reductase/D‑isomer specific 2‑hydroxyacid dehydrogenase family protein
Reads supporting (aligned to +/- strand):  ref base A (79/83);  new base C (6/8);  total (86/91)
Fisher's exact test for biased strand distribution p-value = 7.84e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 177 total aligned reads displayed.

GCAGAGGTTGCGGGTGTGGCCGGTGCGGTTGAAGGGGGAGGCGGTGGAGCTGGCGGTGGCCTGATACCCCTTTGCTGCCTGTGCCGGCCTCTTCGCGGGCAAGCCCGCTCCTACGGAATCTGCACCGCCTTGTAATACGGCGCCGACCCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCAGCACAAGCCGACTAGATAACCAACCCCTGCGGCGTGCGCCCGGCAAACACATCCACCAGGCTGGCCACCACCATCTCGCCCATGCGCGTGCGTGTCTGTACAGTGGCACTGGCG  >  NC_002947/1911588‑1911881
                                                                                                                                                 |                                                                                                                                                    
  tgggcgagatggtggtggccagcctggtggatgtgtttgccgggcgcacgccgcaggggttggttatctagtcggcttgtgctGGCCTCTTCGCGGGCAAGCCCGCTCCTACGGAATCTGCACCGCCTTGTAATACGGCGCCGACCCtgt                                                                                                                                                <  1:4282665‑M1/67‑1 (MQ=255)
  agagGTTGCGGGTGTGGCCGGTGCGGTTGAAGGGGGAGGCGGTGGAGCTGGCGGTGGCCTGATACCCCTTTGCTGCCTGTGCCGGCCTCTTCGCGGGCAAGCCCGCTCCAACCGAAGCTGCACCGCCGTGTAATACGGCGCCCACCCggt                                                                                                                                                >  2:2899426/1‑150 (MQ=39)
   gagGTTGCGGGTGTGGCCGGTGCGGTTGAAGGGGGAGGCGGTGGAGCTGGCGGTGGCCTGATACCCCTTTGCTGCCTGTGCCGGCCTCTTCGCGGGCAAGCCCGCTCCTACGGAATCTGCACCGCCTTGTAATACGGCGCCGACCCtgtg                                                                                                                                               <  1:3954326/150‑1 (MQ=255)
      gTTGCGGGTGTGGCCGGTGCGGTTGAAGGGGGAGGCGGTGGAGCTGGCGGTGGCCTGATACCCCTTTGCTGCCTGTGCCGGCCTCTTCGCGGGCAAGCCCGCTCCTACGGAATCTGCACCGCCTTGTAATACGGCGCCGACCCTGTGGGa                                                                                                                                            <  1:2139108/150‑1 (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                                                                                                                                        <  2:1649943/150‑1 (MQ=255)
          cGGGTGTGGCCGGTGCGGTTGAAGGGGGAGGCGGTGGAGCTGGCGGTGGCCTGATACCCCTTTGCGGCCTGTGCCGGCCTCTTCGCGGGCAAGCCCGCTCCTACGGAATCTGCACCGCCTTGTAATACGGCGCCGACCCTGTGGGAGCgg                                                                                                                                        <  2:4899253/150‑1 (MQ=255)
          atggtggtggccagcctggtggatgtgtttgccgggcgcacgccgcaggggttggttatctagtcggcttgtgctGGCCTCTTCGCGGGCAAGCCCGCTCCTACGGAATCTGCACCGCCTTGTAATACGGCGCCGACCCTGTGGGAGCgg                                                                                                                                        >  2:2404394‑M1/76‑150 (MQ=255)
          atggtggtggccagcctggtggatgtgtttgccgggcgcacgccgcaggggttggttatctagtcggcttgtgctGGCCTCTTCGCGGGCAAGCCCGCTCCTACGGAATCTGCACCGCCTTGTAATACGGCGCCGACCCTGTGGGAGCgg                                                                                                                                        >  2:1806154‑M1/76‑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)
              tgtgGCCGGTGCGGTTGAAGGGGGAGGCGGTGGAGCTGGCGGTGGCCTGATACCCCTTTGCTGCCTGTGCCGGCCACATCGCCGGCCAGCCCCCGCCAAAAGAAAAAGCCCCCCCAGGGAATAAGGCGCCCCCCCCGGGGGAGCGGGcgt                                                                                                                                    >  1:800247/1‑148 (MQ=255)
               gtgGCCGGTGCGGTTGAAGGGGGAGGCGGTGGAGCTGGCGGTGGCCTGATACCCCTTTGCTGCCTGTGCCGGCCTCTTCGCGGGCAAGCCCGCTCCTACGGAATCTGCACCGCCTTGTAATACGGCGCCGACCCTGTGGGAGCGGGCTTg                                                                                                                                   <  2:698164/150‑1 (MQ=255)
               ggtggccagcctggtggatgtgtttgccgggcgcacgccgcaggggttggttatctagtcggcctgtgccGGCCTCTTCGCGGGCAAGCCCGCTCCTACGGAATCTGCACCGCCTTGTAATACGGCGCCGACCCTGTGGGAGCGGGCTTg                                                                                                                                   >  2:4007004‑M1/71‑150 (MQ=255)
                tgGCCGGTGCGGTTGAAGGGGGAGGCGGTGGAGCTGGCGGTGGCCTGATACCCCTTTGCTGCCTGTGCCGGCCTCTTCGCGGGCAAGCCCGCTCCTACGGAATCTGCACCGCCGTGTAATACGGCGCCGACCCTGTGGGAGCGGGCTTGc                                                                                                                                  >  2:4146506/1‑150 (MQ=255)
                  gCCGGTGCGGTTGAAGGGGGAGGCGGTGGAGCTGGCGGTGGCCTGATACCCCTTTGCTGCCTGTGCCGGCCTCTTCGCGGGCCAGCCCCCACCTACGGAATAAGCACCGCCGTGTAATAAGGCGCCGACCCGGTGGGAGCGGGCTTGccc                                                                                                                                >  1:4143782/1‑150 (MQ=255)
                  gCCGGTGCGGTTGAAGGGGGAGGCGGTGGAGCTGGCGGTGGCCTGATACCCCTTTGCTGCCTGTGCCGGCCTCTTCGCGGGCAAGCCCGCTCCTACGGAATCTGCCCCGCCGTGTAATACGGCGCCGACCCGGTGGGAGCGGGCTTGccc                                                                                                                                >  2:3624876/1‑150 (MQ=255)
                  gCCGGTGCGGTTGAAGGGGGAGGCGGTGGAGCTGGCGGTGGCCTGATACCCCTTTGCTGCCTGTGCCGGCCTCTTCGCGGGCAAGCCCGCTCCTACGGAATCTGCACCGCCTTGTAATACGGCGCCGACCCTGTGGGAGCGGGCTTGccc                                                                                                                                <  1:4424102/150‑1 (MQ=255)
                     ggTGCGGTTGAAGGGGGAGGCGGTGGAGCTGGCGGTGGCCTGATACCCCTTTGCTGCCTGTGCCGGCCTCTTCGCGGGCAAGCCCGCACCAACCGAAAAAGCACCGCCGTGTAATACGGCCCCCCCCCAGTGGGCGCGGGCGTGcccccc                                                                                                                             >  1:5076414/1‑147 (MQ=255)
                       tGCGGTTGTTGGGGGAGGGGGTGGAGGTGGCGGTGGCCTGATACCCCTTTGCTGCCTGTGCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACGGAATCTGCACCGCCTTGTAATACGGCGCCGACCCTGTGGGAGCGGGCTTGCCCGCGaa                                                                                                                           <  1:1276197/150‑1 (MQ=39)
                       tGCGGTTGAAGGGGGCGGCGGTGGAGCTGGCGGTGGCCTGATACCCCTTTGCTGCCTGTGCCGGCCTCTTCGCGGGCAAGCCCGCTCCTACGGAATCTGCACCGCCTTGTAATACGGCGCCGACCCTGTGGGAGCGGGCTTGCCCGCGaa                                                                                                                           <  2:3436282/150‑1 (MQ=255)
                       tGCGGTTGAAGGGGGAGGCGGTGTAGCTGGCGGTGGCCTGATACCCCTTTGCTGCCTGTGCCGGCCTCTTCGCGGGCAAGCCCGCTCCTACGGAATCTGCACCGCCTTGTAATACGGCGCCGACCCTGTGGGAGCGGGCTTGCCCGCGaa                                                                                                                           <  2:771494/150‑1 (MQ=255)
                       tGCGGTTGAAGGGGGAGGCGGTGGAGCTGGCGGTGGCCTGATACCCCTTTGCTGCCTGTGCCGGCCTCTTCGCGGGCAAGCCCGCTCCTACGGAATCTGCACCGCCTTGTAATACGGCGCCGACCCTGTGGGAGCGGGCTTGCCCGCGaa                                                                                                                           <  1:2806065/150‑1 (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:2695985/150‑1 (MQ=255)
                              gAAGGGGGAGGCGGTGGAGCTGGCGGTGGCCTGATACCCCTTTGCTGCCTGTGCCGGCCTCTTCGCGGGCAAGCCCGCTCCTACGGAATCTGCACCGCCTTGTATTACGGCGCCGACCCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCa                                                                                                                    <  2:2588200/150‑1 (MQ=255)
                                    gtttgccgggcgcacgccgcaggggttggttatctagtcggcttgtgccGGCCTCTTCGCGGGCAAGCCCGCTCCTACGGAATCTGCACCGCCTTGTAATACGGCGCCGACCCTGTGGGAGCGGGc                                                                                                                                      <  2:2599615‑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)
                                          ggTGGAGCTGGCGGTGGCCTGATACCCCTTTGCTGCCTGTGCCGGCCTCTTCGCGGGCAAGCCCGCTCCTACGGAATCTGCACCGCCTTGTAATACGGCGCCGACCCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCAGCACAAGCCGAc                                                                                                        <  2:4032106/150‑1 (MQ=255)
                                            tGGGGCTGGCGGTGGCCTGATACCCCTTTGCTGCCTGTGCCGGCCTCTTCGCGGGCAAGCCCCCGCCTACGGCATATGCACCGCCGTGTAAAAACGCGCCCACCCCGTGGGAGCGGGCGTGCCCCCGAAGAGGCCCGCACAAGCCCAcaa                                                                                                      >  2:5175034/1‑148 (MQ=255)
                                            ggcgcacgccgcaggggttggttatctagtcggcctgtgccGGCCTCTTCGCGGGCAAGCCCGCTCCTACGGAATCTGCACCGCCTTGTAATACGGCGCCGACCCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCAGCACAAGCCGACTa                                                                                                      >  2:4604013‑M1/42‑150 (MQ=255)
                                                           ccTGATACCCCTTTGCTGCCTGTGCCGGCCTCTTCGCGGGCAAGCCCGCGCCGACGGAATCTGCACCGCCGTGGAATACGGCGCCGACCCAGGGGGAGCGGGCTTGCCCGCGAAGAGGCCCCCACAAGCCGACAAGATAACCAAccccag                                                                                       >  2:3597007/1‑148 (MQ=38)
                                                                   cccTTTGCTGCCTGTGCCGGCCTCTTCGCGGGCAAGCCCCCTCCTACGGAATCAGCACCGCCGTGTAATACGGCGCCGACCCTGTGGGAGCGGGCTTGCCCGCCAAGAGGCCAGCACAAGCCGACAAGATAAACAACCCCAGCGGCGTgc                                                                               >  2:5084828/1‑150 (MQ=38)
                                                                      tttGCTGCCTGTGCCGGCCTCGTCGCGGGCCAGCCCCCGCCGACCGCAACTGCCCCCCCCGGGAAAAACGCGCCCCCCCCGGGGGAGCCGGCGTGCCCCCGCACAGGCCCCCCCAACCCCCCAAGAGAACCCACCCCCGCGGCGGGCGcc                                                                            >  2:1053250/1‑150 (MQ=255)
                                                                       gggctgcCTGTGCCGGCCTCTTCGCGGGCTAGCCCGCTCCTCCGGAATGTGCACCGCCTTGTAATACGGCGCCGACCCGGTGGGAGCGGGCTTGCCCGCGAAGAGGCCAGCACAAGCCGACTAGATAACCAACCCCTGCGGCGTGCGccc                                                                           <  2:4207741/148‑1 (MQ=255)
                                                                        tgctgcCTGTGCCGTCATCTTCGCGGTCAAGGCGGATCCTACTGATTCTGCACCGCCTTGTAATACGGCGGCGACCCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCAGCACAAGCCGACTAGATAACCAACCCCTGCGGCGTGCGCCCg                                                                          <  2:4403035/150‑1 (MQ=38)
                                                                        tgctgcCTGTGCCGGCCTCTTCGCGGGCAAGCCCGCTCCTACGGAATCTGCACCGCCTTGTAATACGGCGCCGACCCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCAGCACAAGCCGACTAGATAACCAACCCCTGCGGCGTGCGCCCg                                                                          <  2:4403169/150‑1 (MQ=255)
                                                                        tgctgcCTGTGCCGGCCGCTGGGCGGCCAAGCCCGCTCCTACGGAATCTGCACCGCCTTGTAATACGGCGCCGACCCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCAGCACAAGCCGACTAGATAACCAACCCCTGCGGCGTGCGCCCg                                                                          <  2:5415649/150‑1 (MQ=255)
                                                                             ccTGGGCCGGCCTCTTCGGGGGCAAGCCCGCACCTCCGGAATCTGCACCGCCTTGTAATCCGGCGCCGACCCTGTGGGCGCGGGCTGGCCCGCGAAGAGGCCAGCACAAGCCGACTAGATAACCAACCCCTGCGGCGTGCGCCCGGCaaa                                                                     <  1:1761245/150‑1 (MQ=255)
                                                                              cTGTGCCGGCCTCTTCGCGGGCAAGCCCGCTCCTACGGAATCTGCACCGCCTTGTAATACGGCGCCGACCCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCAGCACAAGCCGACTAGATAACCAACCCCTGCGGCGTGCGCCCGGCAAac                                                                    <  2:829401/150‑1 (MQ=255)
                                                                                    gGGCCTCTTCGCGGGCAAGCCCGCTCCTACGGTATATGCACCGCCTTGTATTAGGGCGCCGACCGTGGGGGTGCGGGCTGGCCCGCGAAGAGGCCAGCACAAGCCGACTAGATAACCAACCCCTGCGGCGTGCGCCCGGCAATCACATcc                                                              >  2:1276384‑M1/2‑150 (MQ=255)
                                                                                     GGCCTCTTCGCGGGCAAGCCCGCTCCTACGGAATCTGCACCGCCTTGTAATACGGCGCCGACCCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCAGCACAAGCCGACTAGATAACCAACCCCTTCGGCGTGCGCCCGGCAAACACATcca                                                             <  2:3251115‑M1/150‑1 (MQ=255)
                                                                                        CTCTTCGCGGGCAAGCCCGCTCCTACGGAATCTGCACCGCCTTGTAATACGGCGCCGACCCTGTGGGAGCGGGCTTGCCCGCGCAGAGGCCAGCACAAGCCGACTAGATAACCAACCCCTGCGGCGTGCGCCCGGCAAACACATccacca                                                          <  2:5501805‑M1/150‑1 (MQ=255)
                                                                                          CTTCGCGGGCACGCCCGCTCCTACGTAATCTGCACCGCCTTTTATTACGGCGCCGACCCTGTGGGAGGGGGCTTGCCCGCGAAGAGGCCAGCACAAGCCGACTAGATAACCAACCCCTGCGGCGTACGCCCGGCAAACACATCCACCAgg                                                        >  1:2762066‑M1/1‑150 (MQ=255)
                                                                                           TTCGCGGGCAAGCCCGCTCCTACGGAATCTGCACCGCCTTGTAATACGGCGCCGACCCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCAGCACAAGCCGACTAGATAACCAACCCCTGCGGCGTGCGCCCGGCAAACACATCCACCAGGc                                                       <  2:2797522‑M1/150‑1 (MQ=255)
                                                                                                ggccacccccGCCCCCCCGGAACGTGCACCGCTTTTAACTCCGGCGCCGCCCGGGGGGGGGCGGGCGCCCCCGCGAAGAGGCCAGCACAAGCCGACTAGATAACCAACCCCTGCGGCGTGCGCCCGGCAAACACATCCACCAGGCTGGcc                                                  <  1:2364123/143‑1 (MQ=255)
                                                                                                GGGCAAGCCCGCTCCTACGGAATCTGCACCGCCTTGTAATACGGCGCCGACCCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCAGCACAAGCCGACTAGATAACCAACCCCTGCGGCGTGCGCCCGGCAAACACATCCACCAGGCTGGcc                                                  <  2:1520411‑M1/150‑1 (MQ=255)
                                                                                                 GGCAAGCCCGCTCCTCCGGAATCTGCACCGCCTTGTAATCCGGCGCCGACCCTGTGGGCGCGGGCTTGCCCGCGAAGAGGCCAGCACAAGCCGACTAGATAACCAACCCCTGCGGCGTGCGCCCGGCAAACACATCCACCAGGCTGGcca                                                 >  2:3040426‑M1/1‑150 (MQ=255)
                                                                                                 GGCAAGCCCGCTCCTACGGAATCTGCACCGCCTTGTAATACGGCGCCGACCCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCAGCACAAGCCGACTAGATAACCAACCCCTGCGGCGTGCGCCCGGCAAACACATCCACCAGGCTGGcca                                                 <  2:3217852‑M1/150‑1 (MQ=255)
                                                                                                   CAAGCCCGCTCCTACGGAATCTGCACCGCCTTGTAATACGGCGCCGACCCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCAGCACAAGCCGACTAGATAACCAACCCCTGCGGCGTGCGCCCGGCAAACACATCCACCAGGCTGGccacc                                               <  1:2505320‑M1/150‑1 (MQ=255)
                                                                                                   CAAGCCCGCTCCTACGGAATCTGCACCGCCTTGTAATACGGCGCCGACCCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCAGCACAAGCCGACTAGATAACCAACCCCTGCGGCGTGCGCCCGGCAAACACATCCACCAGGCTGGccacc                                               <  1:2505064‑M1/150‑1 (MQ=255)
                                                                                                    AAGCCCGCTCCTACGGAATCTGCACCGCCTTGTAATACGGCGCCGACCCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCAGCACAAGCCGACTAGATAACCAACCCCTGCGGCGTGCGCCCGGCAAACACATCCACCAGGCTGGccacca                                              <  2:5433155‑M1/150‑1 (MQ=255)
                                                                                                       CCCGCTCCTACGGAAACTGCACCGCCTTGAAAACCGGCGCCGACCCCGTGGGAGCGGGCTGGCCCGCGAAGAGGCCAGCACAAGCCGACTAGATAACCAACCCCTGCGGCGTGCGCCCGGCAAACACATCCACCAGGCTGGccaccacca                                           >  1:1144861‑M1/1‑150 (MQ=255)
                                                                                                        CCGCTCCTACGGAATCTGCACCGCCTTGTAATACGGCGCCGACCCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCAGCACAAGCCGACTAGATAACCAACCCCTGCGGCGTGCGCCCGGCAAACACATCCACCAGGCTGGCCACCACCAt                                          <  2:2382646‑M1/150‑1 (MQ=255)
                                                                                                        CCGCTCCTACGGAATCTGCACCGCCTTGTAATACGGCGCCGACCCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCAGCACAAGCCGACTAGATAACCAACCCCTGCGGCGTGCGCCCGGCAAACACATCCACCAGGCTGGCCACCACCAt                                          <  2:2382317‑M1/150‑1 (MQ=255)
                                                                                                         CGCTCCTACGGAATCTGCACCGCCTTGTAATACGGCGCCGACCCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCAGCACAAGCCGACTAGATAACCAACCCCTGCGGCGTGCGCCCGGCAAACACATCCACCAGGCTGGCCACCACCAtc                                         >  1:3217883‑M1/1‑150 (MQ=255)
                                                                                                            TCCTACGGAATCTGCACCGCCTTGTAATACGGCGCCGACCCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCAGCACAAGCCGACTAGATAACCAACCCCTGCGGCGTGCGCCCGGCAAACACATCCACCAGGCTGGCCACCACCATCTCg                                      <  1:4604572‑M1/150‑1 (MQ=255)
                                                                                                            TCCTACGGAATCTGCACCGCCTTGTAATACGGCGCCGACCCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCAGCACAAGCCGACTAGATAACCAACCCCTGCGGCGTGCGCCCGGCAAACACATCCACCAGGCTGGCCACCACCATCTCg                                      <  1:3257540‑M1/150‑1 (MQ=255)
                                                                                                              CTACGGAATCTGCACCGCCTTGTCATACGGCGCCGACCCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCAGCACAAGCCGACTAGATAACCAACCCCTGCGGCGTGCGCCCGGCAAACACATCCACCAGGCTGGCCACCACCATCTCGcc                                    <  1:4507816‑M1/150‑1 (MQ=255)
                                                                                                              CTACGGAATCTGCACCGCCTTGTAATACGGCGCCGACCCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCAGCACAAGCCGACTAGATAACCAACCCCTGCGGCGTGCGCCCGGCAAACACATCCACCAGGCTGGCCACCACCATCTCGcc                                    <  1:1721888‑M1/150‑1 (MQ=255)
                                                                                                               tACGGAATCTGCACCGCCTTGTAATACGGCGCCGACCCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCAGCACAAGCCGACTAGATAACCAACCCCTGCGGCGTGCGCCCGGCAAACACATCCACCAGGCTGGCCACCACCATCTCGccc                                   >  2:1153307/1‑150 (MQ=255)
                                                                                                               tACGGAATCTGCACCGCCTTGTAATACGGCGCCGACCCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCAGCACAAGCCGACTAGATAACCAACCCCTGCGGCGTGCGCCCGGCAAACACATCCACCAGGCTGGCCACCACCATCTCGccc                                   >  2:2057779/1‑150 (MQ=255)
                                                                                                               tACGGAATCTGCACCGCCTTGTAATACGGCGCCGACCCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCAGCACAAGCCGACTAGATAACCAACCCCTGCGGCGTGCGCCCGGCAAACACATCCACCAGGCTGGCCACCACCATCTCGccc                                   >  2:3880704/1‑150 (MQ=255)
                                                                                                                aCGGAATCTGCACCGCCTTGTAATACGGCGCCGACCCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCAGCACAAGCCGACTAGATAACCAACCCCTGCGGCGTGCGCCCGGCAAACACATCCACCAGGCTGGCCACCACCATCTCGCCCa                                  >  2:2436780/1‑150 (MQ=255)
                                                                                                                 cGGAATCTGCACCGCCTTGTAATACGGCGCCGACCCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCAGCACAAGCCGACTAGATAACCAACCCCTGCGGCGTGCGCCCGGCAAACACATCCACCAGGCTGGCCACCACCATCTCGCCCAt                                 >  1:902757/1‑150 (MQ=255)
                                                                                                                 cGGAATCTGCACCGCCTTGTAATACGGCGCCGACCCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCAGCACAAGCCGACTAGATAACCAACCCCTGCGGCGTGCGCCCGGCAAACACATCCACCAGGCTGGCCACCACCATCTCGCCCAt                                 >  2:3768000/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)
                                                                                                                         gCACCGCCTTGTAATACGGCGCCGACCCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCAGCACAAGCCGACTAGATAACCAACCCCTgcg                                                                                     >  1:320255/1‑90 (MQ=255)
                                                                                                                           aCCGCCTTGTAATACGGCGCCGACCCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCAGCACAAGCCGACTAGATAACCAACCCCTGCGGCGTGCGCCCGGCAAACACATCCACCAGGCTGGCCACCACCATCTCGCCCATGCgcgtgcgt                       <  1:5599550/150‑1 (MQ=255)
                                                                                                                           aCCGCCTTGTAATACGGCGCCGACCCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCAGCACAAGCCGACTAGATAACCAACCCCTGCGGCGTGCGCCCGGCAAACACATCCACCAGGCTGGCCACCACCATCTCGCCCATGCgcgtgcgt                       <  1:5599754/150‑1 (MQ=255)
                                                                                                                             cGCCTTGTAATACGGCGCCGACCCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCAGCACAAGCCGACTAGATAACCAACCCCTGCGGCGTGCGCCCGGCAAACACATCCACCAGGCTGGCCACCACCATCTCGCCCATGCGCGTGCgtgt                     >  2:4688514/1‑150 (MQ=255)
                                                                                                                             cGCCTTGTAATACGGCGCCGACCCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCAGCACAAGCCGACTAGATAACCAACCCCTGCGGCGTGCGCCCGGCAAACACATCCACCAGGCTGGCCACCACCATCTCGCCCATGCGCGTGCgtgt                     >  2:2162170/1‑150 (MQ=255)
                                                                                                                             cGCCTTGTAATACGGCGCCGACCCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCAGCACAAGCCGACTAGATAACCAACCCCTGCGGCGTGCGCCCGGCAAACACATCCACCAGGCTGGCCACCACCATCTCGCCCATGCGCGTGCgtgt                     <  2:4507805/150‑1 (MQ=255)
                                                                                                                             cGCCTTGTAATACGGCGCCGACCCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCAGCACAAGCCGACTAGATAACCAACCCCTGCGGCGTGCGCCCGGCAAACACATCCACCAGGCTGGCCACCACCATCTCGCCCATGCGCGTGCgtgt                     >  2:3026659/1‑150 (MQ=255)
                                                                                                                             cGCCTTGTAATACGGCGCCGACCCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCAGCACAAGCCGACTAGATAACCAACCCCTGCGGCGTGCGCCCGGCAAACACATCCACCAGGCTGGCCACCACCATCTCGCCCATGCGCGTGCgtgt                     >  1:486250/1‑150 (MQ=255)
                                                                                                                               ccTTTTAATCGGGCGCCGCCCCTGGGGGCGGGGGCTGGCCCGCGAAGAGGCCAGCACAAGCCGACTAGATAACCAACCCCTGCGGCGTGCGCCCGGCAAACACATCCACCAGGCTGGCCACCACCATCTCGCCCATGCGCGTGCGTGTCt                   <  1:1635676/150‑1 (MQ=255)
                                                                                                                                 ttGTAATACGGCGCCGACCCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCAGCACAAGCCGACTAGATAACCAACCCCTGCGGCGTGCGCCCGGCAAACACATCCACCAGGCTGGCCACCACCATCTCGCCCATGCGCGTGCGTGTCTGt                 >  2:4961366/1‑150 (MQ=255)
                                                                                                                                 ttGTAATACGGCGCCGACCCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCAGCACAAGCCGACTAGATAACCAACCCCTGCGGCGTGCGCCCGGCAAACACATCCACCAGGCTGGCCACCACCATCTCGCCCATGCGCGTGCGTGTCTGt                 >  1:5121828/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)
                                                                                                                                         cGGCGCCGACCCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCAGCACAAGCCGACTAGATAACCAACCCCTGCGGCGTGCGCCCGGCAAACACATCCACCAGGCTGGCCACCACCATCTCGCCCATGCGCGTGCGTGTCTGTACAGTGGc         <  2:5126409/150‑1 (MQ=255)
                                                                                                                                           gcgcCGACCCTGTGGGTGGGGGGTTGCCCGCGAAGTGGCCAGCACAAGCCGACTAGATAACCAACCCCTGCGGCGTGCGCCCGGCAAACTCATCCACCAGGCTGGCCTCCTGGTTCTCGCCCATGCGCGTGCGTGTCTGTACAGTGGCAc       <  2:1717910/150‑1 (MQ=255)
                                                                                                                                               cGACCCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCAGCACAAGCCGACTAGATAACCAACCCCTGCGGCGTGCGCCCGGCAAACACATCCACCAGGCTGGCCACCACCATCTCGCCCATGCGCGTGCGTGTCTGTACAGTGGCACTGgc   >  2:3588150/1‑150 (MQ=255)
                                                                                                                                               cGACCCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCAGCACAAGCCGACTAGATAACCAACCCCTGCGGCGTGCGCCCGGCAAACACATCCACCAGGCTGGCCACCACCATCTCGCCCATGCGCGTGCGTGTCTGTACAGTGGCACTGgc   <  2:1600529/150‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                    
GCAGAGGTTGCGGGTGTGGCCGGTGCGGTTGAAGGGGGAGGCGGTGGAGCTGGCGGTGGCCTGATACCCCTTTGCTGCCTGTGCCGGCCTCTTCGCGGGCAAGCCCGCTCCTACGGAATCTGCACCGCCTTGTAATACGGCGCCGACCCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCAGCACAAGCCGACTAGATAACCAACCCCTGCGGCGTGCGCCCGGCAAACACATCCACCAGGCTGGCCACCACCATCTCGCCCATGCGCGTGCGTGTCTGTACAGTGGCACTGGCG  >  NC_002947/1911588‑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: