Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_000913 | 2,041,159 | G→A | 8.7% | intergenic (+41/‑216) | msrQ → / → zinT | periplasmic protein‑L‑methionine sulfoxide reductase heme binding subunit/metal‑binding protein ZinT |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 2,041,159 | 0 | G | A | 8.7% | 137.7 / 3.3 | 54 | intergenic (+41/‑216) | msrQ/zinT | periplasmic protein‑L‑methionine sulfoxide reductase heme binding subunit/metal‑binding protein ZinT |
Reads supporting (aligned to +/- strand): ref base G (31/18); new base A (2/3); total (33/21) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 3.66e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.22e-01 |
TGGCTGTACTGCTTTTAGCCTTACGGTATAAGAAGTTGCGTTCCCTTTTTAACCGGTTACGCAAGCAGGTTCACAACAAATTATCTGTATAACTTCATATTGCCGACAAAGTACGCGTTATCCCTGCCGGATGTGGCGCGAGCGCCCTATCCGGCATACACAGATCGCATCAGGTTATGTTTCGTTTGTCATCAACCTCAAATCACCAAAATAATTTATTTTATTTCAATAAGATAGATAAGTAGAACTGAGAAAGCCATGCTCTCGTTTCCTAA > NC_000913/2041027‑2041301 | tggctgTACTGCTTTTAGCCTTACGGTATAAGAAGTTGCGTTCCCTTTTTAACCGGTTACGCAAGCAGGTTCACAACAAATTATCTGTATAACTTCATATTGCCGACAAAGTACGCGTTATCCCTGCCGGATGTGGCGCGAgc < 1:151686/143‑1 (MQ=255) gctgTACTGCTTTTAGCCTTACGGTATAAGAAGTTGCGTTCCCTTTTTAACCGGTTACGCAAGCAGGTTCACAACAAATTATCTGTATAACTTCATATTGCCGACAAAGTACGCGTTATCCCTGCCGTATGTGGCGCGAgcgc > 2:376367/1‑143 (MQ=255) gctgTACTGCTTTTAGCCTTACGGTATAAGAAGTTGCGTTCCCTTTTTAACCGGTTACGCAAGCAGGTTCACAACAAACTATCTGAACAACTTCCTAGTGCCGACCAAATACGCGTTATCCCTGCCGCATGCGGCGCGaaccc > 2:53968/1‑139 (MQ=255) gTACTGCTTTTAGCCTTACGGTATAAGAAGTTGCGTTCCCTTTTTAACCGGTTACGCAAGCAGGTTCACAACAAATTATCTGTATAACTTCATATTGCCGACAAAGTACGCGTTATCCCTGCCGGATGTGGCGCGAGCGCCCt > 1:14233/1‑143 (MQ=255) gCTTTTAGCCTTACGGTATAAGAAGTTGCGTTCCCTTTTTAACCGGTTACGCAAGCAGGTTCACAACAAATTATCTGTATAACTTCATATTGCCGACAAAGTACGCGTTATCCCTGCCGGATGTGGCGCGAGCGCCCTATCCg > 1:343370/1‑143 (MQ=255) cTTTTAGCCTTACGGTATAGTAAGTTGCGTTCCCTTTTTAACCGGTTACGCAAGCAGGTTCACAACAAATTATCTGTATAACTTCATATTGCCGACAAAGTACGCGTTATCCCTGCCGTATGTGGCGCGAGCGCCCTATCCgg < 1:448826/143‑1 (MQ=255) tAGTCTTTCGGTATAAGAAGTTGCGTTCTCTTTTTACTCGGTTACGCAAGCGGGTGCACAACAAATTATCTGTATAACTTCATATTGCCGACAAAGTCCGCGTTATCCCTGCCGGATGTGGTGCGAGCCCCCTATCCGGCATa < 2:521021/143‑1 (MQ=255) tAGCCTTACGGTATAAGAAGTTGCGTGCCCTTTTTAACCGGTTACGCAAGCAGGTTCACAACAAATTATCTGTATAACTTCATATTGCCGACAAAGTACGCGTTATCCCTGCCGGATAGGGCGCTCGCGCCACATCCGGCagg < 2:723894/143‑3 (MQ=38) gaGCCTTACGGTATAAGAAGTTGCGTTCCCTTTTTAACCGGAAGCGTATGAGCGTGTACCACCAACAATCTGTATAACTTCATATTGCCGACAAAGTACGCGTTATCCCTGCCGGATAGGGCGCTCGCGCCACATCCGGCagg < 2:363642/142‑3 (MQ=16) agtcTTACTGTATAAGAAGTTGCGTTCCCTTTTTAACCGGTCAGGCAAGCAGGTTCACAACAAATTATCTGTATAACTTCATATTGCCGACAAAGTACGCGTTATCCCTGCCGGATGTGGCGCGAGCGCCCTATCCGGCATac < 1:265402/140‑1 (MQ=255) ccTTACGGTATAGCAAGTTGCGTTCCCTGTTTAACCGGTTACGCAAGCAGGTTCACAACAAATTATCTGTATAACTTCATATTGCCGACAAAGTACGCGTTATCCCTGCCGGATAGGGCGCTCGCGCCACATCCGGCagggat < 2:199765/143‑6 (MQ=37) cgtACGGTATAAGAAGTTGCGTTCCCTTGTTAACCGGTTACGCAAGCAGGTTCACAACAAATTATCTGTATCACTTCATATTGCCGACAAAGTACGCGTTATCCCTGCCGGATGTGGCGCGAGCGCCCTATCCGGCATAcaca < 2:529987/141‑1 (MQ=255) aagaagTTGCGTTCCCTTTTTAACCGGTTACGCAAGCAGGTTCACAACAAATTATCTGTATAACTTCATATTGCCGACAAAGTACGCGTTATCCCTGCCGGATGTGGCGCGAGCGCCCTATCCGGCATACACAGATCGCCTCa > 1:42969/1‑143 (MQ=255) tggtgTTGTGTGCGCTTGTTATCCGGTTACGAAAGCAGGTTCACAACAAATGAGCTGTATAACTTCAGATTGCCGACAAAGTACGCGTTAGCCCTGCCGGATGTTGCGCGAGCGCCCTATCCGGCATACACAGATCGCATCAg < 2:82345/139‑1 (MQ=255) gaagTTGCGTTCCCTTTTTAACCGGTTACGCAAGCAGGTTCACAACACATTATCTGTATAACTTCATATTGCCGACAAAGTACGCGTTATCCCTGCCGGATGTGTCGCGAGCGCCCTATCCGGCATACACAGATCGCCTCAgg > 2:304986/1‑143 (MQ=255) gaagTTGCGTTCCCTTTTTAACCGGTTACGCAAGCAGGTTCACAACAAATTATCTGTATAACTTCATATTGCCGACAAAGTACGCGTTATCCCTGCCGGATGTGGCGCGAGCGCCCTATCCGGCATACACAGATCGCATCAgg > 1:546988/1‑143 (MQ=255) gTTGCGTTCCCTTTTTAACCGGTTACGCAAGCAGGTTCACAACAAATTATCTGTATAACTTCATATTGCCGACAAAGTACGCGTTATCCCTGCCGGATGTGGCGCGAGCGCCCTATCCGGCATACACAGATCGCATCAGGTTa > 1:477614/1‑143 (MQ=255) ttCCCTTTTTAACCGGTTACGCAAGCAGGTTCACAACAAATTATCTGTATAACTTCATATTGCCGACAAAGTACGCGTTATCCCTGCCGGATGTGGCGCGAGCGCCCTCTCCGGCATCCACCCCTCGCCTCATGTCATATTTc > 1:22211/1‑143 (MQ=255) cccTTTTTAACCGGTTACGCAAGCAGGTTCACAACAAATTATCTGTATAACTTCATATTGCCGACAAAGTACGCGTTATCCCTGCCGGATGTGGCGCGAGCGCCCTATCCGGCATACACAGATCGCATCAGGTTATGTTTCGt < 1:14781/143‑1 (MQ=255) ccTTTTTAACCGGTTACGCAGGCAGGTTCACAACAAATTATCTGTATAACTTCATATTGCCGACAAAGTACGCGTTATCCCTGCCGGATGTGGCGCGAGCGCCCTATCCGGCATACACAGATCGCATCAGGTTATGTTTCGtt < 2:299899/143‑1 (MQ=255) ccTTTTTAACCGGTTACGCAAGCAGGTTCACAACAAATTCTCTGTATAACTTCATATTGCCGACAGAGTACGCGTTATCCCTGCCGGATAATGCGCGAGCCCCCCACCCCGCATACTCAGCTAGAATCAGGTGCTGTCTCgat > 2:724171/1‑141 (MQ=255) cTTTTTAACCGGTTACGCAAGCAGGTTCACAACAAATTATCTGTATAACTTCATATTGCCGACAAAGTACGCGTTATCCCTGCCGGATGTGGCGCGAGCGCCCTATCCGGCATACACAGATCGCATCAGGTTATGTTTCGttt > 1:355334/1‑143 (MQ=255) cTTTTTAACCGGTTACGCAAGCAGGTTCACAACAAATTATCTGTATAACTTCATATTGCCGACAAAGTACGCGTTATCCCTGCCGGATGTGGCGCGAGCGCCCTATCCGGCATACACAGATCGCATCAGGTTATGTTTCGtct > 2:735174/1‑141 (MQ=255) cTTGTTAACCGGTTACGCAAGCACGTTCACAACAAATTATCTGTATAACTTCATATTGCCGACAAAGTACGCGTTATCCCTGCCGGATGTGGCGCGAGCGCCCTATCCGGCATACACAGATCGCATCAGGTTATGTTTCGttt < 1:580884/143‑1 (MQ=255) ttttctccGGTTACGCAAGCAGGTTGACAACAAATTATCTGTATACCTTCCTATTGCCGACAAAGTACGCGTTACTCCTGCCGGTTGTGGCGCGAGCGCCCTATCCGGCATACACAGATCGCATCAGGTTATGTTTCGTTTGt < 2:411246/137‑1 (MQ=255) tttAACCGGTTACGCAAGCAGGTTCACAACAAATTATCTGTATAACTTCATATTGCCGACAAAGTACGCGTTATCCCTGCCGGATGTGGCGCGAGCGCCCTATCCGGCATACACCGATCGCATCAGGTTCTGTTTCGTTTGtt > 1:369001/1‑142 (MQ=255) aCCGGTTACGCAAGCAGGTTCACAACAAATTATCTGTATAACTTCATATTGCCGACAAAGTACGCGTTATCCCTGCCGGATGTGGCGCGAGCGCCCTATCCGGCATACACAGATCGCATCAGGTTATGTTTCGTTTGtcatca > 2:186841/1‑143 (MQ=255) ggTTACGCAAGCAGGTTCACAACAAATTATCTGTATAACTTCATATTGCCGACAAAGTACGCGTTATCCCTGCCAGATGTGGCGCGAGCGCCCTATCCGGCATACACAGATCGCATCAGGTTATGTTTCGTTTGTCATCAAcc > 2:96176/1‑143 (MQ=255) gTTACGCAAGCAGGTTCACAACAAATTATCTGTATAACTTCATATTGCCGACAAAGTACGCGTTATCCCTGCCGGATGTGGCGCGAGCGCCCTATCCGGCATACACAGATCGCATCAGGTTATGTTTCGTTTGTCATCAAcgt > 1:326828/1‑141 (MQ=255) gTTACGCAAGCAGGTTCACAACAAATTATCTGTATAACTTCATATTGCCGACAAAGTACGCGTTATCCCTGCCGGATGTGGCGCGAGCGCCCTATCCGGCATACACAGATCGCATCAGGTTATGTTTCGTTTGTCATCAAACt > 2:427485/1‑143 (MQ=255) ttACGCAAGCAGGTTCACAACAAATTATCTGTATAACTTCATATTGCCGACAAAGTACGCGTTATCCCTGCCGGATGTGGCGCGAGCGCCCTATCCGGCATCCCCAGGTCGCCTCGGGTTCTGCTTCGTTTTtcaccatgctg > 1:703687/1‑138 (MQ=255) aCGCAAGCAGGTGCACAACAAATTATCCGTATGACTTCATACTGCCGACGAAGTACGCGTTATCCCTGCCGGATGTGGCGCGAGCGCCCTATCCGGCATACACAGATCGCATCAGGTTATGTTTCGTTTGTCATCAACCTCaa < 2:17997/143‑1 (MQ=255) gCAAGCAGGTTCACAACAAATTATCTGTATAACTTCATATTGCCGACAAAGTACGCGTTATCCCTGCCGGATGTGGCGCGAGCGCCCTATCCGGCCTACCCAGATCGCCTCCGGCTATGTTTCGTCTGTCATCAACCTCAAAt > 1:52270/1‑143 (MQ=255) gCAAGCAGGTTCACAACAAATTATCTGTATAACTTCATATTGCCGACAAAGTACGCGTTATCCCTGCCGGATGTGGCGCGAGCGCCCTATCCGGCATACACAGATCGCATCAGGTTATGTTTCGTGTGTCACCAACCTCaaag > 2:337607/1‑142 (MQ=255) aaccaATTATCTGTATAACTTCATATTGCCGACAAAGTACGAGTTATCCCTGCCGGATGTGGCGCGAGCGCCCTATCCGGCATACACAGATCGCATCAGGTTATGTTTCGTTTGTCATCAACCTCAAATCACCAAAATAattt > 2:720263/1‑143 (MQ=255) acaaTTTATCTGTATAACTTCATATTGCCGACAAAGTACGCGTTATCCCTGCCGGATGTGGCGCGAGCGCCCTCTCCGGCATACACAGCTCGCCTCCGGTTATGTTTCGTTTGTTATCAACCTCAAATGACAAAAATAattaa > 1:277230/1‑141 (MQ=255) aTCGGTATAACTTCATATTGCGGACAAAGTCCGCGTTATCCCTGCCGGATGTGGCGCGAGCGCCCTATCCGGCATACACAGATCGCATCAGGTTATGTTTCGTTTGTCATCAACCTCAAATCACCAAAATAATttattttatt < 2:746242/143‑1 (MQ=255) ttCATATTGCCGACAAAGTACGCGTTATCCCTGCCGGATGTGGCGCGAGCGCCCTATCCGGCATACACAGATCGCATCAGGTTATGTTTCGTTTGTCATCAACCTCAAATCACCAAAATAATTTATTTTATTTCAATAagata > 1:20494/1‑143 (MQ=255) ttCACATTGCCGACAAAGTACGCGTTATCCCGGCCGGATGTGGCGCGAGCGCCCTATCCGGCAGACACAGATCGCATCAGGTTATGTTTCGTTTGTCATCAACCTCAAATCACCAAAATAATTTATTTTATTTCAATAagata < 1:644295/143‑1 (MQ=255) atatTGCCGACAAAGTACGCGTTATCCCTGCCGGATGTGGCGCGAGCGCCCTATCCGGCATACACAGATCGCATCAGGTTATGTTTCGTTTGTCATCAACCTCAAATCACCAAAATAATTTATTTTATTTCAATAagatagat < 1:810843/143‑1 (MQ=255) tatTGCCGACAAAGTACGCGTTATCCCTGCCGGATGTGGCGCGAGCGCCCTATCCGGCATACACAGATCGCATCAGGTTATGTTTCGTTTGTCATCAACCTCAAATCACCAAAATAATTTATTTTATTTGAATAAgatagata > 2:760199/1‑143 (MQ=255) tGCCGACAAAGTACGCGTTATCCCTGCCGGATGTGGCGCGAGCGCCCTATCCGGCATACACAGATCGCATCAGGTTATGTTTCGTTTGTCATCAACCTCAAATCACCAAAATAATTTATTTTATTCCAATAAGATAGATAAGt < 1:389001/143‑1 (MQ=255) gACAAAGTACGCGTTATCCCTGCCGGATGTGGCGCGAGCGCCCTATCCGGCATACACAGATCGCATCAGGGTATGTTTCGTTTGTCATCAACCTCAAATCACCAACATAATTTATTTTATTTCAATAAGATAGATAAGTAGaa > 2:256183/1‑143 (MQ=255) gTACGCGTTATCCCTGCCGGATGTGGCGCGAGCGCCCTATCCGGCATACACAGATCGCATCAGGTTATGTTTCGTTTGTCATCAACCTCAAATCACCAAAATAATTTATTTTATTTCAATAAGATAGATAAGTAGAACTgaga > 1:690181/1‑143 (MQ=255) aCGCGTTATCCCTGCCGGATGTGGCGCGAGCGCCCTATCCGGCATACACAGATCGCATCAGGTTATGTTTCGTTTGTCATCAACCTCAAATCACCAAAATAATTTATTTTATTTCAATAAGATAGATAAGTAGAACTGAGaaa < 2:20494/143‑1 (MQ=255) cgTTATCCCTGCCGGATGTTGCGCGAGCACCCTATCCGGCATCCACCGGTCGCATCAGGGTATGTTTCGTCTGTCATCCACCCCCACTCACCACAACAATTCATTCTATTTCAGTCATATAGATAACTACAACTGATAAAGct > 2:523911/1‑142 (MQ=255) tATCCCTGCCGGATGTGGCGCGAGCGCCCTATCCGGCATACACAGATCGCATCAGGTTATGTTTCGTTTGTCATCAACCTCAAATCACCAAAATAATTTATTTTATTTCAATAAGATAGATAAGTAGAACTGAGAAAGCCATg < 2:343370/143‑1 (MQ=255) ccTGCCGGATGTGGCGCGAGCGCCCTGTCCGGCATACACAGATCGCATCAGGTTATGTTTCGTTTGTCAGCAACCTCAAATCACCAAAATAATTTATTTTATTTCAATAAGATAGATAAGTAGAACTGCGAAAGCCATGctat > 2:678511/1‑141 (MQ=255) ccTGCCGGATGTGGCGCGAGCGCCCTATCCGGCATACACATATCGCATCAGGTTATGTTTCGTTTGTCTTCAAACTAAAATCAACAAAATAATTTATTTTATTTCAATAAGATAGATAAGTAAAACTGAGAAAGCCATGctct > 1:492432/1‑143 (MQ=255) tGCCGGATAGGGCGCTCGCGCCACATCCGGCATACACAGATCGCATCAGGTTATGTTTCGTTTGTCATCAACCTCAAATCACCAAAATAATTTATTTTATTTCAATAAGATAGATAAGTAGAACTGAGAACGCCATGCTCTCg > 2:294333/1‑143 (MQ=255) ccGGATGTGGCGCGAGCGCCCTATCCGGCATACACAGATCGCATCAGGTTATGTTTCGTTTGTCATCAACCTCAAATCACCAAAATAATTTATTTTATTTCCATAAGATAGATAAGTAGAACTGCGAACGCCCTGCTCTCGtc > 2:315843/1‑142 (MQ=255) ccGGATGTGGCGCGAGCGCCCTATCCGGCATACACAGATCGCATCAGGTTATGTTTCGTTTGTCATCAACCTCAAATCACCAAAATAATTTATTTTATTTCAATAAGATAGATAAGTAGAACTGAGAAAGCCATGCTCTCGtt > 2:189122/1‑143 (MQ=255) cGGATGTGGCGCGAGCGCCCTATCCGGCATACACAGATCGCATCAGGTTATGTTTCGTTTGTCATCAACCTCAAATCACCAAAATAATTTATTTTATTTCAATAAGATAGATAAGTAGAACTGAGAAAGCCATGCTCTCGttt > 2:296332/1‑143 (MQ=255) aTGTGGCGCGAGCGCCCTATCCGGCATACACAGATCGCATCAGGTTATGTTTCGTTTGTCATCAACCTCAAATCACCAAAATAATTTATTTTATTTCAATAAGATAGATAAGTAGAACTGAGAAAGCCATGCTCTCGTTTCCt < 1:84634/143‑1 (MQ=255) gtgGCGCGAGCGCCCTATCCGGCATACACAGATCGCATCAGGTTATGTTTCGTTTGTCATCAACCTCAAATCACCAAAATAATTTATTTTATTTCAATAAGATAGATAAGTAGAACTGAGAAAGCCATGCTCTCGTTTCCTaa < 2:375791/143‑1 (MQ=255) | TGGCTGTACTGCTTTTAGCCTTACGGTATAAGAAGTTGCGTTCCCTTTTTAACCGGTTACGCAAGCAGGTTCACAACAAATTATCTGTATAACTTCATATTGCCGACAAAGTACGCGTTATCCCTGCCGGATGTGGCGCGAGCGCCCTATCCGGCATACACAGATCGCATCAGGTTATGTTTCGTTTGTCATCAACCTCAAATCACCAAAATAATTTATTTTATTTCAATAAGATAGATAAGTAGAACTGAGAAAGCCATGCTCTCGTTTCCTAA > NC_000913/2041027‑2041301 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |