Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 2,852,239 | T→C | 11.7% | intergenic (+107/+98) | PP_5531 → / ← PP_2505 | hypothetical protein/GAF domain/GGDEF domain‑containing protein |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 2,852,239 | 0 | T | C | 11.7% | 83.1 / 3.0 | 34 | intergenic (+107/+98) | PP_5531/PP_2505 | hypothetical protein/GAF domain/GGDEF domain‑containing protein |
Reads supporting (aligned to +/- strand): ref base T (10/20); new base C (2/2); total (12/22) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 6.02e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.08e-01 |
TTACGTGCAGCACATGTTCGCCTATCTGCGACAGTGAATCCTGCACGAGGTGGGCGCTAAGCAACTGTGCCGCGCTCGCCACTCGGGGCGGTCGCAGTTTGTTTGTCCAGCGATGTGGACGCCATTCCCGCCCCTGGTTGCAATAAGGGGCGGCGATCAGTTGAAAGCGAGCTGATCAGGCAAATTGGTAACTCAGCCGAATGTGATGCGTACCCTCCCGGATCCGGAGATCTTAACTTCCTCAGGATAACGCTGGCACATAGACAACTTGGTCACGTCCACCATGTTTGGC > NC_002947/2852096‑2852387 | ttACGTGCAGCACATGTTCGCCTATCTGCGACAGTGAATCCTGCACGAGGTGGGCGCTAAGCAACTGTGCCGCGCTCGCCACTCGGGGCGGTCGCAGTTTGTTTGTCCAGCGATGTGGACGCCATTCCCGCCCCTGGTTGCAATAAgggg > 1:252874/1‑150 (MQ=255) cGTGCAGCACATGTTCGCCTATCTGCGACAGTGAATCCTGCACGAGGTGGGCGCTAAGCAACTGTGCCGCGCTCGCCACTCGGGGCGGTCGCAGTTTGTTTGTCCAGCGATGTGGACGCCATTCCCGCCCCTGGTTGCAATAAGGggcgg < 1:16428/150‑1 (MQ=255) gTGCAGCACATGTTCGCCTATCTGCGACAGTGAATCCTGCACGAGGTGGGCGCTAAGCAACTGTGCCGCGCTCGCCCCTCGGGGCGGTCGCAGTTTGTTTGTCCAGCGATGTGGACGCCATTCCCGCCCCTGGTTGCAATAAGGggcggc < 2:252874/150‑1 (MQ=255) gcagcaCATGTTCGCCTATCTGCGACAGTGAATCCTGCACGAGGTGGGCGCTAAGCAACTGTGCCGCGCTCGGCACTCGGGGCGGTCGCAGTTTGTTTGTCCAGCGATGTGGACGCCATTCCCGCCCCTGGTTGCAATAAGGGGCGGCGa < 1:910408/150‑1 (MQ=255) agcaCATGTTCGCCTATCTGCGACAGTGAATCCTGCACGAGGTGGGCGCTAAGCAACTGTGCCGCGCTCGCCACTCGGGGCGGTCGCAGTTTGTTTGTCCAGCGATGTGGACGCCATTCCCGCCCCTTATTGCAACCAGGggcgggaatg < 2:416299/150‑6 (MQ=255) acaTGTTCGCCTATCTGCGACAGTGAATCCTGCACGAGGTGGGCGCTAAGCAACTGTGCCGCGCTCGCCACTCGGGGCGGTCGCAGTTTGTTTGTCCAGCGATGTGGACGCCATTCCCGCCCCTGGTTGCAATAAGGGGCGGCGATCAGt < 2:431101/150‑1 (MQ=255) acaTGTTCGCCTATCTGCGACAGTGAATCCTGCACGAGGGGGGCGCTAAGCAACTGTGCCGCGCTCGCCACTCGGGGCGGTCGCAGTTTGTTTGTCCAGCGATGTGGACGCCATTCCCGCCCCTGGTTGCAATAAGGCGCGGCGATCAGt < 2:431118/150‑1 (MQ=255) caTGTTCGCCTATCTGCGACAGTGAATCCTGCACGAGGTGGGCGCTAAGCAACTGTGCCGCGCTCGCCACTCGGGGCGGTCGCAGTTTGTTTGTCCAGCGATGTGGACGCCATTCCCGCCCCTGGTTGCAATAAGGGGCGGCGATCAGtt < 2:268575/150‑1 (MQ=255) ttCGCCTATCTGCGACAGTGAATCCTGCACGAGGTGGGCGCTAAGCAACTGTGCCGCGCTCGCCACTCGGGGCGGTCGCAGTTTGTTTGTCCAGCGATGTGGACGCCATTCCCGCCCCTGGTTGCAATAAGGGGCGGCGATCAGTTGaaa < 1:57956/150‑1 (MQ=255) cGCCTATCTGCGACAGTGAATCCTGCACGAGGTGGGCGCTAAGCAACTGTGCCGCGCTCGCCACTCGGGGCGGTCGCAGTTTGTTTGTCCAGCGATGTGGACGCCATTCCCGCCCCTGGTTGCAATAAGGGGCGGCGATCAGTTGAAAGc < 2:971313/150‑1 (MQ=255) ccTATCTGCGACAGTGAATCCTGCACGAGGTGGGCGCTAAGCAACTGTGCCGCGCTCGCCACTCGGGGCGGTCGCAGTTTGTTTGTCCAGCGATGTGGACGCCATTCCCGCCCCTGGTTGCAATAAGGGGCGGCGATCAGTTGAAAGCGa < 1:297868/150‑1 (MQ=255) gACAGTGAATCCTGCACGAGGTGGGCGCTAAGCAACTGTGCCGCGCTCGCCACTCGGGGCGGTCGCAGTTTGTTTGTCCAGCGATGTGGACGCCATTCCCGCCCCTGGTTGCAATAAGGGGCGGCGATCAg > 2:868451/1‑131 (MQ=255) gACAGTGAATCCTGCACGAGGTGGGCGCTAAGCAACTGTGCCGCGCTCGCCACTCGGGGCGGTCGCAGTTTGTTTGTCCAGCGATGTGGACGCCATTCCCGCCCCTGGTTGCAATAAGGGGCGGCGATCAg < 1:868451/131‑1 (MQ=255) gACAGTGAATCCTGCACGAGGTGGGCGCTAAGCAACTGTGCCGCGCTCGCCACTCGGGGCGGTCGCAGTTTGTTTGTCCAGCGATGTGGACGCCATTCCCGCCCCTGGTTGCAATAAGGGGCGGCGATCAGTTGAAAGCGAGCTGATCAg > 1:1030201/1‑150 (MQ=255) tGAATCCTGCACGAGGTGGGCGCTAAGCAACTGTGCCGCGCTCGCCACTCGGGGCGGTCGCAGTTTGTTTGTCCAGCGATGTGGACGCCATTCCCGCCCCTGGTTGCAATAAGGGGCGGCGATCAGTTGAAAGCGAGCTGATCAGGCaaa > 2:535293/1‑150 (MQ=255) tAAGCAACTGTGCCGCGCTCGCCACTCGGGGCGGTCGCAGTTTGTTTGTCCAGCGATGTGGACGCCATTCCCGCCCCTGGTTGCAATAAGGGGCGGCGATCAGTTGAAAGCGAGCTGATCAGGCAAATTGGTAACTCAGCCGAATGTGAt > 2:477884/1‑150 (MQ=255) gCCACTCGGGGCGGTCGCAGTTTGTTTGTCCAGCGATGTGGACGCCATTCCCGCCCCTGGTTGCAATAAGGGGCGGCGATCAGTTGAAAGCGAGCTGATCAGGCAAATTGGTAACTCAGCCGAATGTGATGCGTACCCTCCCGGATCCgg > 2:501932/1‑150 (MQ=255) gCGGTCGCAGTTTGTTTGTCCAGCGATGTGGACGCCATTCCCGCCCCTGGTTGCAATAAGGGGCGGCGATCAGTTGAAAGCGAGCTGATCAGGCAAATTGGTAACTCAGCCGAATGTGATGCGTACCCTCCCGGATCCGGAGATCTTAAc < 1:34610/150‑1 (MQ=255) aGTTTGTTTGTCCAGCGATGTGGACGCCATTCCCGCCCCTGGTTGCAATAAGGGGCGGCGATCAGTTGAAAGCGAGCTGATCAGGCAAAtt < 1:1037146/91‑1 (MQ=255) aGTTTGTTTGTCCAGCGATGTGGACGCCATTCCCGCCCCTGGTTGCAATAAGGGGCGGCGATCAGTTGAAAGCGAGCTGATCAGGCAAAtt > 2:1037146/1‑91 (MQ=255) gcttggCCAGCGATGTGGACGCCATTCCCGCCCCTGGTTGCCATAAGGGGCGGCGATCAGTTGAAAGCGAGCTGATCAGGCAAATTGGTAACTCAGCCGAATGTGATGCGTACCCTCCCGGATCCGGAGATCTTAACTTCCTCAGGATaa < 1:68759/148‑1 (MQ=255) gtCCAGCGATGTGGACGCCATTCCCGCCCCTGGTTGCAATAAGGGGCGGCGATCAGTTGAAAGCGAGCTGATCAGGCAAATTGGTAACTCAGCCGAATGTGATGCGTACCCTCCTGGATCCGGAGATCTTAACTTCCTCAGGATAACGCt < 2:1043556/150‑1 (MQ=255) ccAGCGATGTGGACGCCATTCCCGCCCCTGGTTGCAATAAGGGGCGGCGATCAGTTGAAAGCGAGCTGATCAGGCAAATTGGTAACTCAGCCGAATGTGATGCGTACCCTCCCGGATCCGGAGATCTTAACTTCCTCAGGATAACGCTgg < 1:501932/150‑1 (MQ=255) cGATGTGGACGCCATTCCCGCCCCTGGTTGCAATAAGGGGCGGCGATCAGTTGAAAGCGAGCTGATCAGGCAAATTGGTAACTCAGCCGAATGTGATGCGTACCCTCCCGGATCCGGAGATCTTAACTTCCTCAGGATAACGCTGGcaca < 1:535293/150‑1 (MQ=255) gACGCCATTCCCGCCCCTGGTTGCAATAAGGGGCGGCGATCAGTTGAAAGCGAGCTGATCAGGCAAATTGGTAACTCAGCCGAATGTGATGCGTACCCTCCCGGATCCGGAGATCTTAACTTCCTCAGGATAACGCTGGCACATAGACaa < 1:722183/150‑1 (MQ=255) gatcgccGCCCCTTATTGCAACCAGGGGCGGCGATAAGTTGAAAGCGAGCTGATCAGGCAAATTGGTAACTCAGCCGAATGTGATGCGTACCCTCCCGGATCCGGAGATCTTAACTTCCTCAGGATAACGCTGGCACATAGACAACTTgg > 2:419924/6‑150 (MQ=255) gccGCCCCTTATTGCAACCAGGGGCGGCGATCAGTTGAAAGCGAGCTGATCAGGCAAATTGGTAACTCAGCCGAATGTGATGCGTACCCTCCCGGATCCGGAGATCTTAACTTCCTCAGGATAACGCTGGCACATa > 2:589509/2‑136 (MQ=255) gccGCCCCTTATTGCAACCAGGGGCGGCGATCAGTTGAAAGCGAGCTGATCAGGCAAATTGGTAACTCAGCCGAATGTGATGCGTACCCTCCCGGATCCGGAGATCTTAACTTCCTCAGGATAACGCTGGCACATa < 1:589509/135‑1 (MQ=255) gTTGCAATAAGGGGCGGCGATCAGTTGAAAGCGAGCTGATCAGGCAAATTGGTAACTCAGCCGAATGTGATGCGTACCCTCCCGGATCCGGAGATCTTAACTTCCTCAGGATAACGCTGGCACATAGACAACTTGGTCACGTccacca < 1:495033/148‑1 (MQ=255) gTTGCAATAAGGGGCGGCGATCAGTTGAAAGCGAGCTGATCAGGCAAATTGGTAACTCAGCCGAATGTGATGCGTACCCTCCCGGATCCGGAGATCTTAACTTCCTCAGGATAACGCTGGCACATAGACAACTTGGTCACGTccacca > 2:495033/1‑148 (MQ=255) aaTAAGGGGCGGCGATCAGTTGAAAGCGAGCTGATCAGGCAAATTGGTAACTCAGCCGAAtgtg > 2:1024117/1‑64 (MQ=255) aaTAAGGGGCGGCGATCAGTTGAAAGCGAGCTGATCAGGCAAATTGGTAACTCAGCCGAAtgtg < 1:1024117/64‑1 (MQ=255) aaTAAGGGGCGGCGATCAGTTGAAAGCGAGCTGATCAGGCAAATTGGTAACTCAGCCGAATGTGATGCGTACCCTCCCGGATCCGGAGATCTTAACTTCCTCAGGATAACGCTGGCACATAGACAACTTGGTCACGTCCACCATGTTTgg > 1:152904/1‑150 (MQ=255) aTAAGGGGCGGCGATCAGTTGAAAGCGAGCTGATCAGGCAAATTGGTAACTCAGCCGAATGTGATGCGTACCCTCCCGGATCCGGAGATCTTAACTTCCTCAGGATAACGCTGGCACATAGACAACTTGGTCACGTCCACCATGTTTGGc < 2:152904/150‑1 (MQ=255) | TTACGTGCAGCACATGTTCGCCTATCTGCGACAGTGAATCCTGCACGAGGTGGGCGCTAAGCAACTGTGCCGCGCTCGCCACTCGGGGCGGTCGCAGTTTGTTTGTCCAGCGATGTGGACGCCATTCCCGCCCCTGGTTGCAATAAGGGGCGGCGATCAGTTGAAAGCGAGCTGATCAGGCAAATTGGTAACTCAGCCGAATGTGATGCGTACCCTCCCGGATCCGGAGATCTTAACTTCCTCAGGATAACGCTGGCACATAGACAACTTGGTCACGTCCACCATGTTTGGC > NC_002947/2852096‑2852387 |
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: ‑ |