Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 5,131,909 | T→G | 11.5% | L72R (CTG→CGG) | PP_4517 → | hypothetical protein |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 5,131,909 | 0 | T | G | 11.5% | 101.9 / 6.5 | 43 | L72R (CTG→CGG) | PP_4517 | hypothetical protein |
Reads supporting (aligned to +/- strand): ref base T (18/20); new base G (2/3); total (20/23) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.15e-01 |
GCGCTGGACTATGTCGCCGCACACGGCTACCTCACCGCTCTTTCGATCAACGCCGAGGACGTGCCCGAGCGCGAATGGATCGACGCCCTGTTCGCCGAGGAGCCGCACTACGCCAGCGAAGCCCAGCGCACCGAGATCGAG‑GCCA‑CCCTGGTGGCGCTCAAGGCACACATTGCCCGCCAGC‑‑TGGCCAGCGACGAAGAATTCGACCTGCCCTGCGACCTCGACCTGACCGATGAGCCGGACGATTCCGACCTGCGCGGCTGGTGCATCGGCTTCATGGAAGGCGTGTTCCTGCGGGA > NC_002947/5131761‑5132056 | gcgcTGGACTATGTCGCCGCACACGGCTACCTCACCGCTCTTTCGATCAACGCCGAGGACGTGCCCGAGCGCGAATGGATCGACGCCCTGTTCGCCGAGGAGCCGCACTACGCCAGCGAAGCCCAGCGCACCGAGATCGAG‑GCCA‑CCCtg > 2:698061/1‑150 (MQ=255) tGGACTATGTCGCCGCACACGGCTACCTCACCGCTCTTTCGATCAACGCCGAGGACGTGCCCGAGCGCGAATGGATCGACGCCCTGTTCGCCGAGGAGCCGCACTACGCCAGCGAAGCCCAGCGCACCGAGATCGAG‑GCCA‑CCCtggtgg > 1:845872/1‑150 (MQ=255) aCTATGTCGCCGCACACGGCTACCTCACCGCTCTTTCGATCAACGCCGAGGACGTGCCCGAGCGCGAATGGATCGACGCCCTGTTCGCCGAGGAGCCGCACTACGCCAGCGAAGCCCAGCGCACCGAGATCGAG‑GCCA‑CCAGGGTGgcct < 2:82877/150‑3 (MQ=255) cTATGTCGCCGCACACGGCTACCTCACCGCTCTTTCGATCAACGCCGAGGACGTGCCCGAGCGCGAATGGATCGACGCCCTGTTCGCCGAGGAGCCGCACTACGCCAGCGAAGCCCAGCGCACCGAGATCGAG‑GCCA‑CCAGGGTGgcctc < 2:408019/150‑4 (MQ=255) tATGTCGCCGCACACGGCTACCTCACCGCTCTTTCGATCAACGCCGAGGACGTGCCCGAGCGCGAATGGATCGACGCCCTGTTCGCCGAGGAGCCGCACTACGCCAGCGAAGCCCAGCGCACCGAGATCGAG‑GCCA‑CCCTGGTGGCGCTc < 1:263075/150‑1 (MQ=255) gTCGCCGCACACGGCTACCTCACCGCTCTTTCGATCAACGCCGAGGACGTGCCCGAGCGCGAATGGATCGACGCCCTGTTCGCCGAGGAGCCGCACTACGCCAGCGAAGCCCAGCGCACCGAGATCGAG‑GCCA‑CCCTGGTGGCGCTCAAg > 1:528863/1‑150 (MQ=255) gTCGCCGCACACGGCTACCTCACCGCTCTTGCGATCAACGCCGAGGACGTGCCCGAGCGCGAATGGACCGCCGCCCGGTTCGCCGAGGATGGGCACTACGCCAGCGAAGCCCAGCGCACCGAGATCGAG‑GCCA‑CCCTGGTGGCGCTCAAg < 1:104870/150‑1 (MQ=255) gccgcACACGGCTACCTCACCGCTCTTTCGATCAACGCCGAGGACGTGCCCGAGCGCGAATGGATCGACGCCCTGTTCGCCGAGGAGCCGCACTACGCCAGCGAAGCCCAGCGCACCGAGATCGAG‑GCCA‑CCCTGGTGGCGCTCAAGGca < 1:129864/150‑1 (MQ=255) ccgcACACGGCTACCTCACCGCTCTTTCGATCAACGCCGAGGACGTGCCCGAGCGCGAATGGATCGACGCCCTGTTCGCCGAGGAGCCGCACTACGCCAGCGAAGCCCAGCGCACCGAGATCGAG‑GCCA‑CCCTGGTGGCGCTCAAGGcac > 1:702464/1‑150 (MQ=255) acacGGCTACCTCACCGCTCTTTCGATCAACGCCGAGGACGTGCCCGAGCGCGAATGGATCGACGCCCTGTTCGCCGAGGAGCCGCACTACGCCAGCGAAGCCCAGCGCACCGAGATCGAG‑GCCA‑CCCTGGTGGCGCTCAAGGCACACAt > 1:747763/1‑150 (MQ=255) acGGCTACCTCACCGCTCTTTCGATCAACGCCTATTACGTGCCCGAGCGCGAATGGATCGACGCCCTGTTCGCCGAGGAGCCGCACTACGCCAGCGAAGCCCAGCGCACCGAGATCGAG‑GCCA‑CCCTGGTGGCGCTCAAGGCACACATTg < 2:747763/150‑1 (MQ=255) ggCTACCTCACCGCTCTTTCGATCAACGCCGAGGACGTGCCCGAGCGCGAATGGATCGACGCCCTGTTCGCCGAGGAGCCGCACTACGCCAGCGAAGCCCAGCGCACCGAGATCGAG‑GCCA‑CCCTGGTGGCGCTCAAGGCACACATTGcc > 2:627796/1‑150 (MQ=255) cTCACCGCTCTGTCGATCAACGCCGAGGACGTGCCCGAGCGCGAATGGATCGACTCCCTGTTCGCCGAGGAGCCGCCCTACGCCAGCGAAGCCCAGCGCACCGAGACCCAG‑GCCA‑CCCGGGCGGCGCTCAAGGCACACATTGCCCGcccc > 2:393539/1‑148 (MQ=255) gATCAACGCCGAGGACGTGCCCGAGCGCGAATGGATCGACGCCCTGTTCGCCGAGGAGCCGCACTACGCCAGCGAAGCCCAGCGCACCGAGATCGAG‑GCCA‑CCCTGGTGTCCCTCAAGGCACACATTGCCCGCCAGC‑‑TGGCCAGcgacga > 1:649215/1‑150 (MQ=255) cGAGGACGTGCCCGAGCGCGAATGGATCGACGCCCTGTTCGCCGAGGAGCCGCACTACGCCAGCGAAGCCCAGCGCACCGAGATCGAG‑GCCA‑CCCTGGTGGCGCTCAAGGCACACATTGCCCGCCAGC‑‑TGGCCAGCGACGAAGAATTCGa < 1:627796/150‑1 (MQ=255) gAGGACGTGCCCGAGCGCGAATGGATCGACGCCCTGTTCGCCGAGGAGCCGCACTACGCCAGCGAAGCCCAGCGCACCGAGATCGAG‑GCCA‑CCCTGGTGGCGCTCAAGGCACACATTGCCCGCCAGC‑‑TGGCCAGCGACGAAGAATTCGAc < 1:382929/150‑1 (MQ=255) gAATGGATCGACGCCCTGTTCGCCGAGGAGCCGCACTACGCCAGCGAAGCCCAGCGCACCGAGATCGAGCGCCA‑CCAGGGTGgcctcg < 1:247455/88‑5 (MQ=255) gAATGGATCGACGCCCTGTTCGCCGAGGAGCCGCACTACGCCAGCGAAGCCCAGCGCACCGAGATCGAGCGCCA‑CCAGGGTGgcctcg > 2:247455/1‑84 (MQ=255) aTGGATCGACGCCCTGTTCGCCGAGGAGCCGCACTACGCCAGCGAAGCCCAGCGCACCGAGATCGAG‑GCCA‑CCCTGGTGGCGCTCAAGGCACACATTGCCCGCCAGC‑‑TGGCCAGCGAGGACGAATACCACCTGCCCTGcgagctccacct > 2:130497/1‑150 (MQ=255) aTGGATCGACGCCCTGTTCGCCGAGGAGCCGCACTACGCCAGCGAAGCCCAGCGCACCGAGATCGAG‑GCCA‑CCCTGGTGGCGCTCAAGGCACACATTGCCCGCCAGC‑‑TGGCCAGCGACGAAGAATTCGACCTGCCCTGcgacctcgacct < 2:845872/150‑1 (MQ=255) gATCGACGCCCTGTTCGCCGAGGAGCCGCACTACGCCAGCGAAGCCCAGCGCACCGAGATCGAG‑GCCA‑CCCTGCTGGCGCACAAGGCACACGTTGCCCGCCAGC‑‑TGGCCAGCGACGAAGAAGTCGACCTGCCCGTCGTCCTCGTGCTGAc > 2:393849/1‑150 (MQ=255) cGCACTGTTCGCCGAGGAGCCGCACTACGCCAGCGAAGCCCAGCGCACCGAGATCGAG‑GCCA‑CCCTGGTGGCGCTCAAGGCACACATTGCCCGCCAGC‑‑TGGCCAGCGACGAAGAATTCGACCTGCCCTGCGACCTCGACCTGACCGATGa < 1:393539/150‑1 (MQ=255) gTTCTCCGAGGTGCCGCACTACGCCAGCGAAGCCCAGCGCACCGAGATCGAG‑GCCA‑CCCTGGTGGCGCTCAAGGCACACATTGCCCGCCAGC‑‑TGGCCAGCTACGAAGAATTCGACCTGCCCTGCGACCTCGACCTGACCGATGAGCCGGa < 1:130497/150‑1 (MQ=255) gCCGAGGAGCCGCACTACGCCAGCGAAGCCCAGCGCACCGAGATCGAG‑GCCA‑CCCTGGTGGCGCTCAAGGCACACATTGCCCGCCAGC‑‑TGGCCAGCGACGAAGAATTCGACCTGCCCTGCGACCTCGACCTGACCGATGAGCCGGa > 2:174207/1‑146 (MQ=255) gCCGAGGAGCCGCACTACGCCAGCGAAGCCCAGCGCACCGAGATCGAG‑GCCA‑CCCTGGTGGCGCTCAAGGCACACATTGCCCGCCAGC‑‑TGGCCAGCGACGAAGAATTCGACCTGCCCTGCGACCTCGACCTGACCGATGAGCCGGa < 1:174207/146‑1 (MQ=255) cGAGGAGCCGCACTACGCCAGCGAAGCCCAGCGCACCGAGATCGAG‑GCCA‑CCCTGGTGGCGCTCAAGGCACACATTGCCCGCCAGC‑‑TGGCCAGCGACGAAGAATTCGACCTGCCCTGCGACCTCGACCTGACCGATGAGCCGGACGATTc > 1:127422/1‑150 (MQ=255) aggagCCGCACTACGCCAGCGAAGCCCAGCGCACCGAGATCGAG‑GCCA‑CCCTGGTGGCGCTCAAGGCACACATTGCCCGCCAGC‑‑TGGCCAGCGACGAAGAATTCGACCTGCCCTGCGACCTCGACCTGACCGATGAGCCGGACGATTCCg < 1:393849/150‑1 (MQ=255) gagCCGCACTACGCCAGCGAAGCCCAGCGCACCGAGATCGAG‑GCCA‑CCCTGGTGGCGCTCAAGGCACACATTGCCCGCCAGC‑‑TGGCCAGCGACGAAGAATTCGACCTGCCCTGCGACCTCGACCTGACCGATGAGCCGGACGATTCCGAc < 2:528863/150‑1 (MQ=255) ccAGCGAAGCCCAGCGCACCGAGATCGAG‑GCCA‑CCCTGGTGGCGCTCAAGGCACACATTGCCCGCCAGC‑‑TGGCCAGCGACGAAGAATTCGACCTGCCCTGCGACCTCGACCTGACCGATGAGCCGGACGATTCCGACCTGCGCGGCTGGt < 2:200451/150‑1 (MQ=255) gCGAAGCCCAGCGCACCGAGATCGAG‑GCCA‑CCCTGGTGGCGCTCAAGGCACACATTGCCCGCCAGC‑‑TGGCCAGCGACGAAGAATTCGACCTGCCCTGCGACCTCGACCTGACCGATGAGCCGGACGAtt < 1:133371/129‑1 (MQ=255) gCGAAGCCCAGCGCACCGAGATCGAG‑GCCA‑CCCTGGTGGCGCTCAAGGCACACATTGCCCGCCAGC‑‑TGGCCAGCGACGAAGAATTCGACCTGCCCTGCGACCTCGACCTGACCGATGAGCCGGACGAtt > 2:133371/1‑129 (MQ=255) gAAGCCCAGCGCACCGAGATCGAG‑GCCA‑CCCTGGTGGCGCTCAAGGCACACATTGCCCGCCAGC‑‑TGGCCAGCGACgaa < 1:630715/78‑1 (MQ=255) gAAGCCCAGCGCACCGAGATCGAG‑GCCA‑CCCTGGTGGCGCTCAAGGCACACATTGCCCGCCAGC‑‑TGGCCAGCGACgaa > 2:630715/1‑78 (MQ=255) gAAGCCCAGCGCACCGAGATCGAG‑GCCA‑CCCTGGTGGCGCTCAAGGCACACATTGCCCGCCAGC‑‑TGGCCAGCGACGAAGAATTCGACCTGCCCTGCGACCTCGACCTGACCGATGAGCCGGACGATTCCGACCTGCGCGGCTGGTGCATc < 2:127422/150‑1 (MQ=255) aaaGCCCAGCGCACCGCCATCG‑G‑GCCATGCGTGGTGGTGGGGAAGGCACACATTGCCCGCCAGCATTCGCCATC‑‑TTTAGAATTCGACCTGCCCTGCGACCTCGACCTGACCGATGAGCCGGACGATTCCGACCTGCGCGGCTGGTGCATc < 1:787255/149‑1 (MQ=255) cAGCGCACCGAGATCGAG‑GCCA‑CCCTGGTGGCGCTCAAGGCACACATTGCCCGCCAGC‑‑TGGCCAGCGACGAAGAATTCGACCTGCCCTGCGACCTCGACCTGACCGATGAGCCGGACGATTCCGACCTGCGCGGCTGGTGCATCGGCTTc > 2:761073/1‑150 (MQ=255) cGAGATCGAG‑GCCA‑CCCTGGTGGCGCTCAAGGCACACATTGCCCGCCAGC‑‑TGGCCAGCGACGAAGAATTCGACCTGCCCTGCGACCTCGACCTGACCGATGAGCCGGACGATTCCGACCTGCGCGGCTGGTGCATc < 2:250251/136‑1 (MQ=255) cGAGATCGAG‑GCCA‑CCCTGGTGGCGCTCAAGGCACACATTGCCCGCCAGC‑‑TGGCCAGCGACGAAGAATTCGACCTGCCCTGCGACCTCGACCTGACCGATGAGCCGGACGATTCCGACCTGCGCGGCTGGTGCATc > 1:250251/1‑136 (MQ=255) gAG‑GCCA‑CCCTGGTGGCGCTCAAGGCACACATTGCCCGCCAGC‑‑TGGCCAGCGACGAAGAATTCGACCTGCCCTGCGACCTCGACCTGACCGATGCGACCGTCTATTCCGACCTGCGCGGCTGGTGCATCGGCTTCATGGAAGGCGTGTTc > 1:692332/1‑150 (MQ=255) a‑CCCTGGTGGCGCTCAAGGCACACATTGCCCGCCAGC‑‑TGGCCAGCGACGAAGAATTCGACCTGCCCTGCGACCTCGACCTGACCGATGAGCCGGACGATTCCGa > 1:806919/1‑104 (MQ=255) a‑CCCTGGTGGCGCTCAAGGCACACATTGCCCGCCAGC‑‑TGGCCAGCGACGAAGAATTCGACCTGCCCTGCGACCTCGACCTGACCGATGAGCCGGACGATTCCGa < 2:806919/104‑1 (MQ=255) cccTGGTGGCGCTCAAGGCACACATTGCCCGCCAGC‑‑TGGCCAGCGACGAAGAATTCGACCTGCCCTGCGACCTCGACCTGACCGATGAGCCGGACGATTCCGACCTGCGCGGCTGGTGCATCGGCTTCATGGAAGGCGTGTTCCTGCggg > 1:288065/1‑150 (MQ=255) ccTGGTGGCGCTCAAGGCACACATTGCCCGCCAGC‑‑TGGCCAGCGACGACGAATTCGACCTGCCCTGCGACCTCGACCTGACCGATGAGCCGGACGATTCCGACCTGCGCGGCTGGTGCATCGGCTTCATGGAAGGCGTGTTCCTGCGGGa > 1:548486/1‑150 (MQ=255) ccTGGTGGCGCTCAAGGCACACATTGCCCGCCAGC‑‑TGGCCAGCGACGAAGAATTCGACCTGCCCTGCGACCTCGACCTGACCGATGAGCCGGACGATTCCg > 1:75253/1‑101 (MQ=255) ccTGGTGGCGCTCAAGGCACACATTGCCCGCCAGC‑‑TGGCCAGCGACGAAGAATTCGACCTGCCCTGCGACCTCGACCTGACCGATGAGCCGGACGATTCCg < 2:75253/101‑1 (MQ=255) | GCGCTGGACTATGTCGCCGCACACGGCTACCTCACCGCTCTTTCGATCAACGCCGAGGACGTGCCCGAGCGCGAATGGATCGACGCCCTGTTCGCCGAGGAGCCGCACTACGCCAGCGAAGCCCAGCGCACCGAGATCGAG‑GCCA‑CCCTGGTGGCGCTCAAGGCACACATTGCCCGCCAGC‑‑TGGCCAGCGACGAAGAATTCGACCTGCCCTGCGACCTCGACCTGACCGATGAGCCGGACGATTCCGACCTGCGCGGCTGGTGCATCGGCTTCATGGAAGGCGTGTTCCTGCGGGA > NC_002947/5131761‑5132056 |
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: ‑ |