Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 4,414,906 | 0 | T | G | 39.3% | 27.3 / 24.5 | 28 | V211G (GTG→GGG) | aidB | putative acyl‑CoA dehydrogenase AidB |
Reads supporting (aligned to +/- strand): ref base T (0/17); new base G (11/0); total (11/17) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 4.66e-08 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
GGTTCCGATGTTATGAGCAACACCACCCGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCTGGTGGGGCATAAATGGTTTTTCT‑CGGTTCCGCAAAGCGATGCGCATCTGGTGCTGGCGCAGA > NC_000913/4414842‑4414965 | ggTTCCGATGTTATGAGCAACACCACCCGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCtggt < 1:4314695/65‑1 (MQ=255) tATGAGCAACACCACCCGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCTGGGGGGGGGTaaaa > 1:3367901/1‑64 (MQ=255) tATGAGCAACACCACCCGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCTGGGGGGGGATAAAt > 1:115449/1‑65 (MQ=255) tGAGCAACACCACCCGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCTGGGGGGGGATAACTgg > 1:335019/1‑65 (MQ=255) gAGCAACACCACCCGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCTGGGGGGGGGTAAAAggg > 1:2775777/1‑64 (MQ=255) gCAACACCACCCGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCTGGGGGGGGGTAAAAGGgtt > 1:2233390/1‑65 (MQ=255) gCAACACCACCCGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCTGGGGGGGGATAAATGGgtt > 1:1910249/1‑65 (MQ=255) cgacacCACCCGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCTGGTGGGGCATAAATGGtttt < 1:4271119/63‑1 (MQ=255) acacCGCCCGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCTGGTGGGGCATAAATGGTTTTtc < 1:1531374/65‑1 (MQ=255) acacCGCCCGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCTGGTGGGGCATAAATGGTTTTtc < 1:4097663/65‑1 (MQ=255) acacCACCCGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCTGGTGGGGCATAAATGGTTTTtc < 1:431162/65‑1 (MQ=255) accCGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCTGGTGGGGCATAAATGGTTTTTCT‑CGGt < 1:2500718/65‑1 (MQ=255) accCGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCTGGTGGGGCATAAATGGTTTTTCT‑CGGt < 1:194693/65‑1 (MQ=255) accCGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCTGGGGGGGGGTAAATGGTTTTTTT‑TGGt > 1:3580222/1‑65 (MQ=255) accCGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCTGGGGGGGGGGAAAAGGGTTTTTT‑Cggc > 1:4047392/1‑64 (MQ=255) ccGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCTGGGGGGGGGTAAATTGTTTTTTT‑CGGtgc > 1:358943/1‑63 (MQ=255) cGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCTGGGGGGGGGTAAAAGGTTTTTTTCCGGTTc > 1:2630923/1‑65 (MQ=255) tGCAGAGCGTCTGGAAGATGGCTCTTATCGGCTGGTGGGGCATAAATGGTTTTTCT‑CGGTTCCGc < 1:1822537/65‑1 (MQ=255) agagCGTCTGGAAGATGGCTCTTATCGGCTGGGGGGGGGTAAATGGTTTTTTT‑TGGTTTCGCaca > 1:3115799/1‑63 (MQ=255) agagCGTCTGGAAGATGGCTCTTATCGGCTGGGGGGGGGTAAATGGGTTTTCT‑CGGTTCCTCaca > 1:4335036/1‑63 (MQ=255) agCGTCTGGAAGATGGCTCTTATCGGCTGGTGGGGCATAAATGGTTTTTCT‑CGGTTCCGCAAAGc < 1:3041755/65‑1 (MQ=255) agCGTCTGGAAGATGGCTCTTATCGGCTGGTGGGGCATAAATGGTTTTTCT‑CGGTTCCGCAAAGc < 1:4088460/65‑1 (MQ=255) agCGTCTGGAAGATGGCTCTTATCGGCTGGTGGGGCATAAATGGTTTTTCT‑CGGTTCCGCAAAGc < 1:1159878/65‑1 (MQ=255) agCGTCTGGAAGATGGCTCTTATCGGCTGGGGGGGGGTAAATGGTTTTTTT‑CTGTTCCGCACAGc > 1:2608398/1‑65 (MQ=255) gCGTCTGGAAGATGGCTCTTATCGGCTGGGGGGGGATAAATGGGTTTTTT‑TGGGTCCGCAAAACg > 1:680913/1‑65 (MQ=255) tGGAAGATGGCTCTTATCGGCTGGTGGGGCATAAATGGTTTTTCT‑CGGTTCCGCAAAGCGATgcg < 1:3968296/65‑1 (MQ=255) tGGAAGATGGCTCTTATCGGCTGGTGGGGCATAAATGGTTTTTCT‑CGGTTCCGCAAAGCGATgcg < 1:608307/65‑1 (MQ=255) gAAGATGGCTCTTATCGGCTGGTNGGGCATAAATGGTTTTTCT‑CGGTTCCGCAAAGCGATGCGCa < 1:3410941/65‑1 (MQ=255) gAAGATGGCTCTTATCGGCTGGTGGGGCATAAATGGTTTTTCT‑CGGTTCCGCAAAGCGATGCGCa < 1:1900926/65‑1 (MQ=255) gAAGATGGCTCTTATCGGCTGGTGGGGCATAAATGGTTTTTCT‑CGGTTCCGCAAAGCGATGCGCa < 1:1799101/65‑1 (MQ=255) gAAGATGGCTCTTATCGGCTGGTGGGGCATAAATGGTTTTTCT‑CGGTTCCGCAAAGCGATGCCCa < 1:2176829/65‑1 (MQ=255) gAAGATGGCTCTTATCCGCTGGTGGGGCATAAATGGTTTTTCT‑CGGTTCCGCAAAGCGATGCGCa < 1:1797648/65‑1 (MQ=255) gAAGATGGCTCTTACCGGCTGGTGGGGCATAAATGGTTTTTCT‑CGGTTCCGCAAAGCGATGCGCa < 1:2883938/65‑1 (MQ=255) aaGATGGCTCTTATCGGCTGGGGGGGGGTTAATAGTTTTTTT‑TGGTTTCGCCAAACGATgcgcct > 1:3714496/1‑63 (MQ=255) aGATGGCTCTTATCGGCTGGTGGGGCATAAATGGTTTTTCT‑CGGTTCCGCAAAGCGATGCGCATc < 1:2895639/65‑1 (MQ=255) aGATGGCTCTTATCGGCTGGGGGGGGATAAACAGTTTTTTT‑CCGTTCCGTACAGCGCTGCGCATc > 1:3308930/1‑65 (MQ=255) ggCTCTTATCGGCTGGTGGGGCATAAATGGTTTTTCT‑CGGTTCCGCAAAGCGATGCGCATCTGGt < 1:3887560/65‑1 (MQ=255) tctTATCGGCTGGTGGGGCATAAATGGTTTTGCT‑CGGTTCCGCAAAGCGATGCGCATCTGGTGCt < 1:2544018/65‑1 (MQ=255) gCTGGTGGGGCATAAATGGTTTTTCT‑CGGTTCCGCAAAGCGATGCGCATCTGGTGCTGGCGCAGa < 1:3058369/65‑1 (MQ=255) gCTGGTGGGGCATAAATGGTTTTTCT‑CGGTTCCGCAAAGCGATGCGCATCTGGTGCTGGCGCAGa < 1:3824702/65‑1 (MQ=255) gCTGGTGGGGCATAAATGGTTTTTCT‑CGGTTCCGCAAAGCGATGCGCATCTGGTGCTGGCGCAGa < 1:2634210/65‑1 (MQ=255) | GGTTCCGATGTTATGAGCAACACCACCCGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCTGGTGGGGCATAAATGGTTTTTCT‑CGGTTCCGCAAAGCGATGCGCATCTGGTGCTGGCGCAGA > NC_000913/4414842‑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: ‑ |