Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 954,398 | 0 | A | C | 27.9% | 129.1 / 11.2 | 62 | G23G (GGT→GGG) | focA | formate channel FocA |
Reads supporting (aligned to +/- strand): ref base A (22/22); new base C (0/17); total (22/40) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.77e-04 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
TAATCGCCAGATAGAAAGTCTTAAGCGGATGTTTCGTTGCTTTATAGACACCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCA > NC_000913/954349‑954446 | tAATCGCCAGATAGAAAGTCTTAAGCGGATGTTTCGTTGCTTTATAGaca < 1:3926128/50‑1 (MQ=255) aaTCGCCAGATAGAAAGTCTTAAGCGGATGTTTCGTTGCTTTATAGacac > 1:7461891/1‑50 (MQ=255) aTCGCCAGATAGAAAGTCTTAAGCGGATGTTTCGTTGCTTTATAGACAcc > 1:8482549/1‑50 (MQ=255) aTCGCCAGATAGAAAGTCTTAAGCGGATGTTTCGTTGCTTTATAGACAcc < 1:7111254/50‑1 (MQ=255) gCCAGATAGAAAGTCTTAAGCGGATGTTTCTTTGCTTTATAGACACCCg < 1:7039197/49‑1 (MQ=255) gCCAGATAGAAAGTCTTAAGCGGATGTTTCGTTGCTTTATAGACACCCg < 1:4742540/49‑1 (MQ=255) gCCAGATAGAAAGTCTTAAGCGGATGTTTCGTTGCTTTATAGACACCCGc < 1:1166812/50‑1 (MQ=255) gCCAGATAGAAAGTCTTAAGCGGATGTTTCGTTGCTTTATAGACACCCGc < 1:4929559/50‑1 (MQ=255) gCCAGATAGAAAGTCTTAAGCGGATGTTTCGTTGCTTTATAGACACCCGc < 1:4738280/50‑1 (MQ=255) gCCAGATAGAAAGTCTTAAGCGGATGTTTCGTTGCTTTATAGACACCCGc > 1:262063/1‑50 (MQ=255) ccAGATAGAAAGTCTTAAGCGGATGTTTCGTTGCTTTATAGACACCCGcc > 1:2789294/1‑50 (MQ=255) ccAGATAGAAAGTCTTAAGCGGATGTTTCGTTGCTTTATAGACACCCGcc > 1:3717709/1‑50 (MQ=255) cAGATAGAAAGTCTTAAGCGGATGTTTCGTTGCTTTATAGACACCCGCct < 1:6009974/50‑1 (MQ=255) cAGATAGAAAGTCTTAAGCGGATGTTTCGTTGCTTTATAGACACCCGCct > 1:272307/1‑50 (MQ=255) aTAGAAAGTTTTAAGCGGTTGTTTCGTTGCTTTATAGACACCCGCCTCtt < 1:8150792/50‑1 (MQ=255) aTAGAAAGTCTTAAGCGGATGTTTCGTTGCTTTATAGACACCCGCCTCtt < 1:5876936/50‑1 (MQ=255) aTAGAAAGTCTTAAGCGGATGTTTCGTTGCTTTATAGACACCCGCCTCtt < 1:4453758/50‑1 (MQ=255) tAGAAAGTCTTAAGCGGATGTTTCGTTGCTTTATAGACACCCGCCTCTTc > 1:760188/1‑50 (MQ=255) aGAAAGTCTTAAGCGGATGTTTCGTTGCTTTATAGACACCCGCCTCTTCg > 1:3802595/1‑50 (MQ=255) gAAAGTCTTAAGCGGATGTTTCGTTGCTTTATAGACACCCGCCTCTTCgg > 1:34445/1‑50 (MQ=255) aaaGTCTTAAGCGGATGTTTCGTTGCTTTATAGACACCCGCCTCTTCGGc > 1:3342150/1‑50 (MQ=255) aaGTCTTAAGCGGTTGTTTCGTTGTTTTATAGACCCCCGCCTCTTCGGcc < 1:3153945/50‑1 (MQ=255) aGTCTTAAGCGGTTGTTTCGTTGTTTTATAGACCCCCGCCTCTTCGGCCa < 1:3855115/50‑1 (MQ=255) ttAAGCGGATGTTTCGTTGCTTTATAGACACCCGCCTCTTCGGCCACttt > 1:4409336/1‑50 (MQ=255) aaGCGGATGTTTCGTTGCTTTATAGACACCCGCCTCTTCGGCCACTTTgg > 1:5526607/1‑50 (MQ=255) gCGGATTTTTATTTGGTTTATAACCGCCCGCCTCTTCGGCCACTTTGGcc < 1:8308616/50‑1 (MQ=255) gCGGATGTTTCGTTGCTTTATAGACACCCGCCTCATCGGCCACTTTGGc > 1:5081269/1‑49 (MQ=255) gCGGATGTTTCGTTGCTTTATAAACCCCCGCCTCTTCGGCCACTTTGGcc < 1:2908770/50‑1 (MQ=255) cGGATTTTTCTTTGTTTTATAGACACCCGCCTCTTCGGCCACTTTGGCCa < 1:2479066/50‑1 (MQ=255) ggATTTTTGTTTGCTTTATAGACACCCGCCTCTTCGGCCACTTTGGCCAt < 1:2603974/50‑1 (MQ=255) ggATGTTTCGTTGCTTTATAGACACCCGCCTCTTCGGCCACTTTGGCCAt > 1:6452368/1‑50 (MQ=255) gttGTTTCTTTGTTTTAAAGACCCCCGCCTCTTCGGCCACTTTGGCCAtt < 1:7949586/48‑1 (MQ=255) gttGTTTCGTTGTTTTATAGACCCCCGCCTCTTCGGCCACTTTGGCCAtt < 1:2167168/48‑1 (MQ=255) gATGTTTCGTTGTTTAATAGACACCCGCCTCTTCGGCCACTTTGGCCAtt < 1:6651723/50‑1 (MQ=255) gATGTTTCGTTGCTTTATAGACACCCGCCTCTTCGGCCACTTTGGCCAtt > 1:7678275/1‑50 (MQ=255) tGTTTCGTTGCTTTATAGACACCCGCCTCTTCGGCCACTTTGGCCATTg > 1:5923226/1‑49 (MQ=255) gTTTCGTTGTTTTATAGACCCCCGCCTCTTCGGCCACTTTGGCCATTgc < 1:3210325/49‑1 (MQ=255) gTTTCGTTGCTTTATAGACACCCGCCTCTTCGGCCACTTTGGCCATTg > 1:512384/1‑48 (MQ=255) tttCTTTGTTTTATAGACCCCCGCCTCTTCGGCCACTTTGGCCATTgca < 1:2829151/49‑1 (MQ=255) tttCGTTGTTTTATAGCCCCCCGCCTCTTCGGCCACTTTGGCCATTgca < 1:8159872/49‑1 (MQ=255) tttCGTTGTTTTATAGACCCCCGCCTCTTCGGCCACTTTGGCCATTgca < 1:2392192/49‑1 (MQ=255) ttCGTTTTTTTATAGACCCCCGCCTCTTCGGCCACTTTGGCCATTgcag < 1:637411/49‑1 (MQ=255) ttCGTTGTTTTATAGACCCCCGCCTCTTCGGCCACTTTGGCCATTgcag < 1:4695212/49‑1 (MQ=255) ttCGTTGCTTTATAGACACCCGCCTCTTCGGCCACTTTGGCCATTgcag > 1:8166024/1‑49 (MQ=255) tCGTTGCTTTATAGACACCCGCCTCTTCGGCCACTTTGGCCATTgcag > 1:6681533/1‑48 (MQ=255) ggttgttTTAAAGACCCCCGCCTCTTCGGCCACTTTGGCCATTgcagca < 1:2321566/48‑1 (MQ=255) gTTGCTTTATAGACACCCGCCTCTTCGGCCACTTTGGCCATTGcagcag < 1:2814344/49‑1 (MQ=255) ttGCTTTATAGACACCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAgg > 1:7895013/1‑49 (MQ=255) gttttATAGCCACCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGaa < 1:1874881/47‑1 (MQ=255) gttttATAGACCCCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGaa < 1:3306366/47‑1 (MQ=255) gttttATAGACACCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGaa < 1:1584498/47‑1 (MQ=255) gttttATAAACACCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGaa < 1:4286359/47‑1 (MQ=255) ttttATAGACCCCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAg < 1:4090139/48‑1 (MQ=255) cTTTATAGACACCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAg < 1:1263136/49‑1 (MQ=255) cTTTATAGACACCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGt < 1:6246481/50‑1 (MQ=255) tttATAGACACCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGt < 1:4968513/49‑1 (MQ=255) tttAAAGACCCCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGt < 1:2751334/49‑1 (MQ=255) taataGCCCCCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTa < 1:6924464/47‑1 (MQ=255) tataGACCCCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTa < 1:1832042/48‑1 (MQ=255) taGACACCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTaa > 1:5186473/1‑47 (MQ=255) aaGAAACCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTaaaa < 1:7750254/48‑1 (MQ=255) gaacaaCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGa < 1:1025364/48‑1 (MQ=255) aGACACCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTaa > 1:4194047/1‑46 (MQ=255) aGACACCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTaa > 1:293459/1‑46 (MQ=255) gACACCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGAt < 1:6937855/50‑1 (MQ=255) acacCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGATc < 1:2047893/50‑1 (MQ=255) cacCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCa < 1:1677668/50‑1 (MQ=255) cacCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCa < 1:5482419/50‑1 (MQ=255) acCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCa > 1:6711994/1‑49 (MQ=255) | TAATCGCCAGATAGAAAGTCTTAAGCGGATGTTTCGTTGCTTTATAGACACCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCA > NC_000913/954349‑954446 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 18 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |