Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 5,537,058 | A→G | 11.8% | intergenic (+35/+68) | nadE → / ← PP_4870 | NH(3)‑dependent NAD(+) synthetase/azurin |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 5,537,058 | 0 | A | G | 11.8% | 83.3 / 3.8 | 34 | intergenic (+35/+68) | nadE/PP_4870 | NH(3)‑dependent NAD(+) synthetase/azurin |
Reads supporting (aligned to +/- strand): ref base A (11/19); new base G (2/2); total (13/21) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 6.27e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.05e-01 |
CAGATCGATGCCTTCCTGCATGGGCAGCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTCACCACCGGC > NC_002947/5536913‑5537203 | cAGATCGATGCCTTCCTGCATGGGCAGCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGa < 8:150526/149‑1 (MQ=255) aGATCGATGCCTTCCTGCATGGGCAGCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATGCgtc > 3:174699/1‑147 (MQ=255) gATCGATGCCTTCCTGCATGGGCAGCCGGTTGACCAGGCAGCGTTTGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAgg < 5:46342/149‑1 (MQ=255) gATCGATGCCTTCCTGCATGGGCAGCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAgg < 4:33295/149‑1 (MQ=255) aTCGTTGGCTTCCTGCATTGGCTGCCGCTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACGCAGCACTAGCTCGAACTGCCGTTGGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCTCCCGCGAATACGAgg < 8:302447/148‑1 (MQ=255) aTGCCTTCCTGCATGGGCAGCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGaa < 8:224503/149‑1 (MQ=255) ccTGCATGGGCAGCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCAc < 4:282604/147‑1 (MQ=255) aGCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCACCGACGCATTcgcg < 6:74151/149‑1 (MQ=255) gCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATGCGTCGGTGGGTTCACCCTCGTATTCGCgg < 4:174699/149‑1 (MQ=25) gCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCACCGACGCATTCGCgg < 6:129650/149‑1 (MQ=255) gTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAgg < 1:152438/121‑1 (MQ=255) gTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAgg > 2:152438/1‑121 (MQ=255) gACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATGCGTCGGTGGGTTCACCCTCGTATTCGCGGGTGAAc < 5:251839/148‑1 (MQ=21) aCCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAAcc < 2:38990/148‑1 (MQ=255) aCCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCACCGACGCATTAGCGGGTGAAccc > 1:114183/1‑149 (MQ=255) ccACCTACCCAACTACCCACCACACGCGCTACATCCCTTTCTCTCCTTTAAAAGTGGCCTCCAGCCGCGGGTTCACCCGCTAATACTAGGGTTAACCAACCCACGATTTCGCGGGTCAACCCGCTCCTGCAGGGTTCGCCGTCAGCCtt < 6:265069/149‑1 (MQ=255) ccccacgACCCCGCACACGCGCGAACTGCCTTTCGCTCTTTTATAAGGGGCCTCCCGGCGCGGGTTCACCCGCTAATAGTCGGGTCAACCCACCGACGCATTCGCGGTTTAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTAc < 2:114183/143‑1 (MQ=255) cGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTACAGGGTTCGCCGTCa < 1:108040/136‑1 (MQ=255) cGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTAACCCGCGAATACGAGGGTGAAACCACCCACGCATTAGCGGGTGAACCCGCTCCTACAGGGGTCGcccgca > 2:108040/1‑132 (MQ=37) gcaacccccAGCACAAGCGCCAACTGCCGTTCGGCCCTTGATAAGTTGGCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCACCGACGCATTCGGGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACtt < 4:226471/143‑1 (MQ=38) gCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTAACCCGCCAAAACAAAGGTGAACCCACCCACGCATTAGCCGGTGAACCCCCACCAGCAGGGTTAGCCGTAAGCCATTGATTaaa > 2:64084/1‑146 (MQ=255) cccAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTAACCCGCGAATACGAGGGTGAACCCCCCGACGCATTCGCGGGTGAACCCGCTACTGCAGGGTTAGCCGGCAGCCTTTGATTAATTGACGAt > 6:88631/1‑149 (MQ=38) ccAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGGATACGAGGGTGGACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCCGGGGTCGCCGTCAGCCGTTGATTTCTTGACGCTg > 1:13951/1‑149 (MQ=38) acaAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCtt > 1:62734/1‑149 (MQ=255) ccaaCTGCCGTTCGCTCATTTATAAGTGCCCTGCTGGAGCGGGTTCACCCGCTAATACGAGGGTCAACCCACCGACGCATTCGCGGGTTAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTtc < 7:308759/147‑1 (MQ=255) gAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTGACCCGCGAATAAGAGGGTGAACCCACCCACGCATTAGCGGGTGAACCCGCACCGGGAGGGTTCGCCGGCAGCCGTTGATTTCTTGGCGGTGACCGTGCCTTtc > 1:220156/1‑148 (MQ=255) gAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTtca < 7:273654/149‑1 (MQ=255) aCTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTtcatc < 1:64084/149‑1 (MQ=255) cccGCGAATGCGTCGGTGGGTTCACCCTCGTATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTcacc > 8:99823/1‑149 (MQ=25) ccGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCagaa > 4:265853/1‑133 (MQ=255) ccGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCagaa < 3:265853/133‑1 (MQ=255) ccGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTcacca > 5:150440/1‑149 (MQ=255) gcgAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGAt > 8:294469/1‑86 (MQ=255) aaTACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTCACCACCGGc < 2:13951/149‑1 (MQ=255) | CAGATCGATGCCTTCCTGCATGGGCAGCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTCACCACCGGC > NC_002947/5536913‑5537203 |
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: ‑ |