| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 1,527,602 | 0 | G | C | 29.2% | 35.0 / 13.7 | 26 | intergenic (+450/‑300) | yncH/rhsE | IPR020099 family protein/pseudogene, Rhs family |
| Reads supporting (aligned to +/- strand): ref base G (10/7); new base C (0/7); total (10/16) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.88e-02 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.78e-01 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
ACTATTAAGCTAGTGTTTAACGACGTTAAGATGGCAATGTGCAGATCATGTTATCAGGGTGGGGCGGTAGCGTATATTCCTTATACAAACCATCATTCATTTCAATGATTGTTATTGAATTATTTATTGTT > NC_000913/1527539‑1527669 | aCTATTAAGCTAGTGTTTAACGACGTTAAGATGGCAATGTGCAGATCATGTTATCAGGGTGGGGCGGt < 1:89200/68‑1 (MQ=255) gCTAGTGTTTAACGACGTTAAGATGGCAATGTGCAGATCATGTTATCAGGGTGGGGCGGTAGCGtata > 1:145970/1‑68 (MQ=255) gCTAGTGTTTAACGACGTTAAGATGGCAATGTGCAGATCATGTTATCAGGGTGGGGCGGTAGCGtata > 1:48821/1‑68 (MQ=255) gCTAGTGTTTAACGACGTTAAGAAGGCAATGTGCAGATCATGTTATCAGGGTGGGGCGGTAGCGtata > 2:291668/1‑68 (MQ=255) tgtgtTTATCGACGTTAAGATGGCATTGTGCAGATCATGTTATCAGGGTGGGCCGGTAGCGTATATTc < 2:424632/67‑1 (MQ=255) aGTGTTTAACGACGTTAAGATGGCAATGTGCAGATCATGTTATCAGGGTGGGGCGGTAGCGTATATTc < 2:451/68‑1 (MQ=255) aGTGTTTAACGACGTTAAGATGGCAATGTGCAGATCATGTTATCAGGGTGGGCCGGTAGCGTATATTc < 1:105270/68‑1 (MQ=255) aGTGTTAAACGACGTTAAGATGGCAATGTGCAGATCATGTTATCAGGGGGGGCCGGTAGCGTATATTc < 2:729687/68‑1 (MQ=255) ttAACGACGTTAAGATGGCAATGTGCAGATCATGTTATCAGGGTGGGCCGGTAGCGTATATTCCTtat < 1:654377/68‑1 (MQ=255) aaGATGGCAATGTGCAGATCATGTTATCAGGGTGGGGCGGTAGCGTATATTCCTTATACAAACcatca < 2:48821/68‑1 (MQ=255) gATGGCAATGTGCAGATCATGTTATCAGGGTGGGGCGGTAGCGTATATTCCTTATACAAACCAtcatt > 1:600364/1‑68 (MQ=255) tGGCAATGTGCAGATCATGTTATCAGGGTGGGGCGGTAGCGTATATTCCTTATACAAACCAtcattca < 2:748385/68‑1 (MQ=255) aaTGTGCAGATCATGTTATCAGGGTGGGGCGGTAGCGTATATTCCTTATACAAACCATCATTCATTTc > 1:664320/1‑68 (MQ=255) gCAGATCATGTTATCAGGGTGGGCCGGTAGCGTATATTCCTTATACAAACCATCATTCATTTCAATGa < 2:140586/68‑1 (MQ=255) tCATGTTATCAGGGTGGGGCGGTAGCGTATATTCCTTATACAAACCATCATTCATTTCAATGATTGtt > 1:134890/1‑68 (MQ=255) tCATGTTATCAGGGTGGGGCGGTAGCGTATATTCCTTATACAAACCATCATTCATTTCAATGATTGtt < 1:307331/68‑1 (MQ=255) tCATGTTATCAGGGTGGGGCGGTAGCGTATATTCCTTATACAAACCATCATTCATTTCAATGATTGtt > 2:307330/1‑68 (MQ=255) cAGGGTGGGGCGGTAGCGTATATTCCTTATACAAACCATCATTCATTTCAATGa > 1:79470/1‑54 (MQ=255) cAGGGTGGGCCGGTAGCGTATATTCCTTATACAAACCATCATTCATTTCAATGa < 2:79470/54‑1 (MQ=255) aGGGTGGGTCGGTAGCGTATATTCCTTATACAAACCATCATTCATTTCAATGATTGTTATTGAattat < 2:193781/68‑1 (MQ=255) aGGGTGGGTCGGTAGCGTATATTCCTTATACAAACCATCATTCATTTCAATGATTGTTATTGAattat < 1:104710/68‑1 (MQ=255) aGGGTGGGGCGGTAGCGTATATTCCTTATACAAACCATCATTCATTTCAATGATTGTTATTGAattat < 2:664320/68‑1 (MQ=255) aGGGTGGGCCGGTAGCGTATATTCCTTATACAAACCATCATTCATTTCAATGATTGTTATTGAattat < 1:162308/68‑1 (MQ=255) ggTGGGGCGGTAGCGTATATTCCTTATACAAACCATCATTCATTTCAATGATTGTTATTGa > 1:161526/1‑61 (MQ=255) ggTGGGGCGGTAGCGTATATTCCTTATACAAACCATCATTCATTTCAATGATTGTTATTGa < 2:161526/61‑1 (MQ=255) gTGGGGCGGTAGCGTATATTCCTTATACAAACCATCATTCATTTCAATGATTGTTATTGAAttattta > 1:505377/1‑68 (MQ=255) gCGGTAGCGTATATTCCTTATACAAACCATCATTCATTTCAATGATTGTTATTGAATTATTTAttgtt > 2:651189/1‑68 (MQ=255) | ACTATTAAGCTAGTGTTTAACGACGTTAAGATGGCAATGTGCAGATCATGTTATCAGGGTGGGGCGGTAGCGTATATTCCTTATACAAACCATCATTCATTTCAATGATTGTTATTGAATTATTTATTGTT > NC_000913/1527539‑1527669 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |