| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | NC_002947 | 2,852,232 | G→A | 13.8% | intergenic (+100/+105) | 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,232 | 0 | G | A | 13.8% | 65.8 / 4.3 | 29 | intergenic (+100/+105) | PP_5531/PP_2505 | hypothetical protein/GAF domain/GGDEF domain‑containing protein |
| Reads supporting (aligned to +/- strand): ref base G (7/18); new base A (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 = 6.72e-01 | |||||||||||
GGAAAGGCCGATTACGTGCAGCACATGTTCGCCTATCTGCGACAGTGAATCCTGCACGAGGTGGGCGCTAAGCAACTGTGCCGCGCTCGCCACTCGGGGCGGTCGCAGTTTGTTTGTCCAGCGATGTGGACGCCATTCCCGCCCCTGGTTGCAATAAGGGGCGGCGATCAGTTGAAAGCGAGCTGATCAGGCAAATTGGTAACTCAGCCGAATGTGATGCGTACCCTCCCGGATCCGGAGATCTTAACTTCCTCAGGATAACGCTGGCACATAGACAACTTGGTCACGTCCACCA > NC_002947/2852085‑2852379 | 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) gTTGCAATAAGGGGCGGCGATCAGTTGAAAGCGAGCTGATCAGGCAAATTGGTAACTCAGCCGAATGTGATGCGTACCCTCCCGGATCCGGAGATCTTAACTTCCTCAGGATAACGCTGGCACATAGACAACTTGGTCACGTccacca > 2:495033/1‑148 (MQ=255) gTTGCAATAAGGGGCGGCGATCAGTTGAAAGCGAGCTGATCAGGCAAATTGGTAACTCAGCCGAATGTGATGCGTACCCTCCCGGATCCGGAGATCTTAACTTCCTCAGGATAACGCTGGCACATAGACAACTTGGTCACGTccacca < 1:495033/148‑1 (MQ=255) | GGAAAGGCCGATTACGTGCAGCACATGTTCGCCTATCTGCGACAGTGAATCCTGCACGAGGTGGGCGCTAAGCAACTGTGCCGCGCTCGCCACTCGGGGCGGTCGCAGTTTGTTTGTCCAGCGATGTGGACGCCATTCCCGCCCCTGGTTGCAATAAGGGGCGGCGATCAGTTGAAAGCGAGCTGATCAGGCAAATTGGTAACTCAGCCGAATGTGATGCGTACCCTCCCGGATCCGGAGATCTTAACTTCCTCAGGATAACGCTGGCACATAGACAACTTGGTCACGTCCACCA > NC_002947/2852085‑2852379 |
| 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: ‑ |