Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0099745,105,2670AT76.5% ‑1.0 / 17.1 34G195G (GGT→GGA) RPPX_22650benzene 1,2‑dioxygenase
Reads supporting (aligned to +/- strand):  ref base A (3/5);  new base T (14/12);  total (17/17)
Fisher's exact test for biased strand distribution p-value = 6.88e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: E-value score below prediction cutoff.

GTCTTGATCTCGTCACCCGCCTCGCGCTGCTTGCGCTGGTTCTGGGTGGCAGCGTAGTTCCAGTGCACGGAGCTTACGTGGTAGCCGTCGGCGCCGTTTTCGGCGGTGAGCTTCCAGTTGCCTTCGTAGATGTACGAGCTGGCACCGCGCAGCACTTCCAGGCCTTCAGGCGACTGATCGACGATCATGTCGATGATCTTCGCCGACTCGCCCAGGTGCTCGACCAGTGGTTTCACATCGGCATTCAGACTGCCAAACAGGAAGCCCCGGTACGACTCGAATC  >  CP009974/5105124‑5105406
                                                                                                                                               |                                                                                                                                           
gttttGATTTCGTCGCCTGCCTCGCGCTGCTTGCGCTGGTTCTGAGTGGCGGCGTAGTTCCAGTGCACGGAGCTAACATGATAGCCGTCGGCGCCGTTCTCGGCGGTCAGCTTCCAGTTACCTTCGTAGATGTAGGAGCTAGCTCCACgca                                                                                                                                      <  1:5373807/148‑1 (MQ=0)
 ttttGATTTCGTCGCCTGCCTCGCGCTGCTTGCGCTGGTTCTGAGTGGCGGCGTAGTTCCAGTGCACGGAGCTAACATGATAGCCGTCGGCGCCGTTCTCGGCGGTCAGCTTCCAGTTACCTTCGTAGATGTAGGAGCTAGCTCCACgcag                                                                                                                                     <  1:5957066/149‑1 (MQ=0)
      atttcgtcgCCTGCCTCGCGCTGCTTGCGCTGGTTCTGAGTGGCGGCGTAGTTCCAGTGCACGGAGCTAACATGATAGCCGTCGGCGCCGTTCTCGGCGGTCAGCTTCCAGTTACCTTCGTAGATGTAGGAGCTAGCTCCACGCAGCACtt                                                                                                                                <  2:203728/148‑1 (MQ=1)
      atttcgtcgCCTGCCTCGCGCTGCTTGCGCTGGTTCTGAGTGGCGGCGTAGTTCCAGTGCACGGAGCTAACATGATAGCCGTCGGCGCCGTTCTCGGCGGTCAGCTTCCAGTTACCTTCGTAGATGTAGGAGCTAGCTCCACGCAGCACtt                                                                                                                                <  1:7137341/148‑1 (MQ=0)
                            gcgtGCGCGGGTTCTGGGTGGCAGCGTAGTTGCAGTGCGCGGAGCTTACGTGGTAGCCGTCGGCGCCGTTTTCGGCGGTGAGCTTCCAGTTGCCTTCGTAGATGTACGAGCTGGCACCGCGCAGCACTTCCAGGCCTTCAGGCGACTGATc                                                                                                          <  2:3790997/148‑1 (MQ=0)
                                                gcggCGTAGTTCCAGTGCACGGAGCTAACATGATAGCCGTCGGCGCCGTTCTCGGCGGTCAGCTTCCAGTTACCTTCGTAGATGTAGGAGCTAGCTCCACGCAGCACTTCCAGACCGTCCGGCGACTGGTCGACGATCATGTCGATGATCt                                                                                      >  2:3465787/4‑151 (MQ=1)
                                                gcggCGTAGTTCCAGTGCACGGAGCTAACATGATAGCCGTCGGCGCCGTTCTCGGCGGTCAGCTTCCAGTTACCTTCGTAGATGTAGGAGCTAGCTCCACGCAGCACTTCCAGACCGTCCGGCGACTGGTCGACGATCATGTCGATGATCt                                                                                      >  2:1373076/4‑151 (MQ=1)
                                                     ggaggtCCAGTGCACGGAGCTTACGTGGTAGCCGTCGGCGCCGTTTTCGGCGGTGAGCTTCCAGTTGCCTTCGTAGATGTACGAGCTGGCACCGCGCAGCACTTCCAGGCCTTCAGGCGACTGATCGACGATCATGTCGATGATCTTCGcc                                                                                 <  1:3016523/146‑1 (MQ=1)
                                                             agcgCACGGAGCTAACATGATAGCCGTCGGCGCCGTTCTCGGCGGTCAGCTTCCAGTTACCTTCGTAGATGTAGGAGCTAGCTCCACGCAGCACTTCCAGACCGTCCGGCGACTGGTCGACGATCATGTCGATGATCTTCGCCGACTCGcc                                                                         <  1:3910143/148‑1 (MQ=1)
                                                                          taggTGGTAGCCATCGGCGCCGTTTTCGGCGGTGGGGGTCCACTTGCCTTCGTACATGTACGAGCGGGCACGGCGCAGCACTTCCAGGCCTTCAGGCGACTCATCGACGAGCATCTCGATGATCTTCGCCGACTCGACCCGGTGCACGAcc                                                            >  2:3652225/4‑151 (MQ=1)
                                                                           acatgatAGCCGTCGGCGCCGTTCTCGGCGGTCAGCTTCCAGTTACCTTCGTAGATGTAGGAGCTAGCTCCACGCAGCACTTCCAGACCGTCCGGCGACTGGTCGACGATCATGTCGATGATCTTCGCCGACTCGCCCAGGTGCTCGACCa                                                           <  2:1937371/145‑1 (MQ=1)
                                                                           acatgatAGCCGTCGGCGCCGTTATCGGCGGTCATCTCCCAGTTACCTTCGTATATGTATGCGCTAGCTCCACGCTGCACTTCCAGTCCGTTCGGCGACTGGTCGTCGATCATGTCGATGATCTTCGCCGACTCGCCCAGTTGCTCGACCa                                                           >  2:451137/7‑151 (MQ=0)
                                                                            catgatAGCCGTCGGCGCCGTTCTCGGCGGTCAGCTTCCAGTTACCTTCGTAGATGTAGGAGCTAGCTCCACGCAGCACTTCCAGACCGTCCGGCGACTGGTCGACGATCATGTCGATGATCTTCGCCGACTCGCCCAGGTGCTCGACCAg                                                          >  1:1838746/6‑151 (MQ=1)
                                                                            catgatAGCCGTCGGCGCCGTTCTCGGCGGTCAGCTTCCAGTTACCTTCGTAGATGTAGGAGCTAGCTCCACGCAGCACTTCCAGACCGTCCGGCGACTGGTCGACGATCATGTCGATGATCTTCGCCGACTCGCCCAGGTGCTCGACCAg                                                          <  2:3702759/146‑1 (MQ=1)
                                                                            catgatAGCCGTCGGCGCCGTTCTCGGCGGTCAGCTTCCAGTTACCTTCGTAGATGTAGGAGCTAGCTCCACGCAGCACTTCCAGACCGTCCGGCGACTGGTCGACGATCATGTCGATGATCTTCGCCGACTCGCCCAGGTGCTCGACCAg                                                          >  1:6307603/6‑151 (MQ=1)
                                                                             ttgatAGCCGTCGGCGCCGTTCTCGGCGGTCAGCTTCCAGTTACCTTCGTAGATGTAGGAGCTAGCTCCACGCAGCACTTCCAGACCGTCCGGCGACTGGTCGACGATCATGTCGATGATCTTCGCCGACTCGCCCAGGTGCTCGACCAGt                                                         <  2:4808241/147‑1 (MQ=1)
                                                                             atgatAGCCGTCGGCGCCGTTCTCGGCGGTCAGCTTCCAGTTACCTTCGTAGATGTAGGAGCTAGCTCCACGCAGCACTTCCAGACCGTCCGGCGACTGGTCGACGATCATGTCGATGATCTTCGCCGACTCGCCCAGGTGCTCGACCAGt                                                         >  2:586107/5‑151 (MQ=1)
                                                                             atgatAGCCGTCGGCGCCGTTCTCGGCGGTCAGCTTCCAGTTACCTTCGTAGATGTAGGAGCTAGCTCCACGCAGCACTTCCAGACCGTCCGGCGACTGGTCGACGATCATGTCGATGATCTTCGCCGACTCGCCCAGGTGCTCGACCAGt                                                         >  2:3783025/5‑151 (MQ=1)
                                                                             atgatAGCCGTCGGCGCCGTTCGCGGCGGTCAGCTTCCAGTTACCTTCGTAGATGTAGGAGCTAGCTCCACGCAGCACTGCCAGACCGTCCGGCGACTGGTCGACGATCATGTCGATGATCGTCGCCGACTCGCCCAGGTGCTCGACCAGt                                                         <  2:18649/147‑1 (MQ=1)
                                                                                               gtgttCGGCGGTGAGCTTCCAGTTGCCTTCGTAGATGTACGAGCTGGCACCGCGCAGCACTTCCAGGCCTGCAGGCGACTGATCGACGATCATGTCGATGATCTTCGCCGACTCGCCCAGGTGCTCGACCAGCGGCTTCACATCGGCGTTc                                       <  2:5943429/148‑1 (MQ=1)
                                                                                                ttctCGGCGGTCAGCTTCCAGTTACCTTCGTAGATGTAGGAGCTAGCTCCACGCAGCACTTCCAGACCGTCCGGCGACTGGTCGACGATCATGTCGATGATCTTCGCCGACTCGCCCAGGTGCTCGACCAGTGGTTTCACATCGGCATTCa                                      <  2:489067/148‑1 (MQ=1)
                                                                                                ttctCGGCGGTCAGCTTCCAGTTACCTTCGTAGATGTAGGAGCTAGCTCCACGCAGCACTTCCAGACCGTCCGGCGACTGGTCGACGATCATGTCGATGATCTTCGCCGACTCGCCCAGGTGCTCGACCAGTGGTTTCACATCGGCATTCa                                      >  1:3157872/4‑151 (MQ=1)
                                                                                                 ttgcggcggTGAGCGTCCAGTTGCCTTCGTAGATGTACGAGCTGGCACCGCGCAGCACTTCCAGGCCTTCAGGCGACTGATCGACGATCATGTCGATGATCTTCGCCGACTCGCCCAGGTGCTCGACCAGCGGCTTCACATCGGCGTTCAg                                     <  2:4889129/148‑1 (MQ=1)
                                                                                                         gtcaGCTTCCAGTTACCTTCGTAGATGTAGGAGCTAGCTCCACGCAGCACTTCCAGACCGTCCGGCGACTGGTCGACGATCATGTCGATGATCTTCGCCGACTCGCCCAGGTGCTCGACCAGTGGTTTCACATCGGCATTCAGACTGCCaa                             >  1:5396546/4‑151 (MQ=1)
                                                                                                         gtcaGCTTCCAGTTACCTTCGTAGATGTAGGAGCTAGCTCCACGCAGCACTTCCAGACCGTCCGGCGACTGGTCGACGATCATGTCGATGATCTTCGCCGACTCGCCCAGGTGCTCGACCAGTGGTTTCACATCGGCATTCAGACTGCCaa                             >  1:5394936/4‑151 (MQ=1)
                                                                                                         gtcaGCTTCCAGTTACCTTCGTAGATGTAGGAGCTAGCTCCACGCAGCACTTCCAGACCGTCCGGCGACTGGTCGACGATCATGTCGATCATCTTCGCCGACTCGCCCAGGTGCTCGACCCGTGGTTTCACATCGGCATTCAGACTGCCac                             >  1:4824201/4‑150 (MQ=1)
                                                                                                         gtaaGCTTCCGGTTGCCATCGTTGATGTGCGATCTGGCACCGCGCAGCACCTCCAGGTCTTCCGGCGACTGATCGACGATCATGGAGGTGATCTTAGCCCGCGCGCCCAGGTGCTCGCCCAGCGGCTTCGCATTGTCGTTCAGGCTGccga                             <  2:7202088/148‑3 (MQ=0)
                                                                                                                    gTTGCCTTCGTAGATGTACGAGCTGGCACCGCTCAGCACTTCCAGGCCTTCAGGCGACTGATCGACGATCATGTCGATGATCTTCGCCGACTCGCCCAGGTGCTCGACCAGCGGCTTCACATCGGCGTTCAGGCTGCCGAACAAAAAGccc                  >  1:1081403/1‑151 (MQ=1)
                                                                                                                    gTTGCCTTCGTAGATGTACGAGCTGGCACCGCGCAGCACTTCCAGGCCTTCAGGCGACTGATCGACGATCATGTCGATGATCTTCGCCGACTCGCCCAGGTGCTCGACCAGCGGCTTCACATCGGCGTTCAGGCTGCCGAACAAAAAGccc                  >  1:5602863/1‑151 (MQ=1)
                                                                                                                     ttaccTTCGTAGATGTAGGAGCTAGCTCCACGCAGCACTTCCAGACCGTCCGGCGACTGGTCGACGATCATGTCGATGATCTTCGCCGACTCGCCCAGGTGCTCGACCAGTGGTTTCACATCGGCATTCAGACTGCCAAACAGGAAGcccc                 >  2:3610456/4‑151 (MQ=1)
                                                                                                                     ttaccTTCGTAGATGTAGGAGCTAGCTCCACGCAGCACTTCCAGACCGTCCGGCGACTGGTCGACGATCATGTCGATGATCTTCGCCGACTCGCCCAGGTGCTCGACCAGTGGTTTCACATCGGCATTCAGACTGCCAAACAGGAAGcccc                 <  1:5201937/148‑1 (MQ=1)
                                                                                                                                    taggAGCTAGCTCCACGCAGCACTTCCAGACCGTCCGGCGACTGGTCGACGATCATGTCGATGATCTTCGCCGACTCGCCCAGGTGCTCGACCAGTGGTTTCACATCGGCATTCAGACTGCCAAACAGGAAGCCCCGGTACGACTCGAATc  <  1:4872992/148‑1 (MQ=1)
                                                                                                                                    taggAGCTAGCTCCACGCAGCACTTCCAGACCGTCCGGCGACTGGTCGACGATCATGTCGATGATCTTCGCCGACTCGCCCAGGTGCTCGACCAGTGGTTTCACATCGGCATTCAGACTGCCAAACAGGAAGCCCCGGTACGACTCGAATc  >  2:2604406/4‑151 (MQ=1)
                                                                                                                                    taggAGCTAGCTCCACGCAGCACTTCCAGACCGTCCGGCGACTGGTCGACGATCATGTCGATGATCTTCGCCGACTCGCCCAGGTGCTCGACCAGTGGTTTCACATCGGCATTCAGACTGCCAAACAGGAAGCCCCGGTACGACTCGAATc  >  2:900989/4‑151 (MQ=1)
                                                                                                                                               |                                                                                                                                           
GTCTTGATCTCGTCACCCGCCTCGCGCTGCTTGCGCTGGTTCTGGGTGGCAGCGTAGTTCCAGTGCACGGAGCTTACGTGGTAGCCGTCGGCGCCGTTTTCGGCGGTGAGCTTCCAGTTGCCTTCGTAGATGTACGAGCTGGCACCGCGCAGCACTTCCAGGCCTTCAGGCGACTGATCGACGATCATGTCGATGATCTTCGCCGACTCGCCCAGGTGCTCGACCAGTGGTTTCACATCGGCATTCAGACTGCCAAACAGGAAGCCCCGGTACGACTCGAATC  >  CP009974/5105124‑5105406

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