Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 2,738,448 | 0 | A | C | 63.6% | 7.8 / 30.5 | 34 | V84G (GTG→GGG) | trmD | tRNA (guanosine(37)‑N1)‑methyltransferase TrmD |
Reads supporting (aligned to +/- strand): ref base A (9/3); new base C (0/21); total (9/25) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 5.70e-06 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.79e-06 | |||||||||||
Rejected as consensus: E-value score below prediction cutoff. | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
GACCGCACACCAGAATCAATTTTTGATTCGTTGCCAGTTCGCTGACGCCCGCTTGATCAAGCTTGCGTCCCTGTGGTGACAGATAAATCACCTTTGCGCCTTCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAACATCCCCGGTCCGCCGCCG > NZ_CP009273/2738359‑2738537 | gACCGCACACCAGAATCAATTTTTGATTCGTTGCCAGTTCGCTGACGCCCGCTTGATCAAGCTTGCGTCCCTGTGGTGACAGATAAATCa > 2:671316/1‑90 (MQ=255) gACCGCACACCAGAATCAATTTTTGATTCGTTGCCAGTTCGCTGACGCCCGCTTGATCAAGCTTGCGTCCCTGTGGTGACAGATAAATCa > 2:657064/1‑90 (MQ=255) acacCAGAATCAATTTTTGATTCGTTGCCAGTTCGCTGACGCCCGCTTGATCAAGCTTGCGTCCCTGTGGTGACAGATAAATCACCTTTg < 1:460622/90‑1 (MQ=255) aaTTTTTGATTCGTTGCCAGTTCGCTGACGCCCGCTTGATCAAGCTTGCGTCCCTGTGGTGACAGATAAATCACCTTTGCGCCTTCAccc > 1:20419/1‑90 (MQ=255) ttGATTCGTTGCCAGTTCGCTGACCCCCGCTTGATCACGCTTGCGTCCTTGTGGTGAAAGATAAATCCCCTTTGCCCCTCCCCCcgccgc < 1:572751/90‑1 (MQ=255) gCCAGTTCGCTGACGCCCGCTTGATCAAGCTTGCGTCCCTGTGGTGACAGATAAATCCCCTTTGCGCCTTCCCCCGCCGCGGCTTTtgct < 1:603851/90‑1 (MQ=255) gCCAGTTCGCTGACGCCCGCTTGATCAAGCTTGCGTCCCTGTGGTGACAGATAAATCACCTTTGCGCCTTCACCCGCCGCGGCTTTtgct < 2:5483/90‑1 (MQ=255) atttCGTTGACGCCGGCTTGACCAAGCTTGCGGCCCTGTGGTGACAGAAAAATCTCCTTTGCACCTTCCCCCGCCGCGGCTTTTGCTGCa < 2:990815/88‑1 (MQ=255) tCGCTGACGCCCGCTTGATCAAGCTTGCGTCCCTGTGGTGACAGATAAATCACCTTTGCGCCTTCACCCGCCGCGGCTTTTGCTGCATGa > 1:820534/1‑90 (MQ=255) cGCCCGCTTGATCAAGCTTGCGTCCCTGTGGTGACAGATAAATCCCCTTTGCGCCTTCACCCGCCGCGGCTTTTGCTGCATGAATGGCGt < 2:508060/90‑1 (MQ=255) cGCCCCTTTGACAAACCTTCCGCCCCTGTGGTGGCAGATAAACCCCTTTGGCGCCTTCCCCCGCCGCGGCTTTTGCTGCATGAATGGCGt < 2:980039/90‑1 (MQ=255) gCCCGCTTGATCAAGCTTGCGTCCCTGTGGTGACAGATAAATCACCTTTGCGCCTTCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTc > 1:35518/1‑90 (MQ=255) gCCCGCTTGACCAAGCTTGCGTCCCTGTGGTGACAGATAAACCCCCTTTGCGCTTTCCCCCGCCGCGGCTTTTGCTGCATGAATGGCGTc < 2:970051/90‑1 (MQ=255) tttGATCAAGCTTGCGTCCCTGTGGTGACAAAAAAACCCCCTTTGGGCCTTCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCa < 1:1077260/89‑1 (MQ=255) cTTGATCAAGCTTGCGTCCCTGTGGTGACAGATAAATCCCCTTTGCGCCTTCCCCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCa < 1:124085/90‑1 (MQ=255) cTTGATCAAGCTTGCGTCCCTGTGGTGACAGATAAATCCCCTTTGCGCCTTCCCCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCa < 2:890031/90‑1 (MQ=255) cTTGATCAAGCTTGCGTCCCTGTGGTGACAGATAAATCCCCTTTGCGCCTTCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCa < 2:801116/90‑1 (MQ=255) cTTGAAAAAGTTTCCGTCCCTGTGGTGACAGATAAACCCCCTTTGCGCCTTCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCa < 1:379151/90‑1 (MQ=255) ttaaTCAAGCTTGTGTCCCTGTGGTGACAGATAAATCCCCTTTGCGCCTTCCCCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCaa < 2:525144/87‑1 (MQ=255) cTTGCGTCCTTGTGGTGACAGATAAATCCCCTTTGCGCCTTCCCCCGCCGCGGGTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCAc < 1:368139/90‑1 (MQ=255) cTTGCGTCCTTGTGGTGACAGATAAATCCCCTTTGCGCCTTCCCCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCAc < 1:249803/90‑1 (MQ=255) cTTGCGTCCTTGTGGTGACAGAAAAATCCCCTTTGCCCCTTCCCCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCAc < 1:671316/90‑1 (MQ=255) cTTGCGTCCCTGTGGTGACAGATAAATCCCCTTTGCGCCTTCCCCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCAc < 1:50345/90‑1 (MQ=255) cTTGCGTCCCTGTGGTGACAGATAAATCCCCTTTGCGCCTTCCCCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCAc < 2:490181/90‑1 (MQ=255) cTTGCGCCCCTGTGGTGACAGATAAATCCCCTTTGCGCCTTCCCCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCAc < 1:1033734/90‑1 (MQ=255) cTTCCGCCCCTGGGGTGACAGATAAATCCCCTTCGCGCCTTCCCCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCAc < 1:580976/90‑1 (MQ=255) ccTGTGGTGACAGATAAATCACCTTTGCGCCTTCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTa > 2:665549/1‑90 (MQ=255) ccTGTGGTGACAGATAAATCACCTTTGCGCCTTCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTa > 2:678475/1‑90 (MQ=255) ccTGTGGTGACAGATAAATCACCTTTGCGCCTTCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTa > 1:228986/1‑90 (MQ=255) ccTGTGGTGACAGATAAATCACCTTTGCGCCTTCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTa > 1:443812/1‑90 (MQ=255) tgtgGTGACAGATAAATCCCCTTTGCGCCTTCCCCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAAc < 2:1012660/90‑1 (MQ=255) ggtgACAGATAAATCACCTTTGCGCCTTCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAACATc > 1:312890/1‑90 (MQ=255) gACAGATAAACCCCCTTCGCCCCTCCCCCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAACATcccc < 1:683035/90‑1 (MQ=255) cAGATAAATCACCTTTGCGCCTTCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAACATCCCCgg > 1:825319/1‑90 (MQ=255) aaTCCCTTTTGCGCCTCCCCCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAACATCCCCGGTccgcc < 1:403564/90‑1 (MQ=255) tCACCTTTGCGCCTTCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAACATCCCCGGTccgccgc < 1:678475/90‑1 (MQ=255) aCCTTTGCGCCTTCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAACATCCCCGGTccgccgccg > 1:619610/1‑90 (MQ=255) | GACCGCACACCAGAATCAATTTTTGATTCGTTGCCAGTTCGCTGACGCCCGCTTGATCAAGCTTGCGTCCCTGTGGTGACAGATAAATCACCTTTGCGCCTTCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAACATCCCCGGTCCGCCGCCG > NZ_CP009273/2738359‑2738537 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 5 ≤ ATCG/ATCG < 10 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |