Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 1,295,019 | Δ1 bp | 10.3% | intergenic (+29/+41) | slyB → / ← nhaA‑I | outer membrane lipoprotein/Na(+)/H(+) antiporter NhaA |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 1,295,019 | 0 | A | . | 10.3% | 106.2 / 7.3 | 39 | intergenic (+29/+41) | slyB/nhaA‑I | outer membrane lipoprotein/Na(+)/H(+) antiporter NhaA |
Reads supporting (aligned to +/- strand): ref base A (13/22); new base . (2/2); total (15/24) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 6.31e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.87e-01 |
TGGCAGCATGCGCGCCTATGTGCAGGCCGTGCAGGAGAACGAAATCTTCCGTGTTGGCGACCGAGTGCGCATCATGACTGTAGACGGTACCAGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCA‑GTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGCGGCTGGCCATTGCGGTGACGGCGTAGCCGATCACCGCAGCAAAGAACGAACCAGTAAGAATGCCCATGCGGTCCATGCCGGCGTA > NC_002947/1294882‑1295164 | tgGCAGCATGCGCGCCTATGTGCAGGCCGTGCAGGAGAACGAAATCTTCCGTGTTGGTGACCGAGTGCGCATCATGACTGTAGACGGTACCAGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCa > 2:474898/1‑150 (MQ=255) ggCNGCATGCGCGCCTATGTGCAGGCCGTGCAGGAGAACGAAATCTTCCGTGTTGGCGACCGAGTGCGCATCATGACTGTAGACGGTACCAGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCA‑g > 2:701310/1‑150 (MQ=255) agcaTGCGCGCCTATGTGCAGGCCGTGCAGGAGAACGAAATGTTCCGTGTTGGCGACCGAGTGCGCATCATGACTGTAGACGGTACCAGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCA‑GTCg > 1:41964/1‑150 (MQ=255) gcaTGCGCGCCTATGTGCAGGCCGTGCAGGAGAACGAAATCTTCCGTGTTGGCGACCGAGTGCGCATCATGACTGTAGACGGTACCAGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAACCCCCGACATGGCTGCTGCCA‑GTCgg < 2:1123622/150‑1 (MQ=255) caTGCGCGCCTATGTGCAGGCCGTGCAGGAGAACGAAATCTTCCGTGTTGGCGACCGAGTGCGCATCATGACTGTAGACGGTACCAGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGAC‑TGGCAGCAGCCATGTCggg < 2:807678/150‑1 (MQ=255) aTGCGCGCCTATGTGCAGGCCGTGCAGGAGAACGAAATCTTCCGTGTTGGCGACCGAGTGCGCATCATGACTGTAGACGGTACCAGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCA‑GTCgggg > 1:188575/1‑150 (MQ=255) gcgcgcCTATGTGCAGGCCGTGCAGGAGAACGAAATCTTCCGTGTTGGCGACCGAGTGCGCATCATGACTGTAGACGGTACCAGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCA‑GTCGGGGtt < 1:145931/150‑1 (MQ=255) cTANGTGCAGGCCGTGCAGGAGAACGAAATCTTCCGTGTTGGCGACCGAGTGCGCATCATGACTGTAGACGGTACCAGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCA‑GTCGGGGTTTTTTCa > 2:702747/1‑150 (MQ=255) tATGTGCAGGCCGTGCAGGAGAACGAAATGTTCCGTGTTGGCGACCGAGTGCGCATCATGACTGTAGACGGTACCAGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCA‑GTCGGGGTTTTTTCAt < 2:41964/150‑1 (MQ=255) gCAGGCCGTGCAGGAGAACGAAATCTTCCGTGTTGGCGACCGAGTGCGCATCATGACTGTAGACGGTACCAGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCA‑GTCGGGGTTTTTTCATGCCTc > 1:507685/1‑150 (MQ=255) gCAGGCCGTGCAGGAGAACGAAATCTTCCGTGTTGGCGACCGAGTGCGCATCATGACTGTAGACGGTACCAGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCA‑GTCGGGGTTTTTTCATGCCTc > 2:35263/1‑150 (MQ=255) aGGCCGTGCAGGAGAACGAAATCTTCCGTGTTGGCGACCGAGTGCGCATCATGACTGTAGACGGTACCAGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCA‑GTCGGGGTTTTTTCATGCCTCAg < 1:108319/150‑1 (MQ=255) aGGCCGTGCAGGAGAACGAAATCTTCCGTGTTGGCGACCGAGTGCGCATCATGACTGTAGACGGTACCAGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGAC‑TGGCAGCAGCCATGTCGGGGtttt < 2:624351/138‑1 (MQ=38) aGGCCGTGCAGGAGAACGAAATCTTCCGTGTTGGCGACCGAGTGCGCATCATGACTGTAGACGGTACCAGTCGCGCTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGAC‑TGGCAGCAGCCATGTCGGGGtttt > 1:624351/1‑138 (MQ=38) ggAGAACGAAATCTTCCGTGTTGGCGACCGAGTGCGCATCATGACTGTAGACGGTACCAGTCGCGTTTCGCACTGACCGTAGTGAAGTAATAATAACCGCCAACATGACACCTCCCT‑GTCAGGATTTGTTCAGGCCTCAGCCAAGTCAgg < 2:852206/150‑1 (MQ=255) aaaTCTTCCGTGTTGGCGACCGAGTGCGCATCATGACTGTAGACGGTACCAGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCA‑GTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTg < 1:701310/150‑1 (MQ=255) aaaTCTTCCGTGTTGGCGACCGAGTGCGAATCATGACTGTAGACGGTACCAGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACACGGGCGCTTGCA‑GTCGGGTTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTg < 1:1041467/150‑1 (MQ=255) ccGTGTTGGCGACCGAGTGCGCATCATGACTGTAGACGGTACCAGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCA‑GTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGc < 2:771968/150‑1 (MQ=255) gcgcATCATGACTGTAGACGGTACCAGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCA‑GTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGCGGCTGGCCATTGCGGTGa < 1:853249/150‑1 (MQ=255) gcATCATGACTGTAGACGGTACCAGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCA‑GTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGCGGCTGACCATTGCGGTTAct < 2:545450/150‑2 (MQ=255) gACTGTAGACGGTACCAGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCA‑GTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGCGGCTGGCCATTGCGGTGACGGCGTAGc < 2:188575/150‑1 (MQ=255) aGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCA‑GTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGCGGCTGGCCATTGCGGTGACGGCGTAGCCGATCACCGCAGCaaa > 1:649650/1‑150 (MQ=255) tttCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCA‑GTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGCGGCTGGCCATTGCGGTGACGGCGTAGCCGATCACCGCAGCAAAgaacgaa < 2:314872/150‑1 (MQ=255) ttCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCA‑GTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGCGGCTGGCCATTGCGGTGACGGCGTAGCCGATCACCGCAGCAAAgaacgaac < 1:35263/150‑1 (MQ=255) cGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCA‑GTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGCGGCTGGCCATTGCGGTGACGGCGTAGCCGATCACCGCAGCAAAGAACGAACCa < 1:474898/150‑1 (MQ=255) gCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCA‑GTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGCGGCTGGCCATTGCGGTGACGGCGTAGCCGATCACCGCAGCAAAGAACGAACCAg < 1:24699/150‑1 (MQ=255) tGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCA‑GTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGCGGCTGGCCATTGCGGTGACGGCGTAGCCGATCACCgcag > 2:725489/1‑131 (MQ=255) tGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCA‑GTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGCGGCTGGCCATTGCGGTGACGGCGTAGCCGATCACCgcag < 1:725489/131‑1 (MQ=255) gTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCA‑GTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGCGGCTGGCCATTGCGGTGACGGCGTAGCCGATCACCGCAGCAAAGAACGAACCAGTAAGAATGc < 2:649650/150‑1 (MQ=255) taataaAAAACCCCGACATGGCTGCTGCCA‑GTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGCGGCTGGCCATTGCGGTGACGGCGTAGCCGATCAcc < 1:109905/113‑1 (MQ=255) taataaAAAACCCCGACATGGCTGCTGCCA‑GTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGCGGCTGGCCATTGCGGTGACGGCGTAGCCGATCAcc > 2:109905/1‑113 (MQ=255) ataaAAAACCCCGACATGGCTGCTGCCA‑GTCgggg > 2:543343/1‑35 (MQ=255) ataaAAAACCCCGACATGGCTGCTGCCA‑GTCgggg < 1:543343/35‑1 (MQ=255) taaAAAACCCCGACATGGCTGCTGCCA‑GTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGCGGCTGGCCATTGCGGTGACGGCGTAGCCGATCACCGCAGCAAAGAACGAACCAGTAAGAATGCCCATGCGGTCCa < 2:507685/150‑1 (MQ=255) aaaaCCCCGACATGGCTGCTGCCA‑GTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGCGGCTGGCCATTGCGGTGACGGCGTAGCCGATCACCGCAGCAAAGAACGAACCAGTAAGAATGCCCATGCGGTCCATGc > 2:725001/1‑150 (MQ=255) aaCCCCGACATGGCTGCTGCCA‑GTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGCGGCTGGCCATTGCGGTGACGGCGTAGCCGATCACCGCAGCAAAGAACGAACCAGTAAGAATGCCCATGCGGTCCATGCCg < 1:725001/150‑1 (MQ=255) aCCCCGAC‑TGGCAGCAGCCATGTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGCGGCTGGCCATTGCGGTGACGGCGTAGCCGATCACCGCAGCAAAGAACGAACCAGTAAGAATGCCCATGCGGTCCATGCCgg > 1:482410/1‑150 (MQ=255) cGACATGGCTGCTGCCA‑GTCGGGGTTTTTTCATGCCTCAGCCACGCCAGGCAATGCTGGTTTTGCGGCTGGCCATTGCGGTGACGGCGTAGCCGATCACCGCAGCAAAGAACGAACCAGTAAGAATGCCCATGCGGTCCATGCCGGCGTa < 2:760243/150‑1 (MQ=255) cGACATGGCTGCTGCCA‑GTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGCGGCTGGCCATTGCGGTGACGGCGTAGCCGATCACCGCAGCAAAGAACGAACCAGTAAGAATGCCCATGCGGTCCATGCCGGCGTa > 1:40942/1‑150 (MQ=255) | TGGCAGCATGCGCGCCTATGTGCAGGCCGTGCAGGAGAACGAAATCTTCCGTGTTGGCGACCGAGTGCGCATCATGACTGTAGACGGTACCAGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCA‑GTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGCGGCTGGCCATTGCGGTGACGGCGTAGCCGATCACCGCAGCAAAGAACGAACCAGTAAGAATGCCCATGCGGTCCATGCCGGCGTA > NC_002947/1294882‑1295164 |
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: ‑ |