Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 289,256 | 0 | T | C | 54.5% | ‑1.2 / 28.9 | 22 | D125G (GAC→GGC) | yagL | recombinase family protein |
Reads supporting (aligned to +/- strand): ref base T (4/6); new base C (3/9); total (7/15) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 6.52e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
Rejected as consensus: E-value score below prediction cutoff. | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
ATGTCAATCGTATCCTGACGAGCAATCAAAAAAGACAAGCAACTTTCAGCAGAAGGCTCTATCCTTGAAGCTGGATGGTAGCTGTAGATAAACATTTCTTTTTGAAGGCAAAAAAATAAAATATTCTGTAAGTCCTCCACAGTTCTGCCAAGACAGCTTAGATCTGGAATGACAAGCGTGTCCTGCCGGTTCATTTCATGATGAATCAATTGTCTGAGTAAACGTTGCTCTGAAACAAC > NZ_CP009273/289124‑289362 | aTGTCAATCGTATCCTGACGAGCAATCAAAAAAGACAAGCAACTTTCAGCAGAAGGCTCTATCCTTGAAGCTGGATGGTAGCTGTAGATAAACATTTCTTTTTGAAGGCAAAAAAATAAAATATTCTGTAAGTCCTCca < 2:217459/139‑1 (MQ=255) gTATCCTGACGAGCAATCAAAAAAGACAAGCAACTTTCAGCAGAAGGCTCTATCCTTGAAGCTGGATGGTAGCTGTAGATAAACATTTCTTTTTGAAGGCAAAAAAATAAAATATTCTGTAAGCCCTCCACAGTTCTGc > 1:79606/1‑139 (MQ=255) cTGACGAGCAATCAAAAAAGACAAGCAACTTTCAGCAGAAGGCTCTATCCTTGAAGCTGGATGGTAGCTGTAGATAAACATTTCTTTTTGAAGGCAAAAAAATAAAATATTCTGTAAGCCCTCCACAGTTCTGCCAAGa < 2:79606/139‑1 (MQ=255) cTGACGAGCAATCAAAAAAGACAAGCAACTTTCAGCAGAAGGCTCTATCCTTGAAGCTGGATGGTAGCTGTAGATAAACATTTCTTTTTGAAGGCAAAAAAATAAAATATTCTGTAAGCCCTCCACAGTTCTGCCAAGa < 1:356837/139‑1 (MQ=255) aCGAGCAATCAAAAAAGACAAGCAACTTTCAGCAGAAGGCTCTATCCTTGAAGCTGGATGGTAGCTGTAGATAAACATTTCTTTTTGAAGGCAAAAAAATAAAATATTCTGTAAGTCCTCCACAGTTCTGCCAAGACAg > 1:222296/1‑139 (MQ=255) aCGAGCAATCAAAAAAGACAAGCAACTTTCAGCAGAAGGCTCTATCCTTGAAGCTGGATGGTAGCTGTAGATAAACATTTCTTTTTGAAGGCAAAAAAATAAAATATTCTGTAAGCCCTCCACAGTTCTGCCAAGACAg > 1:486300/1‑139 (MQ=255) aTCAAAAAAGACAAGCAACTTTCAGCAGAAGGCTCTATCCTTGAAGCTGGATGGTAGCTGTAGATAAACATTTCTTTTTGAAGGCAAAAAAATAAAATATTCTGTAAGTCCTCCACAGTTCTGCCAAGACAGCt > 2:96419/1‑134 (MQ=255) aaGACAAGCAACTTTCAGCAGAAGGCTCTATCCTTGAAGCTGGATGGTAGCTGTAGATAAACATTTCTTTTTGAAGGCAAAAAAATAAAATATTCTGTAAGTCCTCCACAGTTCTGCCAAGACAGCTTAGATCTGGAAt < 2:222296/139‑1 (MQ=255) aGCAACTTTCAGCAGAAGGCTCTATCCTTGAAGCTGGATGGTAGCTGTAGATAAACATTTCTTTTTGAAGGCAAAAAAATAAAATATTCTGTAAGTCCTCCACAGTTCTGCCAAGACAGCTTAGATCTGGAATGACAAg < 2:279021/139‑1 (MQ=255) gAAGGCTCTATCCTTGAAGCTGGATGGTAGCTGTAGATAAACATTTCTTTTTGAAGGCAAAAAAATAAAATATTCTGTAAGCCCTCCACAGTTCTGCCAAGACAGCTTAGATCTGGAATGACAAGCGTGTCCTGCCGGt < 1:421473/139‑1 (MQ=255) cTTGAAGCTGGATGGTAGCTGTAGATAAACATTTCTTTTTGAAGGCAAAAAAATAAAATATTCTGTAAGTCCTCCACAGTTCTGCCAAGACAGCTTAGATCTGGAATGACAAGCGTGTCCTGCCGGTTCATTTCatga < 1:96419/138‑1 (MQ=255) cTGGATGGTAGCTGTAGATAAACATTTCTTTTTGAAGGCAAAAAAATAAAATATTCTGTAAGCCCTCCACAGTTCTGCCAAGACAGCTTAGATCTGGAATGACAAGCGTGTCCTGCCGGTTCATTTCATGATGAATCaa < 2:264834/139‑1 (MQ=255) aTGGTAGCTGTAGATAAACATTTCTTTTTGAAGGCAAAAAAATAAAATATTCTGTAAGCCCTCCACAGTTCTGCCAAGACAGCTTAGATCTGGAATGACAAGCGTGTCCTGCCGGTTCATTTCATGATGAATCAATTGt < 1:158496/139‑1 (MQ=255) aaCATTTCTTTTTGAAGGCAAAAAAATAAAATATTCTGTAAGCCCTCCACAGTTCTGCCAAGACAGCTTAGATCTGGAATGACAAGCGTGTCCTGCCGGTTCATTTCatgatg < 1:349724/113‑1 (MQ=255) aaCATTTCTTTTTGAAGGCAAAAAAATAAAATATTCTGTAAGCCCTCCACAGTTCTGCCAAGACAGCTTAGATCTGGAATGACAAGCGTGTCCTGCCGGTTCATTTCatgatg > 2:349724/1‑113 (MQ=255) aaCATTTCTTTTTGAAGGCAAAAAAATAAAATATTCTGTAAGCCCTCCACAGTTCTGCCAAGACAGCTTAGATCTGGAATGACAAGCGTGTCCTGCCGGTTCATTTCATGATGAATCAATTGTCTGAGTAAACGTTGct < 2:160247/139‑1 (MQ=255) cATTTCTTTTTGAAGGCAAAAAAATAAAATATTCTGTAAGTCCTCCACAGTTCTGCCAAGACAGCTTAGATCTGGAATGACAAGCGTGTCCTGc > 1:370559/1‑94 (MQ=255) cATTTCTTTTTGAAGGCAAAAAAATAAAATATTCTGTAAGTCCTCCACAGTTCTGCCAAGACAGCTTAGATCTGGAATGACAAGCGTGTCCTGc < 2:370559/94‑1 (MQ=255) tttCTTTTTGAAGGCAAAAAAATAAAATATTCTGTAAGTCCTCCACAGTTCTGCCAAGACAGCTTAGATCTGGAATGACAAGCGTGTCCTGCCGGTTCATTTCATGATGAATCAATTGTCTGAGTAAACGTTGCTCTg > 1:282310/1‑138 (MQ=255) ttGAAGGCAAAAAAATAAAATATTCTGTAAGTCCTCCACAGTTCTGCCAAGACAGCTTAGATCTGGAATGACAAGCGTGTCCTGCCGGTTCATTTCATGATGAATCAATTGTCTGAGTAAACGTTGCTCTGAaacaac < 1:307307/138‑1 (MQ=255) ttGAAGGCAAAAAAATAAAATATTCTGTAAGCCCTCCACAGTTCTGCCAAGACAGCTTAGATCTGGAATGACAAGCGTGTCCTGCCGGTTCATTTCATGATGAATCAATTGTCTGAGTAAACGTTGCTCTGAaacaac < 1:458590/138‑1 (MQ=255) ttGAAGGCAAAAAAATAAAATATTCTGTAAGCCCTCCACAGTTCTGCCAAGACAGCTTAGATCTGGAATGACAAGCGTGTCCTGCCGGTTCATTTCATGATGAATCAATTGTCTGAGTAAACGTTGCTCTGAaacaac < 2:486300/138‑1 (MQ=255) | ATGTCAATCGTATCCTGACGAGCAATCAAAAAAGACAAGCAACTTTCAGCAGAAGGCTCTATCCTTGAAGCTGGATGGTAGCTGTAGATAAACATTTCTTTTTGAAGGCAAAAAAATAAAATATTCTGTAAGTCCTCCACAGTTCTGCCAAGACAGCTTAGATCTGGAATGACAAGCGTGTCCTGCCGGTTCATTTCATGATGAATCAATTGTCTGAGTAAACGTTGCTCTGAAACAAC > NZ_CP009273/289124‑289362 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |