Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,953,314 G→C 11.6% intergenic (‑128/+33) oguA ← / ← PP_2585 8‑oxoguanine deaminase/alpha‑ketoglutaric semialdehyde dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,953,3140GC11.6% 488.0 / 46.4 164intergenic (‑128/+33)oguA/PP_25858‑oxoguanine deaminase/alpha‑ketoglutaric semialdehyde dehydrogenase
Reads supporting (aligned to +/- strand):  ref base G (68/77);  new base C (11/8);  total (79/85)
Fisher's exact test for biased strand distribution p-value = 4.66e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.31e-01

Only 100 of 166 total aligned reads displayed.

GCAGGGTTTGTAGATGACATGGTGCAAGCCTGCGCAATACGCATCGGCGATTACCGATGTCATAGGCGCTCGGTATTGTAAGCAAGCTTTTTACTGATTGTATACAGCTCGAAGTGCCCGTGCCGGCCTCTTGGCGGCTAGAGCCGCGAAGAGGCCGGCACGGGTTGCAGCCATGAAAGATCACACCGCCTCGCGGCTGCTCACACCCTCCAGCAACCGGGCAATCCCCAGCGGGTTAGCGTTTTTCAACGCATCCGGCAGCAGTGCGTCAGGGTAGTTCTGGAAGCACACCGGGCG  >  NC_002947/2953167‑2953463
                                                                                                                                                   |                                                                                                                                                     
gCAGGGTTTGTAGATGACATGGTGCAAGCCTGCGCAATACGCATCGGCGATTACCGATGTCATAGGCGCTCGGTATTGTAAGCAAGCTTTTTACTGATTGTATACAGCTCGAAGTGCCCGTGCCGGCCTCTTGGCGGCTAGAGCCGCGaa                                                                                                                                                     <  1:4760389/150‑1 (MQ=255)
gCAGGGTTTGTAGATGACATGGTGCAAGCCTGCGCAATACGCATCGGCGATTACCGATGTCATAGGCGCTCGGTATTGTAAGCAAGCTTTTTACTGATTGTATACAGCTCGAAGTGCCCGTGCCGGCCTCTTGGCGGCTAGAGCCGCGaa                                                                                                                                                     <  2:1143452/150‑1 (MQ=255)
gCAGGGTTTGTAGATGACATGGTGCAAGCCTGCGCAATACGCATCGGCGATTACCGATGTCATAGGCGCTCGGTATTGTAAGCAAGCTTTTTACTGATTGTATACAGCTCGAAGTGCCCGTGCCGGCCTCTTGGCGGCTAGAGCCGCGaa                                                                                                                                                     <  2:4363127/150‑1 (MQ=255)
gCAGGGTTTGCAGATGACATGGTGCAAGCCTGCGCAATACGCATCGGCGATTACCGATGTCATAGGCGCTCGGTATTGTAAGCAAGCTTTTTACTGATTGTATACAGCTCGAAGTGCCCGTGCCGGCCTCTTGGCGGCTAGAGCCGCGaa                                                                                                                                                     <  1:517538/150‑1 (MQ=255)
 cAGGGTTTGTAGATGACATGGTGCAAGCCTGCGCAATACGCATCGGCGATTACCGATGTCATAGGCGCTCGGTATTGTAAGCAAGCTTTTTACTGATTGTATACAGCTCGAAGTGCCCGTGCCGGCCTCTTGGCGGCTAGAGCCGCGAag                                                                                                                                                    >  1:2618078/1‑150 (MQ=255)
  aGGGTTTGTAGATGACATGGTGCAAGCCTGCGCAATACGCATCGGCGATTACCGATGTCATAGGCGCTCGGTATTGTAAGCAAGCTTTTTACTGATTGTATACAGCTCGAAGTGCCCGTGCCGGCCTCTTGGCGGCTAGAGCCGCGAaga                                                                                                                                                   <  2:4072940/150‑1 (MQ=255)
  aGGGTTTGTAGATGACATGGTGCAAGCCTGCGCAATACGCATCGGCGATTACCGATGTCATAGGCGCTCGGTATTGTAAGCAAGCTTTTTACTGATTGTATACAGCTCGAAGTGCCCGTGCCGGCCTCTTGGCGGCTAGAGCCGCGAaga                                                                                                                                                   <  2:4073032/150‑1 (MQ=255)
  aGGGTTTGTAGATGACATGGTGCAAGCCTGCGCAATACGCATCGGCGATTACCGATGTCATAGGCGCTCGGTATTGTAAGCAAGCTTTTTACTGATTGTATACAGCTCGAAGTGCCCGTGCCGGCCTCTTCGCGGCTCTAGCCGCCAaga                                                                                                                                                   <  2:2637089/150‑1 (MQ=255)
  aGGGTTTGTAGATGACATGGTGCAAGCCTGCGCAATACGCATCGGCGATTACCGATGTCATAGGCGCTCGGTATTGTAAGCAAGCTTTTTACTGATTGTATACAGCTCGAAGTACCCGTGCCGGCCTCTTCGCGGCTCTAGCCGCCAaga                                                                                                                                                   <  2:131117/150‑1 (MQ=39)
   gggTTTGTAGATGACATGGTGCAAGCCTGCGCAATACGCATCGGCGATTACCGATGTCATAGGCGCTCGGTATTGTAAGCAAGCTTTTTACTGATTGTATACAGCTCGAAGTGCCCGTGCCGGCCTCTTGGCGGCTAGAGCCGCGAagag                                                                                                                                                  <  1:1242844/150‑1 (MQ=255)
      tttGTAGATGACATGGTGCAAGCCTGCGCAATACGCATCGGCGATTACCGATGTCATAGGCGCTCGGTATTGTAAGCAAGCTTTTTACTGATTGTATACAGCTCGAAGTGCCCGTGCCGGCCTCTTCGCGGCTCTAGCCGCCAagag                                                                                                                                                  <  2:5021518/147‑1 (MQ=255)
      tttGTAGATGACATGGTGCAAGCCTGCGCAATACGCATCGGCGATTACCGATGTCATAGGCGCTCGGTATTGTAAGCAAGCTTTTTACTGATTGTATACAGCTCGAAGTGCCCGTGCCGGCCTCTTCGCGGCTCTAGCCGCCAagag                                                                                                                                                  >  1:5021518/1‑147 (MQ=255)
       ttGTAGATGACATGGTGCAAGCCTGCGCAATACGCATCGGCGATTACCGATGTCATAGGCGCTCGGTATTGTAAGCAAGCTTTTTACTGATTGTATACAGCTCGAAGTGCCCGTGCCGGCCTCTTGGCGGCTAGAGCCGCGAAGAGGCCg                                                                                                                                              <  1:1878930/150‑1 (MQ=255)
             aTGACATGGTGCAAGCCTGCGCAATACGCATCGGCGATTACCGATGTCATAGGCGCTCGGTATTGTAAGCAAGCTTTTTACTGATTGTATACAGCTCGAAGTGCCCGTGCCGGCCTCTTGGCGGCTAGAGCCGCGAAGAGGCCGGCACgg                                                                                                                                        <  2:1134645/150‑1 (MQ=255)
              tGACATGGTGCAAGCCTGCGCAATACGCATCGGCGATTACCGATGTCATAGGCGCTCGGTATTGTAAGCAAGCTTTTTACTGATTGTATACAGCTCGAAGTGCCCGTGCCGGCCTCTTGGCGGCTAGAGCCGCGAAGAGGCCGGCACggg                                                                                                                                       >  2:2252512/1‑150 (MQ=255)
                aCATGGTGCAAGCCTGCGCAATACGCATCGGCGATTACCGATGTCATAGGCGCTCGGTATTGTAAGCAAGCTTTTTACTGATTGTATACAGCTCGAAGTGCCCGTGCCGGCCTCTTGGCGGCTAGAGCCGCGAAGAGGCCGGCACGGGtt                                                                                                                                     <  1:1883645/150‑1 (MQ=255)
                     gTGCAAGCCTGCGCAATACGCATCGGCGATTACCGATGTCATAGGCGCTCGGTATTGTAAGCAAGCTTTTTACTGATTGTATACAGCTCGAAGTGCCCGTGCCGGCCTCTTGGCGGCTAGAGCCGCGAAGAGGCCGGCACGGGTTGCAGc                                                                                                                                >  2:3725629/1‑150 (MQ=255)
                       gCAAGCCTGCGCAATACGCATCGGCGATTACCGATGTCATAGGCGCTCGGTATTGTAAGCAAGCTTTTTACTGATTGTATACAGCTCGAAGTGCCCGTGCCGGCCTCTTCGCGGCTCTAGCCGCCaa                                                                                                                                                     <  2:4483130/127‑1 (MQ=255)
                       gCAAGCCTGCGCAATACGCATCGGCGATTACCGATGTCATAGGCGCTCGGTATTGTAAGCAAGCTTTTTACTGATTGTATACAGCTCGAAGTGCCCGTGCCGGCCTCTTCGCGGCTCTAGCCGCCaa                                                                                                                                                     >  1:4483130/1‑127 (MQ=255)
                             cTGCGCAATACGCATCGGCGATTACCGATGTCATAGGCGCTCGGTATTGTAAGCAAGCTTTTTACTGATTGTATACAGCTCGAAGTGCCCGTGCCGGCCTCTTGGCGGCTAGAGCCGCGAAGAGGCCGGCACGGGTTGCAGCCATGAAAg                                                                                                                        >  1:3037428/1‑150 (MQ=255)
                              tGCGCAATACGCATCGGCGATTACCGATGTCATAGGCGCTCGGTATTGTAAGCAAGCTTTTTACTGATTGTATACAGCTCGAAGTGCCCGTGCCGGCCTCTTGGCGGCTAGAGCCGCGAAGAGGCCGGCACGGGTTGCAGCCATGAAAGa                                                                                                                       >  2:2631142/1‑150 (MQ=255)
                                cgcAATACGCATCGGCGATTACCGATGTCATAGGCGCTCGGTATTGTAAGCAAGCTTTTTACTGATTGTATACAGCTCGAAGTGCCCGTGCCGGCCTCTTGGCGGCTAGAGCCGCGAAGAGGCCGGCACGGGTTGCAGCCATGAAAGATc                                                                                                                     <  2:86770/150‑1 (MQ=255)
                                       cGCATCGGCGATTACCGATGTCATAGGCGCTCGGTATTGTAAGCAAGCTTTTTACTGATTGTATACAGCTCGAAGTGCCCGTGCCGGCCTCTTGGCGGCTAGAGCCGCGAAGAGGCCGGCACGGGTTGCAGc                                                                                                                                >  1:3814911/1‑132 (MQ=255)
                                             ggCGATTACCGATGTCATAGGCGCTCGGTATTGTAAGCAAGCTTTTTACTGATTGTATACAGCTCGAAGTGCCCGTGCCGGCCTCTTGGCGGCTAGAGCCGCGAAGAGGCCGGCACGGGTTGCAGCCATGAAAGATCACACCGCCTcgcg                                                                                                        >  1:119800/1‑150 (MQ=255)
                                                taTTACCGATGTCATAGTCGCTCGGTATTTTAAGCAAGCTTTTTACTTATTGTATACTGCTCGAAGTGCCCTTGCCGGCCTCTTGGCGGCTAGTGCCGCGAAGAGGCCGGCACGGGTTGCAGCCATGAAAGATCACACCGCCTCGCGgct                                                                                                     <  1:2065638/149‑1 (MQ=255)
                                                gATTACCGATGTCATAGGCGCTCGGTATTGTAAGCAAGCTTTTTACTGATTGTATACAGCTCGAAGTGCCCGTGCCGGCCTCTTGGCGGCTAGAGCCGCGAAGAGGCCGGCACGGGTTGCAGCCATGAAAGATCACACCGCCTCGCGgct                                                                                                     <  1:204172/150‑1 (MQ=255)
                                                   tACCGATGTCATAGGCGCTCGGTATTGTAAGCAAGCTTTTTACTGATTGTATACAGCTCGAAGTGCCCGTGCCGGCCTCTTGGCGGCTAGAGCCGCGAAGAGGCCGGCACGGGTTGCAGCCATGAAAGATCACACCGCCTCGCGgctgct                                                                                                  <  2:3661913/150‑1 (MQ=255)
                                                     ccGATGTCATAGGCGCTCGGTATTGTAAGCAAGCTTTTTACTGATTGTATACAGCTCGAAGTGCCCGTGCCGGCCTCTTGGCGGCTAGAGCCGCGAAGAGGCCGGCACGGGTTGCAGCCATGAAAGATCACACCGCCTCGCGGCTGCTca                                                                                                <  1:4204989/150‑1 (MQ=255)
                                                      cGATGTCATAGGCGCTCGGTATTGTAAGCAAGCTTTTTACTGATTGTATACAGCTCGAAGTGCCCGTGCCGGCCTCTTGGCGGCTAGAGCCGCGAAGAGGCCGGCACGGGTTGCAGCCATGAAAGATCACACCGCCTCGCGGCTGCTcac                                                                                               >  1:3517808/1‑150 (MQ=255)
                                                       gATGTCATAGGCGCTCGGTATTGTAAGCAAGCTTTTTACTGATTGTATACAGCTCGAAGTGCCCGTGCCGGCCTCTTGGCGGCTAGAGCCGCGAAGAGGCCGGCACGGGTTGCAGCCATGAAAGATCACACCGCCTCGCGGCTGCTcaca                                                                                              >  2:5409090/1‑150 (MQ=255)
                                                            cATAGGCGCTCGGTATTGTAAGCAAGCTTTTTACTGATTGTATACAGCTCGAAGTGCCCGTGCCGGCCTCTTGGCGGCTAGAGCCGCGAAGAGGCCGGCACGGGTTGCAGCCATGAAAGATCACACCGCCTCGCGGCTGCTCACACCCTc                                                                                         <  2:4730053/150‑1 (MQ=255)
                                                             aTAGGCGCTCGGTATTGTAAGCAAGCTTTTTACTGATTGTATACAGCTCGAAGTGCCCGTGCCGGCCTCTTGGCGGCTAGAGCCGCGAAGAGGCCGGCACGGGTTGCAGCCATGAAAGATCACACCGCCTCGCGGCTGCTCACACCCTcc                                                                                        >  1:1492942/1‑150 (MQ=255)
                                                              tAGGCGCTCGGTATTGTAAGCAAGCTTTTTACTGATTGTATACAGCTCGAAGTGCCCGTGCCGGCCTCTTGGCGGCTAGAGCCGCGAAGAGGCCGGCACGGGTTGCAGCCATGAAAGATCACACCGCCTCGCGGCTGCTCACACCCTCCa                                                                                       <  1:2955739/150‑1 (MQ=255)
                                                              gaGGCGCTCGGTATTGTAAGCAAGCTTTTTACTGATTGTATACAGCTCGAAGTGCCCGTGCCGGCCTCTTGGCGGCTAGAGCCGCGAAGAGGCCGGCACGGGTTGCAGCCATGAAAGATCACACCGCCTCGCGGCTGCTCACACCCTCCa                                                                                       <  1:4626275/149‑1 (MQ=255)
                                                               aGGCGCTCGGTATTGTAAGCAAGCTTTTTACTGATTGTATACAGCTCGAAGTGCCCGTGCCGGCCTCTTGGCGGCTAGAGCCGCGAAGAGGCCGGCACGGGTTGCAGCCATGAAAGATCACACCGCCTCGCGGCTGCTCACACCCTCCAg                                                                                      <  2:4668473/150‑1 (MQ=255)
                                                                ggCGCTCGGTATTGTAAGCAAGCTTTTTACTGATTGTATACAGCTCGAAGTGCCCGTGCCGGCCTCTTGGCGGCTAGAGCCGCGAAGAGGc                                                                                                                                                <  2:3331205/91‑1 (MQ=255)
                                                                ggCGCTCGGTATTGTAAGCAAGCTTTTGACTGATTGTATACAGCTCGAAGTGCCCGTGCCGGCCTCTTGGCGGCTAGAGCCGCGAAGAGGCCGGCACGGGTTGCAGCCATGAAAGATCACACCGCCTCGCGGCTGCTCACACCCTCCAGc                                                                                     <  2:3408010/150‑1 (MQ=255)
                                                                 gcgcTCGGTATTGTAAGCAAGCTTTTTACTGATTGTATACAGCTCGAAGTGCCCGTGCCGGCCTCTTGGCGGCTAGAGCCGCGAAGAGGCCGGCACGGGTTGCAGCCATGAAAGATCACACCGCCTCGCGGCTGCTCACAcc                                                                                            >  1:2165135/1‑142 (MQ=255)
                                                                   gcTCGGTATTGTAAGCAAGCTTTTTACTGATTGTATACAGCTCGAAGTGCCCGTGCCGGCCTCTTGGCGGCTAGAGCCGCGAAGAGGCCGGCACGGGTTGCAGCCATGAAAGATCACACCGCCTCGCGGCTGCTCACACCCTCCAGCAAc                                                                                  >  1:2779792/1‑150 (MQ=255)
                                                                      cGGTATTGTAAGCAAGCTTTTTACTGATTGTATACAGCTCGAAGTGCCCGTGCCGGCCTCTTGGCGGCTAGAGCCGCGAAGAGGCCGGCACGGGTTGCAGCCATGAAAGATCACACCGCCTCGCGGCTGCTCACACCCTCCAGCAACCgg                                                                               <  1:5083083/150‑1 (MQ=255)
                                                                      cGGTATTGTAAGCAAGCTTTTTACTGATTGTATACAGCTCGAAGTGCCCGTGCCGGCCTCTTGGCGGCTAGAGCCGCGAAGAGGCCGGCACGGGTTGCAGCCATGAAAGATCACACCGCCTCGCGGCTGCTCACACCCTCCAGCAACCgg                                                                               <  2:3273679/150‑1 (MQ=255)
                                                                      cGGTATTGTAAGCAAGCTTTTTACTGATTGTATACAGCTCGAAGTGCCCGTGCCGGCCTCTTGGCGGCTAGAGCCGCGAAGAGGCCGGCACGGGTTGCAGCCATGAAAGATCACACCGCCTCGCGGCTGCTCACACCCTCCAGCAACCgg                                                                               <  1:3826410/150‑1 (MQ=255)
                                                                         tATTGTAAGCAAGCTTTTTACTGATTGTATACAGCTCGAAGTGCCCGTGCCGGCCTCTTGGCGGCTAGAGCCGCGAAGAGGCCGGCACGGGTTGCAGCCATGAAAGATCACACCGCCTCGCGGCTGCTCACACCCTCCAGCAACCGGGCa                                                                            >  2:2835011/1‑150 (MQ=255)
                                                                            tGTAAGCAAGCTTTTTACTGATTGTATACAGCTCGGAGTGCCCGTGCCGGCCTCTTGGCGGCTAGAGCCGCGAAGAGGCCGGCACGGGTTGCAGCCATGAAAGATCACACCGCCTCGCGGCTGCTCACACCCTCCAGCAACCGGGCAATc                                                                         <  2:3520592/150‑1 (MQ=255)
                                                                              tAAGCAAGCTTTTTACTGATTGTATACAGCTCGAAGTGCCCGTGCCGGCCTCTTGGCGGCTAGAGCCGCGAAGAGGCCGGCACGGGTTGCAGCCATGAAAGATCACACCGCCTCGCGGCTGCTCACACCCTCCAGCAACCGGGCAATccc                                                                       >  2:2861311/1‑150 (MQ=255)
                                                                              tAAGCAAGCTTTTTACTGATTGTATACAGCTCGAAGTGCCCGTGCCGGCCTCTTGGCGGCTAGAGCCGCGAAGAGGCCGGCACGGGTTGCAGCCATGAAAGATCACACCGCCTCGCGGCTGCTCACACCCTCCAGCAACCGGGCAATccc                                                                       >  1:2088741/1‑150 (MQ=255)
                                                                              tAAGCAAGCTTTTTACTGATTGTATACAGCTCGAAGTGCCCGTGCCGGCCTCTTGGCGGCTAGAGCCGCGAAGAGGCCGGCACGGGTTGCAGCCATGAAAGATCACACCGCCTCGCGGCTGCTCACACCCTCCAGCAACCGGGCAATccc                                                                       >  1:3750526/1‑150 (MQ=255)
                                                                              tAAGCAAGCTTTTTACTGATTGTATACAGCTCGAAGTGCCCGTGCCGGCCTCTTGGCGGCTAGAGCCGCGAAGAGGCCGGCACGGGTTGCAGCCATGAAAGATCACACCGCCTCGCGGCTGCTCACACCCTCCAGCAACCGGGCAATccc                                                                       >  1:2872824/1‑150 (MQ=255)
                                                                                  caagcTTTTTACTGATTGTATACAGCTCGAAGTGCCCGTGCCGGCCTCTTGGCGGCTAGAGCCGCGAAGAGGCCGGCACGGGTTGCAGCCATGAAAGATCACACCGCCTCGCGGCTGCTCACACCCTCCAGCAACCgg                                                                               >  2:4627818/1‑138 (MQ=255)
                                                                                  caagcTTTTTACTGATTGTATACAGCTCGAAGTGCCCGTGCCGGCCTCTTGGCGGCTAGAGCCGCGAAGAGGCCGGCACGGGTTGCAGCCATGAAAGATCACACCGCCTCGCGGCTGCTCACACCCTCCAGCAACCgg                                                                               <  1:4627818/138‑1 (MQ=255)
                                                                                  caagcTTTTTACTGATTGTATACAGCTCGAAGTGCCCGTGCCGGCCTCTTGGCGGCTAGAGCCGCGAAGAGGCCGGCACGGGTTGCAGCCATGAAAGATCACACCGCCTCGCGGCTGCTCACACCCACCAGCAACCGGGCAATCCCCAGc                                                                   >  1:4782619/1‑150 (MQ=255)
                                                                                        ttttACTGATTGTATACAGCTCGAAGTGCCCGTGCCGGCCTCTTGGCGGCTAGAGCCGCGAAGAGGCCGGCACGGGTTGCAGCCATGAAAGATCACACCGCCTCGCGGCTGCTCACACCCTCCAGCAACCGGGCAATCCCCAGCGGGTTa                                                             <  2:4745248/150‑1 (MQ=255)
                                                                                         tttACTGATTGTATACAGCTCGAAGTGCCCGTGCCGGCCTCTTGGCGGCGAGAGCCGCGAAGAGGCCGGCACGGGTTGCAGCCATGAAAGATCACACCGCCTCGCGGCTGCTCACCCCCTCCAGCAACCGGGCAATCCCCAGCGGGTTAg                                                            >  1:162158/1‑150 (MQ=255)
                                                                                             cTGATTGTATACAGCTCGAAGTGCCCGTGCCGGCCTCTTGGCGGCTAGAGCCGCGAAGAGGCCGGCACGGGTTGCAGCCATGAAAGATCACACCGCCTCGCGGCTGCTCACACCCTCCAGCAACCGGGCAATCCCCAGCGGGTTAGCGtt                                                        >  1:4065587/1‑150 (MQ=255)
                                                                                              tGATTGTATACAGCTCGAAGTGCCCGTGCCGGCCTCTTGGCGGCTAGAGCCGCGAAGAGGCCGGCACGGGTTGCAGCCATGAAAGATCACACCGCCTCGCGGCTGCTCACACCCTCCAGCAACCGGGCAATCCCCAGCGGGTTAGCGttt                                                       >  2:3890651/1‑150 (MQ=255)
                                                                                              tGATTGTATACAGCTCGAAGTGCCCGTGCCGGCCTCTTGGCGGCTAGAGCCGCGAAGAGGCCGGCACGGGTTGCAGCCATGAAAGATCACACCGCCTCGCGGCTGCTCACACCCTCCAGCAACCGGGCAATCCCCAGCGGGTTAGCGttt                                                       >  2:2341616/1‑150 (MQ=255)
                                                                                                aTTGTATACAGCTCGAAGTGCCCGTGCCGGCCTCTTGGCGGCTAGAGCCGCGAAGAGGCCGGCACGGGTTGCAGCCATGAAAGATCACACCGCCTCGCGGCTGCTCACACCCTCCAGCAACCGGGCAATCCCCAGCGGGTTAGCGttttt                                                     >  2:323295/1‑150 (MQ=255)
                                                                                                aTTGTATACAGCTCGAAGTGCCCGTGCCGGCCTCTTGGCGGCTAGAGCCGCGAAGAGGCCGGCACGGGTTGCAGCCATGAAAGATCACACCGCCTCGCGGCTGCTCACACCCTCCAGCAACCGGGCAATCCCCAGCGGGTTAGCGttttt                                                     >  2:2824711/1‑150 (MQ=255)
                                                                                                    tataCAGCTCGAAGTGCCCGTGCCGGCCTCTTGGCGGCTAGAGCCGCGAAGAGGCCGGCACGGGTTGCAGCCATGAAAGATCACACCGCCTCGCGGCTGCTCACACCCTCCAGCAACCGGGCAATCCCCAGCGGGTTa                                                             >  1:5385309/1‑138 (MQ=255)
                                                                                                    tataCAGCTCGAAGTGCCCGTGCCGGCCTCTTGGCGGCTAGAGCCGCGAAGAGGCCGGCACGGGTTGCAGCCATGAAAGATCACACCGCCTCGCGGCTGCTCACACCCTCCAGCAACCGGGCAATCCCCAGCGGGTTa                                                             <  2:5385309/138‑1 (MQ=255)
                                                                                                    tataCAGCTCGAAGTGCCCGTGCCGGCCTCTTGGCGGCTAGAGCCGCGAAGAGGCCGGCACGGGTTGCAGCCATGAAAGATCACACCGCCTCGCGGCTGCTCACACCCTCCAGCAACCGGGCAATCCCCAGCGGGTTAGCGTTTTTCAAc                                                 <  1:5027977/150‑1 (MQ=255)
                                                                                                    tataCAGCTCGAAGTGCCCGTGCCGGCCTCTTGGCGGCTAGAGCCGCGAAGAGGCCGGCACGGGTTGCAGCCATGAAAGATCACACCGCCTCGCGGCTGCTCACACCCTCCAGCAACCGGGCAATCCCCAGCGGGTTAGCGTTTTTCAAc                                                 <  2:3813070/150‑1 (MQ=255)
                                                                                                       aCAGCTCGAAGCGCCCGTGCCGGCCTCTTGGCGGCTAGAGCCGCGAAGAGGCCGGCACGGGTTGCAGCCATGAAAGATCACACCGCCTCGCGGCTGCTCACACCCTCCAGCAACCGGGCAATCCCCAGCGGGTTAGCGTTTTTCAACGCa                                              >  2:375225/1‑150 (MQ=255)
                                                                                                         aGCTCGAAGTGCCCGTGCCGGCCTCTTGGCGGCTAGAGCCGCGAAGAGGCCGGCACGGGTTGCAGCCATGAAAGATCACACCGCCTCGCGGCTGCTCACACCCTCCAGCAACCGGGCAATCCCCAGCGGGTTAGCGTTTTTCAACGCATc                                            >  1:4905986/1‑150 (MQ=255)
                                                                                                            tCGAAGTGCCCGTGCCGGCCTCTTGGCGGCTAGAGCCGCGAAGAGGCCGGCACGGGTTGCAGCCATGAAAGATCACACCGCCTCGCGGCTGCTCACACCCTCCAGCAACCGGGCAATCCCCAGCGGGTTAGCGTTTTTCAACGCATCCgg                                         >  2:2075463/1‑150 (MQ=255)
                                                                                                            tCGAAGTGCCCGTGCCGGCCTCTTGGCGGCTAGAGCCGCGAAGAGGCCGGCACGGGTTGCAGCCATGAAAGATCACACCGCCTCGCGGCTGCTCACACCCTCCAGCAACCGGGCAATCCCCAGCGGGTTAGCGTTTTTCAACGCATCCgg                                         >  2:1961801/1‑150 (MQ=255)
                                                                                                            tCGAAGTGCCCGTGCCGGCCTCTTGGCGGCTAGAGCCGCGAAGAGGCCGGCACGGGTTGCAGCCATGAAAGATCACACCGCCTCGCGGCTGCTCACACCCTCCAGCAACCGGGCAATCCCCAGCGGGTTAGCGTTTTTCAACGCATCCgg                                         <  1:2359707/150‑1 (MQ=255)
                                                                                                             cGAAGTGCCCGTGCCGGCCTCTTGGCGGCTAGAGCCGCGAAGAGGCCGGCACGGGTTGCAGCCATGAAAGATCACACCGCCTCGCGGCTGCTCACACCCTCCAGCAACCGGGCAATCCCCAGCGGGTTAGCGTTTTTCAACGCATCCGgc                                        <  2:3788742/150‑1 (MQ=255)
                                                                                                               aaGTGCCCGTGCCGGCCTCTTGGCGGCTAGAGCCGCGAAGAGGCCGGCACGGGTTGCAGCCATGAAAGATCACACCGCCTCGCGGCTGCTCACACCCTCCAGCAACCGGGCAATCCCCAGCGGGTTAGCGTTTTTCAACGCATCCGgcag                                      >  2:1627014/1‑150 (MQ=255)
                                                                                                                  aacccGTGCCGGCCTCTTCGCGGCTCTAGCCGCCAAGAGGCCGGCACGGGTTGCAGCCATGAAAGATCACACCGCCTCGCGGCTGCTCACACCCTCCAGCAACCGGGCAATCCCCAGCGGGTTAGCGTTTTTCAACGCATCCGGcagcag                                   >  2:1134563/3‑150 (MQ=255)
                                                                                                                   gCCCGTGCCGGCCTCTTGGCGGCTAGAGCCGCGAAGAGGCCGGCACGGGTTGCAGCCATGAAAGATCACACCGCCTCGCGGCTGCTCACACCCTCCAGCAACCGGGCAATCCCCAGCGGGTTAGCGTTTTTCAACGCATCCGGCAGCAGt                                  <  2:3897106/150‑1 (MQ=255)
                                                                                                                   gCCCGTGCCGGCCTCTTGGCGGCTAGAGCCGCGAAGAGGCCGGCACGGGTTGCAGCCATGAAAGATCACACCGCCTCGCGGCTGCTCACACCCTCCAGCAACCGGGCAATCCCCAGCGGGTTAGCGTTTTTCAACGCATCCGGCAGCAGt                                  <  2:2480128/150‑1 (MQ=255)
                                                                                                                    cccGTGCCGGCCTCTTGGCGGCTAGAGCCGCGAAGAGGCCGGCACGGGTTGCAGCCATGAAAGATCACACCGCCTCGCGGCTGCTCACACCCTCCAGCAACCGGGCAATCCCCAGCGGGTTAGCGTTTTTCAACGCATCCGGCAGCAGTg                                 >  2:393374/1‑150 (MQ=255)
                                                                                                                           cGGCCTCTTGGCGGCTAGAGCCGCGAAGAGGCCGGCACGGGTTGCAGCCATGAAAGATCACACCGCCTCGCGgctg                                                                                                    >  1:3304838/1‑76 (MQ=255)
                                                                                                                           cGGCCTCTTGGCGGCTAGAGCCGCGAAGAGGCCGGCACGGGTTGCAGCCATGAAAGATCACACCGCCTCGCGGCTGCTCACACCCTCCAGCAACCGGGCAATCCCCAGCGGGTTAGCGTTTTTCAACGCATCCGGCAGCAGTGCGTCAgg                          >  2:4255458/1‑150 (MQ=255)
                                                                                                                           cGGCCTCTTGGCGGCTAGAGCCGCGAAGAGGCCGGCACGGGTTGCAGCCATGAAAGATCACACCGCCTCGCGGCTGCTCACACCCTCCAGCAACCGGGCAATCCCCAGCGGGTTAGCGTTTTTCAACGCATCCGGCAGCAGTGCGTCAgg                          <  1:2252512/150‑1 (MQ=255)
                                                                                                                           cGGCCTCTTGGCGGCTAGAGCCGCGAAGAGGCCGGCACGGGTTGCAGCCATGAAAGATCACACCGCCTCGCGGCTGCTCACACCCTCCAGCAACCGGGCAATCCCCAGCGGGTTAGCGTTTTTCAACGCATCCGGCAGCAGTGCGTCAgg                          <  2:1836617/150‑1 (MQ=255)
                                                                                                                           cGGCCTCTTCGCGGCTCTAGCCGCCAAGAGGCCGGCACGGGTTGCAGCCCTGAAAGATCACACCGCCTCGCGGCTGCTCACACCCTCCAGCAACCGGGCAATCCCCAGCGGGTTAGCGTTTTTCAACGCATCCGGCAGCAGTGCGTCAgg                          >  1:3395589/1‑150 (MQ=39)
                                                                                                                             gCCTCTTGGCGGCTAGAGCCGCGAAGAGGCCGGCACGGGTTGCAGCCATGAAAGATCACACCGCCTCGCGGCTGCTCACACCCTCCAGCAACCGGGCAATCCCCAGCGGGTTAGCGTTTTTCAACGCATCCGGCAGCAGTGCGTCAGGGt                        <  1:1143611/150‑1 (MQ=255)
                                                                                                                             gCCTCTTGGCGGCTAGAGCCGCGAAGAGGCCGGCACGGGTTGCAGCCATGAAAGATCACACCGCCTCGCGGCTGCTCACACCCTCCAGCAACCGGGCAATCCCCAGCGGGTTAGCGTTTTTCAACGCATCCGGCAGCAGTGCGTCAGGGt                        >  2:5094726/1‑150 (MQ=255)
                                                                                                                               ctctTGGCGGCTAGAGCCGCGAAGAGGCCGGCACGGGTTGCAGCCCTGAAAGATCACACCGCCTCGCGGCTGCGCACAGCCTCCAGCAACCGGGCCCTCCCCCGCGGGTTAGCGTGTTTCAGCGCGTCCGGCGGCAGTGCGGCAGGGtcg                      >  1:4340889/1‑148 (MQ=255)
                                                                                                                               ctctTCGCGGCTCTAGCCGCGAAGAGGCCGGCACGGGTTGCAGCCATGAAAGATCACACCGCCTCGCGGCTGCTCACACCCTCCAGCAACCGGGCAATCCCCAGCGGGTTAGCGTTTTTCAACGCATCCGGCAGCAGTGCGTCAGGGTAg                      >  2:4911662/1‑150 (MQ=255)
                                                                                                                               ctctTCGCGGCTCTAGCCGCCAAGAGGCCGGCACGGGTTGCAGCCATGAAAGATCACACCGCCTCGCGGCTGCTCACACCCTCCAGCAACCGGGCAATCCCCAGCGGGTTAGCGTTTTTCAACGCATCCGGCAGCAGTGCGTCAGGGtcg                      >  2:1969167/1‑148 (MQ=255)
                                                                                                                                tctTGGCGGCTAGAGCCGCGAAGAGGCCGGCACGGGTTGCAGCCATGAAAGATCACACCGCCTCGCGGCTGCTCACACCCTCCAGCAACCGGGCAATCCCCAGCGGGTTAGCGttt                                                       <  2:2375096/116‑1 (MQ=255)
                                                                                                                                tctTGGCGGCTAGAGCCGCGAAGAGGCCGGCACGGGTTGCAGCCATGAAAGATCACACCGCCTCGCGGCTGCTCACACCCTCCAGCAACCGGGCAATCCCCAGCGGGTTAGCGttt                                                       >  1:2375096/1‑116 (MQ=255)
                                                                                                                                tctTGGCGGCTAGAGCCGCGAAGAGGCCGGCACGGGTTGCAGCCATGAAAGATCACACCGCCTCGCGGCTGCTCACACCCTCCAGCAACCGGGCAATCCCCAGCGGGTTAGCGTTTTTCAACGCATCCGGCAGCAGTGCGTCAGGGTCGt                     >  2:2953325/1‑150 (MQ=255)
                                                                                                                                tctTGGCGGCTAGAGCCGCGAAGAGGCCGGCACGGGTTGCAGCCATGAAAGATCACACCGCCTCGCGGCTGCTCACACCCTCCAGCAACCGGGCAATCCCCAGCGGGTTAGCGTTTTTCAACGCATCCGGCAGCAGTGCGTCAGGGTAGt                     >  2:33746/1‑150 (MQ=255)
                                                                                                                                tctTGGCGGCTAGAGCCGCGAAGAGGCCGGCACGGGTTGCAGCCATGAAAGATCACACCGCCTCGCGGCTGCTCACACCCTCCAGCAACCGGGCAATCCCCAGCGGGTTAGCGTTTTTCAACGCATCCGGCAGCAGTGCGTCAGGGTAGt                     >  2:5033907/1‑150 (MQ=255)
                                                                                                                                tctTCGCGGCTCTAGCCGCGAAGAGGCCGGCACGGGTTGCAGTCATGAAAGATCACACCGCCTCGCGGCTGCTCACACCCTCCAGCAACCGGGCAATCCCCAGCGGGTTAGCGTTTTTCAACGCAGCCGGCAGCAGTGCGTCAGGGTAGt                     >  2:1486914/1‑150 (MQ=255)
                                                                                                                                 ctTCGCGGCTCTAGCCGCCAAGAGGCCGGCACGGGTTGCAGCCATGAAAGATCACACCGCCTCGCGGCTGCTCACACCCTCCAGCAACCGGGCAATCCCCAGCGGGTTAGCGTTTTTCa                                                   >  2:2270003/1‑119 (MQ=38)
                                                                                                                                 ctTCGCGGCTCTAGCCGCCAAGAGGCCGGCACGGGTTGCAGCCATGAAAGATCACACCGCCTCGCGGCTGCTCACACCCTCCAGCAACCGGGCAATCCCCAGCGGGTTAGCGTTTTTCa                                                   <  1:2270003/119‑1 (MQ=38)
                                                                                                                                  ttGGCGGCTAGAGCCGCGTAGAGGCCGGCACGGGTTGCAGCCATGAAAGATCACACCTCCTCGCGGCTGCTCACACCCTCCAGCAACCGGGCAATCCCCAGCGGGTTAGCGTTTTTCAACGCATCCGGCAGCAGTGCGTCAGGGTAGTTc                   <  1:3153458/150‑1 (MQ=255)
                                                                                                                                  ttGGCGGCTAGAGCCGCGAAGAGGCCGGCACGGGTTGCAGCCATGAAAGATCACACCGCCTCGCGGCTGCTCACACCCTCCAGCAACCGGGCAATCCCCAGCGGGTTAGCGTTTTTCAACGCATCCGGCAGCAGTGCGTCAGGGTAGTTc                   >  1:795262/1‑150 (MQ=255)
                                                                                                                                  ttGGCGGCTAGAGCCGCGAAGAGGCCGGCACGGGTTGCAGCCATGAAAGATCACACCGCCTCGCGGCTGCTCACACCCTCCAGCAACCGGGCAATCCCCAGCGGGTTAGCGTTTTTCAACGCATCCGGCAGCAGTGCGTCAGGGTAGTTc                   >  1:2338869/1‑150 (MQ=255)
                                                                                                                                     gcggcTAGAGCCGCGAAGAGGCCGGCACGGGTTGCAGCCATGAAAGATCACACCGCCTCGCGGCTGCTCACACCCTCCAGCAACCGGGCAATCCCCAGCGGGTTAGCGTTTTTCAACGCATCCGGCAGCAGTGCGTCAGGGTAGTTCTgg                <  2:2618078/150‑1 (MQ=255)
                                                                                                                                          tAGAGCCGCGAAGAGGCCGGCACGGGTTGCAGCCATGAAAGATCACACCGCCTCGCGGCTGCTCACACCCTCCAGCAACCGGGCAATCCCCAGCGGGTTAGCGTTTTTCAACGCATCCGGCAGCAGTGCGTCAGGGTAGTTCTGGAAGca           >  2:4371617/1‑150 (MQ=255)
                                                                                                                                          gagagCCGCGAAGAGGCCGGCACGGGTTGCAGCCATGAAAGATCACACCGCCTCGCGGCTGCTCACACCCTCCAGCAACCGGGCAATCCCCAGCGGGTTAGCGTTTTTCAACGCATCCGGCAGCAGTGCGTCAGGGTAGTTCTGGAAGca           <  2:3647212/149‑1 (MQ=255)
                                                                                                                                               ccGCGAAGAGGCCGGCACGGGTTGCAGCCATGAAAGATCACACCGCCTCGCGGCTGCTCACACCCTCCAGCAACCGGGCAATCCCCAGCGGGTTAGCGTTTTTCAACGCATCCGGCAGCAGTGCGTCAGGGTAGTTCTGGa               <  2:4821404/141‑1 (MQ=255)
                                                                                                                                               ccGCGAAGAGGCCGGCACGGGTTGCAGCCATGAAAGATCACACCGCCTCGCGGCTGCTCACACCCTCCAGCAACCGGGCAATCCCCAGCGGGTAAGCGTTTTTCAACGCATCCGGCAGCAGTGCGTCAGGGTAGTTCTGGa               >  1:4821404/1‑141 (MQ=255)
                                                                                                                                                   gAAGAGGCCGGCACGGGTTGCAGCCATGAAAGATCACACCGCCTCGCGGCTGCTCACACCCTCCAGCAACCGGGCAATCCCCAGCGGGTTAGCGTTTTTCAACGCATCCGGCAGCAGTGCGTCAGGGTAGTTCTGGAAGCACACCGGgcg  >  2:1354805/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
GCAGGGTTTGTAGATGACATGGTGCAAGCCTGCGCAATACGCATCGGCGATTACCGATGTCATAGGCGCTCGGTATTGTAAGCAAGCTTTTTACTGATTGTATACAGCTCGAAGTGCCCGTGCCGGCCTCTTGGCGGCTAGAGCCGCGAAGAGGCCGGCACGGGTTGCAGCCATGAAAGATCACACCGCCTCGCGGCTGCTCACACCCTCCAGCAACCGGGCAATCCCCAGCGGGTTAGCGTTTTTCAACGCATCCGGCAGCAGTGCGTCAGGGTAGTTCTGGAAGCACACCGGGCG  >  NC_002947/2953167‑2953463

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