Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 1,354,272 | A→C | 12.8% | S371A (TCC→GCC) ‡ | nrdA ← | ribonucleoside‑diphosphate reductase subunit alpha |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 1,354,272 | 0 | A | C | 12.8% | 134.5 / 12.5 | 56 | S371A (TCC→GCC) ‡ | nrdA | ribonucleoside‑diphosphate reductase subunit alpha |
Reads supporting (aligned to +/- strand): ref base A (23/25); new base C (4/3); total (28/28) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 7.05e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.26e-01 |
GTGCCCTTGATGTAAGAGCCCAGTGCGCGCACAGGGGTCCAGTCGTTGCCCAGGCCACCGGCGAATTTCGACAGCATGGCGTTGTCGTGGATCGCGTGGTAGATGCCCGACAGGTCGTCCGGCACGGTGGTCAGGTAGCAGCTGGACAGCTGCGGGCGCAGGGTGCCGGCGTTGAACAGGGTCGGGGTCGAGGCCATGTAGTCGAACGACGACAACAGGTTGTAGAACTCGATCGCACGGGCTTCCTTGTCTTTTTCTTCCAGCGCCAGGCCCATGGCCACACGCATGA > NC_002947/1354127‑1354415 | gTGCCCTTGATGTAAGAGCCCAGTGCGCGCACAGGGGTCCAGTCGTTGCCCAGGCCACCGGCGAATTTCGACAGCATGGCGTTGTCGTGGATCGCGTGGTAGATGCCCGACAGGTCGTCCGGCACGGTGGTCAGGTAGCAGCTGGAcccc < 2:715403/150‑4 (MQ=255) tGCCCTTGATGTAAGAGCCCAGTGCGCGCACAGGGGTCCAGTCGTTGCCCAGGCCACCGGCGAATTTCGACAGCATGGCGTTGTCGTGGATCGCGTGGTAGATGCCCGACAGGTCGTCCGGCACGGTGGTCAGGTAGCAGCTGGACAGCt < 1:1133125/150‑1 (MQ=255) cccTTGATGTAAGAGCCCAGTGCGCGCACAGGGGTCCAGTCGTTGCCCAGGCCACCGGCGAATTTCGACAGCATGGCGTTGTCGTGGATCGCGTGGTAGATGCCCGACAGGTCGTCCGGCACGGTGGTCAGGTAGCAGCTGGACAGCTGc > 1:464276/1‑150 (MQ=255) ccTTGATGTAAGAGCCCAGTGCGCGCACAGGGGTCCAGTCGTTGCCCAGGCCACCGGCGAATTTCGACAGCATGGCGTTGTCGTGGATCGCGTGGTAGATGCCCGACAGGTCGTCCGGCACGGTGGTCAGGTAGCAGCTGGACAGCTGCg < 1:572612/150‑1 (MQ=255) ttGATGTAAGAGCCCAGTGCGCGCACAGGGGTCCAGTCGTTGCCCAGGCCACCTGCGAATTTCGACAGCATGGCGTTGTCGTGGATCGCGTGGTAGATGCCCGACAGGTCGTCCGGCACGGTGGTCAGGTAGCAGCTGTCCAGCTGctac < 2:1213116/150‑4 (MQ=255) ttGATGTAAGAGCCCAGTGCGCGCACAGGGGTCCAGTCGTTGCCCAGGCCACCGGCGAATTTCGACAGCATGGCGTTGTCGTGGATCGCGTGGTAGATGCCCGACAGGTCGTCCGGCACGGTGGTCAGGTAGCAGCTGGACAGCTGCggg > 1:208535/1‑150 (MQ=255) agCCCAGTGCGCGCACAGGGGTCCAGTCGTTGCCCAGGCCACCGGCGAATTTCGACAGCATGGCGTTGTCGTGGATCGCGTGGTAGATGCCCGACAGGTCGTCCGGCACGGTGGTCAGGTAGCAGCTGGACAGCTGCGGGCGCAGGGTGc > 1:221352/1‑150 (MQ=255) gCCCAGTGCGCGCACAGGGGTCCAGTCGTTGCCCAGGCCACCGGCGAATTTCGACAGCATGGCGTTGTCGTGGATCGCGTGGTAGATGCCCGACAGGTCGTCCGGCACGGTGGTCAGGTAGCAGCTGGACAGCTGCGGGCGCAGGGTGcc < 2:475234/150‑1 (MQ=255) cccAGTGCGCGCACAGGGGTCCAGTCGTTGCCCAGGCCACCGGCGAATTTCGACAGCATGGCGTTGTCGTGGATCGCGTGGTAGATGCCCGACAGGTCGTCCGGCACGGTGGTCAGGTAGCAGCTGGa > 2:632771/1‑128 (MQ=255) cccAGTGCGCGCACAGGGGTCCAGTCGTTGCCCAGGCCACCGGCGAATTTCGACAGCATGGCGTTGTCGTGGATCGCGTGGTAGATGCCCGACAGGTCGTCCGGCACGGTGGTCAGGTAGCAGCTGGa < 1:632771/128‑1 (MQ=255) cAGTGCGCGCACAGGGGTCCAGTCGTTGCCCAGGCCACCGGCGAATTTCGACAGCATGGCGTTGTCGTGGATCGCGTGGTAGATGCCCGACAGGTCGTCCGGCACGGTGGTCAGGTAGCAGCTGGACAGCTGCggg > 2:976591/1‑136 (MQ=255) cAGTGCGCGCACAGGGGTCCAGTCGTTGCCCAGGCCACCGGCGAATTTCGACAGCATGGCGTTGTCGTGGATCGCGTGGTAGATGCCCGACAGGTCGTCCGGCACGGTGGTCAGGTAGCAGCTGGACAGCTGCggg < 1:976591/136‑1 (MQ=255) gTGCGCGCACAGGGGTCCAGTCGTTGCCCAGGCCACCGGCGAATTTCGACAGCATGGCGTTGTCGTGGATCGCGTGGTAGATGCCCGACAGGTCGTCCGGCACGGTGGTCAGGTAGCAGCTGGACAGCTGCGGGCGCAGGGTGCCGGCGt < 2:192096/150‑1 (MQ=255) aGGGGTCCAGTCGTTGCCCAGGCCACCGGCGAATTTCGACAGCATGGCGTTGTCGTGGATCGCGTGGTAGATGCCCGACAGGTCGTCCGGCACGGTGGTCAGGTAGCAGCTGGACAGCTGCGGGCGCAGGGTGCCGGCGTTGAACAgggt < 2:221352/150‑1 (MQ=255) aGTCGTTGCCCAGGCCACCGGCGAATTTCGACAGCATGGCGTTGTCGTGGATCGCGTGGTAGATGCCCGACAGGTCGTCCGGCACGGTGGTCAGGTAGCAGCTGGACAGCTGCGGGCGCAGGGTGCCGGCGTTGAACAgggtcggggtcg < 1:1383159/150‑1 (MQ=255) gTCGTTGCCCAGGCCACCGGCGAATTTCGACAGCATGGCGTTGTCGTGGATCGCGTGGTAGATGCCCGACAGGTCGTCCGGCACGGTGGTCAGGTAGCAGCTGGACAGCTGCGGGCGCAGGGTGCCGGCGTTGAACAGGGTCGGGGTCGa > 2:84028/1‑150 (MQ=255) ccAGGCCACCGGCGAATTTCGACAGCATGGCGTTGTCGTGGATCGCGTGGTAGATGCCCGACAGGTCGTCCGGCACGGTGGTCAGGTAGCAGCTGGACAGCTGCGGGCGCAGGGTGCCGGCGTTGAACAGGGTCGGGGTCGAGGCCATGt > 1:391240/1‑150 (MQ=255) cACCGGCGAATTTCGACAGCATGGCGTTGTCGTGGATCGCGTGGTAGATGCCCGACAGGTCGTCCGGCACGGTGGTCAGGTAGCAGCTGGACAGCTGCGGGCGCAGGGTGCCGGCGTTGAACAGGGTCGGGGTCGAGGCCATGTAGTCGa > 2:887211/1‑150 (MQ=255) cACCGGCGAATTTCGACAGCATGGCGTTGTCGTGGATCGCGTGGTAGATGCCCGACAGGTCGTCCGGCACGGTGGTCAGGTAGCAGCTGGACAGCTGCGGGCGCAGGGTGCCGGCGTTGAACAGGGTCGGGGTCGAGGCCATGTAGTCGa > 1:118394/1‑150 (MQ=255) cGGCGAATTTCGACAGCATGGCGTTGTCGTGGATCGCGTGGTAGATGCCCGACAGGTCGTCCGGCACGGTGGTCAGGTAGCAGCTGGACAGCTGCGGGCGCAGGGTGCCGGCGTTGAACAGGGTCGGGGTCGAGGCTATGTAGTCGAacg > 2:33512/1‑150 (MQ=255) cGGCGAATTTCGACAGCATGGCGTTGTCGTGGATCGCGTGGTAGATGCCCGACAGGTCGTCCGGCACGGTGGTCAGGTAGCAGCTGGACAGCTGCGGGCGCAGGGTGCCGGCGTTGAACAGGGTCGCGGTCGAGGCCATGTAGTCGAacg > 2:277247/1‑150 (MQ=255) ggCGAATTTCGACAGCATGGCGTTGTCGTGGATCGCGTGGTAGATGCCCGACAGGTCGTCCGGCACGGTGGTCAGGTAGCAGCTGGACAGCTGCGGGCGCAGGGTGCCGGCGTTGAACAGGGTCGGGGTCGAGGCCATGTAGTCGaa < 2:929714/147‑1 (MQ=255) ggCGAATTTCGACAGCATGGCGTTGTCGTGGATCGCGTGGTAGATGCCCGACAGGTCGTCCGGCACGGTGGTCAGGTAGCAGCTGGACAGCTGCGGGCGCAGGGTGCCGGCGTTGAACAGGGTCGGGGTCGAGGCCATGTAGTCGaa > 1:929714/1‑147 (MQ=255) ggCGAATTTCGACAGCATGGCGTTGTCGTGGATCGCGTGGTAGATGCCCGACAGGTCGTCCGGCACGGTGGTCAGGTAGCAGCTGGACAGCTGCGGGCGCAGGGTGCCGGCGTTGAACAGGGTCGGGGTCGAGGCCATGTAGTCGAacga < 2:208535/150‑1 (MQ=255) gCGAATTTCGACAGCATGGCGTTGTCGTGGATCGCGTGGTAGATGCCCGACAGGTCGTCCGGCACGGTGGTCAGGTAGCAGCTGGACAGCTGCGGGCGCAGGGTGCCGGCGTTGAACAGGGTCGGGGTCGAGGCCATGTAGTCGAacgac < 2:391240/150‑1 (MQ=255) gAATTTCGACAGCATGGCGTTGTCGTGGATCGCGTGGTAGATGCCCGACAGGTCGTCCGGCACGGTGGTCAGGTAGCAGCTGGACAGCTGCGAGCGCAGGGTGCCGGCGTTGAACAGGGTCGGGGTCGAGGCCATGTAGTCGAAcgacga > 1:169027/1‑150 (MQ=255) cAGCATGGCGTTGTCGTGGATCGCGTGGTAGATGCCCGACAGGTCGTCCGGCACGGTGGTCAGGTAGCAGCTGGACAGCTGCGGGCGCAGGGTGCCGGCGTTGAACAGGGTCGGGGTCGAGGCCATGTAGTCGAACGACGACAACAGGtt < 2:464276/150‑1 (MQ=255) cATGGCGTTGTCGTGGATCGCGTGGTAGATGCCCGACAGGTCGTCCGGCACGGTGGTCAGGTAGCAGCTGGACAGCTGCGGGCGCAGGGTGCCGGCGTTGAACAGGGTCGGGGTCGAGGCCATGTAGTCGAACGACGACAACAGGTTGTa < 1:229955/150‑1 (MQ=255) cATGGCGTTGTCGTGGATCGCGTGGTAGATGCCCGACAGGTCGTCCGGCACGGTGGTCAGGTAGCAGCTGGACAGCTGCGGGCGCAGGGTGCCGGCGTTGAACAGGGTCGGGGTCGAGGCCATGTAGTCGAACGACGACAACAGGTTGTa < 1:84028/150‑1 (MQ=255) gCGTTGTCGTGGATCGCGTGGTAGATGCCCGACAGGTCGTCCGGCACGGTGGTCAGGTAGCAGCTGGACAGCTGCGGGCGCAGGGTGCCGGCGTTGAACAGGGTCGGGGTCGAGGCTATGTAGTCGAACGACGACAACAGGTTGTAGAAc < 1:33512/150‑1 (MQ=255) gCGTTGTCGTGGATCGCGTGGTAGATGCCCGACAGGTCGTCCGGCACGGTGGTCAGGTAGCAGCTGGACAGCTGCGGGCGCAGGGTGCCGGCGTTGAACAGGGTCGGGGTCGAGGCCATGTAGTCGAACGACGACAACAGGTTGTAGAAc < 1:58900/150‑1 (MQ=255) tCGTGGATCGCGTGGTAGATGCCCGACAGGTCGTCCGGCACGGTGGTCAGGTAGCAGCTGGACAGCTGCGGGCGCAGGGTGCCGGCGTTGAACAg > 2:1417327/1‑95 (MQ=255) tCGTGGATCGCGTGGTAGATGCCCGACAGGTCGTCCGGCACGGTGGTCAGGTAGCAGCTGGACAGCTGCGGGCGCAGGGTGCCGGCGTTGAACAg < 1:1417327/95‑1 (MQ=255) tCGTGGATCGCGTGGTAGATGCCCGACAGGTCGTCCGGCACGGTGGTCAGGTAGCAGCTGGACAGCTGCGGGCGCAGGGTGCCGGCGTTGAACAGGGTCGGGGTCGAGGCCATGTAGTCGAACGACGACAACAGGTTGTAGAACTCGATc > 2:854769/1‑150 (MQ=255) aTCGCGTGGTAGATGCCCGACAGGTCGTCCGGCACGGTGGTCAGGTAGCAGCTGGACAGCTGCGGGCGCAGGGTGCCGGCGTTGAACAGGGTCGGGGTCGAGGCCATGTAGTCGAACGACGACAACAGGTTGTAGAACTCGATCGCACgg > 1:503994/1‑150 (MQ=255) tCGCGTGGTAGATGCCCGACAGGTCGTCCGGCACGGTGGTCAGGTAGCAGCTGGACAGCTGCGGGCGCAGGGTGCCGGCGTTGAACAGGGTCGGGGTCGAGGCCATGTAGTCGAACGACGACAACAg < 1:121883/127‑1 (MQ=255) tCGCGTGGTAGATGCCCGACAGGTCGTCCGGCACGGTGGTCAGGTAGCAGCTGGACAGCTGCGGGCGCAGGGTGCCGGCGTTGAACAGGGTCGGGGTCGAGGCCATGTAGTCGAACGACGACAACAg > 2:121883/1‑127 (MQ=255) tCGCGTGGTAGATGCCCGACAGGTCGTCCGGCACGGTGGTCAGGTAGCAGATGGACAGCTGCGGGCGCAGGGTGCCGGCGTTGAACAGGGTCGGGGTCGAGGCCATGTAGTCGAACGACGACAACAGGTTGTAGAACTCGATCGCACggg < 1:442679/150‑1 (MQ=255) gcgTGGTAGATGCCCGACAGGTCGTCCGGCACGGTGGTCAGGTAGCAGCTGGACAGCTGCGGGCGCAGGGTGCCGGCGTTGAACAGGGTCGGGGTCGAGGCCATGTAGTCGAACGACGACAACAGGTTGTAGAAc < 1:65430/135‑1 (MQ=255) gcgTGGTAGATGCCCGACAGGTCGTCCGGCACGGTGGTCAGGTAGCAGCTGGACAGCTGCGGGCGCAGGGTGCCGGCGTTGAACAGGGTCGGGGTCGAGGCCATGTAGTCGAACGACGACAACAGGTTGTAGAAc > 2:65430/1‑135 (MQ=255) tAGATGCCCGACAGGTCGTCCGGCACGGTGGTCAGGTAGCAGCTGGACAGCTGCGGGCGCAGGGTGCCGGCGTTGAACAGGGTCGGGGTCGAGGCCATGTAGTCGAACGACGACAACAGGTTGTAGAACTCGATCGCACGGGCTTCCTTg > 2:1176838/1‑150 (MQ=255) gATGCCCGACAGGTCGTCCGGCACGGTGGTCAGGTAGCAGCTGGACAGCTGCGGGCGCAGGGTGCCGGCGTTGAACAGGGTCGGGGTCGAGGCCATGTAGTCGAACGACGACAACAGGTTGTAGAACTCGATCGCACGGGCTTCCTTGTc > 2:345851/1‑150 (MQ=255) gATGCCCGACAGGTCGTCCGGCACGGTGGTCAGGTAGCAGATGGACAGCTGCGGGCGCAGGGTGCCGGCGTTGAACAGGGTCGGGGTCGAGGCCATGTAGTCGAACGACGACAACAGGTTGTAGAACTCGATCGCACGGGCTTCCTTGTc > 1:1358045/1‑150 (MQ=255) aTGCCCGACAGGTCGTCCGGCACGGTGGTCAGGTAGCAGCTGGACAGCTGCGGGCGCAGGGTGCCGGCGTTGAACAGGGTCGGGGTCGAGGCCATGTAGTCGAACGACGACAACAGGTTGTAGAACTCGATCGCACGGGCTTCCTTGTCt < 1:854769/150‑1 (MQ=255) cGACAGGTCGTCCGGCACGGTGGTCAGGTAGCAGCTGGACAGCTGCGGGCGCAGGGTGCCGGCGTTGAACAGGGTCGGGGTCGAGGCCATGTAGTCGAACGACGACAACAGGTTGTAGAACTCGATCGCACGGGCTTCCTTTTCTTTttc > 2:546501/1‑150 (MQ=255) tcgtcCGGCACGGTGGTCAGGTAGCAGCTGGACAGCTGCGGGCGCAGGGTGCCGGCGTTGAACAGGGTCGGGGTCGAGGCCATGTAGTCGAACGACGACAACAGGTTGTAGAACTCGATCGCACGGGCTTCCTTGTCTTTTTCTTCCAgc < 1:345851/150‑1 (MQ=255) tcCGGCACGGTGGTCAGGTAGCAGCTGGACAGCTGCGGGCGCAGGGTGCCGGCGTTGAAc > 1:609981/1‑60 (MQ=255) tcCGGCACGGTGGTCAGGTAGCAGCTGGACAGCTGCGGGCGCAGGGTGCCGGCGTTGAAc < 2:609981/60‑1 (MQ=255) tcCGGCACGGTGGTCAGGTAGCAGCTGGACAGCTGCGGGCGCAGGGTGCCGGCGTTGAACAGGGTCGGGGTCGAGGCCATGTAGTCGAACGACGACAACAGGTTGTAGAACTCGATCGCACGGGCTTCCTTGTCttttt < 1:631130/139‑1 (MQ=255) tcCGGCACGGTGGTCAGGTAGCAGCTGGACAGCTGCGGGCGCAGGGTGCCGGCGTTGAACAGGGTCGGGGTCGAGGCCATGTAGTCGAACGACGACAACAGGTTGTAGAACTCGATCGCACGGGCTTCCTTGTCttttt > 2:631130/1‑139 (MQ=255) ggTAGCAGCTGGACAGCTGCGGGCGCAGGGTGCCGGCGTTGAACAGGGTCGGGGTCGAGGCCATGTAGTCGAACGACGACAACAGGTTGTAGAACTCGATCGCACGGGCTTCCTTGTCTTTTTCTTCCAGCGCCAGGCCCATGGCCacac < 1:1176838/150‑1 (MQ=255) ccgcagcTGTCCAGCTGCGGGCGCAGGGTGCCGGCGTTGAACAGGGTCGGGGTCGAGGCCATGTAGTCGAACGACGACAACAGGTTg < 2:1139687/85‑1 (MQ=255) ccgcagcTGTCCAGCTGCGGGCGCAGGGTGCCGGCGTTGAACAGGGTCGGGGTCGAGGCCATGTAGTCGAACGACGACAACAGGTTg > 1:1139687/3‑87 (MQ=255) ccgcagcTGTCCAGCTGCGGGCGCAGGGTGCCGGCGTTGAACAGGGTCGGGGTCGAGGCCATGTAGTCGAACGACGACAACAGGTTGTAGAACTCGATCGCACGGGCTTCCTTGTCTTTTTCTTCCAGCGCCAGGCCCATgg < 1:1126195/140‑1 (MQ=255) ccgcagcTGTCCAGCTGCGGGCGCAGGGTGCCGGCGTTGAACAGGGTCGGGGTCGAGGCCATGTAGTCGAACGACGACAACAGGTTGTAGAACTCGATCGCACGGGCTTCCTTGTCTTTTTCTTCCAGCGCCAGGCCCATgg > 2:1126195/3‑142 (MQ=255) ccgcagcTGTCCAGCTGCGGGCGCAGGGTGCCGGCGTTGAACAGGGTCGGGGTCGAGGCCATGTAGTCGAACGACGACAACAGGTTGTAGAACTCGATCGCACGGGCTTCCTTGTCTTTTTCTTCCAGCGCCAGGCCCATGGCCACACGc > 1:506267/3‑150 (MQ=255) gcagcTGTCCAGCTGCGGGCGCAGGGTGCCGGCGTTGAACAGGGTCGGGGTCGAGGCCATGTAGTCGAACGACGACAACAGGTTGTAGAACTCGATCGCACGGGCTTCCTTGTCTTTTTCTTCCAGCGCCAGGCCCATGGCCACACGCAt > 2:296328/1‑150 (MQ=255) agcTGGGCAGCTGCGGGCGCAGGGTGCCGGCGTTGAACAGGGTCGGGGTCGAGGCCATGTAGTCGAACGACGACAACAGGTTGTAGAACTCGATCGCACGGGCTTCCTTGTCTTTTTCTTCCAGCGCCAGGCCCATGGCCACACGCATga > 2:708609/1‑150 (MQ=255) | GTGCCCTTGATGTAAGAGCCCAGTGCGCGCACAGGGGTCCAGTCGTTGCCCAGGCCACCGGCGAATTTCGACAGCATGGCGTTGTCGTGGATCGCGTGGTAGATGCCCGACAGGTCGTCCGGCACGGTGGTCAGGTAGCAGCTGGACAGCTGCGGGCGCAGGGTGCCGGCGTTGAACAGGGTCGGGGTCGAGGCCATGTAGTCGAACGACGACAACAGGTTGTAGAACTCGATCGCACGGGCTTCCTTGTCTTTTTCTTCCAGCGCCAGGCCCATGGCCACACGCATGA > NC_002947/1354127‑1354415 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |