Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 1,506,932 | A→G | 24.2% | intergenic (+44/‑73) | petC → / → sspA | ubiquinol‑‑cytochrome c reductase cytochrome c1 subunit/stringent starvation protein A |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 1,506,932 | 0 | A | G | 24.2% | 52.1 / 14.8 | 33 | intergenic (+44/‑73) | petC/sspA | ubiquinol‑‑cytochrome c reductase cytochrome c1 subunit/stringent starvation protein A |
Reads supporting (aligned to +/- strand): ref base A (11/14); new base G (4/4); total (15/18) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.00e-01 |
AACTGGAAAGCCAGCGCATCGGCACCTATGTGTTGCTGTATCTGGCTTTCTTCTTCGTATTCGCCTACTTGCTCAAGCGTGAATACTGGAAGGACGTGCACTGATCACGCAGTAGTTTCTGGTTGTTGAACGCGCCCGGGGTGGCCTAGCAGTATTGCCAGGCCACCCCGGGCGCGTTTTCATTTAAAGTTGCACAAGCATCCGTGCGAGGAGCCGCTACATGGGCGCAACCAACAGGTTAGCCTGCTATTCCGATCCCGCTGATCACTATTCTCATCGGGTTCGCCTGGTTCTG > NC_002947/1506785‑1507079 | agcTGGAAAGCCGGCGCTTCGGCACCTATGTTTTGCTGTATCTGGCTTTCTTCTTCGTATTCGCCTACTTGCTCAAGCGTGAATACTGGGAGGACGTGCACTGATCACGCAGTAGTTTCTGGTTGTTGAACGCGCCCGGGGTGGCGTAGc < 2:1705124/148‑1 (MQ=255) aCTGGAAAGCCAGCGCATCGGCACCTATGTGTTGCTGTATCTGGCTTTCTTCTTCGTATTCGCCTACTTGCTCAAGCGTGAATACTGGAAGGACGTGCACTGATCACGCAGTAGTTTCTGGTTGTTGAACGCGCCCGGGGTGGCCTAGCa < 2:1145018/150‑1 (MQ=255) cTGGAAAGCCAGCGCATCGGCACCTATGTGTTGCTGTATCTGGCTTTCTTCTTCGTATTCGCCTACTTGCTCAAGCGTGAATACTGGAAGGACGTGCACTGATCACGCAGTAGTTTCTGGTTGTTGAACGCGCCCGGGGTGGCCTAGCAg < 1:159343/150‑1 (MQ=255) aGCCAGCGCATCGGCACCTATGTGTTGCTGTATCTGGCTTTCTTCTTCGTATTCGCCTACTTGCTCAAGCGTGAATACTGGAAGGACGTGCACTGATCACGCAGTAGTTTCTGGTTGTTGAACGCGCCCGGGGTGGCCTAGCAGTATTGc < 2:341502/150‑1 (MQ=255) aGCCAGCGCATCGGCACCTATGTGTTGCTGTATCTGGCTTTCTTCTTCGTATTCGCCTACTTGCTCAAGCGTGAATACTGGAAGGACGTGCACTGATCACGCAGTAGTTTCTGGTTGTTGAACGCGCCCGGGGTGGCCTAGCAGTATTGc < 1:865167/150‑1 (MQ=255) ccAGCGCATCGGCACCTATGTGTTGCTGTATCTGGCTTTCTTCTTCGTATTCGCCTACTTGCTCAAGCGTGAATACTGGAAGGACGTGCACTGATCACGCAGTAGTTTCTGGTTGTTGAACGCGCCCGGGGTGGCCTAGCAGTATTGCCa < 1:1234738/150‑1 (MQ=255) aGCGCATCGGCACCTATGTGTTGCTGTATCTGGCTTTCTTCTTCGTATTCGCCTACTTGCTCAAGCGTGAATACTGGGAGGACGTGCACTGATCACGCAGTAGTTTCTGGTTGTTGAACGCGCCCGGGGTGGCCTAGCAGTATTGCCAgg > 1:1498281/1‑150 (MQ=255) gcgcATCGGCACCTATGTGTTGCTGTATCTGGCTTTCTTCTTCGTATTCGCCTACTTGCTCAAGCGTGAATACTGGAAGGACGTGCACTGATCACGCAGTAGTTTCTGGTTGTTGAACGCGCCCGGGGTGGCCTGGCAATACTGCTAGGc < 2:1703288/150‑1 (MQ=255) gcgcATCGGCACCTATGTGTTGCTGTATCTGGCTTTCTTCTTCGTATTCGCCTACTTGCTCAAGCGTGAATACTGGAAGGACGTGCACTGATCACGCAGTAGTTTCTGGTTGTTGAACGCGCCCGGGGTGGCCTAGCAGTATCGCCAGGc < 2:25754/150‑1 (MQ=255) tCGGCACCTATGTGTTGCTGTATCTGGCTTTCTTCTTCGTATTCGCCTACTTGCTCAAGCGTGAATACTGGAAGGACGTGCACTGATCACGCAGTAGTTTCTGGTTGTTGAACGCGCCCGGGGTGGCCTAGCAGTATTGCCAGGCCAccc < 2:1370852/150‑1 (MQ=255) gtTGCTGTATCTGGCTTTCTTCTCCGTATTCGCCTACTTGCTCAAGCGTGAATACTGGAAGGACGTGCACTGATCACGCAGTAGTTTCTGGTTGTTGAACGCGCCCGGGGTGGCCTAGCAGTATTGCCAGGCCACCCCGGGCGCGTTTTc > 2:946580/1‑150 (MQ=255) tATCTGGCTTTCTTCTTCGTATTCGCCTACTTGCTCAAGCGTGAATACTGGAAGGACGTGCACTGATCACGCAGTAGTTTCTGGTTGTTGAACGCGCCCGGGGTGGCCTAGCAGTATTGCCAGGc < 1:279780/125‑1 (MQ=255) tATCTGGCTTTCTTCTTCGTATTCGCCTACTTGCTCAAGCGTGAATACTGGAAGGACGTGCACTGATCACGCAGTAGTTTCTGGTTGTTGAACGCGCCCGGGGTGGCCTAGCAGTATTGCCAGGc > 2:279780/1‑125 (MQ=255) tttCTTCTTCGTATTCGCCTACTTGCTCAAGCGTGAATACTGGAAGGACGTGCACTGATCACGCAGTAGTTTCTGGTTGTTGAACGCGCCCGGGGTGGCCTAGCAGTATTGCCAGGCCACCCCGGGCGCGTTTTCATTTAGAGTTGcaca > 1:1061535/1‑150 (MQ=255) tATTCGCCTACTTGCTCAAGCGTGAATACTGGAAGGACGTGCACTGATCACGCAGTAGTTTCTGGTTGTTGAACGCGCCCGGGGTGGCCTGGCAATac < 2:954258/98‑2 (MQ=255) tATTCGCCTACTTGCTCAAGCGTGAATACTGGAAGGACGTGCACTGATCACGCAGTAGTTTCTGGTTGTTGAACGCGCCCGGGGTGGCCTGGCAATac > 1:954258/1‑97 (MQ=255) cTCAAGCGTGAATACTGGAAGGACGTGCACTGATCACGCAGTAGTTTCTGGTTGTTGAACGCGCCCGGGGTGGCCTAGCAGTATTGCCAGGCCACCCCGGGCGCGTTTTCATTTAAAGTTGCACAAGCATCCGTGCGAGGAGCCGCTACa < 1:946580/150‑1 (MQ=255) aaGCGTGAATACTGGAAGGCCGTGCACTGATCACGCGGGAGTTTCTGGTTGTTGAACGCGCCCGGGGTGGCCTAGCAGTATTGCCAGGCCACCCCGGGCGCGTTTTCATTTAAAGTTGCACAAGCATCCGTGCGAGGAGCCGCTACATgg > 2:1755880/1‑150 (MQ=255) agtagtTTCTGGTTGTTGAACGCGCCCGGGGTGGCCTAGCAGTATTGCCAGGCCACCCCGGGCGCGTTTTCATTTAAAGTTGCACAAGCATCCGTGCGAGGAGCCGCTACATGGGCGcaac > 2:901234/1‑121 (MQ=255) agtagtTTCTGGTTGTTGAACGCGCCCGGGGTGGCCTAGCAGTATTGCCAGGCCACCCCGGGCGCGTTTTCATTTAAAGTTGCACAAGCATCCGTGCGAGGAGCCGCTACATGGGCGcaac < 1:901234/121‑1 (MQ=255) gtagtTTCTGGTTGTTGAACGCGCCCGGGGTGGCCTAGCAGTATTGCCAGGCCACCCCGGGCGCGTTTTCATTTAAAGTTGCACAAGCATCCGTGCGAGGAGCCGCTACATGGGCGCAAACAAAAGATTAGCCTGCTATTCCGATCCCGc > 2:107790/1‑150 (MQ=255) ggTTGTTGAACGCGCCCGGGGTGGCCTAGCAGTATTGCCAGGCCCCCCCGGGCGCGTTTTCATTTAAAGTTGCACAAGCATCCGTGCGAGGAGCCGCTACATGGGCGCAACCAACAGGTTAGCCTGCTATTCCGATCCCGCTGATCACTa < 2:1498281/150‑1 (MQ=255) aCGCGCCCGGGGTGGCCTAGCAGTATTGCCAGGCCACCCCGGGCGCGTTTTCATTTAAAGTTGCACAAGCATCCGTGCGAGGAGCCGCTACATGGGCGCAACCAACAGGTTAGCCTGCTATTCCGATCCCGCTGATCACTATTCTCATCg < 1:107790/150‑1 (MQ=255) gcCCGGGGTGGCCTGGCAATACTGCTAGGCCACCCCGGGCGCGTTTTCATTTAAAGTTGCACAAGCATCCGTGCGAGGAGCCGCTACATGGGCGCAATCAACAGGTTAGCCTGCTATTCCGATCCCGCTGATCACTa > 2:1073849/1‑137 (MQ=37) gcCCGGGGTGGCCTGGCAATACTGCTAGGCCACCCCGGGCGCGTTTTCATTTAAAGTTGCACAAGCATCCGTGCGAGGAGCCGCTACATGGGCGCAATCAACAGGTTAGCCTGCTATTCCGATCCCGCTGATCACTa < 1:1073849/137‑1 (MQ=37) cccGGGGTGGCCTAGCAGTATTGCCAGGCCACCCCGGGCGCGTTTTCATTTAAAGTTGCACAAGCATCCGTGCGAGGAGCCGCTACATGGGCGCAACCAACAg < 1:765886/103‑1 (MQ=255) cccGGGGTGGCCTAGCAGTATTGCCAGGCCACCCCGGGCGCGTTTTCATTTAAAGTTGCACAAGCATCCGTGCGAGGAGCCGCTACATGGGCGCAACCAACAg > 2:765886/1‑103 (MQ=255) ccGGGGTGGCCTGGCAATACTGCTAGGCCACCCCGGGCGCGTTTTCATTTAAAGTTGCACAAGCATCCGTGCGAGGAGCCGCTACATGGGCGCAACCAACAGGTTAGCCTGCTATTCCGATCCCGCTGATCACTATTCTCATCGGGTTCg > 2:176364/1‑150 (MQ=39) ccGGGGTGGCCTAGCAGTATTGCCAGGCCACCCCGGGCGCGTTTTCATTTAAAGTTGCACAAGCATCCGTGCGAGGAGCCGCTACATGGGCGCAACCAACAGGTTAGCCTGCTATTCCGATCCCGCTGATCACTATTCTCATCGGGTTCg > 2:1659485/1‑150 (MQ=255) gTGGCCTGGCAATACTGCTAGGCCACCCCGGGCGCGTTTTCATTTAAAGTTGCACAAGCATCCGTGCGAGGAGCCGCTACATGGGCGCAACCAACa > 2:314863/1‑96 (MQ=25) gTGGCCTGGCAATACTGCTAGGCCACCCCGGGCGCGTTTTCATTTAAAGTTGCACAAGCATCCGTGCGAGGAGCCGCTACATGGGCGCAACCAACa < 1:314863/96‑1 (MQ=25) ggCCTAGCAGTATTGCCAGGCCACCCCGGGCGCGTTTTCATTTAAAGTTGCACAAGCATCCGTGCGAGGAGCCGCTACATGGGCGCAACCAACAGGTTAGCCTGCTATTCCGATCCCGCTGAGCACTATTCTCATCGGGTTCGCCTGGtt > 2:1680015/1‑150 (MQ=255) cTAGCAGTATTGCCAGGCCACCCCGGGCGCGTTTTCATTTAAAGTTGCACAAGCATCCGTGCGAGGAGCCGCTACATGGGCGCAACCAACAGGTTAGCCTGCTATTCCGATCCCGCTGATCACTATTCTCATCGGGTTCGCCTGGTTCTg > 2:1554666/1‑150 (MQ=255) | AACTGGAAAGCCAGCGCATCGGCACCTATGTGTTGCTGTATCTGGCTTTCTTCTTCGTATTCGCCTACTTGCTCAAGCGTGAATACTGGAAGGACGTGCACTGATCACGCAGTAGTTTCTGGTTGTTGAACGCGCCCGGGGTGGCCTAGCAGTATTGCCAGGCCACCCCGGGCGCGTTTTCATTTAAAGTTGCACAAGCATCCGTGCGAGGAGCCGCTACATGGGCGCAACCAACAGGTTAGCCTGCTATTCCGATCCCGCTGATCACTATTCTCATCGGGTTCGCCTGGTTCTG > NC_002947/1506785‑1507079 |
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: ‑ |