Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0099742,738,4170AG76.0% ‑3.6 / 11.2 25E210E (GAA→GAGRPPX_12165benzene 1,2‑dioxygenase
Reads supporting (aligned to +/- strand):  ref base A (3/3);  new base G (9/10);  total (12/13)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.19e-01
Rejected as consensus: E-value score below prediction cutoff.

AATGCCGATGTGAAACCACTGGTCGAGCACCTGGGCGAGTCGGCGAAGATCATCGACATGATCGTCGATCAGTCGCCTGAAGGCCTGGAAGTGCTGCGCGGTGCCAGCTCGTACATCTACGAAGGCAACTGGAAGCTCACCGCCGAAAACGGCGCCGACGGCTACCACGTAAGCTCCGTGCACTGGAACTACGCTGCCACCCAGAACCAGCGCAAGCAGCGCGAGGCGGGTGACGAGATCAAGACCATGAGTGCTGGCGCCTGGGCCAAGC  >  CP009974/2738271‑2738541
                                                                                                                                                  |                                                                                                                            
aacgCCGATGTGAAGCCGCTGGTCGAGCACCTGGGCGAGTCGGCGAAGATCATCGACATGATCGTCGATCAGTCGCCTGAAGGCCTGGAAGTGCTGCGCGGTGCCAGCTCGTACATCTACGAAGGCAACTGGAAGCTCACCGCCGAAAACg                                                                                                                          >  1:4766890/4‑151 (MQ=0)
           gacgccgcTGGTCGAGCACCTGGGCGAGTCGGCGAACATCATCGACATGATCGTCGATCAGTCGCCTGAAGGCCTGGAAGTGCTGCGCGGTGCCAGCTCGTACATCTACGAAGGCAACTGGAAGCTCACCGCCGAAAACGGCGCCGACGGc                                                                                                               <  1:4194278/144‑1 (MQ=0)
            aagccgcTGGTCGAGCACCTGGGCGAGTCGGCGAAGATCATCGACATGATCGTCGATCAGTCGCCTGAAGGCCTGGAAGTGCTGCGCGGTGCCAGCTCGTACATCTACGAAGGCAACTGGAAGCTCACCGCCGAAAACGGCGCCGACGGCt                                                                                                              <  1:1639199/145‑1 (MQ=1)
            aagccgcTGGTCGAGCACCTGGGCGAGTCGGCGAAGATCATCGACATGATCGTCGATCAGTCGCCTGAAGGCCTGGAAGTGCTGCGCGGTGCCAGCTCGTACATCTACGAAGGCAACTGGAAGCTCACCGCCGAAAACGGCGCCGACGGCt                                                                                                              <  1:3092175/145‑1 (MQ=1)
              gccgcTGGTCGAGCACCTGGGCGAGTCGGCGAAGATCATCGACATGATCGTCGATCAGTCGCCTGAAGGCCTGGAAGTGCTGCGCGGTGCCAGCTCGTACATCTACGAAGGCAACTGGAAGCTCACCGCCGAAAACGGCGCCGACGGCTAc                                                                                                            >  2:2087755/5‑151 (MQ=1)
               ccgcTGGTCGAGCACCTGGGCGAGTCGGCGAAGATCATCGACATGATCGTCGATCAGTCGCCTGAAGGCCTGGAAGTGCTGCGCGGTGCCAGCTCGTACATCTACGAAGGCAACTGGAAGCTCACCGCCGAAAACGGCGCCGACGGCTAcc                                                                                                           >  1:3999139/4‑151 (MQ=1)
                         aGCACCTGGGCGAGTCGGCGAAGATCATCGACATGATCGTCGACCAGTCGCCGGACGGTCTGGAAGTGCTGCGTGGAGCTAGCTCCTACATCTACGAAGGTAACTGGAAGCTGACCGCCGAGAACGGCGCCGACGGCTATCATGTTAGCTc                                                                                                 >  1:1682676/1‑151 (MQ=1)
                                                                  gaccAGTCGCCGGACGGTCTGGAAGTGCTGCGTGGAGCTAGCTCCTACATCTACGAAGGTAACTGGAAGCTGACCGCCGAGAACGGCGCCGACGGCTATCATGTTAGCTCCGTGCACTGGAACTACGCCGCCACTCAGAACCAGCGCAagc                                                        >  2:4361186/4‑151 (MQ=1)
                                                                  gaccAGTCGCCGGACGGTCTGGAAGTGCTGCGTGGAGCTAGCTCCTACATCTACGAAGGTAACTGGAAGCTGACCGCCGAGAACGGCGCCGACGGCTATCATGTTAGCTCCGTGCACTGGAACTACGCCGCCACTCAGAACCAGCGCAagc                                                        >  1:2476636/4‑151 (MQ=0)
                                                                  gaccAGTCGCCGGACGGTCTGGAAGTGCTGCGTGGAGCTAGCTCCTACATCTACGAAGGTAACTGGAAGCTGACCGCCGAGAACGGCGCCGACGGCTATCATGTTAGCTCCGTGCACTGGAACTACGCCGCCACTCAGAACCAGCGCAagc                                                        <  1:3838177/148‑1 (MQ=1)
                                                                  gaccAGTCGCCGGACGGTCTGGAAGTGCTGCGTGGAGCTAGCTCCTACATCTACGAAGGTAACTGGAAGCTGACCGCCGAGAACGGCGCCGACGGCTATCATGTTAGCTCCGTGCACTGGAACTACGCCGCCACTCAGAACCAGCGCAagc                                                        >  1:870881/4‑151 (MQ=0)
                                                                  gaccAGTCGCCGGACGGTCTGGAAGTGCTGCGTGGAGCTAGCTCCTACATCTACGAAGGTAACTGGAAGCTGACCGCCGAGAACGGCGCCGACGGCTATCATGTTAGCTCCGTGCACTGGAACTACGCCGCCACTCAGAACCAGCGCAagc                                                        >  1:1466668/4‑151 (MQ=0)
                                                                                    cTGGAAGTGCTGCGTGGAGCTAGCTCCTACATCTACGAAGGTAACTGGAAGCTGACCGCCGAGAACGGCGCCGACGGCTATCATGTTAGCTCCGTGCACTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACg                                      >  1:4531903/1‑151 (MQ=1)
                                                                                                  tggagctaGCTCCTACATCTACGAAGGTAACTGGAAGCTGACCGCCGAGAACGGCGCCGACGGCTATCATGTTAGCTCCGTGCACTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATg                        <  1:2224635/144‑1 (MQ=1)
                                                                                                   ggagctaGCTCCTACATCTACGAAGGTAACTGGAAGCTGACCGCCGAGAACGGCGCCGACGGCTATCATGTTAGCTCCGTGCACTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGa                       >  1:4090766/7‑151 (MQ=1)
                                                                                                     agctaGCTCCTACGTCTACGAAGGTAACTGCAAGCTGACCGCCGAGAACGGCGCCGACGGCTATCATGTTAGCTCCGTGCACTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAgc                     <  2:2485100/147‑2 (MQ=1)
                                                                                                     agctaGCTCCTACATCTACGAAGGTAACTGGAAGCTGACCGCCGAGAACGGCGCCGACGGCTATCATGTTAGCTCCGTGCACTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAgc                     <  1:3911754/147‑2 (MQ=1)
                                                                                                     agctaGCTCCTACATCTACGAAGGTAACTGGAAGCTGACCGCCGAGAACGGCGCCGACGGCTATCATGTTAGCTCCGTGCACTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAgc                     <  1:928391/147‑2 (MQ=1)
                                                                                                     agctaGCTCCTACATCTACGAAGGTAACTGGAAGCTGACCGCCGAGAACGGCGCCGACGGCTATCATGTTAGCTCCGTGCACTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAgc                     >  1:1887559/5‑150 (MQ=1)
                                                                                                      gctaGCTCCTACATCTACGAAGGTAACTGGAAGCTGACCGCCGAGAACGGCGCCGACGGCTATCATGTTAGCTCCGTGCACTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAgcg                    >  1:4934264/4‑149 (MQ=1)
                                                                                                      gctaGCTCCTACATCTACGAAGGTAACTGGAAGCTGACCGCCGAGAACGGCGCCGACGGCTATCATGTTAGCTCCGTGCACTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAgcg                    <  1:172489/148‑3 (MQ=1)
                                                                                                      gctaGCTCCTACAACTACGAAGGTAACTGGAAGCTGACCGCCGAGAACGGCGCCGACGGCTATCATGTTAGCTCCGTGCACTGAAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAgcg                    <  1:5227194/148‑3 (MQ=1)
                                                                                                            tcctACATCTACGAAGGTAACTGGAAGCTGACCGCCGAGAACGGCGCCGACGGCTATCATGTTAGCTCCGTGCACTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGgct              <  1:4932243/148‑2 (MQ=0)
                                                                                                                       ggagggtaaCGGGAAGCTGACCGCCGAGAACGGGGCCGACGGCTATCATGTTAGCTCCGTGCACTGGAACTACGCCGCCACTCAGAACCAGCGCCAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAg   <  2:3519538/144‑1 (MQ=0)
                                                                                                                        gaatgtaaCTGGAAGCTGACCGCCGAGAACGGCGCCGACGGCTATCATGTTAGCTCCGTGCACTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGc  <  2:5068417/145‑1 (MQ=0)
                                                                                                                                                  |                                                                                                                            
AATGCCGATGTGAAACCACTGGTCGAGCACCTGGGCGAGTCGGCGAAGATCATCGACATGATCGTCGATCAGTCGCCTGAAGGCCTGGAAGTGCTGCGCGGTGCCAGCTCGTACATCTACGAAGGCAACTGGAAGCTCACCGCCGAAAACGGCGCCGACGGCTACCACGTAAGCTCCGTGCACTGGAACTACGCTGCCACCCAGAACCAGCGCAAGCAGCGCGAGGCGGGTGACGAGATCAAGACCATGAGTGCTGGCGCCTGGGCCAAGC  >  CP009974/2738271‑2738541

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: