Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 3,818,086 | 0 | T | G | 54.5% | 6.0 / 11.2 | 11 | V71G (GTA→GGA) | trmH | tRNA (guanosine(18)‑2'‑O)‑methyltransferase TrmH |
Reads supporting (aligned to +/- strand): ref base T (2/3); new base G (6/0); total (8/3) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 6.06e-02 | |||||||||||
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. |
CGTTGGCGTACATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGCAGCGGCGGGTAGTAACAGCTGGGTACAGGTGAAAACACACCGCACCATTGGCGATGCCGTCGCTCATCTCAAAGGCCAGGGCATGCAGATTCTGGCAACCCATCTTTCTG > NZ_CP009273/3818000‑3818172 | cGTTGGCGTACATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGCAGCGGCGGGGAGTAACAGCTGGGGACa > 1:90411/1‑90 (MQ=255) aCATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGCAGCGGCGGGTAGTAACAGCTGGGGACAGGGGaaaaa > 1:147122/1‑89 (MQ=255) gTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGCAGCGGCGGGGAGTAAACGCTGGGGACAAGGGGAAAAACAcccc > 2:155438/1‑88 (MQ=255) ggCCTGGTAGCCGCATGCGCACCATGGCTTCGGCAGCGGCGGGTAGTAACAGCTGGGGACAGGGGAAAACACACCGCACCCTTGGCGATg > 1:333178/1‑90 (MQ=255) gTAGCCGCATGCGCACCATGGCTTCGGCAGCGGCGGGGAGGAACAGCTGGGGACAGGGGAAAAAACACCCCACCATTGGCGATGCCGGCg > 1:243679/1‑90 (MQ=255) gcgcACCATGGCTTCGGCAGCGGCGGGTAGTAACAGCTGGGTACAGGTGAAAACACACCGCACCATTGGCGATGCCGTCGCTCATCTCaa < 2:243679/90‑1 (MQ=255) gcgcACCATGGCTTCGGCAGCGGCGGGTAGTAACAGCTGGGTACAGGTGAAAACACACCGCACCATTGGCGATGCCGTCGCTCATCTCaa < 2:30560/90‑1 (MQ=255) tCGGCAGCGGCGGGGAGTAAAAGCTGGGGACAGGTGAAAACACACCGCACCATTGGCGATGCCGGCGCTCATTTCCAAAGCCAGGGCATg > 1:387586/1‑90 (MQ=255) gCTGGGTACAGGTGAAAACACACCGCACCATTGGTGATGCCGTCGCTCATCTCAAAGGCCAGGGCATGCAGATTCTGGCAACCCATCttt < 1:258645/90‑1 (MQ=255) gCTGGGTACAGGTGAAAACACACCGCACCATTGGCGATGCCGTCGCTCATCTCAAAGGCCAGGGCATGCAGATTCTGGCAACCCATCttt > 2:318493/1‑90 (MQ=255) gggTACAGGTGAAAACACACCGCACCATTGGCGATGCCGTCGCTCATCTCAAAGGCCAGGGCATGCAGATTCTGGCAACCCATCTTTCTg > 2:214127/1‑90 (MQ=255) | CGTTGGCGTACATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGCAGCGGCGGGTAGTAACAGCTGGGTACAGGTGAAAACACACCGCACCATTGGCGATGCCGTCGCTCATCTCAAAGGCCAGGGCATGCAGATTCTGGCAACCCATCTTTCTG > NZ_CP009273/3818000‑3818172 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |