| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | NC_002947 | 5,072,314 | A→G | 11.4% | H179H (CAT→CAC) | PP_5689 ← | hypothetical protein |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_002947 | 5,072,314 | 0 | A | G | 11.4% | 88.1 / 3.5 | 35 | H179H (CAT→CAC) | PP_5689 | hypothetical protein |
| Reads supporting (aligned to +/- strand): ref base A (16/15); new base G (2/2); total (18/17) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.67e-01 | |||||||||||
GCGAGTCGGGTAGACCGCATAGATGCCCAGCTCTACCGAACACTACTCGGGCAATATCTCCACCAGATTTCTAAAGAAGCAAACCCATACTCTCGGCCCACCGTCCAAGGCTCAGCTCAAAGCTTTCGTTTGCTTTGCTGCCAAATGCGCCGCAGTCCATAGCTGGTCTGCCAACCACTCCGCTTCAGCAGGACTGAGCACCCACTCTACCCAGCCAAACTTAAGCATAATTGATTCGCCGTAGGGCCTGATAGTCAAATCATAGATGGCCCCTGCTTCGTCCTTGAATCC > NC_002947/5072170‑5072460 | gCGAGTCGGGTAGACCGCATAGATGCCCAGCTCTACCGAACACTACTCGGGCAATATCTCCACCAGATTTCTAAAGAAGCAAACCCATACTCTCGGCCCACCGTCCAAGGCTCAGCTCAAAGCTTTCGTTTGCTTTGCTGCCAAATgcg > 4:225208/1‑149 (MQ=255) gggTAGACCGCATAGATGCCCAGCTCTACCGAACACTACTCGGGCAATATCTCCACCAGATTTCTAAAGAAGCAAACCCATACTCTCGGCCCACCGTCCAAGGCTCAGCTCAAAGCTTTCGTTTGCTTTGCTGCCAAATGCGCCGCAGt < 3:294635/149‑1 (MQ=255) aTGCCCAGCTCTACCGAACACTACTCGGGCAATATCTCCACCAGATTTCTAAAGAAGCAAACCCATACTCTCGGCCCACCGTCCAAGGCTCAGCTCAAAGCTTTCGTTTGCTTTGCTGCCAAATGCGCCGCAGTCCATAGCTGGTCTGc < 2:178374/149‑1 (MQ=255) tGCCCAGCTCTACCGAACACTACTGGGTGAATATCTCGACCAGATTTCTAAAGAAGCAAACCCATAGTCTCGGCCCACCGTCCAAGGCTCTGCTCAAAGCTTTCGTTTGCTTTGCTGCCAAATGCGGCGCAGTACATAGCTTGTCTGcc < 2:187625/149‑1 (MQ=255) ccAGCTCTACCGAACACTACTCGGGCAATATCTCCACCAGATTTCTAAAGAAGCAAACCCATACTCTCGGCCCACCGTCCAAGGCTCAGCTCAAAGCTTTCGTTTGCTTTGCTGCCAAATGCGCCGCAGTCCATAGCTGGTCTGCCAAc < 7:88657/149‑1 (MQ=255) tACCGAACACTACTCGGGCAATATCTCCACCAGATTTCTAAAGAAGCAAACCCATACTCTCGGCCCACCGTCCAAGGCTCAGCTCAAAGCTTTCGTTTGCTTTGCTGCCAAATGCGCCGCAGTCCATAGCTGGTCTGCCAACCACTc < 5:247204/147‑1 (MQ=255) cAATATCTCCACCAGATTTCTAAAGAAGCAAACCCATACTCTCGGCCCACCGTCCAAGGCTCAGCTCAAAGCTTTCGTTTGCTTTGCTGCCAAATGCGCCGCAGTCCATAGCTCGTCTGCCAACCACTCCGCTTCAGCAGGAGTGAGCa > 8:120325/1‑149 (MQ=255) caccaGATTTCTAAAGAAGCAAACCCATACTCTCGGCCCACCGTCCAAGGCTCAGCTCAAAGCTTTCGTTTGCTTTGCTGCCAAATGCGCCGCAGTCCATAGCTGGTCTGCCAACCACTCCGCTTCAGCAGGACTGAGCACCCAc < 6:82845/145‑1 (MQ=255) agCAAACCCATACTCTCGGCCCACCGTCCAAGGCTCAGCTCAAAGCTTTCGTTTGCTTTGCTGCCAAATGCGCCGCAGTCCATAGCTGGTCTGCCAACCACTCCGCTTCAGCAGGACTGAGCACCCACTCTACCCAGCCAAACTTAAg < 3:225208/148‑1 (MQ=255) ccATACTCTCGGCCCACCGTCCAAGGCTCAGCTCAAAGCTTTCGTTTGCTTTGCTGCCAAATGCGCCGCAGTCCATAGCTGGTCTGCCAACCACTCCGCTTCAGCAGGACTGAGCCCCCACTCTACCCAGCCAAACTTAAGCATAAtt > 3:250872/1‑148 (MQ=255) aTACTCTCGGCCCACCGTCCAAGGCTCAGCTCAAAGCTTTCGTTTGCTTTGCTGCCAAATGCGCCGCAGTCCATAGCTGGTCTGCCAACCACTCCGCTTCAGCAGGACTGAGCCCCCACTCTACCCAGCCAAACTTAAGCATAATTGAt > 4:106200/1‑149 (MQ=255) gCCCACCGTCCAAGGCTCAGCTCAAAGCTTTCGTTTGCTTTGCTGCCAAATGCGCCGCAGTCCATAGCTGGTCTGCCAACCACTCCGCTTCAGCAGGACTGAGCACCCACTCTACCCAGCCAAACTTAAGCATAATTGATTCGCCGTAg > 1:70230/1‑149 (MQ=255) cACCGTCCAAGGCTCAGCTCAAAGCTTTCGTTTGCTTTGCTGCCAAGTGCGCCGCAGTCCATAGCTGGTCTGCCAACCACTCCGCTTcagc < 6:68491/91‑1 (MQ=255) cACCGTCCAAGGCTCAGCTCAAAGCTTTCGTTTGCTTTGCTGCCAAGTGCGCCGCAGTCCATAGCTGGTCTGCCAACCACTCCGCTTcagc < 6:63813/91‑1 (MQ=255) cACCGTCCAAGGCTCAGCTCAAAGCTTTCGTTTGCTTTGCTGCCAAGTGCGCCGCAGTCCATAGCTGGTCTGCCAACCACTCCGCTTcagc > 5:68491/1‑91 (MQ=255) cACCGTCCAAGGCTCAGCTCAAAGCTTTCGTTTGCTTTGCTGCCAAGTGCGCCGCAGTCCATAGCTGGTCTGCCAACCACTCCGCTTcagc > 5:63813/1‑91 (MQ=255) cGTCCAAGGCTCAGCTCAAAGCTTTCGTTTGCTTTGCTGCCAAATGCGCCGCAGTCCATAGCTGGTCTGCCAACCACTCCGCTTCAGCAGGACTGAGCACCCACTCTACCCAGCCaaa < 6:222896/118‑1 (MQ=255) cGTCCAAGGCTCAGCTCAAAGCTTTCGTTTGCTTTGCTGCCAAATGCGCCGCAGTCCATAGCTGGTCTGCCAACCACTCCGCTTCAGCAGGACTGAGCACCCACTCTACCCAGCCaaa > 5:222896/1‑118 (MQ=255) cGTCCAAGGCTCAGCTCAAAGCTTTCGTTTGCTTTGCTGCCAAATGCCCCCCAGTCCATAGCTGGTCTGCCAACCACTCCGCTTCAGCAGGACTGAGCACACAATCTAACCAGCCAAACTTAAGCATAATTGATTCGCCGTAGGGcat > 7:10530/1‑146 (MQ=255) agctcaAAGCTTTCGTTTGCTTTGCTGCCAAATGCGCCGCAGTCCATAGCTGGTCTGCCAACCACTCCGCTTCAGCAGGACTGAGCACCCACTCTACCCAGCCAAACTTAAGCATAATTGATTCGCCGTAGGGCCTg < 5:295555/137‑1 (MQ=255) agctcaAAGCTTTCGTTTGCTTTGCTGCCAAATGCGCCGCAGTCCATAGCTGGTCTGCCAACCACTCCGCTTCAGCAGGACTGAGCACCCACTCTACCCAGCCAAACTTAAGCATAATTGATTCGCCGTAGGGCCTg > 6:295555/1‑137 (MQ=255) ctcaAAGCTTTCGTTTGCTTTGCTGCCAAATGCGCCGCAGTCCATAGCTGGTCTGCCAACCACTCCGCTTCAGCAGGACTGAGCACCCACTCTACCCAGCCAAACTTAAGCATAATTGATTCGCCGTAGGGcc < 4:259124/133‑1 (MQ=255) ctcaAAGCTTTCGTTTGCTTTGCTGCCAAATGCGCCGCAGTCCATAGCTGGTCTGCCAACCACTCCGCTTCAGCAGGACTGAGCACCCACTCTACCCAGCCAAACTTAAGCATAATTGATTCGCCGTAGGGcc > 3:259124/1‑133 (MQ=255) tcaAAGCTTTCGTTTGCTTTGCTGCCAAATGCGCCGCAGTCCATAGCTGGTCTGCCAACCACTCCGCTTCAGCAGGACTGAGCACCCACTCTACCCAGCCAAACTTAAGCATAATTGATTCGCCGTAGGGCCTGATAGTCAAATCa > 2:311540/1‑146 (MQ=255) caAAGCTTTCGTTTGCTTTGCTGCCAAATGCGCCGCAGTCCATAGCTGGTCTGCCAACCACTCCGCTTCAGCAGGa > 6:146456/1‑76 (MQ=255) caAAGCTTTCGTTTGCTTTGCTGCCAAATGCGCCGCAGTCCATAGCTGGTCTGCCAACCACTCCGCTTCAGCAGGa < 5:146456/76‑1 (MQ=255) aaaGCTTTCGTTTGCTTTGCTGCCAAATGCGCCGCAGTCCATAGCTGGTCTGCCAACCACTCCGCTTCAGCAGGACTGAGCACCCACTCTACCCAGCCAAACTTAAGCATAATTGATTCGCCGTAGGGCCTGATAGTCAAATCATa > 7:44495/1‑146 (MQ=255) aaaGCTTTCGTTTGCTTTGCTGCCAAATGCGCCGCAGTCCATAGCTGGTCTGCCAACCACTCCGCTTCAGCAGGACTGAGCACCCACTCTACCCAGCCAAACTTAAGCATAATTGATTCGCCGTAGGGCCTGATAGTCAAATCATa < 8:44495/146‑1 (MQ=255) tttCGTTTGCTTTGCTGCCAAATGCGCCGCAGTCCATAGCTGGTCTGCCAACCACTCCGCTTCAGCAGGACTGAGCACCCACTCTACCCAGCCAAACTTAAGCATAATTGATTCGCCGTAGGGCCTGATAGTCAAATCATAGATGGccc > 2:218962/1‑149 (MQ=255) tttgctttgcTGCCAAATGCGCCGCAGTCCATAGCTGGTCTGCCAACCACTCCGCTTCAGCAGGACTGAGCACCCACTCTACCCAGCCAAACTTAAGCATAATTGATTCGCCGTAGGGCCTGATAGTCAAATCATAGATGGCCCCTGCt > 2:228747/1‑149 (MQ=255) tatgctttgctGCCAAATGCGCCGCAGTCCATAGCTGGTCTGCCAACCACTCCGCTTCAGCAGGACTGAGCACCCACTCTACCCAGCCAAACTTAAGCATAATTGATTCGCCGTAGGGCCTGATAGTCAAATCATAGATGGCCCCTGCt < 8:10530/147‑1 (MQ=255) gctttgctGCCAAATGCGCCGCAGTCCATAGCTGGTCTGCCAACCACTCCGCTTCAGCAGGACTGAGCACCCACTCTACCCAGCCAAACTTAAGCATAATTGATTCGCCGTAGGGCCTGATAGTCAAATCATAGATGGCCCCTGCTTc < 3:99796/148‑1 (MQ=255) ctttgctGCCTAATGCGCCGCAGTCCATAGCTTGTCTGCCAACCAGTCCGCTTCAGCAGGACTGAGCACCCACTCTACCCAGCCAAACTTAAGCATAATTGATTCGCCGTAGGGCCTGATAGTGAAATCATAGATGGCGCCTGCTTCg < 4:106477/148‑1 (MQ=255) tgcCAAATGCGCCGCAGTCCATAGCTGGTCTGCCAACCACTCCGCTTCAGCAGGACTGAGCACCCACTCTACCCAGCCAAACTTAAGCATAATTGATTCGCCGTAGGGCCTGATAGTCAAATCATAGATGGCCCCTGCTTCGTCCTTg > 4:158833/1‑148 (MQ=255) ccAAATGCGCCGCAGTCCATAGCTGGTCTGCCAACCACTCCGCTTCAGCAGGACTGAGCACCCACTCTACCCAGCCAAACTTAAGCATAATTGATTCGCCGTAGGGCCTGATAGTCAAATCATAGATGGCCCCTGCTTCGTCGTTGAAt > 8:234885/1‑149 (MQ=255) aaaTGCGCCGCAGTCCATAGCTGGTCTGCCAACCACTCCGCTTCAGCAGGACTGAGCACCCACTCTACCCAGCCAAACTTAAGCATAATTGATTCGCCGTAGGGCCTGATAGTCAAATCATAGATGGCCCCTGCTTCGTCCTTGAATcc > 4:139573/1‑149 (MQ=255) | GCGAGTCGGGTAGACCGCATAGATGCCCAGCTCTACCGAACACTACTCGGGCAATATCTCCACCAGATTTCTAAAGAAGCAAACCCATACTCTCGGCCCACCGTCCAAGGCTCAGCTCAAAGCTTTCGTTTGCTTTGCTGCCAAATGCGCCGCAGTCCATAGCTGGTCTGCCAACCACTCCGCTTCAGCAGGACTGAGCACCCACTCTACCCAGCCAAACTTAAGCATAATTGATTCGCCGTAGGGCCTGATAGTCAAATCATAGATGGCCCCTGCTTCGTCCTTGAATCC > NC_002947/5072170‑5072460 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |