Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | KT2440 | 2,938,239 | 0 | A | C | 52.2% | 10.0 / 20.2 | 23 | K537N (AAA→AAC) | pqiB | Intermembrane transport protein PqiB |
Reads supporting (aligned to +/- strand): ref base A (6/5); new base C (0/12); total (6/17) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 4.58e-03 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.58e-03 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
ATTTCCGAAGATTCGCCGCAGCGCGAGCGAATGGGTGAAACGCTGGACGAGCTGGAGCGCATGTCGCGTTCGCTACGTGACCTTGCCGATTACCTGGGCCGCCACCCCGAGTCGCTGATTCGTGGCCGGCCGAAAGCGGCCGGGTCCGCAGACCTTGAACCCTGAGGAGACAGCCGATGCAACGACTGCGCAATGCATGTGTGGCCAGCCTGCTGGTGTGGCTCGGGGGCTGCAGCAGTACGCCGAACAACTACCACACGTTGGTAC > KT2440/2938105‑2938371 | aTTTCCGAAGATTCGCCGCAGCGCGAGCGAATGGGTGAAACGCTGGACGAGCTGGAGCGCATGTCGCGTTCGCTACGTGACCTTGCCGATTACCTGGGCCGCCACCCCGAGTCGCTGGTTCGTGGCCGGCCGcaa > 1:112601/1‑135 (MQ=255) aTTTCCGAAGATTCGCCGCAGCGCGAGCGAATGGGTGAAACGCTGGACGAGCTGGAGCGCATGTCGCGTTCGCTACGTGACCTTGCCGATTACCTGGGCCGCCACCCCGAGTCGCTGATTCGTGGCCGGCCGaaa > 3:74855/1‑135 (MQ=255) cgccgcAGCGCGAGCGAATGGGTGAAACGCTGGACGAGCTGGAGCGCATGTCGCGTTCGCTACGTGACCTTGCCGATTACCTGGGCCGCCACCCCGAGTCGCTGATTCGTGGCCGGCCGAAAGCGGCCGGGTcc < 8:299083/134‑1 (MQ=255) gcAGCGCGAGCGAATGGGTGAAACGCTGGACGAGCTGGAGCGCATGTCGCGTTCGCTACGTGACCTTGCCGATTACCTGGGCCGCCACCCCGAGTCGCTGATTGGTGGCCGGCCGCAAGCGGCCCGGGGCGCAGa > 5:332982/1‑135 (MQ=255) ccgcgcGAGCTATGGGGTGAAACGCGGGCCGAGCTGGAGCGCTTGTCGCGTGCGCTACGTGACCTTGCCCATTCCCGGGCCCCCCACCCCGAGGCGCTGATGCGGGCCCGGCCGAACGCGGCCGGGTCCGCAGAc < 2:112601/133‑1 (MQ=255) aaCGCGGGACGAGCTGGAGCGCAGGTGGCGTTCGATACGTGACCTTGCCGATTACCTGGGCCGCCACCCCGAGTCGCTGATTCGGGGCCGGCCGAACGCGGCCGGGTCCGCAGACCTTGAACCCTGAGGAGACa < 8:240097/134‑1 (MQ=255) gttgCTGGAGCTAATGTCGGGTGCGCTACGGTACCTGGCCTATTACCTGGGCCCCCCCCCCGTGTCGTTGATGCGGGGCCGGCCCAAAGCGGCCGGGTCCGCAGACCTTGAACCCTGAGGAGACAGCCGATGCaa < 6:332982/132‑1 (MQ=255) cGAGCTGGCGCGCAGGGCGCGTGCGCTCCGTGACCTTGCCGATACCCGGGCCCGCCACCCCGAGGCGCTGATTCGTGCCCGGCCGAACGCGGCCGGGTCCGCAGACCTTGAACCCTGAGGAGACAGCCGATGCaa < 7:284450/135‑1 (MQ=255) gggAGCGCTTGTCGCGTTCGCTACGTGACGTTGCCGATTACCGGGGCCGCCACCCGGAGTCGGTGATTCGGGCCCGGCCTAACGCGGCCGGGTCCGCAGACCTTGAACCCTGAGGAGACAGCCGATGCAACGACt < 3:185805/134‑1 (MQ=255) gtgcGCTCCGTCACCTGGCCGATTCCCGGGGCCGCCACCCCGAGCCGCTGATGCGGGGCCGGCCGAACGCGGCCGGGTCCGCAGACCTTGAACCCTGAGGAGACAGCCGATGCAACGACTGCGCAATGCAtgtgt < 4:155291/132‑1 (MQ=255) tCGCTACGGCACCTTGCCGATAACCGGGGCCCCCCCCCCGAGTCGCTCATGCGTGGCCGGCCGAAAGCGGCCGGGTCCGCAGACCTTGAACCCTGAGGAGACAGCCGATGCAACGACTGCGCAATGCATGTGTgg < 8:201443/135‑1 (MQ=255) gttcccTTGCCGATAACAGGGCCCCCCACCCCGCTGCGCTGATGGGGGGCCGGCCTAACGCGGCCGGGTCCGCAGACCTTGAACCCTGAGGAGACAGCCGATGCAACGACTGCGCAATGCATGTGTGGCCAGc < 8:52597/129‑1 (MQ=255) gggACCTTGCCGATTCCCTGGGCCGCCACCCCGCGGCGGTGTTGCGGGGCCGGCCTAACGCGGCCGGGTCCGCAGACCTTGAACCCTGAGGAGACAGCCGATGCAACGACTGCGCAATGCATGTGTGGCCAGc < 8:46519/131‑1 (MQ=255) aCCTGGCCGATTCCCGGGGCCGCCCCCCGGAGTCGCTGATGGGGGGCGGCCCTAACGCGGCCGGGTCCGCAGACCTTGAACCCTGAGGAGACAGCCGATGCAACGACTGCGCAATGCATGTGTGGCCAGCctgct < 3:236601/135‑1 (MQ=255) cGATTACCTGGGCCGCCACCCCGAGTCGCTGATTCGTGGCCGGCCGAAAGCGGCCGGGTCCGCAGCCCTTGCACCCTGCGGAGACAGCCGATGCAACCACTGCGCCATGCATGTGTGGCCCGCCGGCTGGTGTgg > 4:35882/1‑135 (MQ=255) cGATTACCT‑GGCCGCCACCCCGAGTCGCTGATTAGTGGCCCGCCGAAAGCGGCCGGGTCCGCACACCTTTAAACCTGAGGAGACAGCCGATGCAACGACTGCGCAATGCATTTGTGGCCAGCCGGCTGGTGTGGc > 7:242840/1‑135 (MQ=255) gggcccgccCCCCCGCGGCGCTGATGGGGGGCCGCCCCAACGCGGCCGGGTCCGCAGACCTTGAACCCTGAGGAGACAGCCGATGCAACGACTGCGCAATGCATGTGTGGCCAGCCTGCTGGTGTGGCTCggggg < 5:179573/131‑1 (MQ=255) ccgccACCCCGAGTCGCTGATTCGTGGCCGGCCGAAAGCGGCCGGGTCCGCAGACCTTGAACCCTGAGGAGACAGCCGATGCAACGACTGCGCAATGCATGTGTGGCCAGCCTGCTGGTGTGGCTCGGGGGCt > 7:279486/1‑133 (MQ=255) cgccACCCCGAGTCGCGGATTCGTGGCCGGCCGAAAGCGGCCGGGTCCGCCGCCCGTGTACCCGGTGGAGACACCCGATGGCACGCCTGCGGCATGGATGTGTGGCCCGCCGGCTGGTGTGGGGCggggggggca > 6:103622/1‑130 (MQ=255) cccGAGTCGCTGATTCGTGGCCGGCCGAAAGCGGCCGGGTCCGCAGACCTTGAACCCTGAGGAGACAGCCGCTGCAACGACTGCGCAATGCATGTGTGGCCCGCCTGCTGGTGTGGCTCGGGGGCTGCCGCAGTa > 3:344311/1‑135 (MQ=255) ggAGTCGCTGTTTCGGGGCGGGCCGAACGCGGCCGGGTCCGCAGACCTTGAACCCTGAGGAGACAGCCGATGCAACGACTGCGCAATGCATGTGTGGCCAGCCTGCTGGTGTGGCTCGGGGGCTGCAGCAGTACg < 4:74855/134‑1 (MQ=255) gTCGCGGTTTCGTGCCCGGCCGAAAGCGGCCGGGTCCGCAGACCTTGAACCCTGAGGAGACAGCCGATGCTACGACTGCGCAATGCATGTGTGGCCAGCCTGCTGGTGTGGCTCGGGGGCTGCAGCAGTACGCCg < 7:233563/135‑1 (MQ=255) tgcggggccggccGAAAGCGGCCGGGTCCGCAGACCTTGAACCCTGAGGAGACAGCCGATGCAACGACTGCGCAATGCATGTGTGGCCAGCCTGCTGGTGTGGCTCGGGGGCTGCAGCAGTACGCCGAACAACTa < 7:330465/130‑1 (MQ=255) cggggcccgccAAACGCGGCCGGGTCCGCAGACCTTGAACCCTGAGGAGACAGCCGATGCAACGACTGCGCAATGCATGTGTGGCCAGCCTGCTGGTGTGGCTCGGGGGCTGCAGCAGTACGCCGAACAACTAcc < 7:140874/132‑1 (MQ=255) ggggcgggccGAACGCGGCCGGGTCCGCAGACCTTGAACCCTGAGGAGACAGCCGATGCAACGACTGCGCAATGCATGTGTGGCCAGCCTGCTGGTGTGGCTCGGGGGCTGCAGCAGTACGCCGAACAACTAcc < 4:19915/132‑1 (MQ=255) aaaGCGGCCGGGTCCGCAGACCTTGAACCCTGAGGAGACAGCCGATGCAACGACTGCGCAATGCATGTGTGGCCAGCCTGCTGGTGTGGCTCggggg < 4:31351/97‑1 (MQ=255) aaaGCGGCCGGGTCCGCAGACCTTGAACCCTGAGGAGACAGCCGATGCAACGACTGCGCAATGCATGTGTGGCCAGCCTGCTGGTGTGGCTCggggg > 3:31351/1‑97 (MQ=255) aaaGCGGCCGGGTCCGCAGACCTTGAACCCTGAGGAGACAGCCGATGCAACGACTGCGCAATGCATGTGTGGCCAGCCTGCTGGTGTGGCTCGGGGGCTGCAGCAGTCCGCCGCACAACTACCACACGTTGGTAc > 8:67537/1‑135 (MQ=255) | ATTTCCGAAGATTCGCCGCAGCGCGAGCGAATGGGTGAAACGCTGGACGAGCTGGAGCGCATGTCGCGTTCGCTACGTGACCTTGCCGATTACCTGGGCCGCCACCCCGAGTCGCTGATTCGTGGCCGGCCGAAAGCGGCCGGGTCCGCAGACCTTGAACCCTGAGGAGACAGCCGATGCAACGACTGCGCAATGCATGTGTGGCCAGCCTGCTGGTGTGGCTCGGGGGCTGCAGCAGTACGCCGAACAACTACCACACGTTGGTAC > KT2440/2938105‑2938371 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |