| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 3,948,799 | 1 | . | G | 34.0% | 56.9 / 93.7 | 50 | coding (1171/1521 nt) | yifB | putative magnesium chelatase YifB |
| Reads supporting (aligned to +/- strand): ref base . (23/10); new base G (7/10); total (30/20) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 7.05e-02 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.77e-01 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
AACGGTGGCGCTGCTTTCTCCCGGCACTACCGTTTTACTCAAAATGCC‑GGGGGGTGGTAATGGGATCTCCAGTGAGAGATCGAAGCGGTCGAGAA > NC_000913/3948752‑3948846 | aaCGGTGGCGCTGCTTTCTCCCGGCACTACCGTTTTACTCAAAATGCC‑g < 1:2827769/49‑1 (MQ=255) cGGTGGCGCTGCTTTCTCCCGGCACTACCGTTTTACTCAAAATGCC‑gggg < 1:2151987/50‑1 (MQ=255) cGGTGGCGCTGCTTTCTCCCGGCACTACCGTTTTACTCAAAATGCC‑gggg > 1:3959985/1‑50 (MQ=255) tGGCGCTGCTTTCTCCCGGCACTACCGTTTTACTCAAAATGCC‑ggggggg < 1:3855180/50‑2 (MQ=255) tGGCGCTGCTTTCTCCCGGCACTACCGTTTTACTCAAAATGCC‑GGGGggt > 1:3445827/1‑50 (MQ=255) gcgcTGCTTTCTCCCGGCACTACCGTTTTACTCAAAATGCC‑GGGGggtgg < 1:298644/50‑1 (MQ=255) gcgcTGCTTTCTCCCGGCACTACCGTTTTACTCAAAATGCC‑GGGGggtgg > 1:3538300/1‑50 (MQ=255) gctgctTTCTCCCGGCACTACCGTTTTACTCAAAATGCC‑GGGGggtggt > 1:1060588/1‑49 (MQ=255) gctgctTTCTCCCGGCACTACCGTTTTACTCAAAATGCC‑GGGGggtggt > 1:899703/1‑49 (MQ=255) gctgctTTCTCCCGGCACTACCGTTTTACTCAAAATGCC‑GGGGggtggt > 1:137854/1‑49 (MQ=255) gctgctTTCTCCCGGCACTACCGTTTTACTCAAAATGCC‑GGGGggtggt > 1:673756/1‑49 (MQ=255) gctgctTTCTCCCGGCACTACCGTTTTACTCAAAATGCC‑GGGGGGTGGTa < 1:613755/50‑1 (MQ=255) ctgctTTCTCCCGGCACTACCGTTTTACTCAAAATGCCGGGGGGGTGGTa < 1:3075531/50‑1 (MQ=255) ctgctTTCTCCCGGCACTACCGTTTTACTCAAAATGCC‑GGGGGGTGGTaa > 1:3683962/1‑50 (MQ=255) tgctTTCTCCCGGCACTACCGTTTTACTCAAAATGCC‑GGGGGGTGGTAAt > 1:3665692/1‑50 (MQ=255) gctTTCTCCCGGCACTACCGTTTTACTCAAAATGCCGGGGGGGTGGTAAt > 1:3891067/1‑50 (MQ=255) ctTTCTCCCGGCACTACCGTTTTACTCAAAATGCC‑GGGGGGTGGTAAt < 1:462107/48‑1 (MQ=255) ctTTCTCCCGGCACTACCGTTTTACTCAAAATGCC‑GGGGGGTGGTAATgg < 1:841814/50‑1 (MQ=255) ctTTCTCCCGGCACTACCGTTTTACTCAAAATGCC‑GGGGGGTGGTAATgg < 1:3542833/50‑1 (MQ=255) tttCTCCCGGCACTACCGTTTTACTCAAAATGCC‑GGGGGGTGGTAATgg > 1:815406/1‑49 (MQ=255) tttCTCCCGGCACTACCGTTTTACTCAAAATGCC‑GGGGGGTGGTAATgg > 1:449856/1‑49 (MQ=255) tctcCCGGCACTACCGTTTTACTCAAAATGCCGGGGGGGTGGTAATggg < 1:1452020/49‑1 (MQ=255) tctcCCGGCACTACCGTTTTACTCAAAATGCCGGGGGGGTGGTAATGGGa > 1:760805/1‑50 (MQ=255) ctcCCGGCACTACCGTTTTACTCAAAATGCCGGGGGGGTGGTAATGGGa < 1:2309514/49‑1 (MQ=255) ctcCCGGCACTACCGTTTTACTCAAAATGCC‑GGGGGGTGGTAATGGGAtc > 1:2519915/1‑50 (MQ=255) ccGGCACTACCGTTTTACTCAAAATGCC‑GGGGGGTGGTAATGGGAtct < 1:3439752/48‑1 (MQ=255) cGGCACTACCGTTTTACTCAAAATGCC‑GGGGGGTGGTAATGGGAtctc > 1:4339062/1‑48 (MQ=255) cGGCACTACCGTTTTACTCAAAATGCC‑GGGGGGTGGTAATGGGAtctc > 1:3907820/1‑48 (MQ=255) cGGCACTACCGTTTTACTCAAAATGCC‑GGGGGGTGGTAATGGGATCTCCa < 1:3926770/50‑1 (MQ=255) gCACTACCGTTTTACTCAAAATGCCGGGGGGGTGGTAATGGGATCTCCAg < 1:4140757/50‑1 (MQ=255) cACTACCGTTTTACTCAAAATGCCGGGGGGGTGGTAATGGGATCTCCAGt > 1:898264/1‑50 (MQ=255) aCTACCGTTTTACTCAAAATGCC‑GGGGGGTGGTAATGGGATCTCCAGt > 1:1820995/1‑48 (MQ=255) cTACCGTTTTACTCAAAATGCCGGGGGGGTGGTAATGGGATCTCCAGTga < 1:4120484/50‑1 (MQ=255) cTACCGTTTTACTCAAAATGCC‑GGGGGGTGGTAATGGGATCTCCAGt > 1:1102816/1‑47 (MQ=255) aCCGTTTTACTCAAAATGCCGGGGGGGTGGTAATGGGATCTCCAGTga > 1:704692/1‑48 (MQ=255) aCCGTTTTACTCAAAATGCC‑GGGGGGTGGTAATGGGATCTCCAGTga > 1:2381995/1‑47 (MQ=255) ccGTTTTACTCAAAATGCC‑GGGGGGTGGTAATGGGATCTCCAGTga > 1:4371113/1‑46 (MQ=255) cGTTTTACTCAAAATGCCGGGGGGGTGGTAATGGGATCTCCAGTgagaga < 1:2069866/50‑1 (MQ=255) cGTTTTACTCAAAATGCC‑GGGGGGTGGTAATGGGATCTCCAGTga > 1:3004940/1‑45 (MQ=255) cGTTTTACTCAAAATGCC‑GGGGGGTGGTAATGGGATCTCCAGTGAGAGAt < 1:882815/50‑1 (MQ=255) gTTTTACTCAAAATGCCGGGGGGGTGGTAATGGGATCTCCAGTga > 1:3867798/1‑45 (MQ=255) gTTTTACTCAAAATGCCGGGGGGGTGGTAATGGGATCTCCAGTGAGAGAt < 1:3028388/50‑1 (MQ=255) gTTTTACTCAAAATGCC‑GGGGGGTGGTAATGGGATCTCCAGTga > 1:1689960/1‑44 (MQ=255) ttttACTCAAAATGCCGGGGGGGTGGTAATGGGATCTCCAGTga > 1:3651947/1‑44 (MQ=255) ttACTCAAAATGCC‑GGGGGGTGGTAATGGGATCTCCAGTGAGAGATCg > 1:2721229/1‑48 (MQ=255) ttACTCAAAATGCC‑GGGGGGTGGTAATGGGATCTCCAGTGAGAGATCGa < 1:4357811/49‑1 (MQ=255) aCTCAAAATGCCGGGGGGGTGGTAATGGGATCTCCAGTGAGAGATCGAAg > 1:2221077/1‑50 (MQ=255) aCTCAAAATGCC‑GGGGGGTGGTAATGGGATCTCCAGTGAGAGATCGAAGc > 1:2757279/1‑50 (MQ=255) cTCAAAATGCCGGGGGGGTGGTAATGGGATCTCCAGTGAGAGATCGAAGc < 1:402933/50‑1 (MQ=255) cTCAAAATGCC‑GGGGGGTGGTAATGGGATCTCCAGTGAGAGATCGAAGCg > 1:985015/1‑50 (MQ=255) aaaaTGCC‑GGGGGGTGGTAATGGGATCTCCAGTGAGAGATCGAAGCGGt < 1:4220965/49‑1 (MQ=255) aaaaTGCC‑GGGGGGTGGTAATGGGATCTCCAGTGAGAGATCGAAGCGGTc > 1:4035502/1‑50 (MQ=255) aaaTGCCGGGGGGGTGGTAATGGGATCTCCAGTGAGAGATCGAAGCGGTc < 1:321753/50‑1 (MQ=255) aTGCCGGGGGGGTGGTAATGGGATCTCCAGTGAGAGATCGAAGCGGTCga < 1:1080077/50‑1 (MQ=255) cc‑GGGGGGTGGTAATGGGATCTCCAGTGAGAGATCGAAGCGGTCGAGaa < 1:1629787/49‑1 (MQ=255) | AACGGTGGCGCTGCTTTCTCCCGGCACTACCGTTTTACTCAAAATGCC‑GGGGGGTGGTAATGGGATCTCCAGTGAGAGATCGAAGCGGTCGAGAA > NC_000913/3948752‑3948846 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |