| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 225,465 | 0 | A | C | 25.0% | 35.4 / 15.5 | 36 | intergenic (+8/‑35) | ileV/alaV | tRNA‑Ile/tRNA‑Ala |
| Reads supporting (aligned to +/- strand): ref base A (16/11); new base C (9/0); total (25/11) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 3.41e-02 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.04e-01 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
AAGTCCACTCAGGCCTACCAAATTTGCACGGCAAATTTGAAGAGGTTTTAACTACATGTTATGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCTTTGCACGCAGGAGGTCTGCGGTTCGATCCCGCATAGCTCCACCATCTCTGTAGTGGTTAAATAAAAAATACTTCAGAGTGTAC > NC_000913/225438‑225614 | aaGTCCACTCAGGCCTACCAAATTTGCACGGCAAATTTGAAGAGGTTTTAACTACATGTTATGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCTTTGCACGCAGGAGGTCTGCGGTTCGATCCCGCATAGCTCCACCATCTCTGTAGTgg > 1:969095/1‑150 (MQ=11)aaGTCCACTCAGGCCTACCAAATTTGCACGGCAAATTTGAAGAGGTTTTAACTACATGTTATGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCTTTGCACGCAGGAGGTCTGCGGTTCGATCCCGCATAGCTCCACCATCTCTGTAGTgg < 2:1550129/150‑1 (MQ=11) aGTCCACTCAGGCCTACCAAATTTGCACGGCACATTTGAAGAGGTTTTAACTACATGGTATGGGGCTAAAGCTCAGCTGGGAGCGCGCCTGCTTTGCACGCAGGCGGTCTGCGGTTCGATCCCGCATAGCTCCACCATCGCTGTAGTGGt > 1:354824/1‑150 (MQ=11) aGTCCACTCAGGCCTACCAAATTTGCACGGCAAATTTGAAGAGGTTTTAACTACATGTTATGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCTTTGCACGCAGGAGGTCTGCGGTTCGATCCCGCATAGCTCCACCATCTCTGTAGTGGt < 2:1892314/150‑1 (MQ=31) aGTCCACTCAGGCCTACCAAATTTGCACGGCAAATTTGAAGAGGTTTTAACTACATGTTATGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCTTTGCACGCAGGAGGTCTGCGGTTCGATCCCGCATAGCTCCACCATCTCTGTAGTGGt > 2:1708777/1‑150 (MQ=31) aGTCCACTCAGGCCTACCAAATTTGCACGGCAAATTTGAAGAGGGTTTAACTACATGTTATGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCTTTGCACGCAGGAGGTCTGCGGTTCGATCCCGCATAGCTCCACCATCTCTGTAGTGGt > 2:1617541/1‑150 (MQ=14) gTCCACTCAGGCCTACCAAATTTGCACGGCAAATTTGAAGAGGTTTTAACTACATGTTATGGTGCTATAGCTCAGCTGGGAGAGCGCCTGCTTTGCACGCAGGAGGTCTGCGGTTCGATCCCGCATAGCTCCACCATCTCTGTAGTGGtt < 1:1343879/150‑1 (MQ=17) tCCACTCAGGCCTAACAAATTTGCACGGCAAATTTGAAGAGGTTTTAACTACATGTTATGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCTTTGCACGCAGGAGGTCTGCGGTTCGATCCCGCATAGCTCCACCATCTCTGTAGTGGTta < 2:1889073/150‑1 (MQ=17) ccACTCAGGCCTACCAAATTTGCACGGCAAATTTGAAGAGGTTTTAACTACATGTTATGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCTTTGCACGCAGGAGGTCTGCGGTTCGATCCCGCATAGCTCCACCATCTCTGTAGTGGTtaa > 1:1931185/1‑150 (MQ=32) aCTCAGGCCTACCAAATTTGCACGGCAAATTTGAAGAGGTTTTAACTACATGTTATGGGTCTATAGCTCAGCTGGGAGAGCGCCTGCTTTGCACGCAGGAGGTCTGCGGTTCGATCCCGCATAGCTCCACCATCTCTGTAGTGGTtaaat < 2:1391284/150‑1 (MQ=17) aCTCAGGCCTACCAAATTTGCACGGCAAATTTGAAGAGGTTTTAACTACATGTTATGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCTTTGCACGCAGGAGGTCTGCGGTTCGATCCCGCATAGCTCCACCATCTCTGTAGTGGTtaaat < 2:955395/150‑1 (MQ=32) aCTCAGGCCTACCAAATTTGCACGGCAAATTTGAAGAGGTTTTAACTACATGTTATGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCTTTGCACGCAGGAGGTCTGCGGTTCGATCCCGCATAGCTCCACCATCTCTGTAGTGGTtaaat < 1:1560111/150‑1 (MQ=32) ggCCTTCGAAATTTGCAAGGCAAATTTGAAGATGTTTTAACTACCTGCTTTGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCTTTGCACGCAGGAGGTCTGCGGTTCGATCCCGCATAGCTCCACCATCTCTGAAGTGGTTAAATaaaaa < 2:402946/150‑1 (MQ=11) ggCCTACCAAATTTGCACGGCAAATTTGAAGAGGTTTTAACTACATGTTATGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCTTTGCACGCAGGAGGTCTGCGGTTCGATCCCGCATAGCTCCACCATCTCTGTAGTGGTTAAATaaaaa > 1:435404/1‑150 (MQ=32) gCCTACCAAATTTGCACGGCAAATTTGAAGAGGTTTTAACTACATGTTATGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCTTTGCACGCAGGAGGTCTGCGGTTCGATCCCGCATAGCTCCACCATCTCTGTAGTGGTTAATTaaaaaa > 1:1905343/1‑150 (MQ=17) gCCTACCAAATTTGCACGGCAAATTTGAAGAGGTTTTAACTACATGTTATGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCTTTGCACGCAGGAGGTCTGCGGTTCGATCCCGCATAGCTCCACCATCTCTGTAGTGGTTAAATaaaaaa > 1:2082082/1‑150 (MQ=32) cctcttcAAATTTGCCGTGCAAATTTGAAGAGGTTTTAACTACATGTTATGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCTTTGCACGCAGCAGGTCTGCCGATCGATCCCGCATAGCTCCACCATCTCTGTAGTGGTTAAATaaaaaa > 2:1053427/7‑150 (MQ=11) ccTACCAAATTTGCACGGCAAATTTGAAGAGGTTTTAACTACATGTTATGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCTTTGCACGCAGGAGGTCTGCGGTTCGATCCCGCATAGCTCCACCATCTCTGTAGTGGTTAAATAAAAAAt > 1:563890/1‑150 (MQ=32) cttcAAATTTGCCGTGCAAATTTGAAGAGGTTTTAACTACATGTTATGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCTTTGCACGCAGGAGGTCTGCGGTTCGATCCCGCATAGCTCCACCATCTCTGTAGTGGTTAAATAAAAAATAc > 2:626525/4‑150 (MQ=255) cttcAAATTTGCCGTGCAAATTTGAAGAGGTTTTAACTACATGTTATGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCTTTGCACGCAGGAGGTCTGCGGTTCGATCCCGCATAGCTCCACCATCTCTGTAGTGGTTAAATAAAAAATAc > 2:68590/4‑150 (MQ=255) ttcAAATTTGCCGTGCAAATTTGAAGAGGTTTTAACTACATGTTATGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCTTTGCACGCAGGAGGTCTGCGGTTCGATCCCGCATAGCTCCACCATCTCTGTAGTGGTTAAATAAAAAATACt > 2:1016178/3‑150 (MQ=12) ttcAAATTTGCCGTGCAAATTTGAAGAGGTTTTAACTACATGTTATGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCTTTGCACGCAGGAGGTCTGCGGTTCGATCCCGCATAGCTCCACCATCTCTGTAGTGGTTAAATAAAAAATACt > 1:807775/3‑150 (MQ=12) ttcAAATTTGCCGTGCAAATTTGAAGAGGTTTTAACTACATGTTATGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCTTTGCACGCAGGAGGTCTGCGGTTCGATCCCGCATAGCTCCACCATCTCTGTAGTGGTTAAATAAAAAATACt > 1:925155/3‑150 (MQ=12) ttcAAATTTGCCGTGCAAATTTGAAGAGGTTTTAACTACATGTTATGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCTTTGCACGCAGGAGGTCTGCGGTTCGATCCCGCATAGCTCCACCATCTCTGTAGTGGTTAAATAAAAAATACt > 1:739179/3‑150 (MQ=12) tcAAATTTGCCGTGCAAATTTGAAGAGGTTTTAACTACATGTTATGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCTTTGCACGCAGGAGGTCTGCGGTTCGATCCCGCATAGCTCCACCATCTCTGTAGTGGTTAAATAAAAAATACtt > 1:2139405/2‑150 (MQ=12) ccAAATTTGCACGGCAAATTTGAAGAGGTTTTAACTACATGTTATGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCTTTGCACGCAGGAGGTCTGCGGTTCGATCCCGCATAGCTCCACCATCTCTGTAGTGGTTAAATAAAAAATACtt > 1:1853732/1‑150 (MQ=32) cAAATTTGCACGGCAATTTTGAAGAGTTTTTAACTACATGTTATGGGGCTATAGCTCAGCTGGGAGAGCGCATGCTTTGCACGCAGGAGGTCTGCGGTTCGATCCCGCATAGCTCCACCATCTCTGTAGTGGTTAAATAAAAAATACTTc < 1:848544/150‑1 (MQ=12) cAAATTTGCACGGCAAATTTGAAGTGGTTTTAACTACATGTTATGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCTTTGCACGCAGGAGGTCTGCGGTTCGATCCCGCATTGCTCCACCATCTCTGTAGTGGTTAAATAAAAAATACTTc > 1:1372237/1‑150 (MQ=17) cAAATTTGCACGGCAAATTTGAAGAGGTTTTAACTACATGTTATGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCTTTGCACGCAGGATGTCTTCGGTTCGATCCCGCATAGCTCCACCATCTCTGTAGTGGTTAAATAAAAAATACTTc > 2:1364554/1‑150 (MQ=17) cAAATTTGCACGGCAAATTTGAAGAGGTTTTAACTACATGTTATGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCTTTGCACGCAGGAGGTCTGCGGTTCGATCCCGCATAGCTCCACCATCTCTGTAGTGGTTAAATAAAAAATACTTc > 1:949496/1‑150 (MQ=32) aaaTTTGCACGGCAAATTTGAAGAGGTTTTAACTACATGTTATGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCTTTGCACGCAGGAGGTCTGCGGTTCGATCCCGCATAGCTCCACCATCTCTGTAGTGGTTAAATAAAAAATACTTCa > 1:476446/1‑150 (MQ=32) tttGCACGGCAAATTTGAAGAGGTTTTAACTACATGTTATGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCTTTGCACGCAGGAGGTCTGCGGTTCGATCCCGCATAGCTCCACCATCTCTGTAGTGGTTAAATAAAAAATACTTCagag < 1:1289256/150‑1 (MQ=32) ttGCACGGCAAATTTGAAGAGGTTTTAACTACATGTTATGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCTTTGCACGCAGGAGGTCTGCGGTTCGATCCCGCATAGCTCCACCATCTCTGTAGTGGTTAAATAAAAAATACTTCAGAgt > 2:570610/1‑150 (MQ=32) tGCCCGGCAAATTTGAAGAGGTTTTAACTACATGTTATGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCTTTGCACTCAGGAGGTCTGCGGTTCGATCCCGCATAGCTCCACCATCTCTGTAGTGGTTAAATAAAAAATCCTTCAGAgtg > 1:142712/1‑150 (MQ=12) cACGGCAAATTTGAAGAGGTTTTAACTACATGTTATGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCTTTGCACGCAGGAGGTCTGCGGTTCGATCCCGCATAGCTCCACCAGCTCTGTAGTGGTTAAATAAAAAATACTTCAGAGTGTa < 1:1195106/150‑1 (MQ=17) cACGGCAAATTTGAAGAGGTTTTAACTACATGTTATGGGGCTATAGCTCAGCAGGGAGAGCGCCTGCTTTGCACGCATGATGTCTTCGGTTCGATCCCGCATTGCTACACCAGCTCTGTAGTGGTTAAATAAAAAATACTTCAGAGTGTa > 1:1349534/1‑150 (MQ=11) aCGGCAAATTTGAAGAGGTTTTAACTACATGTTATGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCTTTGCACGCAGGAGGTCTGCGGTTCGATCCCGCATAGCTCCACCATCTCTGTAGTGGTTAAATAAAAAATACTTCAGAGTGTAc < 1:1898822/150‑1 (MQ=32) | AAGTCCACTCAGGCCTACCAAATTTGCACGGCAAATTTGAAGAGGTTTTAACTACATGTTATGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCTTTGCACGCAGGAGGTCTGCGGTTCGATCCCGCATAGCTCCACCATCTCTGTAGTGGTTAAATAAAAAATACTTCAGAGTGTAC > NC_000913/225438‑225614 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |