Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0099745,105,2880TG75.8% ‑0.2 / 11.4 33E188D (GAA→GAC) RPPX_22650benzene 1,2‑dioxygenase
Reads supporting (aligned to +/- strand):  ref base T (3/5);  new base G (16/9);  total (19/14)
Fisher's exact test for biased strand distribution p-value = 2.38e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: E-value score below prediction cutoff.

GCTTGCGCTGGTTCTGGGTGGCAGCGTAGTTCCAGTGCACGGAGCTTACGTGGTAGCCGTCGGCGCCGTTTTCGGCGGTGAGCTTCCAGTTGCCTTCGTAGATGTACGAGCTGGCACCGCGCAGCACTTCCAGGCCTTCAGGCGACTGATCGACGATCATGTCGATGATCTTCGCCGACTCGCCCAGGTGCTCGACCAGTGGTTTCACATCGGCATTCAGACTGCCAAACAGGAAGCCCCGGTACGACTCGAATCGTGCCACCTTG  >  CP009974/5105152‑5105417
                                                                                                                                        |                                                                                                                                 
gcgtGCGCGGGTTCTGGGTGGCAGCGTAGTTGCAGTGCGCGGAGCTTACGTGGTAGCCGTCGGCGCCGTTTTCGGCGGTGAGCTTCCAGTTGCCTTCGTAGATGTACGAGCTGGCACCGCGCAGCACTTCCAGGCCTTCAGGCGACTGATc                                                                                                                     <  2:3790997/148‑1 (MQ=0)
                    gcggCGTAGTTCCAGTGCACGGAGCTAACATGATAGCCGTCGGCGCCGTTCTCGGCGGTCAGCTTCCAGTTACCTTCGTAGATGTAGGAGCTAGCTCCACGCAGCACTTCCAGACCGTCCGGCGACTGGTCGACGATCATGTCGATGATCt                                                                                                 >  2:1373076/4‑151 (MQ=1)
                    gcggCGTAGTTCCAGTGCACGGAGCTAACATGATAGCCGTCGGCGCCGTTCTCGGCGGTCAGCTTCCAGTTACCTTCGTAGATGTAGGAGCTAGCTCCACGCAGCACTTCCAGACCGTCCGGCGACTGGTCGACGATCATGTCGATGATCt                                                                                                 >  2:3465787/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                                                                     >  1:6307603/6‑151 (MQ=1)
                                                catgatAGCCGTCGGCGCCGTTCTCGGCGGTCAGCTTCCAGTTACCTTCGTAGATGTAGGAGCTAGCTCCACGCAGCACTTCCAGACCGTCCGGCGACTGGTCGACGATCATGTCGATGATCTTCGCCGACTCGCCCAGGTGCTCGACCAg                                                                     <  2:3702759/146‑1 (MQ=1)
                                                 ttgatAGCCGTCGGCGCCGTTCTCGGCGGTCAGCTTCCAGTTACCTTCGTAGATGTAGGAGCTAGCTCCACGCAGCACTTCCAGACCGTCCGGCGACTGGTCGACGATCATGTCGATGATCTTCGCCGACTCGCCCAGGTGCTCGACCAGt                                                                    <  2:4808241/147‑1 (MQ=1)
                                                 atgatAGCCGTCGGCGCCGTTCTCGGCGGTCAGCTTCCAGTTACCTTCGTAGATGTAGGAGCTAGCTCCACGCAGCACTTCCAGACCGTCCGGCGACTGGTCGACGATCATGTCGATGATCTTCGCCGACTCGCCCAGGTGCTCGACCAGt                                                                    >  2:3783025/5‑151 (MQ=1)
                                                 atgatAGCCGTCGGCGCCGTTCTCGGCGGTCAGCTTCCAGTTACCTTCGTAGATGTAGGAGCTAGCTCCACGCAGCACTTCCAGACCGTCCGGCGACTGGTCGACGATCATGTCGATGATCTTCGCCGACTCGCCCAGGTGCTCGACCAGt                                                                    >  2:586107/5‑151 (MQ=1)
                                                 atgatAGCCGTCGGCGCCGTTCGCGGCGGTCAGCTTCCAGTTACCTTCGTAGATGTAGGAGCTAGCTCCACGCAGCACTGCCAGACCGTCCGGCGACTGGTCGACGATCATGTCGATGATCGTCGCCGACTCGCCCAGGTGCTCGACCAGt                                                                    <  2:18649/147‑1 (MQ=1)
                                                                   gtgttCGGCGGTGAGCTTCCAGTTGCCTTCGTAGATGTACGAGCTGGCACCGCGCAGCACTTCCAGGCCTGCAGGCGACTGATCGACGATCATGTCGATGATCTTCGCCGACTCGCCCAGGTGCTCGACCAGCGGCTTCACATCGGCGTTc                                                  <  2:5943429/148‑1 (MQ=1)
                                                                    ttctCGGCGGTCAGCTTCCAGTTACCTTCGTAGATGTAGGAGCTAGCTCCACGCAGCACTTCCAGACCGTCCGGCGACTGGTCGACGATCATGTCGATGATCTTCGCCGACTCGCCCAGGTGCTCGACCAGTGGTTTCACATCGGCATTCa                                                 >  1:3157872/4‑151 (MQ=1)
                                                                    ttctCGGCGGTCAGCTTCCAGTTACCTTCGTAGATGTAGGAGCTAGCTCCACGCAGCACTTCCAGACCGTCCGGCGACTGGTCGACGATCATGTCGATGATCTTCGCCGACTCGCCCAGGTGCTCGACCAGTGGTTTCACATCGGCATTCa                                                 <  2:489067/148‑1 (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             >  2:900989/4‑151 (MQ=1)
                                                                                                        taggAGCTAGCTCCACGCAGCACTTCCAGACCGTCCGGCGACTGGTCGACGATCATGTCGATGATCTTCGCCGACTCGCCCAGGTGCTCGACCAGTGGTTTCACATCGGCATTCAGACTGCCAAACAGGAAGCCCCGGTACGACTCGAATc             >  2:2604406/4‑151 (MQ=1)
                                                                                                        taggAGCTAGCTCCACGCAGCACTTCCAGACCGTCCGGCGACTGGTCGACGATCATGTCGATGATCTTCGCCGACTCGCCCAGGTGCTCGACCAGTGGTTTCACATCGGCATTCAGACTGCCAAACAGGAAGCCCCGGTACGACTCGAATc             <  1:4872992/148‑1 (MQ=1)
                                                                                                                  ctccacgcAGCACTTCCAGACCGTCCGGCGACTGGTCGACGATCATGTCGATGATCTTCGCCGACTCGCCCAGGTGCTCGACCAGTGGTTTCACATCGGCATTCAGACTGCCAAACAGGAAGCCCCGGTACGACTCGAATCGTGCCACCtt   <  2:2173939/146‑1 (MQ=1)
                                                                                                                  ctccacgcAGCACTTCCAGACCGTCCGGCGACTGGTCGACGATCATGTCGATGATCTTCGCCGACTCGCCCAGGTGCTCGACCAGTGGTTTCACATCGGCATTCAGACTGCCAAACAGGAAGCCCCGGTACGACTCGAATCGGGCCACCtt   >  1:1509157/6‑151 (MQ=1)
                                                                                                                   tccacgcAGCACTTCCAGACCGTCCGGCGACTGGTCGACGATCATGTCGATGATCTTCGCCGACTCGCCCAGGTGCTCGACCAGTGGTTTCACATCGGCATTCAGACTGCCAAACAGGAAGCCCCGGTACGACTCGAATCGTGCCACCTTg  >  1:589807/5‑151 (MQ=1)
                                                                                                                                        |                                                                                                                                 
GCTTGCGCTGGTTCTGGGTGGCAGCGTAGTTCCAGTGCACGGAGCTTACGTGGTAGCCGTCGGCGCCGTTTTCGGCGGTGAGCTTCCAGTTGCCTTCGTAGATGTACGAGCTGGCACCGCGCAGCACTTCCAGGCCTTCAGGCGACTGATCGACGATCATGTCGATGATCTTCGCCGACTCGCCCAGGTGCTCGACCAGTGGTTTCACATCGGCATTCAGACTGCCAAACAGGAAGCCCCGGTACGACTCGAATCGTGCCACCTTG  >  CP009974/5105152‑5105417

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: ‑