Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | CP009273 | 3,371,135 | G→T | 29.3% | intergenic (‑356/+39) | sspA ← / ← rpsI | stringent starvation protein A, phage P1 late gene activator, RNAP‑associated acid‑resistance protein, inactive glutathione S‑transferase homolog/30S ribosomal subunit protein S9 |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | CP009273 | 3,371,135 | 0 | G | T | 29.3% | 17.2 / 6.2 | 17 | intergenic (‑356/+39) | sspA/rpsI | stringent starvation protein A, phage P1 late gene activator, RNAP‑associated acid‑resistance protein, inactive glutathione S‑transferase homolog/30S ribosomal subunit protein S9 |
Reads supporting (aligned to +/- strand): ref base G (4/8); new base T (3/2); total (7/10) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 5.93e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01 |
GCTTTAGTCTGAAATAGTTGACTTAGTCCCTTATTGGCGATGTGGTTTTTGTTTTACCTGTCTGTCAGGTGGCAGCAAAAAGCAACTTTCCAGTTTTTACGCTGATTCAGATTTTAGCTATAAAAAAACCCGCCGAAGCGGGTTTTTTCGAAAATTGTTTTCTGCCGGAGCAGAAGCCAATTAACGTTTGGAGAACTGCGGACGACGACGTGCTTTACGCAGACCGACTTTCTTACGTTCAACCTGACGAGCGTCACGAGTAACGAAGCCAGCTTTACGCAGTT > CP009273/3370994‑3371277 | gcggttgTCTGAAATAGTTGACTAAGTCCCTTATTGGCGATGTGGTTTTTGTTTTACCGGTCTGTCAGGTGGCAGCAAAAAGCAACTTTCCAGTTTTTACGCTGATTCAGATTTTAGCTATAAAAAAACCCGCCGAAGCGGGTTTTTTc < 1:117992/143‑1 (MQ=255) tAGGCTGAAATAGTTGACTTAGTCCCTTATTGGCGATGTGGTTTTTGTTTTACCTGTCTGT‑AGGTGGCAGCAAAAAGCAACTTTCCAGTTTTTACGCTGATTCAGATTTTAGCTATAAAAAAACCCGCCGAAGCGGGTTTTTTCGaaaa < 2:200896/149‑1 (MQ=255) tCTGAAATAGTTGACTTAGTCCCTTATTGGCGATGTGGTTTTTGTTTTACCTGTCTGTCAGGTGGCAGCAAAAAGCAACTTTCCAGTTTTTACGCTGATTCAGATTTTAGCTATAAAAAAACCCGCCGAAGCGGGTTTTTTCGAAAAtt > 2:14647/1‑149 (MQ=255) gACTTAGTCCCTTATTGGCGATGTGGTTTTTGTTTTACCTGTCTGTCAGGTGGCAGCAAAAAGCAACTTTCCAGTTTTTACGCTGATTCAGATTTTAGCTATAAAAAAACCCGCCGAAGCGGGTTTTTTCGAAAATGGTTTTCTGCCgg > 1:133294/1‑149 (MQ=255) tAGTCCCTTATTGGCGATGTGGTTTTTGTTTTACCTGTCTGTCAGGTGGCAGCAAAAAGCAACTTTCCAGTTTTTACGCTGATTCAGATTTTAGCTATAAAAAAACCCGCCGAAGCGGTTTTTTTCGAAAATTGTTTTCTGCCGGAGCa < 1:102182/149‑1 (MQ=255) ccTTATTGGCGATGTGGTTTTTGTTTTACTTGTCTGTCAGGTGGCAGCAAAAAGCAACTTTCCAGTTTTTACGCTGATTCAGATTTTAGCTATAAAAAAACCCGCCGAAGCGGGTTTTTTCGAAAATTGTTTTCTGCCGGAGCAGAAGc < 2:158079/149‑1 (MQ=255) aTTGGCGATGTGGTGTTTGTTTTACCTGTCTGTCAGGTGGCAGCAAAAGGCAACTTACCAGTTTTTACGCTGATTCAGATTTTAGCTATAAAAAAACCCGCCGAAGCGGTTTTTTTCGATAATTGATTTCTGCCGGAGCAGAGGCCAAt > 2:294231/1‑149 (MQ=255) gTTTTACCTGTCTGTCAGGTGGCAGCAAAAAGCAACTTTCCAGTTTTTACGCTGATTCAGATTTTAGCTATAAAAAAACCCGCCGAAGCGGGTTTTTTCGAAAATTGTTTTCTGCCGGAGCAGAAGCCAATTAACGTTTGGAGAACTGc < 1:14647/149‑1 (MQ=255) gtcAGGTGGCAGCAAAAAGCAACTTTCCAGTTTTTACGCTGATTCAGATTTTAGCTATAAAAAAACCCGCCGAAGCGGGTTTTTTCGAAAATTGTTTTCTGCCGGAGCAGAAGCCAATTAACGTTTGGAGAACTGCGGACGACGACGTg < 2:261934/149‑1 (MQ=255) tcAGGTGGCCGCAAAAAGCAACTTTCCAGTTTTTACGTTGATTCAGATTTTAGCTGTAAAAATACCCGCCGAAGCGGGTTTTTTCGAAAATTGTTTTCTGTCGGAGCCGAAGCCAAGTACCGCTTGGAGAACTGTGGGCGACAAGGTGc > 2:38982/1‑149 (MQ=255) cAGGTGGCAGCAAAAAGCAACTTTCCAGTTTTTACGCTGATTCAGAATTTAGCTATAAAAAAACCCGCCGAAGCGGTTTTTTTCGAAAATTGTTTTCTGCCGGAGCAGAAGCCAATTAACGTTTGGAGAACTGCGGGCGACGACGGGCt > 2:14108/1‑149 (MQ=255) tCCAGTTTTTACGCTGGTTCAGATTTTAGCTATAAAATAACCCGCCGAAGCGGTTTTTTTCGAAAATTGTTTTCTGCCGGAGCAGAAGCCAATTAACGTTTGGAGCACTGCGGGCGACGACGTGCTTTACGCAGACCGACTTTCTTACg < 2:246253/149‑1 (MQ=255) aTTTTAGCTATAAAAAAACCCGCCGAAGCGGGTTTTGTCGAAAATTGTTTTCTGCCGGAGCAGAAGCCAATTAACGTTTGGAGAACTGCGGACGACGACGTGCTTTACGCAGACCGACTTTCTTACGTTCAACCTGACGAGCGTCACGa < 2:9462/149‑1 (MQ=255) tataAAAAAACCCGCCGAAGCGGTTTTTTTCGAAAATTGTTTTCTGCCGGAGCAGAAGCCGATTAACGTTTGGAGAACTGCGGACGACGACGTGCTTTACGCAGACCGACTTTCTTAAGTTCAACCTGACGCGCGGAGCGGGTAACGag > 1:385327/1‑148 (MQ=255) aaaaaaCCCGCCGAAGCGGGTTTTTTCGAAAATTGTTTTCTGCCGGAGCAGAAGCCAATTAACGTTTGGAGAACTGCGGACGACGACGTGCTTTACGCAGACCGACTTTCTTACGTTCAACCTGACGAGCGTCACGAGTAACGAAGCCa > 1:376714/1‑149 (MQ=255) cgcggAAGCGGGTTTTTTCGAAAACTGTTGTCTGCCGGAGCAGAAGCCAATTAACGTTTGGAGAACTGCGGACGACGACGTGCTTTACGCGGACCGACTTTCTTACGTTCAACCTGACGAGCGTCACGAGTAACGAAGCCAGCTTTACg < 2:14560/149‑1 (MQ=255) aaGCGGGTTTTTTCGAAAATTGTTTTCTGCCGGAGCAGAAGCCAATTAACGTTTGGAGAACTGCGGACGACGACGTGCTTTACGCAGACCGACTTTCTTACGTTCAACCTGACGAGCGTCACGAGTAACGAAGCCAGCTTTACGCAGtt < 1:853/149‑1 (MQ=255) | GCTTTAGTCTGAAATAGTTGACTTAGTCCCTTATTGGCGATGTGGTTTTTGTTTTACCTGTCTGTCAGGTGGCAGCAAAAAGCAACTTTCCAGTTTTTACGCTGATTCAGATTTTAGCTATAAAAAAACCCGCCGAAGCGGGTTTTTTCGAAAATTGTTTTCTGCCGGAGCAGAAGCCAATTAACGTTTGGAGAACTGCGGACGACGACGTGCTTTACGCAGACCGACTTTCTTACGTTCAACCTGACGAGCGTCACGAGTAACGAAGCCAGCTTTACGCAGTT > CP009273/3370994‑3371277 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |