Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_017635 | 4,327,004 | 0 | T | G | 51.6% | 35.1 / 47.5 | 31 | intergenic (+77/‑183) | sodA/kdgT | superoxide dismutase [Mn]/2‑keto‑3‑deoxygluconate transporter |
Reads supporting (aligned to +/- strand): ref base T (3/12); new base G (16/0); total (19/12) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 3.22e-06 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.89e-07 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
TATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGTTTTTTTGTGAT‑CAATTTCAAAATAAAAACAATGATCCGATAAAAATAA > NC_017635/4326955‑4327052 | tataGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCggtggt < 2:836400/50‑1 (MQ=255) taGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGtt < 1:1737056/49‑1 (MQ=255) aGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGgttt > 1:1527144/1‑50 (MQ=255) ggCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGtttt < 2:1496109/49‑1 (MQ=255) gCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGtttttt > 1:661651/1‑50 (MQ=255) ccGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGgtttttt > 1:2049103/1‑50 (MQ=255) ccGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGGGGgtttttt > 2:814491/1‑50 (MQ=38) cGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGGTTTTTtg > 1:1179853/1‑50 (MQ=255) gCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGgtttttttt > 2:267736/1‑48 (MQ=39) gCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGgtttttttg > 2:1516999/1‑48 (MQ=39) gCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGTTTTTTtgt < 1:967674/50‑1 (MQ=255) cATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGgtttttttt > 1:971146/1‑47 (MQ=38) cATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGTTTTTTtttg > 1:239887/1‑50 (MQ=255) atatCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGGTTTTTTTTGa > 2:55981/1‑50 (MQ=38) tatCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGTTTTTTTGTGAt < 1:762702/50‑1 (MQ=255) tatCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGGTTTTTTTTGAt > 1:594787/1‑50 (MQ=38) tatCAGCTTAAAAAATGAACCATCGCCAACGGCGGGGGGTTTTTtttggt > 2:1195771/1‑48 (MQ=37) atCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGGTTTTTtttggt > 1:1090833/1‑47 (MQ=38) atCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGGTTTTTTTTGAt > 1:854366/1‑49 (MQ=38) tCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGTTTTTTTGTGAT‑Ca < 2:1849868/50‑1 (MQ=255) cAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGTTTTTTTTTGAT‑Ca > 1:1583695/1‑49 (MQ=255) gCTTAAAAAATGAACCATCGCCAACGGCGGTGGTTTTTTTGTGAT‑CAAtt < 2:2032986/50‑1 (MQ=255) gCTTAAAAAATGAACCATCGCCAACGGCGGTGGGTTTTTTTTGAT‑Caaat > 1:1045848/1‑48 (MQ=38) cTTAAAAAATGAACCATCGCCAACGGCGGTGGGTTTTTTTTGAT‑CAAttt > 1:593377/1‑50 (MQ=38) cTTAAAAAATGAACCATCGCCAACGGCGGTGGGTTTTTTGTTGT‑AAAttt > 1:920349/1‑50 (MQ=37) cTTAAAAAATGAACCATCGCCAACGGCGGTGGGTTTTTTGTGAT‑CAAttt > 2:24671/1‑50 (MQ=255) ttAAAAAATGAACCATCGCCAACGGCGGAGGGTTTTTTTTGATCCAAttt > 2:419058/1‑50 (MQ=25) aaaaaaTGAACCATCGCCAACGGCGGTGGGTTTTTTTTGGT‑CCAtttta > 2:1772922/1‑47 (MQ=25) aaaaaTGAACCATCGCCAACGGCGGTGGGTTTTTTTTGAT‑CAATTTCaa > 1:1425646/1‑49 (MQ=38) aaaaTGAACCATCGCCAACGGCGGTGGTTTTTTTGTGAT‑CAATTTCaaaa < 2:311770/50‑1 (MQ=255) aaTGAACCATCGCCAACGGCGGTGGTTTTTTTGTGAT‑CAATTTCAAAATa < 2:897416/50‑1 (MQ=255) gAACCATCGCCAACGGCGGTGGGTTTTTTTTGAT‑CAATTTTaaaaaaaa > 1:1802582/1‑49 (MQ=25) gAACCATCGCCAACGGCGGTGGGTTTTTTTTGAT‑CAATTTTaaaaaaaa > 2:328195/1‑49 (MQ=25) aaCCATCGCCAACGGCGGTGGTTTTTTTGTGAT‑CAATTTCAAAATaaaa < 2:1086998/49‑1 (MQ=255) aaCCATCGCCAACGGCGGTGGTTTTTTTGTGAT‑CAATTTCAAAATaaaa < 1:1804087/49‑1 (MQ=255) aaCCATCGCCAACGGCGGTGGGTTTTTTTTGAT‑CCATTTtaaaataaaa > 2:1733598/1‑49 (MQ=37) ccATCGCCAACGGCGGTGGTTTTTTTGTGAT‑CAATTTCAAAATAAAAACa < 2:161262/50‑1 (MQ=255) gCCAACGGCGGTGGTTTTTTTGTGAT‑CAATTTCAAAATAAAAACAATGa < 2:1534389/49‑1 (MQ=255) ccAACGGCGGTGGTTTTTTTGTGAT‑CAATTTCAAAATAAAAACAATGATc < 1:615653/50‑1 (MQ=255) cAACGGCGGTGGTTTTTTTGTGAT‑CAATTTCAAAATAAAAACAATGATcc < 2:107979/50‑1 (MQ=255) ggcggTGGTTTTTTTGTGAT‑CAATTTCAAAATAAAAACAATGATCCGat > 1:1066202/1‑49 (MQ=255) ggtggtTTTTTTGTGAT‑CAATTTCAAAATAAAAACAATGATCCGataaaa < 2:894558/50‑1 (MQ=255) gtggtTTTTTTGTGAT‑CAATTTCAAAATAAAAACAATGATCCGataaaaa > 1:362211/1‑50 (MQ=255) tggtTTTTTTGTGAT‑CAATTTCAAAATAAAAACAATGATCCGataaaaa > 1:494964/1‑49 (MQ=255) gtTTTTTTGTGAT‑CAATTTCAAAATAAAAACAATGATCCGataaaaata > 2:1544149/1‑49 (MQ=255) gtTTTTTTGTGAT‑CAATTTCAAAATAAAAACAATGATCCGataaaaata > 2:640146/1‑49 (MQ=255) tttttttGTGAT‑CAATTTCAAAATAAAAACAATGATCCGataaaaataa > 2:1247490/1‑49 (MQ=255) | TATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGTTTTTTTGTGAT‑CAATTTCAAAATAAAAACAATGATCCGATAAAAATAA > NC_017635/4326955‑4327052 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |