| Predicted mutation | ||||||
|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | annotation | gene | description |
| RA | NC_000913 | 2,652,855 | (G)6→7 | coding (362/846 nt) | sseA → | 3‑mercaptopyruvate sulfurtransferase |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 2,652,849 | 1 | . | G | 94.0% | 164.5 / 5.5 | 50 | coding (356/846 nt) | sseA | 3‑mercaptopyruvate sulfurtransferase |
| Reads supporting (aligned to +/- strand): ref base . (2/1); new base G (20/27); total (22/28) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 5.76e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.12e-01 | |||||||||||
GGATAAGCACCTGATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGAT‑TCT‑GGGGG‑GTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAACGCCGCGTTTAAT > NC_000913/2652751‑2652943 | ggATAAGCACCTGATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACATCTGGTGTAGAGAAAGTGTCGAT‑TCT‑g > 2:449419/1‑100 (MQ=255) gATAAGCACCTGATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGAT‑TCT‑gg < 2:1286253/100‑1 (MQ=255) ttccgatctaTTGTCTATGACGAAGGTAATTTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGAT‑TCT‑ggggg < 2:648254/91‑1 (MQ=255) aGCACCTGATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGAT‑TCT‑ggggg > 1:215947/1‑99 (MQ=255) gCACCTGATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGAT‑TCT‑ggggg‑gg > 2:1123936/1‑99 (MQ=255) aCCTGATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGAT‑TCT‑ggggg‑ggtg < 2:720784/100‑4 (MQ=255) tGATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGAT‑TCTGGGGGG‑GTGGAc > 1:872761/1‑100 (MQ=255) tGATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGAT‑TCT‑GGGGG‑GTGGACt < 1:1779831/100‑1 (MQ=255) aTTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGAT‑TCT‑gggggagatcggaa > 1:648250/1‑91 (MQ=255) gTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGAT‑TCTGGGGGG‑GTGGACTTGCa > 2:1116249/1‑100 (MQ=255) gTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGAT‑TCTGGGGGG‑GTGGACTTGCa < 1:944694/100‑1 (MQ=255) tATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGAT‑TCTGGGGGG‑GTGGACTTGCAGGc < 1:213050/100‑1 (MQ=255) tATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGAT‑TCTGGGGGG‑GTGGACTTGCAGGc > 1:1217755/1‑100 (MQ=255) tGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGGGCACCTTTGGTGTAGAGAAAGTGTCGAT‑TCTGGGGGG‑GTGGACTTGCAGGCTg < 2:90641/100‑1 (MQ=255) gACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGAT‑TCTGGGGGG‑GTGGACTTGCAGGCTgg < 2:1347334/100‑1 (MQ=255) gAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGAT‑TCTGGGGGG‑GTGGACTTGCAGGCTGGCAg < 1:1743310/100‑1 (MQ=255) aaGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGAT‑TCTGGGGGG‑GTGGACTTGCAGGCTGGCAgc > 2:639293/1‑100 (MQ=255) gTAATCTTTTGTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGAT‑TCTGGGGGG‑GTGGACTTGCAGGCTGGCAGcgcg < 1:1188594/100‑1 (MQ=255) gTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGAT‑TCTGGGGGG‑GTGGACTTGCAGGCTGGCAGcgcg > 1:1837/1‑100 (MQ=255) tAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGAT‑TCTGGGGGG‑GTGGACTTGCAGGCTGGCAGCGCga > 2:915647/1‑100 (MQ=255) aaTCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGAT‑TCTGGGGGG‑GTGGACTTGCAGGCTGGCAGCGCgat < 1:915641/100‑1 (MQ=255) aTCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGAT‑TCTGGGGGG‑GTGGACTTGCAGGCTGGCAGCGCgatg < 1:587953/100‑1 (MQ=255) cTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGAT‑TCTGGGGGG‑GTGGACTTGCAGGCTGGCAGCGCgatgat > 2:236559/1‑100 (MQ=255) cTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGAT‑TCTGGGGGG‑GTGGACTTGCAGGCTGGCAGCGCgatgat < 2:601717/100‑1 (MQ=255) ttttCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGAT‑TCTGGGGGG‑GTGGACTTGCAGGCTGGCAGCGCGATGATc < 2:406560/100‑1 (MQ=255) ttttCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGAT‑TCTGGGGGG‑GTGGACTTGCAGGCTGGCAGCGCGATGATc < 2:1648001/100‑1 (MQ=255) tttCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGAT‑TCTGGGGGG‑GTGGACTTGCAGGCTGGCAGCGCGATGATct < 1:1332001/100‑1 (MQ=255) tttCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGAT‑TCTGGGGGG‑GTGGACTTGCAGGCTGGCAGCGCGATGATct > 2:1590862/1‑100 (MQ=255) gCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGAT‑TCTGGGGGG‑GTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTa < 1:236558/100‑1 (MQ=255) gCCCCACGAGCATGGTGAATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGAT‑TCTGGGGGG‑GTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGt < 1:1729025/98‑1 (MQ=255) ccccACGAGCATGGTGGATGATGCGCACCTTTGGTGTAGAGCAAGTGTCGAT‑TCT‑GGGGGAGTAAACTTGAAGGCTAGCAGCGCGATGATCTGCTGTTag > 2:497687/1‑100 (MQ=255) nccACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGAT‑TCTGGGGGG‑GTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTaga > 1:281248/2‑100 (MQ=255) aGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGAT‑TCTGGGGGG‑GTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAgg < 2:648041/100‑1 (MQ=255) gCATGGTGGTTGTTGGGCACCTTTGGTGTAGAGAAAGGGTCGAT‑TCTGGGGGG‑GTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAgaagaa < 1:1329391/97‑1 (MQ=255) gCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGAT‑TCTGGGGGG‑GTGGACTTGCAGGCTGGCAGCGCGATGAACTACTGTTAGAAGAAGGt > 2:523238/1‑100 (MQ=255) aTGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGAT‑TCTGGGGGG‑GTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGc < 2:1306911/100‑1 (MQ=255) tggtggATGCTGCGCACCTTTGGAGAAAAGAAAGTATCGATCACG‑GGGGG‑GCGGACTTCCAGGCACGCAGCACGCTGATCTGCTGTTagaagta > 1:586294/1‑92 (MQ=255) gtggATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGAT‑TCTGGGGGG‑GTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGt < 1:1798893/100‑1 (MQ=255) ggATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGAT‑TCTGGGGGG‑GTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTag > 1:283155/1‑100 (MQ=255) aTGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGAT‑TCTGGGGGG‑GTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTagag > 2:980392/1‑100 (MQ=255) tgctgcGCACCTTTGGTGTAGAGAAAGTGTCGAT‑TCTGGGGGG‑GTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGc < 2:180610/100‑1 (MQ=255) gctgcGCCCCTTTGGTGTAGCGAAAGTGTCGAT‑TCTGGGGGG‑GTGGAATTGAAGGCTGGCAGCGCGATGATCAGCGGATAGAAGAAGGTGCGGTagag > 1:1583528/1‑98 (MQ=255) ctgcGCACCTTTGGTGTAGAGAAAGTGTCGAT‑TCTGGGGGG‑GTGGACTTGCAGGCTGGCAGCGCGATAATCTGCTGTTAGAAGAAAGAGCAGTAGAGCTg > 2:359247/1‑100 (MQ=255) ctgcGCACCTTTGGTGTAGAGAAAGTGTCGAT‑TCTGGGGGG‑GTGGACTTGCAGGCTGGCAGCACGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTg > 2:90677/1‑100 (MQ=255) cgcACCTTTGGTGTAGAGAAAGTGTCGAT‑TCTGGGGGG‑GTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCg < 2:838187/100‑1 (MQ=255) cACCTTTGGTGTAGAGAAAGTGTCGAT‑TCTGGGGGG‑GTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGa < 2:1212290/100‑1 (MQ=255) cACCTTTGGTGTAGAGAAAGTGTCGAT‑TCTGGGGGG‑GTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGa < 2:1323935/100‑1 (MQ=255) cTTTGGTGTAGAGAAAGTGTCGAT‑TCTGGGGGG‑GTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAgg > 2:1207860/1‑100 (MQ=255) tttGGTGTAGAGAAAGTGTCGAT‑TCTGGGGGG‑GTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAgg < 2:1461232/99‑1 (MQ=255) gtAGAGAAAGTGTCGAT‑TCTGGGGGG‑GTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGttt < 2:253180/100‑1 (MQ=255) tAGAGAAAGTGTCGAT‑TCTGGGGGG‑GTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTa > 1:1198547/1‑100 (MQ=255) agagAAAGTGTCGAT‑TCTGGGGGG‑GTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTaa > 2:7745/1‑100 (MQ=255) aaaGTGTCGAT‑TCTGGGGGG‑GTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAAcgcc < 2:397016/100‑1 (MQ=255) aaGTGTCGAT‑TCTGGGGGG‑GTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAAcgcc < 2:1184507/99‑1 (MQ=255) cGAT‑TCTGGGGGG‑GTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAACGCCGCGTTTa < 1:827338/100‑1 (MQ=255) gAT‑TCTGGGGGG‑GTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAACGCCACGTTTaa > 2:1765747/1‑100 (MQ=255) aT‑TCTGGGGGG‑GTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGCAAGGAGAGTTTAACGCCGCGTTTAAt > 2:429286/1‑100 (MQ=255) | GGATAAGCACCTGATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGAT‑TCT‑GGGGG‑GTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAACGCCGCGTTTAAT > NC_000913/2652751‑2652943 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |