Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,461,854 T→C 5.7% intergenic (+113/‑68) gapB → / → PP_2150 glyceraldehyde‑3‑phosphate dehydrogenase/thiamine biosynthesis lipoprotein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,461,8540TC 5.7% 657.3 / 18.1 210intergenic (+113/‑68)gapB/PP_2150glyceraldehyde‑3‑phosphate dehydrogenase/thiamine biosynthesis lipoprotein
Reads supporting (aligned to +/- strand):  ref base T (88/108);  new base C (2/10);  total (92/118)
Fisher's exact test for biased strand distribution p-value = 7.25e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.16e-01

Only 100 of 220 total aligned reads displayed.

GGCTGGTGTGAACCCGCCAGCGTTTCCACGCTGATTCGCTTGTAAGGCAATGAAAAAACGGGAACCTTCGGGTTCCCGTTTTTTTTTATTCTGAATGTGTGTTGTCAGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGATTATCGCTGCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCTGGAGCGCCTGGGTGGTCCGA  >  NC_002947/2461708‑2462003
                                                                                                                                                    |                                                                                                                                                     
ggCTGGTGTGAACCCGCCAGCGTTTCCACGCTGATTCGCTTGTAAGGCAATGAAAAAACGGGAACCTTCGGGTTCCCGTTTTTTTTTATTCTGAATGTGTGTTGTCAGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGATTATc                                                                                                                                                    <  1:4516191/150‑1 (MQ=255)
ggCTGGTGTGAACCCGCCAGCGTTTCCACGCTGATTCGCTTGTAAGGCAATGAAAAAACGGGAACCTTCGGGTTCCCGTTTTTTTTTATTCTGAATGTGTGTTGTCAGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGATTATc                                                                                                                                                    <  2:375415/150‑1 (MQ=255)
ggCTGGTGTGAACCCGCCAGCGTTTCCACGCTGATTCGCTTGTAAGGCAATGAAAAAACGGGAACCTTCGGGTTCCCGTTTTTTTTTATTCTGAATGTGTGTTGTCAGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGATTATc                                                                                                                                                    <  1:8087510/150‑1 (MQ=255)
ggCTGGTGTGAACCCGCCAGCGTTTCCACGCTGATTCGCTTGTAAGGCAATGAAAAAACGGGAACCTTCGGGTTCCCGTTTTTTTTTATTCTGAATGTGTGTTGTCAGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGATTATc                                                                                                                                                    <  1:8343460/150‑1 (MQ=255)
   tGGTGTGAACCCGCCAGCGTTTCCACGCTGATTCGCTTGTAAGGCAATGAAAAAACGGGAACCTTCGGGTTCCCGTTTTTTTTTATTCTGAATGTGTGTTGTCAGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGATTATCGCt                                                                                                                                                 <  2:3813350/150‑1 (MQ=255)
    ggTGTGAACCCGCCAGCGTTTCCACGCTGATTCGCTTGTAAGGCAATGAAAAAACGGGAACCTTCGGGTTCCCGTTTTTTTTTATTCTGAATGTGTGTTGTCAGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGATTATCGCTg                                                                                                                                                <  2:2374700/150‑1 (MQ=255)
     gtgtGAACCCGCCAGCGTTTCCACGCTGATTCGCTTGTAAGGCAATGAAAAAACGGGAACCTTCGGGTTCCCGTTTTTTTTTATTCTGAATGTGTGTTGTCAGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGATTATCGCTGc                                                                                                                                               <  2:5171381/150‑1 (MQ=255)
       gtgAACCCGCCAGCGTTTCCACGCTGATTCGCTTGTAAGGCAATGAAAAAACGGGAACCTTCGGGTTCCCGTTTTTTTTTATTCTGAATGTGTGTTGTCAGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGATTATCGCTGccc                                                                                                                                             <  1:840892/150‑1 (MQ=255)
       gtgAACCCGCCAGCGTTTCCACGCTGATTCGCTTGTAAGGCAATGAAAAAACGGGAACCTTCGGGTTCCCGTTTTTTTTTATTCTGAATGTGTGTTGTCAGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGATTATCGCTGccc                                                                                                                                             <  2:862074/150‑1 (MQ=255)
       gtgAACCCGCCAGCGTTTCCACGCTGATTCGCTTGTAAGGCAATGAAAAAACGGGAACCTTCGGGTTCCCGTTTTTTTTTATTCTGAATGTGTGTTGTCAGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGATTATCGCTGccc                                                                                                                                             <  1:662453/150‑1 (MQ=255)
          aaCCCGCCAGCGTTTCCACGCTGATTCGCTTGTAAGGCAATGAAAAAACGGGAACCTTCGGGTTCCCGTTTTTTTTTATTCTGAATGTGTGTTGTCAGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCGACAGGATTATCGCTGccccgg                                                                                                                                          >  2:145600/1‑148 (MQ=255)
           aCCCGCCAGCGTTTCCACGCTGATTCGCTTGTAAGGCAATGAAAAAACGGGAACCTTCGGGTTCCCGTTTTTTTTTATTCTGAATGTGTGTTGTCAGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGATTATCGCTGCCCCTGt                                                                                                                                         <  1:1025613/150‑1 (MQ=255)
           aCCCGCCAGCGTTTCCACGCTGATTCGCTTGTAAGGCAATGAAAAAACGGGAACCTTCGGGTTCCCGTTTTTTTTTATTCTGAATGTGTGTTGTCAGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGATTATCGCTGCCCCTGt                                                                                                                                         <  1:30311/150‑1 (MQ=255)
           aCCCGCCAGCGTTTCCACGCTGATTCGCTTGTAAGGCAATGAAAAAACGGGAACCTTCGGGTTCCCGTTTTTTTTTATTCTGAATGTGTGTTGTCAGAGCCGGCCCTTTCGCGGGTTAACCCGCGCCAACACGCTTATAGCTGccccagg                                                                                                                                         >  1:455375/1‑147 (MQ=255)
            cccGCCAGCGTTTCCACGCTGATTCGCTTTTAAGGCAATGAAAAAACGGGAACCTTCGGGTTCACGTTTTTTTTTATTCTGAATGTGTGTTGTCAGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGATTATCGCTGCCCCTGTa                                                                                                                                        <  1:7234646/150‑1 (MQ=255)
            cccGCCAGCGTTTCCACGCTGATTCGCTTGTAAGGCAATGAAAAAACGGGAACCTTCGGGTTCCCGTTTTTTTTTATTCTGAATGTGTGTTGTCAGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGATTATCGCTGCCCCTGTa                                                                                                                                        <  1:6904843/150‑1 (MQ=255)
             ccGCCAGCGTTTCCACGCTGATTCGCTTGTAAGGCAATGAAAAAACGGGAACCTTCGGGTTCCCGTTTTTTTTTATTCTGAATGTGTGTTGTCAGTGCCGGCCCTTTCGCGGGTTAAACCCCTACAACACGGTTATAGCATCCCCAGTAg                                                                                                                                       >  1:3901228/1‑150 (MQ=255)
             ccGCCAGCGTTTCCACGCTGATTCGCTTGTAAGGCAATGAAAAAACGGGAACCTTCGGGTTCCCGTTTTTTTTTATTCTGAATGTGTGTTGTCAGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGATTATCGCTGCCCCTGTAg                                                                                                                                       <  1:7388541/150‑1 (MQ=255)
             ccGCCAGCGTGTCCACGCTGATTCGCTTGTAAGGCCATGAAAAAACGGGCACCGTCGGGTTCCCGTTTTTTTTTATTCTGAATGTGTGTTGTCAGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGATTATCGCTGCCCCTGTAg                                                                                                                                       <  1:678632/150‑1 (MQ=255)
                ccAGCGTTTCCACGCTGATTCGCTTGTAAGGCAATGAAAAAACGGGAACCTTCGGGTTCCCGTTTTTTTTTATTCTGAATGTGTGTTGTCAGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGATTATCGCTGCCCCGGTAGGAg                                                                                                                                    >  2:3431026/1‑150 (MQ=255)
                ccAGCGTTTCCACGCTGATTCGCTTGTAAGGCAATGAAAAAACGGGAACCTTCGGGTTCCCGTTTTTTTTTATTCTGAATGTGTGTTGTCAGTGCCGGCCCTTTCGCGGGGGAAACCGCGCCGACAGGAGTATCGCTGCCCCCGTAgggg                                                                                                                                    >  1:224868/1‑148 (MQ=255)
                  aGCGTTTCCACGCTTATTCGCTTTTAAGGCAATGAAAAAACGGGAACCTTCGGGTTCCCGTTTTTTTGTTATTCTGAATGTGTGTTGTCAGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGATTATCGCTGCCCCTGTAGGAGc                                                                                                                                   <  2:2733915/150‑1 (MQ=255)
                   gCGTTTCCACGCTGATTCGCTTGTAAGGCAATGAAAAAACGGGAACCTTCGGGTTCCCGTTTTTTTTTATTCTGAATGTGTGTTGTCAGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGATTATCGCTGCCCCTGTAGGAGCgg                                                                                                                                 <  1:2592521/150‑1 (MQ=255)
                   gCGTTTCCACGCTGATTCGCTTGTAAGGCAATGAAAAAACGGGAACCTTCGGGTTCCCGTTTTTTTTTATTCTGAATGTGTGTTGTCAGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGATTATCGCTGCCCCTGTAGGAGCgg                                                                                                                                 <  1:4950298/150‑1 (MQ=255)
                   gCGTTTCCACGCTGATTCGCTTGTAAGGCAATGAAAAAACGGGAACCTTCGGGTTCCCGTTTTTTTTTATTCTGAATGTGTGTTGTCAGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGATTATCGCTGCCCCTGTAGGAGCgg                                                                                                                                 <  2:4348116/150‑1 (MQ=255)
                   gCGTTTCCACGCTGATTCGCTTGTAAGGCAATGAAAAAACGGGAACCTTCGGGTTCCCGTTTTTTTTTATTCTGAATGTGTGTTGTCAGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGATTATCGCTGCCCCGGTAGGAGCgg                                                                                                                                 >  2:2640390/1‑150 (MQ=255)
                    cGTTTCCACGCTGATTCGCTTGTAAGGCAATGAAAAAACGGGAACCTTCGGGTTCCCGTTTTTTTTTTATTCTGAATGTGTGTTGTCAGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGGGCAGCGATAATCCTGTAGGAGCgg                                                                                                                                 <  1:6328386/150‑1 (MQ=25)
                    cGTTTCCACGCTGATTCGCTTGTAAGGCAATGAAAAAACGGGAACCTTCGGGTTCCCGTTTTTTTTTATTCTGAATGTGTGTTGTCAGTGCCGGCCCTTTCGCGGGTTAAACCGCTCCGACAGGATTATCGCAGCCCCAGTTGGAGCggg                                                                                                                                >  2:751846/1‑150 (MQ=39)
                     gTTTCCACGCTGATTCGCTTGTAAGGCAATGAAAAAACGGGAACCTTCGGGTTCCCGTTTTTTTTTATTCTGAATGTGTGTTGTCAGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGATTATCGCTGCCCCTGTAGGAGCGGGt                                                                                                                               <  1:2121706/150‑1 (MQ=255)
                         ccACGCTGATTCGCTTGTAAGGCAATGAAAAAACGGGAACCTTCGGGTTCCCGTTTTTTTTTATTCTGAATGTGTGTTGTCAGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGATTATCGCTGCCCCTGTAGGAGCGGGTTTAc                                                                                                                           <  2:847673/150‑1 (MQ=255)
                         ccACGCTGAGTCGCTTGTAAGGCACTGCAACAACGGGAACCGTCGGGTTCCCGTTTTTTTTTATTCTGAAGGTGTGTTGTCAGGGCCGGCCCTTTCGCGGGTAAACCCGCTCCTCCAGGATTATCGCTGCCCCTGTAGGAGCGGGTTTAc                                                                                                                           <  2:7455426/150‑1 (MQ=38)
                          cACGCTGATTCGCTTGTAAGGCAATGAAAAAACGGGAACCTTCGGGTTCCCGTTTTTTTTTATTCTGAATGTGTGTTGTCAGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGCTGATCGCAGCCCCGGTAGGAGCGGGTTTAcc                                                                                                                          >  2:7013164/1‑150 (MQ=255)
                          cACGCTGATTCGCTTGTAAGGCAATGAAAAAACGGGAACCTTCGGGTTCCCGTTTTTTTTTATTCTGAATGTGTGTTGTCAGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGATTATCGCTGCCCCTGTAGGAGCGGGTTTAcc                                                                                                                          <  1:588941/150‑1 (MQ=255)
                          cACGCTGATTCGCTTGTAAGGCAATGAAAAAACGGGAACCTTCGGGTTCCCGTTTTTTTTATTCTGAATGTGTGTTGTCAGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGATTATAGCTGCCCCTGTAGGAGCGGGTTTAccc                                                                                                                         >  2:7542761/1‑150 (MQ=255)
                             gctgctGCGCTTGTAAGGCAATGAAAAAACGGGAACCTTCGGGTTCCCGTTTTTTTTTATTCTGAATGTGTGTTGTCAGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGATTATCGCTGCCCCTGTAGGAGCGGGTTTACCcgc                                                                                                                       <  1:6477947/150‑1 (MQ=255)
                                gATTCGCTTGTACGGCAATGAAAAAACGGGAACCTTCGGGTTCCCGTTTTTTTTTATTCTGAATGTGTGTTGTCAGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGATTATCGCTGCCCCTGTAGGAGCGGGTTTACCCGCGaa                                                                                                                    <  2:1991929/150‑1 (MQ=255)
                                 aTTCGCTTGTACGGCAATGAAAAAACGGGAACCTTCGGGTTCCCGTTTTTTTTTATTCTGAATGTGTGTTGTCAGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGGGCAGCGATAATCCTGTAGGAGCGGGTTTACCCGCGaaa                                                                                                                   <  2:1611596/150‑2 (MQ=21)
                                  ttCGCTTGTAAGGCAATGAAAAAACGGGACCCTTCGGGTTCCCTTTTTTTTTTATTCTGAATGTGTGTTGTCAGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGATTATCGCTGCCCCGGTAGGAGCGGGTTTACCCGCGAaga                                                                                                                  <  2:5010016/150‑1 (MQ=255)
                                  tgcGCTTGTAAGGCATTGAAAAAACGGGAACCTTCGGGTACCCGTTTTTTTTTATTCTGACTGTGTGTTGTCAGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGATTATCGCTGCCCCTGTTGGTGCGGGTTTACCCGCGAaga                                                                                                                  <  2:6576902/148‑1 (MQ=39)
                                     gCTTGTAAGGCAATGAAAAAACGGGAACCTTCGGGTTCCCGTTTTTTTTTATTCTGAATGTGTGTTGTCAGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGATTATCGCAGCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGc                                                                                                               >  2:5234524/1‑150 (MQ=255)
                                        tGTAAGGCAATGAAAAAACGGGAACCTTCGGGTTCCCGTTTTTTTTTTATTCTGAATGTGTGTTGTCAGTGCCGGCCCTTTCGCGGGTGAAACCCCTACTAAAGGATTATCGCTGGCCCGGGTGGAGCGGGGTTAACCGCGAAGTgggcc                                                                                                             >  1:7721360/1‑147 (MQ=37)
                                          tAAGGCAATGAAAAAACGGGAACCTTCGGGTTCCCGTTTTTTTTTATTCTGAATGTGTGTTGTCAGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGATTATCGCTGCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCa                                                                                                          >  2:8542981/1‑150 (MQ=255)
                                          tAAGGCAATGAAAAAACGGGAACCTTCGGGTTCCCGTTTTTTTTTATTCTGAATGTGTGTTGTCAGTGCCGGCCCTTTCGCGGGTAAAACCGCTCCGACAGGGTTATCGCTGCCCCTGGAGGAGCGGGTTTACCCGCGCACGGGCCAGCa                                                                                                          >  2:5526407/1‑150 (MQ=38)
                                           acggCAATGAAAAAACGGCCCCCTCCGGGTACCCGATTTTTTTTATTCTGAATGTGTGTTGTCAGTGCCGGCCCTTGCGCGGGTAAACCCGCCCCGACAGGATTATCGCTGCCCCTGCAGGAGCGGGTTTACCCGCGAAGAGGCCAGCAc                                                                                                         <  1:6421720/148‑1 (MQ=38)
                                                aaTGAAAAAACGGGAACCTTCGGGTTCCCGTTTTTTTTTTATTCTGAATGTGTGTTGTCAGTGCCGGCCCTTTCGCGGGTAAAACCGCTCCTACAGGATTATCGCTGCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACCGGc                                                                                                     >  2:680597/1‑150 (MQ=255)
                                                 aaaaaaaaaCGGGAACCCGAAGGTTCCCGTTTTTTTTTATTCTGAATGTGTGTTGTCAGTGCCGGCCCTTTCGCGGGTTAAACCGCGCCGACAGGGGTATAGCTGCCCCGGGAGGGGGGGGGTTTCCCCCGGAGCGGCCCGCCCCgggca                                                                                                   >  1:1116782/4‑147 (MQ=255)
                                                   aaaaaaaCGGGAACCTTCGGGTTCCCGTTTTTTTTTATTCTGAATGTGTGTTGTCAGTGCCGGCCCTTTCGCGGGTAAAACCGCTCCTAAAGGCTTATCGCTGCCCCGGTAGGAGCGGGTGTACCCGCGAAGAGGCCAGCACCGGCAAcc                                                                                                 >  1:8014381/2‑150 (MQ=39)
                                                    aaaaaaCGGGAACCTTCGGGTTCCCGTTTTTTTTTATTCTGAATGTGTGTTGTCAGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGATTATCGCTGCCCCTGTAGGAGCgg                                                                                                                                 >  1:2029374/1‑117 (MQ=255)
                                                     aaaaaCGGGAACCTTCGGGTTCCCGTTTTTTTTTATTCTGAATGTGTGTTGTCAGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGATTATCGCTGCCCCGGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACCGGCAACCga                                                                                               >  1:2895297/1‑150 (MQ=255)
                                                      aaaaCGGGAACCTTCGGGTTCCCGTTTTTTTTTATTCTGAATGTGTGTTGTCAGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGATTATCGCTGCCCCt                                                                                                                                           <  2:6941502/105‑1 (MQ=255)
                                                       aaaCGGGAACCTTCGGGTTCCCGTTTTTTTTTATTCTGAATGTGTGTTGTCAGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACCGGATTATCGCTGCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACCGGCAACCgaga                                                                                             <  2:654074/150‑1 (MQ=255)
                                                        aaCGGGAACCCGAAGGTTCCCGTTTTTTTTTATTCTGAATGTGTGTTGTCAGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGATTATCGCTGCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGaa                                                                                            >  2:4247738/1‑150 (MQ=255)
                                                         aCGGGAACCCGAAGGTTCCCGTTTTTTTTTATTCTGAATGTGTGTTGTCAGTGCCGGCCCTTTCGCGGGTAAAACCGCTCCTACAGGATTATCGCTGCCCCGGTAGGAGAGGGTTTAACCCCGAAGAGGCCAGCACCGGCAAccccgaac                                                                                           >  2:2929919/1‑144 (MQ=38)
                                                           gggACCCTTCGGGTTCCCGTTTTTTTTTATTCTGAATGTGTGTTGTCAGTGCCGGCCCTTTCGCGGGTAAACCCGCGACGACAGGATTATCGCTGCCCCGGTCGGCGCGGGGTTACCCGCGAAGAGGCCCGCACCGGCCACCGCGAAccc                                                                                         >  2:184213/1‑150 (MQ=37)
                                                           gggAACCTTCGGGTTCCCGTTTTTTTTTATTCTGAATGTGTGTTGTCAGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGATTATCGCTGCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAAcc                                                                                         >  2:2112629/1‑150 (MQ=255)
                                                              acccTTCGGGTTCCCGATTTTTTTTATTCTGACTGTGTGTGGTCAGGGCCGGCCCTTCCGCGGGAAACCCCGCCCCTACAGGCTTATCGCTGCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTa                                                                                      <  1:4383771/148‑1 (MQ=38)
                                                                       gTTCCCGTTTTTTTTTATTCTGAATGTGTGTTGTCAGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCAACACGAATATCCCAGCCCCCGGACGACCCGGGTTACCCCCCAAGAGGCCCCCACCCGCAACCCAGCAACCCAACAGGgcgcc                                                                             >  2:4459054/1‑149 (MQ=19)
                                                                       gTTCCCGTTTTTTTTTATTCTGAATGTGTGTTGTCAGTGCCGGCCCTTTCGCGGGTAAAACCGCTCCAACAGGATTATAGCAGCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCAACTGTgcgcc                                                                             >  2:4697553/1‑149 (MQ=38)
                                                                               tttttttttATTCTGATTGTGTGTTGTCAGTGCCGGCCGTTTCGGGGGTAAACCGGCTCCTACTGTTTTATCGCTGCCCGTGTTGGAGCGGTTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCt                                                                      <  1:8594631/150‑1 (MQ=255)
                                                                                ttttttttATTCTTATTGTGTGTTTTCTGGGCGGGCCCTTGCGGGGTTAACCCCGCTCCTCCAGGTTTTGCGCTGCCCCTTTAGGTGGGGGTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTa                                                                     <  1:1618216/150‑1 (MQ=255)
                                                                                ttttttttATTCTGATTGTGTGTTGTCAGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGATTATCGCTGCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTa                                                                     <  2:359238/150‑1 (MQ=255)
                                                                                  ttttttttTCTGAATGTGTGTTGTCTGTGCCGGCCCTTTCGCGGTTAAACCCGCTCATACCGGATTATCGCTGCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACt                                                                   <  1:8769180/150‑1 (MQ=39)
                                                                                   tttttATTCTGAATGTGTGTTGTCAGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGATTATCGCTGCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTg                                                                  >  2:7201861/1‑150 (MQ=255)
                                                                                              gAATGTGTGTTGTCAGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGATTATCGCTGCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTctgc                                                        <  2:3228077/150‑1 (MQ=255)
                                                                                                 tgtgtgTTGTCTGTGCCGGCCCTTTCGCGGGTAAACCCGCTCATCCTGTTTTATCGCTGCCCCTGTTGGAGCGGTTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCtgctgc                                                     <  1:3395158/150‑1 (MQ=38)
                                                                                                   tgtgtTGTCAGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGATTATCGCTGCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTa                                                   <  1:5861383/150‑1 (MQ=255)
                                                                                                      gttgtCAGTGCCGGCCCTTTCGCGGGTAAAACCCCTCCGACAGGATTATAGCTGCCCCGGGAGGAGCGGGTTTTCCCGCGAAGAGGCCAGCCCCCGCAACCGAGAAACCCGACAGTGCGCACAGCCCGACTGTTCATTCTGCGGCTACTc                                                >  2:3131677/1‑150 (MQ=255)
                                                                                                           ctgTGCCGGCCCTTTCGGGGGTAAACCCGCTCCTACAGGTTTTTCGCTGCCCGTGTTGGTGGGGTTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGc                                           <  1:7035238/148‑1 (MQ=38)
                                                                                                             gTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGATTATCGCTGCCCCTGTAGGAGCGGTTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCt                                         <  2:1053669/150‑1 (MQ=255)
                                                                                                              tGCCGGCCCTTTCGCGGGTTAACCCCCTCCGACAGGATTATCGCTGCCCCGGTAGGGGCGGGTTTACCCGCGAAGAGGCCAGCCCCGGCAACCCAGAAACCCGACTGTGCGCACACCCCTACTGTTCATTCTGCTGCTACTCACCGccgg                                        >  2:6104919/1‑148 (MQ=38)
                                                                                                                 cGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGATTATCGCTGCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCaa                                     <  2:7665707/150‑1 (MQ=255)
                                                                                                                     ccTTTCGCGGGTAAACCCGCTCCTACAGGATTATCGCTGCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCaa                                 <  2:4535471/150‑1 (MQ=255)
                                                                                                                        ttCGGGGTTAACCCCGCGCCACCCGGATTAGCGCTGCCCCTGTTGGCGGGGTTTCACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGc                              <  1:85940/150‑1 (MQ=255)
                                                                                                                        ttCGCGGGTAAACCCGCTCCTACAGGATTATCGCTGCCCCTGGAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGc                              >  1:3577959/1‑150 (MQ=255)
                                                                                                                              ggTAAACCCGCTCCTACAGGATTATCGCTGCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCAcc                        >  1:465878/1‑150 (MQ=255)
                                                                                                                               gTAAACCCGCTCCTACAGGATTATCGCTGCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCAccc                       >  1:1779832/1‑150 (MQ=255)
                                                                                                                                   aCCCGCTCCTACTGGTTTAGCGCGCCCCCTTTAGGCGGGGTTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCTGGa                   <  1:3054626/150‑1 (MQ=255)
                                                                                                                                   aCCCGCTCCTACAGGATTATCGCTGCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCTGGa                   <  2:3901228/150‑1 (MQ=255)
                                                                                                                                     ccGCTCCTACAGGGGCAGCGATAATCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCTGGAgc                 >  1:2147707/1‑150 (MQ=37)
                                                                                                                                     ccGCTCCTACAGGATTATCGCTGCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCTGGAgc                 >  1:8475246/1‑150 (MQ=255)
                                                                                                                                     ccGCTCCTACAGGATTATCGCTGCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCTGGAgc                 >  1:3095751/1‑150 (MQ=255)
                                                                                                                                        cTCCTACAGGATTATCGCTGCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCTGGAGCGcc              >  2:7839023/1‑150 (MQ=255)
                                                                                                                                        cTCCAACGGTATTCTCGCGGCCCCTGTTGGTGGGGTTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCTGGAGCGcc              <  1:751846/150‑1 (MQ=255)
                                                                                                                                         tCCTACAGGATTATCGCTGCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCTGGa                   >  2:8800355/1‑144 (MQ=255)
                                                                                                                                         tCCTACAGGATTATCGCTGCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCTGGa                   <  1:8800355/144‑1 (MQ=255)
                                                                                                                                         tCCTACAGGATTATCGCTGCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCTGGAGCGCCt             >  2:4889415/1‑150 (MQ=255)
                                                                                                                                           cTACTGGATTATCGCTGCCCCTTTTGGTGCGGTTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCTGGAGCGCCTgg           <  1:3431026/150‑1 (MQ=255)
                                                                                                                                            tACCGGTTTCTCGCGGCCCCGGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCTGGAGCGCCTggg          <  1:2795570/150‑1 (MQ=39)
                                                                                                                                            tACAGGATTATCGCTGCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCTGGAGCGCCTggg          >  2:4051195/1‑150 (MQ=255)
                                                                                                                                             aCAGGATTATCGCTGCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCTGGAGCGCCTGggt         >  2:1103029/1‑150 (MQ=255)
                                                                                                                                              cAGTATTATCGCTGCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCTGGAGCGCCTGggtg        <  1:5670901/150‑1 (MQ=255)
                                                                                                                                                 gATTATCGCTGCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCTGGAGCGCCTGGGTGGTc     >  2:6757813/1‑150 (MQ=255)
                                                                                                                                                 gATTATCGCTGCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCTGGAGCGCCTGGGTGGTc     >  1:3300707/1‑150 (MQ=255)
                                                                                                                                                   ttATCGCTGCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCAGCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCTGGAGCGCCTGGGTGGTCCg   >  2:5482849/1‑150 (MQ=255)
                                                                                                                                                   ttATCGCTGCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCTGGAGCGCCTGGGTGGTCCg   >  1:5186698/1‑150 (MQ=255)
                                                                                                                                                   ttATCGCTGCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCTGGAGCGCCTGGGTGGTCCg   >  1:4619698/1‑150 (MQ=255)
                                                                                                                                                   ttATCGCTGCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCTGGAGCGCCTGGGTGGTCCg   >  2:7310790/1‑150 (MQ=255)
                                                                                                                                                   ttATCGCTGCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCTGGAGCGCCTGGGTGGTCCg   >  2:253319/1‑150 (MQ=255)
                                                                                                                                                   ttATCGCTGCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCTGGAGCGCCTGGGTGGTCCg   >  1:1689242/1‑150 (MQ=255)
                                                                                                                                                    tATCGCTGCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCTGGAGCGCCTGGGTGGTCCGa  <  2:1238399/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
GGCTGGTGTGAACCCGCCAGCGTTTCCACGCTGATTCGCTTGTAAGGCAATGAAAAAACGGGAACCTTCGGGTTCCCGTTTTTTTTTATTCTGAATGTGTGTTGTCAGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGATTATCGCTGCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCTGGAGCGCCTGGGTGGTCCGA  >  NC_002947/2461708‑2462003

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