Predicted mutation | ||||||
---|---|---|---|---|---|---|
evidence | seq id | position | mutation | annotation | gene | description |
RA | NC_000913 | 3,556,033 | A→G | Y273H (TAC→CAC) | rtcA ← | RNA 3'‑terminal phosphate cyclase |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 3,556,033 | 0 | A | G | 98.4% | 216.5 / ‑5.5 | 63 | Y273H (TAC→CAC) | rtcA | RNA 3'‑terminal phosphate cyclase |
Reads supporting (aligned to +/- strand): ref base A (1/0); new base G (29/33); total (30/33) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 4.76e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.92e-01 |
CTCCACCACCGCGATATTGGTCAGCAGATGGCATGAGGGATGGGCGACCGTAAATTCCCCCGCGCCCGCCAGCGCCATCGGTAGCACCAGTTGGTCAGCAAGATATT‑CCCCCACCGCCGCCGTGCTTGCCAGGTAGCGTTTCACCTCTTTCACCAACTGTGCCGCAACCACCTCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACC > NC_000913/3555899‑3556153 | cTCCACCACCGCGATATTGGTCAGCAGATGGCATGAGGGATGGGCGACCGTAAATTCCCCCGCGCCCGCCAGCGCCATCGGTAGCACCAGTTGGTCAGCAAGATATT‑CCCCCACCGCCGCCGTGCTTGCCAGGTGGCGt > 2:557056/1‑139 (MQ=255) cTCCACCACCGCGATATTGGTCAGCAGATGGCATGAGGGATGGGCGACCGTAAATTCCCCCGCGCCCGCCAGCGCCATCGGTAGCACCAGTTGGTCAGCAAGATATT‑CCCCCACCGCCGCCGTGCTTGCCAGGTGGCGt > 2:716378/1‑139 (MQ=255) cTCCACCACCGCGATATTGGTCAGCAGATGGCATGAGGGATGGGCGACCGTAAATTCCCCCGCGCCCGCCAGCGCCATCGGTAGCACCAGTTGGTCAGCAAGATATT‑CCCCCACCGCCACCGTGCTTGCCAGGTGGCGt > 2:314778/1‑139 (MQ=255) ccaccGCGATATTGGTCAGCAGATGGCATGAGGGATGGGCGACCGTAAATTCCCCCGCGCCCGCCAGCGCCATCGGTAGCACCAGTTGGTCAGCAAGATATT‑CCCCCACCGCCGCCGTGCTTGCCAGGTGGCGTTTCAc < 1:557056/139‑1 (MQ=255) ccaccGCGATATTGGTCAGCAGATGGCATGAGGGATGGGCGACCGTAAATTCCCCCGCGCCCGCCAGCGCCATCGGTAGCACCAGTTGGTCAGCAAGATATT‑CCCCCACCGCCGCCGTGCTTGCCAGGTGGCGTTTCAc < 2:1613197/139‑1 (MQ=255) cgATATTGGTCAGCAGATGGCATGAGGGATGGGCGACCGTAAATTCCCCCGCGCCCGCCAGCGCCATCGGTAGCACCAGTTGGTCAGCAAGATATT‑CCCCCACCGCCGCCGTGCTTGCCAGGTGGCGTTTCACCTCttt > 1:509439/1‑139 (MQ=255) cgATATTGGTCAGCAGATGGCATGAGGGATGGGCGACCGTAAATTCCCCCGCGCCCGCCAGCGCCATCGGTAGCACCAGTTGGTCAGCAAGATATT‑CCCCCACCGCCGCCGTGCTTGCCAGGTGGCGTTTCACCTCttt > 2:129779/1‑139 (MQ=255) cgATATTGGTCAGCAGATGGCATGAGGGATGGGCGACCGTAAATTCCCCCGCGCCCGCCAGCGCCATCGGTAGCACCAGTTGGTCAGCAAGATATT‑CCCCCACCGCCGCCGTGCTTGCCAGGTGGCGTTTCACCTCttt > 2:551895/1‑139 (MQ=255) agcagATGGCATGAGGGATGGGCGACCGTAAATTCCCCCGCGCCCGCCAGCGCCATCGGTAGCACCAGTTGGTCAGCAAGATATT‑CCCCCACCGCCGCCGTGCTTGCCAGGTGGCGTTTCACCTCTTTCACCAACtgtg < 2:509439/139‑1 (MQ=255) ggCATGAGTGATGGGCGACCGTATATTCCCCCGCGCCCGCCAGCGCCAGCGGTAGCACCAGTTGGTCAGCAAGATATT‑CCCCCACCGCCGCCGTGCTTGCCAGGTGGCGTTTCACCTCTTTCACCAACTGTGCCGCAac < 1:1014437/139‑1 (MQ=255) ggCATGAGGGATGGGCGACCGTAAATTCCCCCGCGCCCGCCAGCGCCATCGGTAGCACCAGTTGGTCAGCAAGATATT‑CCCCCACCGCCGCCGTGCTTGCCAGGTGGCGTTTCACCTCTTTCACCAACTGTGCCGCAac < 2:611650/139‑1 (MQ=255) ggCATGAGGGATGGGCGACCGTAAATTCCCCCGCGCCCGCCAGCGCCATCGGTAGCACCAGTTGGTCAGCAAGATATT‑CCCCCACCGCCGCCGTGCTTGCCAGGTGGCGTTTCACCTCTTTCACCAACTGTGCCGCAac < 1:551895/139‑1 (MQ=255) gggCGACCGTAAATTCCCCCGCGCCCGCCAGCGCCATCGGTAGCACCAGTTGGTCAGCAAGATATT‑CCCCCACCGCCGCCGTGCTTGCCAGGTGGCGTTTCACCTCTTTCACCAACTGTGCCGCAac < 1:1328888/127‑1 (MQ=255) gggCGACCGTAAATTCCCCCGCGCCCGCCAGCGCCATCGGTAGCACCAGTTGGTCAGCAAGATATT‑CCCCCACCGCCGCCGTGCTTGCCAGGTGGCGTTTCACCTCTTTCACCAACTGTGCCGCAac > 2:1328888/1‑127 (MQ=255) ggCGACCGTAAATTCCCCCGCGCCCGCCAGCGCCATCGGTAGCACCAGTTGGTCAGCAAGATATT‑CCCCCACCGCCGCCGTGCTTGCCAGGTGGCGTTTCACCTCTTTCACCAACTGTGCCGCAACCACCTCGGCACTg < 1:1618883/139‑1 (MQ=255) cGTAAATTCCCCCGCGCCCGCCAGCGCCATCGGTAGCACCAGTTGGTCAGCAAGATATT‑‑CCCCACCGCCGCCGTGCTTGCCAGGTGGCGTTTCACCTCTTTCACCAACTGTGCCGCAACCACCTCGGCACTGACGCGCt > 1:579847/1‑139 (MQ=255) tCCCCCGCGCCCGCCAGCGCCATCGGTAGCACCAGTTGGTCAGCAAGATATT‑‑CCCCACCGCCGCCGTGCTTGCCAGGTGGCGTTTCACCTCTTTCACCAACTGTGCCGCAACCACCTCGGCACTGACGCGCTTTTCAcc < 2:579847/139‑1 (MQ=255) ccccGCGCCCGCCAGCGCCATCGGTAGCACCAGTTGGTCAGCAAGATATT‑CCCCCACCGCCGCCGTGCTTGCCAGGTGGCGTTTCACCTCTTTCACCAAc > 1:1294818/1‑100 (MQ=255) ccccGCGCCCGCCAGCGCCATCGGTAGCACCAGTTGGTCAGCAAGATATT‑CCCCCACCGCCGCCGTGCTTGCCAGGTGGCGTTTCACCTCTTTCACCAAc < 2:1294818/100‑1 (MQ=255) cccGCGCCCGCCAGCGCCATCGGTAGCACCAGTTGGTCAGCAAGATATT‑CCCCCACCGCCGCCGTGCTTGCCAGGTGGCGTTTCACct > 2:1620665/1‑88 (MQ=255) cccGCGCCCGCCAGCGCCATCGGTAGCACCAGTTGGTCAGCAAGATATT‑CCCCCACCGCCGCCGTGCTTGCCAGGTGGCGTTTCACct < 1:1620665/88‑1 (MQ=255) cccGCGCCCGCCAGCGCCATCGGTAGCACCAGTTGGTCAGCAAGATATT‑CCCCCACCGCCGCCGTGCTTGCCAGGTGGCGTTTCACCTCTTTCACCAACTGTGCCGCAACCACCTCGGCACTGACGCGCTTTTCACcga < 1:716378/139‑1 (MQ=255) cgccAGCGCCATCGGTAGCACCAGTTGGTCAGCAAGATATT‑CCCCCACCGCCGCCGTGCTTGCCAGGTGGCGTTTCACCTCTTTCACCAACTGTGcc > 1:1176528/1‑97 (MQ=255) cgccAGCGCCATCGGTAGCACCAGTTGGTCAGCAAGATATT‑CCCCCACCGCCGCCGTGCTTGCCAGGTGGCGTTTCACCTCTTTCACCAACTGTGcc < 2:1176528/97‑1 (MQ=255) aGCGCCATCGGTAGCACCAGTTGGTCAGCAAGATATTCCCCCCACCGCCGCCGTGCTTGCCAGGTGGCGTTTCACCTCTTTCACCAACTGTGCCGCAACCACCTCGGCACTGACGCGCTTTTCACCGACGACAAAAAAg < 2:165204/139‑1 (MQ=255) gcCATCGGTAGCACCAGTTGGTCAGCAAGATATT‑CCCCCACCGCCGCCGTGCTTGCCAGGTGGCGTTTCACCTCTTTCACCAACTGTGCCGCAACCACCTCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGCGtt < 2:695570/139‑1 (MQ=255) ccATCGGTAGCACCAGTTGGTCAGCAAGATATT‑CCCCCACCGCCGCCGTGCTTGCCAGGTGGCGTTTCACCTCTTTCACCAACTGTGCCGCAACCACCTCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGCGTTc < 1:830199/139‑1 (MQ=255) gTTGGTCAGCAAGATATT‑CCCCCACCGCCGCCGTGCTTGCCAGGTGGCGTTTCACCTCTTTCACCAACTGTGCCGCAac > 2:1329320/1‑79 (MQ=255) gTTGGTCAGCAAGATATT‑CCCCCACCGCCGCCGTGCTTGCCAGGTGGCGTTTCACCTCTTTCACCAACTGTGCCGCAac < 1:1329320/79‑1 (MQ=255) ttGGTCAGCAAGATATT‑CCCCCACCGCCGCCGTGCTTGCCAGGTGGCGTTTCACct < 1:1226812/56‑1 (MQ=255) ttGGTCAGCAAGATATT‑CCCCCACCGCCGCCGTGCTTGCCAGGTGGCGTTTCACct > 2:1226812/1‑56 (MQ=255) ttGGTCAGCAAGATATT‑CCCCCACCGCCGCCGTGCTTGCCAGGTGGCGTTTCACCTCTTTCACCAACTGTGCCGCAaccac > 1:1267862/1‑81 (MQ=255) ttGGTCAGCAAGATATT‑CCCCCACCGCCGCCGTGCTTGCCAGGTGGCGTTTCACCTCTTTCACCAACTGTGCCGCAaccac < 2:1267862/81‑1 (MQ=255) aGCAAGATATT‑CCCCCACCGCCGCCGTGCTTGCCAGGTGGCGTTTCACCTCTTTCACCAACtgtg > 1:1302670/1‑65 (MQ=255) aGCAAGATATT‑CCCCCACCGCCGCCGTGCTTGCCAGGTGGCGTTTCACCTCTTTCACCAACtgtg < 2:1302670/65‑1 (MQ=255) cAAGATATT‑CCCCCACCGCCGCCGTGCTTGCCAGGTGGCGTTTCACCTCTTTCACCAACTGTGCCGCAaccac < 2:476947/73‑1 (MQ=255) cAAGATATT‑CCCCCACCGCCGCCGTGCTTGCCAGGTGGCGTTTCACCTCTTTCACCAACTGTGCCGCAaccac > 1:476947/1‑73 (MQ=255) atT‑CCCCCACCGCCGCCGTGCTTGCCAGGTGGCGTTTCAc > 2:844967/1‑40 (MQ=255) atT‑CCCCCACCGCCGCCGTGCTTGCCAGGTGGCGTTTCAc < 1:844967/40‑1 (MQ=255) t‑CCCCCACCGCCGCCGTGCTTGCCAGGTGGCGTTTCACCTCTTTCACCAACTGTGCCGCAACCACCTCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCgaga > 2:189427/1‑139 (MQ=255) ccccACCGCCGCCGTGCTTGCCAGGTGGCGTTTCACCTCTTTCACCAACTGTGCCGCAACCACCTCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACg > 1:501560/1‑139 (MQ=255) ccccACCGCCGCCGTGCTTGCCAGGTGGCGTTTCACCTCTTTCACCAACTGTGCCGCAACCACCTCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACg > 1:1057249/1‑139 (MQ=255) aCCGCCGCCGTGCTTGCCAGGTGGCGTTTCACCTCTTTCACCAACTGTGCCGCAACCACCTCgg < 1:1382290/64‑1 (MQ=255) aCCGCCGCCGTGCTTGCCAGGTGGCGTTTCACCTCTTTCACCAACTGTGCCGCAACCACCTCgg > 2:1382290/1‑64 (MQ=255) ccgccgccgTGCTTGCCAGGTGGCGTTTCACCTCTTTCACCAACTGTGCCGCAaccac < 2:1004565/58‑1 (MQ=255) ccgccgccgTGCTTGCCAGGTGGCGTTTCACCTCTTTCACCAACTGTGCCGCAACCACCt < 2:795819/60‑1 (MQ=255) ccgccgccgTGCTTGCCAGGTGGCGTTTCACCTCTTTCACCAACTGTGCCGCAACCACCt > 1:795819/1‑60 (MQ=255) ccgccgccgTGCTTGCCAGGTGGCGTTTCACCTCTTTCACCAACTGTGCCGCAACCACCTCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTAtt > 1:1536091/1‑139 (MQ=255) ccgccgccgTGCTTGCCAGGTGGCGTTTCACCTCTTTCACCAACTGTGCCGCAACCACCTCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTAtt > 2:1549866/1‑139 (MQ=255) cagccgccgTGCTTGCCAGGTAGCGTTTCACCTCTTTCACCAACTGTGCCGCAaccac > 1:1004565/3‑58 (MQ=255) ccgccgTGCTTGCCAGGTGGCGTTTCTCCTCTTTCACCAACTGTGCCGCAACCACCTCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTAcc < 2:442260/139‑1 (MQ=255) ccgccgTGCTTGCCAGGTGGCGTTTCACCTCTTTCACCAACTGTGCCGCAACCACCTCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGCGtt > 2:681786/1‑96 (MQ=255) ccgccgTGCTTGCCAGGTGGCGTTTCACCTCTTTCACCAACTGTGCCGCAACCACCTCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGCGtt < 1:681786/96‑1 (MQ=255) ccgccgTGCTTGCCAGGTGGCGTTTCACCTCTTTCACCAACTGTGCCGCAACCACCTCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTAcc < 1:82732/139‑1 (MQ=255) ccgccgTGCTTGCCAGGTGGCGTTTCACCTCTTTCACCAACTGTGCCGCAACCACCTCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTAcc < 1:1549866/139‑1 (MQ=255) ccgccgTGCTTGCCAGGTGGCGTTTCACCTCTTTCACCAACTGTGCCGCAACCACCTCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTAcc < 2:1536091/139‑1 (MQ=255) ccgccgTGCTTGCCAGGTGGCGTTTCACCTCTTTCACCAACTGTGCCGCAACCACCTCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTAcc > 2:82732/1‑139 (MQ=255) ccgccgTGCTTGCCAGGTGGCGTTACACCTCTTTCACCAACTGTGCCGCAACCACCTCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTAcc > 2:898400/1‑139 (MQ=255) ccaccgTGCTTGCCAGGTGGCGTTTCACCTCTTTCACCAACTGTGCCGCAACCACCTCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTAcc < 1:314778/136‑1 (MQ=255) cgccgTGCTTGCCAGGTGGCGTTTCACCTCTTTCACCAACTGTGCCGCAACCACCTCGGCACTGACGCGCTTTTCACcga > 1:1112747/1‑80 (MQ=255) cgccgTGCTTGCCAGGTGGCGTTTCACCTCTTTCACCAACTGTGCCGAAACCACCTCGGCACTGACGCGCTTTTCACcga < 2:1112747/80‑1 (MQ=255) gccgTGCTTGCCAGGTGGCGTTTCACCTCTTTCAcc < 2:1619470/36‑1 (MQ=255) gccgTGCTTGCCAGGTGGCGTTTCACCTCTTTCAcc > 1:1619470/1‑36 (MQ=255) | CTCCACCACCGCGATATTGGTCAGCAGATGGCATGAGGGATGGGCGACCGTAAATTCCCCCGCGCCCGCCAGCGCCATCGGTAGCACCAGTTGGTCAGCAAGATATT‑CCCCCACCGCCGCCGTGCTTGCCAGGTAGCGTTTCACCTCTTTCACCAACTGTGCCGCAACCACCTCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACC > NC_000913/3555899‑3556153 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |