Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 5,997,245:1 | +G | 25.0% | intergenic (‑158/+42) | PP_5252 ← / ← PP_5253 | hypothetical protein/arylesterase |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 5,997,245 | 1 | . | G | 25.0% | 36.4 / 18.5 | 20 | intergenic (‑158/+42) | PP_5252/PP_5253 | hypothetical protein/arylesterase |
Reads supporting (aligned to +/- strand): ref base . (9/6); new base G (2/3); total (11/9) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 6.17e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.87e-01 |
GCAGGGGGACTGCGAAGCATGCCGGTTGATGCAAGGCCTGGGCTTGAATGGATGTGCGGTGGGTTGAAGAAGGTTAGCCCTGGTCAAGCTTTTTGCTGCCTGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACA‑ACGTCTGTAGGAGCGGGTTTACCCGCGAATGTTTCCCACAATTACCGCTGAAGGAACGCCAGCAAGTCCTCATTCAGCTGCTGCGCATGAGTCACTGCAAAACCGTGCGGCGCCCCGGCATACACCTTCAGCTCGGCGCCACGAATCAG > NC_002947/5997110‑5997394 | gCAGGGGGACTGCGAAGCATGCCGGTTGATGCAAGGCCTGGGCTTGAATGGATGTGCGGTGGGTTGAAGAAGGTTAGCCCTGGTCAAGCTTTTTGCTGCCTGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGACGT‑TGTAGGAGc < 2:499572/150‑1 (MQ=255) gCAGGGGGACTGCGAAGCATGCCGGTTGATGCAAGGCCTGGGCTTGAATGGATGTGCGGTGGGTTGAAGAAGGTTAGCCCTGGTCAAGCTTTTTGCTGCCTGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACA‑ACGTCTGTAGGAGc < 1:185465/150‑1 (MQ=255) tGCGAAGCATGCCGGTTGATGCAAGGCCTGGGCTTGAATGGATGTGCGGTGGGTTGAAGAAGGTTAGCCCTGGTCAAGCTTTTTGCTGCCTGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGACGT‑TGTAGGAGc < 1:252625/140‑1 (MQ=255) tGCGAAGCATGCCGGTTGATGCAAGGCCTGGGCTTGAATGGATGTGCGGTGGGTTGAAGAAGGTTAGCCCTGGTCAAGCTTTTTGCTGCCTGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGACGT‑TGTAGGAGc > 2:252625/1‑140 (MQ=255) aaGGCCTGGGCTTGAATGGATGTGCGGTGGGTTGAAGAAGGTTAGCCCTGGTCAAGCTTTTTGCTGCCTGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACA‑ACGTCTGTAGGAGCGGGTTTACCCGCGAATGTTTCCCACAATTAcc < 1:62519/150‑1 (MQ=255) tGGATGTGCGGGGGGTTGAAGAAGGCTAGCCCTGGTCAAGCTTTTTGCTGCCTGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACA‑ACGTCTGTAGGAGCGGGTTTACCCGCGAATGTTTCCCACAATTACCGCTGAAGGAACGCCAg < 2:192350/150‑1 (MQ=255) gtgCGGTGGGTTGAAGAAGGTTAGCCCTGGTCAAGCTTTTTGCTGCCTGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACA‑ACGTCTGTAGGAGCGGGTTTACCCGCGAATGTTTCCCACAATTACCGCTGAAGGAACGCCAGCAAGt > 1:145410/1‑150 (MQ=255) tGGTCAAGCTTTTTGCTGCCTGGGCCGNCCTCTTCGCGGGNAANCCCNCTCCTACA‑ACGTCTGTAGAAGCGGGTTTACCCGCGAATGTTTCCCACAATTACCGCTGAAGGAACGCCAGCAAGTCCTCATTCAGCTGCTGCGCATGAGTCa > 1:502270/1‑150 (MQ=255) gTCAAGCTTTTTGCTGCCTGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACA‑ACGTCTGTAGGAGCGGGTTTACCCGCGAATGTTTCCCACAATTACCGCTGAAGGAACGCCAGCAAGTCCTCATTCAGCTGCTGCGCATGAGTCACt > 2:246561/1‑150 (MQ=255) cggccTCTTCGCGGGTAAACCCGCTCCTACA‑ACGTCTGTAGGAGCGGGTTTACCCGCGAATGTTTCCCACAATTACCGCTGAAGGAACGCCAGCAAGTCCTCATTCAGCTGCTGCGCATGAGTCACTGCAAAACCGTGCGGCGCCCCGGc > 2:399992/1‑150 (MQ=255) ttCGCGGGTAAACCCGCTCCTACA‑ACGTCTGTAGGAGCGGGTTTACCCGCGAATGTTTCCCACAATTACCGCTGAAGGAACGCCAGCAAGTCCTCATTCAGCTGCTGCGCATGAGTCACTGCAAAACCGTGCGGCGCCCCGGCATACAcc > 2:244117/1‑150 (MQ=255) gggTAAACCCGCTCCTACA‑ACGTCTGTAGGAGCGGGTTTACCCGCGAATGTTTCCCACAATTACCGCTGAAGGAACGCCAGCAAGTCCTCATTCAGCTGCTGCGCATGAGTCACTGCAAAACCGTGCGGCGCCCCGGCATACACCTTCAg > 1:385337/1‑150 (MQ=255) ggTAAACCCGCTCCTACA‑ACGTCTGTAGGAGCGGGTTTACCCGCGAATGTTTCCCACAATTACCGCTGAAGGAACGCCAg < 2:506868/80‑1 (MQ=255) ggTAAACCCGCTCCTACA‑ACGTCTGTAGGAGCGGGTTTACCCGCGAATGTTTCCCACAATTACCGCTGAAGGAACGCCAg > 1:506868/1‑80 (MQ=255) cccGCTCCTACAGACGT‑TGTAGGAGCGGGTTTACCCGCGAATGTTTCCCACAATTACCGCTGAAGGAACGCCAGCAAGTCCTCATTCAGCTGCTGCGCATGAGTCACTGCAAAACCGTGCGGCGCCCCGGCATACACCTTCAGCTCgg < 2:359775/148‑1 (MQ=255) cccGCTCCTACAGACGT‑TGTAGGAGCGGGTTTACCCGCGAATGTTTCCCACAATTACCGCTGAAGGAACGCCAGCAAGTCCTCATTCAGCTGCTGCGCATGAGTCACTGCAAAACCGTGCGGCGCCCCGGCATACACCTTCAGCTCgg > 1:359775/1‑148 (MQ=255) cccGCTCCTACA‑ACGTCTGTAGGAGCGGGTTTACCCGCGAATTTTTCCCACAATTACCGCTGAAGGAACGCCAGCAAGTCCTCATTCAgctgc < 2:368563/93‑1 (MQ=255) cccGCTCCTACA‑ACGTCTGTAGGAGCGGGTTTACCCGCGAATTTTTCCCACAATTACCGCTGAAGGAACGCCAGCAAGTCCTCATTCAgctgc > 1:368563/1‑93 (MQ=255) tACA‑ACGTCTGTAGGAGCGGGTTTACCCGCGAATGTTTCCCACAATTACCGCTGAAGGAACGCCAGCAAGTCCTCATTCAGCTGCTGCGCATGAGTCACt < 2:1570/100‑1 (MQ=255) tACA‑ACGTCTGTAGGAGCGGGTTTACCCGCGAATGTTTCCCACAATTACCGCTGAAGGAACGCCAGCAAGTCCTCATTCAGCTGCTGCGCATGAGTCACt > 1:1570/1‑100 (MQ=255) a‑aCGTCTGTAGGAGCGGGTTTACCCGCGAATGTTTCCCACAATTACCGCTGAAGGAACGCCAGCAAGTCCTCATTCAGCTGCTGCGCATGAGTCACTGCAAAACCGTGCGGCGCCCCGGCATACACCTTCAGCTCGGCGCCACGAATCAg < 1:399992/150‑1 (MQ=255) | GCAGGGGGACTGCGAAGCATGCCGGTTGATGCAAGGCCTGGGCTTGAATGGATGTGCGGTGGGTTGAAGAAGGTTAGCCCTGGTCAAGCTTTTTGCTGCCTGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACA‑ACGTCTGTAGGAGCGGGTTTACCCGCGAATGTTTCCCACAATTACCGCTGAAGGAACGCCAGCAAGTCCTCATTCAGCTGCTGCGCATGAGTCACTGCAAAACCGTGCGGCGCCCCGGCATACACCTTCAGCTCGGCGCCACGAATCAG > NC_002947/5997110‑5997394 |
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: ‑ |