Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 3,294,154 | 0 | G | T | 54.5% | 3.7 / 14.6 | 11 | T161N (ACC→AAC) | ubiT | SCP2 domain‑containing protein |
Reads supporting (aligned to +/- strand): ref base G (5/0); new base T (0/6); total (5/6) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 2.16e-03 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.28e-03 | |||||||||||
Rejected as consensus: E-value score below prediction cutoff. | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
CCCGGCGCATCAATGGGAATTTCTACTCGAATTAGCATGGTTCACCTACCGATGTCTGTTTGGTTTCAGGCGCGGTTTTCATTCCCGCCTCAACAAAATCCGCCAGTTGCAGCAGCATCATGCGTAAGGCTTTCGGCATTTGCTCC > NZ_CP009273/3294080‑3294225 | cccGGCGCATCAATGGGATTTTCTCCTCGAATTAGCATGGTTCCCCTCCCGATGTCTGTTTGTTTTCGGGCGGGTTTTTCATCCCCGCCt < 2:389452/90‑1 (MQ=255) cGGCGCATCAATGGGAATTTCTACTCGAATTAGCATGGTTCACCTACCGATGTCTGTTTGGTTTCAGGCGCGGTTTTCATTCCCGCCTca > 1:294753/1‑90 (MQ=255) gcgcATCAATGGGAATTTCTACTCGAATTAGCATGGTTCACCTACCGATGTCTGTTTGGTTTCAGGCGCGGTTTTCATTCCCGCCTcaac > 1:255336/1‑90 (MQ=255) cTACCCAAATTAGCATGTTCCCCCTCCCGATGTTTTTTTGTTTTCGGGCGGGTTTTTCTTTCCCGCCTCAACAAAATCCGCCAGTTgcag < 2:434840/90‑1 (MQ=255) gTTCACCTACCGATGTCTGTTTGGTTTCAGGCGCGGTTTTCATTCCCGCCTCAACAAAATCCGCCAGTTGCAGCAGCATCATGCGTAAgg > 2:399555/1‑90 (MQ=255) tACCGATTTCTTTTTGGTTTCAGGCGCGTTTTTCATTCCCGCCTCAACAAAATCCGCCAGTTGCAGCAGCATCATGCGTAAGGCTTTCgg < 1:440807/90‑1 (MQ=255) cGATGTCTGTTTGGTTTCAGGCGCGGTTTTCATTCCCGCCTCAACAAAATCCGCCAGTTGCAGCAGCATCATGCGTAAGGCTTTCGGCAt > 1:44970/1‑90 (MQ=255) tttttGTTTTCAGGCGCGTTTTTCATCCCCGCCTCAACAAAATCCGCCAGTTGCAGCAGCATCATGCGTAAGGCTTTCGGCATTTGCTcc < 1:358770/88‑1 (MQ=255) tgtttgttTTCAGGCGGGTTTTTCATTCCCGCCTCAACAAAATCCGCCAGTTGCAGCAGCATCATGCGTAAGGCTTTCGGCATTTGCTcc < 1:372888/90‑1 (MQ=255) gcgTTTTTCATCCCCGCCTCAACAAAATCCGCCAGTTGCAGCAGCATCATGCGTAAGGCTTTCGGCATTTGCTcc < 2:421020/75‑1 (MQ=255) gcgGTTTTCATTCCCGCCTCAACAAAATCCGCCAGTTGCAGCAGCATCATGCGTAAGGCTTTCGGCATTTGCTcc > 1:421020/1‑75 (MQ=255) | CCCGGCGCATCAATGGGAATTTCTACTCGAATTAGCATGGTTCACCTACCGATGTCTGTTTGGTTTCAGGCGCGGTTTTCATTCCCGCCTCAACAAAATCCGCCAGTTGCAGCAGCATCATGCGTAAGGCTTTCGGCATTTGCTCC > NZ_CP009273/3294080‑3294225 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 6 ≤ ATCG/ATCG < 10 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |