| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | NC_002947 | 1,506,927 | G→A | 17.6% | intergenic (+48/‑69) | 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,927 | 0 | G | A | 17.6% | 72.3 / 8.3 | 34 | intergenic (+48/‑69) | petC/sspA | ubiquinol‑‑cytochrome c reductase cytochrome c1 subunit/stringent starvation protein A |
| Reads supporting (aligned to +/- strand): ref base G (14/14); new base A (4/2); total (18/16) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 6.60e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.17e-01 | |||||||||||
AAGCCAGCGCATCGGCACCTATGTGTTGCTGTATCTGGCTTTCTTCTTCGTATTCGCCTACTTGCTCAAGCGTGAATACTGGAAGGACGTGCACTGATCACGCAGTAGTTTCTGGTTGTTGAACGCGCCCGGGGTGGCCTAGCAGTATTGCCAGGCCACCCCGGG‑CGCGTTTTCATTTAAAGTTGCACAAGCATCCGTGCGAGGAGCCGCTACATGGGCGCAACCAACAGGTTAGCCTGCTATTCCGATCCCGCTGATCACTATTCTCATCGGGTTCGCCTGGTTCTGGCC > NC_002947/1506783‑1507073 | aaGCCAGCGCATCGGCACCTATGTGTTGCTGTATCTGGCTTTCTTCTTCGTATTCGCCTACTTGCTCAAGCGTGAATACTGGAAGGACGTGCACTGATCACGCAGTAGTTTCTGGTTGTTGAACGCGCCCGGGGTGGCCTAGCAGTAtt < 6:268608/149‑1 (MQ=255) cgcATCGTCACGTATGTGTTGCTGTATTTGGCTTTCTTCTTCGTATTCGCCTACTTGCTCTAGAGTGAATGCTGGAAGGACTTGCACTGATCACGCAGTAGTTTCTGGTTGTTGAACGCGCCCGTGGTGGACTAGCAGTATTGCAAGGc < 3:188570/149‑1 (MQ=255) cgcATCGGCACCTATGTGTTGCTGTATCTGGCTTTCTTCTTCGTATTCGCCTACTTGCTCAAGCGTGAATACTGGAAGGACGTGCACTGATCACGCAGTAGTTTCTGGTTGTTGAACGCGCCCGGGGTGGCCTAGCAGTATTGCCAGGc < 7:25938/149‑1 (MQ=255) cgcATCGGCACCTATGTGTTGCTGTATCTGGCTTTCTTCTTCGTATTCGCCTACTTGCTCAAGCGTGAATACTGGAAGGACGTGCACTGATCACGCAGTAGTTTCTGGTTGTTGAACGCGCCCGGGGTGGCCTAGCAGTATTGCCAGGc < 6:109806/149‑1 (MQ=255) gCACCTATGTGTTGCTGTATCTGGCTTTCTTCTTCGTATTCGCCTACTTGCTCAAGCGTGAATACTGGAAGGACGTGCACTGATCACGCAGTAGTTTCTGGTTGTTGAACGCGCCCGGGGTGGCCTAGCAGTATTGCCAGGCCACCCCg > 8:113443/1‑149 (MQ=255) tATGTGTTGCTGTATCTGGCTTTCTTCTTCGTATTCGCCTACTTGCTCAAGCGTGAATACTGGAAGGACGTGCACTGATCACGCAGTAGTTTCTGGTTGTTGAACGCGCCCGGGGTGGCCTAGCAGTATTGCCAGGCCACCCCGGg‑cgc > 7:183825/1‑149 (MQ=255) tgtgTTGCTGTATCTGGCTTTCTTCTTCGTATTCGCCTACTTGCTCATGCGTGAATACTGGAAGGACGTGCACTGATCACGCAGTAGTTTCTGGTTGTTGAACGCGCCCGGGGTGGCCTAGCAGTATTGCCAGGCCACCCCGGG‑CGCGt < 7:113443/149‑1 (MQ=255) ttGCTGTAGCTGGCTTTCTTCTTCGTATTCGCCTACTTGCGCAAGCGGGAATACTGGAAGGACGTGCACTGATCACGCAGTAGTTTCTGGTTGTTGAACGCGCCCGGGGTGGCCTAGCAGTATTGCCAGGCCACCCCGGG‑CGCGTTTTc < 7:181840/149‑1 (MQ=255) tCGCCTACTTGCTCAAGCGTGAATACTGGAAGGACGTGCACTGATCACGCAGTAGTTTCTGGTTGTTGAACGCGCCCGGGGTGGCCTAGCAGTATTGCCAGGCCACCCCGGG‑CGCGTTTTCATTTAAAGTTGCACAAGCATCCGTGc > 2:30922/1‑147 (MQ=255) ccTACTTGCTCAAGCGTGAATACTGGAAGGACGTGCACTGATCACGCAGTAGTTTCTGGTTGTTGAACGCGCCCGGGGTGGCCTAGCAGTATTGCCAGGCCACCCCGGG‑CGCGTTTTCTTTTAAAGTTGCACAAGCATCCGTGCgagga < 1:30922/149‑1 (MQ=255) ccTACTTGCTCAAGCGTGAATACTGGAAGGACGTGCACTGATCACGCAGTAGTTTCTGGTTGTTGAACGCGCCCGGGGTGGCCTAGCAGTATTGCCAGGCCACCCCGGG‑CGCGTTTTCATTTAAAGTTGCACAAGCATCCGTGCgagga < 8:183825/149‑1 (MQ=255) cTACTTGCTCAAGCGTGAATACTGGAAGGACGTGCACTGATCACGCAGTAGTTTCTGGTTGTTGAACGCGCCCGGGGTGGCCTAGCAGTATTGCCAGGCCACCCCGGG‑CGCGTTTTCATTTAAAGTTGCACAAGCATCCGTGCgaggag > 4:139692/1‑149 (MQ=255) aaGCGTGAATACTGGAAGGACGTGCACTGATCACGCAGTAGTTTCTGGTTGTTGAACGCGCCCGGGGTGGCCTGGCAATACTGCTAGGc < 1:1257/89‑1 (MQ=37) aaGCGTGAATACTGGAAGGACGTGCACTGATCACGCAGTAGTTTCTGGTTGTTGAACGCGCCCGGGGTGGCCTGGCAATACTGCTAGGc > 2:1257/1‑89 (MQ=37) ggAAGGACGTGCACTGATCACGCAGTAGTTTCTGGTTGTTGAACGCGCCCGGGGTGGCCTAGCAGTATTGCCAGGCCACCCCGGG‑CGCGTTTTCATTTAAAGTTGCACAAGCATCCGTGCGAGGAGCCGCTACATGGGCGcaaccaa > 7:181702/1‑147 (MQ=255) aTCACGCAGTAGTTTCTGGTTGTTGAACGCGCCCGGGGTGGCCTAGCAGTATTGCCAGGCCACCCCGGG‑CGCGTTTTCATTTAAAGTTGCACAAGCATCCGTGCGAGGAGCCGCTACATGGgc < 4:335282/123‑1 (MQ=255) aTCACGCAGTAGTTTCTGGTTGTTGAACGCGCCCGGGGTGGCCTAGCAGTATTGCCAGGCCACCCCGGG‑CGCGTTTTCATTTAAAGTTGCACAAGCATCCGTGCGAGGAGCCGCTACATGGgc > 3:335282/1‑123 (MQ=255) ggTTGTTGAACGCGCCCGGGGTGGCCTAGCAGTATTGCCAGGCCACCCCGGG‑CGCGTTTTCATTTAAAGTTGCACAAGCATCCGTGCGAGGAGCCGCTACATGGGCGCAACCAACAGGTTAGCCTGCTATTCCGATCCCGCTGATCACt > 3:184814/1‑149 (MQ=255) gttgAACGCGCCCGGGGTGGCCTAGCAGTATTGCCAGGCCACCCCGGG‑CGCGTTTTCATTTAAAGTTGCACAAGCATCCGTGCGAGGAGCCGCTACATGGGCGCAACCAACAGGTTAGCCTGCTATTCCGATCCCGCTGATCACTAtt < 8:181702/148‑1 (MQ=255) cgcgCCAGGGGTGGCCTAGAAGAATTGCCAGGCCACCCCGGG‑CGCCTTTTCATTTAAAGTTGCACAAGCAAACGTGCGAGGCGCCCCTACATGGGCGCAACCAACAGGATAGCCACCTATTCCCATCCCGATGATCACTAATATCATCg > 8:49008/1‑149 (MQ=255) gcgcCCGGGGTGGCCTAGCAGTATTGCCAGGCCACCCCGGG‑CGCGTTTTCATTTAAAGTTGCACAAGCATCCGt < 7:153315/74‑1 (MQ=255) gcgcCCGGGGTGGCCTAGCAGTATTGCCAGGCCACCCCGGG‑CGCGTTTTCATTTAAAGTTGCACAAGCATCCGt > 8:153315/1‑74 (MQ=255) gcCCGGGGTGGCCTAGCAGTATTGCCAGGCCACCCCGGG‑CGCGTTTTCATTTAAAGTTGCACAAGCATCCGTGCGAGGAGCCGCTACATGGGCGCAACCAACAGGTTAGCCTGCTATTCCGATCCCGCTGATCACTATTCTCATCGGGt > 7:112979/1‑149 (MQ=255) cccGGGGTGGCCTAGCAGTATTGCCAGGTCACCCCGGG‑CGCGTTTTCATTTAAAGTTGCACAAGCATCCGTGCGAGGAGCCGCTACATGGGCGCAACCAACAGGTTAGCCTGCTATTCCGATCCCGCTGATCACTATTCTCATCGGGtt > 3:322929/1‑149 (MQ=255) ggTGGCCTAGCAGTATTGCCAGGCCACCCCGGG‑CGCGTTTTCATTTAAAGTTGCACAAGCATCCGTGCGAGGAGCCGCTACATGGGCGCAACCAACAGGTTAGCCTGCTATTCCGATCCCGCTGATCACTATTCTCATCGGGTTCGCCt < 4:184814/149‑1 (MQ=255) gCCTGGCAATACTGCTAGGCCACCCCGGG‑CGCGTTTTCATTTAAAGTTGCACAAGCATCCGTGCGAGGAGCCGCTACATGGGCGCAACCAACAGGTTAGCCTGCTATTCCGATCCCGCTGATCACTAtt > 4:14242/1‑129 (MQ=255) gCCTGGCAATACTGCTAGGCCACCCCGGG‑CGCGTTTTCATTTAAAGTTGCACAAGCATCCGTGCGAGGAGCCGCTACATGGGCGCAACCAACAGGTTAGCCTGCTATTCCGATCCCGCTGATCACTAtt < 3:14242/129‑1 (MQ=255) gCCTAGCAGTATTGCCAGGCCACCCCGGG‑CGCGTTTTCATTTAAAGTTGCACAAGCATCCGTGCGAGGAGCCGCTACATGGGCGCAACCAACAGGTTAGCCTGCTATTCCGATCCCGCTGATCACTATTCTCATCGGGTTCGCCTGGtt > 5:194694/1‑149 (MQ=255) gCCTAGCAGTATTGCCAGGCCACCCCGGG‑CGCGTTTTCATTTAAAGTTGCACAAGCATCCGTGCGAGGAGCCGCTACATGGGCGCAACCAACAGGTTAGCCTGCTATTCCGATCCCGCTGATCACTATTCTCATCGGGTTCGCCTGGt > 8:276062/1‑148 (MQ=255) gCCTAGCAGTATTGCCAGGCCACCCCGGG‑CGCGTTTTCATTTAAAGTTGCACAAGCATCCGTGCGAGGAGCCGCTACATGGGCGCAACCAACAGGTTAGCCTGCTATTCCGATCCCGCTGATCACTATTCTCATCGGGTTCGCCTGGt < 1:345233/148‑1 (MQ=255) ggCAATACTGCTAGGCCACCCCGGGCCGCGTTTTCATTTAAAGTTGCACAAGCATCCGTGGGAGGAGCCGCTACATGGGCGCAACCAACAGGTTAGCCTGCTGTTCCGATCCCGCTGATCACTATTCTCATCGGGTTCGCCTGGTTCTg > 1:185758/2‑149 (MQ=255) gCAGTATTGCCAGGCCACCCCGGG‑CGCGTTTTCATTTAAAGTTGCACAAGCATCCGTGCGAGGAGCCGCTACATGGGCGCAACCAACAGGTTAGCCTGCTATTCCGATCCCGCTGATCACTATTCTCATCGGGTTCGCCTGGTTCTGGc < 8:112979/149‑1 (MQ=255) gCAATACTGCTAGGCCACCCCGGG‑CGCGTTTTCATTTAAAGTTGCACAAGCATCCGTGCGAGGAGCCGCTACATGGGCGCAACCAACAGGTTAGCCTGCTATTCCGATCCCGCTGATCACTATTCTCATCGGGTTCGCCTGGTTCTGGc > 7:275107/1‑149 (MQ=255) aGTATTGCCAGGCCACCCCGGG‑CGCGTTTTCATTTAAAGTTGCACAAGCATCCGTGCGAGGAGCCGCTACATGGGCGCAACCAACAGGTTAGCCTGCTATTCCGATCCCGCTGATCACTATTCTCATCGGGTTCGCCTGGTTCTGGcc > 2:33663/1‑148 (MQ=255) gTATTGCCAGGCCACCCCGGG‑CGCGTTTTCATTTAAAGTTGCACAAGCATCCGTGCGCGGAGCCGCTACATGGGCGCAACCAACAGGTTAGCCTGCTATTCCGATCCCGCTGATCACTATTCTCATCGGGTTCGCCTGGTTCTGGcc > 5:253333/1‑147 (MQ=255) | AAGCCAGCGCATCGGCACCTATGTGTTGCTGTATCTGGCTTTCTTCTTCGTATTCGCCTACTTGCTCAAGCGTGAATACTGGAAGGACGTGCACTGATCACGCAGTAGTTTCTGGTTGTTGAACGCGCCCGGGGTGGCCTAGCAGTATTGCCAGGCCACCCCGGG‑CGCGTTTTCATTTAAAGTTGCACAAGCATCCGTGCGAGGAGCCGCTACATGGGCGCAACCAACAGGTTAGCCTGCTATTCCGATCCCGCTGATCACTATTCTCATCGGGTTCGCCTGGTTCTGGCC > NC_002947/1506783‑1507073 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |