| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | minE | 525,695 | C→T | 100% | P65S (CCG→TCG) | moaD → | molybdopterin synthase, small subunit |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | minE | 525,695 | 0 | C | T | 91.6% | 169.9 / 4.9 | 59 | P65S (CCG→TCG) | moaD | molybdopterin synthase, small subunit |
| Reads supporting (aligned to +/- strand): ref base C (5/0); new base T (54/0); total (59/0) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01 | |||||||||||
| Rejected as polymorphism: Variant not supported by required number of reads on each strand. | |||||||||||
GCAAATTACTGGCTGCCGTCAACCAGACGCTGGTGAGTTTTGACCATCCGCTGACTGACGGCGACGAAGTAGCTTTCTTCCCGCCGGTAACCGGA > minE/525648‑525742 | gCAAATTACTGGCTGCCGTCACCCAGACGCTGGTGAGTTTTGACCATCCGCTGACTGACGGcgacga > 1:1551448/1‑67 (MQ=255)gCAAATTACTGGCTGCCGTCAACCAGACGCTGGTGAGTTTTGACCATCCGCTGACTGACGGcgacga > 1:581767/1‑67 (MQ=255)gCAAATTACTGGCTGCCGTCAACCAGACGCTGGTGAGTTTTGACCATCCGCTGACTGACGGcgacga > 1:313097/1‑67 (MQ=255)gCAAATTACTGGCTGCCGTCAACCAGACGCTGGTGAGTTTTGACCATCCGCTGACTGACGGcgacga > 1:1472449/1‑67 (MQ=255)gCAAATTACTGGCTGCCGTCAACCAGACGCTGGTGAGTTTTGACCATCCGCTGACTGACGGcgacga > 1:1942182/1‑67 (MQ=255) gCTGGTGAGTTTTGACCATTCGCTGACTGACGGcgac > 1:1652737/1‑37 (MQ=255) gCTGGTGAGTTTTGACCATTCGCTGACTGACGGCGACGAAGTGGCTttct > 1:1952049/1‑50 (MQ=39) gCTGGTGAGTTTTGACCATTCGCTGACTGACGGCGACGAAGTGGCTTTCTTCCCGCCGGTa > 1:2197006/1‑61 (MQ=255) gCTGGTGAGTTTTGACCATTCGCTGACTGACGGCGACGAAGTGGCTTTCTTCCCGCCGGTAACCgg > 1:1349821/1‑66 (MQ=255) gCTGGTGAGTTTTGACCATTCGCTGACTGACGGCGACGAAGTGGCTTTCTTCCCGCCGGTAACCgg > 1:438428/1‑66 (MQ=255) gCTGGTGAGTTTTGACCATTCGCTGACTGACGGCGACGAAGTGGCTTTCTTCCCGCCGGTAACCGGa > 1:382424/1‑67 (MQ=255) gCTGGTGAGTTTTGACCATTCGCTGACTGACGGCGACGAAGTGGCTTTCTTCCCGCCGGTAACCGGa > 1:514468/1‑67 (MQ=255) gCTGGTGAGTTTTGACCATTCGCTGACTGACGGCGACGAAGTGGCTTTCTTCCCGCCGGTAACCGGa > 1:2278748/1‑67 (MQ=255) gCTGGTGAGTTTTGACCATTCGCTGACTGACGGCGACGAAGTGGCTTTCTTCCCGCCGGTAACCGGa > 1:363693/1‑67 (MQ=255) gCTGGTGAGTTTTGACCATTCGCTGACTGACGGCGACGAAGTGGCTTTCTTCCCGCCGGTAACCGGa > 1:35120/1‑67 (MQ=255) gCTGGTGAGTTTTGACCATTCGCTGACTGACGGCGACGAAGTGGCTTTCTTCCCGCCGGTAACCGGa > 1:334923/1‑67 (MQ=255) gCTGGTGAGTTTTGACCATTCGCTGACTGACGGCGACGAAGTGGCTTTCTTCCCGCCGGTAACCGGa > 1:247595/1‑67 (MQ=255) gCTGGTGAGTTTTGACCATTCGCTGACTGACGGCGACGAAGTGGCTTTCTTCCCGCCGGTAACCGGa > 1:238012/1‑67 (MQ=255) gCTGGTGAGTTTTGACCATTCGCTGACTGACGGCGACGAAGTGGCTTTCTTCCCGCCGGTAACCGGa > 1:2171824/1‑67 (MQ=255) gCTGGTGAGTTTTGACCATTCGCTGACTGACGGCGACGAAGTGGCTTTCTTCCCGCCGGTAACCGGa > 1:532078/1‑67 (MQ=255) gCTGGTGAGTTTTGACCATTCGCTGACTGACGGCGACGAAGTGGCTTTCTTCCCGCCGGTAACCGGa > 1:627440/1‑67 (MQ=255) gCTGGTGAGTTTTGACCATTCGCTGACTGACGGCGACGAAGTGGCTTTCTTCCCGCCGGTAACCGGa > 1:681850/1‑67 (MQ=255) gCTGGTGAGTTTTGACCATTCGCTGACTGACGGCGACGAAGTGGCTTTCTTCCCGCCGGTAACCGGa > 1:706513/1‑67 (MQ=255) gCTGGTGAGTTTTGACCATTCGCTGACTGACGGCGACGAAGTGGCTTTCTTCCCGCCGGTAACCGGa > 1:743981/1‑67 (MQ=255) gCTGGTGAGTTTTGACCATTCGCTGACTGACGGCGACGAAGTGGCTTTCTTCCCGCCGGTAACCGGa > 1:778838/1‑67 (MQ=255) gCTGGTGAGTTTTGACCATTCGCTGACTGACGGCGACGAAGTGGCTTTCTTCCCGCCGGTAACCGGa > 1:789795/1‑67 (MQ=255) gCTGGTGAGTTTTGACCATTCGCTGACTGACGGCGACGAAGTGGCTTTCTTCCCGCCGGTAACCGGa > 1:796838/1‑67 (MQ=255) gCTGGTGAGTTTTGACCATTCGCTGACTGACGGCGACGAAGTGGCTTTCTTCCCGCCGGTAACCGGa > 1:799899/1‑67 (MQ=255) gCTGGTGAGTTTTGACCATTCGCTGACTGACGGCGACGAAGTGGCTTTCTTCCCGCCGGTAACCGGa > 1:834307/1‑67 (MQ=255) gCTGGTGAGTTTTGACCATTCGCTGACTGACGGCGACGAAGTGGCTTTCTTCCCGCCGGTAACCGGa > 1:841748/1‑67 (MQ=255) gCTGGTGAGTTTTGACCATTCGCTGACTGACGGCGACGAAGTGGCTTTCTTCCCGCCGGTAACCGGa > 1:878334/1‑67 (MQ=255) gCTGGTGAGTTTTGACCATTCGCTGACTGACGGCGACGAAGTGGCTTTCTTCCCGCCGGTAACCGGa > 1:901554/1‑67 (MQ=255) gCTGGTGAGTTTTGACCATTCGCTGACTGACGGCGACGAAGTGGCTTTCTTCCCGCCGGTAACCGGa > 1:93079/1‑67 (MQ=255) gCTGGTGAGTTTTGACCATTCGCTGACTGACGGCGACGAAGTGGCTTTCTTCCCGCCGGTAACCGGa > 1:1456310/1‑67 (MQ=255) gCTGGTGAGTTTTGACCATTCGCTGACTGACGGCGACGAAGTGGCTTTCTTCCCGCCGGTAACCGGa > 1:103100/1‑67 (MQ=255) gCTGGTGAGTTTTGACCATTCGCTGACTGACGGCGACGAAGTGGCTTTCTTCCCGCCGGTAACCGGa > 1:1093394/1‑67 (MQ=255) gCTGGTGAGTTTTGACCATTCGCTGACTGACGGCGACGAAGTGGCTTTCTTCCCGCCGGTAACCGGa > 1:1094911/1‑67 (MQ=255) gCTGGTGAGTTTTGACCATTCGCTGACTGACGGCGACGAAGTGGCTTTCTTCCCGCCGGTAACCGGa > 1:1098840/1‑67 (MQ=255) gCTGGTGAGTTTTGACCATTCGCTGACTGACGGCGACGAAGTGGCTTTCTTCCCGCCGGTAACCGGa > 1:1120447/1‑67 (MQ=255) gCTGGTGAGTTTTGACCATTCGCTGACTGACGGCGACGAAGTGGCTTTCTTCCCGCCGGTAACCGGa > 1:1126192/1‑67 (MQ=255) gCTGGTGAGTTTTGACCATTCGCTGACTGACGGCGACGAAGTGGCTTTCTTCCCGCCGGTAACCGGa > 1:1158416/1‑67 (MQ=255) gCTGGTGAGTTTTGACCATTCGCTGACTGACGGCGACGAAGTGGCTTTCTTCCCGCCGGTAACCGGa > 1:1165349/1‑67 (MQ=255) gCTGGTGAGTTTTGACCATTCGCTGACTGACGGCGACGAAGTGGCTTTCTTCCCGCCGGTAACCGGa > 1:1214960/1‑67 (MQ=255) gCTGGTGAGTTTTGACCATTCGCTGACTGACGGCGACGAAGTGGCTTTCTTCCCGCCGGTAACCGGa > 1:1218728/1‑67 (MQ=255) gCTGGTGAGTTTTGACCATTCGCTGACTGACGGCGACGAAGTGGCTTTCTTCCCGCCGGTAACCGGa > 1:122684/1‑67 (MQ=255) gCTGGTGAGTTTTGACCATTCGCTGACTGACGGCGACGAAGTGGCTTTCTTCCCGCCGGTAACCGGa > 1:1375731/1‑67 (MQ=255) gCTGGTGAGTTTTGACCATTCGCTGACTGACGGCGACGAAGTGGCTTTCTTCCCGCCGGTAACCGGa > 1:227453/1‑67 (MQ=255) gCTGGTGAGTTTTGACCATTCGCTGACTGACGGCGACGAAGTGGCTTTCTTCCCGCCGGTAACCGGa > 1:1597403/1‑67 (MQ=255) gCTGGTGAGTTTTGACCATTCGCTGACTGACGGCGACGAAGTGGCTTTCTTCCCGCCGGTAACCGGa > 1:1716885/1‑67 (MQ=255) gCTGGTGAGTTTTGACCATTCGCTGACTGACGGCGACGAAGTGGCTTTCTTCCCGCCGGTAACCGGa > 1:1846929/1‑67 (MQ=255) gCTGGTGAGTTTTGACCATTCGCTGACTGACGGCGACGAAGTGGCTTTCTTCCCGCCGGTAACCGGa > 1:1961647/1‑67 (MQ=255) gCTGGTGAGTTTTGACCATTCGCTGACTGACGGCGACGAAGTGGCTTTCTTCCCGCCGGTAACCGGa > 1:2000107/1‑67 (MQ=255) gCTGGTGAGTTTTGACCATTCGCTGACTGACGGCGACGAAGTGGCTTTCTTCCCGCCGGTAACCGGa > 1:2079440/1‑67 (MQ=255) gCTGGTGAGTTTTGACCATTCGCTGACTGACGGCGACGAAGTGGCTTTCTTCCCGCCGGTAACCGGa > 1:2091982/1‑67 (MQ=255) gCTGGTGAGTTTTGACCATTCGCTGACTGACGGCGACGAAGTGGCTTTCTTCCCGCCGGTAACCGGa > 1:2106580/1‑67 (MQ=255) gCTGGTGAGTTTTGACCATTCGCTGACTGACGGCGACGAAGTGGCTTTCTTCCCGCCGGTAACCGGa > 1:2164514/1‑67 (MQ=255) gCTGGTGAGTTTTGACCATTCGCTGACTGACGGCGACGAAGTGGCTTTCTTCCCGCCGGTAACCGGa > 1:1008981/1‑67 (MQ=255) gCTGGTGAGTTTTGACCATTCGCTGACTGACGGCGACGAAGTGGCTTTCTTCCCGCCGGTAACCGGa > 1:222/1‑67 (MQ=255) gCTGGTGAGTTTTGACCATTCGCTGACTGACGGCGACGAAGTGGCTTTCTTCCCGCCGGTAACCGGa > 1:226759/1‑67 (MQ=255) | GCAAATTACTGGCTGCCGTCAACCAGACGCTGGTGAGTTTTGACCATCCGCTGACTGACGGCGACGAAGTAGCTTTCTTCCCGCCGGTAACCGGA > minE/525648‑525742 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |