Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 2,546,869 | G→A | 10.9% | intergenic (+25/+15) | PP_2238 → / ← PP_5504 | metallopeptidase/hypothetical protein |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 2,546,869 | 0 | G | A | 10.9% | 137.8 / 7.8 | 55 | intergenic (+25/+15) | PP_2238/PP_5504 | metallopeptidase/hypothetical protein |
Reads supporting (aligned to +/- strand): ref base G (26/23); new base A (2/4); total (28/27) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 4.22e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.41e-01 |
GACCAAGGAGGAGGTGGAGTGGCTGAACGGCTACCACGCCTGCGTACATGAGCGCCTGGCGCCCTTGCTCCAAGGGGATGCGCTAGCCTGGCTGGAGATGCGTACTGCACCGCTGTAACGCAGCGGAAAAAGGGCAGTTTGCGAACTGCCCTTTTCTCTATTTACAGATGATGATCATGCTCCGGCTGGTGTAGCCCGCCGGGTTGATGCCGAACAAATAGTCACCGGCCTCTTCATCGCTATCGCCCGAGCGCGCAATGACCCGGTAGCCGCGCGTGCTGCAG > NC_002947/2546727‑2547010 | gACCAAGGAGGAGGTGGAGTGGCTGAACGGCTACCACGCCTGCGTACATGAGCGCCTGGCGCCCTTGCTCCAAGGGGATGCGCTAGCCTGGCTGGAGATGCGTACTGCACCGCTGTAACGCAGCGGAAAAAGGGCAGTTTGCGAACTGcc < 1:645122/150‑1 (MQ=255) cAAGGAGGAGGTGGAGTGGCTGAACGGCTACCACGCCTGCGTACATGAGCGCCTGGCGCCCTTGCTCCAAGGGGATGCGCTAGCCTGGCTGGAGATGCGTACTGCACCGCTGTAACGCAGCGGAAAAAGGGCAGTTCGCAAACTGCCCtt > 1:1232463/1‑150 (MQ=255) aaGGAGGAGGTGGAGTGGCTGAACGGCTACCACGCCTGCGTACATGAGCGCCTGGCGCCCTTGCTCCAAGGGGATGCGCTAGCCTGGCTGGAGATGCGTACTGCACCGCTGTAACGCAGCGGAAAAAGGGCAGTTCGCAAACTGCCCttt < 1:1921864/150‑1 (MQ=255) ggaggaggTGGAGTGGCTGAACGGCTACCACGCCTGCGTACATGAGCGCCTGGCGCCCTTGCTCCAAGGGGGTGCGCTAGCCTGGCTGGAGATGCGTACTGCACCGCTGTAACGCAGCGGAAAAAGGGCAGTTTGCGAACTGCCCTTTtc > 1:1648600/1‑150 (MQ=255) ggaggaggTGGAGTGGCTGAACGGCTACCACGCCTGCGTACATGAGCGCCTGGCGCCCTTGCTCCAAGGGGATGCGCTAGCCTGGCTGGAGATGCGTACTGCACCGCTGTAACGCAGCGGAAAAAGGGCAGTTTGCGAACTGCCCTTTtc > 1:454487/1‑150 (MQ=255) gaggaggTGGAGTGGCTGAACGGCTACCACGCCTGCGTACATGAGCGCCTGGCGACCTGGCTCCAAGGGGATGCGCTAGCCTGGCTGGAGATGCGTACTGCACCGCTGTAACGCAGCGGAAAAAGGGCAGTTCGCAAACTGCCCtttttc < 2:1232463/150‑3 (MQ=255) ggaggTGGAGTGGCTGAACGGCTACCACGCCTGCGTACATGAGCGCCTGGCGCCCTTGCTCCAAGGGGATGCGCTAGCCTGGCTGGAGATGCGTACTGCACCGCTGTAACGCAGCGGAAAAAGGGCAGTTCGCAAACTGCCCtttttccg < 2:540162/150‑5 (MQ=255) gaggTGGAGTGGCTGAACGGCTACCACGCCTGCGTACATGAGCGCCTGGCGCCCTTGCTCCAAGGGGATGCGCTAGCCTGGCTGGAGATGCGTACTGCACCGCTGTAACGCAGCGGAAAAAGGGCAGTTTGCGAACTGCCCTTTTCTCTa > 2:56529/1‑150 (MQ=255) tGGCTGAACGGCTACCACGCCTGCGTACATGAGCGCCTGGCGCCCTTGCTCCAAGGGGATGCGCTAGCCTGGCTGGAGATGCGTACTGCACCGCTGTAACGCAGCGGAAAAAGGGCAGCTTGCGAACTGCCCTTTTCTCTATTTACAgat < 1:1601882/150‑1 (MQ=255) aaCGGCTACCACGCCTGCGTACATGAGCGCCTGGCGCCCTTGCTCCAAGGGGATGCGCTAGCCCGGCTGGAGATGCGTACTGCACCGCTGTAACGCAGCGGAAAAAGGGCAGTTTGCGAACTGCCCTTTTCTCTATTTACAgatgatgat < 2:468201/150‑1 (MQ=255) aCGGCTACCACGCCTGCGTACATGAGCGCCTGGCGCCCTTGCTCCAAGGGGATGCGCTAGCCTGGCTGGAGATGCGTACTGCACCGCTGTAACGCAGCGGAAAAAGGGCAGTTTGCGAACTGCCCTTTTCTCTATTTACAGATGATGATc < 2:1285986/150‑1 (MQ=255) cTACCACGCCTGCGTACATGAGCGCCTGGCGCCCTTGCTCCAAGGGGATGCGCTAGCCTGGCTGGAGATGCGTACTGCACCGCTGTAACGCAGCGGAAAAAGGGCAGTTTGCGAACTGCCCTTTTCTCTATTTACAGATGATGATCATGc > 1:1740291/1‑150 (MQ=255) ccACGCCTGCGTACATGAGCGCCTGGCGCCCTTGCTCCAAGGGGATGCGCTAGCCTGGCTGGAGATGCGTACTGCACCGCTGTAACGCAGCGGAAAAAGGGCAGTTTGCGAACTGCCCTTTTCTCTATTTACAGATGATGATCATGCTcc > 2:1307820/1‑150 (MQ=255) cACGCCTGCGTACATGAGCGCCTGGCGCCCTTGCTCCAAGGGGATGCGCTAGCCTGGCTGGAGATGCGTACTGCACCGCTGTAACGCAGCGGAAAAAGGGCAGTTTGCGAACTGCCCTTTTCTCTATTTACAGATGATGATc < 1:578336/142‑1 (MQ=255) cACGCCTGCGTACATGAGCGCCTGGCGCCCTTGCTCCAAGGGGATGCGCTAGCCTGGCTGGAGATGCGTACTGCACCGCTGTAACGCAGCGGAAAAAGGGCAGTTTGCGAACTGCCCTTTTCTCTATTTACAGATGATGATc > 2:578336/1‑142 (MQ=255) aCGCCTGCGTACATGAGCGCCTGGCGCCCTTGCTCCAAGGGGATATGCTAGCCTGGCTGGAGATGCGTACTGCACCGCTGTAACGCAGCGGAAAAAGGGCAGTTTGCGAACTGCCCTTTTCTCTATTTACAGATGATGATCATGCTCCgg > 1:452050/1‑150 (MQ=255) gAGCGCCTGGCGCCCTTGCTCCAAGGGGATGCGCTAGCCTGGCTGGAGATGCGTACTGCACCGCTGTAACGCAGCGGAAAAAGGGCAGTTTGCGAACTGCCCTTTTCTCTATTTACAGATGATGATCATGCTCCGGCTGGTGTAGCccgc > 1:672791/1‑150 (MQ=255) cgcCTGGCGCCCTTGCTCCAAGGGGATGCGCTAGCCTGGCTGGAGATGCGTACTGCACCGCTGTAACGCAGCGGAAAAAGGGCAGTTTGCGAACTGCCCTTTTCTCTATTTACAGATGATGATCATGCTCCGGCTGGTGTAGCCCGCCgg < 1:828577/150‑1 (MQ=255) cTGGCGCCCTTGCTCCAAGGGGATGCGCTAGCCTGGCTGGAGATGCGTACTGCACCGCTGTAACGCAGCGGAAAAAGGGCAGTTTGCGAACTGCCCTTTTCTCTATTTACAGATGATGATCATGCTCCGGCTGGTGTAGCCCGCCGGGtt > 1:412771/1‑150 (MQ=255) gcgcCCTTGCTCCAAGGGGATGCGCTAGCCTGGCTGGAGATGCGTACTGCACCGCTGTAACGCAGCGGAAAAAGGGCAGTTTGCGAACTGTCCTTTTCTCTATTTTCAGATTATGATCATGCTCCGGCTGGTGTAGCCCGCCGGGTTGag < 1:1307820/150‑2 (MQ=255) gcgcCCTTGCTCCAAGGGGATGCGCTAGCCTGGCTGGAGATGCGTACTGCACCGCTGTAACGCAGCGGAAAAAGGGCAGTTTGCGAACTGCCCTTTTCTCTATTTACAGATGATGATCATGCTCCGGCTGGTGTAGCCCGCCGGGTTGAt < 1:56529/150‑1 (MQ=255) cgcCCTTGCTCCAAGGGGATGCGCTAGCCTGGCTGGAGATGCGTACTGCACCGCTGTAACGCAGCGGAAAAAGGGCAGTTTGCGAACTGCCCTTTTCTCTATTTACAGATGATGATCATGCTCCGGCTGGTGTAGCCCGCCGGGTTGATg < 1:1117903/150‑1 (MQ=255) gcCCTTGCTCCAAGGGGATGCGCTAGCATGGCTGGAGATGCGTACTGCACCGCTGTAACGCAGCGGAAAAAGGGCAGTTTGCGAACTGCCCTTTTCTCTATTTACAGATGATGATCATGCTCCGGCTGGTGTAGCCCGCCGGGTTGATGc > 1:1843191/1‑150 (MQ=255) cccTTGCTCCAAGGGGATGCGCTAGCCTGGCTGGAGATGCGTACTGCACCGCTGTAACGCAGCGGAAAAAGGGCAGTTTGCGAACTGCCCTTTTCTCTAGTTACAGATGATGATCATGCTCCGGCTGGTGTAGCCCGCCGGGTTGATGcc < 1:1321147/150‑1 (MQ=255) ccTTGCTCCAAGGGGATGCGCTAGCCTGGCTGGAGATGCGTACTGCACCGCTGTAACGCAGCGGAAAAAGGGCAGTTTGCGAACTGCCCTTTTCTCTATTTACAGATGATGATCATGCTCCGGCTGGTGTAGCCCGCCGGGTTGATGCCg > 1:275622/1‑150 (MQ=255) ttGCTCCAACGGGATGCGCTAGCCTGGCTGGAGGTGCGTACTGCACCGCTGTAACGCAGCGGAAAAAGGGCAGTTTGCGAACTGCCCTTTTCTCTATTTACNGATGATGATCATGCTCCGGCTGGTGTAGCCCGCCGGGTTGATGCCGaa < 1:1317833/150‑1 (MQ=255) tGCTCCAAGGGGATGCGCTAGCCTGGCTGGAGATGCGTACTGCACCGCTGTAACGCAGCGGAAAAAGGGCAGTTTGCGAACTGCCCTTTTCTCTATTTACAGATGATGATCATGCTCCGGCTGGTGTAGCCCGCCGGGTTGATGCCGAAc < 2:869150/150‑1 (MQ=255) cTCCAAGGGGATGCGCTAGCCTGGCTGGAGATGCGTACTGCACCGCTGTAACGCAGCGGAAAAAGGGCAGTTTGCGAACTGCCCTTTTCTCTATTTACAGATGATGATCATGCTCCGGCTGGTGTAGCCCGCCGGGTTGATGCCGAACaa > 1:791464/1‑150 (MQ=255) cAAGGGGATGCGCTAGCCTGGCTGGAGATGCGTACTGCACCGCTGTAACGCAGCGGAAAAAGGGCAGTTTGCGAACTGCCCTTTTCTCTATTTACTGATGATGATCATGCTCCGGCTGGTGTAGCCCGCCGGGTTGATGGTGAACAAATa > 1:1751053/1‑150 (MQ=255) cAAGGGGATGCGCTAGCCTGGCTGGAGATGCGTACTGCACCGCTGTAACGCAGCGGAAAAAGGGCAGTTTGCGAACTGCCCTTTTCTCTATTTACAGATGATGATCATGCTCCGGCTGGTGTAGCCCGCCGGGTTGATGCCGAACAAATa < 1:221405/150‑1 (MQ=255) ggggATGCGCTAGCCTGGCTGGAGATGCGTACTGCACCGCTGTAACGCAGCGGAAAAAGGGCAGTTTGCGAACTGCCCTTTTCTCTATTTACAGATGATGATCATGCTCCGGCTGGTGTAGCCCGCCGGGTTGATGCCGAACAAATAGTc > 2:1760065/1‑150 (MQ=255) ggATGCGCTAGCCTGGCTGGAGATGCGTACTGCACCTCTGTAACGCAGCGGAAAAAGGGCAGTTTGCGAACTGCCCTTTTCTCTATTTACAGATGATGATCATGCTCCGGCTGGTGTAGCCCGCCGGGTTGATGCCGAACAAATAGTCAc < 2:1707793/150‑1 (MQ=255) tggctggAGATGCGTACTGCACCGCTGTAACGCAGCGGAAAAAGGGCAGTTTGCGAACTGCCCTTTTCTCTATTTACAGATGATGATCATGCTCCGGCTGGTGTAGCCCGCCGGGTTGATGCCGAACAAATAGTCACCGGCCTCTTCATc < 2:412771/150‑1 (MQ=255) tggctggAGATGCGTACTGCACCGCTGTAACGCAGCGGAAAAAGGGCAGTTTGCGAACTGCCCTTTTCTCTATTTACAGATGATGATCATGCTCCGGCTGGTGTAGCCCGCCGGGTTGATGCCGAACAAATAGTCACCGGCCTCTTCATc > 1:486211/1‑150 (MQ=255) tggAGATGCGTACTGCACCGCTGTCACGCAGCGGAAAAAGGGCAGTTTGCGAAATGCGCTTTTCTCTATTTTTAGATGATGATCATGCTCCGGCTGGTGTAGCCCGCCGGGTTGATGCCGAACAAATAGTCACCGGCCTCTTCATCGCTa > 2:813240/1‑150 (MQ=255) tggAGATGCGTACTGCACCGCTGTAACGCAGCGGAAAAAGGGCAGTTTGCGAACTGCCCTTTTCTCTATTTACAGATGATGATCATGCTCCGGCTGGTGTAGCCCGCCGGGTTGATGCCGAACAAATAGTCACCGGCCTCTTCATCGCTa > 1:1275382/1‑150 (MQ=255) gaTGCGTACTGCACCGCTGTAACGCAGCGGAAAAAGGGCAGTTTGCGAACTGCCCTTTTCTCTATTTACAGATGATGATCATGCTCCGGCTGGTGTAGCCCGCCGGGTTGATGCCGAACAAATAGTCACCGGCCTCTTCATCGCTATCGc > 2:45094/1‑150 (MQ=255) cGTACTGCACCGCTGTAACGCAGCGGAAAAAGGGCAGTTTGCGAACTGCCCTTTTCTCTATTTACAGATGATGATCATGCTCCGGCTGGTGTAGCCCGCCGGGTTGATGCCGAACAAATAGTCACCGGCCTCTTCATCGCTATCGCCCGa < 2:486211/150‑1 (MQ=255) tACTGCACCGCTGTAACGCAGCGGAAAAAGGGCAGTTTGCGAACTGCCCTTTTCTCTATTTACAGATGATGATCATGCTCCGGCTGGTGTAGCCCGCCGGGTTGATGCCGAACAAATAGTCACCGGc > 1:699217/1‑127 (MQ=255) tACTGCACCGCTGTAACGCAGCGGAAAAAGGGCAGTTTGCGAACTGCCCTTTTCTCTATTTACAGATGATGATCATGCTCCGGCTGGTGTAGCCCGCCGGGTTGATGCCGAACAAATAGTCACCGGc < 2:699217/127‑1 (MQ=255) tACTGCACCGCTGTAACGCAGCGGAAAAAGGGCAGTTTGCGAACTGCCCTTTTCTCTATTTACAGATGATGATCATGCTCCGGCTGGTGTAGCCCGCCGGGTTGATGCCGAACAAATAGTCACCGGCCTCTTCATCGCTATCGCCCGAgc > 1:1224436/1‑150 (MQ=255) cACCGCTGTAACGCAGCGGAAAAAGGGCAGTTTGCGAACTGCCCTTTtct < 1:900817/50‑1 (MQ=255) cACCGCTGTAACGCAGCGGAAAAAGGGCAGTTTGCGAACTGCCCTTTtct > 2:900817/1‑50 (MQ=255) cACCGCTGTAACGCAGCGGAAAAAGGGCAGTTTGCGAACTGCCCTTTTCTCTATTTACAgatgatgat > 2:372575/1‑68 (MQ=255) cACCGCTGTAACGCAGCGGAAAAAGGGCAGTTTGCGAACTGCCCTTTTCTCTATTTACAgatgatgat < 1:372575/68‑1 (MQ=255) cGCTGTAACGCAGCGGAAAAAGGGCAGTTTGCGAACTGCCCTTTTCTCTATTTACAGATGATGATCATGCTCCGGCTGGTGTAGCCCGCCGGGTTGATGCCGAACAAATAGTCACCGGCCTCTTCATCGCTATCGCCCGAGCGCGCAATg > 2:944351/1‑150 (MQ=255) gCTGTAACGCAGCGGAAAAAGGGCAGTTTGCGAACTGCCCTTTTCTCTATTTACAGATGATGATCATGCTCCGGCTGGTGTAGCCCGCCGGGTTGATGCCGAACAAATAGTCACCGGCCTCTTCATCGCTATCGCCCGAGCGCGCAATGa < 2:672791/150‑1 (MQ=255) aGCGGAAAAAGGGCAGTTTGCGAACTGCCCTTTTCTCTATTTACAGATGATGATCATGCTCCGGCTGGTGTAGCCCGCCGGGTTGATGCCGAACAAATAGTCACCGGCCTCTTCATCGCTATCGCCCGCGCGCGCAATGACCCGCTAGcc > 2:43295/1‑150 (MQ=255) ctgAAAAGGGGCCTTTTGCGAACTGCCCTTTTCTCTCTTTACAGATGATGATCATGCTCCGCCTGGTGTCGCCCGCCGGGTTGATGCCGAACAAATAGTCGCCGGCTTCTTCATCGCTATCGTCCGAGCGCGCAATGACCCGGTAGCcgc < 1:813240/148‑1 (MQ=255) gAAAAAGGGCAGTTTGCGAACTGCCCTTTTCTCTATTTACAGATGATGATCATGCTCCGGCTGGTGTAGCCCGCCGGGTTGATGCCGAACAAATAGTCACCGGCCTCTTCATCGCTATCGCCCGAGCGCGCAATGACCCGGTAGCcgcgc < 2:1224436/150‑1 (MQ=255) gAAAAAGGGCAGTTTGCGAACTGCCCTTTTCTCTATTTACAGATGATGATCATGCTCCGGCTGGTGTAGCCCGCCGGGTTGATGCCGAACAAATAGTCACCGGCCTCTTCATCGCTATCGCCCGAGCGCGCAATGACCCGGTAGCcgcgc > 1:454574/1‑150 (MQ=255) ggCAGTTCGCAAACTGCCCTTTTCTCTATTTACAGATGATGATCATGCTCCGGCTGGTGTAG‑CCGCCGGGTTGATGCCGAACAAATAGTCACCGGCCTCTTCATCGCTATCGCCCGAgc < 2:322487/119‑1 (MQ=255) ggCAGTTCGCAAACTGCCCTTTTCTCTATTTACAGATGATGATCATGCTCCGGCTGGTGTAG‑CCGCCGGGTTGATGCCGAACAAATAGTCACCGGCCTCTTCATCGCTATCGCCCGAgc > 1:322487/1‑119 (MQ=255) gCAGTTTGCGAACTGCCCTTTTATACCCGTGGAGATGATGATCATGCTCCGGCTGGTGTAGCCCGCCGGGTTGATGCCGAACAAATAGTCACCGGCCTCTTCATCGCTATCGCCCGAGCGCGCAATGACCCGGTAGCCGCGCGTGCTGCa > 1:1507427/1‑150 (MQ=255) cAGTTTGCGAACTGCCCTTTTCTCTATTTACAGATGATGATCATGCTCCGGCTGGTGTAGCCCGCCGGGTTGATGCCGAACAAATAGTCACCGGCCTCTTCATCGCTATCGCCCGAGCGCGCAATGACCCGGTAGCCGCGCGTGCTGCAg < 1:1760065/150‑1 (MQ=255) | GACCAAGGAGGAGGTGGAGTGGCTGAACGGCTACCACGCCTGCGTACATGAGCGCCTGGCGCCCTTGCTCCAAGGGGATGCGCTAGCCTGGCTGGAGATGCGTACTGCACCGCTGTAACGCAGCGGAAAAAGGGCAGTTTGCGAACTGCCCTTTTCTCTATTTACAGATGATGATCATGCTCCGGCTGGTGTAGCCCGCCGGGTTGATGCCGAACAAATAGTCACCGGCCTCTTCATCGCTATCGCCCGAGCGCGCAATGACCCGGTAGCCGCGCGTGCTGCAG > NC_002947/2546727‑2547010 |
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: ‑ |