Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_017635 | 4,327,011 | 0 | G | T | 30.8% | 87.4 / 14.6 | 39 | intergenic (+84/‑176) | sodA/kdgT | superoxide dismutase [Mn]/2‑keto‑3‑deoxygluconate transporter |
Reads supporting (aligned to +/- strand): ref base G (14/13); new base T (12/0); total (26/13) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 2.97e-03 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
CGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGTTTTTTTGTGAT‑CAATTTCAAAATAAAAACAATGATCCGATAAAAATAAAACAGTG > NC_017635/4326962‑4327059 | cGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGGTTTTTtg > 1:1179853/1‑50 (MQ=255) gCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGTTTTTTtgt < 1:967674/50‑1 (MQ=255) 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 > 2:328195/1‑49 (MQ=25) gAACCATCGCCAACGGCGGTGGGTTTTTTTTGAT‑CAATTTTaaaaaaaa > 1:1802582/1‑49 (MQ=25) aaCCATCGCCAACGGCGGTGGTTTTTTTGTGAT‑CAATTTCAAAATaaaa < 1:1804087/49‑1 (MQ=255) aaCCATCGCCAACGGCGGTGGTTTTTTTGTGAT‑CAATTTCAAAATaaaa < 2:1086998/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) ttttttGTGAT‑CAATTTCAAAATAAAAACAATGATCCGataaaaataaaa > 1:380563/1‑50 (MQ=255) ttttttGTGAT‑CAATTTCAAAATAAAAACAATGATCCGataaaaataaaa > 2:2039745/1‑50 (MQ=255) tttttGTGAT‑CAATTTCAAAATAAAAACAATGATCCGataaaaataaaa > 1:590926/1‑49 (MQ=255) ttttGTGAT‑CAATTTCAAAATAAAAACAATGATCCGataaaaataaaa > 1:1783569/1‑48 (MQ=255) tttGTGAT‑CAATTTCAAAATAAAAACAATGATCCGATAAAAATAAAACAg > 2:1374054/1‑50 (MQ=255) ttGTGAT‑CAATTTCAAAATAAAAACAATGATCCGATAAAAATAAAACAgt > 1:500678/1‑50 (MQ=255) ttGTGAT‑CAATTTCAAAATAAAAACAATGATCCGATAAAAATAAAACAg < 1:1472209/49‑1 (MQ=255) tgtgAT‑CAATTTCAAAATAAAAACAATGATCCGATAAAAATAAAACAgt < 1:249869/49‑1 (MQ=255) gtgAT‑CAATTTCAAAATAAAAACAATGATCCGATAAAAATAAAACAgtg < 1:1491347/49‑1 (MQ=255) gtgAT‑CAATTTCAAAATAAAAACAATGATCCGATAAAAATAAAACAgtg < 1:762934/49‑1 (MQ=255) | CGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGTTTTTTTGTGAT‑CAATTTCAAAATAAAAACAATGATCCGATAAAAATAAAACAGTG > NC_017635/4326962‑4327059 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |