| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | NC_000913 | 3,782,912 | A→C | 38.6% | V250G (GTG→GGG) | gpsA ← | glycerol‑3‑phosphate dehydrogenase (NAD+) |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 3,782,912 | 0 | A | C | 38.6% | 22.0 / 8.2 | 30 | V250G (GTG→GGG) | gpsA | glycerol‑3‑phosphate dehydrogenase (NAD+) |
| Reads supporting (aligned to +/- strand): ref base A (10/8); new base C (2/10); total (12/18) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 5.77e-02 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
GTAAAGTCAATGCTGCCTCGCGCGCGTTTTTTCCGCAATAT‑AATACTTGATAAATTTCCTCGGTTATTGGCATTTCAACGCCGAAGCGATGCGCCAGTTCGCGGACTTCTTTCGTATTGCGGTAGCCTTCCACCACCTGACCAATCTTCTCCTGCGCGCTTTGTACATCCATGCCCTGACCGAGCATCATGCCA‑AAACGGCGGT‑TACGCGACTGG‑TTGTCGGTACAGGTAA‑GCACCAGATCGCCAAG‑CCCCGCCATGCCCATAAAGGTGGCAGGGTCGGCACCCAGCGCCGCACCAAGACGCGACATTTCAGCCAGCCCACGGGTGATCAGCGCCGTACGCGCATTCGCACCAAAACCGATACCGTCGGACATCCCCGCACCAATGGCAATAACGTTTTTCACCGCGCCGCCAAGCTGCACGCCAATGAAATCCGGATTGCTGTAAACGCGGAAACTTTTGCCGCAGTGCA > NC_000913/3782679‑3783149 | ttAAAGTCATATCTGCCTCGCGGGCGTTTTTTCCGCAATTT‑AATAATTGATAAATTTCCTCGGTTATTGGCATTTCAACGCCGAAGCGATGCGCCAGTTCGCGGACTTCTTTCGTATTGCGGTAGCCTTCCACCACCTGACCAATCTTCTCCTGCGCGCTTTGTACATCCATGCCCTGACCGAGCATCATGCCA‑AAACGGCGGT‑TACGCGACTGG‑TTGTCGGTACAGGTAA‑GCACCAGa < 2:882799/238‑1 (MQ=255) aaTGCTGCCTCGCGCGCGTTTTTTCCGCAATAT‑AATACTTGATAAATTTCCTCGGTTATTGGCATTTCAACGCCGAAGCGATGCGCCAGTTCGCGGACTTCTTTCGTATTGCGGTAGCCTTCCACCACCTGACCAATCTTCTCCTGCGCGCTTTGTACATCCATGCCCTGACCGAGCATCATGCCA‑AAACGGCGGTCTACCCACCGGG‑TTGGCGG‑ACACGTAA‑GGACCaaaaacccac > 1:1078126/1‑229 (MQ=255) ccTCGCGCGCGTTTTTTCCGCAATAT‑AATACTTGATAAATTTCCTCGGTTATTGGCATTTCAACGCCGAAGCGATGCGCCAGTTCGCGGACTTCTTTCGTATTGCGGTAGCCTTCCACCACCTGACCAATCTTCTCCTGCGCGCTTTGTACATCCATGCCCTGACCGAGCATCATGCCA‑AAACGGCGGT‑TACACGACTGG‑GTGTTGGGACACGGAA‑GGACCaacacccccaacccccc > 1:342672/1‑222 (MQ=255) cttggcggtgTTTTTTCCAAAAAAT‑AATTCTTGAAAAATTGCCTGGGTTATTGGAATTTAAACGCCGAAGCGATGCGCCAGTTCGCGTATTTTTTTCGTATTGGGGTAGCCTTCCCCCCCCTGACCAATCTTCACCTGCGCGTTTTGTACATCCATGCCCTGACCGAGCATCATGCCA‑AAACGGCGGT‑TACGCGACTGG‑TTGTCGGTACAGGTAA‑GCACCAGATCGCCAAG‑CCCCGCCa < 2:69072/230‑1 (MQ=255) ctcgcggctgTTTTTTCCG‑ATTATAAATCTTTGATAATTTTCCTGGGTTATTGGAATTCAAACCCCGACGGGCTGCGCCAGTTCGCGGACTTTTTTCGTATTGCGGTAGCCTTCCCCCCCCTGCCCAACCTCCTCCTGCGCGCTTGGTACATCCATGCCCTGACCGAGCATCATGCCA‑AACCGGCGGT‑TACGCGACTGG‑TTGTCGGTACAGGTAA‑GCCCCAGATCGCCAAG‑CCCCGCCa < 2:869322/230‑1 (MQ=255) caaaatt‑aattCTTGGTAAATTTCCCCGGTTTTTGGCATTCCACCGCGGAGGCGATGCGCAGTTTCGCGCACTTCTTTATTATGGCGGTAGCCTGCCACCACCTGACAAATCTTCCCCTGAGGGTGTTGCACAGCCATGCCCTGACCGAGCATCATGCCA‑AACCGGCGGT‑TACGCGATTGG‑TTGCCGGTACAGGTAA‑GCCCCAGCTCGCCAAG‑CCCCGCCATGCCCATAAAGGTGGCAg < 2:212936/233‑1 (MQ=255) gATAAATTTCCTCGGTTATTGGCATTTCAACGCCGAAGCGATGCGCCAGTTCGCGGACTTCTTTCGTATTGCGGTAGCCTTCCACCACCTGACCAATCTTCTCCTGCGCGCTTTGTACATCCATGCCCTGACCGAGCATCATGCCA‑AAACGGCGGT‑TTCCCGCCTGGTTTGCCGGTA‑AGGTAACACCCCACAGCCCCAAC‑CCCCCCCAGCCCCACAAAGGGGGCAGGG‑CGGGAccccg > 1:1077605/1‑236 (MQ=255) gATAAATTTCCTCGGTTATTGGCATTTCAACGCCGAAGCGATGCGCCAGTTCGCGGACTTCTTTCGTATTGCGGTAGCCTTCCACCACCTGACCAATCTTCTCCTGCGCGCTTTGTACATCCATGCCCTGACCGAGCATCATGCCA‑AAACGGCGGT‑TACGCGACCGG‑GTGTCGGGACAGGTAA‑GACCCAAATCGCCAAA‑CCCCCCCATGCCCCTAAAAGGTGCAGGGTCggggcaccg > 1:679730/1‑231 (MQ=255) gATAAATTTCCTCGGTTATTGGCATTTCAACGCCGAAGCGATGCGCCAGTTCGCGGACTTCTTTCGTATTGCGGTAGCCTTCCACCACCTGACCAATCTTCTCCTGCGCGCTTTGTACATCCATGCCCTGACCGAGCATCATGCCA‑AAACGGCGGT‑TACGCCACAGG‑GTGTCGGTACAGGTAA‑GCACCACAACGCCCAG‑CCCCCCCCCGCCCCAAAAAAGGGGAGGGGCGGCACCCaac > 2:510728/1‑237 (MQ=255) tccccgccagggCGATGCGCCGGTTCGCGGCCTTTTTGCGTTTTGCGGTACCTTTCCCCCCCCTGACCAATTCTCCCCTGCGCTCTTTGTACCTCCAGGCCCTGACCGAGCATCAAGCAA‑AACCGGCGGT‑TCCGCGACTGG‑TTGTCGGTACAGGTAA‑GCCCCAGATCGCCAGG‑CCCCGCCATGCCCATAAAGGTGGCAGGGTCGGCACCCAGCGCCGCACCAAGACGCGACATTTcagc < 2:646471/228‑1 (MQ=255) cGAAGCGATGCCCCAGTTCGCGGACTTCTTTCGTATTGCGGTACCCTTCCCCCACCTGACCAATCTTCTCCTGCGCGCTTTGTACATCCATGCCCTGACCGAGCATCATGCCA‑AAACGGCGGT‑TACGCGACTGG‑TTGTCGGTACAGGTAA‑GCACCAGATCGCCAAG‑CCCCGCCATGCCCATAAAGGTGGCAGGGTCGGCACCCAGCGCCGCACCAAGACGCGACATTTCAGCCAGCCCa < 1:510728/239‑1 (MQ=255) gcggagcccccggttgcggtaCTTCTTTCCTTTTGCGGTTGCATCCCACCCCCTGACCAATCCTCCCGTGGGCCTTTTGTCAAGCTCTGCCCTGACCG‑GGGTCCTGGCACAACCGGCTGT‑TACGCGTCTGG‑TTGTCGGTACAGGTAA‑GCCCAAGATCGCCAAG‑CCCCGCCATGCCCATAAAGGTGGCGGGGGCGGCACCCAGCGCCGCACCAAGACGCGACATTTCAGCCAGCCCACgg < 2:192097/219‑1 (MQ=255) ccTTCCACCCCCTGACCAATCTTCTCCTGCGCGCTTTGTACATCCATGCCCTGACCGACCATCATGCCA‑AAACGGCGGT‑TACGCGACTGG‑TTGTCGGTACAGGTAA‑GCACCAGATCGCCAAG‑CCCCGCCATGCCCATAAAGGTGGCAGGGt < 1:23524/151‑1 (MQ=255) cTGCGCGCTTGGTACATCCATGCCCGGACCGAGCATCATGCCA‑AAACGGGGGT‑TACGCGACTGG‑TTGTCGGTACAGGAAA‑GCACCAGATCGCCAAG‑CCCCGCCATGCCCATAAAGGTGGCAGGGTCGGCACCCAGCGCCGCACCAAGACGCGACATTTCAGCCAGCCCACGGGTGATCAGCGCCGTACGCGCATTCGCACCAAAACCGATACCGTCGGACATCCCCGCACCAATGGCag < 1:662158/239‑2 (MQ=255) cgcTTTGTACATCCATGCCCTGACCGAGCATCATGCCA‑AAACGGCGGT‑TATGCGACTGG‑GTGTCGGTACCGGTAA‑GCACCAAAACGCCCAG‑CCCCCCCATGCCCCCAAAGAAGGGGGGGGCGGCCCCCACCCGCCCCCCCAGAAGCGCCACATCAGCCAGCCCCCCGGGCGTCATCGCCGCACGCGCCCCCGCCCCCAAACCCACACCGCCGGACACCCCCCCCCCACAGGCAaaaaac > 2:345522/1‑237 (MQ=255) tGTACATCCATGCCCTGACCGAGCATCATGCCA‑AAACGGCGGT‑TCCGCGACTGG‑TTGTCGGTACAGGTAA‑GCACCAGATCGCCAAG‑CCCCGCCATGCCCATAAAGGTGGCAGGGTCGGCACCCAGCGCCGCACCAAGACGCGACATTTCAGCCAGCCCACGGGTGAt < 1:294320/167‑1 (MQ=255) cactGCCGTGACCGGGCAACAGGCCA‑AAACGGCGGT‑TCCGCGACTGG‑TTGTCGGTAACGGTAA‑GCACCAGATCGCCAAG‑CCCCGCCATGCCCATAAAGGTGGCAGGGTCGGCACCCAGCGCCGCACCAAGACGCGACATTTCAGCCAGCCCACGGGTGATCAGCGCCGTACGCGCATTCGCACCAAAACCGATACCGTCGGACATCCCCGCACCAATGGCAATAACGTTTTTCACcgcg < 1:1072694/236‑1 (MQ=255) cccTGACCGAGCATCATGCCA‑AAACGGCGGT‑TTCGCGACTGG‑GTGTCGGGACAGGTAA‑GCACCCCATCGCCAAG‑CCCCGCCATTCCCACAAAGgggggcggg > 1:854340/1‑97 (MQ=255) ccTGACCGAGCATCATGCCA‑AAACGGGGGT‑TCCGCGACTGG‑TTGTCGGTACAGGTAA‑GCCCCAGATCGCCAAC‑CCCCGCCATGCCCATAAAGGTGGCAGGGTCGGCACCCAGCGCCGCACCAAGACGCGACATTTCAGCCAGCCCACGGGTGATCAGCGCCGTACGCGCATTCGCACCAAAACCGATACCGTCGGACATCCCCGCACCAATGGCAATAACGTTTTTCACCGCGCCGCCa < 1:246974/239‑1 (MQ=255) ccgggatacatGCAA‑AAGGGGGGGT‑TACGCGACTGG‑TTGTCGGTACAGGTAG‑GCCCCAGTTCGCCAAG‑CCCCGCCATGCCCATAAAGGTGGCAGGGTCGGCACCCAGCGCCGCACCAAGACGCGACATTTCAGCCAGCCCACGGGTGATCAGCGCCGTACGCGCATTCGCACCAAAACCGATACCGTCGGACATCCCCGCACCAATGGCAATAACGTTTTTCACCGCGCCGCCAAGCTg < 1:244848/231‑1 (MQ=255) gggcttcagGCCA‑AAACGGGGGT‑TCCGGGACTGG‑TTGCCGGTAAAGTAAA‑GCCCCAGATCCCAAAG‑CCCCGCCATGCCCATAAAGGTGGCAGGGTCGGCACCCAGCGCCGCACCAAGACGCGACATTTCAGCCAGCCCACGGGTGATCAGCGCCGTACGCGCATTCGCACCAAAACCGATACCGTCGGACATCCCCGCACCAATGGCAATAACGTTTTTCACCGCGCCGCCAAGCTGCa < 1:345522/234‑1 (MQ=255) gcaCTGGTTTGTGGGTACGGGTAG‑CCCCCAGA‑CCCCAAGCCCCCGCCATGCCCATAAAGGTGGCAGGGTCGGCACCCAGCGCCGCACCAAGACGCGACATTTCAGCCAGCCCACGGGTGATCAGCGCCGTACGCGCATTCGCACCAAAACCGATACCGTCGGACATCCCCGCACCAATGGCAATAACGTTTTTCACCGCGCCGCCAAGCTGCACGCCAATGAAATCCGGATTGCTGTaa < 1:445159/237‑1 (MQ=255) aaCAGGTAA‑GCCCCAGATCGCCAAG‑CCCCGCCATGCCCATAAAGGTGGCAGGGTCGGCACCCAGCGCCGCACCAAGACGCGACATTTCAGCCAGCCCACGGGTGATCAGCGCCGTACGCGCATTCGCACCAAAACCGATACCGTCGGACATCCCCGCACCAATGGCAATAACGTTTTTCACCGCGCCGCCAAGCTGCACGCCAATGaa < 1:366370/207‑1 (MQ=255) cAGGTAA‑GCACCAGATCGCCAAG‑CCCCGCCATGCCCATAAAGGTGGCAGGGTCGGCACCCAGCGCCGCACCAAGACGCGACATTTCAGCCAGCCCACGGGTGATCAGCGCCGTACGCGCATTCGCACCAAAACCGATACCGTCGGACATCCCCGCACCAATGGCCATAACGTTTTTCACCCCCCCGCCACGCTCCACGCCCATGCAATCCGGGTTGACTTAAACGCGaaaaacattgcc > 2:717874/1‑231 (MQ=255) cAGGTAA‑GCACCAGATCGCCAAG‑CCCCGCCATGCCCATAAAGGTGGCAGGGTCGGCACCCAGCGCCGCACCAAGACGCGACATTTCAGCCAGCCCACGGGTGATCAGCGCCGTACGCGCATTCGCACCAAAACCGATACCGTCGGACATCCCCGCACCAATGGCAATAACGTTTTTCACCGCGCCGCCAAGCTGCACGCCAATGAAATCCGGATTGCTGTAAACCCGGCAACTTTTGcc > 1:316574/1‑239 (MQ=255) gTAA‑GCCCCGAACCGCAAAG‑CCCCGCCATGCCCATAAAGGTGGCAGGGTCGGCACCCAGCGCCGCACCAAGACGCGACATTTCAGCCAGCCCACGGGTGATCAGCGCCGTACGCGCATTCGCACCAAAACCGATACCGTCGGACATCCCCGCACCAATGGCAATAACGTTTTTCACCGCGCCGCCAAGCTGCACGCCAATGAAATCCGGATTGCTGTAAACGCGGAAACTTTTGCCg < 1:35910/237‑1 (MQ=255) gTAA‑GCACCAGATCGCCAAG‑CCCCGCCATGCCCATAAAGGTGGCAGGGTCGGCACCCAGCGCCGCACCAAGACGCGACATTTCAGCCAGCCCACGGGTGATCAGCGCCGTACGCGCATTCGCACCAAAACCGATACCGTCGGACATCCCCGCACCAATGGCAATAACGTTTTTCACCGCGCCGCCAAGCTGCACGCCAATGAAATCCGGATTGCGCTAAACCAGGCAAATTTTTCCGcc > 1:20391/1‑238 (MQ=255) gaa‑GCACCGGATCCCCAGC‑CCCCGCCATGCCCATAAAGGTGGCAGGGTCGGCACCCAGCGCCGCACCAAGACGCGACATTTCAGCCAGCCCACGGGTGATCAGCGCCGTACGCGCATTCGCACCAAAACCGATACCGTCGGACATCCCCGCACCAATGGCAATAACGTTTTTCACCGCGCCGCCAAGCTGCACGCCAATGAAATCCGGATTGCTGTAAACGCGGAAACTTTTGCCg < 1:1035229/235‑1 (MQ=255) cACCAGATCGCCAAG‑CCCCGCCATGCCCATAAAGGTGGCAGGGTCGGCACCCAGCGCCGCACCAAGACGCGACATTTCAGCCAGCCCACGGGTGATCAGCGCCGTACGCGCATTCGCACCAAAACCGATACCGTCGGCCATCCCCGCACCAATGGCAATAACGTTTTTCACCGCGCCGCCAAGCTGCACGCCAATg > 1:1055246/1‑196 (MQ=255) cACCAGATCGCCAAG‑CCCCGCCATGCCCATAAAGGTGGCAGGGTCGGCACCCAGCGCCGCACCAAGACGCGACATTTCAGCCAGCCCACGGGTGATCAGCGCCGTACGCGCATTCGCACCAAAACCGATACCGTCGGACATCCCCGCACCAATGGCAATAACGTTTTTAACCGCGCCGCCAAGCTGCACGCCAATGAAATCCGGCTTGCTGTAAACGCGGAAACTTTTGCCGCAGTgca > 1:1048000/1‑239 (MQ=255) | GTAAAGTCAATGCTGCCTCGCGCGCGTTTTTTCCGCAATAT‑AATACTTGATAAATTTCCTCGGTTATTGGCATTTCAACGCCGAAGCGATGCGCCAGTTCGCGGACTTCTTTCGTATTGCGGTAGCCTTCCACCACCTGACCAATCTTCTCCTGCGCGCTTTGTACATCCATGCCCTGACCGAGCATCATGCCA‑AAACGGCGGT‑TACGCGACTGG‑TTGTCGGTACAGGTAA‑GCACCAGATCGCCAAG‑CCCCGCCATGCCCATAAAGGTGGCAGGGTCGGCACCCAGCGCCGCACCAAGACGCGACATTTCAGCCAGCCCACGGGTGATCAGCGCCGTACGCGCATTCGCACCAAAACCGATACCGTCGGACATCCCCGCACCAATGGCAATAACGTTTTTCACCGCGCCGCCAAGCTGCACGCCAATGAAATCCGGATTGCTGTAAACGCGGAAACTTTTGCCGCAGTGCA > NC_000913/3782679‑3783149 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 24 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |