Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 1,506,939 | T→C | 23.4% | intergenic (+51/‑66) | 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,939 | 0 | T | C | 23.4% | 92.8 / 22.9 | 51 | intergenic (+51/‑66) | petC/sspA | ubiquinol‑‑cytochrome c reductase cytochrome c1 subunit/stringent starvation protein A |
Reads supporting (aligned to +/- strand): ref base T (21/18); new base C (5/7); total (26/25) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 5.23e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01 |
CAGCGCATCGGCACCTATGTGTTGCTGTATCTGGCTTTCTTCTTCGTATTCGCCTACTTGCTCAAGCGTGAATACTGGAAGGACGTGCACTGATCACGCAGTAGTTTCTGGTTGTTGAACGCGCCCGGGGTGGCCTAGCAGTATTGCCAGGCCACCCCGGGCGCGTTTTCATTTAAAGTTGCACAAGCATCCGTGCGAGGAGCCGCTACATGGGCGCAACCAACAGGTTAGCCTGCTATTCCGATCCCGCTGATCACTATTCTCATCGGGTTCGCCTGGTTCTGGCCGAGAAG > NC_002947/1506796‑1507088 | cAGCGCATCGGCACCTATGTGTTGCTGTATCTGGCTTTCTTCTTCGTATTCGCCTACTTGCTCAAGCGTGAATACTGGAAGGACGTGCACTGATCACGCAGTAGTTTCTGGTTGTTGAACGCGCCCGGGGTGGCCTAGCAGTATTGCCAg < 1:1410816/150‑1 (MQ=255) aGCGCATCGGCACCTATGTGTTGCTGTATCTGGCTTTCTTCTTCGTATTCGCCTACTTGCTCAAGCGTGAATACTGGAAGGACGTGCACTGATCACGCAGTAGTTTCTGGTTGTTGAACGCGCCCGGGGTGGCCTAGCAGTATTGCCAgg < 1:2019878/150‑1 (MQ=255) gcgcATCGGCACCTATGTGTTGCTGTATCTGGCTTTCTTCTTCGTATTCGCCTACTTGCTCAAGCGTGAATACTGGAAGGACGTGCACTGATCACGCAGTAGTTTCTGGTTGTTGAACGCGCCCGGGGTGGCCTGGCAATACTGCTAGGc < 2:94238/150‑1 (MQ=255) gcgcATCGGCACCTATGTGTTGCTGTATCTGGCTTTCTTCTTCGTATTCGCCTACTTGCTCAAGCGTGAATACTGGAAGGACGTGCACTGATCACGCAGTAGTTTCTGGTTGTTGAACGCGCCCGGGGTGGCCTAGCAGTATTGCCAGGc < 1:1881076/150‑1 (MQ=255) cgcATCGGCACCTATGTGTTGCTGTATCTGGCTTTCTTCTTCGTATTCGCCTACTTGCTCAAGCGTGAATACTGGAAGGACGTGCACTGATCACGCAGTAGTTTCTGGTTGTTGAACGCGCCCGGGGTGGCCTAGCAGTATTGCCAGGcc < 2:2059002/150‑1 (MQ=255) aTCGGCACCTATGTGTTGCTGTATCTGGCTTTCTTCTTCGTATTCGCCTACTTGCTCAAGCGTGAATACTGGAAGGACGTGCACTGATCACGCAGTAGTTTCTGGTTGTTGAACGCGCCCGGGGTGGCCTGGCAATACTGCTAGGCcccc < 1:2333303/150‑1 (MQ=39) aCCTATGTGTTGCTGTATCTGGCTTTCTTCTTCGTATTCGCCTACTTGCTCAAGCGTGAATACTGGAAGGACGTGCACTGATCACGCAGTAGTTTCTGGTTGTTGAACGCGCCCGGGGTGGCCTAGCAGTATTGCCAGGCCACCCCGGgc < 1:984209/150‑1 (MQ=255) tgtTGCTGTATCTGGCTTTCTTCTTCGTATTCGCCTACTTGCTCAAGCGTGAATACTGGAAGGACGTGCACTGATCACGCAGTAGTTTCTGGTTGTTGAACGCGCCCGGGGTGGCCTAGCAGTATTGCCAGGCCACCCCg < 1:1996537/140‑1 (MQ=255) tgtTGCTGTATCTGGCTTTCTTCTTCGTATTCGCCTACTTGCTCAAGCGTGAATACTGGAAGGACGTGCACTGATCACGCAGTAGTTTCTGGTTGTTGAACGCGCCCGGGGTGGCCTAGCAGTATTGCCAGGCCACCCCg > 2:1996537/1‑140 (MQ=255) tGCTGTATCTGGCTTTCTTCTTCGTATTCGCCTACTTGCTCAAGCGTGAATACTGGAAGGACGTGCACTGATCACGCAGTAGTTTCTGGTTGTTGAACGCGCCCGGGGTGGCCTAGCAGTATTGCCAGGCCACCCCGGGCGCGTTTTCAt < 1:928123/150‑1 (MQ=255) tGGCTTTCTTCTTCGTATTCGCCTACTTGCTCAAGCGTGAATACTGGAAGGACGTGCACTGATCACGCAGTAGTTTCTGGTTGTTGAACGCGCCCGGGGTGGCCTAGCAGTAtt > 1:2021833/1‑114 (MQ=255) tGGCTTTCTTCTTCGTATTCGCCTACTTGCTCAAGCGTGAATACTGGAAGGACGTGCACTGATCACGCAGTAGTTTCTGGTTGTTGAACGCGCCCGGGGTGGCCTAGCAGTAtt < 2:2021833/114‑1 (MQ=255) tcttcGTATTCGCCTACTTGCTCAAGCGTGAATACTGGAAGGACGTGCACTGATCACGCAGTAGTTTCTGGTTGTTGAACGCGCCCGGGGTGGCCTGGCAATACTGc < 1:1547892/107‑1 (MQ=39) ttcGTATTCGCCTACTTGCTCAAGCATGAATACTGGAAGGACGTGCACTGATCACGCAGTAGTTTCTGGTTGTTGAACGCGCCCGGGGTGGCCTAGCAGTATTGCCAg < 2:1963835/108‑1 (MQ=255) ttcGTATTCGCCTACTTGCTCAAGCATGAATACTGGAAGGACGTGCACTGATCACGCAGTAGTTTCTGGTTGTTGAACGCGCCCGGGGTGGCCTAGCAGTATTGCCAg > 1:1963835/1‑108 (MQ=255) aTTCGCCTACTTGCTCAAGCGTGAATACTGGAAGGACGTGCACTGATCACGCAGTAGTTTCTGGTTGTTGAACGCGCCCGGGGTGGCCTGGCAATACTGCTAGGc > 2:1221975/1‑105 (MQ=37) aTTCGCCTACTTGCTCAAGCGTGAATACTGGAAGGACGTGCACTGATCACGCAGTAGTTTCTGGTTGTTGAACGCGCCCGGGGTGGCCTGGCAATACTGCTAGGc < 1:1221975/105‑1 (MQ=38) ttCGCCTACTTGCTCAAGCGTGAATACTGGAAGGACGTGCACTGATCACGCAGTAGTTTCTGGTTGTTGAACGCGCCCGGGGTGGCCTGGCAATACTGCTAGGc < 1:1848931/104‑1 (MQ=37) ttCGCCTACTTGCTCAAGCGTGAATACTGGAAGGACGTGCACTGATCACGCAGTAGTTTCTGGTTGTTGAACGCGCCCGGGGTGGCCTGGCAATACTGCTAGGc > 2:1848931/1‑104 (MQ=37) cTTGCTCAAGCGTGAATACTGGAAGGACGTGCACTGATCACGCAGTAGTTTCTGGTTGTTGAACGCGCCCGGGGTGGCCTAGCAGTAtt < 1:1403469/89‑1 (MQ=255) cTTGCTCAAGCGTGAATACTGGAAGGACGTGCACTGATCACGCAGTAGTTTCTGGTTGTTGAACGCGCCCGGGGTGGCCTAGCAGTAtt > 2:1403469/1‑89 (MQ=255) tCAAGCGTGAATACTGGAAGGACGTGCACTGATCACGCAGTAGTTTCTGGTTGTTGAACGCGCCCGGGGTGGCCTAGCAGTATTGCCAGGCCACCCCGGGCGCGTTTTCATTTAAAGTTGCACAAGCATCCGTGCGAGGAGCCGCTACAt > 2:572798/1‑150 (MQ=255) aaGCGTGAATACTGGAAGGACGTGCACTGATCACGCCGTAGCTTCTGGTTGTTGAACGCGCCCGGGGTGGCCAGGCAGTATTGACAGGCCACGCCGGGCGCGTTTTCATTTAAAGTTGCACAAGCATCCGTGCGAGGAGCCGCTACATgg > 1:1798581/1‑150 (MQ=255) aaGCGTGAATACTGGAAGGACGTGCACTGATCACGCAGTAGTTTCTGGTTGTTGAACGCGCCCGGGGTGGCCTGGCAATACTGCTAGGc > 2:722233/1‑89 (MQ=25) aaGCGTGAATACTGGAAGGACGTGCACTGATCACGCAGTAGTTTCTGGTTGTTGAACGCGCCCGGGGTGGCCTGGCAATACTGCTAGGc < 1:722233/89‑1 (MQ=25) gAATACTGGAAGGACGTGCACTGATCACGCAGTAGTTTCTGGTTGTTGAACGCGCCCGGGGTGGCCTAGCAGTATTGCCAGGCCACCCCGGGCGCGTTTTCATTTAAAGTTGCACAAGCATCCGTGCGAGGAGCCGCTACATGGGCGcaa > 1:1363534/1‑150 (MQ=255) aaTACTGGAAGGACGTGCACTGATCACGCAGTAGTTTCTGGTTGTTGAACGCGCCCGGGGTGGCCTAGCAGTATTGCCAGGCCACCCCGGGCGCGTTTTCATTTAAAGTTGCACAAGCATCCGTGCGAGGAGCCGCTACATGGGCGcaa < 1:2063388/149‑1 (MQ=255) aaTACTGGAAGGACGTGCACTGATCACGCAGTAGTTTCTGGTTGTTGAACGCGCCCGGGGTGGCCTAGCAGTATTGCCAGGCCACCCCGGGCGCGTTTTCATTTAAAGTTGCACAAGCATCCGTGCGAGGAGCCGCTACATGGGCGcaa > 2:2063388/1‑149 (MQ=255) cTGGAAGGACGTGCACTGATCACGCAGTAGTTTCTGGTTGTTGAACGCGCCCGGGGTGGCCTAGCAGTATTGCCAGGCCACCCCGGGCGCGTTTTCATTTAAAGTTGCACAAGCATCCGTGCGAGGAGCCGCTACATGGGCGcaaccaac > 2:928559/1‑150 (MQ=255) tGGAAGGACGTGCACTGATCACGCAGTAGTTTCTGGTTGTTGAACGCGCCCGGGGTGGCCTAGCAGTAtt < 1:636862/70‑1 (MQ=255) tGGAAGGACGTGCACTGATCACGCAGTAGTTTCTGGTTGTTGAACGCGCCCGGGGTGGCCTAGCAGTAtt > 2:636862/1‑70 (MQ=255) gAAGGACGTGCACTGATCACGCAGTAGTTTCTGGTTGTTGAACGCGCCCGGGGTGGCCTAGCAGTATTGCCAGGCCACCCCGGGCGCGTTTTCATTTAAAGTTGCACAAGCATCCGTGCGAGGAGCCGCTACATGGGCGcaaccaa < 1:1321095/146‑1 (MQ=255) gAAGGACGTGCACTGATCACGCAGTAGTTTCTGGTTGTTGAACGCGCCCGGGGTGGCCTAGCAGTATTGCCAGGCCACCCCGGGCGCGTTTTCATTTAAAGTTGCACAAGCATCCGTGCGAGGAGCCGCTACATGGGCGcaaccaa > 2:1321095/1‑146 (MQ=255) aGGACGTGCACTGATCACGCAGTAGTTTCTGGTTGTTGAACGCGCCCGGGGTGGCCTAGCAGTATTGCCAGGCCACCCCGGGCGCGTTTTCATTTAAAGTTGCACAAGCATCCGTGCGAGGAGCCGCTACATGGGCGcaaccaa > 1:1177704/1‑144 (MQ=255) aGGACGTGCACTGATCACGCAGTAGTTTCTGGTTGTTGAACGCGCCCGGGGTGGCCTAGCAGTATTGCCAGGCCACCCCGGGCGCGTTTTCATTTAAAGTTGCACAAGCATCCGTGCGAGGAGCCGCTACATGGGCGcaaccaa < 2:1177704/144‑1 (MQ=255) cGTGCACTGATCACGCAGTAGTTTCTGGTTGTTGAACGCGCCCGGGGTGGCCTAGCAGTATTGCCAGGCCACCCCGGGCGCGTTTTCATTTAAAGTTGCACAAGCATCCGTGCGAGGAGCCGCTACATGGGCGCAACCAACAGGTTAGcc < 1:572798/150‑1 (MQ=255) tGCACTGATCACGCAGTAGTTTCTGGTTGTTGAACGCGCCCGGGGTGGCCTAGCAGTATTGCCAGGCCACCCCGGGCGCGTTTTCATTTAAAGTTGCACAAGCATCCGTGCGAGGAGCCGCTACATGGGCGCAACCAACAGGTTAGCCTg > 2:1577126/1‑150 (MQ=255) aCTGATCACGCAGTAGTTTCTGGTTGTTGAACGCGCCCGTGGGGGCCTAGCAGTATTGCCGGGCCACCCCGGGCGCGTTTTCATTTAAAGTTGCACAAGCATCCGTGCGAGGAGCCGCTACATGGGCGCAACCAACAGGTGCGCCTGCTa > 2:918475/1‑150 (MQ=255) gATCACGCAGTAGTTTCTGGTTGTTGAACGCGCCCGGGGTGGCCTAGCAGTATTGCCAGGCCACCCCGGGCGCGTTTTCATTTAAAGTTGCACAAGCATCCGTGCGAGGAGCCGCTACATGGGCGCAACCAACAGGTTAGCCTGCTATTc < 1:1577126/150‑1 (MQ=255) gtagtTTCTGGTTGTTGAACGCGCCCGGGGTGGCCTAGCAGTATTGCCAGGCCACCCCGGGCGCGTTTTCATTTAAAGTTGCACAAGCATCCGTGCGAGGAGCCGCTACATGGGCGCAACCAACAGGTTAGCCTGCTATTCCGATCCCGc > 1:812066/1‑150 (MQ=255) tagtTTCTGGTTGTTGAACGCGCCCGGGGTGGCCTAGCAGTATTGCCAGGCCACCCCGGGCGCGTTTTCATTTAAAGTTGCACAAGCATCCGTGCGAGGAGCCGCTACATGGGCGCAACCAACAGGTTAGCCTGCTATTCCGATCCCGCt < 1:928559/150‑1 (MQ=255) tttCTGGTTGTTGAACGCGCCCGGGGTGGCCTAGCAGTATTGCCAg < 1:168967/46‑1 (MQ=255) tttCTGGTTGTTGAACGCGCCCGGGGTGGCCTAGCAGTATTGCCAg > 2:168967/1‑46 (MQ=255) ggTTGTTGAACGCGCCCGGGGTGGCCTAGCAGTATTGCCAGGCCACCCCGGGCGCGTTTTCATTTAAAGTTGCACGAGCATCCGTGCGAGGAGCCGCTACATGGGCGCAACCAACAGGTTAGCCTGCTATTCCGATCCCGCTGATCACTa < 1:1242580/150‑1 (MQ=255) gttgttGAACGCGCCCGGGGTGGCCTAGCAGTATTGCCAGGCAACCCCGGGCGCGTTTTCATTTAAAGTTGCACAAGCATCCGTGCGAGGAGCCGCTACATGGGCGCAACCAACAGGTTAGCCTGCTATTCCGATCCCGCTGATCACTAt > 2:692697/1‑150 (MQ=255) gttgAACGCGCCCGGGGTGGCCTAGCAGTATTGCCAGGCCACCCCGGgc < 1:1357740/49‑1 (MQ=255) gttgAACGCGCCCGGGGTGGCCTAGCAGTATTGCCAGGCCACCCCGGgc > 2:1357740/1‑49 (MQ=255) gcCCGGGGGGGCCTAGCAGTATTGCCAGGCCACCCCGGGCGCGTTTTCATTTAAAGTTGCACAAGCATCCGTGCGAGGAGCCGCTACATGGGCGCAACCAACAGGTTAGCCTGCTATTCCGATCCCGCTGATCACTATTCTCATCGGGtt > 2:228480/1‑150 (MQ=255) tGGCCTAGCAGTATTGCCAGGCCACCCCGGGCGCGTTTTCATTTAAAGTTGCACAAGCATCCGTGCGAGGAGCCGCTACATGGGCGCAACCAACAGGTTAGCCTGCTATTc < 1:1655569/111‑1 (MQ=255) tGGCCTAGCAGTATTGCCAGGCCACCCCGGGCGCGTTTTCATTTAAAGTTGCACAAGCATCCGTGCGAGGAGCCGCTACATGGGCGCAACCAACAGGTTAGCCTGCTATTc > 2:1655569/1‑111 (MQ=255) gCCTGGCAATACTGCTAGGCCACCCGGGGCGCGTTTTCATTTAAAGTTGCACAAGCATCCGTGCGAGGAGCCGCTACATGGGCGCAACCAACAGGTTAGCCTGCTATTCCGATCCCGCTGATCACTATTCTCATCGGGTTCGCCt < 2:136413/145‑1 (MQ=255) gCCTGGCAATACTGCTAGGCCACCCCGGGCGCGTTTTCATTTAAAGTTGCACAAGCATCCGTGCGAGGAGCCGCTACATGGGCGCAACCAACAGGTTAGCCTGCTATTCCGATCCCGCTGATCACTATTCTCATCGGGTTCGCCTGGTTc > 1:258087/1‑150 (MQ=255) ccTGGCAATACTGCTAGGCCACCCCGGGAGCGTTTTCATTTAAAGTTGCACAAGCATCCGTGCGACGAGCCGCTACATGGGCGCAACCAACAGGTTAGCCTGCTATTCCGATCCCGCTGATCACTATTCTCATCGGGTTCGCCTGGTTCt > 1:352972/1‑150 (MQ=255) cTAGCAGTATTGCCAGGCCACCCCGGGCGCGTTTTCATTTAAAGTTGCACAAGCATCCGTGCGAGGAGCCGCTACATGGGCGCAACCAACAGGTTAGCCTGCTATTCCGATCCCGCTGATCACTATTCTCATCGGGTTCGCCTGGTTCTg > 2:1573022/1‑150 (MQ=255) cTAGCAGTATTGCCAGGCCACCCCGGGCGCGTTTTCATTTAAAGTTGCACAAGCATCCGTGCGAGGAGCCGCTACATGGGCGCAACCAACAGGTTAGCCTGCTATTCCGATCCCGCTGATCACTATTCTCATCGGGTTCGCCTGGTTCTg > 2:697525/1‑150 (MQ=255) aGCAGTATTGCCAGGCCACCCCGGGCGCGTTTTCATTTAAAGTTGCACAAGCATCCGTGCGAGGAGCCGCTACATGGGCGCAACCAACAGGTTAGCCTGCTATTCCGATCCCGCTGATCACTATTCTCATCGGGTTCGCCTGGTTCTGGc > 2:366761/1‑150 (MQ=255) aGCAGTATTGCCAGGCCACCCCGGGCGCGTTTTCATTTAAAGTTGCACAAGCATCCGTGCGAGGAGCCGCTACATGGGCGCAACCAACAGGTTAGCCTGCTATTCCAATCCCGCTGATCACTATTCTCATCGGGTTCGCCTGGTTCTGGc > 2:2271508/1‑150 (MQ=255) ttGCCAGGCCACCCCGGGCGCGTTTTCATTTAAAGTTGCACAAGCATCCGTGCGAGGAGCCGCTACATGGGCGCAACCAACAGGTTAGCCTGCTATTCCGATCCCGCTGATCACTATTCTCATCGGGTTCGCCTGGTTCTGGCCGAGAAg > 1:600442/1‑150 (MQ=255) | CAGCGCATCGGCACCTATGTGTTGCTGTATCTGGCTTTCTTCTTCGTATTCGCCTACTTGCTCAAGCGTGAATACTGGAAGGACGTGCACTGATCACGCAGTAGTTTCTGGTTGTTGAACGCGCCCGGGGTGGCCTAGCAGTATTGCCAGGCCACCCCGGGCGCGTTTTCATTTAAAGTTGCACAAGCATCCGTGCGAGGAGCCGCTACATGGGCGCAACCAACAGGTTAGCCTGCTATTCCGATCCCGCTGATCACTATTCTCATCGGGTTCGCCTGGTTCTGGCCGAGAAG > NC_002947/1506796‑1507088 |
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: ‑ |