Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 2,715,349 | 0 | C | A | 66.7% | 33.8 / 40.4 | 39 | intergenic (‑29/‑74) | yfiE/eamB | putative DNA‑binding transcriptional regulator/cysteine and O‑acetylserine exporter |
Reads supporting (aligned to +/- strand): ref base C (5/8); new base A (12/14); total (17/22) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 7.40e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.70e-02 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
AATATGAAAAGTCACCGTCGATTGTGTACAGCACAATTTTTGCGAAGCTCGCAAAAAGGAACCCTCTTCCACCACGGTTTTAAGCGTAATAAAACGGCGCAGATCCATAACCCCAAACCTATCGAAAATATCGAATCTAGAATATAAAAACATTCATTTTTTTAAATGTTCCGTGTCGGGTACTGTCTACCAAAACAGAGGAGATAACAAGTGACACCGACCCTTTTAAGT > NC_000913/2715213‑2715443 | aaTATGAAAAGTCACCGTCGATTGTGTACAGCACAATTTTTGCGAAGCTCGCAAAAAGGAACCCTCTTCCACCACGGTTTTAAGCGTAATAAAACGGCGCAGATCCATAACCCCAAACCTATCGAAAATATCGAATCTAGAATATaaa > 3:123622/1‑148 (MQ=255) aGTCACCGTCGATTGTGTACAGCACAATTTTTGCGAAGCTCGCAAAAAGGAACCCTCTTCCACCACGGTTTTAAGCGTAATAATACGGCGCAGATCCATAACCCCAAACCTATCGAAAATATCGAATATAGAATATAAAAAcattcatt < 6:93625/149‑1 (MQ=255) gtgtACAGCACAATTTTTGCGAAGCTCGCAAAAAGGAACCCTCTTCCACCACGGTTTTAAGCGTAATAAAACGGCGCAGATCCATAACCCCAAACCTATCGAAAATATCGAATCTAGAATATAAAAACATTCATTTTTTTAAATGTTc < 4:40207/148‑1 (MQ=255) tgtACAGCACAATTTTTGCGAAGCTCGCAAAAAGGAACCCTCTTCCACCACGGTTTTAAGCGTAATAAAACGGCGCAGATCCATAACCCCAAACCTATCGAAAATATCGAATCTAGAATATAAAAACATTCATTTTTTTAAATGTTcc < 4:10869/148‑1 (MQ=255) gCACAATTTTTGCGAAGCTCGCAAAAAGGAACCCTCTTCCACCACGGTTTTAAGCGTAATAAAACGGCGCAGATCCATAACCCCAAACCTATCGAAAATATCGAATCTAGAATATAAAAACATTCATTTTTTTAAATGTTCCGTGTCgg < 4:123622/149‑1 (MQ=255) tttGCGAAGCTCGCAAAAAGGAACCCTCTTCCACCACGGTTTTAAGCGTAATAAAACGGCGCAGATCCATAACCCCAAACCTATCGAAAATATCGAATCTAGAATATAAAAACATTCATTTTTTTAAATGTTCCGTGTCGGGTACTGt < 5:66646/148‑1 (MQ=255) cgagcttCGCAAAAAGGAACCCTCTTCCACCACGGTTTTAAGCGTAATAAAACGGCGCAGATCCATAACCCCAAACCTATCGAAAATATCGAATATAGAATATAAAAACATTCATTTTTTTAAATGTTCCGTGTCGGGTACTGTc < 4:50601/139‑1 (MQ=255) cgagcttCGCAAAAAGGAACCCTCTTCCACCACGGTTTTAAGCGTAATAAAACGGCGCAGATCCATAACCCCAAACCTATCGAAAATATCGAATATAGAATATAAAAACATTCATTTTTTTAAATGTTCCGTGTCGGGTACTGTc > 3:50601/7‑145 (MQ=255) cAAAAAGGAACCCTCTTCCACCACGGTTTTAAGCGTAATAAAACGGCGCAGATCCATAACCCCAAACCTATCGAAAATATCGAATATAGAATATAAAAACATTCATTTTTTTAAATGTTCCGTGTCGGGTACTGTCTACCAAAACagag > 8:136011/1‑149 (MQ=255) aaCCCTCTTCCACCACGGTTTTAAGCGTAATAAAACGGCGCAGATCCATAACCCCAAACCTATCGAAAATATCGAATATAGAATATAAAAACATTCATTTTTTTAAATGTTCCGTGTCgg > 8:59039/1‑120 (MQ=255) aaCCCTCTTCCACCACGGTTTTAAGCGTAATAAAACGGCGCAGATCCATAACCCCAAACCTATCGAAAATATCGAATATAGAATATAAAAACATTCATTTTTTTAAATGTTCCGTGTCgg < 7:59039/120‑1 (MQ=255) tctTCCACCACGGTTTTAAGCGTAATAAAACGGCGCAGATCCATAACCCCAAACCTATCGAAAATATCGAATATAGAATATAAAAACATTCATTTTTTTAAATGTTCCGTGTCGGGTACTGTCTACCAAAACAGAGGAGATAACAAGt > 2:23382/1‑148 (MQ=255) accacGGTTTTAAGCGTAATAAAACGGCGCAGATCCATAACCCCAAACCTATCGAAAATATCGAATCTAGAATATAAAAACATTCATTTTTTTaaa > 1:110441/1‑96 (MQ=255) accacGGTTTTAAGCGTAATAAAACGGCGCAGATCCATAACCCCAAACCTATCGAAAATATCGAATCTAGAATATAAAAACATTCATTTTTTTaaa < 2:110441/96‑1 (MQ=255) accacGGTTTTAAGCGTAATAAAACGGCGCAGATCCATAACCCCAAACCTATCGAAAATATCGAATCTAGAATATAAAAACATTCATTTTTTTAAATGTTCCGTGTCGGGTACTGTCTACCAAAACAGAGGAGATAACAAGTGACACCg > 7:29176/1‑149 (MQ=255) accacGGTTTTAAGCGTAATAAAACGGCGCAGATCCATAACCCCAAACCTATCGAAAATATCGAATATAGAATATAAAAACATTCATTTTTTTAAATGTTCCGTGTCGGGTACTGTCTACCAAAACagag > 2:40081/1‑130 (MQ=255) accacGGTTTTAAGCGTAATAAAACGGCGCAGATCCATAACCCCAAACCTATCGAAAATATCGAATATAGAATATAAAAACATTCATTTTTTTAAATGTTCCGTGTCGGGTACTGTCTACCAAAACagag < 1:40081/130‑1 (MQ=255) cacGGTTTTAAGCGTAATAAAACGGCGCAGATCCATAACCCCAAACCTATCGAAAATATCGAATATAGAATATAAAAACATTCATTTTTTTAAATGTTCCGTGTc < 8:55321/105‑1 (MQ=255) cacGGTTTTAAGCGTAATAAAACGGCGCAGATCCATAACCCCAAACCTATCGAAAATATCGAATATAGAATATAAAAACATTCATTTTTTTAAATGTTCCGTGTc > 7:55321/1‑105 (MQ=255) acGGTTTTAAGCGTAATAAAACGGCGCAGATCCATAACCCCAAACCTATCGAAAATATCGAATCTAGAATATAAAAACATTCATTTTTTTAAATGTTCCGTGTCGGGTACTGTCTACCAAAACAGAGGAGATAACAAGTGACACCGAcc > 1:147385/1‑149 (MQ=255) cGGTTTTAAGCGTAATAAAACGGAGCAGATCCATAACCCCAAACCTAACGAAAATATCGAATATAGAATATAAAAACATTCATTTTTTTAAATGTTCCGTGTCGGGTACTGTCTACCAAAACAGAGGAGATAACAAGTGACACCGAccc < 6:104685/149‑1 (MQ=255) ttttAATCGTAATAAAACGGCGCAGATCCATAAACCCAAACCTATCGAAAATATCGAATATAGAATATAAAAACATTCATTTTTTTAAATGTTCCGTGTCGGGTACTGTCTACCAAAACAGAGgaga < 2:10034/127‑1 (MQ=255) ttttAAGCGTAATAAAACGGCGCAGATCCATAACCCCAAACCTATAGAAAATATCGAATATAGAATATAAAAACATTCATTTTTTTAAATGTTCCGTGTCGGGTACTGTCTACCAAAACAGAGgaga > 1:10034/1‑127 (MQ=255) ttAAGCGTAATAAAACGGCGCAGATCCATAACCCCAAACCTATCGAAAATATCGAATCTAGAATATAAAAACATTCATTTTTTTAAATGTTCCGTGTCGGGTACTGTCTACCAAAACAGAGGAGATAACAAGTGACACCGACCCtttt < 2:147385/148‑1 (MQ=255) ttAAGCGTAATAAAACGGCGCAGATCCATAACCCCAAACCTATCGAAAATATCGAATATAGaa > 4:52480/1‑63 (MQ=255) ttAAGCGTAATAAAACGGCGCAGATCCATAACCCCAAACCTATCGAAAATATCGAATATAGaa < 3:52480/63‑1 (MQ=255) aaGCGTAATAAAACGGCGCAGATCCATAACCCCAAACCTATCGAAAATATCGAATCTAGAATATAAAAACATTCATTTTTTTAAATGTTCCGTGTCGGGTACTGTCTACCAAAACAGAGGAGATAACAAGTGACACCGACCCTTTTaa < 8:29176/148‑1 (MQ=255) aaGCGTAATAAAACGGCGCAGATCCATAACCCCAAACCTATCGAAAATATCGAATATAGAATATAAAAACATTCATTTTTTTAAATGTTCCGTGTCgg > 5:110817/1‑98 (MQ=255) aaGCGTAATAAAACGGCGCAGATCCATAACCCCAAACCTATCGAAAATATCGAATATAGAATATAAAAACATTCATTTTTTTAAATGTTCCGTGTCgg < 6:110817/98‑1 (MQ=255) aaGCGTAATAAAACGGCGCAGATCCATAACCCCAAACCTATCGAAAATATCGAATATAGAATATAAAAACATTCATTTTTTTAAATGTTCCGTGTCGGGTACTGTCTACCAAAACAGAGGAGATAACAAGTGACACCGACCCTTTTaa < 7:136011/148‑1 (MQ=255) aGCGTAATAAAACGGCGCAGATCCATAACCCCAAACCTATCGAAAATATCGAATATAGAATATAAAAACATTCATTTTTTTAAATGTTCCGTGTCGGGTACTGTCTACCAAAACAGAGGAGATAACAAGTGACACCGACCCTTTTAAGt < 1:23382/149‑1 (MQ=255) taataaAACGGCGCAGATCCATAACCCCAAACCTATCGAAAATATCGAATATAGAATATAAAAACATTCATTTTTTTAAATGTTCCGTGTCGGGTACTGTCTACCAAAACAGAGGAGATAACAAGt > 2:125007/1‑126 (MQ=255) taataaAACGGCGCAGATCCATAACCCCAAACCTATCGAAAATATCGAATATAGAATATAAAAACATTCATTTTTTTAAATGTTCCGTGTCGGGTACTGTCTACCAAAACAGAGGAGATAACAAGt < 1:125007/126‑1 (MQ=255) ataaAACGGCGCAGATCCATAACCCCAAACCTATCGAAAATATCGAATATAGAATATAAAAACATTCATTTTTTTAAATGTTCCgtg < 5:80053/87‑1 (MQ=255) ataaAACGGCGCAGATCCATAACCCCAAACCTATCGAAAATATCGAATATAGAATATAAAAACATTCATTTTTTTAAATGTTCCgtg > 6:80053/1‑87 (MQ=255) ggCGCAGATCCATAACCCCAAACCTATCGAAAATATCGAATCTAGAATATAAAAACATTCATTTTTTTAAATGTTCCGTGTCgg > 4:103781/1‑84 (MQ=255) ggCGCAGATCCATAACCCCAAACCTATCGAAAATATCGAATCTAGAATATAAAAACATTCATTTTTTTAAATGTTCCGTGTCgg < 3:103781/84‑1 (MQ=255) tCGAAAATATCGAATATAGAATATAAAAACATTCATTTTTTTAAATGTTCCGTGTCGGGTACTGTCTACCAAAACAGAGGAGATAACAAGTGACACCGACCCTTTTAAGt > 7:58022/1‑110 (MQ=255) tCGAAAATATCGAATATAGAATATAAAAACATTCATTTTTTTAAATGTTCCGTGTCGGGTACTGTCTACCAAAACAGAGGAGATAACAAGTGACACCGACCCTTTTAAGt < 8:58022/110‑1 (MQ=255) | AATATGAAAAGTCACCGTCGATTGTGTACAGCACAATTTTTGCGAAGCTCGCAAAAAGGAACCCTCTTCCACCACGGTTTTAAGCGTAATAAAACGGCGCAGATCCATAACCCCAAACCTATCGAAAATATCGAATCTAGAATATAAAAACATTCATTTTTTTAAATGTTCCGTGTCGGGTACTGTCTACCAAAACAGAGGAGATAACAAGTGACACCGACCCTTTTAAGT > NC_000913/2715213‑2715443 |
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: ‑ |