Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 580,498 | 0 | G | A | 56.2% | 6.2 / 21.9 | 16 | V148I (GTT→ATT) | nohD | DLP12 prophage; DNA packaging protein |
Reads supporting (aligned to +/- strand): ref base G (5/2); new base A (1/8); total (6/10) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 3.50e-02 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.83e-01 | |||||||||||
Rejected as consensus: E-value score below prediction cutoff. | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
GGTGCAGCGGCGTTTTCCGGAACTGGAAAACCGACATGTTGATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAAC > NC_000913/580461‑580591 | ggagCAGCGGCGTTTTCCGGAACTGGAAAACCGACATGTTGATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAagca > 2:7263042/4‑82 (MQ=1) gttcAGCGGCGCTTCCCGGAGCTGGAAAACCGACATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAagcag < 2:5523032/79‑1 (MQ=0) gttcAGCGGCGCTTCCCGGAGCTGGAAAACCGACATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAagcag < 2:3488899/79‑1 (MQ=0) gttcAGCGGCGCTTCCCGGAGCTGGAAAACCGACATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAagcag < 2:1990699/79‑1 (MQ=0) tgagcTGGAAAACCGACATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTg < 1:4582667/90‑1 (MQ=1) gatcTGGAAAACCGATATGTTGATTTCCTGAAACGGGATATCATAAAAGCCATGAACAAAGCAGCCGCTCTGGATGAACTGa < 2:766530/79‑1 (MQ=1) tGGAAAACCGACATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATAcc < 2:3200725/82‑1 (MQ=0) tGGAAAACCGACATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATAcc < 2:279524/82‑1 (MQ=0) gACATGTTGATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAatatat < 1:4228662/94‑1 (MQ=31) gACATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAGCAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAatatat > 1:5250851/1‑94 (MQ=2) gcATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATc < 1:6186727/93‑1 (MQ=11) ccATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATc < 1:3111502/93‑1 (MQ=11) aCATGTTGATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAatatatg > 1:1261932/1‑93 (MQ=14) aCATGTTGATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAatatatg > 1:2438108/1‑93 (MQ=14) cATGTTGATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATct > 1:6440450/1‑93 (MQ=17) tGTTGATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGac > 1:6409606/1‑93 (MQ=21) tGTTGATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGaa < 1:1782555/94‑1 (MQ=35) gTTGATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAAc > 1:4103166/1‑94 (MQ=37) gTTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCGAAc < 1:5327094/94‑1 (MQ=11) | GGTGCAGCGGCGTTTTCCGGAACTGGAAAACCGACATGTTGATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAAC > NC_000913/580461‑580591 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 23 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |