Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 4,355,090 | Δ1 bp | 11.8% | coding (246/681 nt) | modB → | molybdate ABC transporter permease ModB |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 4,355,090 | 0 | G | . | 11.8% | 97.1 / 9.0 | 34 | coding (246/681 nt) | modB | molybdate ABC transporter permease ModB |
Reads supporting (aligned to +/- strand): ref base G (16/14); new base . (2/2); total (18/16) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.05e-01 |
GCCTGGTGGCTCGCGCGCACGCGCTCTTGGCTGCGCGGGCCGGTGGGTGCGGTGGTGGCACTGCCGCTGGTACTGCCACCGACGGTGATCGGTTTTTACCTGTTGATCGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAG‑GCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCCTGGTGATTGGCTCGACGGTGTACTCCATGCCGTTCGTGGTGCAACCCCTGCAGAATGCCTTTGGCGCGATCGGCCAGCGCCCCCTGGAAGTCGCCGCTACTCTGCGA > NC_002947/4354941‑4355234 | gCCTGGTGGCTAGAGCGCACGCGCTCTTGGCTGCGCGGGCCGGTGGGTGCGGTGGTGGCACTGCCGCTGGTACTGCCACCTACGGTGATCGGTTTTTACCTGTTGATCGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGgcg > 2:612567/1‑150 (MQ=255) gCTCGCGCGCACGCGCTCTTGGCTGCGCGGGCCGGTGGGTGCGGTGGTGGCACTGCCGCTGGTACTGCCACCTACGGTGATCGGTTTTTACCTGTTGATCGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAG‑GCCt < 1:612567/150‑1 (MQ=255) tCGCGCGCACGCGCTCTTGGCTGCGCGGGCCGGTGGGTGCGGTGGTGGCACTGCCGCTGGTACTGCCACCGACGGTGATCGGTTTTTACCTGTTGATCGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAG‑GCCTgg > 1:360214/1‑150 (MQ=255) cgcACGCGCTCTTGGCTGCGCGGGCCGGTGGGTGCGGTGGTGGCACTGCCGCTGGTACTGCCACCGACGGTGATCGGTTTTTACCTGTTGATCGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGC‑CTAGCGCCTgggtggc < 2:191827/150‑5 (MQ=255) cACGCGCTCTTGGCTGCGCGGGCCGGTGGGTGCGGTGGTGGCACTGCCGCTGGTACTGCCACCGACGGTGATCGGTTTTTACCTGTTGATCGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAG‑GCCTGGGCAgtgt < 2:572693/150‑1 (MQ=255) gcgGGCCGGTGGGTGCGGTGGTGGCACTGCCGCTGGTACTGCCACCGACGGTGATCGGTTTTTACCTGTTGATCGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAG‑GCCTGGGCAGTGTGGTGTTCAGCTTCACAg < 2:262146/150‑1 (MQ=255) ggtgggtgCGGTGGTGGCACTGCCGCTGGTACTGCCACCGACGGTGATCGGTTTTTACCTGTTGATCGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAG‑GCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCCTGGt < 1:51897/150‑1 (MQ=255) gtggCACTGCCGCTGGTACTGCCACCGACGGTGATCGGTTTTTACCTGTTGATCGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAG‑GCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCCTGGTGATTGGCTCGACg < 1:165397/150‑1 (MQ=255) gCACTGCCGCTGGTACTGCCACCGACGGTGATCGGTTTTTACCTGTTGATCGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAG‑GCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCCTGGTGATTGGCTCGACGgtg > 1:433584/1‑150 (MQ=255) gCACTGCCGCTGGTACTGCCACCGACGGTGATCGGTTTTTACCTGTTGATCGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAG‑GCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCCTGGTGATTGGCTCGACGgtg < 2:360214/150‑1 (MQ=255) tGCCACCGACGGTGATCGGTTTTTACCTGTTGATCGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAG‑GCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCCTGGTGATTGGCTCGACGGTGTACTCCATGCCGTTCg > 1:60136/1‑150 (MQ=255) aTCGGTTTTTACCTGTTGATCGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAG‑GCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCCTGGTGATTGGCTCGACGGTGTACTCCATGCCGTTCGTGGTGCAACCCCTg > 2:103109/1‑150 (MQ=255) tCGGTTTTTACCTGTTGATCGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAG‑GCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCCTGGTGATTGGCTCGACGGTGTACTCCATGCCGTTCGTGGTGCAACCCCTgg > 1:234388/1‑149 (MQ=255) tCGGTTTTTACCTGTTGATCGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAG‑GCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCCTGGTGATTGGCTCGACGGTGTACTCCATGCCGTTCGTGGTGCAACCCCTGc > 2:132518/1‑150 (MQ=255) cGGTTTTTACCTGTTGATCGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAG‑GCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCCTGGTGATTGGCTCGACGGTGTACTCCATGCCGTTCGTGGTGCAACCCCTGCa > 1:415740/1‑150 (MQ=255) gTTTTTACCTGTTGATCGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAG‑GCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCCTGGTGATTGGCTCGACGGTGTACTCCATGCCGTTCGTGGTGCAACCCCTGCAGa > 1:422723/1‑150 (MQ=255) ttACCTGTTGATCGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAG‑GCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCCTGGTGATTGGCTCGACGGTGTACTCCATGCCGTTCGTGGTGCAACCCCTGCAGAATGc > 1:456412/1‑150 (MQ=255) cTGTTGATCGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAG‑GCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCCTGGTGATTGGCTCGACGGTGTACTCCATGCCGTTCGTGGTGCAACCCCTGCAGAATGCCttt < 2:60136/150‑1 (MQ=255) aTCGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAG‑GCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCCTGGTGATTGGCTCGACGGTGTACTCCATGCCGTTCGTGGTGCAACCCCTGCAGAATGCCTTTGGcgcg > 2:417077/1‑150 (MQ=255) cGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCACACAGGCGCTAG‑GCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCCTGGTGATTGGCTCGACGGTGTACTCCATGCCGTTCGTGGTGCAACCCCTGCAGAATGCCTTTGGCGCGAt < 1:301157/150‑1 (MQ=255) cccTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAG‑GCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCCTGGTGATTGGCTCGACGGTGTACTCCATGCCGTTCGTGGTGCAACCCCTGCAGAATGCCTTTGGCGCGATCg > 2:396738/1‑150 (MQ=255) cTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAG‑GCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCCTGGTGATTGGCTCGACGGTGTACTCCATGCCGTTCgtg > 2:551284/1‑114 (MQ=255) cTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAG‑GCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCCTGGTGATTGGCTCGACGGTGTACTCCATGCCGTTCgtg < 1:551284/114‑1 (MQ=255) ttGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAG‑GCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCCTGGTGATTGGCTCGACGGTGTACTCCATGCCGTTCGTGGTGCAACCCCTGCAGAATGCCTTTGGCGCGATCGGcc > 1:89381/1‑150 (MQ=255) ccccGCACTGCTGGCTTTGTCATCCCACCCAGGCGCCAG‑GCCTGGGCAGTGTGGTGTTCAGCTTTATAGGCCTGGTGATTGGCTCTAAGGTCTACCCCATGCCTTTCGTGGTGCAACCCCTGCAGAATGCCTTTTGCGCGATCGGCCAgc < 2:234388/149‑1 (MQ=255) cccGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAG‑GCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCCTGGTGATTGGCTCGACGGTGTACTCCATGCCGTTCGTGGTGCAACCCCTGCAGAATGCCTTTGGCGCGATCGGCCAgcg > 1:300582/1‑150 (MQ=255) gCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAG‑GCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCCTGGTGATTGGCTCGACGGTGTACTCCATGCCGTTCGTGGTGCAACCCCTGCAGAATGCCTTTGGCGCGATCGGCCAGCGccc < 2:422723/150‑1 (MQ=255) tCAGGCCACCCAGGCGCTAG‑GCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCCTGGTGATTGGCTCGACGGTGTACTCCATGCCGTTCGTGGTGCAACCCCTGCAGAATGCCTTTGGCGCGATCGGCCAGCGCCCCCTGGAAGTcgccgc < 1:132518/150‑1 (MQ=255) tCAGGCCACCCAGGCGCTAG‑GCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCCTGGTGATTGGCTCGACGATGTACTCCATGCCGTTCGTGGTGCAACCCCTGCAGAATGCCTTTGGCGCGATCGGCCAGCGCCCCCTGGAAGTcgccgc < 2:257942/150‑1 (MQ=255) cacactgcccAGGC‑CTAGCGCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCCTGGTGATTGGCTCGACGGTGTACTCCATGCCGTTCGTGGTGCAACCCCTGCAGAATGCCTTTGGCGCGATCGGCCAGCGCCCCCTGGAAGTCGCCGCt > 2:21039/8‑150 (MQ=255) tgcccAGGC‑CTAGCGCCTGGGTAGTGTGGTGTTCAGCTTCACAGGCCTGGTGATTGGCTCGACGGTGTACTCCATGCCGTTCGTGGTGCAACCCCTGCAGAATGCCTTTg < 1:145651/108‑1 (MQ=255) tgcccAGGC‑CTAGCGCCTGGGTAGTGTGGTGTTCAGCTTCACAGGCCTGGTGATTGGCTCGACGGTGTACTCCATGCCGTTCGTGGTGCAACCCCTGCAGAATGCCTTTg > 2:145651/3‑110 (MQ=255) cACCCAGGCGCTAG‑GCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCCTGGTGATTGGCTCGACGGTGTACTCCATGCCGTTCGTGGTGCAACCCCTGCAGAATGCCTTTGGCGCGATCGGCCAGCGCCCCCTGGAAGTCGCCGCTActct > 1:437113/1‑150 (MQ=255) aCCCAGGTGCTAG‑GCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCCTGGTGATTGGCTCGACGGTGTACTCCATGCCGTTCGTGGTGCAACCCCTGCAGAATGCCTTTGGCGCGATCGGCCAGCGCCCCCTGGAAGTCGCCGCTACTCTg > 1:363685/1‑150 (MQ=255) cAGGCGCTAG‑GCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCCTGGTGATTGGCTCGACGGTGTACTCCATGCCGTTCGTGGTGCAACCCCTGCAGAATGCCTTTGGCGCGATCGGCCAGCGCCCCCTGGAAGTCGCCGCTACTCTGCGa < 2:89381/150‑1 (MQ=255) | GCCTGGTGGCTCGCGCGCACGCGCTCTTGGCTGCGCGGGCCGGTGGGTGCGGTGGTGGCACTGCCGCTGGTACTGCCACCGACGGTGATCGGTTTTTACCTGTTGATCGCCCTTGGCCCGCACGGCTGGCTTGGTCAGGCCACCCAGGCGCTAG‑GCCTGGGCAGTGTGGTGTTCAGCTTCACAGGCCTGGTGATTGGCTCGACGGTGTACTCCATGCCGTTCGTGGTGCAACCCCTGCAGAATGCCTTTGGCGCGATCGGCCAGCGCCCCCTGGAAGTCGCCGCTACTCTGCGA > NC_002947/4354941‑4355234 |
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: ‑ |