Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 3,525,012 | 1 | . | G | 20.0% | 56.7 / 19.0 | 25 | coding (2142/2553 nt) | mrcA | peptidoglycan glycosyltransferase/peptidoglycan DD‑transpeptidase MrcA |
Reads supporting (aligned to +/- strand): ref base . (3/17); new base G (5/0); total (8/17) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.05e-03 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.23e-01 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
ACCGCACGTCATCAACACTCCGCTGGCATTCCTGATTAAGAGTGCTTTGAACACCAATATCTTTGGTGAGCCAGGCTGGCAGGGTACTGGCTGGCGTGCAGGTCGTGATTTGCAGCGTCGCGATATCGGC‑GGGAAAACCGGGACCACTAACAGTTCGAAAGATGCGTGGTTCTCGGGTTACGGTCCGGGCGTTGTGACCTCGGTCTGGATTGGCTTTGATGATCACCGTCGTAATCTCGGTCA > NC_000913/3524883‑3525125 | accgcacGTCATCAACACTCCGCTGGCATTCCTGATTAAGAGTGCTTTGAACACCAATATCTTTGGTGAGCCAGGCTGGCAGGGTACTGGCTGGCGTGCAGGTCGTGATTTGCAGCGTCGCGATATCGGC‑GGGAAAACCg < 1:255045/140‑1 (MQ=255) accgcacGTCATCAACACTCCGCTGGCATTCCTGATTAAGAGTGCTTTGAACACCAATATCTTTGGTGAGCCAGGCTGGCAGGGTACTGGCTGGCGTGCAGGTCGTGATTTGCAGCGTCGCGATATCGGC‑GGGAAAACCg < 2:796981/140‑1 (MQ=255) cgcacGTCATCAACACTCCGCTGGCATTCCTGATTAAGAGTGCTTTGAACACCAATATCTTTGGTGAGCCAGGCTGGCAGGGTACTGGCTGGCGTGCAGGTCGTGATTTGCAGCGTCGCGATATCGGCGGGGAAAACCgg > 2:884074/1‑140 (MQ=255) cgcacGTCATCAACACTCCGCTGGCATTCCTGATTAAGAGTGCTTTGAACACCAATATCTTTGGTGAGCCAGGCTGGCAGGGTACTGGCTGGCGTGCAGGTCGTGATTTGCAGCGTCGCGATATCGGC‑GGGAAAACCggg > 1:521308/1‑140 (MQ=255) cgcacGTCATCAACACTCCGCTGGCATTCCTGATTAAGAGTGCTTTGAACACCAATATCTTTGGTGAGCCAGGCTGGCAGGGTACTGGCTGGCGTGCAGGTCGTGATTTGCAGCGTCGCGATATCGGC‑GGGAAAACCggg < 2:603535/140‑1 (MQ=255) cgcacGCCATCAACACTCCGCTGGGATTCCTGATTAAGAGTGCTTTTAACACCCATATCTTTGGTGAGCCAGGCTGGGAGGGGACTGGCTGGCGGGCAGGGCCTGATTTGCAGCGGCGCGATAACGGC‑GGGaaaaaaggg > 2:18527/1‑135 (MQ=255) cacGTCATCAACACTCCGCTGGCATTCCTGATTAAGAGTGCTTTGAACACCAATATCTTTGGTGAGCCAGGCTGGCAGGGTACTGGCTGGCGTGCAGGTCGTGATTTGCAGCGTCGCGATATCGGCGGGGAAAAACgggg > 2:255045/1‑139 (MQ=255) tcatcaACACTCCGCTGGCATTCCTGATTAAGAGTGCTTTGAACACCAATATCTTTGGTGAGCCAGGCTGGCAGGGTACTGGCTGGCGTGCAGGTCGTGATTTGCAGCGTCGCGATATCGGC‑GGGAAAACCGGGACCACt < 2:275765/140‑1 (MQ=255) tcatcaACACTCCGCTGGCATTCCTGATTAAGAGTGCTTTGAACACCAATATCTTTGGTGAGCCAGGCTGGCAGGGTACTGGCTGGCGTGCAGGTCGTGATTTGCAGCGTCGCGATATCGGC‑GGGAAAACCGGGACCACt < 2:521308/140‑1 (MQ=255) tcaACACTCCGCTGGCATTCCTGATTAAGAGTGCTTTGAACACCAATATCTTTGGTGAGCCAGGCTGGCAGGGTACTGGCTGGCGTGCAGGTCGTGATTTGCAGCGTCGCGATATCGGC‑GGGAAAACCGGGACCACTAAc < 1:378006/140‑1 (MQ=255) gCATTCCTGATTAAGAGTGCTTTGAACACCAATATCTTTGGTGAGCCAGGCTGGCAGGGTACTGGCTGGCGTGCAGGTCGTGATTTGCAGCGTCGCGATATCGGC‑GGGAAAACCGGGACCACTAACAGTTCGAAAGATGc < 1:711703/140‑1 (MQ=255) gagacaggATTAAGAGTGCTTTGAACACCAATATCTTTGGTGAGCCAGGCTGGCAGGGTACTGGCTGGCGTGCAGGTCGTGATTTGCAGCGTCGCGATATCGGC‑GGGAAAACCGGGACCACTAACAGTTCGAAAGATGCg < 1:457964/133‑1 (MQ=255) ccTGATTAAGAGTGCTTTGAACACCAATATCTTTGGTGAGCCAGGCTGGCAGGGTACTGGCTGGCGTGCAGGTCGTGATTTGCAGCGTCGCGATATCGGCGGGGAAAACCGGG‑CCACTAACAATTTGAAAAAAGCGTggg > 2:42709/1‑139 (MQ=255) agTGCTTTGAACACCAATATCTTTGGTGAGCCAGGCTGGCAGGGTACTGGCTGGCGTGCAGGTCGTGATTTGCAGCGTCGCGATATCGGC‑GGGAAAACCGGGACCACTAACAGTTCGAAAGATGCGTGGTTCTCGGGTTa < 1:42709/140‑1 (MQ=255) cTTTGAACACCAATATCTTTGGTGAGCCAGGCTGGCAGGGTACTGGCTGGCGTGCAGGTCGTGATTTGCAGCGTCGCGATATCGGCGGGGAAAACCGGG‑CCCCTAAAAGTTCGAAAGAAGCGTGGGTCTCGGGGTACggg > 1:512427/1‑139 (MQ=255) ggTGAGCCAGGCTGGCAGGGTACTGGCTGGCGTGCAGGTCGTGATTTGCAGCGTCGCGATATCGGCGGGGAAAACCGGG‑CCACCAAAAGTTTGGAAAAAGGGGGGGTCTCGGGGTTCGGGCCGGGGGGTGTGGCCTccgg > 2:737770/1‑137 (MQ=255) ggTGAGCCAGGCTGGCAGGGTACTGGCTGGCGTGCAGGTCGTGATTTGCAGCGTCGCGATATCGGC‑GGGAAAACCGGGACCACTAACAGTTCGAAAGATGCGTGGTTCTCGGGTTACGGTCCGGGCGTTGTGACCTCGGt < 1:370471/140‑1 (MQ=255) aGCCAGGCTGGCAGGGTACTGGCTGGCGTGCAGGTCGTGATTTGCAGCGTCGCGATATCGGC‑GGGGAAAACGGGGGCCCCAACAAATTGAAAAAAGCGGGGGTTTCGGGGTACGGGCCGGGGGGTGTGGCCTCGGGCttg > 2:519161/1‑138 (MQ=255) acagggCAGGGTACTGGCTGGCGTGCAGGTCGTGATTTGCAGCGTCGCGATATCGGC‑GGGAAAACCGGGACCACTAACAGTTCGAAAGATGCGTGGTTCTCGGGTTACGGTCCGGGCGTTGTGACCTCGGTCTGGATTgg < 1:353305/136‑1 (MQ=255) cTGGCAGGGTACTGGCTGGCGTGCAGGTCGTGATTTGCAGCGTCGCGATATCGGC‑GGGAAAACCGGGACCACTAACAGTTCGAAAGATGCGTGGTTCTCGGGTTACGGTCCGGGCGTTGTGACCTCGGTCTGGATTGGCt < 1:884074/140‑1 (MQ=255) cTGGCAGGGTACTGGCTGGCGTGCAGGTCGTGATTTGCAGCGTCGCGATATCGGC‑GGGAAAACCGGGACCACTAACAGTTCGAAAGATGCGTGGTTCTCGGGTTACGGTCCGGGCGTTGTGACCTCGGTCTGGATTGGCt < 1:737770/140‑1 (MQ=255) gCAGGGTACTGGCTGGCGTGCAGGTCGTGATTTGCAGCGTCGCGATATCGGC‑GGGGAAAACCGGACCACCCACAACTTGCAAGAAGGGGGGGTTTCGGGGTTCCGGCCCGGCGGTGTGGCCCCCGGCTGGGTTGGGTTtg > 2:242910/1‑140 (MQ=255) aCTGGCTGGCGTGCAGGTCGTGATTTGCAGCGTCGCGATATCGGC‑GGGTAAACCGGGACCACTAACAGTTCGAAAGATGCGTGGTTCTCGGGTTACGGTCCGGGCGTTGTGACCTCGGTCTGGATTGGCTTTGATGATCa < 1:822938/140‑1 (MQ=255) gTGCAGGTCGTGATTTGCAGCGTCGCGATATCGGC‑GGGAAAACCGGGACCACTAACAGTTCGAAAGATGCGTGGTTCTCGGGTTACGGTCCGGGCGTTGTGACCTCGGTCTGGATTGGCTTTGATGATCACCGTCGTAAt < 1:295698/140‑1 (MQ=255) gCAGGTCGTGATTTGCAGCGTCGCGATATCGGC‑GGGAAAACCGGGACCACTAACAGTTCGAAAGATGCGTGGTTCTCGGGTTACGGTCCGGGCGTTGTGACCTCGGTCTGGATTGGCTTTGATGATCACCGTCGTAAtct < 1:519161/140‑1 (MQ=255) cGTGATTTGCAGCGTCGCGATATCGGC‑GGGAAAACCGGGACCACTAACAGTTCGAAAGATGCGTGGTTCTCGGGTTACGGTCCGGGCGTTGTGACCTCGGTCTGGATTGGCTTTGATGATCACCGTCGTAATCTCGGTCa < 2:512427/140‑1 (MQ=255) | ACCGCACGTCATCAACACTCCGCTGGCATTCCTGATTAAGAGTGCTTTGAACACCAATATCTTTGGTGAGCCAGGCTGGCAGGGTACTGGCTGGCGTGCAGGTCGTGATTTGCAGCGTCGCGATATCGGC‑GGGAAAACCGGGACCACTAACAGTTCGAAAGATGCGTGGTTCTCGGGTTACGGTCCGGGCGTTGTGACCTCGGTCTGGATTGGCTTTGATGATCACCGTCGTAATCTCGGTCA > NC_000913/3524883‑3525125 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |