Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_000913 | 2,028,159 | T→A | 16.9% | intergenic (‑142/+29) | 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,159 | 0 | T | A | 16.9% | 194.4 / 29.5 | 77 | intergenic (‑142/+29) | yedQ/yodC | putative membrane‑anchored diguanylate cyclase/uncharacterized protein |
Reads supporting (aligned to +/- strand): ref base T (27/37); new base A (4/9); total (31/46) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 5.44e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.52e-01 |
CGATAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCTTGA > NC_000913/2028064‑2028258 | cGATAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATc > 1:2710130/1‑100 (MQ=255) gATAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATcc < 2:1871456/100‑1 (MQ=255) tAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCgg < 1:2288658/100‑1 (MQ=255) gCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTA‑‑ACATCAGGCCGGATCACATCCGGCCt < 1:1314623/99‑1 (MQ=37) cGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCt < 2:2092825/100‑1 (MQ=255) cGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCt < 2:1465445/100‑1 (MQ=255) cGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATCACATCCGGCCt < 1:2705276/100‑1 (MQ=255) gTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTg < 2:27236/100‑1 (MQ=255) gTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTg < 1:1606336/100‑1 (MQ=255) gTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATCACATCCGGccn < 1:2581947/99‑2 (MQ=255) gTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATCACATCCGGCCTg < 2:2610131/100‑1 (MQ=255) tACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACATCAGGCCGGATCACATCCGGCCTgg < 1:2674046/100‑2 (MQ=38) tACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGa < 2:2179293/100‑1 (MQ=255) tACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATCACATCCGGCCTgg < 1:2662910/100‑2 (MQ=255) tACCGCAAAACAGCGATAATGTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATCACATCCGGCCTgg < 1:1311595/100‑2 (MQ=255) aCCGCTAACCCGCGATAACCTTATGATGGCATCACAATAAGACGCGACTATAACTCATTGCCCGCATTAGCACACCAGGCCCGATCACATCCGGCCtag < 1:1738703/99‑3 (MQ=255) ccGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATg < 1:2262302/100‑1 (MQ=255) ccGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATg > 2:1208200/1‑100 (MQ=255) gCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGtt > 2:2852745/1‑100 (MQ=255) aaaaCAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTaa < 2:566588/100‑1 (MQ=255) aaaCAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAAt < 2:1581921/100‑1 (MQ=255) aaCAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATg > 2:8027/1‑100 (MQ=255) aaCAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATg < 2:2199636/100‑1 (MQ=255) cAGCGATAATTTTATAATGGTGGCATGATGAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATCACATCCGGCCTGgtgtgctaatg < 2:740457/100‑8 (MQ=25) aGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACAACAGGCCGGATGTGATCCGGCCTGATGTTAATGATa < 1:2701554/100‑1 (MQ=255) gCGAAAATTTTATAATGGTAGCATGAGAAGACGCAGCTATCGTTCATTGGAGGCAGTAGCAGACCAGGCGGGGTGTGATCCGGCCTGATGTTAATGATa > 1:2431096/1‑99 (MQ=255) aTTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTc > 1:87738/1‑100 (MQ=255) tttATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTcac < 2:147059/100‑1 (MQ=255) ataATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACAct > 1:2048960/1‑100 (MQ=255) taATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACActt > 1:508757/1‑100 (MQ=255) aTGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCCGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACActtct < 1:1266032/100‑1 (MQ=255) ggtggCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACttcttc < 2:524279/100‑1 (MQ=255) ggtggCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACttcttc > 2:2440396/1‑100 (MQ=255) gtggCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTcc > 1:1536432/1‑100 (MQ=255) gtggCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTcc < 2:613592/100‑1 (MQ=255) tggCTTTATAAGTCGCGACTTTCATTCATTCCAGACTTTATCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATCATACCCTTTCACACttcttc < 1:474728/98‑1 (MQ=255) ggCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGc < 1:2009363/100‑1 (MQ=255) ggCAGGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGc < 1:2544962/100‑1 (MQ=255) aTAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACg < 1:1926278/100‑1 (MQ=255) aTAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACg < 2:155846/100‑1 (MQ=255) aGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGAct < 2:1317523/100‑1 (MQ=255) gACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGActc < 1:677062/100‑1 (MQ=255) aCGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGActcc < 2:1438251/100‑1 (MQ=255) tcattcatTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCg < 2:852702/100‑1 (MQ=255) tcattcatTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCACCCCCg < 1:286870/100‑1 (MQ=255) cattcattGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCgg > 2:1603082/1‑100 (MQ=255) attcattGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGa > 1:2485233/1‑100 (MQ=255) ttcattGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGaa < 1:448949/100‑1 (MQ=255) ttcattGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGaa < 1:1706686/100‑1 (MQ=255) tcattGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAAc < 2:1997112/100‑1 (MQ=255) cattGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACa > 1:2844083/1‑100 (MQ=255) nttGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCACAGAACGACTCCCCTCCCCCGGAACaa > 2:1278265/2‑100 (MQ=255) ttGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACACTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAg > 2:971934/1‑100 (MQ=255) gCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCt > 2:2318489/1‑100 (MQ=255) tggCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCg < 1:200200/99‑1 (MQ=255) ggCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCgg > 1:153037/1‑100 (MQ=255) ggCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCgg < 1:1603076/100‑1 (MQ=255) ggCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCga > 2:200201/1‑99 (MQ=255) cATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGtt > 1:2408077/1‑100 (MQ=255) aTTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGttt < 1:916119/100‑1 (MQ=255) ttAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCAGGAACAAGCTCGGTTTc > 1:2173563/1‑100 (MQ=255) tAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCGCGGACCAAGCTCGGTTTCa > 1:313805/1‑100 (MQ=255) attaacatcAGGCCGGATCACATCCGGCCTGGTGTTAATGATACCCTTTCACCCTTCTTCCGCAGAACGACTCCAATCCCCCGGGACAAGCTCGGTTTCa > 1:1395607/9‑100 (MQ=21) attaacatcAGGCCGGATCACATCCGGCCTGGTGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCa > 2:2236457/9‑100 (MQ=25) ttaacatcAGGCCGGATCACATCCGGCCTGGTGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCTCCTCCCCCGGAACAAGCTCGGTTTCAt > 1:2701088/8‑100 (MQ=25) gCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATg < 2:2277311/100‑1 (MQ=255) gCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATg > 1:2279952/1‑100 (MQ=255) gCACACCAGGCCGGATGTGATCCGGCCTGATGTTAAGGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATg < 2:2048971/100‑1 (MQ=255) acacCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATg > 1:2124009/1‑98 (MQ=255) acacCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGaa > 2:2475565/1‑100 (MQ=255) catcAGGCCGGATCACATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGaaa > 2:23776/4‑100 (MQ=38) ccAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAg < 2:1215490/100‑1 (MQ=255) aGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCt > 2:1023691/1‑100 (MQ=255) aGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCt > 2:1981149/1‑100 (MQ=255) ggCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCtt > 2:219023/1‑100 (MQ=255) gCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTc < 2:2385824/100‑1 (MQ=255) aTGTGATCCGGCCTGATGTTAACGATACCCCTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCTGAACAAGCTCGGTTTCATGAAAAGCTTc < 1:2719034/95‑1 (MQ=255) tgtgATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGCACGACTCCCCTCCCCCGGAACATGCTCGGTTTCATGACCAGCTTCCCGCtt > 1:2098408/1‑100 (MQ=255) tgtgATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCtt > 2:65498/1‑100 (MQ=255) gtgATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCTTg > 2:2437840/1‑100 (MQ=255) gtgATCCGGCCTGATATGAATGATACCCTGTCACACTTCGTCCGCAGAACGGCTCCCCTCCCCCGGAACAGGCTCTGTTGCGTGAAAAGCTTCCCGCt < 1:1304790/98‑1 (MQ=255) tgATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCTTGa > 2:2231458/1‑100 (MQ=255) | CGATAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCTTGA > NC_000913/2028064‑2028258 |
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: ‑ |