Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | CP009974 | 5,105,300 | 0 | A | G | 70.0% | 1.6 / 34.1 | 50 | D184D (GAT→GAC) | RPPX_22650 | benzene 1,2‑dioxygenase |
Reads supporting (aligned to +/- strand): ref base A (6/9); new base G (21/14); total (27/23) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 2.28e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
Rejected as consensus: E-value score below prediction cutoff. | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
GTTCTGGGTGGCAGCGTAGTTCCAGTGCACGGAGCTTACGTGGTAGCCGTCGGCGCCGTTTTCGGCGGTGAGCTTCCAGTTGCCTTCGTAGATGTACGAGCTGGCACCGCGCAGCACTTCCAGGCCTTCAGGCGACTGATCGACGATCATGTCGATGATCTTCGCCGACTCGCCCAGGTGCTCGACCAGTGGTTTCACATCGGCATTCAGACTGCCAAACAGGAAGCCCCGGTACGACTCGAATCGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTG > CP009974/5105162‑5105450 | gtgcTGAGTGGCGGCGTAGTTCCGGTGCGCGGAGCTAACATGATAGCCGTCGGGGCCGTTCTCGGCGGTCAGCTTCCAGTTGCCTTCGTAGATGTAGGAGCTAGCTCCACGCAGCACTTCCAGACCGTCCGGCGACTGGTCGACGATCATg < 2:3657898/148‑1 (MQ=0) tgagTGGCGGCGTAGTTCCAGTGCACGGAGCTAACATGATAGCCGTCGGCGCCGTTCTCGGCGGTCAGCTTCCAGTTACCTTCGTAGATGTAGGAGCTAGCTCCACGCAGCACTTCCAGACCGTCCGGCGACTGGTCGACGATCATGTCga > 1:175116/4‑151 (MQ=0) tgagTGGCGGCGTAGTTCCAGTGCACGGAGCTAACATGATAGCCGTCGGCGCCGTTCTCGGCGGTCAGCTTCCAGTTACCTTCGTAGATGTAGGAGCTAGCTCCACGCAGCACTTCCAGACCGTCCGGCGACTGGTCGACGATCATGTCga < 2:6378079/148‑1 (MQ=1) tgagTGGCGGCGTAGTTCCAGTGCACGGAGCTAACATGATAGCCGTCGGCGCCGTTCTCGGCGGTCAGCTTCCAGTTACCTTCGTAGATGTAGGAGCTAGCTCCACGCAGCACTTCCAGACCGTCCGGCGACTGGTCGACGATCATGTCga > 2:4350729/4‑151 (MQ=0) tgagTGGCGGCGTAGTTCCAGTGCACGGAGCTAACATGATAGCCGTCGGCGCCGTTCTCGGCGGTCAGCTTCCAGTTACCTTCGTAGATGTAGGAGCTAGCTCCACGCAGCACTTCCAGACCGTCCGGCGACTGGTCGACGATCATGTCga > 2:4350700/4‑151 (MQ=0) gcggCGTAGTTCCAGTGCACGGAGCTAACATGATAGCCGTCGGCGCCGTTCTCGGCGGTCAGCTTCCAGTTACCTTCGTAGATGTAGGAGCTAGCTCCACGCAGCACTTCCAGACCGTCCGGCGACTGGTCGACGATCATGTCGATGATCt > 2:4171102/4‑151 (MQ=0) gcggCGTAGTTCCAGTGCACGGAGCTAACATGATAGCCGTCGGCGCCGTTCTCGGCGGTCAGCTTCCAGTTACCTTCGTAGATGTAGGAGCTAGCTCCACGCAGCACTTCCAGACCGTCCGGCGACTGGTCGACGATCATGTCGATGATCt > 2:1714197/4‑151 (MQ=1) tcgaGTGCACGGAGCTAACATGATAGCCGTCGGCGCCGTTCTCGGCGGTCAGCTTCCAGTTACCTTCGTAGATGTAGGAGCTAGCTCCACGCAGCACTTCCAGACCGTCCGGCGACTGGTCGACGATCATGTCGATGATCTTCGCCGACTc < 2:6065821/148‑1 (MQ=1) aacctgatAGCCGTCGGCGCCGTTCTCGGCGGTCAGCTTCCAGTTACCTTCGTAGATGTAGGAGCTAGCTCCACGCAGCACTTCCAGACCGTCCGGCGACTGGTCGACGATCATGTCGATGATCTTCGCCGACTCGCCCAGGTGCTCGAcc < 1:3634356/144‑1 (MQ=1) aacatgatAGCCGTCGGCGCCGTTCTCGGCGGTCAGCTTCCAGTTACCTTCGTAGATGTAGGAGCTAGCTCCACGCAGCACTTCCAGACCGTCCGGCGACTGGTCGACGATCATGTCGATGATCTTCGCCGACTCGCCCAGGTGCTCGAcc < 2:1696894/144‑1 (MQ=1) aacatgatAGCCGTCGGCGCCGTTCTCGGCGGTCAGCTTCCAGTTACCTTCGTAGATGTAGGAGCTAGCTCCACGCAGCACTTCCAGACCGTCCGGCGACTGGTCGACGATCATGTCGATGATCTTCGCCGACTCGCCCAGGTGCTCGAcc < 2:2466045/144‑1 (MQ=1) cggggTAGCCGTCGGCGCCGTTTTCGGCGGTGAGCTTCCAGTTGCCTTCGTAGATGTACGAGCTGGCACCGCGCAGCACTTCCCGGCCTTCAGGCGACTGATCGACGATCATGTCGATGATCTTCGCCGACGCGCCCAGGTGCTCGACCAg < 1:5941873/148‑1 (MQ=0) catgatAGCCGTCGGCGCCGTTCTCGGCGGTCAGCTTCCAGTTACCTTCGTAGATGTAGGAGCTAGCTCCACGCAGCACTTCCAGACCGTCCGGCGACTGGTCGACGATCATGTCGATGATCTTCGCCGACTCGCCCAGGTGCTCGACCAg < 2:4578957/146‑1 (MQ=1) catgatAGCCGTCGGCGCCGTTCTCGGCGGTCAGCTTCCAGTTACCTTCGTAGATGTAGGAGCTAGCTCCACGCAGCACTTCCAGACCGTCCGGCGACTGGTCGACGATCATGTCGATGATCTTCGCCGACTCGCCCAGGTGCTCGACCAg < 2:636476/146‑1 (MQ=1) catgatAGCCGTCGGCGCCGTTCTCGGCGGTCAGCTTCCAGTTACCTTCGTAGATGTAGGAGCTAGCTCCACGCAGCACTTCCAGACCGTCCGGCGACTGGTCGACGATCATGTCGATGATCTTCGCCGACTCGCCCAGGTGCTCGACCAg > 1:3879467/6‑151 (MQ=1) catgatAGCCGTCGGCGCCGTTCTCGGCGGTCAGCTTCCAGTTACCTTCGTAGATGTAGGAGCTAGCTCCACGCAGCACTTCCAGACCGTCCGGCGACTGGTCGACGATCATGTCGATGATCTTCGCCGACTCGCCCAGGTGCTCGACCAg > 1:6490211/6‑151 (MQ=1) atgatAGCCGTCGGCGCCGTTCTCGGCGGTCAGCTTCCAGTTACCTTCGTAGATGTAGGAGCTAGCTCCACGCAGCACTTCCAGACCGTCCGGCGACTGGTCGACGATCATGTCGATGATCTTCGCCGACTCGCTCAGGTGCTCGACCAGt > 1:3190445/5‑151 (MQ=1) atgatAGCCGTCGGCGCCGTTCTCGGCGGTCAGCTTCCAGTTACCTTCGTAGATGTAGGAGCTAGCTCCACGCAGCACTTCCAGACCGTCCGGCGACTGGTCGACGATCATGTCGATGATCTTCGCCGACTCGCCCAGGTGCTCGACCAGt < 2:1244093/147‑1 (MQ=1) atgatAGCCGTCGGCGCCGTTCTCGGCGGTCAGCTTCCAGTTACCTTCGTAGATGTAGGAGCTAGCTCCACGCAGCACTTCCAGACCGTCCGGCGACTGGTCGACGATCATGTCGATGATCTTCGCCGACTCGCCCAGGTGCTCGACCAGt > 1:4142895/5‑151 (MQ=1) atgatAGCCGTCGGCGCCGTTCTCGGCGGTCAGCTTCCAGTTACCTTCGTAGATGTAGGAGCTAGCTCCACGCAGCACTTCCAGACCCTCCGGCGACTGGGCGACGATCATGTCGATGATCTTCGCCGACTCGCCCAGGTGCTCGACCAGt > 2:2325969/5‑151 (MQ=1) ttctCGGCGGTCAGCTTCCAGTTACCTTCGTAGATGTAGGAGCTAGCTCCACGCAGCACTTCCAGACCGTCCGGCGACTGGTCGACGATCATGTCGATGATCTTCGCCGACTCGCCCAGGTGCTCGACCAGTGGTTTCACATCGGCATTCa < 2:1631504/148‑1 (MQ=1) ttctCGGCGGTCAGCTTCCAGTTACCTTCGTAGATGTAGGAGCTAGCTCCACGCAGCACTTCCAGACCGTCCGGCGACTGGTCGACGATCATGTCGATGATCTTCGCCGACTCGCCCAGGTGCTCGACCAGTGGTTTCACATCGGCATTCa > 1:5137443/4‑151 (MQ=1) gtcaGCTTCCAGTTACCTTCGTAGATGTAGGAGCTAGCTCCACGCAGCACTTCCAGACCGTCCGGCGACTGGTCGACGATCATGTCGATGATCTTCGCCGACTCGCCCAGGTGCTCGACCAGTGGTTTCACATCGGCATTCAGACTGCCaa > 2:5558365/4‑151 (MQ=1) gTTGCCTTCGTAGATGTACGAGCTGGCACCGCGCAGCACTTCCAGGCCTTCAGGCGACTGATCGACGATCATGTCGATGATCTTCGCCGACTCGCCCAGGTGCTCGACCAGCGGCTTCACATCGGCGTTCAGGCTGCCGAACAAAAAGccc < 2:2945293/151‑1 (MQ=1) gTTGCCTTCGTAGATGTACGAGCTGGCACCGCGCAGCACTTCCAGGCCTTCAGGCGACTGATCGACGATCATGTCGATGATCTTCGCCGACTCGCCCAGGTGCTCGACCAGCGGCTTCACATCGGCGTTCAGGCTGCCGAACAAAAAGccc < 1:4083479/151‑1 (MQ=1) ttaccTTCGTAGATGTAGGAGCTAGCTCCACGCAGCACTTCCAGACCGTCCGGCGACTGGTCGACGATCATGTCGATGATCTTCGCCGACTCGCCCAGGTGCTCGACCAGTGGTTTCACATCGGCATTCAGACTGCCAAACAGGAAGcccc > 1:2146450/4‑151 (MQ=1) ttaccTTCGTAGATGTAGGAGCTAGCTCCACGCAGCACTTCCAGACCGTCCGGCGACTGGTCGACGATCATGTCGATGATCTTCGCCGACTCGCCCAGGTGCTCGACCAGTGGTTTCACATCGGCATTCAGACTGCCAAACAGGAAGcccc > 1:2154826/4‑151 (MQ=1) ttaccTTCGTAGATGTAGGAGCTAGCTCCACGCAGCACTTCCAGACCGTCCGGCGACTGGTCGACGATCATGTCGATGATCTTCGCCGACTCGCCCAGGTGCTCGACCAGTGGTTTCACATCGGCATTCAGACTGCCAAACAGGAAGcccc < 2:1357769/148‑1 (MQ=1) ttaccTTCGTAGATGTAGGAGCTAGCTCCACGCAGCACTTCCAGACCGTCCGGCGACTGGTCGACGATCATGTCGATGATCTTCGCCGACTCGCCCAGGTGCTCGACCAGTGGTTTCACATCGGCATTCAGACTGCCAAACAGGAAGcccc > 1:3419843/4‑151 (MQ=1) agctccacgcAGCACTTCCAGACCGTCCGGCGACTGGTCGACGATCATGTCGATGATCTTCGCCGACTCGCCCAGGTGCTCGACCCGTGGTTTCACATCGGCCTTCAGACTGCCGAACAGGAAGCCCCGGTACGACTCGAATCGTgccgcc > 1:5408972/8‑151 (MQ=1) ctccacgcAGCACTTCCAGACCGTCCGGCGACTGGTCGACGATCATGTCGATGATCTTCGCCGACTCGCCCAGGTGCTCGACCAGTGGTTTCACATCGGCATTCAGACTGCCAAACAGGAAGCCCCGCTACGACTCGAATCGTGCCACCtt > 2:1824337/6‑151 (MQ=1) tccacgcAGCACTTCCAGACCGTCCGGCGACTGGTCGACGATCATGTCGATGATCTTCGCCGACTCGCCCAGGTGCTCGACCAGTGGTTTCACATCGGCATTCAGACTGCCAAACAGGAAGCCCCGGTACGACTCGAATCGTGCCACCTTg < 1:4672185/147‑1 (MQ=1) tccacgcAGCACTTCCAGACCGTCCGGCGACTGGTCGACGATCATGTCGATGATCTTCGCCGACTCGCCCAGGTGCTCGACCAGTGGTTTCACATCGGCATTCAGACTGCCAAACAGGAAGCCCCGGTACGACTCGAATCGTGCCACCTTg > 1:1811720/5‑151 (MQ=1) ccacgcAGCACTTCCAGACCGTCCGGCGACTGGTCGACGATCATGTCGATGATCTTCGCCGACTCGCCCAGGTGCTCGACCAGTGGTTTCACATCGGCATTCAGACTGCCAAACAGGAAGCCCCGGTACGACTCGAATCGTGCCACCTTgg > 1:5602860/4‑151 (MQ=1) ccacgcAGCACTTCCAGACCGTCCGGCGACTGGTCGACGATCATGTCGATGATCTTCGCCGACTCGCCCAGGTGCTCGACCAGTGGTTTCACATCGGCATTCAGACTGCCAAACAGGAAGCCCCGGTACGACTCGAATCGTGCCACCTTgg > 1:1830504/4‑151 (MQ=1) cgggcagcaCTTCCAGGCCTTCAGGCGACTGATCGACGAGCATGTCGATGATCGTCGCCGCCTCGCCCAGGTGCTCGACCAGCGGCTTCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGt < 2:5544601/148‑1 (MQ=1) ccactcgcAGGCCTTCAGCCGACTGATCGACGATCATGTCGATGATCTTCGCCGACTCGCCCAGGTGCTCGACCAGCGGCTTCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTc < 2:4718383/144‑1 (MQ=1) accgtccggCGACTGGTCGACGATCATGTCGATGATCTTCGCCGACTCGCCCAGGTGCTCGACCAGTGGTTTCACATCGGCATTCAGACTGCCAAACAGGAAGCCCCGGTACGACTCGAATCGTGCCACCTTGGTCAGGTCGTGGGAGCCg < 2:1737343/144‑1 (MQ=0) cgtccggCGACTGGTCGACGATCATGTCGATGAGCTTCGCCGACTCGCCCAGGTGCTCGACCAGTGGTTTCACATCGGCATTCAGACTGCCAAACAGGAAGCCCCGGTACGACTCGAATCGTGCCACCTTGGTCAGGTCGTGGGAGCCGTc < 2:6164804/146‑1 (MQ=0) tccggCGACTGGTCGACGATCATGTCGATGATCTTCGCCGACTCGCCCAGGTGCTCGACCAGTGGTTTCACATCGGCATTCAGACTGCCAAACAGGAAGCCCCGGTACGACTCGAATCGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGc > 1:362111/4‑151 (MQ=1) cTGATCGACGATCATGTCGATGATCTTCGCCGACTCGCCCAGGTGCTCGACCAGCGGCTTCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAg > 2:1349807/1‑151 (MQ=11) cTGATCGACGATCATGTCGATGATCTTCGCCGACTCGCCCAGGTGCTCGACCAGCGGCTTCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAg > 1:3938430/1‑151 (MQ=11) cTGATCGACGATCATGTCGATGATCTTCGCCGACTCGCCCAGGTGCTCGACCAGCGGCTTCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAg > 1:4310440/1‑151 (MQ=11) cTGATCGACGATCATGTCGATGATCTTCGCCGACTCGCCCAGGTGCTCGACCAGCGGCTTCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGCCGCAGTTGAAg > 2:3637124/1‑151 (MQ=11) tGATCGACGATCATGTCGATGATCTTCGCCGGCTCGCCCAGGTGCTCGACCAGCGGCTTCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGc < 1:2386899/151‑1 (MQ=11) tGATCGACGATCATGTCGATGATCTTCGCCGACTCGCCCAGGTGCTCGACCAGCGGCTTCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGc < 1:2386822/151‑1 (MQ=11) tGATCGACGATCATGTCGATGATCTTCGCCGACTCGCCCAGGTGCTCGACCAGCGGCTTCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGc > 1:4449343/1‑151 (MQ=11) gATCGACGATCATGTCGATGATCTTCGCCGACTCGCCCAGGTGCTCGACCAGCGGCTTCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCt > 2:5339482/1‑151 (MQ=11) gATCGACGATCATGTCGATGATCTTCGCCGACTCGCCCAGGTGCTCGACCAGCGGCTTCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCt < 1:2942891/151‑1 (MQ=11) gATCGACGATCATGTCGATGATCTTCGCCGACTCGCCCAGGTGCTCGACCAGCGGCTTCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCt < 1:1547156/151‑1 (MQ=11) aTCGACGATCATGTCGATGATCTTCGCCGACTCGCCCAGGTGCTCGACCAGCGGCTTCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTg < 2:5498016/151‑1 (MQ=11) aTCGACGATCATGTCGATGATCTTCGCCGACTCGCCCAGGTGCTCGACCAGCGGCTTCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTg < 2:4447356/151‑1 (MQ=11) | GTTCTGGGTGGCAGCGTAGTTCCAGTGCACGGAGCTTACGTGGTAGCCGTCGGCGCCGTTTTCGGCGGTGAGCTTCCAGTTGCCTTCGTAGATGTACGAGCTGGCACCGCGCAGCACTTCCAGGCCTTCAGGCGACTGATCGACGATCATGTCGATGATCTTCGCCGACTCGCCCAGGTGCTCGACCAGTGGTTTCACATCGGCATTCAGACTGCCAAACAGGAAGCCCCGGTACGACTCGAATCGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTG > CP009974/5105162‑5105450 |
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: ‑ |