| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NZ_CP009273 | 2,738,448 | 0 | A | C | 70.8% | ‑5.8 / 15.5 | 24 | V84G (GTG→GGG) | trmD | tRNA (guanosine(37)‑N1)‑methyltransferase TrmD |
| Reads supporting (aligned to +/- strand): ref base A (6/1); new base C (0/17); total (6/18) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 5.20e-05 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.93e-05 | |||||||||||
| Rejected as consensus: E-value score below prediction cutoff. | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
GACCGCACACCAGAATCAATTTTTGATTCGTTGCCAGTTCGCTGACGCCCGCTTGATCAAGCTTGCGTCCCTGTGGTGACAGATAAATCACCTTTGCGCCTTCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAACATCCCCGGTCCGCCGCCG > NZ_CP009273/2738359‑2738537 | gACCGCACACCAGAATCAATTTTTGATTCGTTGCCAGTTCGCTGACGCCCGCTTGATCAAGCTTGCGTCCCTGTGGTGACAGATAAATCa > 2:23624/1‑90 (MQ=255)gACCGCACACCAGAATCAATTTTTGATTCGTTGCCAGTTCGCTGACGCCCGCTTGATCAAGCTTGCGTCCCTGTGGTGACAGATAAATCa > 2:660976/1‑90 (MQ=255) ccGCACACCAGAATCAATTTTTGATTCGTTGCCAGTTCGCTGACGCCCGCTTGATCAAGCTTGCGTCCCTGTGGTGACAGATAAATCAcc > 2:670750/1‑90 (MQ=255) gCACACCAGAATCAATTTTTGATTCGTTGCCAGTTCGCTGACGCCCGCTTGATCAAGCTTGCGTCCCTGTGGTGACAGATAAATCACCtt > 2:386140/1‑90 (MQ=255) acacCAGAATCAATTTTTGATTCGTTGCCAGTTCGCTGACGCCCGCTTGATCAAGCTTGCGTCCCTGTGGTGACAGATAAATCACCTTTg < 1:663561/90‑1 (MQ=255) gTTGCCGGTTCGGGGACGCCCCCTTGACCAAGCTTGCGTCCCTGTGGTGACAGATAAATCCCCTTTGCTCCTTCGCCCGCCGCGGCtttt < 1:650732/90‑1 (MQ=255) gTTGCCAGTTCGCTGACGCCCGCTTGATCAAGCTTGCGTCCCTGTGGTGACAGATAAATCACCTTTGCGCCTTCACCCGCCGCGGCtttt > 2:84856/1‑90 (MQ=255) gCCAGTTCGCTGACGCCCGCTTGATCAAGCTTGCGTCCCTGTGGTGACAGATAAATCCCCTTTGCGCCTTCCCCCGCCGCGGCTTTtgct < 1:238337/90‑1 (MQ=255) gCCAGTTCGCTGACGCCCGCTTGATCAAGCTTGCGTCCCTGTGGTGACAGAAAAATCCCCTTTGCGCCTTCCCCCGCCGCGGCTTTtgct < 2:537887/90‑1 (MQ=255) tCGCTGACGCCCGCTTGATCAAGCTTGCGTCCCTGTGGTGACAGATAAATCACCTTTGCGCCTTCACCCGCCGCGGCTTTTGCTGCATGa > 1:275193/1‑90 (MQ=255) gCTGCCCCCCGCTTGATCAAGCTTGCGCCCCTGGGGTGACAGAAAAACCCCCTTTGCCCCTTCCCCCGCCGCGGCTTTTGCTGCATGAAt < 1:197175/90‑1 (MQ=255) cccccGCTTGCCCAATCTTGCGTCCCTGTGGGGAAACATAAATCCCCTTCGGGCTTTCTCCCGCCGCGGCTTTTGCTGCATGAATGGCGt < 2:58837/88‑1 (MQ=255) gCCCGCTTGATCAAGCTTGCGTCCCTGTGGTGACAGATAAATCCCCTTTGCGCCTTCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTc < 1:666175/90‑1 (MQ=255) gttGCGGCCCGGTGGAAAAATAAAAATCCCCTTGCCGCCCTCCCCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCAc < 2:41118/89‑1 (MQ=255) ctggcgcCCTTGTGGGGACAAAAAAACCCCCTTTGCGCCTTCCCCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCAc < 2:244948/87‑1 (MQ=255) cTTGCGTCCCTGTGGTGACAGAAAAATCCCCTTTGCCCCTTCCCCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCAc < 2:536249/90‑1 (MQ=255) cTTGCGCCCCTGTGGTGACAGAAAAACCCCCTTTGCGCTTTCCCCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCAc < 2:196224/90‑1 (MQ=255) tGCGTCCCTGTGGTGACAGATAAATCCCCTTTGCGCCTTCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACca < 1:30374/90‑1 (MQ=255) tCCCTGTGGTGACAGATAAATCCCCTTTGCGCCTTCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACcatcat < 1:614672/90‑1 (MQ=255) ccccTGTGGTGACAGATAAATCCCCTTTGCCCCTTCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACcatcat < 1:386140/89‑1 (MQ=255) cttgtgGTGACAGATAAATCCCCTTTGCCCCTTCCCCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTa < 2:275193/88‑1 (MQ=255) tgtgGTGACAGATAAATCCCCTTTGCCCCTTCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAAc < 1:23624/90‑1 (MQ=255) ggtgACAGATAAATCACCTTTGCGCCTTCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAACATc > 2:82935/1‑90 (MQ=255) gACAGATAAATCACCTTTGCGCCTTCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAACATcccc > 1:456003/1‑90 (MQ=255) gATAAATCCCCTTTGCGCCTTCCCCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAACATCCCCGGTc < 1:660976/90‑1 (MQ=255) aaTCCCCTTTGCGCCTCCCCCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAACATCCCCGGTccgcc < 1:84856/90‑1 (MQ=255) aCCTTTGCGCCTTCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAACATCCCCGGTccgccgccg > 1:324461/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 < 12 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |