| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | ACB122 | 5,009,861 | 0 | A | . | 26.3% | 44.5 / 14.9 | 19 | intergenic (‑156/+78) | CJ019_04705/spuC_5 | hypothetical protein/Putrescine‑‑pyruvate aminotransferase |
| Reads supporting (aligned to +/- strand): ref base A (10/4); new base . (2/3); total (12/7) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 3.05e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.60e-01 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
ATGAGATGAATAAGCCCGCACCGGTCTCCGGCGCGGGCTTTTACTTTCCGGGCCTTGGACAGTATCTGACCCGAAAGTACGTCCTTGTGAACTGATGACAGCGTCTTGGCTTTGATTGTGATAGCACTAC‑AAAAAAAA‑‑AAAGCCCGCCTAAAGGCGGGCTTAAAACGAAGAGGGTAATTCAGCTGGCATCATTGCCAACGTTGATGGCTCAGACCTTCAGCTCATCAAGGACGCTCAGCACAGCTTCGTGTGC > ACB122/5009731‑5009983 | atgagatgaATAAGCCCGCACCGGTCTCCGGCGCGGGCTTTTACTTTCCGGGCCTTGGACAGTATCTGACCCGAAAGTACGTCCTTGTGAACTGATGACAGCGTCTTGGCTTTGATTGTGATAGCACTAC‑aaaa > 4:50566/1‑134 (MQ=255) attaaTACGCCCGCACCGGTCTCCGGCGCGGGCTTTTACTTTCCGGGCCTTGGACATTATCTGACCCGAAAGTACGTCCTTGTGAACTTATGACAGCGTCTTGGCTTTGATTGTGATAGCACTAC‑aaaaaaaa‑‑aa < 6:172566/132‑1 (MQ=255) cGCACCGGTCTCCGGCGCGGGCTTTTACTTTCCGGGCCTTGGACAGTATCTGACCCGAAAGTACGTCCTTGTGAACTGATGACAGCGTCTTGGCTTTGATTGTGTTTGCACTAC‑AAAAAAAA‑‑AAAGCCCGCCt > 5:114996/1‑133 (MQ=255) cGGTCTCCGGCGCGGGCTTTTACTTTCCGGGCCTTGGACAGTATCTGACCCGAAAGTACGTCCTTGTGAACTGATGACAGCGTCTTGGCTTTGATTGTGCTAGCACTAC‑AAAAAAAA‑‑AAAGCCCGCCTAACGGCg > 3:60546/1‑135 (MQ=255) tGACCCGAAAGTACGTCCTTGTGAACTGATGACAGCGTCTTGGATTTGATTGTGATAGCACTAC‑AAAAAAAA‑‑AAAGCCCGCCTAAAGGCGGGCTTAAAACGAAGAGGGTAATTCAGCTGGCATCATTGCCAACg < 8:139674/134‑1 (MQ=255) gACCCGAAAGTACGTCCTTGTGAACTGATGACAGCGTCTTGGCTTTGATTGTGATAGCACTAC‑AAAAAAAA‑‑AAAGCCCGCCTAAAGGCGGGCTTAAAACGAAGAGGGTAATTCAGCTGGCATCATTGCCAACGtt > 6:204777/1‑135 (MQ=255) taatgTAGGTCCTTGTGAAGTTCTGACTGCGGCTTGGCTTTGTTTGTGTTTCCCCTAC‑AAAAAAAA‑‑AAAGCCCGCCTAAAGGCGGGCTTAAAACGAAGAGGGTAATTCAGCTGGCATCATTGCCAACGTTGATgg < 3:50566/131‑1 (MQ=255) caacgaaCGTCTTTGTGACCTGATGCCAGCGTCTGGGCTTTGATTGTAATAGAACAAA‑AAAAAAAA‑‑AAAGCCCGCCTAAAGGCGGGCTTAAAACGAAGAGGGTAATTCAGCTGGCATCATTGCCAACGTTGATgg < 3:129847/129‑1 (MQ=255) aGTACGTCCTTGTGAACTGATGACAGCGTCTTGGCTTTGATTGTGATAGCACTAC‑AAAAAAAA‑‑AAAGCCCGCCTAAAGGCGGGCTTAAAACGAAGAGGGTAATTCAGCTGGCATCATTGCCAACGTTGATGGc > 2:132220/1‑133 (MQ=255) aGTACGTCCTTGTGAACTGAAGACAGCGTCTTGGCTTTGATTGTGATAGCACTAC‑‑AAAAAAA‑‑AAATCCAGCCTAAAGGCGGGCTTAAAACGAACAGGGTATTTCAGCTGGCATCATTGCCATCGTTGATGGc > 7:188424/1‑132 (MQ=255) tCCTTGTGAACTGATGACAGCGTCTTGGCTTTGATTGTGATAGCACTAC‑‑AAAAAAA‑‑AAAGCCCGCCTAAAGGCGGGCTTAAAACGAAGAGGGTAATTCAGCTGGCATCATTGCCAACGTTGATGGCTCAGACCtt < 5:35577/135‑1 (MQ=255) gAACTGATGACAGCGTCTTGGCTTTGATTGTGATAGCACTAC‑AAAAAAAA‑‑AAAGCCCGCCTAAAGGCGGGCTTAAAACGAAGAGGGTAATTCAGCTGGCATCATTGCCAACGTTGATGGCTCAGACCTTCAGCtc > 5:169835/1‑135 (MQ=255) aaCTGATGACAGCGTCTTGGCTTTGATTGTGATAGCACTAC‑AAAAAAAA‑‑AAAGCCCGCCTAAAGGCGGGCTTAAAACGAAGAGGGTAATTCAGCTGGCATCATTGCCAACGTTGATGGCTCAGCCCTTCAGCtca > 5:105241/1‑135 (MQ=255) gaCAGCGTCTTGGCTTTGATTGTGATAGCACTAC‑‑AAAAAAA‑‑AAAGCCCGCCTAAAGGCGGGCTTAAAACGAAGAGGGTAATTCAGCTGGCATCATTGCCAACGTTGATGGCTCAGACCTTCAGCTCATCAAGGAc < 3:100254/135‑1 (MQ=255) aCAGCGTCTTGGCTTTGATTGTGATAGCACTAC‑AAAAAAAA‑‑AAAGCCCGCCTAAAGGCGGGCTTAAAAAGAAGAGGGTAATTCAGCTGGCATCATTGCCAACGTTGATGGCTCAGACCTTCAGCACATCAAGGAc > 3:198640/1‑135 (MQ=255) aGCGGCTTGGCTTTTTTTGTTATACAACTAC‑AAAAAAAAATAAAGCCCGCCTAAAGGCGGGCTTAAAACGAAGAGGGTAATTCAGCTGGCATCATTGCCAACGTTGATGGCTCAGACCTTCAGCTCATCAAGGAc < 1:138355/135‑1 (MQ=255) gCGTCTTGGCTTTGATTGTGATAGCACTACAAAAAAAAA‑‑AAAGCCCCCCAAAAGGGGGGGTTAAAACGACGAGGGTATTTAAGCTGGCAGCATTGCCAACGTTGATGGCTCAGACATGCAGCGCCTCAGGCACGc > 6:125044/1‑135 (MQ=255) gTCTTGGCTTTGATTGTGATAGCACTAC‑AAAAAAAA‑‑AAAGCCCGCCTAAAGGCGGGCTTAAAACGAAGAGGGTAATTCAGCTGGCATCATTGCCAACGTTGATGGCTCAGACCTTCAGCTCATCAAGGACGCTCa > 8:20017/1‑135 (MQ=255) ggatttGATTGTGATAGAACTAC‑AAAAAAAA‑‑AAATCCCCCCTAAATGCTGCCTTAAAAAGAATAGGGAAATTAAGATGTCATCATTGAAAAAGATGATGGCTAAGAACTTCATCTAATCAAGGACGCTCAGcaca > 6:200350/4‑135 (MQ=255) atttttatatCACTAC‑‑AAAAAAA‑‑AAAGCCCGCCTAAAGGCGGGATTAAAACGAAGAGGGTAATTCAGCTGGCATCATTGCCAACGTTGATGGCTCAGACCTTCAGCTCATCAAGGACGCTCAGCACATCTTCgt < 8:135328/128‑1 (MQ=255) gtgATAGCACTAC‑‑AAAAAAA‑‑AAAGCCCGCCTAAAGGCGGGCTTAAAACGAAGAGGGTAATTCAGCTGGCATAATTGCCAACGTTGATGGGTCAGACCTTCAGCTCATCAAGGACGCTCAGCAAAGCTACGtgtgg > 8:30768/1‑134 (MQ=255) | ATGAGATGAATAAGCCCGCACCGGTCTCCGGCGCGGGCTTTTACTTTCCGGGCCTTGGACAGTATCTGACCCGAAAGTACGTCCTTGTGAACTGATGACAGCGTCTTGGCTTTGATTGTGATAGCACTAC‑AAAAAAAA‑‑AAAGCCCGCCTAAAGGCGGGCTTAAAACGAAGAGGGTAATTCAGCTGGCATCATTGCCAACGTTGATGGCTCAGACCTTCAGCTCATCAAGGACGCTCAGCACAGCTTCGTGTGC > ACB122/5009731‑5009983 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |