Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 1,527,602 | 0 | G | C | 28.1% | 67.2 / 16.7 | 32 | intergenic (+450/‑300) | yncH/rhsE | IPR020099 family protein/pseudogene, Rhs family |
Reads supporting (aligned to +/- strand): ref base G (17/6); new base C (0/9); total (17/15) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.78e-04 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
TTAAGCTAGTGTTTAACGACGTTAAGATGGCAATGTGCAGATCATGTTATCAGGGTGGGGCGGTAGCGTATATTCCTTATACAAACCATCATTCATTTCAATGATTGTTATTGAATTATTTATTGTT > NC_000913/1527543‑1527669 | ttAAGCTTGTGTTTAACGACGTTAAGATGGCAATGTGCAGATCATGTTATCAGGGTGGGGCGGTAGCg < 2:374455/68‑1 (MQ=255) ttAAGCTAGTGTTTAACGACGTTAAGATGGCAATGTGCAGATCATGTTATCAGGGTGGGCCGGTAGCg < 2:473821/68‑1 (MQ=255) tAAGCTAGTGTTTAACGACGTTAAGATGGCAATGTGCAGATCATGTTATCAGGGTGGGGCGGTAGCGt > 1:804212/1‑68 (MQ=255) tAAGCTAGTGTTTAACGACGTTAAGATGGCAATGTGCAGATCATGTTATCAGGGTGGGCCGGTAGCGt < 1:533402/68‑1 (MQ=255) gCTAGTGTTTAACGACGTTAAGATGGCAATGTGCAGATCATGTTATCAGGGTGGGGCGGTAGCg > 2:777209/1‑64 (MQ=255) gCTAGTGTTTAACGACGTTAAGATGGCAATGTGCAGATCATGTTATCAGGGTGGGGCGGTAGCGtata > 2:105579/1‑68 (MQ=255) gCTAGTGTTTAACGACGTTAAGATGGCAATGTGCAGATCATGTTATCAGGGTGGGGCGGTAGCGtata > 2:1115748/1‑68 (MQ=255) gCTAGTGTTTAACGACGTTAAGATGGCAATGTGCAGATCATGTTATCAGGGTGGGGCGGTAGCGtata > 2:206348/1‑68 (MQ=255) gCTAGTGTTTAACGACGTTAAGATGGCAATGTGCAGATCATGTTATCAGGGTGGGGCGGTAGCGtata > 2:202114/1‑68 (MQ=255) gCTAGTGTTTAACGACGTTAAGATGGCAATGTGCAGATCATGTTATCAGGGTGGGCCGGTAGCg < 1:777209/64‑1 (MQ=255) gCTAGTGTTTAACGACGTTAAGATGACAATGTGCAGATCATGTTATCAGGGTGGGGCGGTAGCGtata > 1:816697/1‑68 (MQ=255) tAGTGTTTAACGACGTTAAGATGGCAATGTGCAGATCATGTTATCAGGGTGGGGCGGTAGCGTATAtt > 2:1030557/1‑68 (MQ=255) tAGTGTTTAACGACGTTAAGATGGCAATGTGCAGATCATGTTATCAGGGTGGGGCGGTAGCGTATAtt > 2:63758/1‑68 (MQ=255) tgtgtTTAACTACTTTAATATGGCAATGTGCATATCATGTTATCTGGTTGGTCCGGTAGCGTATATTc < 1:206348/67‑1 (MQ=255) aGTGTTTAACGACGTTAAGATGGCAATGTGCAGATCATGTTATCAGGGTGGGGCTGTAGCGTATATTc < 2:804212/68‑1 (MQ=255) aGTGTTTAACGACGTTAAGATGGCAATGTGCAGATCATGTTATCAGGGTGGGCCGGTAGCGTATATTc < 1:1115750/68‑1 (MQ=255) aGTGTTTAACGACGTTAAGATGGCAATGTGCAGATCATGTTATCAGGGTGGGCCGGTAGCGTATATTc < 1:998429/68‑1 (MQ=255) gacgTTAAGATGGCAATGTGCAGATCATGTTATCAGGGTGGGGCGGTAGCGTATATTCCTTATACaaa > 1:173225/1‑68 (MQ=255) cgTTAAGATGGCAATGTGCAGATCATGTTATCAGGGTGGGGCGGTAGCGTATATTCCTTATACAAAcc > 2:837888/1‑68 (MQ=255) cgTTAAGATGGCAATGTGCAGATCATGTTATCAGGGTGGGGCGGTAGCGTATATTCCTTATACAAAcc < 1:202114/68‑1 (MQ=255) cgTTAAGATGGCAATGTGCAGATCATGTTATCAGGGTGGGGCGGTAGCGTATATTCCTTATACAAAcc < 1:105579/68‑1 (MQ=255) tAAGATGGCAATGTGCAGATCATGTTATCAGGGTGGGGCGGTAGCGTATATTCCTTATACAAACcatc < 1:837888/68‑1 (MQ=255) gATGGCAATGTGCAGATCATGTTATCAGGGTGGGGCGGTAGCGTATATTCCTTATACAAACCAtcatt > 2:1056774/1‑68 (MQ=255) tgCAGATCATGTTATCAGGGTGGGGCGGTAGCGTATATTCCTTATACAAACCATCATTCATTTCAATg > 1:1062173/1‑68 (MQ=255) gCAGATCATGTTATCAGGGTGGGGCGGTAGCGTATATTCCTTATACAAACCATCATTCATTTCAATGa > 1:364256/1‑68 (MQ=255) aGATCATGTTATCAGGGTGGGCCGGTAGCGTATATTCCTTATACAAACCATCATTCATTTCAATGAtt < 2:518746/68‑1 (MQ=255) gATCATGTTATCAGGGTGGGGCGGTAGCGTATATTCCTTATACAAACCATCATTCATTTCAATGATTg > 1:1114707/1‑68 (MQ=255) tATCAGGGTGGGGCGGTAGCGTATATTCCTTATACAAACCATCATTCATTTCAATGATTGTTATTGaa > 1:332241/1‑68 (MQ=255) aTCAGGGTGGGGCGGTAGCGTATATTCCTTATACAAACCATCATTCATTTCAATGATTGTTATTGAat > 2:393150/1‑68 (MQ=255) aTCAGGGTGGGGCGGTAGCGTATATTCCTTATACAAACCATCATTCATTTCAATGATTGTTATTGAat < 1:1056775/68‑1 (MQ=255) aGGGTGGGCCGGTAGCGTATATTCCTTATACAAACCATCATTCATTTCAATGATTGTTATTGAattat < 2:149698/68‑1 (MQ=255) aGGGTGGGCCGGTAGCGTATATTCCTTATACAAACCATCATTCATTTCAATGATTGTTATTGAattat < 1:1030558/68‑1 (MQ=255) gCGGTAGCGTATATTCCTTATACAAACCATCATTCATTTCAATGATTGTTTTTGAATTATTTAttgtt > 1:710222/1‑68 (MQ=255) gCGGTAGCGTATATTCCTAATACACACCATCAATCATTTCAATGATTGTTATTGAATTATTTAttgtt > 2:1122304/1‑68 (MQ=255) | TTAAGCTAGTGTTTAACGACGTTAAGATGGCAATGTGCAGATCATGTTATCAGGGTGGGGCGGTAGCGTATATTCCTTATACAAACCATCATTCATTTCAATGATTGTTATTGAATTATTTATTGTT > NC_000913/1527543‑1527669 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |