Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 3,416,904 | C→A | 8.3% | R127L (CGC→CTC) | PP_3029 ← | DNA cytosine methyltransferase family protein |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 3,416,904 | 0 | C | A | 8.3% | 190.0 / 4.8 | 60 | R127L (CGC→CTC) | PP_3029 | DNA cytosine methyltransferase family protein |
Reads supporting (aligned to +/- strand): ref base C (27/28); new base A (2/3); total (29/31) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.53e-01 |
GGATGATCGGCAGGCAAGCCAGTGGTGAGGCAGGCGACAAATGCCTTGAAGGTGCGACCCATTTCCGACTTGATCGGCTTGCCTTCGTCGTTCAGCGGCCCCCAGGCTTGGAACTCCTCAACGTTCTCCATCAGCAGCAGGCGTGGGCGGGTGGCAAACACCCAGTGAACCACAACCCATGGCAGCGAGCGCACTGCCTTACTGCGTGGGGCGCCGCCCTTGGCCTTGGAGAAGTGGCGGCAGTCAGGGGACGCCCACAAGATGCCCACTGGCTGGCCGCCGGTGGCTTCCAGCGG > NC_002947/3416757‑3417052 | ggATGATCGGCAGGCAAGCCAGTGGTGAGGCAGGCGACAAATGCCTTGAAGGTGCGACCCATTTCCGACTTGATCGGCTTGCCTTCGTCGTTCAGCGGCCCCCAGGCTTGGAACTCCTCAACGTTCTCCATCAGCAGCAGGCGTGGGCgg < 1:196102/150‑1 (MQ=255) gatgatCGGCAGGCAAGCCAGTGGTGAGGCAGGCGACAAATGCCTTGAAGGTGCGACCCATTTCCGACTTGATCGGCTTGCCTTCGTCGTTCAGCGGCCCCCAGGCTTGGAACTCCTCAACGTTCTCCATCAGCAGCAGGCGTGGGCggg > 2:531213/1‑150 (MQ=255) atgatCGGCAGGCAAGCCAGTGGTGAGGCAGGCGACAAATGCCTTGAAGGTGCGACCCATTTCCGACTTGATCGGCTTGCCTTCGTCGTTCAGCGGCCCCCAGGCTTGGAACTCCTCAACGTTCTCCATCAGCAGCAGGCGTGGGCGGGt < 2:555465/150‑1 (MQ=255) atgatCGGCAGGCAAGCCAGTGGTGAGGCAGGCGACAAATGCCTTGAAGGTGCGACCCATTTCCGACTTGATCGGCTTGCCTTCGTCGTTCAGCGGCCCCCAGGCTTGGAACTCCTCAACGTTCTCCATCAGCAGCAGGCGTGGGAGGGt < 2:1338460/150‑1 (MQ=255) tgatCGGCAGGCAAGCCAGTGGTGAGGCAGGCGACAAATGCCTTGAAGGTGCGACCCATTTCCGACTGGATCGGCTTGCCTTCGTCGTTCAGCGGCCCCCAGGCTTGGACCTCCTCAACGTTCTCCATCAGCAGCAGGCGTGGGCGGGTg < 2:47628/150‑1 (MQ=255) ggcaggcaAGCCAGTGGTGAGGCAGGCGACAAATGCCTTGAAGGTGCGACCCATTTCCGACTTGATCGGCTTGCCTTCGTCGTTCAGCGGCCCCCAGGCTTGGAACTCCTCAACGTTCTCCATCAGCAGCAGGCGTGGGCGGGTGGCaaa > 1:474788/1‑150 (MQ=255) ggcaggcaAGCCAGTGGTGAGGCAGGCGACAAATGCCTTGAAGGTGCGACCCATTTCCGACTTGATCGGCTTGCCTTCGTCGTTCAGCGGCCCCCAGGCTTGGAACTCCTCAACGTTCTCCATCAGCAGCAGGCGTGGGCGGGTGGCaaa < 2:215203/150‑1 (MQ=255) caggcGACAAATGCCTTGAAGGTGCGACCCATTTCCGACTTGATCGGCTTGCCTTCGTCGTTCAGCGGCCCCCAGGCTTGGAACTCCTCAACGTTCTCCATCAGCAGCAGGCGTGGGCGGGTGGCAAACACCCAGTGAACCACAACCCAt < 2:2174176/150‑1 (MQ=255) aaaTGCCTTGAAGGTGCGACCCATTTCCGACTTGATCGGCTTGCCTTCGTCGTTCAGCGGCCCCCAGGCTTGGAACTCCTCAACGTTCTCCATCAGCAGCAGGCGTGGGCGGGTGGCAAaca < 2:1896178/122‑1 (MQ=255) aaaTGCCTTGAAGGTGCGACCCATTTCCGACTTGATCGGCTTGCCTTCGTCGTTCAGCGGCCCCCAGGCTTGGAACTCCTCAACGTTCTCCATCAGCAGCAGGCGTGGGCGGGTGGCAAaca > 1:1896178/1‑122 (MQ=255) aTGCCTTGAAGGTGCGACCCATTTCCGACTTGATCCGCTTGCCTTCGTCGTTCAGCGGCCCCCAGGCTTGGAACTCCTCAACGTTCTCCATCAGCAGCAGGCGTGGGCGGGTGGCAAACACCCAGTGAACCACAACCCATGGCagcgagc > 1:564530/1‑150 (MQ=255) tGCCTTGAAGGTGCGACCCATTTCCGACTTGCTCGGACGGCCTTCGTCGTTCAGCGGCCCCCAGGCTTGGAACTCCTCAACGTTCTCCATCAGCAGCAGGCGTGGGCGGGTGGCAAACACCCAGTGAACCACAACCCATGGCagcgagcg < 2:564530/150‑1 (MQ=255) gAAGGTGCGACCCATTTCCGACTTGATCGGCTTGCCTTCGTCGTTCAGCGGCCCCCAGGCTTGGAACTCCTCAACGTTCTCCATCAGCAGCAGGCGTGGGCGGGTGGCAAACACCCAGTGAACCACAACCCATgg < 2:2305825/135‑1 (MQ=255) gAAGGTGCGACCCATTTCCGACTTGATCGGCTTGCCTTCGTCGTTCAGCGGCCCCCAGGCTTGGAACTCCTCAACGTTCTCCATCAGCAGCAGGCGTGGGCGGGTGGCAAACACCCAGTGAACCACAACCCATgg > 1:2305825/1‑135 (MQ=255) gTGCGACCCATTTCCGACTTGATCGGCTTGCCTTCGTCGTTCAGCGGCCCCCAGGCTTGGAACTCCTCAACGTTCTCCATCAGCAGCAGGCGTGGGCGGGTGGCAAACACCCAGTGAACCACAACCCATGGCAGTGCGCTCGCTGccatg < 1:481022/150‑4 (MQ=255) gACCCATTTCCGACTTGATCGGCTTGCCTTCGTCGTTCAGCGGCCCCCAGGCTTGGAACTCCTCAACGTTCTCCATCAGCAGCAGGCGTGGGCGGGTGGCAAACACCCAGTGAACCACAACCCATGGCAGCGAGCGCACTGCCTTACTGc > 1:1912719/1‑150 (MQ=255) ccATTTCCGACTTGATCGGCTTGCCTTCGTCGTTCAGCGGCCCCCAGGCTTGGAACTCCTCAACGTTCTCCATCAGCAGCAGGCGTGGGCGGGTGGCAAACACCCAGTGAACCACAACCCATGGCAGCGAGCGCACTGCCTTACTGCGTg > 1:2058267/1‑150 (MQ=255) cGACTTGATCGGCTTGCCTTCGTCGTTCAGCGGCCCCCAGGCTTGGAACTCCTCAACGTTCTCCATCAGCAGCAGGCGTGGGCGGGTGGCAAACACCCAGTGAACCACAACCCATGGCAGCGAGCGCACTGCCTTACTGCGTGGGGcgcc < 1:1298693/150‑1 (MQ=255) gACTTGATCGGCTTGCCTTCGTCGTTCAGCGGCCCCCAGGCTTGGAACTCCTCAACGTTCTCCATCAGCAGCAGGCGTGGGCGGGTGGCAAACACCCAg > 2:871196/1‑99 (MQ=255) gACTTGATCGGCTTGCCTTCGTCGTTCAGCGGCCCCCAGGCTTGGAACTCCTCAACGTTCTCCATCAGCAGCAGGCGTGGGCGGGTGGCAAACACCCAg < 1:871196/99‑1 (MQ=255) aCTTGATCGGCTTGCCGGCGTCGTTCAGCGGCCCCCAGGCTTGGAACTCCTCAACGTTCTCCATCAGCAGCAGGCGTGGGCGGGTGGCAAACACCCAGTGAACCACAACCCATGGCAGCGAGCGCACTGCCTTACTGCGTGGGGcgccgc > 1:683181/1‑150 (MQ=255) ttGATCGGCTTGCCTTCGTCGTTCAGCGGCCCCCAGGCTTGGAACTCCTCAACGTTCTCCATCAGCAGCAGGCGTGGGCGGGTGGCAAACACCCAGTGAACCACAACCCATGGCAGCGAGCGCACTGCCTTACTGCGTGGGGCGCCGccc > 1:2413765/1‑150 (MQ=255) aTCGGCTCGCCTTCGTCGTTCAGCGGCCCCCAGGCTTGGAACTCCTCAACGTTCTCCATCAGCAGCAGGCGTGGGCGGGTGGCAAACACCCAGTGAACCACAACCCATGGCAGCGAGCGCACTGCCTTACTGCGTGGGGCGCCGCccttg < 1:275779/150‑1 (MQ=255) tCGGCTTGCCTTCGTCGTTCAGCGGCCCCCAGGCTTGGAACTCCTCAACGTTCTCCATCAGCAGCAGGCGTGGGCGGGTGGCAAACACCCAGTGAACCACAACCCATGGCAGCGAGCGCACTGCCTTACTGCGTGGGGCGCCGCccttgg > 1:592307/1‑150 (MQ=255) cTTGCCTTCGTCGTTCAGCGGCCCCCAGGCTTGGAACTCCTCAACGTTCTCCATCAGCAGCAGGCGTGGGCGGGTGGCAAACACCCAGTGAACCACAACCCATGGCAGCGAGCGCACTGCCTTACTGCGTggg > 1:2429881/1‑133 (MQ=255) cTTGCCTTCGTCGTTCAGCGGCCCCCAGGCTTGGAACTCCTCAACGTTCTCCATCAGCAGCAGGCGTGGGCGGGTGGCAAACACCCAGTGAACCACAACCCATGGCAGCGAGCGCACTGCCTTACTGCGTggg < 2:2429881/133‑1 (MQ=255) ttGCCTTCGTCGTTCAGCGGCCCCCAGGCTTGGAACTCCTCAACGTTCTCCATCAGCAGCAGGCGTGGGCGGGTGGCAAACACCCAGTGaa > 1:2222076/1‑91 (MQ=255) ttGCCTTCGTCGTTCAGCGGCCCCCAGGCTTGGAACTCCTCAACGTTCTCCATCAGCAGCAGGCGTGGGCGGGTGGCAAACACCCAGTGaa < 2:2222076/91‑1 (MQ=255) ccTTCGTCGTTCAGCGGCCCCCAGGCTTGGAACTCCTCAACGTTCTCCATCAGCAGCAGGCGTGGGAGGGTGGCAAACACCCAGTGAACca > 1:235978/1‑91 (MQ=255) ccTTCGTCGTTCAGCGGCCCCCAGGCTTGGAACTCCTCAACGTTCTCCATCAGCAGCAGGCGTGGGAGGGTGGCAAACACCCAGTGAACca < 2:235978/91‑1 (MQ=255) cTTCGTCGTTCAGCGGCCCCCAGGCTTGGAACTCCTCAACGTTCTCCATCAGCAGCAGGCGTGGGAGGGTGGCAAACACCCAGTGAACCACAACCCATGGCAGCGAGCGCACTGCCTTACTGCGTGGGGCGCCGCCCTTGGCCTTGgaga > 1:1988709/1‑150 (MQ=255) gtcgtTCAGCGGCCCCCAGGCTTGGAACTCCTCAACGTTCTCCATCAGCAGCAGGCGTGGGCGGGTGGCAAACACCCAGTGAACCACAACCCATGGCAGCGAGCGCACTGCCTTACTGCGTGGGGCGCCGCCCTTGGCCTTGGAGAAg < 1:1392252/148‑1 (MQ=255) gtcgtTCAGCGGCCCCCAGGCTTGGAACTCCTCAACGTTCTCCATCAGCAGCAGGCGTGGGCGGGTGGCAAACACCCAGTGAACCACAACCCATGGCAGCGAGCGCACTGCCTTACTGCGTGGGGCGCCGCCCTTGGCCTTGGAGAAg > 2:1392252/1‑148 (MQ=255) cccccAGGCTTGGAACTCCTCAACGTTCTCCATCAGCAGCAGGCGTGGGCGGGTGGCAAACACCCAGTGAACCACAACCCATGGCAGCGAGCGCACTGCCTTACTGCGTGGGGCGCCGCCCTTGGCCTTGGAGAAGTGGCGGCAGTCAgg < 1:37566/150‑1 (MQ=255) cccccAGGCTTGGAACTCCTCAACGTTCTCCATCAGCAGCAGGCGTGGGCGGGTGGCAAACACCCAGTGAACCACAACCCATGGCAGCGAGCGCACTGCCTTACTGCGTGGGGCGCCGCCCTTGGCCTTGGAGAAGTGGCGGCAGTCAgg > 2:1613998/1‑150 (MQ=255) ccccAGGCTTGGAACTCCTCAACGTTCTCCATCAGCAGCAGGCGTGGGCGGGTGGCAAACACCCAGTGAACCACAACCCATGGCAGCGAGCGCACTGCCTTACTGCGTGGGGCGCCGGCCTTGGCCTTGGAGAAGTGGCGGCAGTCCggg < 1:83415/150‑1 (MQ=255) ccccAGGCTTGGAACTCCTCAACGTTCTCCATCAGCAGCAGGCGTGGGCGGGTGGCAAACACCCAGTGAACCACAACCCATGGCAGCGAGCGCACTGCCTTACTGCGTGGGGCGCCGCCCTTGGCCTTGGAGAAGt > 1:1716221/1‑136 (MQ=255) ccccAGGCTTGGAACTCCTCAACGTTCTCCATCAGCAGCAGGCGTGGGCGGGTGGCAAACACCCAGTGAACCACAACCCATGGCAGCGAGCGCACTGCCTTACTGCGTGGGGCGCCGCCCTTGGCCTTGGAGAAGt < 2:1716221/136‑1 (MQ=255) ccAGGCTTGGAACTCCTCAACGTTCTCCATCAGCAGCAGGCGTGGGCGGGTGGCAAACACCCAGTGAACCACAACCCATGGCAGCGAGCGCACTGCCTTACTGCGTGGGGCGCCGCCCTTGGCCTTGGAGAAGTGGCGGCAg < 2:1679484/142‑1 (MQ=255) ccAGGCTTGGAACTCCTCAACGTTCTCCATCAGCAGCAGGCGTGGGCGGGTGGCAAACACCCAGTGAACCACAACCCATGGCAGCGAGCGCACTGCCTTACTGCGTGGGGCGCCGCCCTTGGCCTTGGAGAAGTGGCGGCAg > 1:1679484/1‑142 (MQ=255) cAGGCTTGGAACTCCTCAACGTTCTCCATCAGCAGCAGGCGTGGGCGGGTGGCAAACACCCAGTGAACCACAACCCATGGCAGCGAGCGCACTGCCTTACTGCGTGGGGCGCCGCCCTTGGCCTTGGAGAAGTGGCGGCAGTCAGGGGAc > 1:1513229/1‑150 (MQ=255) ttGGAACTCCTCAACGTTCTCCATCAGCAGCAGGCGTGGGCGGGTGGCAAACACCCAGTGAACCACAACCCATGGCAGCGAGCGCACTGCCTTACTGCGTGGGGCGCCGCCCTTGGCCTTGGAGAAGTGGCGGCAGTCAGGGGACGCCca > 1:1610435/1‑150 (MQ=255) ggAACTCCTCAACGTTCTCCATCAGCAGCAGGCGTGGGCGGGTGGCAAACACCCAGTGAACCACAACCCATGGCAGCGAGCGCACTGCCTTACTGCGTGGGGCGCCGCCCTTGGCCTTGGAGAAGTGGCGGCAGTCAGGGGACGCCcaca > 2:2076695/1‑150 (MQ=255) tcctcAACGTTCTCCATCAGCAGCAGGCGTGGGCGGGTGGCAAACACCCAGTGAACCACAACCCATGGCAGCGAGCGCACTGCCTTACTGCGTGGGGCGCCGCCCTTGGCCTTGGAGAAGTGGCGGCAGTCAGGGGACGCCCACAAGATg > 2:1994355/1‑150 (MQ=255) ctcAACGTTCTCCATCAGCAGCAGGCGTGGGCGGGTGGCAAACACCCAGTGAACCACAACCCATGGCAGCGAGCGCACTGCCTTACTGCGTGGGGCGCCGCCCTTGGCCTTGGAGAAGTGGCGGCAGTCAGGGGACGCCCACAAGATGcc > 1:1214811/1‑150 (MQ=255) cAACGTTCTCCATCAGCAGCAGGCGTGGGCGGGTGGCAAACACCCAGTGAACCACAACCCATGGCAGCGAGCGCACTGCCTTACTGCGTGGGGCGCCGCCCTTGGCCTTGGAGAAGTGGCGGCAGTCAGGGGACGCCCACAAGATGCCCa < 2:2058267/150‑1 (MQ=255) cAACGTCCGCCATCAGCAGCAGGCGTGGGCGGGTGGCAAACACCCAGTGAACCACAACCCATGGCAGCGAGCGCACTGCCTTACTGCGTGGGGCGCCGCCCTTGGCCTTGGAGAAGTGGCGGCAGTCAGGGGACGCCCACAAGATGCCCa < 1:339325/150‑1 (MQ=255) aaCGTTCTCCATCAGCAGCAGGCGTGGGCGGGTGGCAAACACCCAGTGAACCACAACCCATGGCAGCGAGCGCACTGCCTTACTGCGTGGGGCGCCGCCCTTGGCCTTGGAGAAGTgg < 2:1033724/118‑1 (MQ=255) aaCGTTCTCCATCAGCAGCAGGCGTGGGCGGGTGGCAAACACCCAGTGAACCACAACCCATGGCAGCGAGCGCACTGCCTTACTGCGTGGGGCGCCGCCCTTGGCCTTGGAGAAGTgg > 1:1033724/1‑118 (MQ=255) cGTTCTCCATCAGCAGCAGGCGTGGGCGGGTGGCAAACACCCAGTGAACCACAACCCATGGCAGCGAGCGCACTGCGTTACTGCGTGGGGCGCCGCCCTTGGCCTTGGAGAAGTGGCGGCAGTCAGGGGACGCCCACAAGATGCCCActg < 1:975201/150‑1 (MQ=255) ttCTCCATCAGCAGCAGGCGTGGGCGGGTGGCAAACACCCAGTGAACCACAACCCATGGCAGCGAGCGCACTGCCTTACTGCGTGGGGCGCCGCCCTTGGCCTTGGAGAAGTGGCGGCAGTCAGGGGACGCCCACAAGATGCCCActggc > 1:226394/1‑150 (MQ=255) ttCTCCATCAGCAGCAGGCGTGGGCGGGTGGCAAACACCCAGTGAACCACAACCCATGGCAGCGAGCGCACTGCCTTACTGCGTGGGGCGCCGCCCTTGGCCTTGGAGAAGTGGCGGCAGTCAGGGGACGCCCACAAGATGCCCActggc > 1:1469784/1‑150 (MQ=255) ttCTCCATCAGCAGCAGGCGTGGGCGGGTGGCAAACACCCAGTGAACCACAACCCATGGCAGCGAGCGCACTGCCTTACTGCGTGGGGCGCCGCCCTTGGCCTTGGAGAAGTGGCGGCAGTCAGGGGACGCCCACAAGATGCCCActggc > 1:1290302/1‑150 (MQ=255) ttCTCCATCAGCAGCAGGCGTGGGAGGGTGGCAAACACCCAGTGAACCACAACCCATGGCAGCGAGCGCACTGCCTTACTGCGTGGGGCGCCGCCCTTGGCCTTGGAGAAGTGGCGGCAGTCAGGGGACGCCCACAAGATGCCCActggc < 2:1988709/150‑1 (MQ=255) ctcCATCAGCAGCAGGCGTGGGCGGGTGGCAAACACCCAGTGAACCACAACCCATGGCAGCGAGCGCACTGCCTTACTGCGTGGGGCGCCGCCCTTGGCCTTGGAGAAGTGGCGGCAGTCAGGGGACGCCCACAAGATGCCCActggctg < 1:531213/150‑1 (MQ=255) ccATCAGCAGCAGGCGTGGGCGGGTGGCAAACACCCAGTGAACCACAACCCATGGCAGCGAGCGCACTGCCTTACTGCGTGGGGCGCCGCCCTTGGCCTTGGAGAAGTGGCGGCAGTTAGGGGACGCCCACAAGATGCCCACtggctggc < 2:1349822/150‑1 (MQ=255) cATCAGCAGCAGGCGTGGGCGGGTGGCAAACACCCAGTGAACCACAACCCATGGCAGCGAGCGCACTGCCTTACTGCGTGGGGCGCCGCCCTTGGCCTTGGAGAAGTGGCGGCAGTCAGGGGACGCCCACAAGATGCCCACTGGCTGgcc > 2:2306303/1‑150 (MQ=255) agcagGCGTGGGCGGGTGGCAAACACCCAGTGAACCACAACCCATGGCAGCGAGCGCACTACCTTAGACGATTCGGCGCCGCCCTTGGCCTTGGAGAAGTGGCGGCAGTCAGGGGACGCCCACAAGATGCCCACTGGCTGGCCGCCGGTg < 1:1614451/150‑1 (MQ=255) gCGGGTGGCAAACACCCAGTGAACCACAACCCATGGCAGCGAGCGCACTGCCTTACTGCGTGGGGCGCCGCCCTTGGCCTTGGAGAAGTGGCGGCAGTCAGGGGACGCCCACAAGATGCCCACTGGCTGGCCGCCGGTGGCTTCCAGCgg < 1:1753585/150‑1 (MQ=255) gCGGGTGGCAAACACCCAGTGAACCACAACCCATGGCAGCGAGCGCACTGCCTTACTGCGTGGGGCGCCGCCCTTGGCCTTGGAGAAGTGGCGGCAGTCAGGGGACGCCCACAAGATGCCCACTGGCTGGCCGCAGGTGGCTTCGAGCgg < 2:592307/150‑1 (MQ=255) | GGATGATCGGCAGGCAAGCCAGTGGTGAGGCAGGCGACAAATGCCTTGAAGGTGCGACCCATTTCCGACTTGATCGGCTTGCCTTCGTCGTTCAGCGGCCCCCAGGCTTGGAACTCCTCAACGTTCTCCATCAGCAGCAGGCGTGGGCGGGTGGCAAACACCCAGTGAACCACAACCCATGGCAGCGAGCGCACTGCCTTACTGCGTGGGGCGCCGCCCTTGGCCTTGGAGAAGTGGCGGCAGTCAGGGGACGCCCACAAGATGCCCACTGGCTGGCCGCCGGTGGCTTCCAGCGG > NC_002947/3416757‑3417052 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |