Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 2,852,231 | G→T | 13.8% | intergenic (+99/+106) | 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,231 | 0 | G | T | 13.8% | 73.6 / 3.4 | 29 | intergenic (+99/+106) | PP_5531/PP_2505 | hypothetical protein/GAF domain/GGDEF domain‑containing protein |
Reads supporting (aligned to +/- strand): ref base G (7/18); new base T (2/2); total (9/20) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 5.68e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.38e-01 |
CTCGGAAAGGCCGATTACGTGCAGCACATGTTCGCCTATCTGCGACAGTGAATCCTGCACGAGGTGGGCGCTAAGCAACTGTGCCGCGCTCGCCACTCGGGGCGGTCGCAGTTTGTTTGTCCAGCGATGTGGACGCCATTCCCGCCCCTGGTTGCAATAAGGGGCGGCGATCAGTTGAAAGCGAGCTGATCAGGCAAATTGGTAACTCAGCCGAATGTGATGCGTACCCTCCCGGATCCGGAGATCTTAACTTCCTCAGGATAACGCTGGCACATAGACAACTTGG > NC_002947/2852082‑2852367 | cTCGGAAAGGCCGATTACGTGCAGCACATGTTCGCCTATCTGCGACAGTGAATCCTGCACGAGGTGGGCGCTAAGCAACTGTGCCGCGCTCGCCACTCGGGGCGGTCGCAGTTTGTTTGTCCAGCGATGTGGACGCCATTCCCGCCCCTg > 1:268575/1‑150 (MQ=255) ggAAAGGCCGATTACGTGCAGCACATGTTCGCCTATCTGCGACAGTGAATCCTGCACGAGGTGGGCGCGAAGCAACTGTGCCGCGCTCGCCACTCGGGGCGGTCGCAGTTTGTTTGTCCAGCGATGTGGACGCCATTCCCGCCCCTGGtt < 2:723944/150‑1 (MQ=255) 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 < 1:868451/131‑1 (MQ=255) gACAGTGAATCCTGCACGAGGTGGGCGCTAAGCAACTGTGCCGCGCTCGCCACTCGGGGCGGTCGCAGTTTGTTTGTCCAGCGATGTGGACGCCATTCCCGCCCCTGGTTGCAATAAGGGGCGGCGATCAg > 2:868451/1‑131 (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 > 2:1037146/1‑91 (MQ=255) aGTTTGTTTGTCCAGCGATGTGGACGCCATTCCCGCCCCTGGTTGCAATAAGGGGCGGCGATCAGTTGAAAGCGAGCTGATCAGGCAAAtt < 1:1037146/91‑1 (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 < 1:589509/135‑1 (MQ=255) gccGCCCCTTATTGCAACCAGGGGCGGCGATCAGTTGAAAGCGAGCTGATCAGGCAAATTGGTAACTCAGCCGAATGTGATGCGTACCCTCCCGGATCCGGAGATCTTAACTTCCTCAGGATAACGCTGGCACATa > 2:589509/2‑136 (MQ=255) | CTCGGAAAGGCCGATTACGTGCAGCACATGTTCGCCTATCTGCGACAGTGAATCCTGCACGAGGTGGGCGCTAAGCAACTGTGCCGCGCTCGCCACTCGGGGCGGTCGCAGTTTGTTTGTCCAGCGATGTGGACGCCATTCCCGCCCCTGGTTGCAATAAGGGGCGGCGATCAGTTGAAAGCGAGCTGATCAGGCAAATTGGTAACTCAGCCGAATGTGATGCGTACCCTCCCGGATCCGGAGATCTTAACTTCCTCAGGATAACGCTGGCACATAGACAACTTGG > NC_002947/2852082‑2852367 |
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: ‑ |