| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 3,948,799 | 1 | . | G | 29.7% | 103.3 / 104.3 | 64 | coding (1171/1521 nt) | yifB | putative magnesium chelatase YifB |
| Reads supporting (aligned to +/- strand): ref base . (22/23); new base G (10/9); total (32/32) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.14e-02 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
TAACGGTGGCGCTGCTTTCTCCCGGCACTACCGTTTTACTCAAAATGC‑C‑‑GGGGGGTGGTAATGGGATCTCCAGTGAGAGATCGAAGCGGTCGAGAAAG > NC_000913/3948751‑3948848 | tAACGGTGGCGCTGCTTTCTCCCGGCACTACCGTTTTACTCAAAATGc‑c < 1:1173130/49‑1 (MQ=255) aaCGGTGGCGCTGCTTTCTCCCGGCACTACCGTTTTACTCAAAATGC‑C‑‑gg > 1:5484532/1‑50 (MQ=255) aaCGGTGGCGCTGCTTTCTCCCGGCACTACCGTTTTACTCAAAATGC‑C‑‑g < 1:5293812/49‑1 (MQ=255) aCGGTGGCGCTGCTTTCTCCCGGCACTACCGTTTTACTCAAAATGC‑C‑‑gg < 1:1618828/49‑1 (MQ=255) ggTGGCGCTGCTTTCTCCCGGCACTACCGTTTTACTCAAAATGC‑C‑‑ggggg > 1:4192837/1‑50 (MQ=255) ggTGGCGCTGCTTTCTCCCGGCACTACCGTTTTACTCAAAATGC‑C‑‑ggggg > 1:4781527/1‑50 (MQ=255) tGGCGCTGCTTTCTCCCGGCACTACCGTTTTACTCAAAATGC‑C‑‑ggggggg > 1:2281633/1‑49 (MQ=255) tGGCGCTGCTTTCTCCCGGCACTACCGTTTTACTCAAAATGC‑C‑‑GGGGggt < 1:3321915/50‑1 (MQ=255) tGGCGCTGCTTTCTCCCGGCACTACCGTTTTACTCAAAATGC‑C‑‑GGGGggt < 1:3272285/50‑1 (MQ=255) gcgcTGCTTTCTCCCGGCACTACCGTTTTACTCAAAATGC‑C‑‑gggggggt < 1:1771788/49‑3 (MQ=255) gcgcTGCTTTCTCCCGGCACTACCGTTTTACTCAAAATGC‑C‑‑GGGGggtgg < 1:455087/50‑1 (MQ=255) gcgcTGCTTTCTCCCGGCACTACCGTTTTACTCAAAATGC‑C‑‑GGGGggtgg < 1:2576092/50‑1 (MQ=255) gcgcTGCTTTCTCCCGGCACTACCGTTTTACTCAAAATGC‑C‑‑GGGGggtgg > 1:4921976/1‑50 (MQ=255) cgcTGCTTTCTCCCGGCACTACCGTTTTACTCAAAATGC‑C‑‑GGGGggtggt < 1:760252/50‑1 (MQ=255) cgcTGCTTTCTCCCGGCACTACCGTTTTACTCAAAATGC‑C‑‑GGGGggtgg > 1:3017749/1‑49 (MQ=255) tgctTTCTCCCGGCACTACCGTTTTACTCAAAATGC‑C‑‑GGGGGGTGGTAAt > 1:5138596/1‑50 (MQ=255) gctTTCTCCCGGCACTACCGTTTTACTCAAAATGC‑C‑‑GGGGGGTGGTAATg > 1:4214607/1‑50 (MQ=255) ctTTCTCCCGGCACTACCGTTTTACTCAAAATGC‑CG‑GGGGGGTGGTAAt < 1:3143027/49‑1 (MQ=255) ctTTCTCCCGGCACTACCGTTTTACTCAAAATGC‑CG‑GGGGGGTGGTAATg > 1:4325484/1‑50 (MQ=255) tttCTCCCGGCACTACCGTTTTACTCAAAATGC‑C‑‑GGGGGGTGGTAATgg > 1:201443/1‑49 (MQ=255) tttCTCCCGGCACTACCGTTTTACTCAAAATGC‑C‑‑GGGGGGTGGTAATgg > 1:4214446/1‑49 (MQ=255) tctcCCGGCACTACCGTTTTACTCAAAATGC‑C‑‑GGGGGGTGGTAATGGGa < 1:2655355/49‑1 (MQ=255) tcCCGGCACTACCGTTTTACTCAAAATGC‑CG‑GGGGGGTGGTAATGGGAtc > 1:1372659/1‑50 (MQ=255) cccGGCACTACCGTTTTACTCAAAATGC‑C‑‑GGGGGGTGGTAATGGGAtctc > 1:198709/1‑50 (MQ=255) cccGGCACTACCGTTTTACTCAAAATGC‑C‑‑GGGGGGTGGTAATGGGAtctc > 1:2790321/1‑50 (MQ=255) ccGGCACTACCGTTTTACTCAAAATGC‑CG‑GGGGGGTGGTAATGGGAtctc > 1:4363201/1‑50 (MQ=255) ccGGCACTACCGTTTTACTCAAAATGC‑CG‑GGGGGGTGGTAATGGGAtctc > 1:3450916/1‑50 (MQ=255) ccGGCACTACCGTTTTACTCAAAATGC‑CG‑GGGGGGTGGTAATGGGAt < 1:2638923/47‑1 (MQ=255) ccGGCACTACCGTTTTACTCAAAATGC‑CG‑GGGGGGTGGTAATGGGAt < 1:1571810/47‑1 (MQ=255) ccGGCACTACCGTTTTACTCAAAATGC‑C‑‑GGGGGGTGGTAATGGGAtct < 1:4384599/48‑1 (MQ=255) ccGGCACTACCGTTTTACTCAAAATGC‑C‑‑GGGGGGTGGTAATGGGAtc < 1:3557615/47‑1 (MQ=255) cGGCACTACCGTTTTACTCAAAATGC‑C‑‑GGGGGGTGGTAATGGGAtctc > 1:5096653/1‑48 (MQ=255) ggCACTACCGTTTTACTCAAAATGC‑CG‑GGGGGGTGGTAATGGGAtctc > 1:3938791/1‑48 (MQ=255) ggCACTACCGTTTTACTCAAAATGC‑C‑‑GGGGGGTGGTAATGGGATCTCCAg > 1:5458368/1‑50 (MQ=255) cACTACCGTTTTACTCAAAATGC‑CG‑GGGGGGTGGTAATGGGATCTCCAGt < 1:4826253/50‑1 (MQ=255) cACTACCGTTTTACTCAAAATGC‑C‑‑GGGGGGTGGTAATGGGATCTCCAGTg < 1:4284866/50‑1 (MQ=255) aCTACCGTTTTACTCAAAATGC‑CG‑GGGGGGTGGTAATGGGATCTCCAGTg < 1:4528842/50‑1 (MQ=255) aCTACCGTTTTACTCAAAATGC‑C‑‑GGGGGGTGGTAATGGGATCTCCAGt > 1:1583157/1‑48 (MQ=255) cTACCGTTTTACTCAAAATGC‑C‑‑GGGGGGTGGTAATGGGATCTCCAGTgag < 1:4579051/50‑1 (MQ=255) tACCGTTTTACTCAAAATGC‑CG‑GGGGGGTGGTAATGGGATCTCCAGt > 1:5161018/1‑47 (MQ=255) tACCGTTTTACTCAAAATGC‑CG‑GGGGGGTGGTAATGGGATCTCCAGt > 1:1846851/1‑47 (MQ=255) aCCGTTTTACTCAAAATGC‑C‑‑GGGGGGTGGTAATGGGATCTCCAGTga > 1:1072126/1‑47 (MQ=255) aCCGTTTTACTCAAAATGC‑C‑‑GGGGGGTGGTAATGGGATCTCCAGTga > 1:3186488/1‑47 (MQ=255) ccGTTTTACTCAAAATGCAG‑‑GGGGGGTGGTAATGGGATCTCCAGTgagag < 1:5116531/50‑1 (MQ=255) ccGTTTTACTCAAAATGC‑CG‑GGGGGGTGGTAATGGGATCTCCAGTgag < 1:3476185/48‑1 (MQ=255) ccGTTTTACTCAAAATGC‑C‑‑GGGGGGTGGTAATGGGATCTCCAGTgagaga < 1:4589839/50‑1 (MQ=255) ccGTTTTACTCAAAATGC‑C‑‑GGGGGGTGGTAATGGGATCTCCAGTgagag < 1:3345378/49‑1 (MQ=255) ccGTTTTACTCAAAATGC‑C‑‑GGGGGGTGGTAATGGGATCTCCAGTga > 1:1697889/1‑46 (MQ=255) cGTTTTACTCAAAATGC‑CG‑GGGGGGTGGTAATGGGATCTCCAGTga > 1:3552772/1‑46 (MQ=255) cGTTTTACTCAAAATGC‑C‑‑GGGGGGTGGTAATGGGATCTCCAGTGAGAGAt < 1:933954/50‑1 (MQ=255) cGTTTTACTCAAAATGC‑C‑‑GGGGGGTGGTAATGGGATCTCCAGTGAGAGAt < 1:4843247/50‑1 (MQ=255) gTTTTACTCAAAATGC‑CGGGGGGGGTGGTAATGGGATCTCCAGTga > 1:3380288/1‑46 (MQ=255) gTTTTACTCAAAATGC‑C‑‑GGGGGGTGGTAATGGGATCTCCAGTga > 1:512295/1‑44 (MQ=255) gTTTTACTCAAAATGC‑C‑‑GGGGGGTGGTAATGGGATCTCCAGTGAGAGATc < 1:3281283/50‑1 (MQ=255) ttACTCAAAATGC‑C‑‑GGGGGGTGGTAATGGGATCTCCAGTGAGAGATCg > 1:84840/1‑48 (MQ=255) aCTCAAAATGC‑CG‑GGGGGGTGGTAATGGGATCTCCAGTGAGAGATCGaa < 1:840970/49‑1 (MQ=255) aCTCAAAATGC‑CG‑GGGGGGTGGTAATGGGATCTCCAGTGAGAGATCGAAg > 1:4072121/1‑50 (MQ=255) cTCAAAATGC‑CG‑GGGGGGTGGTAATGGGATCTCCAGTGAGAGATCGAAg < 1:319433/49‑1 (MQ=255) cTCAAAATGC‑CG‑GGGGGGTGGTAATGGGATCTCCAGTGAGAGATCGAAGc < 1:230097/50‑1 (MQ=255) cTCAAAATGC‑C‑‑GGGGGGTGGTAATGGGATCTCCAGTGAGAGATCGAAGc < 1:859515/49‑1 (MQ=255) cTCAAAATGC‑C‑‑GGGGGGTGGTAATGGGATCTCCAGTGAGAGATCGAAGc < 1:1887395/49‑1 (MQ=255) cTCAAAATGC‑C‑‑GGGGGGTGGTAATGGGATCTCCAGTGAGAGATCGAAGCg > 1:3386682/1‑50 (MQ=255) tCAAAATGC‑C‑‑GGGGGGTGGTAATGGGATCTCCAGTGAGAGATCGAAGCgg > 1:4252550/1‑50 (MQ=255) cAAAATGC‑C‑‑GGGGGGTGGTAATGGGATCTCCAGTGAGAGATCGAAGCGGt > 1:5103257/1‑50 (MQ=255) aaaaTGC‑C‑‑GGGGGGTGGTAATGGGATCTCCAGTGAGAGATCGAAGCGGt < 1:5147756/49‑1 (MQ=255) aaaaTGC‑C‑‑GGGGGGTGGTAATGGGATCTCCAGTGAGAGATCGAAGCGGTc > 1:794856/1‑50 (MQ=255) aaaTGC‑C‑‑GGGGGGTGGTAATGGGATCTCCAGTGAGAGATCGAAGCGGTc > 1:2057428/1‑49 (MQ=255) aaaTGC‑C‑‑GGGGGGTGGTAATGGGATCTCCAGTGAGAGATCGAAGCGGTCg < 1:4691403/50‑1 (MQ=255) aaaTGC‑C‑‑GGGGGGTGGTAATGGGATCTCCAGTGAGAGATCGAAGCGGTCg < 1:3442802/50‑1 (MQ=255) aTGC‑C‑‑GGGGGGTGGTAATGGGATCTCCAGTGAGAGATCGAAGCGGTCgag < 1:5535001/50‑1 (MQ=255) aTGC‑C‑‑GGGGGGTGGTAATGGGATCTCCAGTGAGAGATCGAAGCGGTCgag < 1:5475544/50‑1 (MQ=255) gC‑C‑‑GGGGGGTGGTAATGGGATCTCCAGTGAGAGATCGAAGCGGTCgaga > 1:1617406/1‑49 (MQ=255) c‑‑GGGGGGTGGTAATGGGATCTCCAGTGAGAGATCGAAGCGGTCGAGAAAg < 1:2016469/50‑1 (MQ=255) c‑‑GGGGGGTGGTAATGGGATCTCCAGTGAGAGATCGAAGCGGTCGAGAAAg < 1:250358/50‑1 (MQ=255) | TAACGGTGGCGCTGCTTTCTCCCGGCACTACCGTTTTACTCAAAATGC‑C‑‑GGGGGGTGGTAATGGGATCTCCAGTGAGAGATCGAAGCGGTCGAGAAAG > NC_000913/3948751‑3948848 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |