Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 4,097,638 | C→G | 10.0% | G28R (GGG→CGG) | PP_3605 ← | LysR family transcriptional regulator |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 4,097,638 | 0 | C | G | 10.0% | 115.0 / 4.9 | 40 | G28R (GGG→CGG) | PP_3605 | LysR family transcriptional regulator |
Reads supporting (aligned to +/- strand): ref base C (18/18); new base G (2/2); total (20/20) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.70e-01 |
TCGAGAACCTTGGCTACCAGCCGCTCGCCCTCATCGGTAAAGCGCAGCTGCCGCGTACTGCGCTCCACCAGGCGAATGCCCAGCCGCGTCTCTAAAGCACTCAGGCGTTTCGACAGCACCGTAGGGTGGCGCTCCAGTTGCCGCCCGGCGGCAGCAAACGAGCCCTGCTCATACAAGGCCAGCAGGGTGGCAAGTTCATCAGCGCGACGACTGTCGAAGGTATCCATGGCAAGGCTCGGCTCAGGGCAGTGAAATCAGGGGGCGGATCAGGCAGTCAGCACCAGCGCGCCCT > NC_002947/4097493‑4097784 | tCGAGAACCTTGGCTACCAGCCGCTCGCCCTCATCGGTAAAGCGCAGCTGCCGCGTACTGCGCTCCACCAGGCGAATGCCCAGCCGCGTCTCTAAAGGACTCGGGCGTTTCGACAGCACCGTAGGGTGGCGCTCCAGTTGCCGCCcggcg > 2:1051300/1‑150 (MQ=255) agaACCTTGGCTACCAGCCGCTCGCCCTCATCGGTAAAGCGCAGCTGCCGCGTACTGCGCTCCACCAGGCGAATGCCCAGCCGCGTCTCTAAAGCACTCAGGCGTTTCGACAGCACCGTAGGGTGGCGCTCCAGTTGCCGCCCGGCGgca < 1:898914/150‑1 (MQ=255) gaACCTTGGCTACCAGCCGCTCGCCCTCATCGGTAAAGCGCAGCTGCCGCGTACTGCGCTCCACCAGGCGAATGCCCAGCCGCGTCTCTAAAGCACTCAGGCGTTTCGACAGCACCGTAGGGTGGCGCTCCAGTTGCCGCCGGGCGgcaa < 2:1040815/150‑2 (MQ=255) tACCAGCCGCTCGCCCTCATCGGTAAAGCGCAGCTGCCGCGTACTGCGCTCCACCAGGCGAATGCCCAGCCGCGTCTCTAAAGCACTCAGGCGTTTCGACAGCACCGTAGGGTGGCGCTCCAGTTGCCGCCCGGCGGCAGCAAACGAGcc > 2:177411/1‑150 (MQ=255) gCCCTCATCGGTAAAGCGCAGCTGCCGCGTACTGCGCTCCACCAGGCGAATGCCCAGCCGCGTCTCTAAAGCACTCAGGCGTTTCGACAGCACCGTAGGGTGGCGCTCCAGTTGCCGCCCGGCGGCAGCAAACGAGCCCTGCTCATACaa < 1:1052430/150‑1 (MQ=255) ccTCATCGGTAAAGCGCAGCTGCCGCGTACTGCGCTCCACCAGGCGAATGCCCAGCCGCGTCTCTAAAGCACTCAGGCGTTTCGACAGCACCGTAGGGTGGCGCTCCAGTTGCCGCCCGGCGGCAGCAAACGAGCCCTGCTCATACAAgg > 1:403297/1‑150 (MQ=255) gTAAAGCGCAGCTGCCGCGTACTGCGCTCCACCAGGCGAATGCCCAGCCGCGTCTCTAAAGCACTCAGGCGTTTCGACAGCACCGTAGGGTGGCGCTCCAGTTGCCGCCCGGCGGCAGCAAACGAGCCCTGCTCATACAAGGCCAGCAgg > 2:1082924/1‑150 (MQ=255) aaaGCGCAGCTGCCGCGTACTGCGCTCCACCAGGCGAATGCCCAGCCGCGTCTCTAAAGCACTCAGGCGTTTCGACAGCACCGTAGGGTGGCGCTCCAGTTGCCGCCCGGCGGCAGCAAACGAGCCCTGCTCATACAAGGCCAGCAGGGt < 2:454217/150‑1 (MQ=255) aaGCGCAGCTGCCGCGTACTGCGCTCCACCAGGCGAATGCCCAGCCGCGTCTCTAAAGCACTCAGGCGTTTCGACAGCACCGTAGGGTGGCGCTCCAGTTGCCGCCCGGCGGCAGCAAACTAGCCCTGCTCATACAAGGCCAGCAGGGTg < 2:89561/150‑1 (MQ=255) aaGCGCAGCTGCCGCGTACTGCGCTCCACCAGGCGAATGCCCAGCCGCGTCTCTAAAGCACTCAGGCGTTTCGACAGCACCGTAGGGTGGCGCTCCAGTTGCCGCCCGGCGGCAGCAAACGAGCCCTGCTCATACAAGGCCAGCAGGGTg < 2:333061/150‑1 (MQ=255) tGCGCTCCACCAGGCGAATGCCCAGCCGCGTCTCTAAAGCACTCAGGCGTTTCAAAAGCCTCGTAGGGTGGCGCGCCATGTGCCGCCCGGCGGCAGCAAACGAGCCCTGCTCATACAAGGCCAGCAGGGTGGCAAGTTCATCAGCGcgac > 1:136128/1‑150 (MQ=255) cTCCACCAGGCGAATGCCCAGCCGCGTCTCTAAAGCACTCAGGCGTTTCGACAGCACCGTAGGGTGGCGCTCCAGTTGCCGCCCGGCGGCAGCAAACGAGCCCTGCTCATACAAGGCCAGCAGGGTGGCAAGTTCATCAGCGCGACGACt < 1:696438/150‑1 (MQ=255) cTCCACCAGGCGAATGCCCAGCCGCGTCTCTAAAGCACTCAGGCGTTTCGACAGCACCGTAGGGTGGCGCTCCAGTTGCCGCCCGGCGGCAGCAAACGAGCCCTGCTCATACAAGGCCAGCAGGGTGGCAAGTTCATCAGCGCGACGACt > 2:696438/1‑150 (MQ=255) tCCACCAGGCGAATGCCCAGCCGCGTCTCTAAAGCACTCAGGCGTTTCGACAGCACCGTAGGGTGGCGCTCCAGTTGCCGCCCGGCGGCAGCAAACGAGCCCTGCTCATACAAGGCCAGCAGGGTGGCAAGTTCATCAGCGCGACGACTg < 1:729342/150‑1 (MQ=255) accaGGCGAATGCCCAGCCGCGTCTCTAAAGCACTCAGGCGTTTCGACAGCACCGTAGGGTGGCGCTCCAGTTGCCGCCGGGCGgcaactg < 1:718654/91‑5 (MQ=255) accaGGCGAATGCCCAGCCGCGTCTCTAAAGCACTCAGGCGTTTCGACAGCACCGTAGGGTGGCGCTCCAGTTGCCGCCGGGCGgcaactg > 2:718654/1‑87 (MQ=255) accaGGCGAATGCCCAGCCGCGTCTCTAAAGCACTCAGGCGTTTCGACAGCACCGTAGGGTGGCGCTCCAGTTGCCGCCCGGCGGCAGCAAACGAGCCCTGCTCATACAAGGCCAGCAGGGTGGCAAGTTCATCAGCGCGACGACTGTCg > 1:251533/1‑150 (MQ=255) ccaGGCGAATGCCCAGCCGCGTCTCTAAAGCACTCCAAAATTTCGACAGCACCGTAGGGTGGCGCTCCAGTTGCCGCCCGGCGGCAGCAAACGAGCCCTGCTCATACAAGGCCAGCAGGGTGGCAAGTTCATCAGCGCGACGACTGTCGa < 1:1081875/150‑1 (MQ=255) aGGCGAATGCCCAGCCGCGTCTCTAAAGCACTCAGGCGTTTCGACAGCACCGTAGGGTGGCGCTCCAGTTGCCGCCCGGCGGCAGCAAACGAGCCCTGCTCATACAAGGCCAGCAGGGTGGCAAGTTCATCAGCGCGACGACTGTCGAAg < 1:177411/150‑1 (MQ=255) aGGCGAATGCCCAGCCGCGTCTCTAAAGCACTCAGGCGTTTCGACAGCACCGTAGGGTGGCGCTCCAGTTGCCGCCCGGCGGCAGCAAACGAGCCCTGCTCATACAAGGCCAGCAGGGTGGCAAGTTCATCAGCGCGACGACTGTCGAAg < 2:309998/150‑1 (MQ=255) ggCGAATGCCCAGCCGCGTCTCTAAAGCACTCAGGCGTTTCGACAGCACCGTAGGGTGGCGCTCCAGTTGCCGCCCGGCGGCAGCAAACGAGCCCTGCTCATACAAGGCCAGCAGGGTGGCAAGTTCATCAGCGCGACGACTGTCGAAgg > 1:146282/1‑150 (MQ=255) gCGAATGCCCAGCCGCGTCTCTAAAGCACTCAGGCGTTTCGACAGCACCGTAGGGTGGCGCTCCAGTTGCCGCCCGGCGGCAGCAAACGAGCCCTGCTCATACAAGGCCAGCAGGGTGGCAAGTTCATCAGCGCGACGACTGTCGAAGGt > 1:819299/1‑150 (MQ=255) aTGCCCAGCCGCGTCTCTAAAGCACTCAGGCGTTTCGACAGCACCGTAGGGTGGCGCTCCAGTTGCCGCCCGGCGGCAGCAAACGAGCCCTGCTCATACAAGGCCAGCAGGGTGGCAAGTTCATCAGCGCGACGACTGTCGAAGGTATcc < 2:403297/150‑1 (MQ=255) ccAGCCGCGTCTCTAAAGCACGCAGGCGTTTCGACAGCACCGTAGGGTGGCGCTCCAGTTGCCGCCCGGCGGCAGCAAACGAGCCCTGCTCATACAAGGCTAGCAGGGTGGCAAGTTCATCAGCGCGACGACTGTCGAAGGTATCCATgg > 2:572240/1‑150 (MQ=255) cgcgTCTCTAAAGCACTCAGGCGTTTCGACAGCACCGTAGGGTGGCGCTCCAGTTGCCGCCCGGCGGCAGCAAACGAGCCCTGCTCATACAAGGCCAGCAGGGTGGCAAGTTCATCAGCGCGACGACTGTCGAAGGTATCCATGGCAAgg < 2:146282/150‑1 (MQ=255) gcgTCTCTTAAGCACTCAGGCGTTTCGACAGGACCGTAGGGTGGNACAANAGTTGCAGCCCGGCGGCAGCAAGCAAGCCCTGCTCATACAAGGCAAGAAGGTTGGCAAGTTCATCAGCGCGACGACTGTCGCAGGTATCCGTgggcggga > 1:248516/1‑144 (MQ=255) gcgTCTCTAAAGCACTCAGGCGTTTCGACAGCACCGTAGGGTGGCGCTCCAGTTGCCGCCCGGCGGCAGCAAACGAGCCCTGCTCATACAAGGCCAGCAGGGTGGCAAGTTCATCAGCGCGACGACTGTCGAAGGTATAACTGGCAAggc < 2:103500/150‑1 (MQ=255) aGCACTCAGGCGTTTCGACAGCACCGTAGGGTGGCGCTCCAGTTGCCGCCCGGCGGCAGCAAACGAGCCCTGCTCATACAAGGCCAACAGGGTGGCAAGTTCATCAGCGCGACGACTGTCGAAGGTATCCATGGCAAGGCTCGGCTCAgg > 2:112448/1‑150 (MQ=255) tCAGGCGTTTCGACAGCACCGTAGGGTGGCGCTCCAGTTGCCGCCCGGCGGCAGCAAACGAGCCCTGCTCATACAAGGCCAGCAGGGTGGCAAGTTCATCAGCGCGACGACTGTCGAAGGTATCCATGGCAAGGCTCGGCTCAGGGCAGt > 2:815331/1‑150 (MQ=255) gTTTCGACAGCACCGTAGGGTGGCGCTCCAGTTGCCGCCCGGCGGCAGCAAACGAGCCCTGCTCATACAAGGCCAGCAGGGTGGCAAGTTCATCAGCGCGACGACTGTCGAAGGTATCCATGGCAAGGCTCGGCTCAGGGCAGTGAAATc > 1:948780/1‑150 (MQ=255) cGTAGGGTGGCGCTCCAGTTGCCGCCCGGCGGCAGCAAACGAGCCCTGCTCATACAAGGCCAGCAGGGTGGCAAGTTCATCAGCGCGACGACTGTCGAAGGTATCCATGGCAAGGCTCGGCTCAGGGCAGTCACATCAGCAAGCGGAtga < 2:15880/150‑3 (MQ=255) gTAGGGTGGCGCTCCAGTTGCCGCCCGGCGGCAGCAAACGAGCCCTGCTCATACAAGGCCAGCAGGGTGGCAAGTTCATCAGCGCGACGACTGTCGAAGGTATCCATGGCAAGGCTCGGCTCAGGGCAGTGAAATCAGGGAGCGGATCAg > 1:597001/1‑150 (MQ=255) tAGGGTGGCGCTCCAGTTGCCGCCCGGCGGCAGCAAACGAGCCCTGCTCATACAAGGCCAGCAGGGTGGCAAGTTCATCAGCGCGACGACTGTCGAAGGTATCCATGGCAAGGCTCGGCTCAGGGCAGTGAAATCAGGGGTCGGATCAgg < 1:1082924/150‑1 (MQ=255) ggTGGCGCTCCAGTTGCCGCCCGGCGGCAGCAAACGAGCCCTGCTCATACAAGGCCAGCAGGGTGGCAAGTTCATCAGCGCGACGACTGTCGAAGGTATCCATGGCAAGGCTCGGCTCAGGGCAGTGAAATCAGGGGGCGGATCAGGCa < 1:881991/149‑1 (MQ=255) ggTGGCGCTCCAGTTGCCGCCCGGCGGCAGCAAACGAGCCCTGCTCATACAAGGCCAGCAGGGTGGCAAGTTCATCAGCGCGACGACTGTCGAAGGTATCCATGGCAAGGCTCGGCTCAGGGCAGTGAAATCAGGGGGCGGATCAGGCa > 2:881991/1‑149 (MQ=255) gTGGCGCTCCAGTTGCCGCCCGGCGGCAGCAAACGAGCCCTGCTCATACAAGGCCAGCAGGGTGGCAAGTTCATCAGCGCGACGACTGTCGAAGGTATCCATGGCAAGGCTCGGCTCAGGGCAGTGAAATCAGGGGGCGGATCAGGCAGt < 2:473995/150‑1 (MQ=255) tCCAGTTGCCGCCCGGCGGCAGCAAACGAGCCCTGCTCATACAAGGCCAGCAGGGTGGCAAGTTCATCAGCGCGACGACTGTCGAAGGTATCCATGGCAAGGCTCGGCTCAGGGCAGTGAAATCAGGGGGCGGATCAGGCAGTCAGCAcc > 2:601113/1‑150 (MQ=255) gtttgctGCCGCCGGGCGGCAACAAACGAGCCCTGCTCATACAAGGCCAGCAGGGTGGCAAGTTCATCAGCGCGACGACTGTCGAAGGTATCCATGGCAAGGCTCGGCTCAGGGCAGTGAAATCAGGGGGCGGATCAGGCAGTCAGCAcc > 1:670613/7‑150 (MQ=255) gTTGCCGCCCGGCGGCAGCAAACGAGCCCTGCTCATACAAGGCCAGCAGGGTGGCAAGTTCATCAGCGCGACGACTGTCGAAGGTATCCATGGCAAGGCTCGGCTCAGGGCAGTGAAATCAGGGGGCGGATCAGGCAGTCAGCACCAgcg > 2:1051532/1‑150 (MQ=255) gccCGGCGGCAGCAAACGAGCCCTGCTCATACAAGGCCAGCAGGGTGGCAAGTTCATCAGCGCGACGACTGTCGAAGGTATCCATGGCAAGGCTCGGCTCAGGGCAGTGAAATCAGGGGGCGGATCAGGCAGTCAGCACCAGCGCGCCCt < 2:434312/150‑1 (MQ=255) gccCGGCGGCAGCAAACGAGCCCTGCTCATACAAGGCCAGCAGGGTGGCAAGTTCATCAGCGCGACGACTGTCGAAGGTATCCATGGCAAGGCTCGGCTCAGGGCAGTGAAATCAGGGGGCGGATCAGGCAGTCAGCACCAGCGCGCCCt > 2:708068/1‑150 (MQ=255) | TCGAGAACCTTGGCTACCAGCCGCTCGCCCTCATCGGTAAAGCGCAGCTGCCGCGTACTGCGCTCCACCAGGCGAATGCCCAGCCGCGTCTCTAAAGCACTCAGGCGTTTCGACAGCACCGTAGGGTGGCGCTCCAGTTGCCGCCCGGCGGCAGCAAACGAGCCCTGCTCATACAAGGCCAGCAGGGTGGCAAGTTCATCAGCGCGACGACTGTCGAAGGTATCCATGGCAAGGCTCGGCTCAGGGCAGTGAAATCAGGGGGCGGATCAGGCAGTCAGCACCAGCGCGCCCT > NC_002947/4097493‑4097784 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |