Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 2,738,448 | 0 | A | C | 51.8% | 8.5 / 25.2 | 27 | V84G (GTG→GGG) | trmD | tRNA (guanosine(37)‑N1)‑methyltransferase TrmD |
Reads supporting (aligned to +/- strand): ref base A (11/2); new base C (0/14); total (11/16) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 5.98e-06 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.03e-05 | |||||||||||
Rejected as consensus: E-value score below prediction cutoff. | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
GACCGCACACCAGAATCAATTTTTGATTCGTTGCCAGTTCGCTGACGCCCGCTTGATCAAGCTTGCGTCCCTGTGGTGACAGATAAATCACCTTTGCGCCTTCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAACATCCCCGGTC > NZ_CP009273/2738359‑2738529 | gACCGCACACCAGAATCAATTTTTGATTCGTTGCCAGTTCGCTGACGCCCGCTTGATCAAGCTTGCGTCCCTGTGGTGACAGATAAATCa > 2:1731966/1‑90 (MQ=255) gACCGCACACCAGAATCAATTTTTGATTCGTTGCCAGTTCGCTGACGCCCGCTTGATCAAGCTTGCGTCCCTGTGGTGACAGATAAATCa > 2:1656306/1‑90 (MQ=255) ccGCACACCAGAATCAATTTTTGATTCGTTGCCAGTTCGCTGACGCCCGCTTGATCAAGCTTGCGTCCCTGTGGTGACAGATAAATCAcc > 1:160976/1‑90 (MQ=255) aaTTTTTGATTCGTTGCCAGTTCGCTGACGCCCGCTTGATCAAGCTTGCGTCCCTGTGGTGACAGATAAATCACCTTTGCGCCTTCAccc > 2:1743929/1‑90 (MQ=255) tttGATTCGTTGCCAGTTCGCTGACGCCCGCTTGATCAAGCTTGCGTCCCTGTGGTGACAGATAAATCCCCTTTGCGCCTTCACccgccg < 1:408549/90‑1 (MQ=255) ttGATTCGTTGCCAGTTCGCTGACGCCCGCTTGATCAAGCTTGCGTCCCTGTGGTGACAGATAAATCACCTTTGCGCCTTCACCcgccgc < 1:1743929/90‑1 (MQ=255) gCCAGTTCGCTGACGCCCGTTTGATCAAGCTTGCGTCCCTGTGGTGACAGATAAATCACCTTTGCGCCTTCACCCGCCGCGGCTTTtgct < 1:1180975/90‑1 (MQ=255) cccAGTTCGCTGACGCCCGCTTGATCAAGCTTGCGTCCTTGTGGTGAAAGATAAATCCCCTTTGCGCCTTCCCCCGCCGCGGCTTTtgct < 2:840489/89‑1 (MQ=255) cAGTTCGCTGACGCCCGCTTGATCAAGCTTGCGTCCCTGTGGTGACAGATAAATCCCCTTTGCGCCTTCACCCGCCGCGGCTTTtgctgc < 1:1402182/90‑1 (MQ=255) gCCCGCTTGATCAAGCTTGCGTCCCTGTGGTGACAGATAAATCa < 1:1551467/44‑1 (MQ=255) gCCCGCTTGATCAAGCTTGCGTCCCTGTGGTGACAGATAAATCa > 2:1551467/1‑44 (MQ=255) gCCCGCTTGATCAAGCTTGCGTCCCTGTGGTGACAGATAAATCACCTTTGCGCCTTCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTc > 1:550933/1‑90 (MQ=255) cTTGATCACGCTGGCGTCCCTGTGGTGACAAAAAAATCCCCTTTGCCCCTCCCCCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCa < 2:1221085/90‑1 (MQ=255) cTTGACCAAGCTTGCGTCCCTGTGGTGACAGATAAATCCCCTTTGCGCCTTCCCCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCa < 2:902462/90‑1 (MQ=255) cAAGCTTGCGTCCCTGTGGTGATAGATAAATCCCCTTTGCGCCTTCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGtt < 1:1662388/90‑1 (MQ=255) cAAGCTTGCGTCCCTGTGGTGACAGATAAATCACCTTTGCGCCTTCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGtt > 2:394074/1‑90 (MQ=255) aaGCTTGCGTCCCTGTGGTGACAGATAAATCCCCTTTGCGCCTTCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTg < 1:1173552/90‑1 (MQ=255) aGCTTTCGTCCCTGTGGTGACAGATAAATCACCTTTGCGCCTTCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGc > 2:1735238/1‑90 (MQ=255) ctggCGTCCCTGTGGTGACAAAAAAATCCCCTTTGCGCCTTCCCCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCAc < 2:160976/87‑1 (MQ=255) cTTGCGTCCCTGTGGTGACAGATAAATCCCCTTTGCGCTTTCCCCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCAc < 2:578781/90‑1 (MQ=255) aggTCCTTGTGGTGACAGATAAACCCCCTTTGCGCCTTCCCCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACcat < 2:1595890/88‑1 (MQ=255) tCCCTGTGGTGACAGATAAATCCCCTTTGCCCCTTCCCCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACcatcat < 1:116316/90‑1 (MQ=255) ggtgACAGATAAATCCCCTTTGCGCCTTCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACcatcat < 2:906263/83‑1 (MQ=255) ggtgACAGATAAATCACCTTTGCGCCTTCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACcatcat > 1:906263/1‑83 (MQ=255) tgACAGATAAATCACCTTTGCGCCTTCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAACATcct > 1:1595890/1‑89 (MQ=255) gACAGATAAATCCCCTTTGCGCCTTCACCCGCCGCGGCTTTTGCTGCATg < 2:1580094/50‑1 (MQ=255) gACAGATAAATCACCTTTGCGCCTTCACCCGCCGCGGCTTTTGCTGCATg > 1:1580094/1‑50 (MQ=255) gACAGATAAATCACCTTTGCGCCTTCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAACATcccc > 1:110375/1‑90 (MQ=255) cAGATAAATCCCCTTTGCGCCTTCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCACGCAAGGGTTGCACCATCATTAACATCCCCgg < 1:1653594/90‑1 (MQ=255) cAGATAAATCACCTTTGCGCCTTCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAACATCCCCgg > 2:574744/1‑90 (MQ=255) gATAAATCACCTTTGCGCCTTCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAACATCCCCGGTc > 1:250210/1‑90 (MQ=255) | GACCGCACACCAGAATCAATTTTTGATTCGTTGCCAGTTCGCTGACGCCCGCTTGATCAAGCTTGCGTCCCTGTGGTGACAGATAAATCACCTTTGCGCCTTCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAACATCCCCGGTC > NZ_CP009273/2738359‑2738529 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 17 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |