| Predicted mutation | ||||||
|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | annotation | gene | description |
| RA | NC_000913 | 3,518,683 | G→A | P128S (CCG→TCG) | aroK ← | shikimate kinase 1 |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 3,518,683 | 0 | G | A | 91.7% | 22.8 / ‑4.1 | 13 | P128S (CCG→TCG) | aroK | shikimate kinase 1 |
| Reads supporting (aligned to +/- strand): ref base G (1/0); new base A (5/6); total (7/6) | |||||||||||
| 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.41e-01 | |||||||||||
AATTGGGTAACTACGTTCCCCGAGAGTAACGACAATCCTCTCCATAACGCGACATCCACCTTAATTACTGTACCCGCAGACGAGTGTATATAAAGCCAGAATTAGTTGCTTTCCAGCATGTGAATAATCTGGTTTGCAACCACTTTAGCGCTTTGATCATCAGTACGAATGGTCACGTCGGCAATCTCTTCATACAGCGGATTGCGTTCATTGGCCAACGCTTCCAGAACTTCACGCGGCGGTGTTTCAACGTGCAGCAACGGGCGTTTTTTATCACGCTGCGTGCGTGCAAGTTGCTTTTCGATGGTCGTTTCAAGATAAACGACAACGCCACGAGCGGAAAGACGGTTACGCGTTTCACGGGATTTCACAGAGCCGCCGCCAGTAGCCAGCACAATACCCTGTTTCTCGGTCAACTCATTGATGACCTTTTCTTCGCGATCGCGGAAGCCTTCTTCGCCTTCTAAATCGAAAACCCAG > NC_000913/3518442‑3518921 | aaTTGGGTAACTACGTTCCCCGAGAGTAACGACAATCCTCTCCATAACGCGACATCCACCTTAATTACTGTACCCGCAGACGAGTGTATATAAAGCCAGAATTAGTTGCTTTCCAGCATGTTAATAATCTGGTTTGCATCCACTTTAGCGCTTTGCTCATCAGTCCGAATGGTCACGTCGGCAATCTCTTCATACAGCGGATTGCGTTCATTGTCCACCGCTTTCAGAACTTCACGCGGCGATGTTTCaan < 1:90515/251‑2 (MQ=255) ggTAACTACGTTCCCCGAGAGTAACGACACTCCTCTCAATAACGCGACATCCACCTTACTTACTGTACCCGCCGACGAGTGTATATAAAGCAAGAATTAGTTGCTTTCCAGCATGTGAATAATCTGGTTTGCAACCACTTTAGCGCTTTGATCATCAGTACGAATGGTCACGTCGGCAATCTCTTCATACAGCGGATTGCGTTCATTGGCCAACGCTTCCAGAACTTCACGCGGCGGTGTTTCAACGTgc > 1:448126/1‑250 (MQ=255) ccGAGAGTAACGACAATCCTCTCCATAACGCGACATCCACCTTAATTACTGTACCCGCAGACGAGTGTATATAAAGCCAGAATTAGTTGCTTTCCAGCATGTGAATAATCTGGTTTGCAACCACTTTAGCGCTTTGAGCATCAGGACGAATGGTCACGTCGGCAATCTCTTCATACATCGGATTGCGTTCATTGGGCAACGCTTCGAGAACTTCACGCGGCGATGTTTCAACGTGCAGCAACggn < 1:45860/245‑2 (MQ=255) ataAAGCCAGAATTAGTTGCTTTCCAGCATGTGAATAATCTGGTTTGCAACCTCTTTAGCGATTTGATCATCAGTACGAATGGTCACGTCGGCAATCTCTTCATACAGCGGATTGCGTTCATTGGCCAACGCTTCCAGAACTTCACGCGGCGATGTTTCAACGTGCAGCAACGGGCGTTTTTTAt > 1:230154/1‑185 (MQ=255) tttCCAGCATGGGAATAATCTGGTTTGCAACCACTTTCGCCCTTTGATCCTCGATACGAACAGTCAAGTCAGCAATCCCTTCATACAGCGAATTCCCATCCTTGGCCAACACTTCCAGAACATCACGCGGCGATGTTTCAACGTGCAGCAACGGGCGTTTTTTATCACGCTGCGTGCGTGCAAGATGCTTTTCGATGGTCGGTTCAAGATAAACGACAACCCCACGAGCGGAAACACGGTTACGCGTTTc > 1:246117/1‑250 (MQ=255) gAATAATCTGGTTTGCAACCACTTTAGCGCTTTGATCATCAGTACGAATGGTCACGGCGGCAATCTCTTCATACAGCGGATTGCGTTCATTGGCCAACGCTTCCAGAACTTCACGCGGCGATGTTTCAACGTGCAGCAACGGGCGTTTTTTATCACGCTGCGTGCGTGCAAGTTGCTTTTCGATGGTCGTTTCAAGATAAACTACAACGCCACGAGCGGAAAGACGGTTAcgcgn < 1:46711/235‑2 (MQ=255) ggTTTGCAACCACGTTAGCGCTTTGATCATCAGTACGAAGGGTCACGTCGGAAATCTCTTCATACAGCGGATTGCGTTCATTGGCCAACGCTTCCGGAACTTCACGCGGCGATTTTTCAACTTGCATCAACGGTCGTTTTTTATCACGCTGCGTGCGTGCAATTTGCGTTTCGATTGTCGTTTCAAGATAAACGACAAGGCCACGAGCGGAAAGACGGTTACGCGTTTCACGGGATTTCACAGAgccgct < 1:66922/250‑2 (MQ=255) aCCACTTTAGCGCTTTGATAATCAGTACGAATGGTCACGTCGGCAATATCTTCATACAGCGGATTGCGTTCATTGGCCAACGCTTCAAGAACTTCACGCGGCGATGTTTCAACGTGCAGCAACGGGCGTTTTTTATCACGATGCGTGCGTGCAAGTTGCTTTTCGATGGTCGTTTCAAGATAAACGACAACGCCACGAGCGGAAAGACGGTTACGAGTTTAACGGGATTTcac > 1:168000/1‑233 (MQ=255) ttgttaCGTCGGCAATCGCTGCATACAGCTGATTGCGTTCAGTGGCCAACGCTTCCAGAACCTCACGCTGCGATGGTGCAACGTGCAGCAACTGGTGTTTTTTATCATGCTGCGTGCGTACGAGTTGCTTTTCGAAGGTCGTTACAAGATAAACGACAACTCCACGAGCGGAAAGACGGTTACGCGTTTCACGTGATTTCACAGAgccg < 1:505595/204‑1 (MQ=255) ggCCAAAGCTTCCAGAACTTCACGCGGCGATGATTCACCGTGCAGCAACGGGCGTTTTTTATCACGCTGCGTGCGTGCAAGTTGCTTTTCGATGGTCGTTTCAAGATAAACGACACCGCCACGAGCGGAAAGACGGTTACGCGTTTCACGGGATTTCACAGAGCCGCCGCCAGTAGCCAGCACAATACCCTGTTTCTCGGTCAACTCATTGATGACCTTTTCTTCGCGATCGCGGAAGCCTTCTTCGCCtt > 1:186497/1‑251 (MQ=255) naGAACTTCTCGCGGCGATGTTTCACCGTGCAGCAACGAGCGTTTTTTATCACGCTGCGTGCGTGCACGTTGCTTTTCGATGGTCGTTTCAAGATAACCGACAACGCAACGCGCGGTAAGACGGTTACGCGTTTCACGGGATTTCACCGAGCCGCCGCCAGTCGCCAGCACACTACCCTGTTTCTCGTTCAACTCCTTGATGCACTTTTCTTCGCGATcgc > 1:105971/2‑221 (MQ=255) cAGAANTTNNCGNNGAGNTGTTTCAACGTGCAGCAACGGGCGTTTTTTATCACNc > 2:263924/1‑55 (MQ=38) cAGAACTTCACGCGGCGATGTTTCAACGTGCAGCAACGGGCGTTTTTTATCACGc < 1:263924/55‑1 (MQ=255) aGAACTNCNANCNGCGCTGTTTCAACGTGNAGCAACGGGCGTTTTTTAtnn > 2:471205/1‑49 (MQ=25) tcggcgaTGTTTCAACGTGCAGCAACGGGCGTTTTGTATCACGCTGCGTTCGTGCAAGTTGCTTTTCGGTGGTCGTTGCAAGATAAACGACAACGCCACGAGCGGAAAGACGGTTACGCTTTTCACGGGATTTCACAGAGCCGCCGCCATTAGCCAGCACAATACCCTGTTTCTCGGTCAACTCAGTTATGACCTTTTCTTCGCGATCGCGGAAGCCTTCTTCGCCTTCTAAATCGAAAACCCAg < 1:500309/244‑1 (MQ=255) | AATTGGGTAACTACGTTCCCCGAGAGTAACGACAATCCTCTCCATAACGCGACATCCACCTTAATTACTGTACCCGCAGACGAGTGTATATAAAGCCAGAATTAGTTGCTTTCCAGCATGTGAATAATCTGGTTTGCAACCACTTTAGCGCTTTGATCATCAGTACGAATGGTCACGTCGGCAATCTCTTCATACAGCGGATTGCGTTCATTGGCCAACGCTTCCAGAACTTCACGCGGCGGTGTTTCAACGTGCAGCAACGGGCGTTTTTTATCACGCTGCGTGCGTGCAAGTTGCTTTTCGATGGTCGTTTCAAGATAAACGACAACGCCACGAGCGGAAAGACGGTTACGCGTTTCACGGGATTTCACAGAGCCGCCGCCAGTAGCCAGCACAATACCCTGTTTCTCGGTCAACTCATTGATGACCTTTTCTTCGCGATCGCGGAAGCCTTCTTCGCCTTCTAAATCGAAAACCCAG > NC_000913/3518442‑3518921 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 24 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |