Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,461,853 T→G 6.3% intergenic (+112/‑69) 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,8530TG 6.3% 641.8 / 23.3 210intergenic (+112/‑69)gapB/PP_2150glyceraldehyde‑3‑phosphate dehydrogenase/thiamine biosynthesis lipoprotein
Reads supporting (aligned to +/- strand):  ref base T (84/108);  new base G (4/10);  total (92/118)
Fisher's exact test for biased strand distribution p-value = 4.02e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.90e-01

Only 100 of 219 total aligned reads displayed.

GGCTGGTGTGAACCCGCCAGCGTTTCCACGCTGATTCGCTTGTAAGGCAATGAAAAAACGGGAACCTTCGGGTTCCCGTTTTTTTTTATTCTGAATGTGTGTTGTCAGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGATTATCGCTGCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCTGGAGCGCCTGGGTGGTCCG  >  NC_002947/2461708‑2462002
                                                                                                                                                    |                                                                                                                                                     
ggCTGGTGTGAACCCGCCAGCGTTTCCACGCTGATTCGCTTGTAAGGCAATGAAAAAACGGGAACCTTCGGGTTCCCGTTTTTTTTTATTCTGAATGTGTGTTGTCAGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGATTATc                                                                                                                                                   <  1:8343460/150‑1 (MQ=255)
ggCTGGTGTGAACCCGCCAGCGTTTCCACGCTGATTCGCTTGTAAGGCAATGAAAAAACGGGAACCTTCGGGTTCCCGTTTTTTTTTATTCTGAATGTGTGTTGTCAGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGATTATc                                                                                                                                                   <  2:375415/150‑1 (MQ=255)
   tGGTGTGAACCCGCCAGCGTTTCCACGCTGATTCGCTTGTAAGGCAATGAAAAAACGGGAACCTTCGGGTTCCCGTTTTTTTTTATTCTGAATGTGTGTTGTCAGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGATTATCGCt                                                                                                                                                <  2:3813350/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:842027/150‑1 (MQ=255)
       gtgAACCCGCCAGCGTTTCCACGCTGATTCGCTTGTAAGGCAATGAAAAAACGGGAACCTTCGGGTTCCCGTTTTTTTTTATTCTGAATGTGTGTTGTCAGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGATTATCGCTGccc                                                                                                                                            <  1:662453/150‑1 (MQ=255)
         gAACCCGCCAGCGTTTCCACGCTGATTCGCTTGTAAGGCAATGAAAAAACGGGAACCTTCGGGTTCCCGTTTTTTTTTATTCTGAATGTGTGTTGTCAGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGATTATCGCTGCCCCt                                                                                                                                          <  1:6878155/150‑1 (MQ=255)
          aaCCCGCCAGCGTTTCCACGCTGATTCGCTTGTAAGGCAATGAAAAAACGGGAACCTTCGGGTTCCCGTTTTTTTTTATTCTGAATGTGTGTTGTCAGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGGTTATCGCTGccccgg                                                                                                                                         >  1:6378820/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)
             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)
                  aGCGTTTCCACGCTGATTCGCTTGTAAGGCAATGAAAAAACGGGAACCTTCGGGTTCCCGTTTTTTTTTATTCTGAATGTGTGTTGTCAGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGATTATCGCTGCCCCGGTAGGAGCg                                                                                                                                 >  2:1185346/1‑150 (MQ=255)
                   gCGTTTCCACGCTGATTCGCTTGTAAGGCAATGAAAAAACGGGAACCTTCGGGTTCCCGTTTTTTTTTATTCTGGATGTGTGTTGTCAGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGATTATCGCTGCCCCTGTAGGAGCgg                                                                                                                                <  1:4762639/150‑1 (MQ=255)
                   gCGTTTCCACGCTGATTCGCTTGTAAGGCAATGAAAAAACGGGAACCTTCGGGTTCCCGTTTTTTTTTATTCTGAATGTGTGTTGTCAGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGATTATCGCTGCCCCTGTAGGAGCgg                                                                                                                                <  2:4348116/150‑1 (MQ=255)
                   gCGTTTCCACGCTGATTCGCTTGTAAGGCAATGAAAAAACGGGAACCTTCGGGTTCCCGTTTTTTTTTATTCTGAATGTGTGTTGTCAGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGATTATCGCTGCCCCTGTAGGAGCgg                                                                                                                                <  1:2592521/150‑1 (MQ=255)
                   gCGTTTCCACGCTGATTCGCTTGTAAGGCAATGAAAAAACGGGAACCTTCGGGTTCCCGTTTTTTTTTATTCTGAATGTGTGTTGTCAGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGATTATCGCTGCCCCGGTAGGAGCgg                                                                                                                                >  2:2640390/1‑150 (MQ=255)
                    cGTTTCCACGCTGATTCGCTTGTAAGGCAATGAAAAAACGGGAACCTTCGGGTTCCCGTTTTTTTTTTATTCTGAATGTGTGTTGTCAGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGGGCAGCGATAATCCTGTAGGAGCgg                                                                                                                                <  1:6328012/150‑1 (MQ=25)
                    cGTTTCCACGCTGATTCGCTTGTAAGGCAATGAAAAAACGGGAACCTTCGGGTTCCCGTTTTTTTTTATTCTGAATGTGTGTTGTCAGTGCCGGCCCTTTCGCGGGTTAAACCGCTCCGACAGGATTATCGCAGCCCCAGTTGGAGCggg                                                                                                                               >  2:751846/1‑150 (MQ=39)
                     gTTTCCACGCTGATTCGCTTGTAAGGCAATGAAAAAACGGGAACCTTCGGGTTCCCGTTTTTTTTTATTCTGAATGTGTGTTGTCAGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGATTATCGCTGCCCCTGTAGGAGCGGGt                                                                                                                              <  1:2121706/150‑1 (MQ=255)
                         ccACGCTGATTCGCTTGTAAGGCAATGAAAAAACGGGCACCTTCGGGGTCCCGTTTTTTTTTATTCTGAATGTGTGTTGTCAGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGATTATCGCTGCCCCTGTAGGAGCGGGTTTAc                                                                                                                          <  2:5239031/150‑1 (MQ=255)
                         ccACGCTGATTCGCTTGTAAGGCAATGAAAAAACGGGAACCTTCGGGTTCCCGTTTTTTTTTATTCTGAATGTGTGTTGTCAGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGGGCAGCGATAATCCTGTAGGAGCGGGTTTAc                                                                                                                          <  1:3230289/150‑1 (MQ=25)
                         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)
                                 aTTCGCTTGTAAGGAAATGAAAAAACGGTAACCTTCGGTTTCCCTTTTTTTTTTATTCTGAATGTGTGTTGTCAGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGGGCAGCGATAATCCTGTAGGAGCGGGTTTACCCGCGaaa                                                                                                                  <  2:1611594/150‑2 (MQ=18)
                                  ttCGCTTGTAAGGCAATGAAAAAACGGGACCCTTCGGGTTCCCTTTTTTTTTTATTCTGAATGTGTGTTGTCAGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGATTATCGCTGCCCCGGTAGGAGCGGGTTTACCCGCGAaga                                                                                                                 <  2:5010016/150‑1 (MQ=255)
                                  tgcGCTTGTAAGGCATTGAAAAAACGGGAACCTTCGGGTACCCGTTTTTTTTTATTCTGACTGTGTGTTGTCAGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGATTATCGCTGCCCCTGTTGGTGCGGGTTTACCCGCGAaga                                                                                                                 <  2:6576902/148‑1 (MQ=39)
                                          tAAGGCAATGAAAAAACGGGAACCTTCGGGTTCCCGTTTTTTTTTATTCTGAATGTGTGTTGTCAGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGATTATCGCTGCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCa                                                                                                         >  2:8542981/1‑150 (MQ=255)
                                          tAAGGCAATGAAAAAACGGGAACCTTCGGGTTCCCGTTTTTTTTTATTCTGAATGTGTGTTGTCAGTGCCGGCCCTTTCGCGGGTAAAACCGCTCCGACAGGGTTATCGCTGCCCCTGGAGGAGCGGGTTTACCCGCGCACGGGCCAGCa                                                                                                         >  2:5526407/1‑150 (MQ=38)
                                              gCAATTAAAAAACGGGAACCTTCGGGTTCCCTTTTTTTTTTTATTCTTATTGTTTGTTGTCTGTGCCGGCCCTTTCGCGGTTAAACCCGCTCCTACAGGATTATCGCTGCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACCg                                                                                                      <  2:332793/150‑1 (MQ=37)
                                                aaTGAAAAAACGGGAACCTTCGGGTTCCCGTTTTTTTTTTATTCTGAATGTGTGTTGTCAGTGCCGGCCCTTTCGCGGGTAAAACCGCTCCTACAGGATTATCGCTGCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACCGGc                                                                                                    >  2:680597/1‑150 (MQ=255)
                                                 aaaaaaaaaCGGGAACCCGAAGGTTCCCGTTTTTTTTTATTCTGAATGTGTGTTGTCAGTGCCGGCCCTTTCGCGGGTTAAACCGCGCCGACAGGGGTATAGCTGCCCCGGGAGGGGGGGGGTTTCCCCCGGAGCGGCCCGCCCCgggca                                                                                                  >  1:1116782/4‑147 (MQ=255)
                                                     aaaaaCGGGAACCTTCGGGTTCCCGTTTTTTTTTATTCTGAATGTGTGTTGTCAGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGATTATCGCTGCCCCGGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACCGGCAACCga                                                                                              >  1:2895297/1‑150 (MQ=255)
                                                      aaaaCGGGAACCTTCGGGTTCCCGTTTTTTTTTATTCTGAATGTGTGTTGTCAGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGATTATCGCTGCCCCt                                                                                                                                          >  1:6941502/1‑105 (MQ=255)
                                                      aaaaCGGGAACCTTCGGGTTCCCGTTTTTTTTTATTCTGAATGTGTGTTGTCAGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGATTATCGCTGCCCCt                                                                                                                                          <  2:6941502/105‑1 (MQ=255)
                                                       aaaCGGGAACCTTCGGGTTCCCGTTTTTTTTTATTCTGAATGTGTGTTGTCAGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACCGGATTATCGCTGCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACCGGCAACCgaga                                                                                            <  2:654074/150‑1 (MQ=255)
                                                       aaaCGGGAACCTTCGGGTTCCCGTTTTTTTTTATTCTGAATGTGTGTTGTCAGTGCCGGCCCTTTCGCGGGGAAAACCCCGCCTACAGGATTATAGCTGCCCCAGTAGGTGGGGGTTTACCCCGGCCGAGGGCACCACCCGCAAccccga                                                                                            >  1:4467219/1‑146 (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)
                                                                       gTTCCCGTTTTTTTTTATTCTGAATGTGTGTTGTCAGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCAACACGAATATCCCAGCCCCCGGACGACCCGGGTTACCCCCCAAGAGGCCCCCACCCGCAACCCAGCAACCCAACAGGgcgcc                                                                            >  2:4459054/1‑149 (MQ=19)
                                                                       gTTCCCGTTTTTTTTTATTCTGAATGTGTGTTGTCAGTGCCGGCCCTTTCGCGGGTAAAACCGCTCCAACAGGATTATAGCAGCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCAACTGTgcgcc                                                                            >  2:4697553/1‑149 (MQ=38)
                                                                                 ttttttttATTCTTATTGTGTGTTTTCTGGGCGGGCCCTTGCGGGGTTAACCCCGCTCCTCCAGGTTTTGCGCTGCCCCTTTAGGTGGGGGTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTa                                                                    <  1:1618216/150‑1 (MQ=255)
                                                                                 ttttttttATTCTGATTGTGTGTTGTCAGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGATTATCGCTGCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTa                                                                    <  2:359238/150‑1 (MQ=255)
                                                                                   ttttttttTCTGAATGTGTGTTGTCTGTGCCGGCCCTTTCGCGGTTAAACCCGCTCATACCGGATTATCGCTGCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACt                                                                  <  1:8769180/150‑1 (MQ=39)
                                                                                   ttttttATTCTGAATGTGTGTTGTCAGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACCGGATTATCGCTGCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACt                                                                  <  2:1614218/150‑1 (MQ=255)
                                                                                   ttttttATACTGAAGGTGGGTTGTCAGTGCCGGCCATTACGCGGGTAAACCCGCTCCAACCGTTTTATCGCTGCCCCGGTAGGCGCGGTTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACt                                                                  <  1:4338145/150‑1 (MQ=37)
                                                                                    tttttATTCTGAATGTGTGTTGTCAGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGATTATCGCTGCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTg                                                                 >  2:7201861/1‑150 (MQ=255)
                                                                                                  tgtgtgTTGTCTGTGCCGGCCCTTTCGCGGGTAAACCCGCTCATCCTGTTTTATCGCTGCCCCTGTTGGAGCGGTTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCtgctgc                                                    <  1:3395158/150‑1 (MQ=38)
                                                                                                    tgtgtTGTCAGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGATTATCGCTGCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTa                                                  <  1:5861383/150‑1 (MQ=255)
                                                                                                       gttgtCAGTGCCGGCCCTTTCGCGGGTAAAACCCCTCCGACAGGATTATAGCTGCCCCGGGAGGAGCGGGTTTTCCCGCGAAGAGGCCAGCCCCCGCAACCGAGAAACCCGACAGTGCGCACAGCCCGACTGTTCATTCTGCGGCTACTc                                               >  2:3131677/1‑150 (MQ=255)
                                                                                                        ttgtCAGTGCCGGCCCTTTCGCGGGTAAAACCCCTCCTACAGGATTATCGCTGCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTAATGTTCATTCTGCTGCTACTCa                                              >  2:2582498/1‑150 (MQ=255)
                                                                                                              gTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGTTTTATCGCTGCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCt                                        <  2:1053422/150‑1 (MQ=255)
                                                                                                               tGCCGGCCCTTTCGCGGGTTAACCCCCTCCGACAGGATTATCGCTGCCCCGGTAGGGGCGGGTTTACCCGCGAAGAGGCCAGCCCCGGCAACCCAGAAACCCGACTGTGCGCACACCCCTACTGTTCATTCTGCTGCTACTCACCGccgg                                       >  2:6104919/1‑148 (MQ=38)
                                                                                                                  cGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGATTATCGCTGCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCaa                                    <  2:7665707/150‑1 (MQ=255)
                                                                                                                      ccTTTCGCGGTTAAACCCGCTCCTACAGGATTATCGCTGCCCCTGTAGTAGCGGGTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCaa                                <  2:2452953/150‑1 (MQ=255)
                                                                                                                        tttCGCGGGTAAACCCCCTCCAACACGAATATCGCATCCCCAGTAGGAGCCGGTTTACCCGCGAAGAGGCCAGCACCGGCAAACGAGAAACCCAACAGTGCGCACAGCCCGACTGTTCATTCTGCTGCTACTCACCCCCTGCAACCaacg                              >  2:107274/1‑148 (MQ=37)
                                                                                                                         ttCGCGGGTAAACCCGCTCCTACAGGATTATCGCTGCCCCTGGAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGc                             >  1:3577959/1‑150 (MQ=255)
                                                                                                                             cgGGCAACCCCGCTCCTCCTCGATTATCGCCGCCCCTGTAGGAGCGGTTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCa                         <  1:1958739/150‑1 (MQ=39)
                                                                                                                               ggTAAACCCGCTCCTACAGGATTATCGCTGCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCAcc                       >  1:465878/1‑150 (MQ=255)
                                                                                                                                gTAAACCCGCTCCTACAGGATTATCGCTGCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCAccc                      >  1:1779832/1‑150 (MQ=255)
                                                                                                                                    aCCCGCTCCTACTGGTTTAGCGCGCCCCCTTTAGGCGGGGTTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCTGGa                  <  1:3054626/150‑1 (MQ=255)
                                                                                                                                    aCCCGCTCCTACAGGATTATCGCTGCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCTGGa                  <  1:7857525/150‑1 (MQ=255)
                                                                                                                                     cccGCTCCTACAGGATTATCGCTGCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCTGGAg                 >  1:8025103/1‑150 (MQ=255)
                                                                                                                                      ccGCTCCTACAGGGGCAGCGATAATCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCTGGAgc                >  1:2147707/1‑150 (MQ=37)
                                                                                                                                      ccGCTCCTACAGGATTATCGCTGCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCTGGAgc                >  1:3095751/1‑150 (MQ=255)
                                                                                                                                         cTCCTACCGGATTATCGCGGCCCCTGTCGGTGGGGTTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCTGGAGCGcc             <  1:859436/150‑1 (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)
                                                                                                                                             tACCGGATTATCGCGGCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTCCTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCTGGAGCGCCTggg         <  1:2795552/150‑1 (MQ=255)
                                                                                                                                             tACAGGATTATCGCTGCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCTGGAGCGCCTggg         >  1:5141439/1‑150 (MQ=255)
                                                                                                                                             tACAGGATTATCGCTGCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCTGGAGCGCCTggg         >  2:4051195/1‑150 (MQ=255)
                                                                                                                                             tACAGGATTATCGCTGCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCTGGAGCGCCTggg         >  2:3961548/1‑150 (MQ=255)
                                                                                                                                               cAGTATTATCGCTGCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCTGGAGCGCCTGggtg       <  1:5670901/150‑1 (MQ=255)
                                                                                                                                                aGGATTATCGCTGCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCAGGCGCTCCAGggtgg      <  1:8537992/150‑1 (MQ=255)
                                                                                                                                                  gATTATCGCTGCCCCTGTTGGAGCGGTTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGCACCACCCAGGCGCTCCAGggtgggg    <  1:7677746/150‑3 (MQ=255)
                                                                                                                                                  gATTATCGCTGCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCTGGAGCGCCTGGGTGGTc    >  2:6757813/1‑150 (MQ=255)
                                                                                                                                                   aTTATCGCTGCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCTGGAGCGCCTGGGTGGTcc   <  2:455375/150‑1 (MQ=255)
                                                                                                                                                    ttATCGCTGCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCAGCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCTGGAGCGCCTGGGTGGTCCg  >  2:5482849/1‑150 (MQ=255)
                                                                                                                                                    ttATCGCTGCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCTGGAGCGCCTGGGTGGTCCg  >  1:4619698/1‑150 (MQ=255)
                                                                                                                                                    ttATCGCTGCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCTGGAGCGCCTGGGTGGTCCg  >  2:7310790/1‑150 (MQ=255)
                                                                                                                                                    ttATCGCTGCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCTGGAGCGCCTGGGTGGTCCg  >  1:1689242/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
GGCTGGTGTGAACCCGCCAGCGTTTCCACGCTGATTCGCTTGTAAGGCAATGAAAAAACGGGAACCTTCGGGTTCCCGTTTTTTTTTATTCTGAATGTGTGTTGTCAGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGATTATCGCTGCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCTGGAGCGCCTGGGTGGTCCG  >  NC_002947/2461708‑2462002

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: