Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_000913 | 3,782,913 | A→C | 37.8% | 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,913 | 0 | A | C | 37.8% | 21.1 / 7.2 | 28 | V250G (GTG→GGG) | gpsA | glycerol‑3‑phosphate dehydrogenase (NAD+) |
Reads supporting (aligned to +/- strand): ref base A (9/8); new base C (2/9); total (11/17) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.15e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 |
GTAAAGTCAATGCTGCCTCGCGCGCGTTTTTTCCGCAATAT‑AATACTTGATAAATTTCCTCGGTTATTGGCATTTCAACGCCGAAGCGATGCGCCAGTTCGCGGACTTCTTTCGTATTGCGGTAGCCTTCCACCACCTGACCAATCTTCTCCTGCGCGCTTTGTACATCCATGCCCTGACCGAGCATCATGCCAAAACGGCGGT‑TACGCGACTGG‑TTGTCGGTACAGGTAA‑GCACCAGATCGCCAAG‑CCCCGCCATGCCCATAAAGGTGGCAGGGTCGGCACCCAGCGCCGCACCAAGACGCGACATTTCAGCCAGCCCACGGGTGATCAGCGCCGTACGCGCATTCGCACCAAAACCGATACCGTCGGACATCCCCGCACCAATGGCAATAACGTTTTTCACCGCGCCGCCAAGCTGCACGCCAATGAAATCCGGATTGCTGTAAACGCGGAAACTTTTGCCGCAGTGCA > NC_000913/3782680‑3783150 | ttAAAGTCATATCTGCCTCGCGGGCGTTTTTTCCGCAATTT‑AATAATTGATAAATTTCCTCGGTTATTGGCATTTCAACGCCGAAGCGATGCGCCAGTTCGCGGACTTCTTTCGTATTGCGGTAGCCTTCCACCACCTGACCAATCTTCTCCTGCGCGCTTTGTACATCCATGCCCTGACCGAGCATCATGCCAAAACGGCGGT‑TACGCGACTGG‑TTGTCGGTACAGGTAA‑GCACCAGa < 2:882799/238‑1 (MQ=255) aaTGCTGCCTCGCGCGCGTTTTTTCCGCAATAT‑AATACTTGATAAATTTCCTCGGTTATTGGCATTTCAACGCCGAAGCGATGCGCCAGTTCGCGGACTTCTTTCGTATTGCGGTAGCCTTCCACCACCTGACCAATCTTCTCCTGCGCGCTTTGTACATCCATGCCCTGACCGAGCATCATGCCAAAACGGCGGTCTACCCACCGGG‑TTGGCGG‑ACACGTAA‑GGACCaaaaacccac > 1:1078126/1‑229 (MQ=255) cttggcggtgTTTTTTCCAAAAAAT‑AATTCTTGAAAAATTGCCTGGGTTATTGGAATTTAAACGCCGAAGCGATGCGCCAGTTCGCGTATTTTTTTCGTATTGGGGTAGCCTTCCCCCCCCTGACCAATCTTCACCTGCGCGTTTTGTACATCCATGCCCTGACCGAGCATCATGCCAAAACGGCGGT‑TACGCGACTGG‑TTGTCGGTACAGGTAA‑GCACCAGATCGCCAAG‑CCCCGCCa < 2:69072/230‑1 (MQ=255) ctcgcggctgTTTTTTCCG‑ATTATAAATCTTTGATAATTTTCCTGGGTTATTGGAATTCAAACCCCGACGGGCTGCGCCAGTTCGCGGACTTTTTTCGTATTGCGGTAGCCTTCCCCCCCCTGCCCAACCTCCTCCTGCGCGCTTGGTACATCCATGCCCTGACCGAGCATCATGCCAAACCGGCGGT‑TACGCGACTGG‑TTGTCGGTACAGGTAA‑GCCCCAGATCGCCAAG‑CCCCGCCa < 2:869322/230‑1 (MQ=255) caaaatt‑aattCTTGGTAAATTTCCCCGGTTTTTGGCATTCCACCGCGGAGGCGATGCGCAGTTTCGCGCACTTCTTTATTATGGCGGTAGCCTGCCACCACCTGACAAATCTTCCCCTGAGGGTGTTGCACAGCCATGCCCTGACCGAGCATCATGCCAAACCGGCGGT‑TACGCGATTGG‑TTGCCGGTACAGGTAA‑GCCCCAGCTCGCCAAG‑CCCCGCCATGCCCATAAAGGTGGCAg < 2:212936/233‑1 (MQ=255) gATAAATTTCCTCGGTTATTGGCATTTCAACGCCGAAGCGATGCGCCAGTTCGCGGACTTCTTTCGTATTGCGGTAGCCTTCCACCACCTGACCAATCTTCTCCTGCGCGCTTTGTACATCCATGCCCTGACCGAGCATCATGCCAAAACGGCGGT‑TTCCCGCCTGGTTTGCCGGTA‑AGGTAACACCCCACAGCCCCAAC‑CCCCCCCAGCCCCACAAAGGGGGCAGGG‑CGGGAccccg > 1:1077605/1‑236 (MQ=255) gATAAATTTCCTCGGTTATTGGCATTTCAACGCCGAAGCGATGCGCCAGTTCGCGGACTTCTTTCGTATTGCGGTAGCCTTCCACCACCTGACCAATCTTCTCCTGCGCGCTTTGTACATCCATGCCCTGACCGAGCATCATGCCAAAACGGCGGT‑TACGCGACCGG‑GTGTCGGGACAGGTAA‑GACCCAAATCGCCAAA‑CCCCCCCATGCCCCTAAAAGGTGCAGGGTCggggcaccg > 1:679730/1‑231 (MQ=255) gATAAATTTCCTCGGTTATTGGCATTTCAACGCCGAAGCGATGCGCCAGTTCGCGGACTTCTTTCGTATTGCGGTAGCCTTCCACCACCTGACCAATCTTCTCCTGCGCGCTTTGTACATCCATGCCCTGACCGAGCATCATGCCAAAACGGCGGT‑TACGCCACAGG‑GTGTCGGTACAGGTAA‑GCACCACAACGCCCAG‑CCCCCCCCCGCCCCAAAAAAGGGGAGGGGCGGCACCCaac > 2:510728/1‑237 (MQ=255) tccccgccagggCGATGCGCCGGTTCGCGGCCTTTTTGCGTTTTGCGGTACCTTTCCCCCCCCTGACCAATTCTCCCCTGCGCTCTTTGTACCTCCAGGCCCTGACCGAGCATCAAGCAAAACCGGCGGT‑TCCGCGACTGG‑TTGTCGGTACAGGTAA‑GCCCCAGATCGCCAGG‑CCCCGCCATGCCCATAAAGGTGGCAGGGTCGGCACCCAGCGCCGCACCAAGACGCGACATTTcagc < 2:646471/228‑1 (MQ=255) cGAAGCGATGCCCCAGTTCGCGGACTTCTTTCGTATTGCGGTACCCTTCCCCCACCTGACCAATCTTCTCCTGCGCGCTTTGTACATCCATGCCCTGACCGAGCATCATGCCAAAACGGCGGT‑TACGCGACTGG‑TTGTCGGTACAGGTAA‑GCACCAGATCGCCAAG‑CCCCGCCATGCCCATAAAGGTGGCAGGGTCGGCACCCAGCGCCGCACCAAGACGCGACATTTCAGCCAGCCCa < 1:510728/239‑1 (MQ=255) ccTTCCACCCCCTGACCAATCTTCTCCTGCGCGCTTTGTACATCCATGCCCTGACCGACCATCATGCCAAAACGGCGGT‑TACGCGACTGG‑TTGTCGGTACAGGTAA‑GCACCAGATCGCCAAG‑CCCCGCCATGCCCATAAAGGTGGCAGGGt < 1:23524/151‑1 (MQ=255) cTGCGCGCTTGGTACATCCATGCCCGGACCGAGCATCATGCCAAAACGGGGGT‑TACGCGACTGG‑TTGTCGGTACAGGAAA‑GCACCAGATCGCCAAG‑CCCCGCCATGCCCATAAAGGTGGCAGGGTCGGCACCCAGCGCCGCACCAAGACGCGACATTTCAGCCAGCCCACGGGTGATCAGCGCCGTACGCGCATTCGCACCAAAACCGATACCGTCGGACATCCCCGCACCAATGGCag < 1:662158/239‑2 (MQ=255) cgcTTTGTACATCCATGCCCTGACCGAGCATCATGCCAAAACGGCGGT‑TATGCGACTGG‑GTGTCGGTACCGGTAA‑GCACCAAAACGCCCAG‑CCCCCCCATGCCCCCAAAGAAGGGGGGGGCGGCCCCCACCCGCCCCCCCAGAAGCGCCACATCAGCCAGCCCCCCGGGCGTCATCGCCGCACGCGCCCCCGCCCCCAAACCCACACCGCCGGACACCCCCCCCCCACAGGCAaaaaac > 2:345522/1‑237 (MQ=255) tGTACATCCATGCCCTGACCGAGCATCATGCCAAAACGGCGGT‑TCCGCGACTGG‑TTGTCGGTACAGGTAA‑GCACCAGATCGCCAAG‑CCCCGCCATGCCCATAAAGGTGGCAGGGTCGGCACCCAGCGCCGCACCAAGACGCGACATTTCAGCCAGCCCACGGGTGAt < 1:294320/167‑1 (MQ=255) cactGCCGTGACCGGGCAACAGGCCAAAACGGCGGT‑TCCGCGACTGG‑TTGTCGGTAACGGTAA‑GCACCAGATCGCCAAG‑CCCCGCCATGCCCATAAAGGTGGCAGGGTCGGCACCCAGCGCCGCACCAAGACGCGACATTTCAGCCAGCCCACGGGTGATCAGCGCCGTACGCGCATTCGCACCAAAACCGATACCGTCGGACATCCCCGCACCAATGGCAATAACGTTTTTCACcgcg < 1:1072694/236‑1 (MQ=255) cccTGACCGAGCATCATGCCAAAACGGCGGT‑TTCGCGACTGG‑GTGTCGGGACAGGTAA‑GCACCCCATCGCCAAG‑CCCCGCCATTCCCACAAAGgggggcggg > 1:854340/1‑97 (MQ=255) ccTGACCGAGCATCATGCCAAAACGGGGGT‑TCCGCGACTGG‑TTGTCGGTACAGGTAA‑GCCCCAGATCGCCAAC‑CCCCGCCATGCCCATAAAGGTGGCAGGGTCGGCACCCAGCGCCGCACCAAGACGCGACATTTCAGCCAGCCCACGGGTGATCAGCGCCGTACGCGCATTCGCACCAAAACCGATACCGTCGGACATCCCCGCACCAATGGCAATAACGTTTTTCACCGCGCCGCCa < 1:246974/239‑1 (MQ=255) ccgggatacatGCAAAAGGGGGGGT‑TACGCGACTGG‑TTGTCGGTACAGGTAG‑GCCCCAGTTCGCCAAG‑CCCCGCCATGCCCATAAAGGTGGCAGGGTCGGCACCCAGCGCCGCACCAAGACGCGACATTTCAGCCAGCCCACGGGTGATCAGCGCCGTACGCGCATTCGCACCAAAACCGATACCGTCGGACATCCCCGCACCAATGGCAATAACGTTTTTCACCGCGCCGCCAAGCTg < 1:244848/231‑1 (MQ=255) gggcttcagGCCAAAACGGGGGT‑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‑AATACTTGATAAATTTCCTCGGTTATTGGCATTTCAACGCCGAAGCGATGCGCCAGTTCGCGGACTTCTTTCGTATTGCGGTAGCCTTCCACCACCTGACCAATCTTCTCCTGCGCGCTTTGTACATCCATGCCCTGACCGAGCATCATGCCAAAACGGCGGT‑TACGCGACTGG‑TTGTCGGTACAGGTAA‑GCACCAGATCGCCAAG‑CCCCGCCATGCCCATAAAGGTGGCAGGGTCGGCACCCAGCGCCGCACCAAGACGCGACATTTCAGCCAGCCCACGGGTGATCAGCGCCGTACGCGCATTCGCACCAAAACCGATACCGTCGGACATCCCCGCACCAATGGCAATAACGTTTTTCACCGCGCCGCCAAGCTGCACGCCAATGAAATCCGGATTGCTGTAAACGCGGAAACTTTTGCCGCAGTGCA > NC_000913/3782680‑3783150 |
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: ‑ |