Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 3,556,207 | 0 | T | C | 30.8% | 61.3 / 14.1 | 39 | S215G (AGT→GGT) ‡ | rtcA | RNA 3'‑terminal phosphate cyclase |
Reads supporting (aligned to +/- strand): ref base T (25/2); new base C (0/12); total (25/14) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 2.33e-08 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
GCACTGACGCGCTTTTCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAA‑A‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGCGAGACTTCC > NC_000913/3556074‑3556343 | gCACTGACGCTTTTTTCACCGACGAAAAAAAAGCGTTCGGTGATTTTTTCACTTTCGACTTCAAGCGAGACGGTATTCCCCGGCCCCTGGCCGCGCGGCAGGTTATAAATATTCTGTTCATGCAGGAAAAA‑C‑CCCCCCGc < 2:1663911/140‑1 (MQ=255) cTGACGCGCTTTTCACCGACGACAAAAAAGCGTCCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAA‑A‑CTCCCCGCCAg > 2:906107/1‑140 (MQ=255) aCGCGCTTTTCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAA‑A‑CTCCCCGCCAgtgt > 2:157367/1‑140 (MQ=255) aCGCGCTTTTCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAA‑A‑CTCCCCGCCAgtgt > 1:1886633/1‑140 (MQ=255) cACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAA‑A‑CTCCCCGCCAGTGTAGCGATTTCa > 1:754345/1‑140 (MQ=255) ccGACGAAAAAAACGCGTTGGGTGATTTTTTCTCTTTCGACTTAAAGCGAGCCGTTATTCCCGGGCCCGGGGTCGCGCGGCAGTTTAGAAATTTTTTTTTCATGCAGGGAAAA‑C‑CCCCCCGCCAGTGTAGCGATTTCACg < 2:1100678/140‑1 (MQ=255) ccGACGAAAAAAAAGCGTTGGGTGATTTTTTCCCTTTCGACTTCAAGCGAGACGGTATTCCCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGAAAAA‑C‑CCCCCCGCCAGTGTAGCGATTTCACg < 2:754345/140‑1 (MQ=255) ccGACAAAAAAAAGGCGTGGGGTTATTTTTTCTTTTTCGTTTTAAAGCGAGCCGTTATTCCCCGCCCCCGGGCCGCGGGGGGGGTTATAAATTTTCTGTTCATGCAGGAAAAA‑C‑CCCCCCGCCAGTGTAGCGATTTCACg < 2:433010/140‑1 (MQ=255) aaGCTTTCGGTGATATTTTCTTTTCCAAATTCAAGGGAGACGGTATCCCCCGGCCCCTGGCCGCGGGGCAGGTTATAAATTTTCTGTTCATGCAGGAAAAA‑T‑CCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATg < 2:478092/140‑1 (MQ=255) ctttCGGTGATTTTTTCTTTTTCGACTTAAAGCGAGAGGGTATTCCCCGCCCCCTGGTCGCGCGGCAGTTTATAAATTTTCTTTTCATGAAGGAAAAA‑C‑CCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGgcg < 2:2090108/138‑1 (MQ=255) ttgggTGATTTTTTCTCTTTCGACTTCAAGCGAGACGGTATTCCCCGGCCCCTGGTCGCGCGGCAGTTTATAAATTTTCTGTTCATGCAGGAAAAA‑C‑CCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGcgcg < 2:1260734/137‑1 (MQ=255) ttCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAA‑A‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGcgcg > 1:2090108/1‑140 (MQ=255) tttttttttCCTTTTCGCTTTAAAGGGAGACGGTTTTCCCCGGCCCCTGGTCGCGCGGGGGGTTATAAATTTTCTGTTCATGCGGGAAAAA‑C‑CCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGcaca < 2:1481734/135‑1 (MQ=255) ttattttttCTTTTTCGACTTCAAGGGAGACGGTATTCCCCGGCCCGTGGTCGCGCGGCAGGTTATGAATTTTCTGTTCAGGCAGGAAAAA‑C‑CCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGcaca < 2:1886633/135‑1 (MQ=255) aCTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAA‑A‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATa > 1:1481734/1‑140 (MQ=255) aCTTTCGACTTCAAGCGAGACGGGATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAA‑A‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATa > 2:1175562/1‑140 (MQ=255) tttCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAA‑A‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGa > 1:1953416/1‑140 (MQ=255) gccTTAAAGCGAGACGGTTTTCCCCGGCCCCGGGCCGCGCGGCAGTTTATAAATTTTCTGTTCATGAAGGAAAAA‑C‑CCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACtt < 2:1953416/138‑1 (MQ=255) aaaGCGAGCCGGTATCCCCCGGCCCCGGGCCGCGCGGCAGTTTATGAATTTTCTTTTCATGAAGGAAAAA‑A‑CCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCa < 2:2186867/139‑1 (MQ=255) gggaCGTTATCCCCCGGCCCCTGGCCGCGCGGAGGTTTATAAATTTTCTGTTCATGAAGGAAAAA‑C‑CCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCAt < 2:1765754/138‑1 (MQ=255) gagaCGTTATTCCCCGGCCCCTGGTCGCGCGGCAGGTTATGAATTTTCTGTTCATGCAGGAAAAA‑C‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCAt < 1:781248/140‑1 (MQ=255) gaCGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAA‑A‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCt > 1:841846/1‑140 (MQ=255) tATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAA‑A‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGcaca > 1:1011598/1‑140 (MQ=255) ccccTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAA‑A‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCcgc > 2:581725/1‑140 (MQ=255) ccccTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAA‑A‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCcgc > 2:1320941/1‑140 (MQ=255) ccgcgcgGCAGGTTATAAATATTCTGTTCATGCAGGAAAAA‑C‑CCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGcc < 1:157367/139‑1 (MQ=255) cgcgGCAGGTTATGAATATTCTGTTCATGCAGGGAAAA‑A‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAg > 1:972162/1‑140 (MQ=255) cgcgGCAGGTTATGAATATTCTGTTCATGCAGGGAAAA‑A‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAg > 1:2290065/1‑140 (MQ=255) cgcgGCAGGTTATGAATATTCTGTTCATGCAGGGAAAA‑A‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAg > 1:1945858/1‑140 (MQ=255) cgcgGCAGGTTATGAATATTCTGTTCATGCAGGGAAAA‑A‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAg > 1:1727876/1‑140 (MQ=255) gCAGGTTATGAATATTCTGTTCATGCAGGGAAAA‑A‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGc > 1:1232932/1‑140 (MQ=255) gagttaaTGAATATCTTGTAAATG‑AGGGAAAATC‑CCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCa < 2:264705/134‑1 (MQ=255) gaGGTTATGAATTTTCTGTTCAGGCAGGAAAAA‑C‑CCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCa < 2:1463125/139‑1 (MQ=255) aTGAATATTCTGTTCATGCAGGAAAAA‑C‑CCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGgtgt < 2:1011598/140‑1 (MQ=255) ttCATGCAGGGAAAA‑A‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCa > 1:708974/1‑140 (MQ=255) ttCATGCAGGGAAAA‑A‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCa > 2:329866/1‑140 (MQ=255) ttCAGGGAGGAAAAA‑C‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCa < 2:418363/140‑1 (MQ=255) cATGCAGGGAAAA‑A‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACg > 1:1524510/1‑140 (MQ=255) gtggaGGGAAAA‑ACCTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACg < 2:1271307/136‑1 (MQ=255) aTGCAGGGAAAA‑A‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACgg > 1:205406/1‑140 (MQ=255) tGCAGGGAAAA‑A‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACggg > 1:2257562/1‑140 (MQ=255) aGGGAAAA‑A‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGCga > 2:903914/1‑140 (MQ=255) aGGAAAAA‑C‑CCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGCga < 2:972162/140‑1 (MQ=255) ggAAAA‑A‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGCgaga > 2:380999/1‑140 (MQ=255) a‑a‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGCGAGACTTcc < 1:329866/140‑1 (MQ=255) a‑a‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGCGAGACTTcc > 1:1841242/1‑140 (MQ=255) | GCACTGACGCGCTTTTCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAA‑A‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGCGAGACTTCC > NC_000913/3556074‑3556343 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 24 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |