Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_000913 | 2,028,158 | G→C | 16.2% | intergenic (‑141/+30) | yedQ ← / ← yodC | putative membrane‑anchored diguanylate cyclase/uncharacterized protein |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 2,028,158 | 0 | G | C | 16.2% | 201.9 / 32.6 | 74 | intergenic (‑141/+30) | yedQ/yodC | putative membrane‑anchored diguanylate cyclase/uncharacterized protein |
Reads supporting (aligned to +/- strand): ref base G (25/37); new base C (3/9); total (28/46) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 5.17e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.41e-01 |
CCTGGCGATAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCTTG > NC_000913/2028059‑2028257 | ccTGGCGATAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGAtg > 1:1999792/1‑100 (MQ=255) ccTGGCGATAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGAtg > 1:8045/1‑100 (MQ=255) ccTGGCGATAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATCATAATACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGAtg < 1:401655/100‑1 (MQ=255) gCGATAGCGTACCGCAAAACAGCGATAATTTTATAAGGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGAt < 2:1210626/100‑1 (MQ=255) cGATAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATc > 2:57624/1‑100 (MQ=255) cGATAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATc > 1:827032/1‑100 (MQ=255) cGATAGCGTACCGCAAAACAGCGATAATTTTATAATGGCGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATc > 2:1482500/1‑100 (MQ=255) gATAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATcc < 2:231042/100‑1 (MQ=255) tggcatACCGCAAAACGGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCgg < 2:1977422/95‑1 (MQ=255) aGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGc < 1:2017947/100‑1 (MQ=255) gCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGcc < 1:943250/100‑1 (MQ=255) cGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACATCAGGCCGGATCACATCCGGCCt < 1:2507517/100‑1 (MQ=38) cGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATCACATCCGGCCt < 1:649895/100‑1 (MQ=255) cGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATCACATCCGGCCt < 1:1308699/100‑1 (MQ=255) tACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGa < 2:838847/100‑1 (MQ=255) tACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATCACATCCGGCCTgg < 1:1473175/100‑2 (MQ=255) tACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATCACATCCGGCCTgg < 2:136372/100‑2 (MQ=255) tACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATCACATCCGGCCTgg < 1:2103449/100‑2 (MQ=255) tACCGCAAAACAGCGATAATTGTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGa < 2:919857/100‑1 (MQ=255) aCCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGAt < 2:1212842/100‑1 (MQ=255) aCCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGAt < 1:2162698/100‑1 (MQ=255) ccGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATg > 2:1251475/1‑100 (MQ=255) ccGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATg < 1:2112595/100‑1 (MQ=255) ccGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATg < 2:851382/100‑1 (MQ=255) ccGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATCACATCCGGCCTGGTg < 1:1468016/100‑1 (MQ=38) ccGCAAAACAGCGATAATGTTATAATGGTGGCACGATAAGACGCGACTATCATCCATTGCAGGCATTAGCACACCAGGCCGGATCACATCCGGCCTgg < 2:861522/98‑2 (MQ=255) cGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGt < 1:801867/100‑1 (MQ=255) gCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATCACATCCGGCCTGgtgtg < 2:635616/100‑2 (MQ=38) cAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACGAGGGCGGATGTGATCCGGCCTGATGTTa > 1:2597619/1‑100 (MQ=255) aaaaCAACGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTaa < 1:788839/100‑1 (MQ=255) aaaCAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAAt < 1:427031/100‑1 (MQ=255) aaaCAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAAt < 1:427006/100‑1 (MQ=255) aaCAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATg < 2:609471/100‑1 (MQ=255) aCAGCGATAATTTTATAATGGTGGCATGATAACACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAAt > 2:2119954/1‑98 (MQ=255) aGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGTACACCAGGGCGGATGTGTTCCGGCCTGATGTTAATGAt < 2:1000477/99‑1 (MQ=255) ttttATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTca > 2:1047799/1‑100 (MQ=255) tttATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTcac < 1:2527211/100‑1 (MQ=255) tttATAATGGTGGCATGATAAGACGCGACTATCATACATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTcac < 2:158826/100‑1 (MQ=255) ttATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTcaca < 1:783755/100‑1 (MQ=255) ttATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTcaca > 2:1678494/1‑100 (MQ=255) ataATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACAct > 2:2303315/1‑100 (MQ=255) taATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACActt < 2:29017/100‑1 (MQ=255) aaggtggCATGACAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACCcttc < 2:2339982/97‑1 (MQ=255) ggtggCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACttcttc < 1:948504/100‑1 (MQ=255) gtggCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTcc < 1:513427/100‑1 (MQ=255) ggCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGc < 1:1045035/100‑1 (MQ=255) tcATTATATGATGCGACTATCATTCATTGCAGGCATTAGCACATCAGGCCGTATGTGATCCTTCCTGATCTTAATGTTACCTTTTCCCACTTCTTCCGc < 2:367772/98‑1 (MQ=255) cATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAg < 2:459317/100‑1 (MQ=255) gATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAAc < 2:1546983/100‑1 (MQ=255) tAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGa > 2:1090334/1‑100 (MQ=255) aGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGAct > 2:1557527/1‑100 (MQ=255) aCGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGActcc < 1:1150972/100‑1 (MQ=255) cgcgACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGActccc > 2:1129960/1‑100 (MQ=255) cgcgACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGActcc < 2:1121516/99‑1 (MQ=255) ttcattGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGaa < 1:1423426/100‑1 (MQ=255) tcattGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAAc < 1:1176937/100‑1 (MQ=255) gCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCt > 1:665385/1‑100 (MQ=255) gCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTCCCGCAGAACGACTCCCCTCCCCCGGAACAAGCt > 2:298275/1‑100 (MQ=255) cAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTc > 2:47325/1‑100 (MQ=255) aGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCg > 1:1851740/1‑100 (MQ=255) aGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCg < 2:1056806/100‑1 (MQ=255) ggCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCgg > 1:852457/1‑100 (MQ=255) cATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGtt > 1:1606158/1‑100 (MQ=255) cATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGtt < 2:999630/100‑1 (MQ=255) cattaacatcAGGCCGGATCACATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAGCGACTCCCCTCCCCCGGAACAAGCTCGGttt > 2:1828752/10‑99 (MQ=25) attaacatcAGGCCGGATCACATCCGGCCTGGTGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCa > 2:372688/9‑100 (MQ=25) aGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCAt < 1:563293/100‑1 (MQ=255) aGCACACCACGCCGGATGTGATCCGCCCGTATGTTAATGATCCCCTTTCACACTCCTTCCGCAGACCGGCTCCCCTCCCCCGGAACAAGCTCGGTTTCa < 2:1715415/99‑1 (MQ=255) cacaCCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGa > 1:1453765/1‑100 (MQ=255) cAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGc > 1:2225066/1‑100 (MQ=255) aGGCCGGATCACATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCt > 2:2143815/1‑100 (MQ=255) ggCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCtt > 1:229249/1‑100 (MQ=255) cGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTccc > 1:572372/1‑100 (MQ=255) cGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTccc > 2:2202056/1‑100 (MQ=255) ggATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCg > 1:660386/1‑100 (MQ=255) gATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCTCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGc < 1:1839806/100‑1 (MQ=255) gATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGc > 2:1591785/1‑100 (MQ=255) tgtgATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCtt > 2:17996/1‑100 (MQ=255) gtgATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCTTg < 2:1550780/100‑1 (MQ=255) | CCTGGCGATAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCTTG > NC_000913/2028059‑2028257 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |