| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 225,467 | 1 | . | T | 25.0% | 20.0 / 18.2 | 20 | intergenic (+10/‑33) | ileV/alaV | tRNA‑Ile/tRNA‑Ala |
| Reads supporting (aligned to +/- strand): ref base . (7/8); new base T (5/0); total (12/8) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 5.47e-02 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.87e-01 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
CTTTGCAGTGCTCACACAGATTGTCTGATGAAAATGAGCAGTAAAACCTCTACAGGCTTGTAGCTCAGGTGGTTAGAGCGCACCCCTGATAAGGGTGAGGTCGGTGGTTCAAGTCCACTCAGGCCTACCAAATTTGCACG‑GCAAATTTGAAGAGGTTTTAACTACATGTTATGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCTTTGCACGCAGGAGGTCTGCGGTTCGATCCCGCATAGCTCCACCATCTCTGTAGTGGTTAAATAAAAAATACTTCAGAGTGTACCT > NC_000913/225328‑225616 | ctgtGCAGTGCTCACACAGATTGTCTGATGAAAATGAGCAGTAAAACCTCTACAGGCTTGTAGCTCAGGTGGTTAGAGCGCACCCCTGATAAGGGTGAGGTCGGTGGTTCAAGTCCACTCAGGCCTACCAAATTTGCCGT‑GCAAATTTgg < 1:301729/147‑2 (MQ=0) aaGTCCACTCAGGCCTACCAAATTTGCACG‑GCAAATTTGAAGAGGTTTTAACTACATGTTATGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCTTTGCACGCAGGAGGTCTGCGGTTCGATCCCGCATAGCTCCACCATCTCTGTAGTgg > 1:1201424/1‑150 (MQ=11) aGTCCACTCAGGGCTACCGAATTTGCACG‑GCAAATTTGGAGAGGTTTTAACTACATGTTATGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCTTTGCACGCAGGAGGTCTGCGGTTCGATCCCGCATAGCTCCACCATCTCTGTAGTGGt < 2:1226389/150‑1 (MQ=9) gTCCACTCAGGCCTACCAAATTTGCACG‑GCAAATTTGAAGAGGTTTTAACTACATGTTATGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCTTTGCACGCAGGAGGTCTGCGGTTCGATCCCGCATAGCTCCACCATCTCTGTAGTGGtt < 1:9251/150‑1 (MQ=32) aGGCCTACCAAATTTGCACG‑GCAAATTTGAAGAGGTTTTAACTACATGTTATGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCTTTGCACGCAGGAGGTCTGCGGTTAGATCCCGCATAGCTCCACCATCTCTGTAGTGGTTAAATaaaa < 1:517745/150‑1 (MQ=17) ctcttcAAATTTGCCGT‑GCAAATTTGAAGAGGTTTTAACTACATGTTATGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCTTTGCACGCAGGAGGTCTGCGGTTCGATCCCGCATAGCTCCACCATCTCTGTAGTGGTTAAATAAAAAAt > 2:1055305/6‑150 (MQ=11) ctcttcAAATTTGC‑CGTGCAAATTTGGAGAGGTTTTAACTACATGTTATGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCTTTGCACGCAGGAGGTCTGCGGTTCGATCCCGCATAGCTCCACCATCTCTGTAGTGGTTAAATAAAAAAt > 3:88840/6‑150 (MQ=11) ccTACCAAATTTGCACG‑GCAAATTTGAAGAGGTTTTAACTACATGTTATGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCTTTGCACGCAGGAGGTCTGCGGTTCGATCCCGCATAGCTCCACCATCTCTGTAGTGGTTACATAAAAAAt > 4:90069/1‑150 (MQ=17) cTACCAAATTTGCACG‑GCAAATTTGAAGAGGTTTTAACTACATGTTATGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCTTTGCACGCAGGAGGTCTGCGGTTCGATCCCGCATAGCTCCACCATCTCTGTAGTGGTTAAATAAAAAATa < 3:488040/150‑1 (MQ=32) ttcAAATTTGC‑CGTGCAAATTTGAAGAGGTTTTAACTACATGTTATGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCTTTGCACGCAGGAGGTCTGCGGTTCGATCCCGCATAGCTCCACCATCTCTGTAGTGGTTAAATAAAAAATACt > 2:958520/3‑150 (MQ=255) tcAAATTTGC‑CGTGCAAATTTGAAGAGGTTTTAACTACATGTTATGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCTTTGCACGCAGGAGGTCTGCGGTTCGATCCCGCATAGCTCCACCATCTCTGTAGTGGTTAAATAAAAAATACtt > 2:1202488/2‑150 (MQ=255) cAAATTTGCACG‑GCAAATTTGAAGAGGTTTTAACTACATGTTATGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCTTTGCACGCAGGAGGTCTGCGGTTCGATCCCGCATAGCTCCACCATCTCTGTAGTGGTTAAATACAAAATACTTc > 2:781771/1‑150 (MQ=17) cAAATTTGCACG‑GCAAATTTGAAGAGGTTTTAACTACATGTTATGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCTTTGCACGCAGGAGGTCTGCGGTTCGATCCCGCATAGCTCCACCATCTCTGTAGTGGTTAAATAAAAAATACTTc > 2:107300/1‑150 (MQ=32) cAAATTTGC‑CGTGCAAATTTGAAGAGGTTTTAACTACATGTTATGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCTTTGCACGCAGGAGGTCTGCGGTTCGATCCCGCATAGCTCCACCATCTCTGTAGTGGTTAAATAAAAAATACTTc > 2:148367/1‑150 (MQ=12) cAAATTTGC‑CGTGCAAATTTGAAGAGGTTTTAACTACATGTTATGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCTTTGCACGCAGGAGGCCTGCGGTTCGATCCCGCATCGCTCCACCATCTCTGTAGTGGTTAAATAAAAAATACTTc > 1:57318/1‑150 (MQ=11) aaaTTTGCACG‑GCAAATTTGAAGAGGTTTTAACTACATGTTATGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCTTTGCACGCAGGAGGTCTGCGGTTCGATCCCGCATAGCTCCACCATCTCTGTAGTGGTTAAATAAAAAATACTTCa > 1:223456/1‑150 (MQ=32) aTTTGCACG‑GCAAATTTGAAGAGGTTTTAACTACATGTTATGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCTTTGCACGCAGGAGGTCTGCGGTTCGATCCCGCATAGCTCCACCATCTCTGTAGTGGTTAAATAAAAAATACTTCaga < 1:1193786/150‑1 (MQ=32) ttGCACG‑GCAAATTTGAAGAGGTTTTAACTACATGTTATGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCTTTGCACGCAGGAGGTCTGCGGTTCGATCCCGCATAGCTCCACCATCTCTGTAGTGGTTAAATAAAAAATACTTCAGAgt > 2:1066852/1‑150 (MQ=32) aCG‑GCAAATTTGAAGAGGTTTTAACTACATGTTATGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCTTTGCACGCAGGAGGTCTGCGGTTCGATCCCGCATAGCTCCACCATCTCTGTAGTGGTTAAATAAAAAATACTTCAGAGTGTAc < 2:173828/150‑1 (MQ=32) aCG‑GCAAATTTGAAGAGGTTTTAACTACATGTTATGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCTTTGCACGCAGGAGGTCTGCGGTTCGATCCCGCATAGCTCCACCATCTCTGTAGTGGTTAAATAAAAAATACTTCAGAGTGTAc < 2:173934/150‑1 (MQ=32) g‑gCCAATTTGAAGAGGTTTTAACTACATGTTATGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCTTTGCACGCAGGAGGTCTGCGGTTCGATCCCGCATAGCTCCACCATCTCTGTAGTGGTTAAATAAAAAATACTTCAGAGTGTACCt < 3:437119/150‑1 (MQ=17) | CTTTGCAGTGCTCACACAGATTGTCTGATGAAAATGAGCAGTAAAACCTCTACAGGCTTGTAGCTCAGGTGGTTAGAGCGCACCCCTGATAAGGGTGAGGTCGGTGGTTCAAGTCCACTCAGGCCTACCAAATTTGCACG‑GCAAATTTGAAGAGGTTTTAACTACATGTTATGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCTTTGCACGCAGGAGGTCTGCGGTTCGATCCCGCATAGCTCCACCATCTCTGTAGTGGTTAAATAAAAAATACTTCAGAGTGTACCT > NC_000913/225328‑225616 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |