Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 2,738,448 | 0 | A | C | 54.2% | 8.5 / 23.0 | 24 | V84G (GTG→GGG) | trmD | tRNA (guanosine(37)‑N1)‑methyltransferase TrmD |
Reads supporting (aligned to +/- strand): ref base A (6/5); new base C (0/13); total (6/18) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 3.43e-03 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.02e-03 | |||||||||||
Rejected as consensus: E-value score below prediction cutoff. | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
GACCGCACACCAGAATCAATTTTTGATTCGTTGCCAGTTCGCTGACGCCCGCTTGATCAAGCTTGCGTCCCTGTGGTGACAGATAAATCACCTTTGCGCCTTCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAACATCCCCGGTCCGCC > NZ_CP009273/2738359‑2738533 | gACCGCACACCAGAATCAATTTTTGATTCGTTGCCAGTTCGCTGACGCCCGCTTGATCAAGCTTGCGTCCCTGTGGTGACAGATAAATCa < 1:380963/90‑1 (MQ=255) ccGCACACCAGAATCAATTTTTGATTCGTTGCCAGTTCGCTGACGCCCGCTTGATCAAGCTTGCGTCCCTGTGGTGACAGATAAATCAcc > 1:274050/1‑90 (MQ=255) aTTTTTGATTCGTTGCCAGTTCGCTGACGCCCGCTTGATCAAGCTTGCGTCCCTGTGGTGACAGATAAATCa < 1:1099255/72‑1 (MQ=255) aTTTTTGATTCGTTGCCAGTTCGCTGACGCCCGCTTGATCAAGCTTGCGTCCCTGTGGTGACAGATAAATCa > 2:1099255/1‑72 (MQ=255) tttGATTCGTTGCCAGTTCGCTGACGCCCGCTTGATCAAGCTTGCGTCCCTGTGGTGACAGATAAATCACCTTTGCGCCTTCACccgccg < 2:406993/90‑1 (MQ=255) tGATTCGTTGCCAGTTCGCTGACGCCCGCTTGATCAAGCTTGCGTCCCTGTGGTGACAGATAAATCACCTTTGCGCCTTCACCCGCcgcg > 2:491053/1‑90 (MQ=255) tGACGCCCGCTTGATCAAGCTTGCGTCCCTGTGGTGACAGATAAATCCCCTTTGCGCCTTCACCCGCCGCGGCTTTTGCTGCATGAATgg < 2:1060750/90‑1 (MQ=255) cTTGATCAAGCTTGCGTCCCTGTGGTGACAGAAAAATCCCCTTTGCGCCTTCCCCCGCCGCGGCTTTTGCTGCATGAATGGCATCCCGCa < 1:464508/90‑1 (MQ=255) tGATCAAGCTTGCGTCCCTGTGGTGACAGATAAATCCCCTTTGCGCCTTCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAg < 2:1046633/90‑1 (MQ=255) tGATCAAGCTTGCGCCCCTGTGGTGACAGAAAAATCCCCTTTGCGCCTTCCCCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAg < 1:720545/90‑1 (MQ=255) cTTGCGTCCCTGTGGTGACAGATAAATCCCCTTTGCGCCTTCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCAc < 1:237718/90‑1 (MQ=255) cTTGCGTCCCTGTGGTGACAGATAAATCACCTTTGCGCCTTCCCCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCAc < 2:917901/90‑1 (MQ=255) cTTGCGTCCCTGTGGTGACAGATAAATCACCTTTGCGCCTTCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCAc < 2:455235/90‑1 (MQ=255) cTTGCGTCCCTGTGGTGAAAGAAAAATCCCCTTTGCGCCTTCCCCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCAc < 1:102443/90‑1 (MQ=255) cTTGCGCCCCGGTGGTGACAGATAAATCCCCTTTGCCCCTTCCCCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCAc < 1:491053/90‑1 (MQ=255) ttGCGTCCTTGTGGTGACAGAAAAACCCCCTTTGCCCCTTCCCCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCAcc < 1:718817/90‑1 (MQ=255) ttGCGTCCCTGTGGTGACAGATAAATCCCCTTTGCGCCTTCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCAcc < 2:614854/90‑1 (MQ=255) gCGTCCCTGTGGTGACAGATAAATCACCTTTGCGCCTTCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACcat > 1:688436/1‑90 (MQ=255) tCCCTGTGGTGACAGATAAATCACCTTTGCGCCTTCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACcatcat < 2:467541/90‑1 (MQ=255) tCCCTGTGGTGACAAAAAAATCCCCTTTCCGCCTCCCCCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACcatcat < 2:763932/90‑1 (MQ=255) cgtgtgTTGAAAGATAAACCCCCTTTGCGCCTTCTCCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTa < 2:674748/88‑1 (MQ=255) cTGTGGTGACAGATAAATCACCTTTGCGCCTTCACCCGCCGCGGCTTTTGCTGCATGAATgg > 1:75269/1‑62 (MQ=255) cTGTGGTGACAGATAAATCACCTTTGCGCCTTCACCCGCCGCGGCTTTTGCTGCATGAATgg < 2:75269/62‑1 (MQ=255) ggtgACAGATAAATCACCTTTGCGCCTTCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAACATc > 2:88199/1‑90 (MQ=255) gACAGATAAATCACCTTTGCGCCTTCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAACATcccc > 2:768420/1‑90 (MQ=255) cAGATAAATCCCCTTTGCGCCTTCCCCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAACATCCCCgg < 2:274050/90‑1 (MQ=255) aaTCCCCTTTGCCCCTTCCCCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAACATCCCCGGTccgcc < 2:359377/90‑1 (MQ=255) | GACCGCACACCAGAATCAATTTTTGATTCGTTGCCAGTTCGCTGACGCCCGCTTGATCAAGCTTGCGTCCCTGTGGTGACAGATAAATCACCTTTGCGCCTTCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAACATCCCCGGTCCGCC > NZ_CP009273/2738359‑2738533 |
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: ‑ |