Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 5,255,730 | G→T | 14.3% | G71V (GGC→GTC) | folM → | bifunctional dihydrofolate reductase/dihydromonapterin reductase |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 5,255,730 | 0 | G | T | 14.3% | 69.8 / 3.8 | 28 | G71V (GGC→GTC) | folM | bifunctional dihydrofolate reductase/dihydromonapterin reductase |
Reads supporting (aligned to +/- strand): ref base G (13/11); new base T (2/2); total (15/13) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.51e-01 |
AGCCAGCCAGCGCGTCGGGCTGGCCCTGGCCCTCGAACTGGCACAGGCCGGCCATACGGTGGTCAGTGCCAGCCGCAGCATCCAGCCGCAGGCAGCCCACCCCAACATCGTGCAGTTTCAGGCCGACCTGTGCCAGGCGGGCGACCGCCAGGCGCTGATCGATCACCTGCACGATCGTTATGACGGACTACGCGCCATCATCCACAACGCCTCGCTATGGCTCGACGACGGCCTCGACAACCTCGAGACCATGTTCCGCCTGCACGTCGAAGCGCCCTACCACCTCAACC > NC_002947/5255590‑5255879 | agccagccAGCGCGTCGGGCTGGCCCTGGCCCTCGAACTGGCACAGGCCGGCCATACGGTGGTCAGTGCCAGCCGCAGCATCCAGCCGCAGGCAGTCCACCCCAACATCGTGCAGTTTCAGGCCGACCTGTGCCAGGCGGGCGACCGCCa < 1:351941/150‑1 (MQ=255) gggccagcCCGACGCGCTGGCCCTGGCCCTCGAACTGGCACAGGCCGGCCATACGGTGGTCAGTGCCAGCCGCAGCATCCAGCCGCAGGCAGCCCACCCCAACATCGTGCAGTTTCAGGCCGACCTGTGCCAGGCGGGCGACCGCCAGgc > 2:714517/3‑150 (MQ=255) gccagcGCGTCGGGCTGGCCCTGGCCCTCGAACTGGCACAGGCCGGCCATACGGTGGTCAGTGCCAGCCGCAGCATCCAGCCGCAGGCAGCCCACCCCAACATCGTGCAGTTTCAGGCCGACCTGTGCCAGGCGGGCGACCGCCAGgcgc < 1:631534/150‑1 (MQ=255) gcgTCGGGCTGGCCCTGGCCCTCGAACTGGCACAGGCCGGCCATACGGTGGTCAGTGCCAGCCGCAGCATCCAGCCGCAGGCAGCCCACCCCAACATCGTGCAGTTTCAGGCCGACCTGTGCCAGGCGGGCGACCGCCAGGCGCTgatcg < 1:617743/150‑1 (MQ=255) gggCTGGCCCTGGCCCTCGAACTGGCACAGGCCGGCCATACGGTGGTCAGTGCCAGCCGCAGCATCCAGCCGCAGGCAGCCCACCCCAACATCGTGCAGTTTCAGGCCGACCTGTGCCAGGCGGGCGACCGCCAGGCGCTGATCGATCAc < 2:2612/150‑1 (MQ=255) cctggccctCGAACTGGCACAGGCCGGCCATACGGTGGTCAGTGCCAGCCGCAGCATCCAGCCGCAGGCAGCCCACCCCAACATCGTGCAGTTTCAGGCCGACCTGTGCCAGGCGGGCGACCGCCAGGCGCTGATCGATCACCTGCCCGa > 1:72885/1‑150 (MQ=255) cTGGCACAGGCCGGCCATACGGTGGTCAGTGCCAGCCGCCGCATCCAGCCTTATTCAGCCCACCCCAACATCGTGCAGTTTCAGGCCGACCTGGGCCAGTTGGGCGTTCGCCAGGCGCTGATCGATCACCTGCACGATCGTTATGACGGa < 1:618974/150‑1 (MQ=255) gCACAGGCCGGCCATACGGTGGTCAGTGCCAGCCGCAGCATCCAGCCGCAGGCAGCCCACCCCAACATCGTGCAGTTTCAGGCCGACCTGTGCCAGGCGGGCGACCGCCAGGCGCTGATCGATCACCTGCACGATCGTTATGACGGACTa > 1:104078/1‑150 (MQ=255) aGGCCGGCCATACGGTGGTCAGTGCCAGCCGCAGCATCCAGCCGCAGGCAGCCCACCCCAACATCGTGCAGTTTCAGGCCGACCTGTGCCAGGCGGTcgcccgcctg > 2:621000/1‑105 (MQ=255) aGGCCGGCCATACGGTGGTCAGTGCCAGCCGCAGCATCCAGCCGCAGGCAGCCCACCCCAACATCGTGCAGTTTCAGGCCGACCTGTGCCAGGCGGTcgcccgcctg < 1:621000/107‑3 (MQ=255) gccggccATACGGTGGTCAGTGCCAGCCGCAGCATCCAGCCGCAGGCAGCCCACCCCAACATCGTGCAGTTTCAGGCCGACCTGTGCCAGGCGGGCGACCGCCAGGCGCTGATCGATCACCTGCACGa < 1:764459/128‑1 (MQ=255) gccggccATACGGTGGTCAGTGCCAGCCGCAGCATCCAGCCGCAGGCAGCCCACCCCAACATCGTGCAGTTTCAGGCCGACCTGTGCCAGGCGGGCGACCGCCAGGCGCTGATCGATCACCTGCACGa > 2:764459/1‑128 (MQ=255) gtggtCAGTGCCAGCCGCAGCATCCAGCCGCAGGCAGCCCACCCCAACATCGTGCAGTTTCAGGCCGACCTGTGCCAGGCGGGCGACCGCCAGGCGCTGATCGATCACCTGCACGATCGTTATGACGGACTACGCGCCATCATCCACAAc < 2:682755/150‑1 (MQ=255) gtCAGTGCCAGCCGCAGCATCCAGCAGCGTGCAGCCCACCCCAACATCGTGCAGTTTCAGGCCGACCTGTGCCAGGCGGGCGACCGCCAGGCGCTGATCGATCACCTGCACGATCGTTATGACGGACTACGCGCCATCATCCACAACGcc > 1:743874/1‑150 (MQ=255) tCAGTGCCAGCCGCAGCATCCAGCCGCAGGCAGCCCACCCCAACATCGTGCAGTTTCAGGCCGACCTGTGCCAggcgggc < 1:479315/80‑1 (MQ=255) tCAGTGCCAGCCGCAGCATCCAGCCGCAGGCAGCCCACCCCAACATCGTGCAGTTTCAGGCCGACCTGTGCCAggcgggc > 2:479315/1‑80 (MQ=255) cAGCCGCAGCATCCAGCCGCAGGCAGCCCACCCCAACATCGTGCAGTTTCAGGCCGACCTGTGCCTGGCGGTCGCCCGCCTGgcac < 1:679189/86‑3 (MQ=37) cAGCCGCAGCATCCAGCCGCAGGCAGCCCACCCCAACATCGTGCAGTTTCAGGCCGACCTGTGCCTGGCGGTCGCCCGCCTGgcac > 2:679189/1‑84 (MQ=37) tCCAGCCGCAGGCAGCCCACCCCAACATCGTGCAGTTTCAGGCCGACCTGTGCCAGGCGGGCGACCGCCAGGCGCTGATCGATCACCTGCACGATCGTTATGACGGACTACGCGCCATCATCCACAACGCCTCGCTATGGCTCGacgacg > 2:738296/1‑150 (MQ=255) aGCCGCAGGCAGCCCACCCCAACATCGTGCAGTTTCAGGCCGACCTGTGCCAGGCGGGCGACCGCCAGGCGCTGATCGATCACCTGCACGATCGTTATGACGGACTACGCGCCATCATCCACAACGCCTCGCTATGGCTCGACGACGGcc > 2:378705/1‑150 (MQ=255) gcagCCCACCCCAACATCGTGCAGTTTCAGGCCGACCTGTGCCAGGCGGGCGACCGCTAGGCGCTGATCGATCACCTGCACGATCGTTATGACGGACTACGCGCCATCATCCACAGCGCCTCGCTATGGCTCGACGACGGCCTCGACAAc > 2:264929/1‑150 (MQ=255) gcagCCCACCCCAACATCGTGCAGTTTCAGGCCGACCTGTGCCAGGCGGGCGACCGCCAGGCGCTGATCGATCACCTGCACGATCGTTATGACGGACTACGCGCCATCATCCACAACGCCTCGCTATGGCTCGACGACGGCCTCGACAAc > 1:176389/1‑150 (MQ=255) gcagCCCACCCCAACATCGTGCAGTTTCAGGCCGACCTGTGCCAGGCGGGCGACCGCCAGGCGCTGATCGATCACCTGCACGATCGTTATGACGGACTACGCGCCATCATCCACAACGCCTCGCTATGGCTCGACGACGGCCTCGACAAc > 2:263917/1‑150 (MQ=255) gcagCCCACCCCAACATCGTGCAGTTTCAGGCCGACCTGTGCCAGGCGGGCGACCGCCAGGCGCTGATCGATCACCTGCACGATCGTTATGACGGACTACGCGCCATCATCCACAACGCCTCGCTATGGCTCGACGACGGCCTCGACAAc < 2:52151/150‑1 (MQ=255) ccACCCCAACATCGTGCAGTTTCAGGCCGACCTGTGCCAGGCCGGCGACCGCCAGGCGCTGATCGATCACCTGCACGATCGTTATGACGGACTACGCGCCATCATCCACAACGCCTCGCTATGGCTCGACGACGGCCTCGACAACCTCga < 2:254938/150‑1 (MQ=255) ccAACATCGTGCAGTTTCAGGCCGACCTGTGCCAGGCGGGCGACCGCCAGGCGCTGATCGATCACCTGCACGATCGTTATGACGGACTACGCGCCATCATCCACAACGCCTCGCTATGGCTCGACGACGGCCTCGACAACCTCGAGACCa < 2:104078/150‑1 (MQ=255) aGTTTCAGGCCGACCTGTGCCAGGCGGGCGACCGCCAGGCGCTGATCGATCACCTGCACGATCGTTATGACGGACTACGCGCCATCATCCACAACGCCTCGCTATGGCTCGACGACGGCCTCGACAACCTCGAGACCATGTTCCGCCTGc > 1:656192/1‑150 (MQ=255) tttCAGGCCGACCTGTGCCAGGCGGGCGACCGCCAGGCGCTGATCGATCACCTGCACGATCGTTATGACGGACTACGCGCCATCATCCACAACGCCTCGCTATGGCTCGACGACGGCCTCGACAACCTc < 2:16247/129‑1 (MQ=255) tttCAGGCCGACCTGTGCCAGGCGGGCGACCGCCAGGCGCTGATCGATCACCTGCACGATCGTTATGACGGACTACGCGCCATCATCCACAACGCCTCGCTATGGCTCGACGACGGCCTCGACAACCTc > 1:16247/1‑129 (MQ=255) aGGCGGGCGACCGCCAGGCGCTGATCGATCACCTGCACGATCGTTATGACGGACTACGCGCCATCATCCACAACGCCTCGCTATGGCTCGACGACGGCCTCGACAACCTCGAGACCATGTTCCGCCTGCACGTCGAAGCGCCCTaccacc > 1:759608/1‑150 (MQ=255) ggcgACCGCCAGGCGCTGATCGATCACCTGCACGATCGTTATGACGGACTACGCGCCATCATCCACAACGCCTCGCTATGGCTCGACGACGGCCTCGACAACCTCGAGACCATGTTCCGCCTGCACGTCGAAGCGCCCTACCACCTCAAc < 2:656192/150‑1 (MQ=255) ggcgACCGCCAGGCGCTGATCGATCACCTGCACGATCGTTATGACGGACTACGCGCCATCATCCACAACGCCTCGCTATGGCTCGACGACGGCCTCGACAACCTCGAGACCATGTTCCGCCTGCACGTCGAAGCGCCCTACCACCTCAAc > 1:224186/1‑150 (MQ=255) gcgACCGCTAGGCGCTGATCGATCACCTGCACGATCGTTATGACGGACTACGCGCCATCATCCACAACGCCTCGCTATGGCTCGACGACGGCCTCGACAACCTCGAGACCATGTCCCGCCTGCACGTCGAAGCGCCCTACCACCTCAAcc < 1:264929/150‑1 (MQ=255) | AGCCAGCCAGCGCGTCGGGCTGGCCCTGGCCCTCGAACTGGCACAGGCCGGCCATACGGTGGTCAGTGCCAGCCGCAGCATCCAGCCGCAGGCAGCCCACCCCAACATCGTGCAGTTTCAGGCCGACCTGTGCCAGGCGGGCGACCGCCAGGCGCTGATCGATCACCTGCACGATCGTTATGACGGACTACGCGCCATCATCCACAACGCCTCGCTATGGCTCGACGACGGCCTCGACAACCTCGAGACCATGTTCCGCCTGCACGTCGAAGCGCCCTACCACCTCAACC > NC_002947/5255590‑5255879 |
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: ‑ |