Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 1,547,033 | T→A | 9.3% | W124R (TGG→AGG) | PP_1357 → | hypothetical protein |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 1,547,033 | 0 | T | A | 9.3% | 112.3 / 3.6 | 43 | W124R (TGG→AGG) | PP_1357 | hypothetical protein |
Reads supporting (aligned to +/- strand): ref base T (17/22); new base A (2/2); total (19/24) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.93e-01 |
GGACGCCTACAGCGGAGAGAAAGCCAAGTCGCTCAAGGCGGCAGAGGCAGGCAAGGTGACCTTGGCCAATGGCGATTCACCGAAGACGGTGGCGCTGGCCAGCATCGAGCAGATCATGAAGCCCAAGCCGCTGGTCGAGGATTTCGTGTGGAAAGGTAATGTCGACGTGGCGCTGGACTACAAGCGCGCCGAGAACGACACCGACGATTATGACGTGGGCTTCAAGACCACTGCCCGCCATGGCCGCTGGCGGCACAACGCCGAGGGTGAGTACAACCGCGAGACCAAGGACGACG > NC_002947/1546885‑1547180 | ggACGCCTACAGCGGAGAGAAAGCCAAGTCGCTCAAGGCGGCAGAGGCAGGCAAGGTGACCTTGGCCAATGGCGATTCACCGAAGACGGTGGCGCTGGCCAGCATCGAGCAGATCATGAAGCCCAAGCCGCTGGTCGAGGATTTCgtgt < 8:132039/149‑1 (MQ=255) gACGCCTACAGCGGAGAGAAAGCCAAGTCGCTCAAGGCGGCAGAGGCAGGCAAGGTGACCTTGGCCAATGGCGATTCACCGAAGACGGTGGCGCTGGCCAGCATCGAGCAGATCATGAAGCCCAAGCCGCTGGTCGAGGATTTCGtgtg > 2:41223/1‑149 (MQ=255) cGCCTACAGCGGAGAGAAAGCCAAGTCGCTCAAGGCGGCAGAGGCAGGCAAGATGACCTTGGCCAATGGCGAATCACCGAAGACGGTGGCGCTGGCCAGCATCGAGCAGATCATTAAGCCCAAGCCGCTGGTCGAGGATTTCgtgttg > 8:297392/1‑146 (MQ=255) cTACAGCGGAGAGAAAGCCAAGTCGCTCAAGGCGGCAGAGGCAGGCAAGGTGACCTTGGCCAATGGCGATTCACCGAAGACGGTGGCGCTGGCCAGCATCGAGCAGATCATGAAGCCCAAGCCGCTGGTCGAGGATTTCGTGTGGAAAg < 8:102715/149‑1 (MQ=255) gCCAAGTCGCTCAAGGCGGCAGAGGCAGGCAAGGTGACCTTGGCCAATGGCGATTCACCGAAGACGGTGGCGCTGGCCAGCATCGAGCAGATCATGAAGCCCAAGCCGCTGGTCGAGGATTTCGTGTGGAAAGGTAATGTCGACGTGgc > 1:84082/1‑149 (MQ=255) aaGTCGCTCAAGGCGGCAGAGGCAGGCAAGGTGACCTTGGCCAATGGCGATTCACCGAAGACGGTGGCGCTGGCCAGCATCGAGCAGATCATGAAGCCCAAGCCGCTGGTCGAGGATTTCGTGTGGAAAGGTAATGTCGACGTGGCGCt < 5:99717/149‑1 (MQ=255) aaGTCGCTCAAGGCGGCAGAGGCAGGCAAGGTGACCTTGGCCAATGGCGATTCACCGAAGACGGTGGCGCTGGCCAGCATCGAGCAGATCATGAAGCCCAAGCCGCTGGTCGAGGATTTCGTGTGGAAAGGTAATGTCGACGTGGCGCt < 5:174723/149‑1 (MQ=255) aaGTCGCTCAAGGCGGCAGAGGCAGGCAAGGTGACCTTGGCCAATGGCGATTCACCGAAGACGGTGGCGCTGGCCAGCATCGAGCAGATCATGAAGCCCAAGCCGCTGGTCGAGGATTTCGTGTGGAAAGGTAATGTCGACGTGGCGCt < 4:184131/149‑1 (MQ=255) cTCAAGGCGGCAGAGGCAGGCAAGGTGACCTTGGCCAATGGCGATTCACCGAAGACGGTGGCGCTGGCCAGCATCGAGCAGATCATGAAGCCCAAGCCGCTGGTCGAGGATTTCGTGTGGAAAGGTAATGTCGACGTGGCGCTGGACTa > 2:297955/1‑149 (MQ=255) aaGGCGGCAGAGGCAGGCAAGGTGACCTTGGCCAATGGCGATTCACCGAAGACGGTGGCGCTGGCCAGCATCGAGCAGATCATGAAGCCCAAGCCGCTGGTCGAGGATTTCGTGTGGAAAGGTAATGTCGACGTGGCGCTGGACTAc > 3:275713/1‑147 (MQ=255) aaGGCGGCAGAGGCAGGCAAGGTGACCTTGGCCAATGGCGATTCACCGAAGACGGTGGCGCTGGCCAGCATCGAGCAGATCATGAAGCCCAAGCCGCTGGTCGAGGATTTCGTGTGGAAAGGTAATGTCGACGTGGCGCTGGACTACaa < 2:221534/149‑1 (MQ=255) aaGGCGGCAGAGGCAGGCAAGGTGACCTTGGCCAATGGCGATTCACCGAAGACGGTGGCGCTGGCCAGCATCGAGCAGATCATGAAGCCCAAGCCGCTGGTCGAGGATTTCGTGTGGAAAGGTAATGTCGACGTGGCGCTGGACTACaa < 2:216045/149‑1 (MQ=255) ggcggcAGAGGCAGGCAAGGTGACCTTGGCCAATGGCGATTCACCGAAGACGGTGGCGCTGGCCAGCATCGAGCAGATCATGAAGCCCAAGCCGCTGGTCGAGGATTTCGTGAGGAAAGGTAATGTCGACGTGGCGCTGGACTACAAgc > 3:364409/1‑149 (MQ=255) ggcggcAGAGGCAGGCAAGGTGACCTTGGCCAATGGCGATTCACCGAAGACGGTGGCGCTGGCCAGCATCGAGCAGATCATGAAGCCCAAGCCGCTGGTCGAGGATTTCGTGAGGAAAGGTAATGTCGACGTGGCGCTGGACTACAAgc > 3:358524/1‑149 (MQ=255) gcggcAGAGGCAGGCAAGGTGACCTTGGCCAATGGCGATTCACCGAAGACGGTGGCGCTGGCCAGCATCGAGCAGATCATGAAGCCCAAGCCGCTGGTCGAGGATTTCGTGTGGAAAGGTAATGTCGACGTGGCGCTGGACTACAAgcg > 8:531/1‑149 (MQ=255) ggcAGAGGCAGGCAAGGTGACCTTGGCCAATGGCGATTCACCGAAGACGGTGGCGCTGGCCAGCATCGAGCAGATCATGAAGCCCAAGCCGCTGGTCGAGGATTTCGTGTGGAAAGGTAATGTCGACGTGGCGCt < 2:360105/135‑1 (MQ=255) ggcAGAGGCAGGCAAGGTGACCTTGGCCAATGGCGATTCACCGAAGACGGTGGCGCTGGCCAGCATCGAGCAGATCATGAAGCCCAAGCCGCTGGTCGAGGATTTCGTGTGGAAAGGTAATGTCGACGTGGCGCt > 1:360105/1‑135 (MQ=255) cAGAGGCAGGCAAGGTGACCTTGGCCAATGGCGATTCACCGAAGACGGTGGCGCTGGCCAGCATCGAGCAGATCATGAAGCCCAAGCCGCTGGTCGAGGATTTCGTGAGGAAAGGTAATGTCGACGTGGCGCTGGACTACAAGCGCGcc < 4:364409/149‑1 (MQ=255) cAGAGGCAGGCAAGGTGACCTTGGCCAATGGCGATTCACCGAAGACGGTGGCGCTGGCCAGCATCGAGCAGATCATGAAGCCCAAGCCGCTGGTCGAGGATTTCGTGAGGAAAGGTAATGTCGACGTGGCGCTGGACTACAAGCGCGcc < 4:358524/149‑1 (MQ=255) agagGCAGGCAAGGTGACCTTGGCCAATGGCGATTCACCGAAGACGGTGGCGCTGGCCAGCATCGAGCAGATCATGAAGCCCAAGCCGCTGGTCGAGGATTTCGTGTGGAAAGGTAATGTCGACGTGGCGCTGGACTACAAGCGCGcc > 4:188118/1‑148 (MQ=255) gcaggcaAGGTGACCTTGGCCAATGGCGATTCACCGAAGACGGTGGCGCTGGCCAGCATCGAGCAGATCATGAAGCCCAAGCCGCTGGTCGAGGATTTCGTGTGGAAAGGTAATGTCGACGTGGCGCTGGACTACAAGCGCGCCGAGaa > 7:96717/1‑149 (MQ=255) caggcaAGGTGACCTTGGCCAATGGCGATTCACCGAAGACGGTGGCGCTGGCCAGCATCGAGCAGATCATGAAGCCCAAGCCGCTGGTCGAGGATTTCGTGTGGAAAGGTAATGTCGACGTGGCGCTGGACTACAAGCGCGCCGAGAAc < 3:228870/149‑1 (MQ=255) aggcaAGGTGACCTTGGCCAATGGCGATTCACCGAAGACGGTGGCGCTGGCCAGCATCGAGCAGATCATGAAGCCCAAGCCGCTGGTCGAGGATTTCGTGTGGAAAGGTAATGTCGACGTGGCGCTGGACTACAAGCGCGCCGAGAACg < 1:41223/149‑1 (MQ=255) caAGGTGACCTTGGCCAATGGCGATTCACCGAAGACGGTGGCGCTGGCCAGCATCGAGCAGATCATGAAGCCCAAGCCGCTGGTCGAGGATTTCGTGTGGAAAGGTAATGTCGACGTGGCGCTGGACTACAAGCGCGCCGAGAACGaca < 1:59923/149‑1 (MQ=255) aaGGTGACCTTGGCCAATGGCGATTCACCGAAGACGGTGGCGCTGGCCAGCATCGAGCAGATCATGAAGCCCAAGCCGCTGGTCGAGGATTTCGTGTGGAAAGGTAATGTCGACGTGGCGCTGGACTACAAGCGCGCCGAGAACGacac > 5:144210/1‑149 (MQ=255) cTTGGCCAATGGCGATTCACCGAAGACGGTGGCGCTGGCCAGCATCGAGCAGATCATGAAGCCCAAGCCGCTGGTCGAGGATTTCGTGTGGAAAGGTAATGTCGACGTGGCGCTGGACTACAAGCGCGCCGAGAACGACACCGACGAtt < 3:188118/149‑1 (MQ=255) cTTGGCCAATGGCGATTCACCGAAGACGGTGGCGCTGGCCAGCATCGAGCAGATCATGAAGCCCAAGCCGCTGGTCGAGGATTTCGTGTGGAAAGGTAATGTCGACGTGGCGCTGGACTACAAGCGCGCCGAGAACGACACCGACGAt > 1:114271/1‑148 (MQ=255) aTGGCGATTCACCGAAGACGGTGGCGCTGGCCAGCATCGAGCAGATCATGAAGCCCAAGCCGCTGGTCGAGGATTTCGTGTGGAAAGGTAATGTCGACGTGGCGCTGGACTACAAGCGCGCCGAGAACGACACCGACGATTATGACGTg < 7:531/149‑1 (MQ=255) aTGGCGATTCACCGAAGACGGTGGCGCTGGCCAGCATCGAGCAGATCATGAAGCCCAAGCCGCTGGTCGAGGATTTCGTGTGGAAAGGTAATGTCGACGTGGCGCTGGACTACAAGCGCGCCGAGAACGACACCGACGATTATGACGTg < 4:369697/149‑1 (MQ=255) cGATTCACCGAAGACGGTGGCGCTGGCCAGCATCGAGCAGATCATGAAGCCCAAGCCGCTGGTCGAGGATTTCGTGTGGAAAGGTAATGTCGACGTGGCGCTGGACTACAAGCGCGCCgag < 7:223014/121‑1 (MQ=255) cGATTCACCGAAGACGGTGGCGCTGGCCAGCATCGAGCAGATCATGAAGCCCAAGCCGCTGGTCGAGGATTTCGTGTGGAAAGGTAATGTCGACGTGGCGCTGGACTACAAGCGCGCCgag > 8:223014/1‑121 (MQ=255) aCCGAAGACGGTGGCGCTGGCCAGCATCGAGCAGATCATGAAGCCCAAGCCGCTGGTCGAGGATTTCGTGTGGAAAGGTAATGTCGACGTGGCGCTGGACTACAAGCGCGCCGAGAACGACACCGACGATTATGACGTGGGCTTCaa > 6:274403/1‑147 (MQ=255) aCCGAAGACGGTGGCGCTGGCCAGCATCGAGCAGATCATGAAGCCCAAGCCGCTGGTCGAGGATTTCGTGTGGAAAGGTAATGTCGACGTGGCGCTGGACTACAAGCGCGCCGAGAACGACACCGACGATTATGACGTGGGCTTCaa < 5:274403/147‑1 (MQ=255) aCCGAAGACGGTGGCGCTGGCCAGCATCGAGCAGATCATGAAGCCCAAGCCGCTGGTCGAGGATTTCGTGTGGAAAGGTAATGTCGACGTGGCGCTGGACTACAAGCGCGCCGAGAACGACACCGACGATTATGACGTGGGCTTCAAGa < 5:208466/149‑1 (MQ=255) ggTGGCGCTGGCCACCATCGAGCAGATCATGAAGCCCAAGCCGCTGGTCGAGGATTTCGTGTGGAAAGGTAATGTCGACGTGGCGCTGGACTACCAGCGCGGCGAGAACGACACCGACGGTTATGACGTGGGCTTCAAGCCCACTGccc > 8:160378/1‑149 (MQ=255) gcTGGCCAGCATAGAGCAGATCATGAAGCCCAAGCCGCTGGTCGAGGAAATCGTGTGGAAAGGTAATCTCGACCTGCCGCTGGACTACAAGCGCGCCGAGAACGACACAGACGATTAAGACGTGGGCTTCAAGACCACTGCCCGCCATg > 5:86090/1‑149 (MQ=255) cAGCATCGAGCAGATCATGAAGCCCAAGCCGCTGGTCGAGGATTTCGTGTGGAAAGGTAATGTCGACGTGGCGCTGGACTACAAGCGCGCCGAGAACGACACCGACGATTATGACGTGGGCTTCAAGACCACTGCCCGCCATGGCCGCt > 1:1862/1‑149 (MQ=255) cATCGAGCAGATCATGAAGCCCAAGCCGCTGGTCGAGGATTTCGTGTGGAAAGGTAATGTCGACGTGGCGCTGGACTACAAGCGCGCCGAGAACGACACCGACGATTATGACGTGGGCTTCAAGACCACTGCCCGCCATGGCCGCTggc < 8:334844/149‑1 (MQ=255) aGCAGATCATGAAGCCCAAGCCGCTGGTCGAGGATTTCGTGTGGAAAGGTAATGTCGACGTGGCGCTGGACTACAAGCGCGCCGAGAACGACACCGACGATTATGACGTGGGCTTCAAGACCACTGCCCGCCATGGCCGCTGGCGGca > 4:175447/1‑148 (MQ=255) tGAAGCCCAAGCCGCTGGTCGAGGATTTCGTGTGGAAAGGTAATGTCGACGTGGCGCTGGACTACAAGCGCGCCGAGAACGACACcga > 4:268461/1‑88 (MQ=255) tGAAGCCCAAGCCGCTGGTCGAGGATTTCGTGTGGAAAGGTAATGTCGACGTGGCGCTGGACTACAAGCGCGCCGAGAACGACACcga < 3:268461/88‑1 (MQ=255) aGCCCAAGCCGCTGGTCGAGGATTTCGTGTGGAAAGGTAATGTCGACGTGGCGCTGGACTACAAGCGCGCCGAGAACGACACCGACGATTATGACGTGGGCTTCAAGACCACTGCCCGCCATGGCCGCTGGCGGCACAACGCCGAGGGt < 2:114271/149‑1 (MQ=255) cAAGCCGCTGGTCGAGGATTTCGTGTGGAAAGGTAATGTCGACGTGGCGCTGGACTACAAGCGCGCCGAGAACGACACCGACGATTATGACGTGGGCTTCAAGACCACTGCCCGCCATGGCCGCTGGCGGCACAACGCCGAGGGTGAGt < 8:96717/149‑1 (MQ=255) gCTGGTCGAGGACCTCCTGTGGAAAGGTAAGGCCTCCCTGGCCCCGGACAACAAGCGCGCCGCCAACCACACCGACGATTATGACGCTGGCTTCAAGTCCACTGCCCGCCATGGCCGCTCGCGGCACAACGCCGAGGGTGAGTACAccc < 4:275713/149‑1 (MQ=255) tCGAGGATTTCGTGTGGAAAGGTAATGTCGACGTGGCGCTGGACTACAAGCGCGCCGAGAACGACACCGACGATTATGACGTGGGCTTCAAGACCACTGCCCGCCATGGCCGCTGGCGGCACAACGCCGAGGGTGAGTACAACcgc > 7:232982/1‑146 (MQ=255) aTTTCGTGTGGAAAGGTAATGTCGACGTGGCGCTGGACTACAAGCGCGCCGAGAACGACACCGACGATTATGACGTGGGCTTCAAGACCACTGCCCGCCATGGCCGCTGGCGGCACAACGCCGAGGGTGAGTACAACCGCGAGACCAAg < 7:297392/149‑1 (MQ=255) tgtgGAAAGGTAATGTCGACGTGGCGCTGGACTACAAGCGCGCCGAGAACGACACCGACGATTATGACGTGGGCTTCAAGACCACTGCCCGCCATGGCCGCTGGCGGCACAACGCCGAGGGTGAGTACAACCGCGAGACCAAGgacgac > 4:149752/1‑149 (MQ=255) tgtgGAAAGGTAATGTCGACGTGGCGCTGGACTACAAGCGCGCCGAGAACGACACCGACGATTATGACGTGGGCTTCAAGACCACTGCCCGCCATGGCCGCTGGCGGCACAACGCCGAGGGTGAGTACAACCGCGAGACCAAGgacgac < 3:249697/149‑1 (MQ=255) gtgGAAAGGTAATGTCGACGTGGCGCTGGACTACAAGCGCGCCGAGAACGACACCGACGATTATGACGTGGGCTTCAAGACCACTGCCCGCCATGGCCGCTGGCGGCACAACGCCGAGGGTGAGTACAACCGCGAGACCAAGGacgacg > 2:353832/1‑149 (MQ=255) tgGAAAGGTAATGTCGACGTGGCGCTGGACTACAAGCGCGCCGAGAACGACACCGACGATTATGACGTGGGCTTCAAGACCACTGCCCGCCATGGCCGCTGGCGGCACAAc > 2:368779/1‑111 (MQ=255) tgGAAAGGTAATGTCGACGTGGCGCTGGACTACAAGCGCGCCGAGAACGACACCGACGATTATGACGTGGGCTTCAAGACCACTGCCCGCCATGGCCGCTGGCGGCACAAc < 1:368779/111‑1 (MQ=255) | GGACGCCTACAGCGGAGAGAAAGCCAAGTCGCTCAAGGCGGCAGAGGCAGGCAAGGTGACCTTGGCCAATGGCGATTCACCGAAGACGGTGGCGCTGGCCAGCATCGAGCAGATCATGAAGCCCAAGCCGCTGGTCGAGGATTTCGTGTGGAAAGGTAATGTCGACGTGGCGCTGGACTACAAGCGCGCCGAGAACGACACCGACGATTATGACGTGGGCTTCAAGACCACTGCCCGCCATGGCCGCTGGCGGCACAACGCCGAGGGTGAGTACAACCGCGAGACCAAGGACGACG > NC_002947/1546885‑1547180 |
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: ‑ |