Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0099745,105,2580CG78.8% ‑0.3 / 14.3 33S198S (TCG→TCC) RPPX_22650benzene 1,2‑dioxygenase
Reads supporting (aligned to +/- strand):  ref base C (2/5);  new base G (17/9);  total (19/14)
Fisher's exact test for biased strand distribution p-value = 1.06e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: E-value score below prediction cutoff.

GCCAGCACTCATGGTCTTGATCTCGTCACCCGCCTCGCGCTGCTTGCGCTGGTTCTGGGTGGCAGCGTAGTTCCAGTGCACGGAGCTTACGTGGTAGCCGTCGGCGCCGTTTTCGGCGGTGAGCTTCCAGTTGCCTTCGTAGATGTACGAGCTGGCACCGCGCAGCACTTCCAGGCCTTCAGGCGACTGATCGACGATCATGTCGATGATCTTCGCCGACTCGCCCAGGTGCTCGACCAGTGGTTTCACATCGGCATTCAGACTGCCAAACAGGAAGCCCCGGTA  >  CP009974/5105111‑5105395
                                                                                                                                                   |                                                                                                                                         
gccggcgcTCATGGTTTTGATTTCGTCGCCTGCCTCGCGCTGCTTGCGCTGGTTCTGAGTGGCGGCGTAGTTCCAGTGCACGGAGCTAACATGATAGCCGTCGGCGCCGTTCTCGGCGGTCAGCTTCCAGTTACCTTCGTAGATGTAGGAg                                                                                                                                        >  1:1627764/8‑151 (MQ=1)
  cggcgcTCATGGTTTTGATTTCGTCGCCTGCCTCGCGCTGCTTGCGCTGGTTCTGAGTGGCGGCGTAGTTCCAGTGCACGGAGCTAACATGATAGCCGTCGGCGCCGTTCTCGGCGGTCAGCTTCCAGTTACCTTCGTAGATGTAGGAGCt                                                                                                                                      >  2:80121/6‑151 (MQ=1)
  cggcgcTCATGGTTTTGATTTCGTCGCCTGCCTCGCGCTGCTTGCGCTGGTTCTGAGTGGCGGCGTAGTTCCAGTGCACGGAGCTAACATGATAGCCGTCGGCGCCGTTCTCGGCGGTCAGCTTCCAGTTACCTTCGTAGATGTAGGAGCt                                                                                                                                      >  2:3684421/6‑151 (MQ=1)
  cggcgcTCATGGTTTTGATTTCGTCGCCTGCCTCGCGCTGCTTGCGCTGGTTCTGAGTGGCGGCGTAGTTCCAGTGCACGGAGCTAACATGATAGCCGTCGGCGCCGTTCTCGGCGGTCAGCTTCCAGTTACCTTCGTAGATGTAGGAGCt                                                                                                                                      >  1:4076414/6‑151 (MQ=1)
  cggcgcTCATGGTTTTGATTTCGTCGCCTGCCTCGCGCTGCTTGCGCTGGTTCTGAGTGGCGGCGTAGTTCCAGTGCACGGAGCTAACATGATAGCCGTCGGCGCCGTTCTCGGCGGTCAGCTTCCAGTTACCTTCGTAGATGTAGGAGCt                                                                                                                                      <  2:1530013/146‑1 (MQ=1)
              ttttGATTTCGTCGCCTGCCTCGCGCTGCTTGCGCTGGTTCTGAGTGGCGGCGTAGTTCCAGTGCACGGAGCTAACATGATAGCCGTCGGCGCCGTTCTCGGCGGTCAGCTTCCAGTTACCTTCGTAGATGTAGGAGCTAGCTCCACgcag                                                                                                                          <  2:2052055/149‑1 (MQ=1)
              ttttGATTTCGTCGCCTGCCTCGCGCTGCTTGCGCTGGTTCTGAGTGGCGGCGTAGTTCCAGTGCACGGAGCTAACATGATAGCCGTCGGCGCCGTTCTCGGCGGTCAGCTTCCAGTTACCTTCGTAGATGTAGGAGCTAGCTCCACgcag                                                                                                                          >  1:435692/3‑151 (MQ=0)
              ttttGATTTCGTCGCCTGCCTCGCGCTGCTTGCGCTGGTTCTGAGTGGCGGCGTAGTTCCAGTGCACGGAGCTAACATGATAGCCGTCGGCGCCGTTCTCGGCGGTCAGCTTCCAGTTACCTTCGTAGATGTAGGAGCTAGCTCCACgcag                                                                                                                          <  1:4645983/149‑1 (MQ=0)
                   atttcgtcgCCTGCCTCGCGGTGCTTGCGCTGGTTCTGAGTGGCGGCGTAGTTCCAGTGCACGGAGCTAACATGATAGCCGTCGGCGCCGTTCTCGGCGGTCAGCTTCCAGTTACCTTCGTAGATGTAGGAGCTAGCTCCACGCAGCACtt                                                                                                                     <  2:1633285/148‑1 (MQ=1)
                   atttcgtcgCCTGCCTCGCGCTGCTTGCGCTGGTTCTGAGTGGCGGCGTAGTTCCAGTGCACGGAGCTAACATGATAGCCGTCGGCGCCGTTCTCGGCGGTCAGCTTCCAGTTACCTTCGTAGATGTAGGAGCTAGCTCCACGCAGGACtt                                                                                                                     >  1:4410420/4‑151 (MQ=1)
                   atttcgtcgCCTGCCTCGCGCTGCTTGCGCTGGTTCTGAGTGGCGGCGTAGTTCCAGTGCACGGAGCTAACATGATAGCCGTCGGCGCCGTTCTCGGCGGTCAGCTTCCAGTTACCTTCGTAGATGTAGGAGCTAGCTCCACGCAGCACtt                                                                                                                     >  2:4904446/4‑151 (MQ=1)
                   atttcgtcgCCTGCCTCGCGCTGCTTGCGCTGGTTCTGAGTGGCGGCGTAGTTCCAGTGCACGGAGCTAACATGATAGCCGTCGGCGCCGTTCTCGGCGGTCAGCTTCCAGTTACCTTCGTAGATGTAGGAGCTAGCTCCACGCAGCACtt                                                                                                                     >  2:269696/4‑151 (MQ=1)
                   atttcgtcgCCTGCCTCGCGCTGCTTGCGCTGGTTCTGAGTGGCGGCGTAGTTCCAGTGCACGGAGCTAACATGATAGCCGTCGGCGCCGTTCTCGGCGGTCAGCTTCCAGTTACCTTCGTAGATGTAGGAGCTAGCTCCACGCAGCACtt                                                                                                                     >  1:3239041/4‑151 (MQ=0)
                     ctcgcgaCCCGCCCCGCGCTGCGTGCGCCGGCTCTGGGTGGCAGCGTAGTTCCAGTGCACGGGGCTTACGGGGTCGCCGTCGGCGCCGTTTTCGGCGGTGAGCTGCCAGTTGCCTTCGTAGATGTACGAGCTGGCACCGCGCAGCACTTcc                                                                                                                   <  2:3322828/145‑1 (MQ=1)
                                                             gcggCGTAGTTCCAGTGCACGGAGCTAACATGATAGCCGTCGGCGCCGTTCTCGGCGGTCAGCTTCCAGTTACCTTCGTAGATGTAGGAGCTAGCTCCACGCAGCACTTCCAGACCGTCCGGCGACTGGTCGACGATCATGTCGATGATCt                                                                           >  1:139127/4‑151 (MQ=1)
                                                                                 ggggCTTACGGGGTAGCCGTCGGCGCCGTTTTCGGCGGTGAGCTTCCAGTTGCCTTCGTAGATGTACGAGCTGGCACCGCGCAGCACTTCCAGGCCTTCAGGCGACTGATCGACGATCATGTCGATGATCTTCGCCGACTCGCCCAGGTGc                                                       <  2:3680425/148‑1 (MQ=1)
                                                                                       aacatgatAGCCGTCGGCGCCGTTCTCGGCGGTCAGCTTCCAGTTACCTTCGTAGATGTAGGAGCTAGCTCCACGCAGCACTTCCAGACCGTCCGGCGACTGGTCGACGATCATGTCGATGATCTTCGCCGACTCGCCCAGGTGCTCGAcc                                                 <  2:4147227/144‑1 (MQ=1)
                                                                                          atgatAGCGGTCGGCGCCGTTCTCGGCGGTCAGCTTCCAGTTACCTTCGTAGATGTAGGAGCTAGCTCCACGCAGCACTTCCAGACCGTCCGGCGACTGGTCGACGATCATGTCGATGATCTTCGCCGACTCGCCCAGGTGCTCGACCAGt                                              <  2:667921/147‑1 (MQ=1)
                                                                                          atgatAGCCGTCGGCGCCGTTCTCGGCGGTCAGCTTCCAGTTACCTTCGTAGATGTAGGAGCTAGCTCCACGCAGCACTTCCAGACCGTCCGGCGACTGGTCGACGATCATGTCGATGATCTTCGCCGACTCGCCCAGGTGCTCGACCAGt                                              >  2:13944/5‑151 (MQ=1)
                                                                                          atgatAGCCGTCGGCGCCGTTCTCGGCGGTCAGCTTCCAGTTACCTTCGTAGATGTAGGAGCTAGCTCCACGCAGCACTTCCAGACCGTCCGGCGACTGGTCGACGATCATGTCGATGATCTTCGCCGACTCGCCCAGGTGCTCGACCAGt                                              <  1:2636810/147‑1 (MQ=1)
                                                                                           tgatAGCCGTCGGCGCCGTTCTCGGCGGTCAGCTTCCAGTTACCTTCGTAGATGTAGGAGCTAGCTCCACGCAGCACTTCCAGACCGTCCGGCGACTGGTCGACGATCATGTCGATGATCTTCGCCGACTCGCCCAGGTGCTCGACCAGTg                                             >  2:1432336/4‑151 (MQ=1)
                                                                                                 ccGTCGGCGCCGTTTTCGGCGGTGAGCTTCCAGTTGCCTTCGTAGATGTACGAGCTGGCACCGCGCAGCACTTCCAGGCCTTCAGGCGACTGATCGACGATCATGTCGATGATCTTCGCCGACTCGCCCAGGTGCTCGACCAGCGGCTTca                                       >  1:5075057/1‑151 (MQ=0)
                                                                                                 ccGTCGGCGCCGTTTTCGGCGGTGAGCTTCCAGTTGCCTTCGTAGATGTACGAGCTGGCACCGCGCAGCACTTCCAGGCCTTCAGGCGACTGATCGACGATCATGTCGATGATCTTCGCCGACTCGCCCAGGTGCTCGACCAGCGGCTTca                                       >  1:4544586/1‑151 (MQ=0)
                                                                                                             ttctCGGCGATCAGCTTCCAGTTACCTTCGTAGATGTAGGAGCTAGCTCCACGCAGCACTTCCAGACCGTCCGGCGACTGGTCGACGATCATGTCGATGATCTTCGCCGACTCGCCCAGGTGCTCGACCAGTGGTTTCACATCGGCATTCa                           >  2:2585471/4‑151 (MQ=1)
                                                                                                                      gtcaGCTTCCAGTTACCTTCGTAGATGTAGGAGCTAGCTCCACGCAGCACTTCCAGACCGTCCGGCGACTGGTCGACGATCATGTCGATGATCTTCGCCGACTCGCCCAGGTGCTCGACCAGTGGTTTCACATCGGCATTCAGACTGCCaa                  >  2:1220759/4‑151 (MQ=1)
                                                                                                                      gtcaGCTTCCAGTTACCTTCGTAGATGTAGGAGCTAGCTCCACGCAGCACTTCCAGACCGTCCGGCGACTGGTCGACGATCATGTCGATGATCTTCGCCGACTCGCCCAGGTGCTCGACCAGTGGTTTCACATCGGCATTCAGACTGCCaa                  >  1:150966/4‑151 (MQ=1)
                                                                                                                      gtcaGCTTCCAGTTACCTTCGTAGATGTAGGAGCTAGCTCCACGCAGCACTTCCAGACCGTCCGGCGACTGGTCGACGATCATGTCGATGATCTTCGCCGACTCGCCCAGGTGCTCGACCAGTGGTTTCACATCGGCATTCAGACTGCCaa                  >  1:3043802/4‑151 (MQ=1)
                                                                                                                      gtcaGCTTCCAGTTACCTTCGTAGATGTAGGAGCTAGCTCCACGCAGCACTTCCAGACCGTCCGGCGACTGGTCGACGATCATGTCGATGATCTTCGCCGACTCGCCCAGGTGCTCGACCAGTGGTTTCACATCGGCATTCAGACTGCCaa                  >  1:3043435/4‑151 (MQ=1)
                                                                                                                                 gTTGCCTTCGTAGATGTACGAGCTGGCACCGCGCAGCACTTCCAGGCCTTCAGGCGACTGATCGACGATCATGTCGATGATCTTCGCCGACTCGCCCAGGTGCTCGACCAGCGGCTTCACATCGGCGTTCAGGCTGCCGAACAAAAAGccc       <  1:884125/151‑1 (MQ=1)
                                                                                                                                 gTTGCCTTCGTAGATGTACGAGCTGGCACCGCGCAGCACTTCCAGGCCTTCAGGCGACTGATCGACGATCATGTCGATGATCTTCGCCGACTCGCCCAGGTGCTCGACCAGCGGCTTCACATCGGCGTTCAGGCTGCCGAACAAAAAGccc       <  1:3049541/151‑1 (MQ=1)
                                                                                                                                  ttaccTTCGTAGATGTAGGAGCTAGCTCCACGCAGCACTTCCAGACCGTCCGGCGACTGGTCGACGATCATGTCGATGATCTTCGCCGACTCGCCCAGGTGCTCGACCAGTGGTTTCACATCGGCATTCAGACTGCCAAACAGGAAGcccc      <  1:3706044/148‑1 (MQ=1)
                                                                                                                                  ttaccTTCGTAGATGTAGGAGCTAGCTCCACGCAGCACTTCCAGACCGTCCGGCGACTGGTCGACGATCATGTCGATGATCTTCGCCGACTCGCCCAGGTGCTCGACCAGTGGTTTCACATCGGCATTCAGACTGCCAAACAGGAAGcccc      <  1:1110405/148‑1 (MQ=1)
                                                                                                                                      ctgcGTAGGTGTACGAGCTGGCACCGCGCAGCACTTCCAGGCCTTCAGGCGACTGATCGACGATCATGTCGATGATCTTCGCCGACTCGCCCAGGTGCTCGACCAGCGGCTTCACATCGGCGTTCAGGCTGCAGAACAAAAAGCCCCGGTa  <  2:2022719/148‑1 (MQ=1)
                                                                                                                                                   |                                                                                                                                         
GCCAGCACTCATGGTCTTGATCTCGTCACCCGCCTCGCGCTGCTTGCGCTGGTTCTGGGTGGCAGCGTAGTTCCAGTGCACGGAGCTTACGTGGTAGCCGTCGGCGCCGTTTTCGGCGGTGAGCTTCCAGTTGCCTTCGTAGATGTACGAGCTGGCACCGCGCAGCACTTCCAGGCCTTCAGGCGACTGATCGACGATCATGTCGATGATCTTCGCCGACTCGCCCAGGTGCTCGACCAGTGGTTTCACATCGGCATTCAGACTGCCAAACAGGAAGCCCCGGTA  >  CP009974/5105111‑5105395

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