Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 1,691,825 | C→A | 28.2% | intergenic (+50/‑121) | PP_1487 → / → PP_1488 | hypothetical protein/methyl‑accepting chemotaxis transducer |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 1,691,825 | 0 | C | A | 28.2% | 53.1 / 17.3 | 35 | intergenic (+50/‑121) | PP_1487/PP_1488 | hypothetical protein/methyl‑accepting chemotaxis transducer |
Reads supporting (aligned to +/- strand): ref base C (12/13); new base A (4/6); total (16/19) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 7.23e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.50e-01 |
GGAGCGTGCCTACCTGGACGAGCTGGCCCGCCAGCTGCGGCTGGACCCGGGCTTGCGCCAGGAACTGGAAAGCCAGGTGCGGTTGGCTGCAGCGTAGTGGCGGTGCTTGCCACGGCCTATTGAGCGCCTGTGAGATCGAGCGCCGCCCGTGCGGCGCTCCGTCTCACAGGCGCAAAATTCCTGACGACAATCACATCAATACCGTAAGACCTCTGGGCTATACTTCGGCGATTTTTCCCGCTCGTTGTGAATTCCTGAGGGCTGAATGTGAAGAACTGGACCTTGCGCCAACGGA > NC_002947/1691679‑1691973 | ggAGCGTGCCTACCTGGACGAGCTGGCCCGCCAGCTGCGGCTGGACCCGGGCTTGCGCCAGGAACTGGAAAGCCAGGTGCGGTTGGCTGCAGCGTAGTGGCGGTGCTTGCCACGGCCTATTGAGCGCCTGTGAGACGGAGCGCCGCACgg > 2:476565/1‑149 (MQ=255) aGCGTGCCTACCTGGACGAGCTGGCCCGCCAGCTGCGGCTGGACCCGGGCTTGCGCCAGGAACTGGAAAGCCAGGTGCGGTTGGCTGCAGCGTAGTGGCGGTGCTTGCCACGGCCTATTGAGCGCCTGTGAGATCGAGCGCCGCCCGTgc < 2:220710/150‑1 (MQ=255) gCGTGCCTACCTGGACGAGCTGGCCCGCCAGCTGCGGCTGGACCCGGGCTTGCGCCAGGAACTGGAAAGCCAGGTGCGGTTGGCTGCAGCGTAGTGGCGGTGCTTGCCACGGCCTATTGAGCGCCTGTGAGACGGAGCGCCGCACGGgcg > 2:684652/1‑150 (MQ=255) gCCTACCTGGACGAGCTGGCCCGCCAGCTGCGGCTGGACCCGGGCTTGCGCCAGGAACTGGAAAGCCAGGTGCGGTTGGCTGCAGCGTAGTGGCGGTGCTTGCCACGGCCTATTGAGCGCCTGTGAGATCGAGCGCCGCCCGTGCGgcgc < 2:197928/150‑1 (MQ=255) gtACCTGGACGAGCTGGCCCGCCAGCTGCGGCTGGACCCGGGCTTGCGCCAGGAACTGGAAAGCCAGGTGCGGTTGGCTGCAGCGTAGTGGCGGTGCTTGCCACGGCCGATTGAGCGCCTGTGAGATCGAGCGCCGCCCGTGCGGCGCTc > 2:545734/2‑150 (MQ=255) cTACCTGGACGAGCTGGCCCGCCAGCTGCGGCTGGACCCGGGCTTGCGCCAGGAACTGGAAAGCCAGGTGCGGTTGGCTGCAGCGTAGTGGCGGTGCTTGCCACGGCCTATTGAGCGCCTGTGAGATCGAGCGCCGCCCGTGCGGCGCTc > 2:269291/1‑150 (MQ=255) ccTGGACGAGCTGGCCCGCCAGCTGCGGCTGGACCCGGGCTTGCGCCAGGAACTGGAAAGCCAGGTGCGGTTGGCTGCAGCGTAGTGGCGGTGCTTGCCACGGCCTATTGAGCGCCTGTGAGACGGAGCGCCGCACGGGCGGCGCTcgat < 1:476565/150‑4 (MQ=39) cTGGACGAGCTGGCCCGCCAGCTGCGGCTGGACCCGGGCTTGCGCCAGGAACTGGAAAGCCAGGTGCGGTTGGCTGCAGCGTAGTGGCGGTGCTTGCCACGGCCTATTGAGCGCCTGTGAGACGGAGCGCCGCACGGGCGGCGCTcgatc < 1:684652/150‑5 (MQ=38) gagaTGGCCTGCCAGCTGCGGCTGGACCCGGGCTTGCGCCAGGAACTGGAAAGCCAGGTGCGGTTGGCTGCAGCGTAGTGGCGGTGCTTGCCACGGCCTATTCAGCGCCTGTGAGATCGAGCGCCGCCCGTGCGGCGCTCCGTCTCACAg < 2:43668/150‑1 (MQ=255) gAGCTGGCCCGCCAGCTGCGGCTGGACCCGGGCTTGCGCCAGGAACTGGAAAGCCAGGTGCGGTTGGCTGCAGCGTAGTGGCGGTGCTTGCCACGGCCTATTGAGCGCCTGTGAGACGGAGCGCCGCACGGGCGGCGCTCGATCTCACAg < 1:333955/150‑1 (MQ=38) gAGCTGGCCCGCCAGCTGCGGATGGACCTGCGCTTGCTCCAGGAAGGGGAAGGCCAGGTGCGGTGGGCTGTAGCGTAGTGGCGGTGCTTGGCACGGCCTATTGAGCGCCTGTGAGATCGAGCGCCGCCCGTGCGGCGCTCCGTCTCACAg < 2:372563/150‑1 (MQ=38) aGCTGGCCCGCCAGCTGCGGCTGGACCCGGGCTTGCGCCAGGAACTGGAAAGCCAGGTGCGGTTGGCTGCAGCGTAGTGGCGGTGCTTGCCACGGCCTATTGAGCGCCTGTGAGACGGAGCGCCGCACGGGCGGCGCTCGATCTCACAgg < 1:21994/150‑1 (MQ=37) tGGCCCGCCAGCTGCGGCTGGACCCGGGCTTGCGCCAGGAACTGGAAAGCCAGGTGCGGTTGGCTGCAGCGTAGTGGCGGTGCTTGCCACGGCCTATTGAGCGCCTGTGAGATAGAGCGCCGCCCGTGCGGCGCTCCGTCTCACAGgcgc > 1:205207/1‑150 (MQ=255) cAGCTGCGGCTGGACCCGGGCTTGCGCCAGGAACTGGAAAGCCAGGTGCGGTTGGCTGCAGCGTAGTGGCGGTGCTTGCCACGGCCTATTGAGCGCCTGTGAGATCGAGCGCCGCCCGTGCGGCGCTCCGTCTCACAGGCGCAAAATTcc > 1:183120/1‑150 (MQ=255) aGCTGCGGCTGGACCCGGGCTTGCGCCAGGAACTGGAAAGCCAGGTGCGGTTGGCTGCAGCGTAGTGGCGGTGCTTGCCACGGCCTATTGAGCGCCTGTGAGATCGAGCGCCGCCCGTGCGGCGCTCCGTCTCACAGGCGCAAAATTCCt < 2:296594/150‑1 (MQ=255) aGCTGCGGCTGGACCCGGGCTTGCGCCAGGAACTGGAAAGCCAGGTGCGGTTGGCTGCAGCGTAGTGGCGGTGCTTGCCACGGCCTATTGAGCGCCTGTGAGATCGAGCGCCGCCCGTGCGGCGCTCCGGCTCACAGGCGCAAAATTCCt < 1:269291/150‑1 (MQ=255) gACCCGGGCTTGCGCCAGGAACTGGAAAGCCAGGTGCGGTTGGCTGCAGCGTAGTGGCGGTGCTTGCCACGGCCTATTGAGCGCCTGTGAGATCGAGCGCCGCCCGTGCGGCGCTCCGTCTca < 2:195395/123‑1 (MQ=255) gACCCGGGCTTGCGCCAGGAACTGGAAAGCCAGGTGCGGTTGGCTGCAGCGTAGTGGCGGTGCTTGCCACGGCCTATTGAGCGCCTGTGAGATCGAGCGCCGCCCGTGCGGCGCTCCGTCTca > 1:195395/1‑123 (MQ=255) cccGGGCTTGCGCCAGGAACTGGAAAGCCAGGTGCGGTTGGCTGCAGCGTAGTGGCGGTGCTTGCCACGGCCTATTGAGCGCCTGTGAGATCGAGCGCCGCCCGTGCGGCGCTCCGTCTCACAGGCGCAAAATTCCTGACGACAATcaca < 2:183120/150‑1 (MQ=255) cccGGGCTTGCGCCAGGAACTGGAAAGCCAGGTGCGGTTGGCTGCAGCGTAGTGGCGGTGCTTGCCACGGCCTATTGAGCGCCTGTGAGATCGAGCGCCGCCCGTGCGGCGCTCCGTCTCACAGGCGCAAAATTCCTGACGACAATcaca < 2:235372/150‑1 (MQ=255) tGGAAAGCCAGGTGCGGTTGGCTGCAGCGTAGTGGCGGTGCTTGCCACGGCCTATTGAGCGCCTGTGAGATCGAGCGCCGCCCGTGCGGCGCTCCGTCTCACAGGCGCAAAATTCCTGACGACAATCACATCAATACCGTAAGACCTCTg < 2:436164/150‑1 (MQ=255) tGGAAAGCCAGGTGCGGTTGGCTGCAGCGTAGTGGCGGTGCTTGCCACGGCCTATTGAGCGCCTGTGAGACGGAGCGCCGCACGGGCGGCGCTCGATCTCACAg > 1:672991/1‑104 (MQ=21) tGGAAAGCCAGGTGCGGTTGGCTGCAGCGTAGTGGCGGTGCTTGCCACGGCCTATTGAGCGCCTGTGAGACGGAGCGCCGCACGGGCGGCGCTCGATCTCACAg < 2:672991/104‑1 (MQ=21) ccAGGTGCGGTTGGCTGCAGCGTAGTGGCGGTGCTTGCCACGGCCTATTGAGCGCCTGTGAGATCGAGCGCCGCCCGTGCGGCGCTCCGTCTCACAGGCGCAAAATTCCTGACGACAATCACATCAATACCGTAAGACCTCTGGGCtata > 1:480869/1‑150 (MQ=255) ggCTGCAGCGTAGTGGCGTTGCTTGCCACGGCCTATTGAGCGCCTGTGAGATCGAGCGCCGCCCGTGCGGCGCTCCGTCTCACAGGCGCAAAATTCCTGACGACAATCACATCAATACCGTAAGACCTCTGGGCTATACTTCGGCGAttt < 1:213778/150‑1 (MQ=255) tGCAGCGTAGTGGCGGTGCTTGCCACGGCCTATTGAGCGCCTGTGAGATCAAGCGCCGCCCGTGCGGCGCTCCGTCACACAGGCGCAAAATTCCTGACGACAATCACATCAATACCGTAAGACCTCTTGGCTATACTTCGGCGAtttttt > 2:609047/1‑149 (MQ=255) gCAGCGTAGTGGCGGTGCTTGCCACGGCCTATTGAGCGCCTGTGAGATCGAGCGCCGCCCGTGCGGCCCTCCGTCTCACAGGCGCAAAATTCCTGACGACAATCACATCAATACCGTAAGACATCTGGGCTATACTTCGGCGATTTTTcc > 2:352343/1‑150 (MQ=255) ccTATTGAGCGCCTGTGAGATCGAGCGCCGCCCGTGCGGCGCTCCGTCTCACAGGCGCAAAATTCCTGACGACAATcac < 1:169158/79‑1 (MQ=255) ccTATTGAGCGCCTGTGAGATCGAGCGCCGCCCGTGCGGCGCTCCGTCTCACAGGCGCAAAATTCCTGACGACAATcac > 2:169158/1‑79 (MQ=255) ccTATTGAGCGCCTGTGAGATCGAGCGCCGCCCGTGCGGCGCTCCGTCTCACAGGCGCAAAATTCCTGACGACAATCACATCAATACCGTAAGACCTCTGGGCTATACTTCGGCGATTTTTCCCGCTCGTTGTGAATTCCTGAGGGCTGa > 1:793227/1‑150 (MQ=255) ccTGTGAGATCGAGCGCCGCCCGTGCGGCGCTCCGTCTCACAGGCGCAAAATTCCTGACGACAATcac < 2:193752/68‑1 (MQ=255) ccTGTGAGATCGAGCGCCGCCCGTGCGGCGCTCCGTCTCACAGGCGCAAAATTCCTGACGACAATcac > 1:193752/1‑68 (MQ=255) gagaTCGAGCGCCGCCCGTGCGGCGCTCCGTCTCACAGGCGCAAAATTCCTGACGACAATCACATCAATACCGTAAGACCTCTGGGCTATACTTCGGCGATTTTTCCCGCTCGTTGTGAATTCCTGAGGGCTGAATGTGAAGAACTGGAc > 1:472269/1‑150 (MQ=255) gacggAGCGCCGCACGGGCGGCGCTCGATCTCACAGGCGCAAAATTCCTGACGACAATCACATCAATACCGTAAGACCTCTGGGCTATACTTCGGCGATTTTTCCCGCTCGTTGTGAATTCCTGAGGGCTGAATGTGAAGAACt < 2:163874/140‑1 (MQ=38) gacggAGCGCCGCACGGGCGGCGCTCGATCTCACAGGCGCAAAATTCCTGACGACAATCACATCAATACCGTAAGACCTCTGGGCTATACTTCGGCGATTTTTCCCGCTCGTTGTGAATTCCTGAGGGCTGAATGTGAAGAACt > 1:163874/5‑144 (MQ=38) gccCGTGCGGCGCTCCGTCTCACAGGCGCAAAATTCCTGACGACAATCACATCACTAAAGTAAGACCTCTGGGCGAGAGGTCGGCGATTTTTCCCGCTCGTTGTGAATTCCTGAGGGCTGAATGTGAAGAACTGGACCTTGCGCCAACgg > 1:695973/1‑150 (MQ=255) cccGTGCGGCGCTCCGTCTCACAGGCGCAAAATTCCTGACGACAATCACATCAATACCGTAAGACCTCTGGGCTATACTTCGGCGATTTTTCCCGCTCGTTGTGAATTCCTGAGGGCTGAATGTGAAGAACTGGACCTTGCGCCAACGGa > 1:302860/1‑150 (MQ=255) | GGAGCGTGCCTACCTGGACGAGCTGGCCCGCCAGCTGCGGCTGGACCCGGGCTTGCGCCAGGAACTGGAAAGCCAGGTGCGGTTGGCTGCAGCGTAGTGGCGGTGCTTGCCACGGCCTATTGAGCGCCTGTGAGATCGAGCGCCGCCCGTGCGGCGCTCCGTCTCACAGGCGCAAAATTCCTGACGACAATCACATCAATACCGTAAGACCTCTGGGCTATACTTCGGCGATTTTTCCCGCTCGTTGTGAATTCCTGAGGGCTGAATGTGAAGAACTGGACCTTGCGCCAACGGA > NC_002947/1691679‑1691973 |
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: ‑ |