Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | CP002416 | 1,820,227 | 0 | G | T | 27.6% | 121.1 / 13.5 | 58 | V23F (GTT→TTT) | Clo1313_1561 | hypothetical protein |
Reads supporting (aligned to +/- strand): ref base G (12/30); new base T (16/0); total (28/30) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 3.80e-07 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
ACAAGCTTAACAAAGGTCTTTTCCGTACAAATCTCCCTTGAAAATTATATATTTAAAACCAAAAGTAAACCAGGAGTGGTAACATGAAAAAGTTTGCAGCATTGCTGATTTTTACGGTCGTTTTTCTGCTTGTCGGATGCGGCGGTTTTGTT‑CCTC‑AAAATG‑CAGGAAAAGATACAGCACTGCCAACCGGCACCGGTGTCAATAAACCTGATACAGCAGCAAAAACAGCTGAGCCCATGACTCCTGCTGTTTCGGAGGAAAGCAGGATAAAATCATATAATCTGGGTGATATTGTTGAG > CP002416/1820078‑1820376 | aCAAGCTTAACAAAGGTCTTTTCCGTACAAATCTCCCTTGAAAATTATATATTTAAAACCAAAAGTAAACCAGGAGTGGTAACATGAAAAAGTTTGCAGCATTGCTGATTTTTACGGTCGTTTTTCTGCTTGTCGGATGCGGCGGTTTTg < 1:196718/150‑1 (MQ=255) aaaGGTCTTTTCCGTACAAATCTCCCTTGAAAATTATATATTTAAAACCAAAAGTAAACCAGGAGTGGTAACATGAAAAAGTTTGCAGCATTGCTGATTTTTACGGTCGTTTTTCTGCTTGTCGGATGCGGCGGTTTTGTT‑CCTC‑AAAAt < 1:1304144/150‑1 (MQ=255) aGGTCTTTTCCGTACAAATCTCCCTTGAAAATTATATATTTAAAACCAAAAGTAAACCAGGAGTGGTAACATGAAAAAGTTTGCAGCATTGCTGATTTTTACGGTCGTTTTTCTGCTTGTCGGATGCGGCGGTTTTGTT‑CCTC‑AAAATG‑c < 1:631669/150‑1 (MQ=255) gTCTTTTCCGTACAAATCTCCCTTGAAAATTATATATTTAAAACCAAAAGTAAACCAGGAGTGGTAACATGAAAAAGTTTGCAGCATTGCTGATTTTTACGGTCGTTTTTCTGCTTGTCGGATGCGGCGGTTTTGTT‑CCTC‑AAAATG‑CAg < 1:1758804/150‑1 (MQ=255) gTCTTTTCCGTACAAATCTCCCTTGAAAATTATATATTTAAAACCAAAAGTAAACCAGGAGTGGTAACATGAAAAAGTTTGCAGCATTGCTGATTTTTACGGTCGTTTTTCTGCTTGTCGGATGCGGCGGTTTTGTT‑CCTC‑AAAATG‑CAg < 1:886116/150‑1 (MQ=255) tCTTTTCCGTACAAATCTCCCTTGAAAATTATATATTTAAAACCAAAAGTAAACCAGGAGTGGTAACATGAAAAAGTTTGCAGCATTGCTGATTTTTACGGTCGTTTTTCTGCTTGTCGGATGCGGCGGTTTTGTT‑CCTC‑AAAATG‑CAgg < 1:320846/150‑1 (MQ=255) tttCCGTACAAATCTCCCTTGAAAATTATATATTTAAAACCAAAAGTAAACCAGGAGTGGTAACATGAAAAAGTTTGCAGCATTGCTGATTTTTACGGTCGTTTTTCTGCTTGTCGGATGCGGCGGTTTTGTT‑CCTC‑AAAATG‑CAGGaaa < 1:426192/150‑1 (MQ=255) tttCCGTACAAATCTCCCTTGAAAATTATATATTTAAAACCAAAAGTAAACCAGGAGTGGTAACATGAAAAAGTTTGCAGCATTGCTGATTTTTACGGTCGTTTTTCTGCTTGTCGGATGCGGCGGTTTTGTT‑CCTC‑AAAATG‑CAGGaaa < 1:426170/150‑1 (MQ=255) aCAAATCTCCCTTGAAAATTATATATTTAAAACCAAAAGTAAACCAGGAGTGGTAACATGAAAAAGTTTGCAGCATTGCTGATTTTTACGGTCGTTTTTCTGCTTGTCGGATGCGGCGGTTTTTTT‑TCTC‑AAAATG‑TAGGgaaa > 1:1111076/1‑144 (MQ=255) cAAATCTCCCTTGAAAATTATATATTTAAAACCAAAAGTAAACCAGGAGTGGTAACATGAAAAAGTTTGCAGCATTGCTGATTTTTACGGTCGTTTTTCTGCTTGTCGGATGCGGCGGTTTTGTT‑CCTC‑AAAATG‑CAGGAAAAGATACAg < 1:153372/150‑1 (MQ=255) ctcCCTTGAAAATTATATATTTAAAACCAAAAGTAAACCAGGAGTGGTAACATGAAAAAGTTTGCAGCATTGCTGATTTTTACGGTCGTTTTTCTGCTTGTCGGATGCGGCGGTTTTTTT‑CCTC‑AAAATG‑CAGGAAAAGATACa > 1:1373557/1‑144 (MQ=255) ttGAAAATTATATATTTAAAACCAAAAGTAAACCAGGAGTGGTAACATGAAAAAGTTTGCAGCATTGCTGATTTTTACGGTCGTTTTTCTGCTTGTCGGATGCGGCGGTTTTGTTTCCCC‑AAAA‑G‑CAGGAAAAGaaaaa > 1:474548/1‑135 (MQ=255) tGAAAATTATATATTTAAAACCAAAAGTAAACCAGGAGTGGTAACATGAAAAAGTTTGCAGCATTGCTGATTTTTACGGTCGTTTTTCTGCTTGTCGGATGCGGCGGTTTTGGT‑CCTCCAAAAGG‑AAGGaaaaat > 1:2137793/1‑133 (MQ=255) aaTTATATATTTAAAACCAAAAGTAAACCAGGAGTGGTAACATGAAAAAGTTTGCAGCATTGCTGATTTTTACGGTCGTTTTTCTGCTTGTCGGATGCGGCGGTTTTGTT‑TCTC‑AAAAAG‑CAGGgaaa > 1:365596/1‑128 (MQ=255) aTTATATATTTAAAACCAAAAGTAAACCAGGAGTGGTAACATGAAAAAGTTTGCAGCATTGCTGATTTTTACGGTCGTTTTTCTGCTTGTCGGATGCGGCGGTTTTGTT‑CCTC‑AAAATG‑CAGGAAAAGATACAGCACTGCCAACCGGCAc < 1:726569/150‑1 (MQ=255) ttATATATTTAAAACCAAAAGTAAACCAGGAGTGGTAACATGAAAAAGTTTGCAGCATTGCTGATTTTTACGGTCGTTTTTCTGCTTGTCGGATGCGGCGGTTTTTTT‑TCTC‑AAAATG‑CAGGAAAAGa > 1:953911/1‑128 (MQ=255) ttATATATTTAAAACCAAAAGTAAACCAGGAGTGGTAACATGAAAAAGTTTGCAGCATTGCTGATTTTTACGGTCGTTTTTCTGCTTGTCGGATGCGGCGGTTTTGTT‑CCTC‑AAAATG‑CAGGAAAAGATACAGCACTGCCAACCGGCAcc < 1:536132/150‑1 (MQ=255) tatatTTAAAACCAAAAGTAAACCAGGAGTGGTAACATGAAAAAGTTTGCAGCATTGCTGATTTTTACGGTCGTTTTTCTGCTTGTCGGATGCGGCGGTTTTTTT‑TCTC‑CAAATG‑CAGGgaaa > 1:1311431/1‑123 (MQ=255) tatatTTAAAACCAAAAGTAAACCAGGAGTGGTAACATGAAAAAGTTTGCAGCATTGCTGATTTTTACGGTCGTTTTTCTGCTTGTCGGATGCGGCGGATTTGTT‑CCCC‑aa > 1:20443/1‑111 (MQ=255) atTTAAAACCAAAAGTAAACCAGGAGTGGTAACATGAAAAAGTTTGCAGCATTGCTGATTTTTACGGTCGTTTTTCTGCTTGTCGGATGCGGCGGTTTTGTT‑CCTC‑AAAAAG‑CAGGGAAAGATACa > 1:671053/1‑126 (MQ=255) aaaCCAAAAGTAAACCAGGAGTGGTAACATGAAAAAGTTTGCAGCATTGCTGATTTTTACGGTCGTTTTTCTGCTTGTCGGATGCGGCGGTTTTGTT‑TCTC‑AAAAAG‑CAGGAAAAGATACa > 1:1743129/1‑121 (MQ=255) ccAAAAGTAAACCAGGAGTGGTAACATGAAAAAGTTTGCAGCATTGCTGATTTTTACGGTCGTTTTTCTGCTTGTCGGATGCGGCGGTTTTTTT‑TCTC‑AAAAAG‑CAGGgaaaaa > 1:342371/1‑112 (MQ=255) aaaaGTAAACCAGGAGTGGTAACATGAAAAAGTTTGCAGCATTGCTGATTTTTACGGTCGTTTTTCTGCTTGTCGGATGCGGCGGTTTTTTT‑TCTC‑CAAAAG‑CAGGaaaa > 1:1677536/1‑110 (MQ=255) aaaGTAAACCAGGAGTGGTAACATGAAAAAGTTTGCAGCATTGCTGATTTTTACGGTCGTTTTTCTGCTTGTCGGATGCGGCGGTTTTTTT‑TCTCCAAAAT‑‑CAGGaaaa > 1:1180688/1‑109 (MQ=255) aaaGTAAACCAGGAGTGGTAACATGAAAAAGTTTGCAGCATTGCTGATTTTTACGGTCGTTTTTCTGCTTGTCGGATGCGGCGGTTTTGTT‑CCTC‑AAAATG‑CAGGAAAAGATACAGCACTGCCAACCGGCACCGGTGTCAATAAACCTGa < 1:1719005/150‑1 (MQ=255) aGGAGTGGTAACATGAAAAAGTTTGCAGCATTGCTGATTTTTACGGTCGTTTTTCTGCTTGTCGGATGCGGCGGTTTTTTT‑TCTC‑AAAATG‑CAGGAAAAAATACAGCCCTGcc > 1:501079/1‑113 (MQ=255) aGTGGTAACATGAAAAAGTTTGCAGCATTGCTGATTTTTACGGTCGTTTTTCTGCTTGTCGGATGCGGCGGTTTTGTT‑CCTC‑AAAATG‑CAGGAAAAGATACAGCACTGCCAACCGGCACCGGTGTCAATAAACCTGATACAGCAGCaaaa < 1:1726588/150‑1 (MQ=255) tGGTAACATGAAAAAGTTTGCAGCATTGCTGATTTTTACGGTCGTTTTTCTGCTTGTCGGATGCGGCGGTTTTGTT‑CCTC‑AAAATG‑CAGGAAAAGATACAGCACTGCCAACCGGCACCGGTGTCAATAAACCTGATACAGCAGCAAAAAc < 1:121589/150‑1 (MQ=255) cATGAAAAAGTTTGCAGCATTGCTGATTTTTACGGTCGTTTTTCTGCTTGTCGGATGCGGCGGTTTTGTT‑TCCC‑AAAATG‑GAGGAAAAAATAAAACACTGCCaa > 1:1946437/1‑104 (MQ=255) tGAAAAAGTTTGCAGCATTGCTGATTTTTACGGTCGTTTTTCTGCTTGTCGGATGCGGCGGCTTTTTT‑TCTC‑CAAAAG‑CAGGGAAAAATa > 1:334261/1‑90 (MQ=255) aaaaGTTTGCAGCATTGCTGATTTTTACGGTCGTTTTTCTGCTTGTCGGATGCGGCGGCTTTGGT‑cc > 1:2017354/1‑67 (MQ=255) aaaGTTTGCAGCATTGCTGATTTTTACGGTCGTTTTTCTGCTTGTCGGATGCGGCGGTTTTTTT‑TCCC‑CAAATG‑CAAGAAAAGATaaa > 1:1312232/1‑86 (MQ=255) aaaGTTTGCAGCATTGCTGATTTTTACGGTCGTTTTTCTGCTTGTCGGATGCGGCGGTTTTGTT‑CCTC‑AAAATG‑CAGGAAAAGATACAGCACTGCCAACCGGCACCGGTGTCAATAAACCTGATACAGCAGCAAAAACAGCTGAGCCCAt < 1:448558/150‑1 (MQ=255) gTTTGCAGCATTGCTGATTTTTACGGTCGTTTTTCTGCTTGTCGGATGCGGCGGTTTTGTT‑CCTC‑AAAATG‑CAGGAAAAGATACAGCACTGCCAACCGGCACCGGTGTCAATAAACCTGATACAGCAGCAAAAACAGCTGAGCCCATGAc < 1:19730/150‑1 (MQ=255) gcagcaTTGCTGATTTTTACGGTCGTTTTTCTGCTTGTCGGATGCGGCGGTTTTGTT‑CCTC‑AAAATG‑CAGGAAAAGATACAGCACTGCCAACCGGCACCGGTGTCAATAAACCTGATACAGCAGCAAAAACAGCTGAGCCCATGACTCct < 1:438308/150‑1 (MQ=255) gcagcaTTGCTGATTTTTACGGTCGTTTTTCTGCTTGTCGGATGCGGCGGTTTTGTT‑CCTC‑AAAATG‑CAGGAAAAGATACAGCACTGCCAACCGGCACCGGTGTCAATAAACCTGATACAGCAGCAAAAACAGCTGAGCCCATGACTCct < 1:1660861/150‑1 (MQ=255) caTTGCTGATTTTTACGGTCGTTTTTCTGCTTGTCGGATGCGGCGGTTTTTTT‑TCTC‑CAAATGTCGGGAAAA‑AAAAAGCACTGCCaaac > 1:1184270/1‑87 (MQ=255) ttGCTGATTTTTACGGTCGTTTTTCTGCTTGTCGGATGCGGCGGTTTTGTT‑CCTC‑AAAATG‑CAGGAAAAGATACAGCACTGCCAACCGGCACCGGTGTCAATAAACCTGATACAGCAGCAAAAACAGCTGAGCCCATGACTCCTGCTGtt < 1:256921/150‑1 (MQ=255) tGCTGATTTTTACGGTCGTTTTTCTGCTTGTCGGATGCGGCGGTTTTGTT‑CCTC‑AAAATG‑CAGGAAAAGATACAGCACTGCCAACCGGCACCGGTGTCAATAAACCTGATACAGCAGCAAAAACAGCTGAGCCCATGACTCCTGCTGttt < 1:21270/150‑1 (MQ=255) gCTGATTTTTACGGTCGTTTTTCTGCTTGTCGGATGCGGCGGTTTTGTT‑CCTC‑AAAATG‑CAGGAAAAGATACAGCACTGCCAACCGGCACCGGTGTCAATAAACCTGATACAGCAGCAAAAACAGCTGAGCCCATGACTCCTGCTGTTTc < 1:883715/150‑1 (MQ=255) gATTTTTACGGTCGTTTTTCTGCTTGTCGGATGCGGCGGTTTTTTT‑TCCC‑AAAAAG‑CACGAAAAGATACAGc > 1:1966763/1‑72 (MQ=255) ttttACGGTCGTTTTTCTGCTTGTCGGATGCGGCGGTTTTGTT‑CCTC‑AAAATG‑CAGGAAAAGATACAGCACTGCCAACCGGCACCGGTGTCAATAAACCTGATACAGCAGCAAAAACAGCTGAGCCCATGACTCCTGCTGTTTCggagga < 1:1700400/150‑1 (MQ=255) tttACGGTCGTTTTTCTGCTTGTCGGATGCGGCGGTTTTGTT‑CCTC‑AAAATG‑CAGGAAAAGATACAGCACTGCCAACCGGCACCGGTGTCAATAAACCTGATACAGCAGCAAAAACAGCTGAGCCCATGACTCCTGCTGTTTCGGAGGaa < 1:376183/150‑1 (MQ=255) tCGTTTTTCTGCTTGTCGGATGCGGCGGTTTTGTT‑CCTC‑AAAATG‑CAGGAAAAAAAACAACACTGCCAACCGGCACCGGGGTCAATAAAACTGATACAGCAGCAAAAAAAACTGAGccc > 1:703409/1‑119 (MQ=255) cGTTTTTCTGCTTGTCGGATGCGGCGGTTTTTTT‑CCTC‑AAAAAG‑CAGGGAAAAATAAAGCACTGCCCACCCGCACCCGTGTCAAAAAACCTGATACAACAGCaaaaa > 1:905373/1‑107 (MQ=255) tttCTGCTTGTCGGATGCGGCGGTTTTTTT‑TCTC‑AAAAAG‑CAGGAAAAAATACAACACTGCCCACCCGCACCCGGGTCAATAAAACTGATACAACAGCAAAAAAAGCTGAGCCCATTACTCCTGCTGTTTCGGGGGGAAGCAGGATaaaa > 1:1181439/1‑150 (MQ=255) ttCTGCTTGTCGGATGCGGCGGTTTTGTT‑CCTC‑AAAATG‑CAGGAAAAGATACAGCACTGCCAACCGGCACCGGTGTCAATAAACCTGATACAGCAGCAAAAACAGCTGAGCCCATGACTCCTGCTGTTTCGGAGGAAAGCAGGATAAAAt < 1:987878/150‑1 (MQ=255) gCTTGTCGGATGCGGCGGTTTTTTT‑TCTC‑CAAAAG‑CCGGTAAAAAAACAGCACTGCCCAcc > 1:1634218/1‑61 (MQ=255) cTTGTCGGATGCGGCGGTTTTTTT‑CTTC‑CAAATG‑CTTGAAAAGATTCATCACTGccca > 1:394579/1‑56 (MQ=255) cTTGTCGGATGCGGCGGTTTTGTT‑CCTC‑AAAATG‑CAGGAAAAGATACAGCACTGCCAACCGGCACCGGTGTCAATAAACCTGATACAGCAGCAAAAACAGCTGAGCCCATGACTCCTGCTGTTTCGGAGGAAAGCAGGATAAAATCatat < 1:2154205/150‑1 (MQ=255) ttGTCGGATGCGGCGGTTTTGTT‑CCTC‑AAAATG‑CAGGAAAAGATACAGCACTGCCAACCGGCACCGGTGTCAATAAACCTGATACAGCAGCAAAAACAGCTGAGCCCATGACTCCTGCTGTTTCGGAGGAAAGCAGGATAAAATCAtata < 1:498406/150‑1 (MQ=255) tGTCGGATGCGGCGGTTTTGTT‑CCTC‑AAAATG‑CAGGAAAAGATACAGCACTGCCAACCGGCACCGGTGTCAATAAACCTGATACAGCAGCAAAAACAGCTGAGCCCATGACTCCTGCTGTTTCGGAGGAAAGCAGGATAAAATCATATaa < 1:1807981/150‑1 (MQ=255) ggTTTTGTT‑CCTC‑AAAATG‑CAGGAAAAGATACAGCACTGCCAACCGGCACCGGTGTCAATAAACCTGATACAGCAGCAAAAACAGCTGAGCCCATGACTCCTGCTGTTTCGGAGGAAAGCAGGATAAAATCATATAATCTGGGTGATAtt < 1:669053/150‑1 (MQ=255) ggTTTTGTT‑CCTC‑AAAATG‑CAGGAAAAGATACAGCACTGCCAACCGGCACCGGTGTCAATAAACCTGATACAGCAGCAAAAACAGCTGAGCCCATGACTCCTGCTGTTTCGGAGGAAAGCAGGATAAAATCATATAATCTGGGTGATAtt > 1:1445173/1‑150 (MQ=255) ttttGTT‑CCTC‑AAAATG‑CAGGAAAAGATACAGCACTGCCAACCTGCACCGGTGTCAATAAACCTGATACAGCAGCAAAAATAGCTGAGTCCATGACTCCTGCTGTTTCGGAGGAAAGCAGGATAAAACCATATAATCTGGGTGATAgtgt > 1:2092074/1‑150 (MQ=255) ttttGTT‑CCTC‑AAAATG‑CAGGAAAAGATACAGCACTGCCAACCGGCACCGGTGTCAATAAACCTGATACAGCAGCAAAAACAGCTGAGCCCATGACTCCTGCTGTTTCGGAGGAAAGCAGGATAAAATCATATAATCTGGGTGATAttgt < 1:21201/150‑1 (MQ=255) ttGTT‑CCTC‑AAAATG‑CAGGAAAAGATACAGCACTGCCAACCGGCACCGGTGTCAATAAACCTGATACAGCAGCAAAAACAGCTGAGCCCATGACTCCTGCTGTTTCGGAGGAAAGCAGGATAAAATCATATAATCTGGGTGATAttgttg < 1:2215622/150‑1 (MQ=255) ttGTT‑CCTC‑AAAATG‑CAGGAAAAGATACAGCACTGCCAACCGGCACCGGTGTCAATAAACCTGATACAGCAGCAAAAACAGCTGAGCCCATGACTCCTGCTGTTTCGGAGGAAAGCAGGATAAAATCATATAATCTGGGTGATAttgttg > 1:228170/1‑150 (MQ=255) tGTT‑CCTC‑AAAATG‑CAGGAAAAGATACAGCACTGCCAACCGGCACCGGTGTCAATAAACCTGATACAGCAGCAAAAACAGCTGAGCCCATGACTCCTGCTGTTTCGGAGGAAAGCAGGATAAAATCATATAATCTGGGTGATATTGTTGa < 1:334026/150‑1 (MQ=255) gTT‑CCTC‑AAAATG‑CAGGAAAAGATACAGCACTGCCAACCGGCACCGGTGTCAATAAACCTGATACAGCAGCAAAAACAGCTGAGCCCATGACTCCTGCTGTTTCGGAGGAAAGCAGGATAAAATCATATAATCTGGGTGATATTGTTGAg < 1:2037327/150‑1 (MQ=255) gTT‑CCTC‑AAAATG‑CAGGAAAAGATACAGCACTGCCAACCGGCACCGGTGTCAATAAACCTGATACAGCAGCAAAAACAGCTGAGCCCATGACTCCTGCTGTTTCGGAGGAAAGCAGGATAAAATCATATAATCTGGGTGATATTGTTGAg > 1:1312155/1‑150 (MQ=255) gTT‑CCTC‑AAAATG‑CAGGAAAAGATACAGCACTGCCAACCGGCACCGGTGTCAATAAACCTGATACAGCAGCAAAAACAGCTGAGCCCATGACTCCTGCTGTTTCGGAGGAAAGCAGGATAAAATCATATAATCTGGGTGATATTGTTGAg < 1:1311432/150‑1 (MQ=255) gTT‑CCTC‑AAAATG‑CAGGAAAAGATACAGCACTGCCAACCGGCACCGGTGTCAATAAACCTGATACAGCAGCAAAAACAGCTGAGCCCATGACTCCTGCTGTTTCGGAGGAAAGCAGGATAAAATCATATAATCTGGGTGATATTGTTGAg < 1:339007/150‑1 (MQ=255) | ACAAGCTTAACAAAGGTCTTTTCCGTACAAATCTCCCTTGAAAATTATATATTTAAAACCAAAAGTAAACCAGGAGTGGTAACATGAAAAAGTTTGCAGCATTGCTGATTTTTACGGTCGTTTTTCTGCTTGTCGGATGCGGCGGTTTTGTT‑CCTC‑AAAATG‑CAGGAAAAGATACAGCACTGCCAACCGGCACCGGTGTCAATAAACCTGATACAGCAGCAAAAACAGCTGAGCCCATGACTCCTGCTGTTTCGGAGGAAAGCAGGATAAAATCATATAATCTGGGTGATATTGTTGAG > CP002416/1820078‑1820376 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |