Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 3,203,095 | 0 | T | C | 55.6% | 5.5 / 12.8 | 18 | S270G (AGT→GGT) | tsaD | tRNA (adenosine(37)‑N6)‑threonylcarbamoyltransferase complex transferase subunit TsaD |
Reads supporting (aligned to +/- strand): ref base T (4/4); new base C (0/10); total (4/14) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 2.29e-02 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.74e-03 | |||||||||||
Rejected as consensus: E-value score below prediction cutoff. | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
TACAAAATTCCGGACGCGCGTAGAACACTTCGCCGCGGCGTTTTTTCATCATTTCAGCCAGCTTCGCCCGTAACGTACGGTTAGCA‑CTCACGCCGCCCGCCATGACCAGTCGCTTAAAGCCCGTCTGATCCAGCGCCCGCTTGCACTTAATCATCAGCGTATCGACCACCGCATCTT > NZ_CP009273/3203008‑3203184 | tACAAAATTCCGGACGCGCGTAGAACACTTCGCCGCGGCGTTTTTTCATCATTTCAGCCAGCTTCGCCCGTAACGTACGGTTAGCA‑CTCa < 1:839861/90‑1 (MQ=255) tACAAAATTCCGGACGCGCGTAGAACACTTCGCCGCGGCGTTTTTTCATCATTTCAGCCAGCTTCGCCCGTAACGTACGGTTAGCA‑CTCa > 2:448267/1‑90 (MQ=255) aCAAAATTCCGGACGCGCGTAGAACACTTCGCCGCGGCGTTTTTTCATCATTTCAGCCAGCTTCGCCCGTAACGTACGGTTAGCA‑CTCAc < 1:1471751/90‑1 (MQ=255) cgcgTAGAACACTTCGCCGCGGCGTTTTTTCATCATTTCAGCCAGCTTCGCCCGTAACGTACGGTTAGCA‑CTCACGCCGCCCGCCATGAc > 1:304955/1‑90 (MQ=255) gcgTACACCCCTCCGCCGCGCCTTTTTTTCATCTTTCCCGCCGGCTCCCCCCGTAACGTCCGGTTACCC‑CCCCCGCCGCCCGCCATGAcc < 1:143695/90‑1 (MQ=255) aGAACCCTTCCCCGGGGTTTTTTTTCATCATTTCAGCCAGCTTCGCCCGTAACGTAGGGTTAGCA‑CTCCCCCCGCCCGCCATGACCAGTc < 2:1178434/90‑1 (MQ=255) aGAACCCTCCGCCGCGGCTTTTTTTCATCTTTTCAGCCAGCTTCGCCCGTAACGTAGGGTTAGCC‑CCCCCGCCGCCCGCCATGACCAGTc < 2:1357925/90‑1 (MQ=255) tccgccgcGGCTTTTTTTCATTTTTCCACCCGGCTTCCCCCGTATCGTACGTTTACCC‑CCCACCCCGCCCGCCATGACCAGTCGCTTaaa < 1:705075/88‑1 (MQ=255) tCGCCGCGGCGTTTTTTCATCATTTCAGCCAGCTTCGCCCGTAACGTACGGTTAGCA‑CTCACGCCGCCCGCCATGACCAGTCGCTTAAAg > 1:1357925/1‑90 (MQ=255) cccgcGGCTTTTTTTCATCATTCCAGCCAGCTCCGCCCGTAACGTACGGTTACCA‑CCCCCGCCGCCCGCCATGACCAGTCGCTTAAAGcc < 2:1181481/89‑1 (MQ=255) cccgcGATGTTTTTCCCTCTTTCCGCCCAGCTCCGCCCGTAACGCCCGATTAGCC‑CCCACGCCGCCCGCCATGCCCAGTCGCTTAAAGcc < 2:864403/89‑1 (MQ=255) tCAGCCAGCTCCGCCCGTAACGTCCGGTTAGCC‑CCCCCGCCGCCCGCCATGACCAGTCGCTTAAAGCCCGTCTGATCCAGCGCCCGCTTg < 2:651675/90‑1 (MQ=255) agccggcTCCGCCCGTAACGTA‑GGTTAAAACCCCCCCCCCCCCGCCATGACCAGTCGCTTAAAGCCCGTCTGATCCAGCGCCCGCTTGCa < 1:220483/90‑1 (MQ=255) ccagcTCCCCCCGTACCGTCCGTTTACCC‑CCCCCCCCCCCCGCCATGACCAGTCGCTTAAAGCCCGTCTGATCCAGCGCCCGCTTGCACt < 1:894588/90‑1 (MQ=255) agcTCCCCCCGAAACGTCCGGTTACCC‑CCCCCCCCCCCCGCCATGACCAGTCGCTTAAAGCCCGTCTGATCCAGCGCCCGCTTGCACTTa < 1:1476916/90‑1 (MQ=255) tccGCCCGAAACGTCCGGTTAGCC‑CCCCCGCCGCCCGCCATGACCAGTCGCTTAAAGCCCGTCTGATCCAGCGCCCGCTTGCACTTAatc < 2:902010/88‑1 (MQ=255) gTACGGTTAGCA‑CTCACGCCGCCCGCCATGACCAGTc < 2:438429/37‑1 (MQ=255) gTACGGTTAGCA‑CTCACGCCGCCCGCCATGACCAGTc > 1:438429/1‑37 (MQ=255) tCACGCCGCCCGCCATGACCAGTCGCTTAAAGCCCGTCTGATCCAGCGCCCGCTTGCACTTAATCATCAGCGTATCGACCACCGCATCtt > 2:110713/1‑90 (MQ=255) | TACAAAATTCCGGACGCGCGTAGAACACTTCGCCGCGGCGTTTTTTCATCATTTCAGCCAGCTTCGCCCGTAACGTACGGTTAGCA‑CTCACGCCGCCCGCCATGACCAGTCGCTTAAAGCCCGTCTGATCCAGCGCCCGCTTGCACTTAATCATCAGCGTATCGACCACCGCATCTT > NZ_CP009273/3203008‑3203184 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 5 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 18 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |