| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 4,037,291 | 0 | T | C | 48.1% | 14.4 / 37.3 | 27 | noncoding (32/76 nt) | alaT | tRNA‑Ala |
| Reads supporting (aligned to +/- strand): ref base T (7/7); new base C (11/2); total (18/9) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.03e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.38e-01 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
GTTAGAGCGCACCCCTGATAAGGGTGAGGTCGGTGGTTCAAGTCCACTCAGGCCTACCAAATTTGCACGGCAAATTTGAAGAGGTTTTAACTACATGTTATGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCTTTGCACGCAGGAGGTCTGCGGTTCGATCCCGCATAGCTCCACCATCTCTGTAGTGATTAAATAAAAAATACTTCAGAGTGTA > NC_000913/4037159‑4037373 | aggtctgtagctcagttggttagagcgcacccctgataagggtgaggtcggcagttcaaatctgcccagacctaccattacatggttcagccgtagaatacGGGGCCATAGCTCAGCTGGGAGAGCGCCTGCTTTGCACGc > 2:2109384‑M1/102‑141 (MQ=255)aggtctgtagctcagttggttagagcgcacccctgataagggtgaggtcggcagttcaaatctgcccagacctaccattacatggttcagccgtagaatacGGGGCCATAGCTCAGCTGGGAGAGCGCCTGCTTTGCACGc > 2:1945640‑M1/102‑141 (MQ=255)aggtctgtagctcagttggttagagcgcacccctgataagggtgaggtcggcagttcaaatctgcccagacctaccattacatggttcagccgtagaatacGGGGCCATAGCTCAGCTGGGAGAGCGCCTGCTTTGCACGc > 2:1945636‑M1/102‑141 (MQ=255)aggtctgtagctcagttggttagagcgcacccctgataagggtgaggtcggcagttcaaatctgcccagacctaccattacatggttcagccgtagaatacGGGGCCATAGCTCAGCTGGGAGAGCGCCTGCTTTGCACGc > 2:1944634‑M1/102‑141 (MQ=255) ttagagcgcacccctgataagggtgaggtcggcagttcaaatctgcccagacctaccattacatggttcagccgtagaatacGGGGCCATAGCTCAGCTGGGAGAGCGCCTGCCTTGCACGCAGGAGGTCAGCGGTTCGAt > 1:685297‑M1/83‑141 (MQ=255) gagcgcacccctgataagggtgaggtcggcagttcaaatctgcccagacctaccattacatggttcagccgtagaatacGGGGCCATAGCTCAGCTGGGAGAGCGCCTGCCTTGCACGCAGGAGGTCAGCGGTTCGATccc < 1:2019618‑M1/62‑1 (MQ=255) gagcgcacccctgataagggtgaggtcggcagttcaaatctgcccagacctaccattacatggttcagccgtagaatacGGGGCCATAGCTCAGCTGGGAGAGCGCCTGCCTTGCACGCAGGAGGTCAGCGGTTCGATccc < 1:1098393‑M1/62‑1 (MQ=255) ccctgataagggtgaggtcggcagttcaaatctgcccagacctaccattacatggttcagccgtagaatacGGGGCCATAGCTCAGCTGGGAGAGCGCCTGCCTTGCACGCAGGAGGTCAGCGGTTCGATCCCGCATAGCt > 2:1439943‑M1/72‑141 (MQ=255) gaggtcggcagttcaaatctgcccagacctaccattacatggttcagccgtagaatacGGGGCCATAGCTCAGCTGGGAGAGCGCCTGCCTTGCACGCAGGAGGTCAGCGGTTCGATCCCGcttggctccaccactttcgc > 2:793747‑M1/59‑134 (MQ=255) gttcaaatctgcccagacctaccattacatggttcagccgtagaatacGGGGCCATAGCTCAGCTGGGAGAGCGCCTGCCTTGCACGCAGGAGGTCTGCGGTTCGATCCCGCATAGCTCCACCATCTCTGTAGTGATtaaa > 2:1093341‑M1/49‑141 (MQ=255) gttcaaatctgcccagacctaccattacatggttcagccgtagaatacGGGGCCATAGCTCAGCTGGGAGAGCGCCTGCCTTGCACGCAGGAGGTCTGCGGTTCGATCCCGCATAGCTCCACCATCTCTGTAGTGATtaaa > 2:1093246‑M1/49‑141 (MQ=255) gttcaaatctgcccagacctaccattacatggttcagccgtagaatacGGGGCCATAGCTCAGCTGGGAGAGCGCCTGCCTTGCACGCAGGAGGTCTGCGGTTCGATCCCGCATAGCTCCACCATCTCTGTAGTGATtaaa > 2:1093207‑M1/49‑141 (MQ=255) gttcaaatctgcccagacctaccattacatggttcagccgtagaatacGGGGCCATAGCTCAGCTGGGAGAGCGCCTGCCTTGCACGCAGGAGGTCTGCGGTTCGATCCCGCATAGCTCCACCATCTCTGTAGTGATtaaa > 2:1093202‑M1/49‑141 (MQ=255) gttcaaatctgcccagacctaccattacatggttcagccgtagaatacGGGGCCATAGCTCAGCTGGGAGAGCGCCTGCCTTGCACGCAGGAGGTCTGCGGTTCGATCCCGCATAGCTCCACCATCTCTGTAGTGATtaaa > 2:1093201‑M1/49‑141 (MQ=255) caaatctgcccagacctaccattacatggttcagccgtagaatacGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCTTTGCACGCAGGAGGTCTGCGGTTCGATCCCGCATAGCTCCACCATCTCTGTAGTGATtaaataa > 2:356099‑M1/46‑141 (MQ=255) caaatctgcccagacctaccattacatggttcagccgtagaatacGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCTTTGCACGCAGGAGGTCTGCGGTTCGATCCCGCATAGCTCCACCATCTCTGTAGTGATtaaataa > 2:2152620‑M1/46‑141 (MQ=255) tctgcccagacctaccattacatggttcagccgtagaatacGGGGCCATAGCTCAGCTGGGAGAGCGCCTGCCTTGCACGCAGGAGGTCAGCGGTTCGATCCCGCTTGGCTCCACCATCTCTGTAGTGATTAAATaaaaaa > 2:36977‑M1/42‑141 (MQ=255) tctgcccagacctaccattacatggttcagccgtagaatacGGGGCCATAGCTCAGCTGGGAGAGCGCCTGCCTTGCACGCAGGAGGTCAGCGGTTCGATCCCGCTTGGCTCCACCATCTCTGTAGTGATTAAATaaaaaa > 2:1013856‑M1/42‑141 (MQ=255) tctgcccagacctaccattacatggttcagccgtagaatacGGGGCCATAGCTCAGCTGGGAGAGCGCCTGCCTTGCACGCAGGAGGTCAGCGGTTCGATCCCGCTTGGCTCCACCATCTCTGTAGTGATTAAATaaaaaa > 2:231712‑M1/42‑141 (MQ=255) atgtGCACGGCAAATTTGAAGAGGTTTTAACTACATGTTATGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCTTTGCACGCAGGAGGTCTGCGGTTCGATCCCGCATAGCTCCACCATCTCTGTAGTGATTAAATaaaaaa < 2:1237302/138‑1 (MQ=1) atgtGCACGGCAAATTTGAAGAGGTTTTAACTACATGTTATGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCTTTGCACGCAGGAGGTCTGCGGTTCGATCCCGCATAGCTCCACCATCTCTGTAGTGATTAAATaaaaaa < 2:496222/138‑1 (MQ=1) atgtGCACGGCAAATTTGAAGAGGTTTTAACTACATGTTATGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCTTTGCACGCAGGAGGTCTGCGGTTCGATCCCGCATAGCTCCACCATCTCTGTAGTGATTAAATaaaaaa < 2:1783632/138‑1 (MQ=1) atgtGCACGGCAAATTTGAAGAGGTTTTAACTACATGTTATGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCTTTGCACGCAGGAGGTCTGCGGTTCGATCCCGCATAGCTCCACCATCTCTGTAGTGATTAAATaaaaaa < 2:1678113/138‑1 (MQ=1) ctgcccagacctaccattacatggttcagccgtagaatacGGGGCCATAGCTCAGCTGGGAGAGCGCCCGCTTTGCACGCAGGAGGTCTGCGGTTCGATCCCGCATAGCTCCACCATCTCTGTAGTGATTAAATAAAAAAt > 2:734623‑M1/41‑141 (MQ=255) ttggAAGAGGTTTTAACTACATGTTATGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCTTTGCACGCAGGAGGTCTGCGGTTCGATCCCGCATAGCTCCACCATCTCTGTAGTGATTAAATAAAAAATACTTCAGAGTGTa < 2:1481037/138‑1 (MQ=1) ttagAAGAGGATTTAACTACATGTTATGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCTTTGCACGCAGGAGGTCAGCGGTTCGATCCCGCTTGGCTCCACCATCTCTGTAGTGATTAAATCAAAAATACTTCAGAGTGTa < 2:428442/138‑1 (MQ=1) tcagacgagGTTTTAACTACATGTTATGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCTGTGCACGCAGGAGGTCTGCGGTTCGATCCCGCATAGCTCCACCATCTCTGTAGTGATTAAATAAAAAATACTTCAGAGTGTa < 1:489901/135‑1 (MQ=0) | GTTAGAGCGCACCCCTGATAAGGGTGAGGTCGGTGGTTCAAGTCCACTCAGGCCTACCAAATTTGCACGGCAAATTTGAAGAGGTTTTAACTACATGTTATGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCTTTGCACGCAGGAGGTCTGCGGTTCGATCCCGCATAGCTCCACCATCTCTGTAGTGATTAAATAAAAAATACTTCAGAGTGTA > NC_000913/4037159‑4037373 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |