Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 4,512,506 | C→A | 25.0% | G9V (GGT→GTT) | lolA ← | outer membrane lipoprotein carrier protein |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 4,512,506 | 0 | C | A | 25.0% | 34.2 / 8.3 | 20 | G9V (GGT→GTT) | lolA | outer membrane lipoprotein carrier protein |
Reads supporting (aligned to +/- strand): ref base C (9/6); new base A (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 = 8.75e-01 |
TTCCAGCAACTGGGTCAGGCGTTGTACGTCTTGCTCACCGGCATAAGCAGTAACCGAGCCCAGGGTCAGGGCAGAAACCAACAGCATGCGAATCGCGCGCATGGGAATCCTCATTGAGCGTTGAAAAGGCCTGGCGCCACCGTTGGCGCCAGGCAGGGTGTTCATCAGTCGCGCGGGCCGCCCGGGGCTATCACTTCCCGCGAGCCGTTGCTGTTCATCGGGGTGACAACGCCGGCCATCTCCATCGACTCGATCATGCGCGCGGCACGGTTGTAGC > NC_002947/4512367‑4512643 | ttCCAGCAACTGGGTCAGGCGTTGTACGTCTTGCTCACCGGCATAAGCAGTAACCGAGCCCAGGGTCAGGGCAGAAACCAACAGCATGCGAATCGCGCGCATGGGAATCCTCATTGAGCGTTGAAAAGGCCTGGCGCCACCGTTGGCGcc < 1:147837/150‑1 (MQ=255) ccAGCAACTGGGTCAGGCGTTGTACGTCTTGCTCACCGGCATAAGCAGTAACCGAGCCCAGGGTCAGGGCAGAAACCAACAGCATGCGAATCGCGCGCATGGGAATCCTCATTGAGCGTTGAAAAGGCCTGGCGCCACCGTTGGCGCcag > 2:39555/1‑150 (MQ=255) ccAGCAACTGGGTCAGGCGTTGTACGTCTTGCTCACCGGCATAAGCAGTAACCGAGCCCAGGGTCAGGGCAGAAACCAACAGCATGCGAATCGCGCGCATGGGAATCCTCATTGAGCGTTGAAAAGGCCTGGCGCCAACGGTGGCGCcag < 2:123176/150‑1 (MQ=255) tACGTCTTGCTCACCGGCATAAGCAGTAACCGAGCCCAGGGTCAGGGCAGAAACCAACAGCATGCGAATCGCGCGCATGGGAATCCTCATTGAGCGTTGAAAAGGCCTGGCGCCACCGTTGGCGCCAGGCAGGGTGTTCATCAGTcgcgc < 1:39555/150‑1 (MQ=255) aCGTCTTGCTCACCGGCATAAGCAGTAACCGAGCCCAGGGTCAGGGCAGAAACCAACAGCATGCGAATCGCGGGCATGGGAATCCTCATTGAGCGTTGAAAAGGCCTGGCGCCACCGTTGGCGCCAGGCAGGGTGTTCATCAGTcgcgcg < 2:172500/150‑1 (MQ=255) gTCTTGCTCACCGGCATAAGCAGTAACCGAGCCCAGGGTCAGGGCAGAAACCAACAGCATGCGAATCGCGCGCATGGGAATCCTCATTGAGCGTTGAAAAGGCCTGGCGCCACCGTTGGCGCCAGGCAGGGTGTTCATCAGTCGCGCggg > 1:281666/1‑150 (MQ=255) aCCGGCATAAGCAGTAACCGAGCCCAGGGTGAGGGCAGAAACCAACAGCATGCGAATCGCGCGCATGGGAATCCTCATTGAGCGTTGAAAAGGCCTGGCGCCACCGTTGGCGCCAGGCAGGGTGTTCATCAGGCGCGCGGGCCGCCCggg > 1:53727/1‑150 (MQ=255) tAACCGAGCCCAGGGTCAGGGCAGAAACCAACAGCATGCGAATCGCGCGCATGGGAATCCTCATTGAGCGTTGAAAAGGCCTGGCGCCACCGTTGGCGCCAGGCAGGGTGTTCATCAGTCGCGCGGGCCGCCCGGGGCTATCACTTCCcg > 2:420597/1‑150 (MQ=255) ccGAGCCCAGGGTCAGGGCAGAAACCAACAGCATGCGAATCGCGCGCATGGGAATCCTCATTGAGCGTTGAAAAGGCCTGGCGCCAACGGTGGCGCcag > 2:66806/1‑99 (MQ=255) ccGAGCCCAGGGTCAGGGCAGAAACCAACAGCATGCGAATCGCGCGCATGGGAATCCTCATTGAGCGTTGAAAAGGCCTGGCGCCAACGGTGGCGCcag < 1:66806/99‑1 (MQ=255) cGAGCCCAGGGTCAGGGCAGAAACCAACAGCATGCGAATCGCGCGCATGGGAATCCTCATTGAGCGTTGAAAAGGCCTGGCGCCACCGTTGGCGCCAGGCAGGGTGTTCATCAGTCGCGCGGGCCGCCCGGGGCTATCACTTCCCGCGAg > 2:245435/1‑150 (MQ=255) tCAGGGCAGAAACCAACAGCATGCGAATCGCGCGCATGGGAATCCTCATTGAGCGTTGAAAAGGCCTGGCGCCAACGGTGGCGCcag < 1:331543/87‑1 (MQ=255) tCAGGGCAGAAACCAACAGCATGCGAATCGCGCGCATGGGAATCCTCATTGAGCGTTGAAAAGGCCTGGCGCCAACGGTGGCGCcag > 2:331543/1‑87 (MQ=255) gCAGAAACCAACAGCATGCGAATCGCGCGCATGGGAATACTCATTGAGCGTTGAAAAGGCCTGGGGCCACCGTTGGCGCCAGGCAGGGCGTTCATCAGTCGCGCGGGCCGCCCGGGGCTATCACTTCCCGCGAGCCATTGCTGTTCATCg > 2:141729/1‑150 (MQ=255) tCATTGAGCGTTGAAAAGGCCTGGCGCCACCGTTGGCGCCAGGCAGGGTGTTCATCAGTCGCGCGGGCCGCCCGGGGCTATCACTTCCCGCGAGCCGTTGCTGTTCATCGGGGTGACAACGCCGGCCATCTCCATCGACTCGATCATgcg > 1:320228/1‑150 (MQ=255) gAAAAGGCCTGGCGCCACCGTTGGCGCCAGGCAGGGTGTTCATCAGTCGCGCGGGCCGCCCGGGGCTATCACTTCCCGCGAGCCGTTGCTGTTCATCGGGGTGACAACGCCGGCCATCTCCATCGACTCGATCATGCGCGCGGCACGGtt < 2:144408/150‑1 (MQ=255) gAAAAGGCCTGGCGCCACCGTTGGCGCCAGGCAGGGTGTTCATCAGTCGCGCGGGCCGCCCGGGGCTATCACTTCCCGCGAGCCGTTGCTGTTCATCGGGGTGACAACGCCGGCCATCTCCATCGACTCGATCATGCGCGCGGCACGGtt > 1:153096/1‑150 (MQ=255) ggCCTGGCGCCACCGTTGGCGCCAGGCAGGGTGTTCATCAGTCGCGCGGGCCGCCCGGGGCTATCACTTCCCGCGAGCCGTTGCTGTTCATCGGGGTGACAACGCCGGCCATCTCCATCGACTCGATCATGCGCGCGGCACGGTTGTAGc < 2:320228/150‑1 (MQ=255) gcgcCACCGTTGGCGCCAGGCAGGGTGTTCATCAGTCGCGCGGGCCGCCCGGGGCTATCACTTCCCGCGAGCCGTTGCTGTTCATCGGGGt > 2:386150/1‑91 (MQ=255) gcgcCACCGTTGGCGCCAGGCAGGGTGTTCATCAGTCGCGCGGGCCGCCCGGGGCTATCACTTCCCGCGAGCCGTTGCTGTTCATCGGGGt < 1:386150/91‑1 (MQ=255) | TTCCAGCAACTGGGTCAGGCGTTGTACGTCTTGCTCACCGGCATAAGCAGTAACCGAGCCCAGGGTCAGGGCAGAAACCAACAGCATGCGAATCGCGCGCATGGGAATCCTCATTGAGCGTTGAAAAGGCCTGGCGCCACCGTTGGCGCCAGGCAGGGTGTTCATCAGTCGCGCGGGCCGCCCGGGGCTATCACTTCCCGCGAGCCGTTGCTGTTCATCGGGGTGACAACGCCGGCCATCTCCATCGACTCGATCATGCGCGCGGCACGGTTGTAGC > NC_002947/4512367‑4512643 |
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: ‑ |