Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 884,407 | 0 | C | G | 21.7% | 49.3 / 10.7 | 23 | I245M (ATC→ATG) | mdfA | multidrug efflux system protein |
Reads supporting (aligned to +/- strand): ref base C (8/10); new base G (3/2); total (11/12) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 6.40e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
ACGCGTATAGGCGAGAAACTGTCACTGAAAGAACTCGGTCGTGACTATAAGCTGGTGCTGAAGAACGGCCGCTTTGTGGCGGGGGCGCTGGCGCTGGGATTCGTTAGTCTGCCGTTGCTGGCGTGGATCGCCCAGTCGCCGATTATCATCATTACCGGCGAGCAGTTGAGCAGCTATGAATATGGCTTGCTGCAAGTGCCTATTTTCGGGGCGTTAATTGCGGGTAACTTGCTGTTAGCGCGTCTGACCTCGCGCCGCACCGTACGTTCGCTGATTATTATGGGCGGCTGGCCGA > NC_000913/884261‑884555 | aCGCGTATAGGCGAGAAACTGTCACTGAAAGAACTCGGTCGTGACTATAAGCTGGTGCTGAAGAACGGCCGCTTTGTGGCGGGGGCGCTGGCGCTGGGATTCGTTAGTCTGCCGTTGCTGGCGTGGATCGCCCAGTCGCCGATTatcat > 1:28727/1‑149 (MQ=255) tCACTGAAAGAACTCGGTCGTGACTATAAGCTGGTGCTGAAGAACGGCCGCTTTGTGGCGGGGGCGCTGGCGCTGGGATTCGTTAGTCTGCCGTTGCTGGCGTGGATCGCCCAGTCGCCGATTATCATCATTACCGGCGAGCAGTTGa < 2:28727/148‑1 (MQ=255) tGAAAGAACTCGGTCGTGACTATAAGCTGGTGCTGAAGAACGGCCGCTTTGTGGCGGGGGCGCTGGCGCTGGGATTCGTTAGTCTGCCGTTGCTGGCGTGGATCGCCCAGTCGCCGATTATCATCATTACCGGCGAGCAGTTGagcagc > 8:141099/1‑149 (MQ=255) aaaGAACTCGGTCGTGACTATAAGCTGGTGCTGAAGAACGGCCGCTTTGTGGCGGGGGCGCTGGCGCTGGGATTCGTTAGTCTGCCGTTGCTGGCGTGGATCGCCCAGTCGCCGATTatcat < 1:195345/122‑1 (MQ=255) aaaGAACTCGGTCGTGACTATAAGCTGGTGCTGAAGAACGGCCGCTTTGTGGCGGGGGCGCTGGCGCTGGGATTCGTTAGTCTGCCGTTGCTGGCGTGGATCGCCCAGTCGCCGATTatcat > 2:195345/1‑122 (MQ=255) aaaGAACTCGGTCGTGACTATAAGCTGGTGCTGAAGAACGGCCGCTTTGTGGCGGGGGCGCTGGCGCTGGGATTCGTTAGTCTGCCGTTGCTGGCGTGGATCGCCCAGTCGCCGATTATCATCATTACCGGCGa < 8:96280/134‑1 (MQ=255) aaaGAACTCGGTCGTGACTATAAGCTGGTGCTGAAGAACGGCCGCTTTGTGGCGGGGGCGCTGGCGCTGGGATTCGTTAGTCTGCCGTTGCTGGCGTGGATCGCCCAGTCGCCGATTATCATCATTACCGGCGa > 7:96280/1‑134 (MQ=255) aaGAACTCGGTCGTGACTATAAGCTGGTGCTGAAGAACGGCCGCTTTGTGGCGGGGGCGCTGGCGCTGGGATTCGTTAGTCTGCCGTTGCTGGCGTGGATCGCCCAGTCGCCGATTATCATCATTACCGGCGAGCAGTTGAGCAGCTAt < 7:141099/149‑1 (MQ=255) gCTGGTGCTGAAGAACGGCCGCTTTGTGGCGGGGGCGCTGGCGCTGGGATTCGTTAGTCTGCCGTTGCTGGCGTGGATCGCCCAGTCGCCGATTATCATCATTACCGGCGAGCAGTTGAGCAGCTATGAATATGGCTTGCTGCAAGTGc > 5:169177/1‑149 (MQ=255) tGGTGCTGAAGAACGGCCGCTTTGTGGCGGGGGCGCTGGCGCTGGGATTCGTTAGTCTGCCGTTGCTGGCGTGGATCGCCCAGTCGCCGATTATCATCATTACCGGCGAGCAGTTGagc < 2:99429/119‑1 (MQ=255) tGGTGCTGAAGAACGGCCGCTTTGTGGCGGGGGCGCTGGCGCTGGGATTCGTTAGTCTGCCGTTGCTGGCGTGGATCGCCCAGTCGCCGATTATCATCATTACCGGCGAGCAGTTGagc > 1:99429/1‑119 (MQ=255) tGAAGAACGGCCGCTTTGTGGCGGGGGCGCTGGCGCTGGGATTCGTTAGTCTGCCGTTGCTGGCGTGGATCGCCCAGTCGCCGATTATCATCATTACCg > 1:114662/1‑99 (MQ=255) tGAAGAACGGCCGCTTTGTGGCGGGGGCGCTGGCGCTGGGATTCGTTAGTCTGCCGTTGCTGGCGTGGATCGCCCAGTCGCCGATTATCATCATTACCg < 2:114662/99‑1 (MQ=255) gaaCGGCCGCTTTGTGGCGGGGGCGCTGGCGCTGGGATTCGTTAGTCTGCCGTTGCTGGCGTGGATCGCCCAGTCGCCGATTATCATCATTACCGGCGAGCAGTTGAGCAGCTATGAATATGGCTTGCTGCAAg > 8:124004/1‑134 (MQ=255) gaaCGGCCGCTTTGTGGCGGGGGCGCTGGCGCTGGGATTCGTTAGTCTGCCGTTGCTCGCGTGGATCGCCCAGTCGCCGATTATCATCATTACCGGCGAGCAGTTGAGCAGCTATGAATATGGCTTGCTGCAAg < 7:124004/134‑1 (MQ=255) ccGCTTTGTGGCGGGGGCGCTGGCGCTGGGATTCGTTAGTCTGCCGTTGCTGGCGTGGATCGCCCAGTCGCCGATTATCATCATTACCGGCGAGCAGTTGAGCAGCTATGAATATGGCTTGCTGCAAGTGCCTATTTTCGGGGCGTTaa < 6:169177/149‑1 (MQ=255) cTTTGTGGCGGGGGCGCTGGCGCTGGGATTCGTTAGTCTGCCGTTGCTGGCGTGGATCGCCCAGTCGCCGATTATCATCATTACCGGCGAGCAGTTGAGCAGCTATGAATATGGCTTGCTGCAAGTGCCTATTTTCGGGGCGTTAATTg < 2:45550/149‑1 (MQ=255) aTTCGTTAGTCTGCCGTTGCTGGCGTGGATCGCCCAGTCGCCGATTATCATCATTACCGGCGAGCAGTTGAGCAGCTATGAATATGGCTTGCTGCAAGTGCCTATTTTCg > 2:173403/1‑110 (MQ=255) aTTCGTTAGTCTGCCGTTGCTGGCGTGGATCGCCCAGTCGCCGATTATCATCATTACCGGCGAGCAGTTGAGCAGCTATGAATATGGCTTGCTGCAAGTGCCTATTTTCg < 1:173403/110‑1 (MQ=255) tgctCGCCGGTAATGATGATAATCGGCGAGCAGTTGAGCAGCTATGAATATGGCTTGCTGCAAGTGCCTATTTTCGGGGCGTTAATTGCGGGTAACTTGCTGTTAGCGCGTCTGACGTCGCGCCGCACCGTACGTTCGCt < 8:20613/137‑1 (MQ=255) tgctCGCCGGTAATGATGATAATCGGCGAGCAGTTGAGCAGCTATGAATATGGCTTGCTGCAAGTGCCTATTTTCGGGGCGTTAATTGCGGGTAACTTGCTGTTAGCGCGTCTGACCTCGCGCCGCACCGTACGTTCGCt > 7:20613/4‑140 (MQ=255) tgctCGCCGGTAATGATGATAATCGGCGAGCAGTTGAGCAGCTATGAATATGGCTTGCTGCAAGTGCCTATTTTCGGGGCGTTAATTGCGGGTAACTTGCTGTTAGCGCGTCTGACCTCGCGCCGCACCGTACGTTCGCTGattat > 2:147763/4‑146 (MQ=255) tgctCGCCGGTAATGATGATAATCGGCGAGCAGTTGAGCAGCTATGAATATGGCTTGCTGCAAGTGCCTATTTTCGGGGCGTTAATTGCGGGTAACTTGCTGTTAGCGCGTCTGACCTCGCGCCGCACCGTACGTTCGCTGattat < 1:147763/143‑1 (MQ=255) tgctCGCCGGTAATGATGATAATCGGCGAGCAGTTGAGCAGCTATGAATATGGCTTGCTGCAAGTGCCTATTTTCGGGGCGTTAATTGCGGGTAACTTGCTGTTAGCGCGTCTGACCTCGCGCCGCACCGTACGTTCGCTGATtattat > 2:93568/4‑149 (MQ=255) aTTATCATCATTACCGGCGAGCAGTTGAGCAGCTATGAATATGGCTTGCTGCAAGTGCCTATTTTCGGGGCGTTAATTGCGGGTAACTTGCTGTTATCGAGTCTGACCTCGCGCCGCACCGTACGTTCGCTGattatt > 3:22908/1‑138 (MQ=255) aTTATCATCATTACCGGCGAGCAGTTGAGCAGCTATGAATATGGCTTGCTGCAAGTGCCTATTTTCGGGGCGTTAATTGCGGGTAACTTGCTGTTAGCGCGTCTGACCTCGCGCCGCACCGTACGTTCGCTGattatt < 4:22908/138‑1 (MQ=255) catcatTACCGGCGAGCAGTTGAGCAGCTATGAATATGGCTTGCTGCAAGTGCCTATTTTCGGGGCGTTAATTGCGGGTAACTTGCTGTTAGCGCGTCTGACCTCGCGCCGCACCGTACGTTCGCTGATTATTATGGGCGGCTGGCCGa > 5:76461/1‑149 (MQ=255) | ACGCGTATAGGCGAGAAACTGTCACTGAAAGAACTCGGTCGTGACTATAAGCTGGTGCTGAAGAACGGCCGCTTTGTGGCGGGGGCGCTGGCGCTGGGATTCGTTAGTCTGCCGTTGCTGGCGTGGATCGCCCAGTCGCCGATTATCATCATTACCGGCGAGCAGTTGAGCAGCTATGAATATGGCTTGCTGCAAGTGCCTATTTTCGGGGCGTTAATTGCGGGTAACTTGCTGTTAGCGCGTCTGACCTCGCGCCGCACCGTACGTTCGCTGATTATTATGGGCGGCTGGCCGA > NC_000913/884261‑884555 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |