Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 4,287,891 | C→A | 21.1% | R153L (CGT→CTT) | PP_3759 ← | CheB‑like methylesterase |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 4,287,891 | 0 | C | A | 21.1% | 33.2 / 4.6 | 19 | R153L (CGT→CTT) | PP_3759 | CheB‑like methylesterase |
Reads supporting (aligned to +/- strand): ref base C (5/10); new base A (2/2); total (7/12) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 6.03e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.38e-01 |
TCAGCAGTTTGGCGGTGGTGTGGCTTAGCATGCGCTGTATTCCAGGGTCGCGAGCAACTGCCCGATACCGCTCAGAGTAAGGATATGGTCGGGCTGGTGCAGGGCCAGGGCGGCCTCCGGCATCAGGGCAACCTGTGCGTCACGAGGGTCCTGGACGATCGTGCGCCCCCCATTGTTCTTGATGTAGGCAAGGCCCCTGGCGCCGTCTTCGTTGGCCCCGGTGAGCAGCACACCCAGCAGGCCATCGCCATAGGCATCGGCCGCCGAGATGAACAGGAAATCGATC > NC_002947/4287749‑4288034 | tCAGCAGTTTGGCGGTGGTGTGGCTTAGCATGCGCTGTATTCCAGGGTCGCGAGCAACTGCCCGATACCGCTCAGAGTAAGGATATGGTCGGGCTGGTGCAGGGCCAGGGCGGCCTCCGGCATCAGGGCAACCTGTGCGTCACGAGGGTc > 1:44986/1‑150 (MQ=255) gtggtgTGGCTTAGCATGCGCTGTATTCCAGGGTCGCGAGCAACTGCCCGATACCGCTCAGAGTAAGGATATGGTCGGGCTGGTGCAGGGCCAGGGCGGCCTCCGGCATCAGGGCAACCTGTGCGTCACGAGGGTCCTGGACGATCGTgc < 2:144305/150‑1 (MQ=255) tGCGCTGTATTCCAGGGTCGCGAGCAACTGCCCGATACCGCTCAGAGTAAGGATATGGTCGGGCTGGTGCAGGGCCAGGGCGGCCTCCGGCATCAGGGCAACCTGTGCGTCACGAGGGTCCTGGACGATCGTGCGCCCCCCATTGTTCtt > 2:367141/1‑150 (MQ=255) cgcTGTATTCCAGGGTCGCGAGCAACTGCCCGATACCGCTCAGAGTAAGGATATGGTCGGGCTGGTGCAGGGCCAGGGCGGCCTCCGGCATCAGGGCAACCTGTGCGTCACGAGGGTCCTGGACGATCGTGCGCCCCCCATTGTTCTTGa < 1:170708/150‑1 (MQ=255) tGTATTCCAGGGTCGCGAGCAACTGCCCGATACCGCTCAGAGTAAGGATATGGTCGGGCTGGTGCAGGGCCAGGGCGGCCTCCGGCATCAGGGCAACCTGTGCGTCACGAGGGTCCTGGACGATCGTGCGCCCCCCATTGTTCTTGATGt < 2:244177/150‑1 (MQ=255) aGGGTCGCGAGCAACTGCCCGATACCGCTCAGAGTAAGGATATGGTCGGGCTGGTGCAGGGCCAGGGCGGCCTCCGGCATCAGGGCAACCTGTGCGTCACGAGGGTCCTGGACGATCGTGCGCCCCCCATTGTTCTTGATGTAGGCAAgg < 2:44986/150‑1 (MQ=255) gTCGCGAGCAACTGCCCGATACCGCTCAGAGTAAGGATATGGTCGGGCTGGTGCAGGGCCAGGGCGGCCTCCGGCATCAGGGCAACCTGTGCGTCACGAGGGTCCTGGACGATCGTGCGCCCCCCATTGTTCTTGATGTAGGCAAGGccc < 1:30255/150‑1 (MQ=255) cccGATACCGCTCAGAGTAAGGATATGGTCGGGCTGGTGCAGGGCCAGGGCGGCCTCCGGCATCAGGGCAACCTGTGCGTCACGAGGGTCCTGGACGATCGTGCGCCCCCCATTGTTCTTGATGTAGGCAAGGCCCCTGGCGCCGTCTTc > 1:362525/1‑150 (MQ=255) gATACCGCTCAGAGTAAGGATATGGTCGGGCTGGTGCAGGGCCAGGGCGGCCTCCGGCATCAGGGCAACCTGTGCGTCACGAGGGTCCTGGACGATCGTGCGCCCCCCATTGTTCTTGATGTAGGCAAGGCCCCTGGCGCCGTCTTCGtt > 2:325820/1‑150 (MQ=255) gATACCGCTCAGAGTAAGGATATGGTCGGGCTGGTGCAGGGCCAGGGCGGCCTCCGGCATCAGGGCAACCTGTGCGTCACGAGGGTCCTGGACGATCGTGCGCCCCCCATTGTTCTTGATGTAGGCAAGGCCCCTGGCGCCGTCTTCGtt < 2:420656/150‑1 (MQ=255) gggCTGGTGCAGGGCCAGGGCGGCCTCCGGCATCAGGGCAACCTGTGCGTCACGAGGGTCCTGGACGATCGTGCGCCCCCCATTGTTCTTGATGTAGGCAAGGCCCCTGGCGCCGTCTTCGTTGGCCCCGGTGAGCAGCACACCcagcag < 2:2611/150‑1 (MQ=255) tGGTGCAGGGCCAGGGCGGCCTCCGGCATCAGGGCAACCTGTGCGTCAAGAGGGTCCTGGACGATCGTGCGCCCCCCATTGTTCTTGATGTAGGCAAGGCCCCTGGCGCCGTCTTCGTTGGCCCCGGTGAGCAGCACACCCAGCAGGCCa > 1:374432/1‑150 (MQ=255) tGGTGCAGGGCCAGGGCGGCCTCCGGCATCAGGGCAACCTGTGCGTCAAGAGGGTCCTGGACGATCGTGCGCCCCCCATTGTTCTTGATGTAGGCAAGGCCCCTGGCGCCGTCTTCGTTGGCCCCGGTGAGCAGCACACCCAGCAGGCCa > 1:374422/1‑150 (MQ=255) gccagggcGGCCTCCGGCATCAGGGCAACCTGTGCGTCACGAGGGTCCTGGACGATCGTGCGCCCCCCATTGTTCTTGATGTAGGCAAGGCCCCTGGCGCCGTCTTCGTTGGCCCCGGTGAGCAGCACACCCAGCAGGCCATCGCCATAg < 1:325820/150‑1 (MQ=255) cagggcGGCCTCCGGCATCAGGGCAACCTGTGCGTCACGAGGGTCCTGGACGATCGTGCGCCCCCCATTGTTCTTGATGTAGGCAAGGCCCCTGGCGCCGTCTTCGTTGGCCCCGGTGAGCAGCACACCCAGCAGGCCATCGCCATAGGc > 1:94034/1‑150 (MQ=255) tCCGGCATCAGGGCAACCTGTGCGTCACGAGGGTCCTGGACGATCGTGCGCCCCCCATTGTTCTTGATGTAGGCAAGGCCCCTGGCGCCGTCTTCGTTGGCCCCGGTGAGCAGCACACCCAGCAGGCCATCGCCATAGGCATCGgccgcc < 2:412592/150‑1 (MQ=255) gCATCAGGGCAACCTGTGCGTCACGAGGGTCCTGGACGATCGTGCGCCCCCCATTGTTCTTGATGTAGGCAAGGCCCCTGGCGCCGTCTTCGTTGGCCCCGGTGAGCAGCACACCCAGCAGGCCATCGCCATAGGCATCGGCCGCCgaga < 1:261536/150‑1 (MQ=255) gCGTCAAGAGGGTCCTGGACGATCGTGCGCCCCCCATTGTTCTTGATGTAGGCAAGGCCCCTGGCGCCGTCTTCGTTGGCCCCGGTGAGCAGCACACCCAGCAGGCCATCGCCATAGGCATCGGCCGCCGAGATGAACAGGAAatcgatc < 2:374422/150‑1 (MQ=255) gCGTCAAGAGGGTCCTGGACGATCGTGCGCCCCCCATTGTTCTTGATGTAGGCAAGGCCCCTGGCGCCGTCTTCGTTGGCCCCGGTGAGCAGCACACCCAGCAGGCCATCGCCATAGGCATCGGCCGCCGAGATGAACAGGAAatcgatc < 2:374432/150‑1 (MQ=255) | TCAGCAGTTTGGCGGTGGTGTGGCTTAGCATGCGCTGTATTCCAGGGTCGCGAGCAACTGCCCGATACCGCTCAGAGTAAGGATATGGTCGGGCTGGTGCAGGGCCAGGGCGGCCTCCGGCATCAGGGCAACCTGTGCGTCACGAGGGTCCTGGACGATCGTGCGCCCCCCATTGTTCTTGATGTAGGCAAGGCCCCTGGCGCCGTCTTCGTTGGCCCCGGTGAGCAGCACACCCAGCAGGCCATCGCCATAGGCATCGGCCGCCGAGATGAACAGGAAATCGATC > NC_002947/4287749‑4288034 |
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: ‑ |