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.6% | 88.4 / 20.4 | 42 | intergenic (+450/‑300) | yncH/rhsE | IPR020099 family protein/pseudogene, Rhs family |
Reads supporting (aligned to +/- strand): ref base G (24/6); new base C (0/12); total (24/18) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.68e-06 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
CCACTATTAAGCTAGTGTTTAACGACGTTAAGATGGCAATGTGCAGATCATGTTATCAGGGTGGGGCGGTAGCGTATATTCCTTATACAAACCATCATTCATTTCAATGATTGTTATTGAATTATTTATTGTT > NC_000913/1527537‑1527669 | ccACTATTAAGCTAGTGTTTAACGACGTTAAGATGGCAATGTGCAGATCATGTTATCAGGGTGGGGCg > 2:506172/1‑68 (MQ=255) aCTATTAAGCTAGTGTTTAACGACGTTAAGATGGCAATGTGCAGATCATGTTATCAGGGTGGGCCGGt < 2:491627/68‑1 (MQ=255) ttAAGCTAGTGTTTAACGACGTTAAGATGGCAATGTGCAGATCATGTTATCAGGGTGGGCCGGTAGCg < 1:868916/68‑1 (MQ=255) taaaGCTAGTGTTTAACGACGTTAAGATGGCAATGTGCAGATCATGTTATCAGGGTGGGCCGGTAGCg < 1:642946/66‑1 (MQ=255) tAAGCTAGTGTTTAACGACGTTAAGATGGCAATGTGCAGATCATGTTATCAGGGTGGGGCGGTAGCGt > 2:206595/1‑68 (MQ=255) gCTAGTGTTTAACGACGTTAAGATGGCAATGTGCAGATCATGTTATCAGGGTGGGGCGGTAGCGtata > 1:266048/1‑68 (MQ=255) gCTAGTGTTTAACGACGTTAAGATGGCAATGTGCAGATCATGTTATCAGGGTGGGGCGGTAGCGtata > 2:481641/1‑68 (MQ=255) gCTAGTGTTTAACGACGTTAAGATGGCAATGTGCAGATCATGTTATCAGGGTGGGGCGGTAGCGtata > 1:44814/1‑68 (MQ=255) gCTAGTGTTTAACGACGTTAAGATGGCAATGTGCAGATCATGTTATCAGGGTGGGGCGGTAGCGtata > 1:470024/1‑68 (MQ=255) gCTAGTGTTTAACGACGTTAAGATGGCAATGTGCAGATCATGTTATCAGGGTGGGGCGGTAGCGtata > 2:717741/1‑68 (MQ=255) gCTAGTGTTTAACGACGTTAAGATGGCAATGTGCAGATCATGTTATCAGGGTGGGGCGGTAGCGtata > 1:548186/1‑68 (MQ=255) gCTAGTGTTTAACGACGTTAAGATGGCAATGTGCAGATCATGTTATCAGGGTGGGGCGGTAGCGtata > 1:574563/1‑68 (MQ=255) gCTAGTGTTTAACGACGTTAAGATGGCAATGTGCAGATCATGTTATCAGGGTGGGGCGGTAGCGtata > 2:175027/1‑68 (MQ=255) gCTAGTGTTTAACGACGTTAAGATGGCAATGTGCAGATCATGTTATCAGGGTGGGGCGGTAGCGtata > 1:724184/1‑68 (MQ=255) tAGTGTTTAACGACGTTAAGATGGCAATGTGCAGATCATGTTATCAGGGTGGGGCGGTAGCGTATAtt > 1:206097/1‑68 (MQ=255) aGTGTTTAACGACGTTAAGATGGCAATGTGCAGATCATGTTATCAGGGTGGGGCGGTAGCGTATATTc < 2:44814/68‑1 (MQ=255) aGTGTTTAACGACGTTAAGATGGCAATGTGCAGATCATGTTATCAGGGTGGGGCGGTAGCGTATATTc < 1:742292/68‑1 (MQ=255) aGTGTTTAACGACGTTAAGATGGCAATGTGCAGATCATGTTATCAGGGTGGGCCGGTAGCGTATATTc < 2:590659/68‑1 (MQ=255) aGTGTTTAACGACGTTAAGATGGCAATGTGCAGATCATGTTATCAGGGTGGGCCGGCACCGTATATTc < 2:649403/68‑1 (MQ=255) aGTGTTTAACGACGTTAAGATGGCAATGTGCAGACCATGTTATCAGGGTGGGCCGGTAGCGTATATTc < 2:903420/68‑1 (MQ=255) aGTGGTTAACGACGTTAAGATGGCAACGTGCAGATCATGTTATCAGGGTGGGCCGGTAGCGTATATTc < 2:526848/68‑1 (MQ=255) tgtTTAACGACGTTAAGATGGCAATGTGCAGATCATGTTATCAGGGTGGGGCGGTAGCGTATATTCCt > 2:269749/1‑68 (MQ=255) tttAACGACGTTAAGATGGCAATGTGCAGATCATGTTATCAGGGTGGGCCGGAAGCGTATATTCCTta < 1:426248/68‑1 (MQ=255) ttAACGACGTTAAGATGGCAATGTGCAGATCATGTTATCAGGGTGGGGCGGTAGCGTATATTCCTtat < 2:724182/68‑1 (MQ=255) cgTTAAGATGGCAATGTGCAGATCATGTTATCAGGGTGGGCCGGTAGCGTATATTCCTTATACAAAcc < 1:206595/68‑1 (MQ=255) cgTTAAGATGGCAATGTGCAGATCATGTTATCAGGGTGGGCCGGTAGCGTACATTCCTTATACAAAcc < 2:548771/68‑1 (MQ=255) tgCAGATCATGTTATCAGGGTGGGGCGGTAGCGTATATTCCTTATACAAACCATCATTCATTTCAATg > 1:95958/1‑68 (MQ=255) gCAGATCATGTTATCAGGGTGGGGCGGTAGCGTATATTCCTTATACACACCATCATACAATTCAATGa > 2:456564/1‑68 (MQ=255) ttATCAGGGTGGGCCGGTAGCGTATATTCCTTATACAAACCATCATTCATTTCAATGATTGTTATTGa < 2:574563/68‑1 (MQ=255) aTCAGGGTGGGGCGGTAGCGTATATTCCTTATACAAACCATCATTCATTTCAATGATTGTTATTGa > 2:298920/1‑66 (MQ=255) aTCAGGGTGGGGCGGTAGCGTATATTCCTTATACAAACCATCATTCATTTCAATGATTGTTATTGa < 1:298920/66‑1 (MQ=255) aTCAGGGTGGGGCGGTAGCGTATATTCCTTATACAAACCATCATTCATTTCAATGATTGTTATTGAat > 2:114898/1‑68 (MQ=255) aTCAGGGTGGGGCGGTAGCGTATATTCCTTATACAAACCATCATTCATTTCAATGATTGATATTGAat > 2:865513/1‑68 (MQ=255) tCAGGGTGGGGCGGTAGCGTATATTCCTTATACAAACCATCATTCATTTCAATGATTGTTTTTGAatt > 2:903536/1‑68 (MQ=255) tCAGGGTGGGGCGGTAGCGTATATTCCTTATACAAACCATCATTCATTTCAATGATTGATATTGAatt > 1:592001/1‑68 (MQ=255) tCAGGGTGGGCCGGTAGCGTATATTCCTTATACAAACCATCATTCATTTCAATGATTGTTATTGAatt < 1:480152/68‑1 (MQ=255) aGGGTGGGGCGGTAGCGTATATTCCTTATACAAACCATCATTCATTTCAATGATTGTTATTGAattat < 1:376112/68‑1 (MQ=255) gggTGGGGCGGTAGCGTATATTCCTTATACAAACCATCATTCATTTCAATGATTGTTATTGAattatt > 1:887747/1‑68 (MQ=255) ggTGGGGCGGTAGCGTATATTCCTTATACAAACCATCATTCATTTCAATGATTGTTATTGAAttattt > 1:294209/1‑68 (MQ=255) gTGGGGCGGTAGCGTATATTCCTTATACAAACCATCATTCATTTCAATGATTGTTATTGa > 1:929056/1‑60 (MQ=255) gTGGGGCGGTAGCGTATATTCCTTATACAAACCATCATTCATTTCAATGATTGTTATTGa < 2:929053/60‑1 (MQ=255) gTGGGGCGGTAGCGTATATTCCTTATACAAACCATCATTCATTTCAATGATTGTTATTGAAttattta > 1:690592/1‑68 (MQ=255) gCGGTAGCGTATATTCCTTATACAAACCATCATTCATTTCAATGATTGTTATTGAATTATTTAttgtt > 2:272398/1‑68 (MQ=255) | CCACTATTAAGCTAGTGTTTAACGACGTTAAGATGGCAATGTGCAGATCATGTTATCAGGGTGGGGCGGTAGCGTATATTCCTTATACAAACCATCATTCATTTCAATGATTGTTATTGAATTATTTATTGTT > NC_000913/1527537‑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: ‑ |