Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 3,105,669 | G→C | 12.0% | L123V (CTG→GTG) | PP_2716 ← | arsenate reductase |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 3,105,669 | 0 | G | C | 12.0% | 111.4 / 7.7 | 42 | L123V (CTG→GTG) | PP_2716 | arsenate reductase |
Reads supporting (aligned to +/- strand): ref base G (17/19); new base C (2/3); total (19/23) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.65e-01 |
GAGGTCGAGGTTAGGCAGTTGCTCGGACATGTTTTCTTCCTGTTCAAAGCGAGCCGATGCGATCCAGCGCGTGCTTGAGCTGCTCGCGGTCGAGGTGATCGAACGGCAGCGCGAAGAAGGCCTGGCAGCGAGACTCGATACGGGCCAGTGTGGCCCGGAAAGCCGCATCGACCGTAGCTTCGTCACCCACGACATCGGAGGGGTCTTCCAGCCCCCAATGCGACTTCAACGCAGGGCCGAAATACACCGGGCAGGCTTCGCCCGCTGCCTTGTCGCAGACGGTGATGACGATGT > NC_002947/3105522‑3105815 | gAGGTCGAGGTTAGGCAGTTGCTCGGACATGTTTTCTTCCTGTTCAAAGCGAGCCGATGCGATCCAGCGCGTGCTTGAGCTGCTCGCGGTCGAGGTGATCGAACGGCAGCGCGAAGAAGGCCTGGCAGCGAGACTCGATACGGGCCAgt > 4:168629/1‑149 (MQ=255) aGGTCGAGGTTAGGCAGTTGCTCGGACATGTTTTCTTCCTGTTCAAAGCGAGCCGATGCGATCCAGCGCGTGCTTGAGCTGCTCGCGGTCGAGGTGATCGAACGGCAGCGCGAAGTAGGCCTGGCAGCGAGACTCGATACGGGCCAgtg < 1:260556/149‑1 (MQ=255) tCGAGGTTAGGCAGTTGCTCGGACATGTTTTCTTCCTGTTCAAAGCGAGCCGATGCGATCCAGCGCGTGCCTGAGCTGCTCGCGGTCGAGGTGATCGAACGGCAGCGCGAAGAAGGCCTGGCAGCGAGACTCGATACGGGCCAGTGTgg > 7:377329/1‑149 (MQ=255) ttAGGCAGTTGCTCGGACATGTTTTCTTCCTGTTCAAAGCGAGCCGATGCGATCCAGCGCGTGCTTGAGCTGCTCGCGGTCGAGGTGATCGAAAGGCAGCGCGAAGAAGGCCTGGCAGCGAGCCTCGATACGGGCCAGTGTGGCCCgg > 8:375738/1‑148 (MQ=255) gCAGTTGCTCGGACATGTTTTCTTCCTGTTCAAAGCGAGCCGATGCGATCCAGCGCGTGCTTGAGCTGCTCGCGGTCGAGGTGATCGAACGGCAGCGCGAAGAAGGCCTGGCAGCGAGACTCGATACGGGCCACACTGGCCCgtatcga < 6:186208/149‑7 (MQ=255) gCAGTTGCTCGGACATGTTTTCTTCCTGTTCAAAGCGAGCCGATGCGATCCAGCGCGTGCTTGAGCTGCTCGCGGTCGAGGTGATCGAACGGCAGCGCGAAGAAGGCCTGGCAGCGAGACTCGATACGGGCCACACTGGCCCgtatcga < 4:198204/149‑7 (MQ=255) gCTCGGACATGTTTTCTTCCTGTTCAAAGCGAGCCGATGCGATCCAGCGCGTGCTTGAGCTGCTCGCGGTCGAGGTGATCGAACGGCAGCGCGAAGAAGGCCTGGCAGCGAGACTCGATACGGGCCAGTGTGGCCCGGAAAGCCGCAt < 7:206368/148‑1 (MQ=255) cTCGGACATGTTTTCTTCCTGTTCAAAGCGAGCCGATGCGATCCAGCGCGTGCTTGAGCTGCTCGCGGTCGAGGTGATCGAACGGCAGCGCGAAGAAGGCCTGGCAGCGAGACTCGATACGGGCCAGTGTGGCCCGGAAAGCCGCATc < 3:366470/148‑1 (MQ=255) tCGGACATGTTTTCTTCCTGTTCAAAGCGAGCCGATGCGATCCAGCGCGTGCTTGAGCTGCTCGCGGTCGAGGTGATCGAACGGCAGCGCGAAGAAGGCCTGGCAGCGAGACTCGATACGGGCCAGTGTGGCCCGGAAAGCCGCATCGa < 8:296911/149‑1 (MQ=255) ggACATGTTTTCTTCCTGTTCAAAGCGAGCCGATGCGATCCAGCGCGTGCTTGAGCTGCTCGCGGTCGAGGTGATCGAACGGCAGCGCGAAGAAGGCCTGGCAGCGAGACTCGATACGGGCCAGTGTGGCCCGGAAAGCCGCATCGAcc > 2:313039/1‑149 (MQ=255) gaaaTGTTTTCTTCCTGTTCAAAGCGAGCCGATGCGATCCAGCGCGTGCTTGAGCTGCTCGCGGTCGAGGTGATCGAACGGCAGCGCGAAGAAGGCCTGGCAGCGAGACTAGATACGGGCCAGTGTGGCCCGGAAAGCCTCATCGACCg > 8:188792/4‑149 (MQ=255) gACATGTTTTCTTCCTGTTCAAAGCGAGCCGATGCGATCCAGCGCGTGCTTGAGCTGCTCGCGGTCGAGGTGATCGAACGGCAGCGCGAAGAAGGCCTGGCAGCGAGACTCGATACGGGCCAGTGTGGCCCGGAAAGCCGCATCGACCg > 8:194195/1‑149 (MQ=255) tcttcCTGTTCAAAGCGAGCCGATGCGATCCAGCGCGTGCTTGAGCTGCTCGCGGTCGAGGTGATCGAACGGCAGCGCGAAGAAGGCCTGGCAGCGAGACTCGATACGGGCCAGTGTGGCCCg < 7:192150/123‑1 (MQ=255) tcttcCTGTTCAAAGCGAGCCGATGCGATCCAGCGCGTGCTTGAGCTGCTCGCGGTCGAGGTGATCGAACGGCAGCGCGAAGAAGGCCTGGCAGCGAGACTCGATACGGGCCAGTGTGGCCCg > 8:192150/1‑123 (MQ=255) ttCAAAGCGAGCCGATGCGATCCAGCGCGTGCTTGAGCTGCTCGCGGTCGAGGTGATCGAACGGCAGCGCGAAGAAGGCCTGGCAGCGAGACTCGATACGGGCCAGTGTGGCCCGGAAAGCCGCATCGACCGTAGCTTCGTCACCCAc > 5:323021/1‑148 (MQ=255) ttCAAAGCGAGCCGATGCGATCCAGCGCGTGCTTGAGCTGCTCGCGGTCGAGGTGATCGAACGGCAGCGCGAAGAAGGCCTGGCAGCGAGACTCGATACGGGCCAGTGTGGCCCGGAAAGCCGCATCGACCGTAGCTTCGTCACCCAc > 6:66212/1‑148 (MQ=255) cAAAGCGAGCCGATGCGATCCAGCGCGTGATTGGGCTGCTCGCGGTCGAGGTGATCGAACGGCAGCGCGAAGAAGGCCTGGCAGCGAGACTCGATACGGGCCAGTGTGGCCCGGAAAGCCGCATCGACCGTAGCTTCGTCACCCACGa > 1:235335/1‑148 (MQ=255) aGCGAGCCGATGCGATCCAGCGCGTGCTTGAGCTGCTCGCGGTCGAGGTGATCGAACGGCAGCGCGAAGAAGGCCTGGCAGCGAGACTCGATACGGGCCAGTGTGGCCCGGAAAGCCGCATCGACCGTAGCTTCGTCACCCACGACATc > 3:105607/1‑149 (MQ=255) ccGATGCGATCCAGCGCGTGCTTGAGCTGCTCGCGGTCGAGGTGATCGAACGGCAGCGCGAAGAAGGCCTGGCAGCGAGACTCGATACGGGCCAGTGTGGCCCGGAAAGCCGCATCGACCGTAGCTTCGTCACCCACGACATCGGAggg < 3:168629/149‑1 (MQ=255) cGATCCAGCGCGTGCTTGAGCTGCTCGCGGTCGAGGTGATCGAACGGCAGCGCGAAGAAGGCCTGGCAGCGAGACTCGATACGGGCCAGTGTGGCCCGGAAAGCCGCATCGACCGTAGCTTCGTCACCCACGACATCGGAGGGGTCTTc < 1:313039/149‑1 (MQ=255) gATCCAGCGCGTGCTTGAGCTGCTCGCGGTCGAGGTGATCGAACGGCAGCGCGAAGAAGGCCTGGCAGCGAGACTCGATACGGGCCAGTGTGGCCCGGAAAGCCGCATCGACCGTAGCTTCGTCACCCACGACATCGGAGGGGTCTTc > 6:409550/1‑148 (MQ=255) aTCCAGCGCGTGCTTGAGCTGCTCGCGGTCGAGGTGATCGAACGGCAGCGCGAAGAAGGCCTGGCAGCGAGACTCGATACGGGCCAGTGTGGCCCGGAAAGCCGCATCGACCGTAGCTTCGTCACCCACGACATCGGAGGGGTCTTCCa > 3:101657/1‑149 (MQ=255) gcgcGTGCTTGAGCTGCTCGCGGTCGAGGTGATCGAACGGCAGCGCGAAGAAGGCCTGGCAGCGAGACTCGATACGGGCCAGTGTGGCCCGGAAAGCCGCATCGACCGTAGCTTCGTCACCCACGACATCGGAGGGGTCTTCCAGcccc < 6:66507/149‑1 (MQ=255) gcgTGCTTGAGCTGCTCGCGGTCGAGGTGATCGAACGGCAGCGCGAAGAAGGCCTGGCAGCGAGACTCGATACGGGCCAGTGTGGCCCGGAAAGCCGCATCGACCGTAGCTTCGTCACCCACGACATCGGAGGGGTCTTCCAGCCCCCa < 5:66212/149‑1 (MQ=255) cgTGCTTGAGCTGCTCGCGGTCGAGGTGATCGAACGGCAGCGCGAAGAAGGCCTGGCAGCGAGACTCGATACGGGCCAGTGTGGCCCGGAAAGCCGCATCGACCGTAGCTTCGTCACCCACGACATCGGAGGGGTCTTCCAGCCCCCaa < 2:80394/149‑1 (MQ=255) gAGCTGCTCGCGGTCGAGGTGATCGAACGGCAGCGCGAAGAAGGCCTGGCAGCGAGACTCGATACGGGCCAGTGTGGCCCGGAAAGCCGCATCGACCGTAGCTTCGTCACCCACGACATCGGAGGGGTCTTCCAGCCCCCAATGCGa < 5:409550/147‑1 (MQ=255) gctgctCGCGGTCGAGGTGATCGAACGGCAGCGCGAAGAAGGCCTGGCAGCGAGACTCGATACGGGCCAGTGTGGCCCGGAAAGCCGCATCGACCGTAGCTTCGTCACCCACGACATCGGAGGGGTCTTCCAGCCCCCAATGCGACTTc < 4:105607/149‑1 (MQ=255) ctCGCGGTCGAGGTGATCGAACGGCAGCGCGAAGAAGGCCTGGCAGCGAGACTCGATACGGGCCAGTGTGGCCCGGAAAGCCGCATCGACCGTAGCTTCGTCACCCACGACATCGGAGGGGTCTTCCAGCCCCCAATGCGACTTCAAc < 8:377329/148‑1 (MQ=255) gTCGAGGTGATCGAACGGCAGCGCGAAGAAGGCCTGGCAGCGAGACTCGATACGGGCCAGTGTGGCCCGGAAAGCCGCATCGACCGTAGCTTCGTCACCGACGACATCGGAGGTGTCTTCCAGCCCCCAATGCGACTTCAACGCAGGGc < 7:375738/149‑1 (MQ=255) gAGGTGATCGAACGGCAGCGCGTAGAAGGCCTGGCAGCGAGACTCGATACGGGCCAGTGTGGCCCGGAAAGCCGCATCGACCGTAGCTTCGTCACCCACGACATCGGAGGGGTCTTCCAGCCCCCAATGCGACTTCAACGCAGGGCCg < 6:37188/148‑1 (MQ=255) gTGATCGAACGGCAGCGCGAAGAAGGCCTGGCAGCGAGACTCGATACGGGCCAGTGTGGCCCGGAAAGCCGCATCGACCGTAGCTTCGTCACCCACGACATCGGAGGGGTCTTCCAGCCCCCAATGCGACTTCAACGCAGGGCCGaaa < 4:178412/148‑1 (MQ=255) gTGATCGAACGGCAGCGCGAAGAAGGCCTGGCAGCGAGACTCGATACGGGCCAGTGTGGCCCGGAAAGCCGCATCGACCGTAGCTTCGTCACCCACGACATCGGAGGGGTCTTCCAGCCCCCAATGCGACTTCAACGCAGGGCCGAAAt < 7:194195/149‑1 (MQ=255) tGATCGAACGGCAGCGCGAAGAAGGCCTGGCAGCGAGACTCGATACGGGCCATTGTGGCCCGGAAAGCCGCATCGACCGTAGCTTCGTCAACCACGACATCGGAGGTGTCTTCCAGCCCCCAATGCGACTTCAACGCAGGGCCGAAAt < 7:188792/148‑1 (MQ=255) cAGCGCGAAGAAGGCCTGGCAGCGAGACTCGATACGGGCCAGTGTGGCCCGGAAAGCCGCATCGACCGTAGCTTCGTCACCCACGACATCGGAGGGGTCTTCCAGCCCCCAATGc > 7:171089/1‑115 (MQ=255) cAGCGCGAAGAAGGCCTGGCAGCGAGACTCGATACGGGCCAGTGTGGCCCGGAAAGCCGCATCGACCGTAGCTTCGTCACCCACGACATCGGAGGGGTCTTCCAGCCCCCAATGc < 8:171089/115‑1 (MQ=255) gaagaaGGCCTGGCAGCGAGACTCGATACGGGCCAGTGTGGCCCGGAAAGCCGCATCGACCGTAGCTTCGTCACCCACGACATCGGAGGGGTCTTCCAGCCCCCAATGCGACTTCAACGCAGGGCCGAAATACACCGGGCAGGCTTCGc > 8:345665/1‑149 (MQ=255) gCAGCGAGACTCGATACGGGCCAGTGTGGCCCGGAAAGCCGCATCGACCGTAGCTTCGTCACCCACGACATCGGAGGGGTCTTCCAGCCCCCAATGCGACTTCAACGCAGGGCCGAAATACACCGGGCAGGCTTCGCCCGCTGCCTTGt > 2:165019/1‑149 (MQ=255) tCGATACGGGCCAGTGTGGCCCGGAAAGCCGCATCGACCGTAGCTTCGTCACCCACGACATCGGAGGGGTCTTCCAGCCCCCAATGCGACTTCAACGCAGGGCCGAAATACACCGGGCAGGCTTCGCCCGCTGCCTTGTc < 7:320576/140‑1 (MQ=255) tCGATACGGGCCAGTGTGGCCCGGAAAGCCGCATCGACCGTAGCTTCGTCACCCACGACATCGGAGGGGTCTTCCAGCCCCCAATGCGACTTCAACGCAGGGCCGAAATACACCGGGCAGGCTTCGCCCGCTGCCTTGTc > 8:320576/1‑140 (MQ=255) tttccGGGCCACACTGGCCCGTAAAGCCGCATCGACCGTAGCTTCGTCACCCACGACATCGGAGGGGTCTTCCAGCCCCCAATGCGACTTCAACGCAGGGcc > 8:185429/5‑102 (MQ=255) tttccGGGCCACACTGGCCCGTAAAGCCGCATCGACCGTAGCTTCGTCACCCACGACATCGGAGGGGTCTTCCAGCCCCCAATGCGACTTCAACGCAGGGcc < 7:185429/98‑1 (MQ=255) tttccGGGCCACACTGGCCCGGAAAGCCGCATCGACCGTAGCTTCGTCACCCACGACATCGGAGGGGTCTTCCAGCCCCCAATGCGACTTCAACGCAGGGCCGAAATACACCGGGCAGGCTTCGCCCGCTGCCTTGTCGCAGACGGtga > 1:124909/5‑149 (MQ=255) gATACGGGCCAGTGTGGCCCGGAAAGCCGCATCGACCGTAGCTTCGTCACCCACGACATCGGAGGGGTCTTCCAGCCCCCAATGCGACTTCAACGCAGGGCCGAAATACACCGGGCAGGCTTCGCCCGCTGCCTTGTCGCAGACGGtg < 8:78552/148‑1 (MQ=255) cAGTGTGGCCCGGAAAGCCGCATCGACCGTAGCTTCGTCACCCACGACATCGGAGGGGTCTTCCAGCCCCCAATGCGACTTCAACGCAGGGCCGAAATACACCGGGCAGGCTTCGCCCGCTGCCTTGTCGCAGACGGTGATGACGATGt > 7:136781/1‑149 (MQ=255) | GAGGTCGAGGTTAGGCAGTTGCTCGGACATGTTTTCTTCCTGTTCAAAGCGAGCCGATGCGATCCAGCGCGTGCTTGAGCTGCTCGCGGTCGAGGTGATCGAACGGCAGCGCGAAGAAGGCCTGGCAGCGAGACTCGATACGGGCCAGTGTGGCCCGGAAAGCCGCATCGACCGTAGCTTCGTCACCCACGACATCGGAGGGGTCTTCCAGCCCCCAATGCGACTTCAACGCAGGGCCGAAATACACCGGGCAGGCTTCGCCCGCTGCCTTGTCGCAGACGGTGATGACGATGT > NC_002947/3105522‑3105815 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |