Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 2,428,361 | 0 | T | G | 52.6% | 5.3 / 15.4 | 19 | Y252S (TAC→TCC) | truA | tRNA pseudouridine(38‑40) synthase TruA |
Reads supporting (aligned to +/- strand): ref base T (1/8); new base G (10/0); total (11/8) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.19e-04 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.65e-04 | |||||||||||
Rejected as consensus: E-value score below prediction cutoff. | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
TGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCC > NZ_CP009273/2428278‑2428448 | tGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGGTTTTGGAGATCATACCCGGCAGGGGAATccc > 2:1560561/1‑89 (MQ=255) aTGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACcgc < 2:633127/90‑1 (MQ=255) tCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGGAGATCATACCGGGCAGGGGAATCCACCGCGAcc > 2:1643863/1‑90 (MQ=255) tCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGGAGATCATACCGGGCAGGGGAAACCACCCCGAcc > 1:1399492/1‑90 (MQ=255) tCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGAcc < 2:1131291/90‑1 (MQ=255) cATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATAc < 1:1643863/90‑1 (MQ=255) ttattCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAAATCATACCCGGCAGGGGAATCCACCGCGACCAGATACaa > 1:1114490/1‑89 (MQ=255) gTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGGTTTTGGAGATCATACCCGGCAGGGGAATACACCCCGACCAGATACAACCCTTccc > 2:167068/1‑89 (MQ=255) tAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGGAATCCACCGCGACCAGATACAACCCTTCCCCt > 1:1228890/1‑90 (MQ=255) gTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGGAGATCATACCCGGCAGGGGAATCCACCCCGACCAGATACAGCCCCTCCCCttt > 1:1540892/1‑90 (MQ=255) gCGGGCCCATTGGCGGTTTTTGAAGAACATAACCGGCAGGGGAAATCACCCCGACCAGATACAACCCCTCCGCTTTTTCCGGTCGTgccg > 1:845415/1‑90 (MQ=255) gCGGGCCCATTGGCGGTTTTTGAAAATCATAACGGGCAGGGGAAACCACCCCCCCCCGATACAGCCCCTCCCCTTTTTCCGTTGCTgccc > 2:1637455/1‑89 (MQ=255) gCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGGCAGGGGAATCCACCGCGACCAGATACAGCCCTTCCCCTTTTTCCGTTTCTgccc > 1:1483787/1‑89 (MQ=255) gCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAg < 1:1727515/90‑1 (MQ=255) gATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCa < 2:1483787/90‑1 (MQ=255) gATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCa < 2:845415/90‑1 (MQ=255) aTCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAg > 1:1587011/1‑90 (MQ=255) cAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTa < 2:1399492/90‑1 (MQ=255) ggTAATCCCCCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATcc < 2:1631090/90‑1 (MQ=255) | TGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCC > NZ_CP009273/2428278‑2428448 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 5 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 26 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |