Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | CP009273 | 2,978,089 | A→C | 18.4% | M288R (ATG→AGG) | yqeF ← | short chain acyltransferase |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | CP009273 | 2,978,089 | 0 | A | C | 18.4% | 127.8 / 3.9 | 58 | M288R (ATG→AGG) | yqeF | short chain acyltransferase |
Reads supporting (aligned to +/- strand): ref base A (12/34); new base C (9/2); total (22/36) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.07e-03 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 |
GACCCGACGCTCATCCCACTCAAGCATCTTGCCAACCGAAAGCGCCTGTGCAGCAAACGCTTCATTAGCCTCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATACACCGGCG‑CAATTCCCATCAATGCCGGATCTACACCAACGCTGGCAAATGCGCGAATGCGGGCCAGCACGGGTAAATTCAACGCTCGTGCTTTGGCTTCGCTCATCATCATTACCGCAGCTGCGCCATCGTTTATGGATGATGCATTACCCGCTGTCACAGAAC > CP009273/2977941‑2978236 | gacCCGACGCTCATCCCACTCAAGCATCTTGCCAACCGAAAGCGCCTGTGCAGCAAACGCTTCATTAGCCTCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATACACCCGCC‑CCAATCCCAt > 1:475451/1‑150 (MQ=255) acCCGACGCTCATCCCACTCAAGCATCTTGCCAACCGAAAGCGCCTGTGCAGCAAACGCTTCATTAGCCTCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATACACCGGCG‑CAATTCCCAt < 2:792424/149‑1 (MQ=255) ccGACGCTCATCCCACTCAAGCATCTTGCCAACCGAAAGCGCCTGTGCAGCAAACGCTTCATTAGCCTCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATACACCGGCG‑CAATTCCCATCaa < 1:570818/150‑1 (MQ=255) ccGACGCTCATCCCACTCAAGCATCTTGCCAACCGAAAGCGCCTGTGCAGCAAACGCTTCATTAGCCTCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATACACCGGCG‑CAATTCCCATCaa < 1:414219/150‑1 (MQ=255) gCTCATCCCACTCAAGCATCTTGCCAACCGAAAGCGCCTGTGCAGCAAACGCTTCATTAGCCTCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATACACCGGCG‑CAATTCCCATCAATGCCg < 1:279965/150‑1 (MQ=255) aTCCCACTCAAGCATCTTGCCAACCGAAAGCGCCTGTGCAGCAAACGCTTCATTAGCCTCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATACACCGGCG‑CAATTCCCATCAATGCCGGAt < 1:420942/149‑1 (MQ=255) aTCCCACTCAAGCATCTTGCCAACCGAAAGCGCCTGTGCAGCAAACGCTTCATTAGCCTCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGGCGCCTACACCCGCG‑CAAATCCCATCAAATCCGGATc > 1:182490/1‑150 (MQ=255) ccACTCAAGCATCTTGCCAACCGAAAGCGCCTGTGCAGCAAACGCTTCATTAGCCTCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATACACCGGCG‑CAATTCCCATCAATGCCGGATCTac < 1:742396/150‑1 (MQ=255) cACTCAAGCATCTTGCCAACCGAAAGCGCCTGTGCAGCAAACGCTTCATTAGCCTCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATATAACGGCG‑CAATTTCCCTCAATGCCGGATCTaca > 1:548814/1‑150 (MQ=255) aaGCATCTTGCCAACCGAAAGCGCCTGTGCAGCAAACGCTTCATTAGCCTCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATACTTCGGAG‑TAATTCCGAGCCATGCGAGATCTACACCa < 2:43493/148‑1 (MQ=255) cATCTTGCCAACCGAAAGCGCCTGTGCAGCAAACGCTTCATTAGCCTCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATACACCGGCG‑CAATTCCCATCAATGCCGGATCTACACCAACGCt < 1:933839/150‑1 (MQ=255) cATCTTGCCAACCGAAAGCGCCTGTGCAGCAAACGCTTCATTAGCCTCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGGCGCATACACCCGCG‑CAATTTCCCTCCATGCCCGGTCTACACCCACCCt > 2:322159/1‑150 (MQ=255) gAAAGCGCCTGTGCAGCAAACGCTTCATTAGCCTCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATACACCCGCG‑CAATTCCCATCAATGCCGGAGCTACACCAACCCCGGCAAATGcgcg > 1:1033983/1‑149 (MQ=255) ccTGTGCAGGAAGCGCTTCATTAGCGTCGATAAGGTCGACTTCAGCCAACTGCCGGCCTACACGCTCCAGGCAACGGCGGGTCGCGTTCTCCTTAA‑TGATTCGGCGCACACCGAGATCTACACCAACGCTGGCAAATGCGCGAATGCggg < 2:123962/150‑1 (MQ=255) gtgCAGCAAACGCTTCATTAGCCTCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCCTACACCCCCG‑CAATTTCCATCCATGCCGGCTCTACACCAACCCCCGCAAATGCGCGAAAGCGGGc > 2:125864/1‑148 (MQ=255) tgCAGCAAACGCTTCATTAGCCTCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATACACCGGCG‑CAATTCCCATCAATGCCGGATCTACACCAACGCTGGCAAATGCGCGAATGCGGGCCa < 1:612791/149‑1 (MQ=255) gcaAACGCTTCATTAGCCTCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATACACCGGCG‑CAATTCCCATCAATGCCGGATCTACACCAACGCTGGCAAATGCGCGAATGCGGGCCAGCACg < 2:789637/150‑1 (MQ=255) aaaCGCTTCATTAGCCTCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATACACCGGCG‑CAATTCCCATCAATGCCGGATCTACACCAACGCTGGCAAATGcgcg < 1:558885/132‑1 (MQ=255) aaCGCTTCATTAGCCTCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATACACCGGCG‑CAATTCCCATCAATGCCGGATCTACACCAACGCTGGCAAATGCGCGAATGCGGGCCAGCACGGGt < 1:792153/150‑1 (MQ=255) aaCGCTTCATTAGCCTCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGGTGCAAACACCCGCG‑CAAATATCAATCATATCCGCTCTACACCAACGCTGGCCAATTCCCCAATGCGGGGCAGCACcgg > 2:428389/1‑149 (MQ=255) aCGCTTCATTAGCCTCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGGCGCCTATCCCCGCC‑CCATTCCCCTCAACTCCCGACCTACACCAAACCCGGCCAAAGCGCGCAAACGGGCCAGCCCgggg > 2:721788/1‑148 (MQ=255) ccTCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATACACCGGCG‑CAATTTCCCTCAAAGCCCGGTCTACACCAACGCTGGCAAATGCGCGAATGCGGGCCCGCCCCGGGAAATTCAACGccc > 1:238573/1‑148 (MQ=255) aTAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGGCGCCTACAACCCCGCCAATTCCC‑TCAACTCCGGATCTACACCCACCCCGGCCAAAGCGCGAAAGCGGGCCCGCACCGGGAAAATCAACACTcccgct > 2:522675/1‑145 (MQ=255) tAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATACACCGGCG‑CAATTCCCATCAATGCCGGATCTACACCAACGCTGGCAAATGCGCGAATGCGGGCCAGCACGGGTAAATTCAACGCTCGTGCtt < 1:287137/150‑1 (MQ=255) aaGATCGACTTCAGCCAGCTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATACACCCGCG‑CAAATCCCCTCAATGCCGGATCTACACCAACGCTGGCAAATGCGCGAATGCGGGCCAGCACGGGTAAATTCAACGCCCGTGCttt > 1:533404/1‑150 (MQ=255) gACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATACACCGGCG‑CAATTCCCATCAATGCCGGATCTACACCAACGCTGGCAAATGCGCGAATGCGGGCCAGCACGGGTAAATTCAACGCTCGTGCTTTGGCtt < 1:264082/149‑1 (MQ=255) gACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATACACCGGCG‑CAATTCCCATCAATGCCGGATCTACACCAACGCTGGCAAATGCGCGAATGCGGGCCAGCACGGGTAAATTCAACGCTCGTGCTTTGGCTTc < 1:199323/150‑1 (MQ=255) ttCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATACACCCGCG‑CCATTCCCCTCAATGCCGGGTCTACACCCACCCCGGCAAAAGCGCGAAAGCGGGCCAGCCCCGGGAAAATCAACGCCCCTGCCTTTGTTTCCCt > 2:422458/1‑150 (MQ=255) cAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATTCATCTGCG‑CAATTCCCAGCAATGCCAGATCTACACCAACGCTGGCAAATGCGCGAATGCGGGCCAGCACGGGTAAATTCAACGCTCGTGCTTTGGCTTCGCt < 2:571992/148‑1 (MQ=255) aGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATACACCGGCG‑CAATTCCCATCAATGCCGGATCTACACCAACGCTGGCAAATGCGCGAATGCGGGCCAGCACGGGTAAATTCAACGCTCGTGCTTTGGCTTCGCtcat < 1:740657/150‑1 (MQ=255) aGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATACACCCGCG‑CCATTCCCATCAATGCCGGATCTACACCCACCCTGGCAAAAGCGCGAATGCGGGCCAGCACCGGGAAAATCAAAGCTCGTGCTTTGGCTTCGCtcat > 1:607016/1‑150 (MQ=255) ccAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATACACCGGCG‑CAATTCCCATCAATGCCGGATCTACACCAACGCTGGCAAATGCGCGAATGCGGGCCAGCACGGGTAAATTCAACGCTCGTGCTTTGGCTTCGCtcatca < 1:17989/150‑1 (MQ=255) cAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCCTACACCCGCG‑CCAATCCCCCCAAAGCCCGGTCTACACCAAAGCCGGCCAAAACGCGAATGCGGGCCCGCCCCGGGAAATTCAACGCTCCTGCTTTTGCTTCGCtcatcat > 2:981634/1‑150 (MQ=255) aaCTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGGATACACCCGCG‑CCACTCCCCTCCATGCCGGGTCTACACCAACGCCGGCAAAAGCGCGGAAGCGGGCCCGCACCGGGGAAATCAACGCCCGTGCTTTTGCTTCTCtcatcct > 2:530809/1‑147 (MQ=255) aaCTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATACTTCGGCA‑CAATTCCCCTCAATGCCGGATCTACACCAACGCTGGCAAATGCGCGAATGCGGGCCAGCACGGGTAAATTCAACGCTCGTGCTTTGGCTTCGCtcatcatc < 2:475451/150‑1 (MQ=255) cTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATACACCCGCC‑CCATTTCCCTCAATATCCGGTCTACACCAACCCTGGCCAATATGCGAATGCGGGGCCGCACGGGGAAAATCAAAGCTCGTGCTTTTGCTTTGCTcctcatcat > 2:307685/1‑150 (MQ=255) ccAGCCTACACGCTCCAGGCAACGGCGGGTCGCATACACCGGCG‑CAATTCCCATCAATGCCGGATCTACACCAACGCTGGCAAATGCGCGAATGCGGGCCAGCACGGGTAAATTCAACGCTCGTGCTTTGGCTTCGCTCATCATCAtt < 1:89989/148‑1 (MQ=255) cAGCCTACACGCTCCAGGCAACGGCGGGTCGCATACACCGGCG‑CAATTCCCATCAATGCCGGATCTACACCAACGCTGGCAAATGCGCGAATGCGGGCCAGCACGGGTAAATTCAACGCTCGTGCTTTGGCTTCGCTCATCATCATTAcc < 1:37693/150‑1 (MQ=255) aGCCTACACGCTCCAGGCAACGGCGGGTCGCATACACCGGCG‑CAATTCCCATCAATGCCGGATCTACACCAACGCTGGCAAATGCGCGAATGCGGGCCAGCACGGGTAAATTCAACGCTCGTGCTTTGGCTTCGCTCATCATCATTACCg < 1:1028377/150‑1 (MQ=255) cacGCTCCAGGCAACGGCGGGTCGCATACACCGGCG‑CAATTCCCATCAATGCCGGATCTACACCAACGCTGGCAAATGCGCGAATGCGGGCCAGCACGGGTAAATTCAACGCTCGTGCTTTGGCTTCGCTCATCATCATTACCGCAGCTg < 1:322159/150‑1 (MQ=255) cGCTCCAGGCAACGGCGGGTCACATACACCGGCG‑CAATTCCCATCAATGCCGGATCTACACCAACGCTGGCAAATGCGCGAATGCGGGCCAGCACGGGTAAATTCAACGCTCGTGCTTTGGCTTCGCTCATCATCATTACCGCAGCTgcg < 1:766919/150‑1 (MQ=255) cAGGCAACGGCGGGTCGCATACACCGGCG‑CAATTCCCATCAATGCCGGATCTACACCAACGCTGGCAAATGCGCGAATGCGGGCCAGCACGGGTAAATTCAACGCTCGTGCTTTGGCTTCGCTCATCATCATTACCGCAGCTGCGCCATc < 2:254178/150‑1 (MQ=255) gcggGTCGCATACACCGGCG‑CAATTCCCATCAATGCCGGATCTACACCAACGCTGGCAAATGCGCGAATGCGGGCCAGCACGGGTAAATTCAACGCTCGTGCTTTGGCTTCGCTCATCATCATTACCGCAGCTGCGCCATCGTTTATGga > 2:770223/1‑150 (MQ=255) cggGTCGCATACTTCGGCG‑CAATCCGGATCAATGCGGGATCTACACCAACGCTGGCAAATGCGCGAATGCGGGCCAGCACGGGTAAATTCAACGCTCGTGCTTTGGCTTCGCTCATCATCATTACCGCAGCTGCGCCATCGTTTATGgat < 2:853155/150‑1 (MQ=255) gggTCGCATACACCGGCG‑CAATTCCCATCAATGCCGGATCTACACCAACGCTGGCAAATGCGCGAATGCGGGCCAGCACGGGTAAATTCAACGCTCGTGCTTTGGCTTCGCTCATCATCATTACCGCAGCTGCGCCATCGTTTATGgat > 1:354754/1‑149 (MQ=255) tCGCATACACCGGCG‑CAATTCCCATCAATGCCGGATCTACACCAACGCTGGCAAATGCGCGAATGCGGGCCAGCACGGGTAAATTCAACGCTCGTGCTTTGGCTTCGCTCATCATCATTACCGCAGCTGCGCCATCGTTTATGGatgatg < 1:307685/150‑1 (MQ=255) cATACACCGGCG‑CAATTCCCATCAATGCCGGATCTACACCAACGCTGGCAAATGCGCGAATGCGGGCCAGCACGGGTAAATTCAACGCTCGTGCTTTGGCTTCGCTCATCATCATTACCGCAGCTGCGCCATCGTTTATGGATGATGcgt < 1:478986/150‑3 (MQ=255) ccGGCG‑CAATTCCCATCAATGCCGGATCTACACCAACGCTGGCAAATGCGCGAATGCGGGCCAGCACGGGTAAATTCAACGCTCGTGCTTTGGCTTCGCTCATCATCATTACCGCAGCTGCGCCATCGTTTATGGATGATGCATTACCCg < 1:833867/150‑1 (MQ=255) cGGCG‑CAATTCCCATCAATGCCGGATCTACACCAACGCTGGCAAATGCGCGAATGCGGGCCAGCACGGGTAAATTCAACGCTCGTGCTTTGGCTTCGCTCATCATCATTACCGCAGCTGCGCCATCGTTTATGGATGATGCATTACCCGc > 1:393206/1‑150 (MQ=255) g‑cAATTCCCATCAATGCCGGATCTACACCAACGCTGGCAAATGCGCGAATGCGGGCCAGCACGGGTAAATTCAACGCTCGTGCTTTGGCTTCGCTCATCATCATTACCGCAGCTGCGCCATCGTTTATGGATGATGCATTACCCGCTGt > 2:603141/1‑149 (MQ=255) g‑cAATTCCCATCAATGCCGGATCTACACCAACGCTGGCAAATGCGCGAATGCGGGCCAGCACGGGTAAATTCAACGCTCGTGCTTTGGCTTCGCTCATCATCATTACCGCAGCTGCGCCATCGTTTATGGATGATGCATTACCCGCTGt > 2:115234/1‑149 (MQ=255) aaTTCCCATCAATGCCGGATCTACACCAACGCTGGCAAATGCGCGAATGCGGGCCAGCACGGGTAAATTCAACGCTCGTGCTTTGGCTTCGCTCATCATCATTACCGCAGCTGCGCCATCGTTTATGGATGATGCATTACCCGCTGTcac < 1:522675/150‑1 (MQ=255) aaTTCCCATCAATGCCGGATCTACACCAACGCTGGCAAATGCGCGAATGCGGGCCAGCACGGGTAAATTCAACGCTCGTGCTTTGGCTTCGCTCATCATCATTACCGCAGCTGCGCCATCGTTTATGGATGATGCATTACCCGCTGTca > 1:732816/1‑149 (MQ=255) aTTCCCATCAATGCCGGATCTACACCAACGCTGGCAAATGCGCGAATGCGGGCCAGCACGGGTAAATTCAACGCTCGTGCTTTGGCTTCGCTCATCATCATTACCGCAGCTGCGCCATCGTTTATGGATGATGCATTACCCGCTGTcaca < 1:13882/150‑1 (MQ=255) ttCCCATCAATGCCGGATCTACACCAACGCTGGCAAATGCGCGAATGCGGGCCAGCACGGGTAAATTCAACGCTCGTGCTTTGGCTTCGCTCATCATCATTACCGCAGCTGCGCCATCGTTTATGGATGATGCATTACCCGCTGTCACAg < 2:556389/150‑1 (MQ=255) atCCCATCAAAGAAGGAACTAAACAAACGCAGGCAAAAGCGCGAATGCGGGCCAGAAAGGGTAAAAACAACGAAAGAGATATGGATTCGCTAAAAAAAAATAAAGAAGCTGCGCCATCGTTTATGGATGATGCATTACCCGCTGTCACAg < 1:886382/149‑1 (MQ=255) ccATCAATGCCGGATCTACACCAACGCTGGCAAATGCGCGAATGCGGGCCAGCACGGGTAAATTCAACGCTCGTGCTTTGGCTTCGCTCATCATCATTACCGCAGCTGCGCCATCGTTTATGGATGATGCATTACCCGCTGTCACAGaa < 1:530809/149‑1 (MQ=255) ccATCAATGCCGGATCTACACCAACGCTGGCAAATGCGCGAATGCGGGCCAGCACGGGTAAATTCAACGCTCGTGCTTTGGCTTCGCTCATCATCATTACCGCAGCTGCGCCATCGTTTATGGATGATGCATTACCCGCTGTCACAGAAc < 1:125864/150‑1 (MQ=255) | GACCCGACGCTCATCCCACTCAAGCATCTTGCCAACCGAAAGCGCCTGTGCAGCAAACGCTTCATTAGCCTCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATACACCGGCG‑CAATTCCCATCAATGCCGGATCTACACCAACGCTGGCAAATGCGCGAATGCGGGCCAGCACGGGTAAATTCAACGCTCGTGCTTTGGCTTCGCTCATCATCATTACCGCAGCTGCGCCATCGTTTATGGATGATGCATTACCCGCTGTCACAGAAC > CP009273/2977941‑2978236 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |