Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 5,840,670:1 | +A | 10.8% | intergenic (‑34/+2) | sseA ← / ← PP_5119 | 3‑mercaptopyruvate sulfurtransferase/TetR family transcriptional regulator |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 5,840,670 | 1 | . | A | 10.8% | 125.5 / 11.5 | 37 | intergenic (‑34/+2) | sseA/PP_5119 | 3‑mercaptopyruvate sulfurtransferase/TetR family transcriptional regulator |
Reads supporting (aligned to +/- strand): ref base . (15/18); new base A (2/2); total (17/20) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.43e-01 |
TGCGCTGGCCATAGTCCACATCCTCGAGGGCAAAGCGGCAGTCGAGGATCACCAACTTGGGCGAGCCCAGACGCTCGGCCAACTGCTGCGGGGTGATCAATTGCGCAAGGGGCATGACAATCTCCTGATCGATGGCACAAGAGCCCTGG‑GCTAGAGCTCCAGGGCCTGGTTGAACGGTACATGGAATTCCTGACACAGCGCTTCTACGGCAGTACGGGCATTGTCGGTGACGAAACCGAGCTCCAGCACCAGCACCTGATAGACGCCGCGCTTGAAGGCTTCTTCACCCAA > NC_002947/5840522‑5840812 | tGCGCTGGCCATAGTCCACATCCTCGAGGGCAAAGCGGCAGTCGAGGATCACCAACTTGGGCGAGCCCAGACGCTCGGCCAACTGCTGCGGGGTGATCAATTGCGCAAGGGGCATGACAATCTCCTGATCGATGGCACAAGAGCCCTgg‑g < 2:845826/150‑1 (MQ=255) tGGCCATAGTCCACATCCTCGAGGGCAAAGCGGCAGTCGAGGATCACCAACTTGGGCGAGCCCAGACGCTCGGCCAACTGCTGCGGGGTGATCAATTGCGCAAGGGGCATGACAATCTCCTGATCGATGGCACAAGAGACCTGG‑GCTaga > 1:177849/1‑150 (MQ=255) ccACATCCTCGAGGGCAAAGCGGCAGTCGAGGATCACCAACTTGGGCGAGCCCAGACGCTCGGCCAACTGCTGCGGGGTGATCAATTGCGCAAGGGGCATGACAATCTCCTGATCGATGGCACAAGAGCCCTGG‑GCTAGAGCTccagggc > 2:456616‑M2/1‑143 (MQ=255) cacaTCCTCGAGGGCAAAGCGGCAGTCGAGGATCACCAACTTGGGCGAGCCCAGACGCTCGGCCAACTGCTGCGGGGTGATCAATTGCGCAAGGGGCATGACAATCTCCTGATCGATGGCACAAGCGCCCTGG‑GCTCGAGCTccggcgcc < 2:607191‑M2/150‑9 (MQ=255) gggCAAAGCGGCAGTCGAGGATCACCAACTTGGGCGAGCCCAGACGCTCGGCCAACTGCTGCGGGGTGATCAATTGCGCAAGGGGCATGACAATCTCCTGATCGATGGCACAAGAGCCCTGG‑GCTAGAGCTCCAGGGCCTgg > 1:443574/1‑142 (MQ=255) gggCAAAGCGGCAGTCGAGGATCACCAACTTGGGCGAGCCCAGACGCTCGGCCAACTGCTGCGGGGTGATCAATTGCGCAAGGGGCATGACAATCTCCTGATCGATGGCACAAGAGCCCTGG‑GCTAGAGCTCCAGGGCCTgg < 2:443574/142‑1 (MQ=255) cAAAGCGGCAGTCGAGGATCACCAACTTGGGCGAGCCCAGACGCTCGGCAAACTGCTGCGGGGTGATCAATTGCGCAAGGGGCATGACAATCTCCTGATCGATGGCACAAGAGCCCTGG‑GCTAGAGCTCCAGGGCCTGGTTGAACGGTAc < 1:612705/150‑1 (MQ=255) aaGCGGCAGTCGAGGATCACCAACTTGGGCGAGCCCAGACGCTCGGCCAACTGCTGCGGGGTGATCAATTGCGCAAGGGGCATGACAATCTCCTGATCGATGGCACAAGAGCCCTGG‑GCTAGAGCTCCAGGGCCTGGTTGAACGGTACAt < 2:332304/150‑1 (MQ=255) aGCGGCAGTCGAGGATCACCAACTTGGGCGAGCCCAGACGCTCGGCCAACTGCTGCGGGGTGATCAATTGCGCAAGGGGCATGACAATCTCCTGATCGATGGCACAAGAGCCCTGG‑GCTAGAGCTCCAGGGCCTGGTTGAACGGTACATg < 2:452423/150‑1 (MQ=255) gAGGATCACCAACTTGGGCGAGCCCAGACGCTCGGCCAACTGCTGCGGGGTGATCAATTGCGCAAGGGGCATGACAATCTCCTGATCGATGGCACAAGAGCCCTGG‑GCTAGAGCTCCAGGGCCTGGTTGAACGGTACATGGAATTCCTGa < 1:745050/150‑1 (MQ=255) aGGATCACCAACTTGGGCGAGCCCAGACGCTCGGCCAACTGCTGCGGGGTGATCAATTGCGCAAGGGGCATGACAATCTCCTGATCGATGGCACAAG‑GCCCTGGAGCTCTAGCcagggctctt > 1:637492‑M2/1‑114 (MQ=255) aGGATCACCAACTTGGGCGAGCCCAGACGCTCGGCCAACTGCTGCGGGGTGATCAATTGCGCAAGGGGCATGACAATCTCCTGATCGATGGCACAAG‑GCCCTGGAGCTCTAGCcagggctctt < 2:637492‑M2/124‑11 (MQ=255) ccAACTTGGGCGAGCCCAGACGCTCGGCCAACTGCTGCGGGGTGATCAATTGCGCAAGGGGCATGACAATCTCCTGATCGATGGCACAAGAGCCCTGGAGCTCTAGCcagggctcttgtgccatcgatcaggagattg > 2:77326‑M2/1‑108 (MQ=255) ccAACTTGGGCGAGCCCAGACGCTCGGCCAACTGCTGCGGGGTGATCAATTGCGCAAGGGGCATGACAATCTCCTGATCGATGGCACAAGAGCCCTGGAGCTCTAGCcagggctcttgtgccatcgatcaggagattg < 1:77326‑M2/139‑32 (MQ=255) cGAGCCCAGACGCTCGGCCAACTGCTGCGGGGTGATCAATTGCGCAAGGGGCATGACAATCTCCTGATCGATGGCACAAGAGCCCTGG‑GCTAGAGCTCCAGGGCCTGGTTGAACGTTACATGGAATTCCTGACACAGCTCTTCTACGGCa < 2:177849/150‑1 (MQ=255) gCCCAGACGCTCGGCCAACTGCTGCGGGGTGATCAATTGCGCAAGGGGCATGACAATCTCCTGATCGATGGCACAAGAGCCCTGG‑GCTAGAGCTCCAGGGCCTGGTTGAACGGTACATGGAATTCCTGACACAGCGCTTCTACGGCAGTa > 2:676667/1‑150 (MQ=255) gACGCTCGGCCAACTGCTGCGGGGTGATCAATTGCGCAAGGGGCATGACAATCTCCTGATCGATGGCACAAGAGCCCTGG‑GCTAGAGCTCCAGGGCCTGGTTGAACGGTACATGGAATTCCTGACACAGCGCTTCTACGGCAGTACGGGc > 1:576268/1‑150 (MQ=255) aaCTGCTGCGGGGTGATCAATTGCGCAAGGGGCATGACAATCTCCTGATCGATGGCACAAGAGCCCTGGagctctagcccagggctc < 2:513019‑M1/87‑19 (MQ=255) aaCTGCTGCGGGGTGATCAATTGCGCAAGGGGCATGACAATCTCCTGATCGATGGCACAAGAGCCCTGGagctctagcccagggctc > 1:513019‑M1/1‑69 (MQ=255) aTTGCGCAAGGGGCATGACAATCTCCTGATCGATGGCACAAGAGCCCTGG‑GCTAGAGCTCCAGGGCCTGGTTGAACGGTACATGGAATTCCTGACACAGCGCTTCTACGGCAGTACGGGCATTGTCGGTGACGAAACCGAGCTccagcac < 1:208420/150‑1 (MQ=255) aTTGCGCAAGGGGCATGACAATCTCCTGATCGATGGCACAAGAGCCCTGG‑GCTAGAGCTCCAGGGCCTGGTTGAACGGTACATGGAATTCCTGACACAGCGATTTTACGGCAGTACGGGCATTGTCGGTGACGAAACCGAGCTccagcac > 1:520386/1‑150 (MQ=255) ttGCGCAAGGGGCATGACAATCTCCTGATCGATGGCACAAGAGCCCTGG‑GCTAGAGCTCCAGGGCCTGGTTGAACGGTACATGGAATTCCTGACACAGCGCTTCTACGGCAGTACGGGCATTGTCGGTGACGAAACCGAGCTccagcacc > 2:578924/1‑150 (MQ=255) gcAAGGGGCATGACAATCTCCTGATCGATGGCACAAGAGCCCTGG‑GCTAGAGCTCCAGGGCCTGGTTGAACGGTACATGGAATTCCTGACACAGCGCTTCTACGGCAGTACGGGCATTGTCGGTGACGAAACCGAGCTccagcaccagca > 2:647785/1‑150 (MQ=255) ggggCATGACAATCTCCTGATCGATGGCACAAGAGCCCTGG‑GCTAGAGCTCCAGGGCCTGGTTGAACGGTACATGGAATTCCTGACACAGCGCTTCTACGGCAGTACGGGCATTGTCGGTGa > 2:684183/1‑122 (MQ=255) ggggCATGACAATCTCCTGATCGATGGCACAAGAGCCCTGG‑GCTAGAGCTCCAGGGCCTGGTTGAACGGTACATGGAATTCCTGACACAGCGCTTCTACGGCAGTACGGGCATTGTCGGTGa < 1:684183/122‑1 (MQ=255) aTGACAATCTCCTGATCGATGGCACAAGAGCCCTGG‑GCTAGAGCTccagg < 2:713140‑M2/50‑6 (MQ=255) aTGACAATCTCCTGATCGATGGCACAAGAGCCCTGG‑GCTAGAGCTccagg > 1:713140‑M2/1‑45 (MQ=255) gACAATCTCCTGATCGATGGCACAAGAGCCCTGG‑GCTAGAGCTCCAGGGCCTGGTTGAACGGTACATGGAATTCCTGAGTGGCTTCTTCTACGGCAGTACGGGCATTGTCGGTGACGAAACCGAGCTCCAGGACCAGCACCTGATAGAcg > 1:53464/1‑150 (MQ=255) cAATCTCCTGATCGATGGCACAAGAGCCCTGGgctagagctccagggctcttgtgccatcgatcaggagattgtcatgccccttgcgcaattgatcaccccgcagcagttggccgagcgtctgggctcgcccaagttgg < 1:77326‑M1/139‑108 (MQ=255) cAATCTCCTGATCGATGGCACAAGAGCCCTGGgctagagctccagggctcttgtgccatcgatcaggagattgtcatgccccttgcgcaattgatcaccccgcagcagttggccgagcgtctgggctcgcccaagttgg > 2:77326‑M1/1‑32 (MQ=255) aTCGATGGCACAAGAGCCCTGG‑GCTAGAGCTCCAGGGCCTGGTTGAACGGTACATGGAATTCCTGACACAGCGCTTCTACGGCAGTACgg < 1:843384/90‑1 (MQ=255) aTCGATGGCACAAGAGCCCTGG‑GCTAGAGCTCCAGGGCCTGGTTGAACGGTACATGGAATTCCTGACACAGCGCTTCTACGGCAGTACgg > 2:843384/1‑90 (MQ=255) aTGGCACAAGAGCCAAGG‑GCTAGAGGGCCAGGGCCTGGTTGAACGGTACATGGAATTCCTGACACACCCATCCGACGGCAGTACGGGCATTGTCGGTGACGAAACCGAGCTCCAGCACCAGCACCTGATAGACGCCGCGCTTGAAGGctt < 1:15419/150‑1 (MQ=255) acaAGAGCCCTGG‑GCTAGAGCTCCAGGGCCTGGTTGAACGGTACATGGAATTCCTGACACAGCGCTTCTACGGCAGTACGGGCATTGTCGGTGACGAAACCGAGCTCCAGCACCAGCACCTGATAGACGCCGCGCTTGAAGGCTTCTTCa > 2:303348/1‑150 (MQ=255) acaAGAGCCCTGG‑GCTAGAGCTCCAGGGCCTGGTTGAACGGTACATGGAATTCCTGACACAGCGCTTCTACGGCAGTACGGGCATTGTCGGTGACGAAACCGAGCTCCAGCACCAGCACCTGATAGACGCCGCGCTTGAAGGCTTCTTCa < 1:676667/150‑1 (MQ=255) aaGAGCCCTGGgctagagctccagggccttgtgccatcgatcaggagattgtcatgccccttgcgcaattgatcaccccgcagcagttggccgagcgtctgggctcgcccaagttggtgatcct < 2:637492‑M1/124‑114 (MQ=255) aaGAGCCCTGGgctagagctccagggccttgtgccatcgatcaggagattgtcatgccccttgcgcaattgatcaccccgcagcagttggccgagcgtctgggctcgcccaagttggtgatcct > 1:637492‑M1/1‑11 (MQ=255) aaGAGCCCTGG‑GCTAGAGCTCCAGGGCCTGGTTGAACGGTACATGGAATTCCTGACACAGCGCTTCTACGGCAGTACGGGCATTGTCGGTGACGAAACCGAGCTCCAGCACCAGCACCTGATAGACGCCGCGCTTGAAGGCTTCTTCAcc > 1:32136/1‑150 (MQ=255) gagCCCTGG‑GCTAGAGCTccagggctcttgtgccatcgatcaggagattgtcatgccccttgcgcaattgatcaccccgcagcagtt < 2:513019‑M2/87‑70 (MQ=255) gagCCCTGG‑GCTAGAGCTccagggctcttgtgccatcgatcaggagattgtcatgccccttgcgcaattgatcaccccgcagcagtt > 1:513019‑M2/1‑18 (MQ=255) ggcgccggagctcgagcccagggcgcttgtgccatcgatcaggagattgtcatgccccttgcgcaattgatcaccccgcagcagttggccgagcgtctgggctcgcccaagttggtgatcctcgactgccgctttgccctcgaggatgtg < 2:607191‑M1/144‑143 (MQ=255) gCCCTGGagctctagcccagggctcttgtgccatcgatcaggagattgtcatgccccttgcgcaattgatcaccccgcagcagttggccgagcgtctgggctcgcccaagttggtgatcctcgactgccgctttgccctcgaggatgtgg > 2:456616‑M1/1‑7 (MQ=255) ccTGGagctctagcccagggctcttgtgccatcgatcaggagattgtcat < 2:713140‑M1/50‑46 (MQ=255) ccTGGagctctagcccagggctcttgtgccatcgatcaggagattgtcat > 1:713140‑M1/1‑5 (MQ=255) | TGCGCTGGCCATAGTCCACATCCTCGAGGGCAAAGCGGCAGTCGAGGATCACCAACTTGGGCGAGCCCAGACGCTCGGCCAACTGCTGCGGGGTGATCAATTGCGCAAGGGGCATGACAATCTCCTGATCGATGGCACAAGAGCCCTGG‑GCTAGAGCTCCAGGGCCTGGTTGAACGGTACATGGAATTCCTGACACAGCGCTTCTACGGCAGTACGGGCATTGTCGGTGACGAAACCGAGCTCCAGCACCAGCACCTGATAGACGCCGCGCTTGAAGGCTTCTTCACCCAA > NC_002947/5840522‑5840812 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |