Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 2,738,448 | 0 | A | C | 61.5% | 4.5 / 10.9 | 13 | V84G (GTG→GGG) | trmD | tRNA (guanosine(37)‑N1)‑methyltransferase TrmD |
Reads supporting (aligned to +/- strand): ref base A (3/2); new base C (0/8); total (3/10) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 3.50e-02 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.13e-03 | |||||||||||
Rejected as consensus: E-value score below prediction cutoff. | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
GACCGCACACCAGAATCAATTTTTGATTCGTTGCCAGTTCGCTGACGCCCGCTTGATCAAGCTTGCGTCCCTGTGGTGACAGATAAATCACCTTTGCGCCTTCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAACATCCCCGGTCCGCCGCCG > NZ_CP009273/2738359‑2738537 | gACCGCACACCAGAATCAATTTTTGATTCGTTGCCAGTTCGCTGACGCCCGCTTGATCAAGCTTGCGTCCCTGTGGTGACAGATAAATCa > 1:414546/1‑90 (MQ=255) cacCAGAATCAATTTTTGATTCGTTGCCAGTTCGCTGACGCCCGCTTGATCAAGCTTGCGTCCCTGTGGTGACAGATAAATCACCTTTgc < 2:236984/90‑1 (MQ=255) ttttGATTCGTTGCCAGTTCGCTGACGCCCGCTTGATCAAGCTTGCGTCCCTGTGGTGACAGATAAATCCCCTTTGCGCCTTCACccgcc < 1:71540/90‑1 (MQ=255) cGTTGCCAGTTCGCTGACGCCCGCTTGATCAAGCTTGCGTCCCTGTGGTGACAGATAAATCACCTTTGCGCCTTCACCCGCCGCGGCttt > 2:508867/1‑90 (MQ=255) gCCAGTTCGCTGACGCCCGCTTGATCAAGCTTGCGTCCCTGTGGTGACAGATAAATCACCTTTGCGCCTTCCCCCGCCGCGGCTTTtgct < 2:414546/90‑1 (MQ=255) gcgaaCGCCCCTTTAATCAAGTTGGCCCCCGTGTGGTCACAGAAAAACCCCTTTTCCGCCTCCCCCCGCCGCGGCTTTTGCTGCATGAAt < 2:507043/86‑1 (MQ=255) cTTGATCAAGCTTGCGCCCTTGTGGTGACAGATAAACCCCCTTTGCGCCTTCCCCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCa < 2:504308/90‑1 (MQ=255) tGATCAAGCTTGCGTCCCTGTGGTGACAGATAAATCCCCTTTGCGCCTTCCCCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAg < 1:595503/90‑1 (MQ=255) cTTGCGTCCCTGTGGTGACAGATAAATCCCCTTTGCGCCTTCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCAc < 1:508197/90‑1 (MQ=255) cTTGCGCCCCTGTGGTAACAGAAAAACCCCCTTTCCCCCTCCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCAc < 2:588629/90‑1 (MQ=255) tGCGTCCCTGTGGTGACAGATAAATCCCCTTTGCGCCTTCCCCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACca < 1:374603/90‑1 (MQ=255) gACAGATAAATCCCCTTTGCGCCTTCCCCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAACATcccc < 1:137298/90‑1 (MQ=255) gACAGATAAATCACCTTTGCGCCTTCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAACATcccc > 2:531409/1‑90 (MQ=255) cAGATAAATCACCTTTGCGCCTTCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAACATCCCCgg > 1:327695/1‑90 (MQ=255) aCCTTTGCGCCTTCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAACATCCCCGGTccgccgccg > 2:350548/1‑90 (MQ=255) | GACCGCACACCAGAATCAATTTTTGATTCGTTGCCAGTTCGCTGACGCCCGCTTGATCAAGCTTGCGTCCCTGTGGTGACAGATAAATCACCTTTGCGCCTTCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAACATCCCCGGTCCGCCGCCG > NZ_CP009273/2738359‑2738537 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |