Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 4,414,906 | 0 | T | G | 40.9% | 13.0 / 23.7 | 22 | V211G (GTG→GGG) | aidB | putative acyl‑CoA dehydrogenase AidB |
Reads supporting (aligned to +/- strand): ref base T (0/13); new base G (9/0); total (9/13) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 2.01e-06 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
TCCGATGTTATGAGCAACACCACCCGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCTGGTGGGGCATAAA‑TGGTTTTTCTC‑GGTTCCGCAAAGCGATGCGCATCTGGTGCTGGCGCAGA > NC_000913/4414845‑4414965 | tCCGATGTTATGAGCAACACCACCCGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCTgggggg > 1:3757982/1‑65 (MQ=255) gTTATGAGCAACACCACCCGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCTGGTGGGGCATaa < 1:1056703/65‑1 (MQ=255) tGAGCAACACCACCCGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCTGGTGGGGCATAAA‑Tgg < 1:703232/65‑1 (MQ=255) gAGCAGCGCCACCCGTGCAGAGCGTGTGGAAGATGGCTGTTATCGGCTGGTGGGGCATAAA‑TGGt < 1:648349/65‑1 (MQ=255) gAGCAACGCCACCCGTGCAGATCGTCTGGAAGATGGCTCTTATCGGCTGGTGGGGCATAAA‑TGGt < 1:3692204/65‑1 (MQ=255) gCAACACCACCCGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCTGGGGGGGGGTAAA‑AGGgtt > 1:2333291/1‑65 (MQ=255) accaccCGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCTGGGGGGGGGGAAATTGGttttttt > 1:2913801/1‑63 (MQ=255) accCGTNCAGAGCGTCTGGAAGATGGCTCTTATCGGCTGGGGGGGGATAAA‑AGGGTTTTTTT‑GGt > 1:2211062/1‑65 (MQ=255) accCGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCTGGTGGGGCGTAAA‑TGGTTTTTCTC‑GGt < 1:2796312/65‑1 (MQ=255) ccGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCTGGGGGGGGGGAAA‑TGGTTTTTTTC‑GGTTc > 1:2663617/1‑65 (MQ=255) agCGTCTGGAAGATGGCTCTTATCGGCTGGGGGGGGGTAAA‑AAGGTTTTTTT‑GGTTTCGCCAAGc > 1:2183551/1‑65 (MQ=255) agCGTCTGGAAGATGGCTCTGATCGGCTGGTGGGGCATAAA‑TGGTTTTTCTC‑GGTTCCGCAAAGc < 1:1584331/65‑1 (MQ=255) gCGTCTGGAAGACGGCTCTTATCGGCTGGGGGGGGGTAAA‑TGGGTTTTTTCTGGTCCCCCaaacg > 1:2292683/1‑63 (MQ=255) cGTCTGGAAGATGGCTCTTATCGGCTGGGGGGGGGTAAA‑AGGGTTTTTTC‑GGGTCCGCCAAACGa > 1:3667294/1‑65 (MQ=255) gTCTGGAAGATGGCTCTTATCGGCTGGTGGGGCATAAA‑TGGTTTTTCTC‑GGTTCCGCAAAGCGAt < 1:1624090/65‑1 (MQ=255) tGGAAGATGGCTCTTATCGGCTGGTGGGGCATAAA‑TGGTTTTTCTC‑GGTTCCGCAAAGCGATgcg < 1:793562/65‑1 (MQ=255) tGGAAGATGGCTCTTATCGGCTGGTGGGGCATAAA‑TGGTTTTTCTC‑GGTTCCGCAAAGCGATgcg < 1:1041603/65‑1 (MQ=255) tGGAAGATGGCTCTTATCGGCTGGTGGGGCATAAA‑TGGTTTTTCTC‑GGTTCCGCAAAGCGATgcg < 1:1592419/65‑1 (MQ=255) gAAGATGGCTCTTATCGGCTGGTGGGGCATAAA‑TGGTTTTTCTC‑GGTTCCGCAAAGCGATGCGCa < 1:334503/65‑1 (MQ=255) aGATGGCTCTTATCGGCTGGTGGGGCATAAA‑TGGTTTTTCTC‑GGTTCCGCAAAGCGATGCGCATc < 1:2684352/65‑1 (MQ=255) aGATGGCTCTTATCGGCTGGTGGGGCATAAA‑TGGTTTTTCTC‑GGTTCCGCAAAGCGATGCGCATc < 1:221479/65‑1 (MQ=255) aTGGCTCTTATCGGCTGGGGGGGGGTAAA‑TGGGTTTTTTT‑GGTTCCGCAAAACGATGCGCATCTg > 1:1090132/1‑65 (MQ=255) aTGGCTCTTATCGGCTGGGGGGGGATAAA‑TGGGTTTTTTC‑CGTTCCGCCAAGCGATGCGCCTCTg > 1:2400168/1‑65 (MQ=255) tctTATCGGCTGGGGGGGGGTAAA‑AGGTTTTTCTC‑CGTTCCCCAAAACGAAGCGCATCTGGTGcc > 1:3573477/1‑64 (MQ=255) ttATCGGCTGGTGGGGCATAAA‑TGGTTTTTCTC‑GGTTCCGCAAAGCGATGCGCATCTGGTGCTgg < 1:3280218/65‑1 (MQ=255) gCTGGTGGGGCATAAA‑TGGTTTTTCTC‑GGTTCCGCAAAGCGATGCGCATCTGGTGCTGGCGCAGa < 1:2980497/65‑1 (MQ=255) gCTGGTGGGGCATAAA‑TGGTTTTTCTC‑GGTTCCGCAAAGCGATGCGCATCTGGTGCTGGCGCAGa < 1:3556237/65‑1 (MQ=255) gCTAGTGGGGCATAAA‑TGGTTTTTCTG‑GGTTCCGCAAAGCGAGGCGCATCTGGTGCTGGCGCAGa < 1:2524999/65‑1 (MQ=255) | TCCGATGTTATGAGCAACACCACCCGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCTGGTGGGGCATAAA‑TGGTTTTTCTC‑GGTTCCGCAAAGCGATGCGCATCTGGTGCTGGCGCAGA > NC_000913/4414845‑4414965 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 20 ≤ ATCG/ATCG < 28 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |