Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | CP009273 | 602,762 | C→A | 12.2% | P26P (CCG→CCT) | ybdK ← | weak gamma‑glutamyl:cysteine ligase |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | CP009273 | 602,762 | 0 | C | A | 12.2% | 88.2 / 4.4 | 33 | P26P (CCG→CCT) | ybdK | weak gamma‑glutamyl:cysteine ligase |
Reads supporting (aligned to +/- strand): ref base C (15/14); new base A (2/2); total (17/16) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.86e-01 |
CGTTCGTTATCGCATACCTCCTGACGCTGCCATTTCTGAAACGGGTGCGTGCCACCGCCGCAAATTTCCAGATGATGGTCTGTGGCTGCCTGCAATACGACTTTCTGCATCGCTGAAAACTGCCCGGCAGCCTGGTTGATATCACGGCAAACATCCGTCGCCAGCTCCAGCATACTTTCGGTGATATCGTGCTTTACCTCTCCGGCCGTGATCTTATTTTTAACCGCGTCAATCAGCATTGAAGAGTCCTGGCTTAAGTCATAGCCCGGCGGATTAACCACCTGCATTTCCAGTTCAATACCGAGGGTAAAAGGTTCAGAAACATGAAAATCGGGTAATGGCATAGGTTTCTCTTATGTTGGCGTTTTCTATTCAGTATAGAAGTCGGAGCGGCTGGGCGAGATGCGGAAGTTCTGGAATGTTTCTTTTTTTTGGTGATGGTGAACTGATGGTGCCTGAAGCAATTTGGCTACTTTTGCAATGTGACAAGTTATGGCAGTGCTGACTGGTGGCGAAGAATTTTGACGATTGAGTCATG > CP009273/602496‑603033 | cgttcgttATCGCATACCGCCTGACGCTGCCATATCTGAAACGGGTGCGTGCCACCGCCGTAAATTTCCAGAGGATGGTCTGTGGCTGCCGGCACTACGACAATCTGCATCGCTGAAAATTGCCTGGCAGCCAGGTACATATCACGGCGAACATCCGTCGCAATCTCCAGCATACTGTAGTTGATCTCGGGCTTACCCTCACCAGCCGTCCACTAATTTTTAACCGCGACACTCAGCATTTAAGAGTCCTGGCATACGTCATANGCCGGCTGact > 2:71142/1‑273 (MQ=255) aTACCTCCTTACGCTGCCATTTCTGAAACGGGTGCGTGCCACCGCTGCAAATTTCCAGATGATGGTCTGTGGCTGCCTGCAATACGACTTTCTGCATCGCTGAAAACTGCCCGGCAGCCTGGTTGATATCACGGCAAACATCCGTCGCCAGCTCCAGCATACTTTCGGTGATATCGTGCTTTACCTCTCCGGCCGTGATCTTATTTTTAACCGCGTCAATCAGCATTGAAGAGTCCTGGCTTAAGTCATAGCCCGGCGGATTAACCACCTGCAtt < 1:3932/275‑1 (MQ=255) aCCTCCTGACTCTGTCTTTTCTGAAACGTGTGCGTGCCACCGCCGCAAATTTCCAGATGATGGTCTGTGGCTGCCTGCAATACGACTCTCTNCATCGCTGAAATCTGCCCGGCAGCCTGGTTGATATCACGGCAAACATCCGTCGCCAGCTCCAGCATACTTTCGGTGAGATCGTGCTTTACCTCTCCGGCCGTGATCTTATTTTTAACCGCGTCAATCAGCATTGAAGAGTCCTGGCTTAAGTCGTAGCACGCCGGATTAACCACCTGCATTTc < 1:510020/275‑1 (MQ=255) aCCTCCTGACGCGGCCATTTCTGAAACGGGTGCGTGCCACCGCCGCAAATTTCCAGATGATTGTCTGTGGCTGCCTGCAATACGTCTTTCTGCATCGCTGAAAACTGCCCGGCAGCCTGGTTGATATCACGGCAAACATCCGTCGCCAGCTCCAGCATACTTTCGGTGATATCGTGCTTTACCTCTCCGGCCGTGATCTTATTTTTAACCGCGTCAATCAGCATTGAAGAGTCCTGGCTTAAGTCATAGCCCGGCGGATTAACCACCTGCATTTc < 1:233802/275‑1 (MQ=255) ctGACGCTGCCATTTCTGAAACGGGTGCGTGCCACCGCCGCAAATTTCCAGATGATGGTCTGTGGCTGCCTGCAATACGACTTTCTGCATCGCTGAAAACTGCCCGGCAGCCTGGTTGATATCACGGCAAACATCCGTCGCCAGCTCCAGCATACTTTCGGTGATATCGTGCTTTACCTCTCCGGCCGTGATCTTATTTTTAACAGCGTCAATCAGCATTGAAGAGTCCTGGCTTAAGTCCTAGCCCGGCAGAATAACCCCCTGCATTTCCAggt > 2:233800/1‑273 (MQ=255) gttgGTCTGTGGCTGCCTGCAATACGACTTTCTGCATCGCTGAAAACTGCCCGGCAGCCTGGTTGATATCACGGCAAACTTCCGTCGCCAGCTCCAGCATACTTTCGGTGATATCGTGCTTTACCTCTCCGGCCGTGATCTTATTTTTAACCGCGTCAATCAGCATTGAAGAGTCCTGGCTTAAGTCATAGCCCGGCGGATTAACCACCTGCATTTCCAGTTCAATACCGAGGGTAAAAGGTTCAGAAACATGAAAATCGGGTAATGGCATAGGt < 1:735536/273‑1 (MQ=255) tgGTCTGTGGCTGCCTGCAATACGACTTTCTGCATCGCTGAAAACTGACCGGCAGCCTGGTTGATATCACGGCAAACATCCGTCGCCAGCTCCAGCATACTTTCGGTGATATCGTGCTTTACCTCTCCGGCCGTGATCTTATTTTTAACCGCGTCAATCAGCATTGAAGAGTCCTGGCTTAAGTCATAGCCAGGCAGATTAACCACCTGCATTTCCAGTTCAATACCGAGGGTAAAAGGTTCAGAAACATGAAAATCGGGAAATGGCATAGGttt > 1:687076/1‑275 (MQ=255) ctgcAATACTACTTTCTGCATCGCTGAAAACTGCCCGGCAGCCTGGTTGATATCACGGCAAACATCCGTCGCCAGCTCGAGCATACTTTCGGTGATATCGTGCTTTACCTCTCCGGCCGTGATCTTATTTTTAACCGCGTCAATCAGCATTGAAGAGTCCTGGCTTAAGTCATAGCCCGGCGGATTAACCACCTGCATTTCCAGTTCAATACCGAGGGTAAAAGGTTCAGAAACATGAAAATCGGGTAATGGCATAGGTTTCTCTTATGTTGGc < 1:731002/274‑1 (MQ=255) gaaaaCGCCTTTCTGCATCGCTGAAAACGGCCCGGCAGCTTGGTTGATAGCACGGCAAACATCCGTCGCCAGCGCCAGCATACTTTCGGTGATATCGTGCTTTACCTCGCCGGCCGTGATCTTATTTTTAACCGCGTCAATCAGCATTGAAGAGTCCTGGCTTAAGTCATAGCCAGGCGGATTAACCACCTGCATTTCCAGTTCAATACCGAGGGTAAAAGGTTCAGAAACATGAAAATCGGGTAATGGCATAGGTTTCTCTTATGTTGGCGtt < 2:687062/270‑1 (MQ=255) atCACGGCAAACATCCGTCGCCAGCTCCAGCATACTTTCGGTGATATCGTGCTTTACCTCTCCGGCCGTGATCTTATTTTTAACCGCGTCAATCAGCATTGAAGAGTCCTGGCTTAAGTCATAGCCCGGCGGATTAac > 1:102102/1‑138 (MQ=255) atCACGGCAAACATCCGTCGCCAGCTCCAGCATACTTTCGGTGATATCGTGCTTTACCTCTCCGGCCGTGATCTTATTTTTAACCGCGTCAATCAGCATTGAAGAGTCCTGGCTTAAGTCATAGCCCGGCGGATTAACCACCTGCATTTCCANTTCACTACCGAGGGTAAAAGGTTCAGAAACATGAAAATCGGGTAATGGCATAGGTTTCTCTTATGTTGGCg > 1:615003/1‑224 (MQ=255) cGGCAAACATCCGTCGCCAGCTCCAGCATACTTTCGGTGATATCGTGCTTTACCTCTCCGGCCGTGATCTTATTTTTAACCGCGTCAATCAGCATTGAAGAGTCCTGGCTTAAGTCATAGCCCGGCGGATTAaccacc < 2:102101/138‑1 (MQ=255) cGGCAAACATCCGTCGCCAGCTCCAGCATACTTTCGGTGATATCGTGCTTTACCTCTCCGGCCGTGATCTTATTTTTAACCGCGTCAATCAGCATTGAAGAGTCCTGGCTTAAGTCATAGCCCGGCGGATTAACCACCTGCAtt > 2:718049/1‑144 (MQ=255) ggCAAACATCCGTCGCCAGCTCCAGCATACTTTCGGTGATATCGTGCTTTACCTCTCCGGCCGTGATCTTATTTTTAACCGCGTCAATCAGCATTGAAGAGTCCTGGCTTAAGTCATAGCCCGGCGGATTAac > 2:564330/1‑133 (MQ=255) aaaCATCCGTCGCCAGCTCCAGCATACTTTCGGTGATATCGTGCTTTACCTCTCCGGCCGTGATCTTATTTTTAACCGCGTCAATCAGCATTGAAGAGTCCTGGCTTAAGTCATAGCCCGGCGGATTAACCACTTGCATTTCCa < 1:718063/144‑1 (MQ=255) aaCATCCGTCGCCAGCTCCAGCATACTTTCGGTGATATCGTGCTTTACCTCTCCGGCCGTGATCTTATTTTTAACCGCGTCAATCAGCATTGAAGAGTCCTGGCTTAAGTCATAGCCCGGCGGATTAaccacc < 1:564342/133‑1 (MQ=255) ccTCTCCGGCCGTGATCTTATTTTTAACCGCGTCAATCAGCATTGAAGAGTCCTGGCTTAAGTCATAGCCCGGCGGATTAACCACCTGCATTTCCAGTTCAATACCGAGGGTAAACGGTTCAGAAACATGAAAATCGGGTAAt > 2:589499/1‑143 (MQ=255) tcCGGCCGTGATCTTATTTTTAACCGCGTCAATCAGCATTGAAGAGTCCTGGCTTAAGTCATAGCCCGGCGGATTAACCACCTGCATTTCCAGTTCAATACCGAGGGTAAAAGGTTCAGAAACATGAAAATCGGGTAATGGCa < 1:589512/143‑1 (MQ=255) cGGCCGTGATCTTATTTTTAACCGCGTCAATCAGCATTGAAGAGTCCTGGCTTAAGTCATAGCCCGGCGGATTAACAACCTGCATTTCAAGTCCAAAACCGAGGGTAAAAGGTTCGGAAACATGAAAATCGGGTACGGGAATAGGTTGTTCTTAGGTTGGCGTTTTATATTCAGAATAGAAGTCGGAGCGGATGGGGGAGAGGAGGAGGTACGGGAAGGATCCTTGTTTATGGGGCTGGTGCACTGGTGTGGCGTGAAGACATTTTTGTACtttt > 2:364619/1‑275 (MQ=255) aTTTTTAACCGCGTCAATCAGCATTGAAGAGTCCTGGCTTAAGTCATAGCCCGGCGGATTAACCACCTGCATTTCCAGTTCAATACCGAGGGTAAAAGGTTCAGAAACATGAAAATCGGGTAATGGCATAGGTTTCTCTTATGTTGACGTTTTCTATTNAGTATAGAAGTCGGAGCGGCTGGGCGAGATGCGGAAGTTCTGGAATGTTTCTTTTTTTTGGTGATGGTGAACTGATGGTGCCTGAAGCAATTTGGCTACTTTTGCAATGTGACAAg > 1:663811/1‑275 (MQ=255) tGAAGAGTCCTGGCTTAAGTCATAGCCAGGCGGATTAACCACCTGCATTTCCAGTTCAATACCGAGGGTAAAAGGTTCAGAAACATGAAAATCGGGTAATGGCATAGGTTTCTCTTATGTTGGCGTTTTCTATTCAGTATAGAAGTCGGAGCGGCTGGGCGAGATGCGGAAGTTCTGGAATGTTTCTTTTTTTTGGTGATGGTGAACTGATgg > 1:690215/1‑213 (MQ=255) gagTCCTGGCTTAAGTCATAGCCAGGCGGATTAACCACCTGCATTTCCAGTTCAATACCGAGGGTAAAAGGTTCAGAAACATGAAAATCGGGTAATGGCATAGGTTTCTCTTATGTTGGCGTTTTCTATTCAGTATAGAAGTCGGAGCGGCTGGGCGAGATGCGGAAGTTCTGGAATGTTTCTTTTTTTTGGTGATGGTGAACTGATGGTGcc < 2:690201/213‑1 (MQ=255) tCCTGGCTTAAGTCATAGCCCGGCGGATTAACCACCTGCATTTCCAGTTCAATACCGAGGGTAAAAGGTTCAGAAACATGAAAATCGGGTAATGGCATAGGTTTCTCTTATGTTGGCGTTTTc > 2:538686/1‑123 (MQ=255) cTGGCTTAAGTCATAGCCCGGCGGATTAACCACCTGCATTTCCAGTTCAATACCGAGGGTAAAAGGTTCAGAAACATGAAAATCGGGTAATGGCATAGGTTTCTCTTATGt > 1:352379/1‑111 (MQ=255) ggCTTAAGTCATAGCCCGGCGGATTAACCACCTGCATTTCCAGTTCAATACCGAGGGTAAAAGGTTCAGAAACATGAAAATCGGGTAATGGCATAGGTTTCTCTTATGTTGGCGTTTTCTAtt < 1:538696/123‑1 (MQ=255) gCTTAAGTCATAGCCCGGCGGATTAACCACCTGCATTTCCAGTTCAATACCGAGGGTAAAAGGTTCAGAAACATGAAAATCGGGTAATGGCATAGGTTTCTCTTATGTTGGCGTTTTCTATTCAGTATAGAAGTCGGAGCGGCTGGGCGAGATGCGGAAGTTCTGGAATGTTTCTTTTTTTTGTTGATGGTGAAc > 2:147567/1‑195 (MQ=255) cTTAAGTCATAGCCCGGCGGATTAACCACCTGCATTTCCAGTTCAATACCGAGGGTAAAAGGTTCAGAAACATGAAAATCGGGTAATGGCATAGGTTTCTCTTATGTTGGc < 2:352373/111‑1 (MQ=255) aaGTCATAGCCCGGCGGATTAACCACCTGCATTTCCAGTTCAATACCGAGGGTAAAAGGTTCAGAAACATGAAAATCGGGTAATGGCATAGGTTTCTCTTATGTTGGCGTTTTCTATTCAGTATAGAAGTCGGAGCGGCTGGGCGAGATGCGGAAGTTCTGGAATGTTTCTTTTTTTTGTTGATGGTGAACTGAt < 1:147568/195‑1 (MQ=255) ggTCATAGCCCGGCGGATTAACCACCTTCATTTCCAGTTCAATACCGAGGGTAAAAGGTTCAGAAACATGAAAATCGGGTAATGGCATAGGTTTCTCTTATGTTGGCGTTTTCTATTCAGTATAGAAGTCGGAGCGGCTGGGCGAGATGCGGAAGTTCTGGAATGTTTCTTTTTTTTGGTGATGGTGAACTGATGGTGCCTGAAGCAATTTGGCTACTTTTGCAATGTGACAAGTTATGGCAGTGCTGACTGGTGGCGAAGAATTTTGACGAtt < 1:463325/273‑1 (MQ=255) aGTCATAGCCCGGCGGATTAACCACCTGCATTTCCAGTTCAATACCGAGGGTAAAAGGTTCAGAAACATGAAAATCGGGTAATGGCATAGGTTTCTCTTATGTTGGCGTTTTCTATTCAGTATAGAAGTCGGAGCGGCTGGGCGAGATGCGGAAGTTCTGGAATGTTTCTTTTTTTTGGTGATGGTGAACTGATGGTGCCTGAAGCAATTTGGCTACTTTTGCAATGTGAAAAGTTATGGCAGTGCTGAAGGGTGGGGAAGAATTTTGACGaggt > 2:202405/1‑272 (MQ=255) tAGCCCGGCGGATTAACCACCTGCATTTCCAGTTCAATACCGAGGGTAAAAGGTTCAGAAACATGa > 2:175505/1‑66 (MQ=255) tAGCCCGGCGGATTAACCACCTGCATTTCCAGTTCAATACCGAGGGTAAAAGGTTCAGAAACATGAAAATCGGGTAATGGCATAGGTTTCTCTTATGTTGGCGTTTTCTATTCAGTATAGAAGTCGGAg > 2:652356/1‑129 (MQ=255) gCCCGTCGGATTAACCTCCTGCATTTCCAGTTCAATACCGAGGGTAAAAGTTTCAGAATCATGAAAATCGGGTAATGGCATAGGTTTCTCTTATGTTGGCGTTTTCTATTCAGTATAGAAGTCGGAGCGGCTGGGCGAGATGCGGAAGTTCTGGAATGTTTCTTTTTTTTGGTGATGGTGAACTGATGGTGCCTGAAGCAATTTGGCTACTTTTGCAATGTGACAAGTTATGGCAGTGCTGACTGGTGGCGAAGAATTTTGACGATTGAGTCATg < 1:202407/275‑1 (MQ=255) ccGGCGGATTAACCACCTGCATTTCCAGTTCAATACCGAGGGTAAAAGGTTCAGAAACATGAAAAt < 1:175507/66‑1 (MQ=255) | CGTTCGTTATCGCATACCTCCTGACGCTGCCATTTCTGAAACGGGTGCGTGCCACCGCCGCAAATTTCCAGATGATGGTCTGTGGCTGCCTGCAATACGACTTTCTGCATCGCTGAAAACTGCCCGGCAGCCTGGTTGATATCACGGCAAACATCCGTCGCCAGCTCCAGCATACTTTCGGTGATATCGTGCTTTACCTCTCCGGCCGTGATCTTATTTTTAACCGCGTCAATCAGCATTGAAGAGTCCTGGCTTAAGTCATAGCCCGGCGGATTAACCACCTGCATTTCCAGTTCAATACCGAGGGTAAAAGGTTCAGAAACATGAAAATCGGGTAATGGCATAGGTTTCTCTTATGTTGGCGTTTTCTATTCAGTATAGAAGTCGGAGCGGCTGGGCGAGATGCGGAAGTTCTGGAATGTTTCTTTTTTTTGGTGATGGTGAACTGATGGTGCCTGAAGCAATTTGGCTACTTTTGCAATGTGACAAGTTATGGCAGTGCTGACTGGTGGCGAAGAATTTTGACGATTGAGTCATG > CP009273/602496‑603033 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 10 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |