Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 85,077 T→G 8.5% intergenic (+55/+75) 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,0770TG 8.5% 696.9 / 39.1 236intergenic (+55/+75)aroE/PP_0075shikimate 5‑dehydrogenase/choline sulfate transporter
Reads supporting (aligned to +/- strand):  ref base T (97/115);  new base G (9/12);  total (106/130)
Fisher's exact test for biased strand distribution p-value = 8.23e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.74e-01

Only 100 of 238 total aligned reads displayed.

AGCAGGCAGCTGAAGCCTTCTTTATCTGGCGTGGCGTGAGGCCGGATACGGCGCCGGTGTTGGCTGAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTACGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGACCCAGAATCCCAACCTCATTCCTCAAACCGGATCGGGCACAACTCCACCCCTTCCAGCTTCTGCAACTCCTCGATCACCTGCGGCCTGGC  >  NC_002947/84928‑85226
                                                                                                                                                     |                                                                                                                                                       
     gcagCTGAAGCCTTCTTTATCTGGCGTGGCGTGAGGCCGGATACGGCGCCGGTGTTGGCTGAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTAcgtc                                                                                                                                                   >  2:4376518/1‑150 (MQ=255)
       agCTGAAGCCTTCTTTATCTGGCGTGGCGTGAGGCCGGATACGGCGCCGGTGTTGGCTGAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTAcgtcgt                                                                                                                                                 >  1:5240521/1‑150 (MQ=255)
       agCTGAAGCCTTCTTTATCTGGCGTGGCGTGAGGCCGGATACGGCGCCGGTGTTGGCTGAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTAcgtcgt                                                                                                                                                 >  2:4839876/1‑150 (MQ=255)
       agCTGAAGCCTTCTTTATCTGGCGTGGCGTGAGGCCGGATACGGCGCCGGTGTTGGCTGAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAACGCCGTTACCACAAAATAactcgt                                                                                                                                                 >  1:936818/1‑150 (MQ=255)
        gCTGAAGCCTTCTTTATCTGGCGTGGCGTGAGGCCGGATACGGCGCCGGTGTTGGCTGAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTACGTCGTa                                                                                                                                                <  1:1771367/150‑1 (MQ=255)
        gCTGAAGCCTTCTTTATCTGGCGTGGCGTGAGGCCGGATACGGCGCCGGTGTTGGCTGAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTACGTCGTa                                                                                                                                                <  1:5623090/150‑1 (MQ=255)
          tGAAGCCTTCTTTATCTGGCGTGGCGTGAGGCCGGATACGGCGCCGGTGTTGGCTGAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAATTAAGTCGTACg                                                                                                                                              >  1:782211/1‑150 (MQ=255)
          tGAAGCCTTCTTTATCTGGCGTGGCGTGAGGCCGGATACGGCGCCGGTGTTGGCTGAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTACCGCGTAcc                                                                                                                                              >  1:4102939/1‑149 (MQ=255)
           gAAGCCTTCTTTATCTGGCGTGGCGTGAGGCCGGATACGGCGCCGGTGTTGGCTGAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTACGTCGTACGc                                                                                                                                             <  1:3391580/150‑1 (MQ=255)
           gAAGCCTTCTTTATCTGGCGTGGCGTGAGGCCGGATACGGCGCCGGTGTTGGCTGAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTACGTCGTACGc                                                                                                                                             >  2:1281170/1‑150 (MQ=255)
                 ttcttTATCTGGCGTGGCGTGAGGCCGGATACGGCGCCGGTGTTGGCTGAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGTCGCAAGGGCGTTGCCACAAAGGACGTCGTACGCCGAGcc                                                                                                                                       >  1:394916/1‑150 (MQ=255)
                 ttcttTATCTGGCGTGGCGTGAGGCCGGATACGGCGCCGGTGTTGGCTGAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCATGGCGGGTCCCAGAAAGTGGGTGGTACGGCGATcg                                                                                                                                       >  1:3580062/1‑149 (MQ=255)
                 ttcttTATCTGGCGTGGCGTGAGGCCGGATACGGCGCCGGTGTTGGCTGAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTACGTCGTACGCCGATcc                                                                                                                                       >  2:2922267/1‑150 (MQ=255)
                    tttATCTGGCGTGGCGTGAGGCCGGATACGGCGCCGGTGTTGGCTGAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTACGTCGTACGCCGATCctct                                                                                                                                    <  2:3499633/150‑1 (MQ=255)
                      tATCTGGCGTGGCGTGAGGCCGGATACGGCGCCGGTGTTGGCTGAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTACGTCGTACGCCGATCCTCtgt                                                                                                                                  <  2:1743425/150‑1 (MQ=255)
                      tATCTGGCGTGGCGTGAGGCCGGATACGGCGCCGGTGTTGGCTGAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTACGTCGTACGCCGATCCTCtgt                                                                                                                                  <  2:3661902/150‑1 (MQ=255)
                      tATCTGGCGTGGCGTGAGGCCGGATACGGCGCCGGTGTTGGCTGAACTGCGCCGGCAACCCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTACGTCGTACGCCGATCCTCtgt                                                                                                                                  <  2:1378391/150‑1 (MQ=255)
                       agcgggcgtggcgaGAGGCCGGATACGGCGCCGGTGTTGGCTGAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTACGTCGTACGCCGATCCTCtgtg                                                                                                                                 <  2:5500879/146‑1 (MQ=255)
                          tggcgtggcgTGAGGCCGGATACGGCGCCGGTGTTGGCTGAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTACGTCGTACGCCGATCCTCTGTGGGa                                                                                                                              <  1:426430/150‑1 (MQ=255)
                            gcgtggcgtgAGGCCGGATACGGCGCCGGTGTTGGCTGAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTACGTCGTACGCCGATCCTCTGTGGGAGc                                                                                                                            <  1:3154165/150‑1 (MQ=255)
                            gcgtggcgtgAGGCCGGATACGGCGCCGGTGTTGGCTGAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGCTCCCACAGAGGATCGGCGTACGACGTACTTTGTGGGaac                                                                                                                            <  1:2285662/150‑3 (MQ=25)
                             cgtggcgtgAGGCCGGATACGGCGCCGGTGTTGGCTGAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTTGCAGCCCATCGCGACGCAAGGCCGCTCCCACAGAGGATCGGCGTACGACGTACTTTGTGGGAACgg                                                                                                                          <  2:273559/150‑1 (MQ=17)
                             cgtggcgtgAGGCCGGATACGGCGCCGGTGTTGGCTGAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTTGCAGCCCATCGCGACGCAAGGCCGCTCCCACAGAGGATCGGCGTACGACGTACTTTGTGGGAACgg                                                                                                                          <  1:3663128/150‑1 (MQ=17)
                               tggcgtgAGGCCGGATACGGCGCCGGTGTTGGCTGAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCCATCGCGACGCAAGGCCGCTCCCACAGAGGATCGGCGTACGACGTACTTTGTGGGAACGGc                                                                                                                         <  1:3101926/150‑1 (MQ=21)
                                 gcgtgAGGCCGGATACGGCGCCGGTGTTGGCTGAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGCTCCCACAGAGGATCGGCGTACGACGTACTTTGTGGGAACGGCCt                                                                                                                       <  2:2173981/150‑1 (MQ=21)
                                  cgtgAGGCCGGATACGGCGCCGGTGTTGGCTGAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTACGTCGTACGCCGATCCTCTGTGGGAGCGGCCtt                                                                                                                      <  2:1436106/150‑1 (MQ=255)
                                   gtgAGGCCGGATACGGCGCCGGTGTTGGCTGAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTACGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTg                                                                                                                     <  2:2436565/150‑1 (MQ=255)
                                     gAGGCCGGATACGGCGCCGGTGTTGGCTGAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTACGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCg                                                                                                                   <  1:456855/150‑1 (MQ=255)
                                      aGGCCGGATACGGCGCCGGTGTTGGCTGAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTACGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGt                                                                                                                  <  1:2880433/150‑1 (MQ=255)
                                      aGGCCGGATACGGCGCCGGTGTTGGCTGAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTACGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGt                                                                                                                  >  1:834637/1‑150 (MQ=255)
                                         ccGGATACGGCGCCGGTGTTGGCTGAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTACGTCTTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTcgc                                                                                                               <  1:2931977/150‑1 (MQ=255)
                                         ccGGATACGGCGCCGGTGTTGGCTGAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTACGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTcgc                                                                                                               <  1:5023263/150‑1 (MQ=255)
                                            gATACGGCGCCGGTGTTGGCTGAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTACGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGAt                                                                                                            <  1:373130/150‑1 (MQ=255)
                                               aCGGCGCCGGTGTTGGCTGACCTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTACGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATggg                                                                                                        <  2:1514575/150‑1 (MQ=255)
                                                 ggCGCCGGTGTTGGCTGAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTACGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCt                                                                                                      >  2:1894553/1‑150 (MQ=255)
                                                  gcgcCGGTGTTGGCTGAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTACGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTg                                                                                                     >  2:3274160/1‑150 (MQ=255)
                                                      cGGTGTTGGCTGAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTACGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCaaa                                                                                                 <  2:1789715/150‑1 (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)
                                                                 gAACTGCGCCGGCAACTCGCTCGGGGCTGCGCGGGCGGGGGCCGCTTGGCGGCCCTACCGCGACGCAAGGCCGTGCCCACAACGTACGTCGTCCGCCGATCCTCTGTGGGCGCGGCCTTGCGCCGCGATGGGCTGCAAAGCAGACCCAgg                                                                                      <  1:3274160/150‑2 (MQ=255)
                                                                  aaCTGCGCCGGCAACGCGCTCGGGGCTGAGCGGGCCGGGGCCGCTGGGCGGCCCTATCGCGCCGCAAGGCCGTTCCCACAAAGTACGTCGTACGCCGATCCACTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGACCCAGaa                                                                                     <  2:2998748/150‑1 (MQ=38)
                                                                      gcgcCGGCAACTCGCTGGGGGCTGTGGTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAATGTACGTCGTACGCCGATCCTGTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGACCCAGAATccc                                                                                 <  1:4547785/150‑1 (MQ=38)
                                                                      gcgcCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTACGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGACCCAGAATccc                                                                                 <  1:3677112/150‑1 (MQ=255)
                                                                       cgcCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTACGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGACCCAGAATCCCa                                                                                >  2:4916003/1‑150 (MQ=255)
                                                                         ccGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCACGGCCGTTCCCACAAAGTACGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGACCCAGAATCCCAAc                                                                              <  2:923933/150‑1 (MQ=255)
                                                                          cGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTACGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGACCCAGAATCCCAAcc                                                                             >  2:650612/1‑150 (MQ=255)
                                                                          cGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTACGTCGTACGCCGATCCTCTGTGGGAGCGGCCGTGCGTCGCGATGGGCTGCACAGCAGACCCAGAAACCCAAcc                                                                             >  2:3063838/1‑150 (MQ=255)
                                                                                ctcgctcgGGGCTGAGCTGGCGGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTACGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGACCCAGAATCCCAACCTCATTc                                                                       <  2:5203492/150‑1 (MQ=255)
                                                                                 tcgctcgGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTACGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGACCCAGAATCCCAACCTCATTcc                                                                      >  2:3230552/1‑150 (MQ=255)
                                                                                  cgctcgGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTACGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGACCCAGAATCCCAACCTCATTCCt                                                                     <  2:3129407/150‑1 (MQ=255)
                                                                                    ctcgGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTACGTCGTACGCCGATCCTCTGTGGGAGCGGCCATGCGTCGCGATGGGCTGCAAAGCAGACCCAGAATCCCAACCTCATTCCTCa                                                                   <  2:1962651/150‑1 (MQ=255)
                                                                                    ctccGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTACGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGACCCAGAATCCCAACCTCATTCCTCa                                                                   >  2:5462715/1‑150 (MQ=255)
                                                                                     tcgGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTACGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGACCCAGAATCCCAACCTCATTCCTCaa                                                                  >  2:3690310/1‑150 (MQ=255)
                                                                                      gggggGTGTGGTGGCGGGGGCCGTTGGGCGGCCCTATCGCGCGGATTGGCCGTTCCCAAAAAGTACGTGGTACGCCGATCCTCTGTGGGAGGGGCTTTGCGTCGCGATGGGCTGCAAAGCAGACCCAGAATCCCAACCTCATTCCTCaaa                                                                 <  1:2916436/149‑1 (MQ=255)
                                                                                         ggCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTACGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGACCCAGAATCCCAACCTCATTCCTCAAACCg                                                              <  2:1026048/150‑1 (MQ=255)
                                                                                          gCTGAGCGGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTACGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGACCCAGAATCCCAACCTCATTCCTCAAACCgg                                                             <  2:3934969/150‑1 (MQ=255)
                                                                                             gAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTACGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGACCCAGAATCCCAACCTCATTCCTCAAACCGGATc                                                          <  1:2273708/150‑1 (MQ=255)
                                                                                               gCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTACGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGTGTCGCGATGGGCTGCAAAGCAGACCCAGAATCCCAACCTCATTCCTCAAACCGGATCgg                                                        <  2:453242/150‑1 (MQ=255)
                                                                                                gtGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCACGGCCGTTCCCACAAAGTACGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGACCCAGAATCCCAACCTCATTCCTCAAACCGGATCggg                                                       <  1:4007980/149‑1 (MQ=255)
                                                                                                cTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCATGGCCGTTCCCACAAAGTACGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGACCCAGAATCCCAACCTCATTCCTCAAACCGGATCggg                                                       <  1:2154936/150‑1 (MQ=255)
                                                                                                 tGGCGGGGGCCGCTTGGCGGCCCTATCGCGACGAAAGGCCGTTCCCACAAAGTACGCCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGACCCAGAATCCCAACCTCATTCCTCAAACCGGATCGGGc                                                      <  1:1768374/150‑1 (MQ=255)
                                                                                                  ggCGGGGGCCGTTTGGCGGCCCTATCGCGCCGAATGGCCGTTCCCAAAACTGACTCCGTCCGCCTATCCTGTGTGGTAGGGGCCTTGCGCCGCGATGGGCTGCAAAGCAGACCCAGAATCCCAACCTCATTCCTCAAACCGGATCGGGca                                                     <  1:101912/150‑1 (MQ=255)
                                                                                                  cgatagggCCGCCAAGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTACGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGACCCAGAATCCCAACCTCATTCCTCAAACCGGATCGGGca                                                     >  2:2361251/6‑150 (MQ=255)
                                                                                                    ccGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTACGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGACCCAGAATCCCAACCTCATTCCGCAAACCGGATCGGGcaca                                                   >  2:402398/1‑150 (MQ=255)
                                                                                                           cGCTTGGCGGCCCTATCGCGACGCAAGGCCATTCCCACAAAGTACGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGACCCAGAATCCCAACCTCATTCCTCAAACCGGATCGGGCACAACTCCAc                                            <  1:4931942/150‑1 (MQ=255)
                                                                                                            gCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTACGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGACCCAGAATCCCAACCTCATTCCTCAAACCGGATCGGGCACAACTCCAcc                                           >  2:3553350/1‑150 (MQ=255)
                                                                                                                 gcggcACTATCGCGACGCAAGGCCGTTCCCACAAAGTACGTCGTACGCCGATCCTCTGTGGGACCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGACCCAGAATCCCAACCACATTCCTCAAACCGGATCGGGCACAACTCCACCCCTTc                                      <  2:387583/150‑1 (MQ=255)
                                                                                                                   ggcCCTATGGCGCCGCAAGGCGGTTCCCACAATGTAGGTCGTACGCCGATCCTGTGTGGGAGCGGCCTTGGGTCGCGATGGGCTGCAAAGCAGACCCAGAATCCCAACCTCATTCCTCAAACCGGATCGGGCACAACTCCACCCCTTCCa                                    <  1:4376518/150‑1 (MQ=38)
                                                                                                                   ggcCCTATCGCGACGCAAGGCCTTTCCCAAAAAGTACGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGACCCAGAATCCCAACCTCATTCCTCAAACCGGATCGGGCACAACTCCACCCCTTCCa                                    <  1:1311020/150‑1 (MQ=255)
                                                                                                                   ggcCCTATCGCGACGCAAGGCCGTTCCCACAAATTACGTCGTACGACGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGACCCAGAATCCCAACCTCATTCCTCAAACCGGATCGGGCACAACTCCACCCCTTCCa                                    <  2:2087535/150‑1 (MQ=255)
                                                                                                                    gcCCTTTCGCCACGCATGGCCGTTCCCAAAACGTACGTGGTACGCCGATCCTCTGTGGGTGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGACCCAGAATCCCAACCTCATTCCTCAAACCGGATCGGGCACAACTCCACCCCTTCCAg                                   <  1:3523966/150‑1 (MQ=39)
                                                                                                                     ccatagcgcgACGCAAGGCCGTTCCAACAACGTCGGCCGTCCGCCGAGCATCGGTGGGTGGGGCCTTGCGTCGCGATGGGCTGCAAAGCAGACCCAGAATCCCAACCTCATTCCTCAAACCGGATCGGGCACAACTCCACCCCTTCCAGc                                  <  1:3323889/144‑1 (MQ=255)
                                                                                                                        taccgctACGCAAGGCCGTTCCCACAACCTACGTCGTACGCCTAACCTCGGGGGGTGGGGCCTTGCGCCGCGATGGGCTGCAAAGCAGACCCAGAATCCCAACCTCATTCCTCAAACCGGATCGGGCACAACTCCACCCCTTCCAGCTTc                               <  1:125787/147‑1 (MQ=255)
                                                                                                                        tATCGCGACGCACGGCCGTTCCCACAAAGTACGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGACCCAGAATCCCAACCTCATTCCTCAAACCGGATCGGGCACAACTCCACCCCTTACAGCTTc                               <  2:5269364/150‑1 (MQ=255)
                                                                                                                         aTCGCGACGCAAGGCCGTTCCCACAAAGTACGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGACCCAGAATACCAACCACATTCCACAAACCGGATCGGGCACAACGCCACCCCATCCAGCATCt                              >  1:2661264/1‑150 (MQ=39)
                                                                                                                           cgcgACGCAAGGCCGTTCCCACAAAGTACGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGACCCAGAATCCCAACCTCATTCCACAAACCGGATCGGGCACAACTCCACCCCTTCCAGCTTCTGc                            >  2:33377/1‑150 (MQ=255)
                                                                                                                            gcgACGCAAGGCCGTTCCCACAAAGTACGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGACCCAGAATCCCAACCTCATTCCTCAAACCGGATCGGGCACAACTCCACCCCTTCCAGCTTCTGCa                           >  2:676361/1‑150 (MQ=255)
                                                                                                                              gACGCAAGGCCGCTCCCACAGAGGATCGGCGTACGACGTACTTTGTGGGAACGGCCTTGCGTCGCGATGGGCTGCAAAGCAGACCCAGAATCCCAACCTCATTCCTCAAACCGGATCGGGCACAACTCCACCCCTTCCAGCTTCTGCAAc                         >  2:1704739/1‑150 (MQ=21)
                                                                                                                               aCGCAAGGCCGTTCCCACAAAGTACGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGACCCAGAATCCCAACCTCATTCCTCAAACCGGATCGGGCACAACTCCACCCCTTCCAGCTTCTGCAAct                        >  1:3465979/1‑150 (MQ=255)
                                                                                                                                 gCAAGGCCGTTCCCACAAAGTACGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGACCCAGAATCCCAACCTCATTCCTCAAACCGGATCGGGCACAACTCCACCCCTTCCAGCTTCTGCAActcc                      >  2:3018139/1‑150 (MQ=255)
                                                                                                                                 gCAAGGCCGTTCCCACAAAGTACGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGACCCAGAATCCCAACCTCATTCCTCAAACCGGATCGGGCACAACTCCACCCCTTCCAGCTTCTGCAActcc                      >  2:4234043/1‑150 (MQ=255)
                                                                                                                                 gCAAGGCCGCTCCCACAGAGGATCGGCGTACGACGTACTTTGTGGGAACGGCCTTGCGTCGCGATGGGCTGCAAAGCAGACCCAGAATCCCAACCTCATTCCTCAAACCGGATCGGGCACAACTCCACCCCTTCCAGCTTCTGCAActcc                      >  1:4426709/1‑150 (MQ=21)
                                                                                                                                    aGGCCGTTCCCACAAAGTACGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGACCCAGAATCCCAACCTCATTCCTCAAACCGGATCGGGCACAACTCCACCCCTTCCAGCTTCTGCAACTCCTCg                   <  2:5418786/150‑1 (MQ=255)
                                                                                                                                    aGGCCGCTCCCACAGAGGATCGGCGTACGACGTACTTTGTGGGAACGGCCTTGCGTCGCGATAGGGCTGCAAAGCAGACCCAGAATCCCAACCTCATTCCTCAAACCGGATCGGGCACAACTCCACCCCTTCCAGCTTCTGCAActcctc                    >  2:4446593/1‑150 (MQ=18)
                                                                                                                                     ggCCGTTCCCACAAAGTACGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGACCCAGAATCCCAACCTCATTCCTCAAACCGGATCGGGCACAACTCCACCCCTTCCAGCTTCTGCAACTCCTCGa                  <  2:1816866/150‑1 (MQ=255)
                                                                                                                                     ggCCGTTCCCACAAAGTACGTCGTACGCCGATCCTCTGTGGGAGCGGAATTGCGTCGCGATTGGCTGCAAAGAAGACCCAGAATACCAACCTCATTCCTCAAACAGGATCGGGCACAACTCCGCACCTTCAAGCTTCTGCAAATCCACGa                  >  2:1861122/1‑150 (MQ=255)
                                                                                                                                       ccGTTCCCACAAAGTACGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGACCCAGAATCCCAACCTCATTCCTCAAACCGGATCGGGCACAACTCCACCCCTTCCAGCTTCTGCAACTCCTCGATc                >  1:809378/1‑150 (MQ=255)
                                                                                                                                         gTTCCCACAAAGTACGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGACCCAGAATCCCAACCTCATTCCTCAAACCGGATCGGGCACAACTCCACCCCTTCCAGCTTCTGCAACTCCTCGATCAc              <  2:3410341/150‑1 (MQ=255)
                                                                                                                                          ttCCCACAAAGTACGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGACCCAGAATCCCAACCTCATTCCTCAAACCGGATCGGGCACAACTCCACCCCTTCCAGCTTCTGCAACTCCTCGATCAcc             <  1:1827381/150‑1 (MQ=255)
                                                                                                                                             ccACAAAGTACTTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGACCCAGAATCCCAACCTCAGTCCTCAAACCGGATCGGGCACAACTCCACCCCTTCCAGCTTCTGCAACTCCTCGATCACCTGc          <  2:4838594/150‑1 (MQ=255)
                                                                                                                                             ccACAAAGTACGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGACCCAGAATCCCAACCTCATTCCTCAAACCGGATCGGGCACAACTCCACCCCTTCCAGCTTCTGCTACTCCTCGATCACCTGc          >  2:4659356/1‑150 (MQ=255)
                                                                                                                                             ccACAAAGTACGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGACCCAGAATCCCAACCTCATTCCTCAAACCGGATCGGGCACAACTCCACCCCTTCCAGCTTCTGCAACTCCTCGATCACCTGc          >  2:3568677/1‑150 (MQ=255)
                                                                                                                                              cacaAAGTACGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGACCCAGAATCCCAACCTCATTCCTCAAACCGGATCGGGCACAACTCCACCCCTTCCAGCTTCTGCAACTCCTCGATCACCTGCg         >  2:2215236/1‑150 (MQ=255)
                                                                                                                                              cacaAAGTACGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGACCCAGAATCCCAACCTCATTCCTCAAACCGGATCGGGCACAACTCCACCCCTTCCAGCTTCTGCAACTCCTCGATCACCTGCg         >  1:2186462/1‑150 (MQ=255)
                                                                                                                                               acaAAGTACGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGACCCAGAATCCCAACCGCATTCCTCAAACCGGATCGGGCACAACTCCACCCCTTCCAGCTTCTGCAACTCCTCGATCACCTGCgg        >  1:5537087/1‑150 (MQ=255)
                                                                                                                                                caAAGTACGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGACCCAGAATCCCAACCTCATTCCTCAAACCGGATCGGGCACAACTCCACCCCTTCCAGCTTCTGCAACTCCTCGATCACCTGCGGc       <  2:2186462/150‑1 (MQ=255)
                                                                                                                                                  aaGTACGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGACCCAGAATCCCAACCTCATTCCTCAAACCGGATCGGGCACAACTCCACCCCTTCCAGCTTCTGCAACTCCTCGATCACCTGCGGCCt     <  2:5620987/150‑1 (MQ=255)
                                                                                                                                                   aGTACGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGACCCAGAATCCCAACCTCATTCCTCAAACCGGATCGGGCACAACTCCACCCCTTCCAGCTTCTGCAACTCCTCGATCACCTGCGGCCTg    <  2:3694628/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                       
AGCAGGCAGCTGAAGCCTTCTTTATCTGGCGTGGCGTGAGGCCGGATACGGCGCCGGTGTTGGCTGAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTACGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGACCCAGAATCCCAACCTCATTCCTCAAACCGGATCGGGCACAACTCCACCCCTTCCAGCTTCTGCAACTCCTCGATCACCTGCGGCCTGGC  >  NC_002947/84928‑85226

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: