Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 4,159,110 | C→T | 13.3% | intergenic (+14/‑195) | PP_3661 → / → PP_3662 | membrane protein/hypothetical protein |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 4,159,110 | 0 | C | T | 13.3% | 75.3 / 3.8 | 30 | intergenic (+14/‑195) | PP_3661/PP_3662 | membrane protein/hypothetical protein |
Reads supporting (aligned to +/- strand): ref base C (10/16); new base T (2/2); total (12/18) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.12e-01 |
CATGGAGACCCGCTTCAGCCAGATTCGCCAGGCTGGGCCACGGGCGCTGGCCACCGGTGCGATCCTCAACCTGTGGCTGGTTGGCGGGGGGCTGGCGATTACCCTGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTCGCGGGTAAACCCGCTCCTACAGGTGGGTGCTGCCCTGAAGACGGCAAGGCCCCTGTAGGAGCGGGTTCACCCGCGAAGCAGGCAGCTCGGTCTCAAGGTTTTTCACAGCCCCCCTCGTTGCCGTCCACATTTCT > NC_002947/4158967‑4159253 | cATGGAGACCCGCTTCAGCCAGATTCGCCAGGCTGGGCCACGGGAGATAGCGACCTGTGCGATCCTCAACCTGTGGCTGGTTGGCGGGGGGCTGGCGATTACCCTGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGgctgc < 1:31253/150‑1 (MQ=255) ccGCTTCAGCCAGATTCGCCAGGCTGGGCCACGGGCGCTGGCCACCGGTGCGATCCTCAACCTGTGGCTGGTTGGCGGGGGGCTGGCGATTACCCTGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTCGCgg < 1:274731/148‑1 (MQ=255) ccGCTGCAGCCAGATTCGCCAGGCTGGGCCACGGGCGCTGGCCACCGGTGCGATCCTCAACCTGTGGCTGGTTGGCGGGGGGCTGGCGATTACCCTGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTCGCgg > 2:274731/1‑148 (MQ=255) tCAGCCAGATTCGCCAGGCTGGGCCACGGGCGCTGGCCACCGGTGCGATCCTCAACCTGTGGCTGGTTGGCGGGGGGCTGGCGATTACCCTGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTCGCGGGTAAAcc < 2:163880/150‑1 (MQ=255) aGCCAGATTCGCCAGGCTGGGCCACGGGCGCTGGCCACCGGTGCGATCCTCAACCTGTGGCTGGTTGGCGGGGGGCTGGCGATTACCCTGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTCGCGGGTAAACCCg > 1:405229/1‑150 (MQ=255) gCCAGATTCGCCAGGCTGGGCCACGGGCGCTGGCCACCGGTGCGATCCTCAACCTGTGGCTGGTTGGCGGGGGGCTGGCGATTACCCTGGGGGTGCAGAAGCTGTTGGGCTGGCTGGAGACCGAGTCGGCTGCTTCGCGGGTAAACCCGc < 1:459788/150‑1 (MQ=255) tCGCCAGGCTGGGCCACGGGCGCTGGCCACCGGTGCGATCCTCAACCTGTGGCTGGTTGGCGGGGGGCTGGCGATTACCCTGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTCGCGGGTAAACCCGCTCCTACa < 2:19294/150‑1 (MQ=255) gCCAGGCTGGGCCACGGGCGCTGGCCACCGGTGCGATCCTCAACCTGTGGCTGGTTGGCGGGGGGCTGGCGATTACCCTGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTCGCGGGTAAACCCGCTCCTACAgg > 2:239984/1‑150 (MQ=255) gCTGGGCCACGGGCGCTGGCCACCGGTGCGATCCTCAACCTGTGGCTGGTTGGCGGGGGGCTGGCGATTACCCTGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTCGCGGGTAAACCCGATCCTACAggtgggt < 2:85685/150‑1 (MQ=255) cTGGCCACCGGTGCGATCCTCAACCTGTGGCTGGTTGGCGGGGGGCTGGCGATTACCCTGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTCGCGGGTAAACCCGCTCCTACAGGTGGGTGCTGCCCTGAAGACg < 1:464937/150‑1 (MQ=255) gATCCTCAACCTGTGGCTGGTTGGCGGGGGGCTGGCGATTACCCTGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTCGCGGGTAAACCCGCTCCTACAGGTGGGTGCTGCCCTGAAGACGGCAAGGccc > 2:41378/1‑145 (MQ=255) gATCCTCAACCTGTGGCTGGTTGGCGGGGGGCTGGCGATTACCCTGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTCGCGGGTAAACCCGCTCCTACAGGTGGGTGCTGCCCTGAAGACGGCAAGGccc < 1:41378/145‑1 (MQ=255) ccTGTGGCTGGTTGGCGGGGGGCTGGCGATTACCCTGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTCGCGGGTAAACCCGCTCCTACAGGTGGGTGCTGCCCTGAAGACGGCAAGGCCCCTGTAGGa < 2:111298/144‑1 (MQ=255) ccTGTGGCTGGTTGGCGGGGGGCTGGCGATTACCCTGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTCGCGGGTAAACCCGCTCCTACAGGTGGGTGCTGCCCTGAAGACGGCAAGGCCCCTGTAGGa > 1:111298/1‑144 (MQ=255) gCTGGTTGGCGGGGGGCTGGCGATTACCCTGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTCGCGGGTAAACCCGCTCCTACAGGTGGGTGCTGCCCTGAAGACGGCAAGGCCCCTGTAGGAGCGGGTTCAccc > 2:330951/1‑150 (MQ=255) tGGTTGGCGGGGGGCTGGAGATTAACATCGGGGTGCAGAAGCTGTTTTGCTGACTGGAGACCGAGTCGGCTGCTTCGCGGGTAAACCCGCTCCTACAGGTGGGTGCTGCCCTGAAGACGGCAAGGCCCCTGTAGGAGCGGGTTCACCcgc < 2:387123/150‑1 (MQ=38) gCGATTCCCCTGGGGGTGCAGAAGCTGTTGTGCTGATAGGAGAGCGAGTCGGCTGTTTCGCGGGTGAACCCGCTCCTACAGGTGGGTGCTGCCCTCAAGACGGCAAGGTCCCTGTAGGAGCGGGTTCACCCGCGAAGCAGGCAGCTCGGt < 2:163103/150‑1 (MQ=255) aTTACCCTGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTCGCGGGTAAACCCGCTCCTACAGGTGGGTGCTGCCCTGAAGACGGCAAGGCCCCTGTAGGAGCggg > 2:161134/1‑121 (MQ=255) aTTACCCTGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTCGCGGGTAAACCCGCTCCTACAGGTGGGTGCTGCCCTGAAGACGGCAAGGCCCCTGTAGGAGCggg < 1:161134/121‑1 (MQ=255) ttACCCTGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTCGCGGGTAAACCCGCTCCTACAGGTGGGTGCTGCCCTGAAGACGGCAAGGCCCCTGTAGGAGCGGGTTCACCCGCGAAGCAGGCAGCTCGGTCTCa < 2:508081/150‑1 (MQ=255) gAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTCGCGGGTAAACCCGCTCCTACAGGTGGGTGCTGCCCTGAAGACGGCAAGGCCCCTGTAGGAGCGGGTTCACCCGCGAAGCAGGCAGCTCGGTCTCAAGGTTTTTCACAGccc > 1:158675/1‑150 (MQ=255) tGGGCTGACTGGAGACCGAGTCGGCTGCTTCGCGGGTAAACCCGCTCCTACAGGTGGGTGCTGCCCTGAAGACGGCAAGGCCCCTGTAGGAGCGGGTTCACCCGCGAAGCAGGCAGCTCGGTCTCAAGGTTTTTCACAGCCCCCCTCGtt < 2:391819/150‑1 (MQ=255) gtcTGACTGGAGACCGAGTCGGCTGCTTCGCGGGTAAACCAGCTCATACCGGCGGGTGCTGCCCTGAAGACGTCAAGGCCACTGTCGAAACGGATTCATTCTTTAAGCAGGCAGCTCGGTCTCAAGGTTTTTCACAGCCCCCCTCGTTga > 1:391819/3‑149 (MQ=255) ttgagaCCGAGCTGCCTGCTTCGCGGGTGAACCCGCTCCTACAGGTGGGCGCTGCCCTGAAGACGGCAAGGCCCCTGTAGGa > 2:191310/3‑82 (MQ=17) ttgagaCCGAGCTGCCTGCTTCGCGGGTGAACCCGCTCCTACAGGTGGGCGCTGCCCTGAAGACGGCAAGGCCCCTGTAGGa > 2:191311/3‑82 (MQ=17) ttgagaCCGAGCTGCCTGCTTCGCGGGTGAACCCGCTCCTACAGGTGGGCGCTGCCCTGAAGACGGCAAGGCCCCTGTAGGa < 1:191311/80‑1 (MQ=18) ttgagaCCGAGCTGCCTGCTTCGCGGGTGAACCCGCTCCTACAGGTGGGCGCTGCCCTGAAGACGGCAAGGCCCCTGTAGGa < 1:191310/80‑1 (MQ=17) ccGAGTCGGCTGCTTCGCGGGTAAACCCGCTCCTACAGGTGGGTGCTGCCCTGAAGACGGCAAGGCCCCTGTAGGAGCGGGTTCACCCGCGAAGCAGGCAGCTCGGTCTCAAGGTTTTTCACAGCCCCCCTCGTTGCCGTCCACATTtct < 2:405229/150‑1 (MQ=255) aGTCGGCTGCTTCGCGGGTAAACCCGCTCCTACAGGTGGGTGCTGCCCTGAAGACGGCAAGGCCCCTGTAGGAGCGGGTTCACCCGCGAAGCAGGCAGCTCGGTCTCAAGGTTTTTCACAGCCCCCCTCGTTg < 1:76112/133‑1 (MQ=255) aGTCGGCTGCTTCGCGGGTAAACCCGCTCCTACAGGTGGGTGCTGCCCTGAAGACGGCAAGGCCCCTGTAGGAGCGGGTTCACCCGCGAAGCAGGCAGCTCGGTCTCAAGGTTTTTCACAGCCCCCCTCGTTg > 2:76112/1‑133 (MQ=255) | CATGGAGACCCGCTTCAGCCAGATTCGCCAGGCTGGGCCACGGGCGCTGGCCACCGGTGCGATCCTCAACCTGTGGCTGGTTGGCGGGGGGCTGGCGATTACCCTGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTCGCGGGTAAACCCGCTCCTACAGGTGGGTGCTGCCCTGAAGACGGCAAGGCCCCTGTAGGAGCGGGTTCACCCGCGAAGCAGGCAGCTCGGTCTCAAGGTTTTTCACAGCCCCCCTCGTTGCCGTCCACATTTCT > NC_002947/4158967‑4159253 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |