Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 3,556,207 | 0 | T | C | 35.7% | 49.1 / 23.2 | 42 | S215G (AGT→GGT) ‡ | rtcA | RNA 3'‑terminal phosphate cyclase |
Reads supporting (aligned to +/- strand): ref base T (25/2); new base C (0/15); total (25/17) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.38e-09 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
TCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGCGAG > NC_000913/3556071‑3556337 | ccGGCATTGACGCGTTTTCCCCCGACGAAAAAAAGGCTTTCGGTGATTTTTTCCCTTTCGACTACAAGCGAGACGGTTTTCCCCGCCCCCGGGTCGCGGGGCAGTTTATGAATATTCTGTTCATGAAGGAAAAAC‑CTccc < 2:878662/139‑1 (MQ=255) cTTTTCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGa > 1:1172236/1‑140 (MQ=255) ttCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGAttt > 2:1116700/1‑140 (MQ=255) ccGACAAAAAAAAAGCGTGCGTTATTTTTTTCTTTTTCGACTTGAGGGGGGCGGTTATTCCCCGCCCCCTGGCGGCGGGGCAGTTTATAAATTTTCTGTTCAGGCGGGAAAAAC‑CCCCCCGCCAGTGTAGCGATTTCACg < 2:653124/140‑1 (MQ=255) gacgacAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCt > 1:749707/1‑140 (MQ=255) aggCGTTCGGTTATTTTTTCACTTCCGCCTTCAAGCGAGACGGTATTCCCCGGCCCCTGGTCGGGCGGCAGTTTATGAATTTTCTGTTCATGCAGGAAAAAC‑CCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATg < 2:664768/138‑1 (MQ=255) aaGTGTTAGGTGATTTTTTCTTTTTCGGCTTCAGGAGGGAGGGTTTTCCCCGGCCCCTGGGGGCGCGGCGGGTTAGGAATTTTCTTTTCATG‑AGGGAAAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATg < 2:749707/140‑1 (MQ=255) ttCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGcgcg > 1:983276/1‑140 (MQ=255) ggtatttttttCCCTTCCGATTTAAAGCGAGACGTTTTCCCCCGCCCCCTGGTCGCGGGGCAGTTTATAAATTTTCTTTTCATGAAGGAAAAAC‑CCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGca < 2:953177/133‑1 (MQ=255) ggTGATTTTTCCCCTTCCGACTTCAAGCGAGAGGTTATTCCCCGGCCCCTGGCCGCGCGGCAGTTTATGAATATTCTGTCCATGCAGGAAAAAC‑CCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGca < 2:1197880/140‑1 (MQ=255) tGATTTTTCCACTTTCGACTTCAAGCGAGACGGTATTCCCCGGCCCCTGGCCGCGCGGCAGGTTATAAATATTCTGTTCATGCAGGAAAAAC‑CCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGcaca < 2:980702/140‑1 (MQ=255) tGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGcaca > 1:430131/1‑140 (MQ=255) tGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGcaca > 1:1394456/1‑140 (MQ=255) gATTTTTTCACTTTCGACTTAAAGCGAGACGGTATCCCCCGGCCCCTGGCCGCGCGGCAGTTTATAAATTTTCTGTTCATGCGGGAAAAAC‑CCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCacac < 2:1394456/140‑1 (MQ=255) aattttttCACTTTCGATTTAAAGCGAGACGGTATTCCCCGGCCCCGGGTCGCGCGGCAGTTTATAAATATTCTGTTCATGAAGGGAAAAC‑CCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCacac < 2:430131/136‑1 (MQ=255) tatTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCa > 1:736556/1‑140 (MQ=255) tCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTaa > 1:675282/1‑140 (MQ=255) aCTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATa > 1:1009411/1‑140 (MQ=255) aCTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATa > 1:606042/1‑140 (MQ=255) aCTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATa > 1:456139/1‑140 (MQ=255) aCTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATa > 2:127274/1‑140 (MQ=255) cTTCAAGCGAGCCGTTATCCCCCGCCCCCTGGTCGCGCGGAAGTTTATAAATATTCTGTTCAGGCAGGAAAAAC‑CCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTtct < 2:1172236/140‑1 (MQ=255) aaaGGGAGACGTTATTCCCCGCCCCCGGGTCGCGGGGCAGTTTATAAATTTTTTTTTCAGGCAGGAAAAAC‑CCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCa < 2:675282/139‑1 (MQ=255) gaCGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCt > 1:669461/1‑140 (MQ=255) cggttATTCCCCGGCCCCTGGTCGCGCGGCAGGTTATAAATTTTCTGTTCAGGCAGGAAAAAC‑CCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTg < 2:983276/136‑1 (MQ=255) ggTATTCCCCGCCCCCTGGCCGCGCGGCAGTTTATAAATTTTCTGTTCATGCAGGAAAAAC‑CCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGca < 2:456139/140‑1 (MQ=255) cccGCCCCCTGGTCGCGCGGAAGGTGATAAATATTCTGTTCAGGAGGGAAAAAA‑CCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGtt < 1:1296158/140‑1 (MQ=255) ccccTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCcgc > 2:1130736/1‑140 (MQ=255) ccccTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCcgc > 1:1143698/1‑140 (MQ=255) ccccTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCcgc > 1:294691/1‑140 (MQ=255) cgcgGCGGTTTATAAATATTCTGTTCAGGAAGGAAAAAC‑CCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAg < 1:35201/140‑1 (MQ=255) cgcgGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAg > 2:257419/1‑140 (MQ=255) ggCAGGTTATGAATATTCTGTTCATGCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTg > 1:176054/1‑140 (MQ=255) gCAGTTTATAAATTTTCTTTTAATGAAGGAAAAAC‑CCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGc < 2:475375/140‑1 (MQ=255) ggTTATAAATATTCTGTTCATGAAGGAAAAAC‑CCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAg < 2:1188039/140‑1 (MQ=255) ttCTGTTCATGCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGa > 1:1188039/1‑140 (MQ=255) tCTGTTCATGCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGAt > 2:1367289/1‑140 (MQ=255) gTTCATGCAGGGAAAAC‑CCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCaa < 1:1083283/124‑1 (MQ=255) gTTCATGCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCaa > 2:1083283/1‑124 (MQ=255) gTTCATGAAGGAAAAAC‑CCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGcc < 2:1143698/140‑1 (MQ=255) ttCATGCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCa > 1:20914/1‑140 (MQ=255) aTGCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACgg > 1:794271/1‑140 (MQ=255) gaaGGAAAAAC‑CCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGc < 2:20914/138‑1 (MQ=255) gCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGc > 2:627970/1‑140 (MQ=255) aGGAAAAAC‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGCga < 2:197992/140‑1 (MQ=255) gggAAAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGCga < 2:669461/140‑1 (MQ=255) gggAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGCgag > 1:780072/1‑140 (MQ=255) | TCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGCGAG > NC_000913/3556071‑3556337 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |