Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 85,093 C→A 7.2% intergenic (+71/+59) aroE → / ← PP_0075 shikimate 5‑dehydrogenase/choline sulfate transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00294785,0930CA 7.2% 747.8 / 33.0 232intergenic (+71/+59)aroE/PP_0075shikimate 5‑dehydrogenase/choline sulfate transporter
Reads supporting (aligned to +/- strand):  ref base C (85/130);  new base A (8/9);  total (93/139)
Fisher's exact test for biased strand distribution p-value = 6.11e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.87e-01

Only 100 of 242 total aligned reads displayed.

CTTCTTTATCTGGCGTGGCGTGAGGCCGGATACGGCGCCGGTGTTGGCTGAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTACGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGACCCAGAATCCCAACCTCATTCCTCAAACCGGATCGGGCACAACTCCACCCCTTCCAGCTTCTGCAACTCCTCGATCACCTGCGGCCTGGCCCGGTGCAAGGTCAAA  >  NC_002947/84944‑85242
                                                                                                                                                      |                                                                                                                                                     
 ttcttTATCTGGCGTGGCGTGAGGCCGGATACGGCGCCGGTGTTGGCTGAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGTCGCAAGGGCGTTGCCACAAAGGACGTCGTACGCCGAGcc                                                                                                                                                      >  1:394916/1‑150 (MQ=255)
 ttcttTATCTGGCGTGGCGTGAGGCCGGATACGGCGCCGGTGTTGGCTGAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCATGGCGGGTCCCAGAAAGTGGGTGGTACGGCGATcg                                                                                                                                                      >  1:3580062/1‑149 (MQ=255)
   cttTATCTGGCGTGGCGTGAGGCCGGATACGGCGCCGGTGTTGGCTGAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTACGTCGTACGCCGATCctc                                                                                                                                                    <  1:3215936/150‑1 (MQ=255)
    tttATCTGGCGTGGCGTGAGGCCGGATACGGCGCCGGTGTTGGCTGAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTACGTCGTACGCCGATCctct                                                                                                                                                   <  2:3499633/150‑1 (MQ=255)
      tATCTGGCGTGGCGTGAGGCCGGATACGGCGCCGGTGTTGGCTGAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTACGTCGTACGCCGATCCTCtgt                                                                                                                                                 <  2:5591667/150‑1 (MQ=255)
      tATCTGGCGTGGCGTGAGGCCGGATACGGCGCCGGTGTTGGCTGAACTGCGCCGGCAACCCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTACGTCGTACGCCGATCCTCtgt                                                                                                                                                 <  2:1378391/150‑1 (MQ=255)
      gaTCTGGCGTGGCGTGAGGCCGGATACGGCGCCGGTGTTGGCTGAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTACGTCGTACGCCGATCCTCtgt                                                                                                                                                 <  2:3662046/149‑1 (MQ=255)
       agcgggcgtggcgaGAGGCCGGATACGGCGCCGGTGTTGGCTGAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTACGTCGTACGCCGATCCTCtgtg                                                                                                                                                <  2:5500879/146‑1 (MQ=255)
            gcgtggcgtgAGGCCGGATACGGCGCCGGTGTTGGCTGAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTACGTCGTACGCCGATCCTCTGTGGGAGc                                                                                                                                           <  1:3154165/150‑1 (MQ=255)
            gcgtggcgtgAGGCCGGATACGGCGCCGGTGTTGGCTGAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAACTACGACGTACGCCCAACCTCAGTGGGAGc                                                                                                                                           >  1:4838594/1‑150 (MQ=39)
            gcgtggcgtgAGGCCGGATACGGCGCCGGTGTTGGCTGAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGCTCCCACAGAGGATCGGCGTACGACGTACTTTGTGGGaac                                                                                                                                           <  1:2285662/150‑3 (MQ=25)
            gcgtggcgtgAGGCCGGATACGGCGCCGGTGTTGGCTGAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCCATCGCGACGCAAGGCCGCTCCCACAGAGGATCGGCGTACGACGTACTTTGTGGGaac                                                                                                                                           <  1:84540/150‑3 (MQ=21)
             cgtggcgtgAGGCCGGATACGGCGCCGGTGTTGGCTGAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTTGCAGCCCATCGCGACGCAAGGCCGCTCCCACAGAGGATCGGCGTACGACGTACTTTGTGGGAACgg                                                                                                                                         <  1:3663264/150‑1 (MQ=17)
             cgtggcgtgAGGCCGGATACGGCGCCGGTGTTGGCTGAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTTGCAGCCCATCGCGACGCAAGGCCGCTCCCACAGAGGATCGGCGTACGACGTACTTTGTGGGAACgg                                                                                                                                         <  2:273559/150‑1 (MQ=17)
               tggcgtgAGGCCGGATACGGCGCCGGTGTTGGCTGAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTACGTCGTACGCCGATCCTCTGTGGGAGCGGc                                                                                                                                        <  1:1799386/150‑1 (MQ=255)
               tggcgtgAGGCCGGATACGGCGCCGGTGTTGGCTGAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCCATCGCGACGCAAGGCCGCTCCCACAGAGGATCGGCGTACGACGTACTTTGTGGGAACGGc                                                                                                                                        <  1:3101926/150‑1 (MQ=21)
                  cgtgAGGCCGGATACGGCGCCGGTGTTGGCTGAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTACGTCGTACGCCGATCCTCTGTGGGAGCGGCCtt                                                                                                                                     <  2:1436106/150‑1 (MQ=255)
                   gtgAGGCCGGATACGGCGCCGGTGTTGGCTGAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTACGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTg                                                                                                                                    <  2:2436565/150‑1 (MQ=255)
                     gAGGCCGGATACGGAGCCGGTGGTGGCTGAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCTGCCATATCGCGACGCAAGTCCGTTCACACAAAGGACGTCGTCCGCCGATCCTCGGTGGGAGCGGCCTTGCg                                                                                                                                  <  1:456745/150‑1 (MQ=255)
                      aGGCCGGATACGGCGCCGGTGTTGGCTGAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTACGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGt                                                                                                                                 <  1:2880433/150‑1 (MQ=255)
                            gATACGGCGCCGGTGTTGGCTGAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTACGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGAt                                                                                                                           <  1:372800/150‑1 (MQ=255)
                             aTACGGCGCCGGTGTTGGCTGAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTACGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATg                                                                                                                          >  2:927300/1‑150 (MQ=255)
                               aCGGCGCCGGTGTTGGCTGACCTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTACGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATggg                                                                                                                        <  2:1514575/150‑1 (MQ=255)
                                 ggCGCCGGTGTTGGCTGAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTACGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCt                                                                                                                      >  2:1785258/1‑150 (MQ=255)
                                 ggCGCCGGTGTTGGCTGAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTACGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCt                                                                                                                      >  2:1894553/1‑150 (MQ=255)
                                  gcgcCGGTGTTGGCTGAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTACGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTg                                                                                                                     >  2:3274160/1‑150 (MQ=255)
                                    gcCGGTGTTGGCTGAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTACGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCa                                                                                                                   >  2:5226684/1‑150 (MQ=255)
                                        gtgtTGGCTGACCTGCGCCGGAACCGCGCGCGGGGCTGCGCTGGCCGGGGCCGCTTGGCGGCCCTACCGCGACGCAAGGCCGTTCCCACAATGTACGTCGTACGCCGAGCCTCTGTGGGCGCGGCCTTGCGTCGCGATGGGCTGCAAAGc                                                                                                               <  2:8355/150‑1 (MQ=37)
                                        gtgtGGGCTTAACTGCGCCGGCAACTCGCGCGGGGCTGAGCTGGCGGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTCCCCAAAAAGTACGTCTTACGCCGACCCTCGGTGGGAGCGGCCTTGCGCCGCGATGGGCTGCAAAGc                                                                                                               <  1:377190/150‑1 (MQ=255)
                                         tgtTGGCTTAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCGGGGGCCGATTGGCGGCCCTATCGCGACGCACGGCCGTTCCCACAAAGTACGTCGTACGCCGATCCTCGGGGGGAGCGGCCTTGCGACGCGATGGGCTGCAAAGCa                                                                                                              <  1:848646/150‑1 (MQ=255)
                                             ggCTGAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGGCGGTCCCACAAAATAAGTAGGAACCCCATACAAAGTTGGTGCGGCCGTTCGTCGCGATGGGCTGCaaacaagaac                                                                                                          >  1:4696009/1‑143 (MQ=19)
                                                tGAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTACGGCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGACCCAg                                                                                                       >  2:2999862/1‑150 (MQ=255)
                                                  aaCTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGGTCCCACAAAGTACGTCGTACGCCGATCCTCTGTGGGAGCGGCCGTGCGTCGCGATGGGCTGCAAAGCAGACCCAGaa                                                                                                     >  1:2172389/1‑150 (MQ=255)
                                                  aaCTGCGCCGGCAACGCGCTCGGGGCTGAGCGGGCCGGGGCCGCTGGGCGGCCCTATCGCGCCGCAAGGCCGTTCCCACAAAGTACGTCGTACGCCGATCCACTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGACCCAGaa                                                                                                     <  2:2998748/150‑1 (MQ=38)
                                                    cggcgcCGGCACCTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTACGACGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGACGCGATGGGCTGCAAAGCAGACCCAGAATc                                                                                                   <  2:3557567/148‑1 (MQ=255)
                                                    cTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTACGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGTGTCGCGATGGGCTGCAAAGCAGACCCAGAATc                                                                                                   <  1:4525653/150‑1 (MQ=255)
                                                       cgcCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTACGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGACCCAGAATCCCa                                                                                                >  2:4916003/1‑150 (MQ=255)
                                                          cGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTACGTCGTACGCCGATCCTCTGTGGGAGCGGCCGTGCGTCGCGATGGGCTGCACAGCAGACCCAGAAACCCAAcc                                                                                             >  2:3063838/1‑150 (MQ=255)
                                                                ctcgctcgGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTACGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGACCCAGAATCCCAACCTCATTc                                                                                       <  2:2759814/150‑1 (MQ=255)
                                                                ctcccacgGGGCTCAGCTGCCCGGGCCCGCTTGCCGCCCCAAACGCCACCCAAGCCCGTACCCACAAACAACCACCAACCCCCAACCCCTGTGGCACCCCCCTACCGACGCGATGGGCTGCAAAGCAGACCCAGAATCCCAACCTCATTc                                                                                       <  2:1899842/144‑1 (MQ=255)
                                                                 tcgctcgGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTACGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGACCCAGAATCCCAACCTCATTcc                                                                                      >  2:3230548/1‑150 (MQ=255)
                                                                  cgctcgGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTACGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGACCCAGAATCCCAACCTCATTCCt                                                                                     <  2:3129407/150‑1 (MQ=255)
                                                                    ctccGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTACGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGACCCAGAATCCCAACCTCATTCCTCa                                                                                   >  2:5462715/1‑150 (MQ=255)
                                                                     tcgGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTACGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGACCCAGAATCCCAACCTCATTCCTCaa                                                                                  >  2:3690310/1‑150 (MQ=255)
                                                                         ggCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTACGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGACCCAGAATCCCAACCTCATTCCTCAAACCg                                                                              <  2:1026048/150‑1 (MQ=255)
                                                                          ggtgcggtGGCGGGGGCCGCGTGGCGCCCCTAGCGCGACGAACGGCCGTTCCCACACAGTAGGTCGTACGCCGATCCTCTGGGGGAGGGGCCTTGCGTCGCGATGGGCTGCAAAGCAGACCCAGAATCCCAACCTCATTCCTCAAACCgg                                                                             <  1:3390184/143‑1 (MQ=255)
                                                                                gtGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCACGGCCGTTCCCACAAAGTACGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGACCCAGAATCCCAACCTCATTCCTCAAACCGGATCggg                                                                       <  1:4007980/149‑1 (MQ=255)
                                                                                 tGGCGGGGGCCGCTTGGCGGCCCTATCGCGACGAAAGGCCGTTCCCACAAAGTACGCCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGACCCAGAATCCCAACCTCATTCCTCAAACCGGATCGGGc                                                                      <  1:1768374/150‑1 (MQ=255)
                                                                                  ggCGGGGGCCGTTTGGCGGCCCTATCGCGCCGAATGGCCGTTCCCAAAACTGACTCCGTCCGCCTATCCTGTGTGGTAGGGGCCTTGCGCCGCGATGGGCTGCAAAGCAGACCCAGAATCCCAACCTCATTCCTCAAACCGGATCGGGca                                                                     <  1:101912/150‑1 (MQ=255)
                                                                                  ggCGGGGGCCGCTTGGCGGCCCTATCGGCACTCACGGCGGTCCCCCAAACGTACGCCGAACGCCGACCCGCTGTGGGCGGGGCCTTGCGCCGCGATGGGCTGCAAAGCAGACCCAGAATCCCAACCTCATTCCTCAAACCGGATCGGGca                                                                     <  1:675711/150‑1 (MQ=255)
                                                                                  cgatagggCCGCCAAGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTACGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGACCCAGAATCCCAACCTCATTCCTCAAACCGGATCGGGca                                                                     >  2:2361251/6‑150 (MQ=255)
                                                                                   gatagggCCGCCAAGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTACGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGACCCAGAATCCCAACCTCATTCCTCAAACCGGATCGGGcac                                                                    >  2:377940/5‑150 (MQ=255)
                                                                                    ccGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTACGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGACCCAGAATCCCAACCTCATTCCGCAAACCGGATCGGGcaca                                                                   >  2:402398/1‑150 (MQ=255)
                                                                                           cGCTTGGCGGCCCTATCGCGACGCAAGGCCATTCCCACAAAGTACGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGACCCAGAATCCCAACCTCATTCCTCAAACCGGATCGGGCACAACTCCAc                                                            <  1:2100527/150‑1 (MQ=255)
                                                                                            gCTTGGCGGCCCTATCGCGACGCAACGCCGTTCCCACAAAGTACGTCGAACGCCGATCCTCAGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGACCAAGAATCCCAACCTCATTCCTCAAACAGGATCGGGCACAACTCCAcc                                                           >  2:3881036/1‑150 (MQ=255)
                                                                                                ggcggcCCTATCGCGACGCAAGGCCGTTCCCACAAAGTACGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGACCCAGAATCCCAACCTCATTCCTCAAACCGGATCGGGCACAACTCCACCCCtt                                                       <  2:799234/150‑1 (MQ=255)
                                                                                                   ggcCCTATCGCGACGCAAGGCCTTTCCCAAAAAGTACGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGACCCAGAATCCCAACCTCATTCCTCAAACCGGATCGGGCACAACTCCACCCCTTCCa                                                    <  1:1311020/150‑1 (MQ=255)
                                                                                                   ggcCCTATCGCGACGCAAGGCCGTTCCCACAAATTACGTCGTACGACGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGACCCAGAATCCCAACCTCATTCCTCAAACCGGATCGGGCACAACTCCACCCCTTCCa                                                    <  2:2087535/150‑1 (MQ=255)
                                                                                                    gcCCTTTCGCCACGCATGGCCGTTCCCAAAACGTACGTGGTACGCCGATCCTCTGTGGGTGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGACCCAGAATCCCAACCTCATTCCTCAAACCGGATCGGGCACAACTCCACCCCTTCCAg                                                   <  1:3523966/150‑1 (MQ=39)
                                                                                                    gcCCTATCGCGACGCATGGCCTTTCCCACAAAGTACGTCGTCCGCCGATCCTCTGTGGGAGCGGCGTTGCGTCGCGATGGGCTGCAAAGCAGACCCAGAATCCCAACCTCATTCCTCAAACCGGATCGGGCACAACTCCACCCCTTCCAg                                                   <  1:2922267/150‑1 (MQ=255)
                                                                                                     ccatagcgcgACGCAAGGCCGTTCCAACAACGTCGGCCGTCCGCCGAGCATCGGTGGGTGGGGCCTTGCGTCGCGATGGGCTGCAAAGCAGACCCAGAATCCCAACCTCATTCCTCAAACCGGATCGGGCACAACTCCACCCCTTCCAGc                                                  <  1:3323889/144‑1 (MQ=255)
                                                                                                        taccgctACGCAAGGCCGTTCCCACAACCTACGTCGTACGCCTAACCTCGGGGGGTGGGGCCTTGCGCCGCGATGGGCTGCAAAGCAGACCCAGAATCCCAACCTCATTCCTCAAACCGGATCGGGCACAACTCCACCCCTTCCAGCTTc                                               <  1:125787/147‑1 (MQ=255)
                                                                                                        tATCGCGACGCAAGGCCGTTCCCACAAAGTACGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGACCCAGAATCCCAACCTCATTCCTCAAACCGGATCGGGCACAACTCCACCCCTTCCAGCTTc                                               <  1:5131825/150‑1 (MQ=255)
                                                                                                        tATCGCGACGCAAGGCCGTTCCCACAAAGTACGTCGTACGCCCATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGACCCAGAATCCCAACCTCATTCCTCAAACCGGATCGGGCACAACTCCACCCCTTCCAGCtta                                               >  1:432109/1‑149 (MQ=255)
                                                                                                           cgcgACGCAAGGCCGTTCCCACAAAGTACGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGACCCAGAATCCCAACCTCATTCCACAAACCGGATCGGGCACAACTCCACCCCTTCCAGCTTCTGc                                            >  2:33377/1‑150 (MQ=255)
                                                                                                            gcgACGCAAGGCCGTTCCCACAAAGTACGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGACCCAGAATCCCAACCTCATTCCTCAAACCGGATCGGGCACAACTCCACCCCTTCCAGCTTCTGCa                                           >  2:1821150/1‑150 (MQ=255)
                                                                                                            gcgACGCAAGGCCGTTCCCACAAAGTACGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGACCCAGAATCCCAACCTCATTCCTCAAACCGGATCGGGCACAACTCCACCCCTTCCAGCTTCTGCa                                           >  2:676361/1‑150 (MQ=255)
                                                                                                             cgACGCAAGGCCGTTCCCACAAAGTACGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGACCCAGAATCCCAACCTCATTCCTCAAACCGGATCGGGCACAACTCCACCCCTTCCAGCTTCTGCaa                                          >  2:4437988/1‑150 (MQ=255)
                                                                                                              gACGCAAGGCCGTTCCCACAAAGTAGGTCGTACGCCGATCCTCGGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGACCCAGAATCCCAACCTCATTCCTCAAACCGGATCGGGCACAATTCCACCCCTTCCAGCTTCTGCAAc                                         <  2:1972073/150‑1 (MQ=255)
                                                                                                              gACGCAAGGCCGCTCCCACAGAGGATCGGCGTACGACGTACTTTGTGGGAACGGCCTTGCGTCGCGATGGGCTGCAAAGCAGACCCAGAATCCCAACCTCATTCCTCAAACCGGATCGGGCACAACTCCACCCCTTCCAGCTTCTGCAAc                                         >  2:1704782/1‑150 (MQ=21)
                                                                                                               aCGCAAGGCCGTTCCCACAAAGTACGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGACCCAGAATCCCAACCTCATTCCTCAAACCGGATCGGGCACAACTCCACCCCTTCCAGCTTCTGCAAct                                        >  1:3465979/1‑150 (MQ=255)
                                                                                                                 gCACGGCCGTGCCCACAACTCCCGTCGTACGCCCACCCTCTGTGGGAGCGGCCTCCCGACGCGATGGGCTGCAAAGCAGACCCAGAATCCCAACCTCATTCCTCAACCCGGATCGGGCACAACTCCACCCCTTCCAGCTTCTGCAActcc                                      <  1:2215212/150‑1 (MQ=255)
                                                                                                                 gCAAGGCCGTTCCCACAAAGTACGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGACCCAGAATCCCAACCTCATTCCTCAAACCGGATCGGGCACAACTCCACCCCTTCCAGCTTCTGCAActcc                                      >  2:4234043/1‑150 (MQ=255)
                                                                                                                 gCAAGGCCGCTCCCACAGAGGATCGGCGTACGACGTACTTTGTGGGAACGGCCTTGCGTCGCGATGGGCTGCAAAGCAGACCCAGAATCCCAACCTCATTCCTCAAACCGGATCGGGCACAACTCCACCCCTTCCAGCTTCTGCAActcc                                      >  1:4426507/1‑150 (MQ=21)
                                                                                                                    aGGCCGTTCCCACAAAGTACGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGACCCAGAATCCCAACCTCATTCCTCAAACCGGATCGGGCACAACTCCACCCCTTCCAGCTTCTGCAACTCCTCg                                   <  2:5418786/150‑1 (MQ=255)
                                                                                                                     ggCCGTTCCCACAAAGTACGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGACCCAGAATCCCAACCTCATTCCTCAAACCGGATCGGGCACAACTCCACCCCTTCCAGCTTCTGCAACTCCTCGa                                  <  2:1816857/150‑1 (MQ=255)
                                                                                                                      gCCGTTCCCACAAAGTACGTCGTACGCCGATCCGCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGACCCAGAATCCCAACCTCAGTCCTCAAACCGGATCGGGCACAACTCCACCCCTTCCAGCTTCTGCAACTCCTCGag                                 >  2:3600431/1‑149 (MQ=255)
                                                                                                                       ccGTTCCCACAAAGTACGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGACCCAGAATCCCAACCTCATTCCTCAAACCGGATCGGGCACAACTCCACCCCTTCCAGCTTCTGCAACTCCTCGATc                                >  1:406986/1‑150 (MQ=255)
                                                                                                                       ccGTTCCCACAAAGTACGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGACCCAGAATCCCAACCTCATTCCTCAAACCGGATCGGGCACAACTCCACCCCTTCCAGCTTCTGCAACTCCTCGATc                                >  1:809378/1‑150 (MQ=255)
                                                                                                                         gTTCCCACAAAGTACGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGACCCAGAATCCCAACCTCATTCCTCAAACCGGATCGGGCACAACTCCACCCCTTCCAGCTTCTGCAACTCCTCGATCAc                              <  2:3410341/150‑1 (MQ=255)
                                                                                                                           tCCCACAGAGGATCGGCGTACGACGTACTTTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGACCCAGAATCCCAACCTCATTCCTCAAACCGGATCGGGCACAACTCCACCCCTTCCAGCTTCTGCAACTCCTCGATCACCt                            >  1:5008888/1‑150 (MQ=25)
                                                                                                                             ccACAAAGTACGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGACCCAGAATCCCAACCTCATTCCTCAAACCGGATCGGGCACAACTCCACCCCTTCCAGCTTCTGCTACTCCTCGATCACCTGc                          >  2:4659356/1‑150 (MQ=255)
                                                                                                                              cacaAAGTACGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGACCCAGAATCCCAACCTCATTCCTCAATCCGGATCGGGCACAACTCCACCCCTTCCAGCTTCTGCAACTCCTCGATCACCTGCg                         >  2:590037/1‑150 (MQ=255)
                                                                                                                               acaAAGTACGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGTAAAGCAGACCCAGAATCCCAACCTCATTCCTCAAACCGGATCGGGCACAACTCCACCCCTTCCAGCTTCTGCACCTCCTCGATCACCTGCgg                        >  2:4495600/1‑150 (MQ=255)
                                                                                                                               acaAAGTACGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGACCCAGAATCCCAACCTCATTCCTCAAACCGGATCGGGCACAACTCCACCCCTTCCAGCTTCTGCAACTCCTCGATCACCTGCgg                        >  2:4768808/1‑150 (MQ=255)
                                                                                                                                caAAGTACGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGACCCAGAATCCCAACCTCATTCCTCAAACCGGATCGGGCACAACTCCACCCCTTCCAGCTTCTGCAACTCCTCGATCACCTGCGGc                       <  2:2186462/150‑1 (MQ=255)
                                                                                                                                        cgtcgtACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGACCCAGAATACCAACCTCATTCCTCAAACCGGATCGGGCACAACTCCACCCCTTCCAGCTTCTGCAACTCCTCGATCACCTGCGGCCTGGccc                >  1:662858/1‑150 (MQ=255)
                                                                                                                                         gtcgtACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGACCCAGCATCCCAACCTCATTCCTCAAACCGGATCGGGCACAACTCCACCCCTTCCAGCTTCTGCAACTCCTCGATCACCTGCGGCCTGGCCCg               >  1:4238210/1‑150 (MQ=255)
                                                                                                                                           cgtACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGACCCAGAATCCCAACCTCATTCCTCAAACCGGATCGGGCACAACTCCACCCCTTCCAGCTTCTGCAACTCCTCGATCACCTGCGGCCTGGCCCGGt             <  2:5043326/150‑1 (MQ=255)
                                                                                                                                               cGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGACCCAGAATCCCAACCTCATTCCTCAAACCGGATCGGGCACAACTCCACCCCTTCCAGCTTCTGCAACTCCTCGATCACCTGCGGCCTGGCCCGGTGCaa         <  2:5143621/150‑1 (MQ=255)
                                                                                                                                                gCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGACCCAGAATCCCAACCTCATTCCTCAAACCGGATCGGGCACAACTCCACCCCTTCCAGCTTCTGCAACTCCTCGATCACCTGCGGCCTGGCCCGGTGCAAg        <  2:936818/150‑1 (MQ=255)
                                                                                                                                                gCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGACCCAGAATCCCAACCTCATTCCTCAAACCGGATCGGGCACAACTCCACCCCTTCCAGCTTCTGCAACTCCTCGATCACCTGCGGCCTGGCCCGGTGCAAg        >  2:2147414/1‑150 (MQ=255)
                                                                                                                                                gCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGACCCAGAATCCCAACCTCATTCCTCAAACCGGATCGGGCACAACTCCACCCCTTCCAGCTTCTGCAACTCCTCGATCACCTGCGGCCTGGCCCGGTGCAAg        >  2:4977338/1‑150 (MQ=255)
                                                                                                                                                  cGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGACCCAGAATCCCAACCTCATTCCTCAAACCGGATCGGGCACAACTCCACCCCTTCCAGCTTCTGCAACTCCTCGATCACCTGCGGCCTGGCCCGGTGCAAGGt      >  2:5295257/1‑150 (MQ=255)
                                                                                                                                                  cGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGACCCAGAATCCCAACCTCATTCCTCAAACCGGATCGGGCACAACTCCACCCCTTCCAGCTTCTGCAACTCCTCGATCACCTGCGGCCTGGCCCGGTGCAAGGt      >  2:1537333/1‑150 (MQ=255)
                                                                                                                                                    aTCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGACCCAGAATCCCAACCTCATTCCTCAAACCGGATCGGGCACAACTCCACCCCTTCCAGCTTCTGCAACTCCTCGATCACCTGCGGCCTGGCCCGGTGCAAGGTCa    >  2:4537847/1‑150 (MQ=255)
                                                                                                                                                    aTCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGACCCAGAATCCCAACCTCATTCCTCAAACCGGATCGGGCACAACTCCACCCCTTCCAGCTTCTGCAACTCCTCGATCACCTGCGGCCTGGCCCGGTGCAAGGTCa    <  2:3701253/150‑1 (MQ=255)
                                                                                                                                                    aTCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGACCCAGAATCCCAACCTCATTCCTCAAACCGGATCGGGCACAACTCCACCCCTTCCAGCTTCTGCAACTCCTCGATCACCTGCGGCCTGGCCCGGTGCAAGGTCa    <  2:3449792/150‑1 (MQ=255)
                                                                                                                                                     tCCTCTGTGGGAGGGGCCTTGCGTCGCGATGGGCTGCAAAGCAGACCCAGAATCCCAACCTCATTCCTCAAACCGGATCGGGCACAACTCCACCCCTTCCAGCTTCTGCAACTCCTCGATCACCTGCGGCCTGGCCCGGTGCAAGGTCaa   <  1:246993/150‑1 (MQ=255)
                                                                                                                                                      ccTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGACCCAGAATCCCAACCTCATTCCTCAAACCGGATCGGGCACAACTCCACCCCTTCCAGCTTCTGCAACTCCTCGATCACCTGCGGCCTGGCCCGGTGCAAGGTCaaa  <  1:5379126/150‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                                     
CTTCTTTATCTGGCGTGGCGTGAGGCCGGATACGGCGCCGGTGTTGGCTGAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTACGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGACCCAGAATCCCAACCTCATTCCTCAAACCGGATCGGGCACAACTCCACCCCTTCCAGCTTCTGCAACTCCTCGATCACCTGCGGCCTGGCCCGGTGCAAGGTCAAA  >  NC_002947/84944‑85242

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