Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 2,852,440 | T→A | 15.8% | Q312L (CAG→CTG) | PP_2505 ← | 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,440 | 0 | T | A | 15.8% | 94.0 / 11.2 | 38 | Q312L (CAG→CTG) | PP_2505 | GAF domain/GGDEF domain‑containing protein |
Reads supporting (aligned to +/- strand): ref base T (17/15); new base A (3/3); total (20/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 = 7.01e-01 |
TACCCTCCCGGATCCGGAGATCTTAACTTCCTCAGGATAACGCTGGCACATAGACAACTTGGTCACGTCCACCATGTTTGGCCGCATAGAGACCGCCGTCAGCGCGCTTTATTGCATCATCAATCTGCTCACCCTGGGTGAGCAGGGCGCACCCAATGCTCATGGTCATCTTCCGAGGGCCGCCGGGTAACAGAAACGGCGCGCTGCCGATTGACTGACGCAAGCGCTCCGCAATCTGCGCCAATTCTTGGCTATCTTTAACCCGTAGCAAGGCTACGAATTC > NC_002947/2852306‑2852588 | tACCCTCCCGGATCCGGAGATCTTAACTTCCTCAGGATAACGCTGGCACATAGACAACTTGGTCACGTCCACCATGTTTGGCCGCATAGAGACCGCCGTCAGCGCGCTTTATTGCATCATCAATCTGCTCACCCTGGGTGAGCAGGgcgc < 2:1030201/150‑1 (MQ=255) ccGGATCCGGAGATCTTAACTTCCTCAGGATAACGCTGGCACATAGACAACTTGGTCACGTCCACCATGTTTGGCCGCATAGAGACCGCCGTCAGCGCGCTTTATTGCATCATCAATCTGCTCACCCTGGGTGAGCAGGGCGCACCCAAt > 1:343835/1‑150 (MQ=255) cGGATCCGGAGATCTTAACTTCCTCAGGATAACGCTGGCACATAGACAACTTGGTCACGTCCACCATGTTTGGCCGCATAGAGACCGCCGTCAGCGCGCTTTATTGCATCATCAATCTGCTCACCCTGGGTGAGCAGGGCGCACCCAATg > 2:1120034/1‑150 (MQ=255) cGGATCCGGAGATCTTAACTTCCTCAGGATAACGCTGGCACATAGACAACTTGGTCACGTCCACCATGTTTGGCCGCATAGAGACCGCCGTCAGCGCGCTTTATTGCATCATCAATCTGCTCACCCTGGGTGAGCAGGGCGCACCCAATg < 1:477884/150‑1 (MQ=255) ggAGATCTTAACTTCCTCAGGATAACGCTGGCACATAGACAACTTGGTCACGTCCACCATGTTTGGCCGCATAGAGACCGCCGTCAGCGCGCTTTATTGCATCATCAATCTGCTCACCCTGGGTGAGCAGGGCGCACCCAATGCTCATgg > 2:447914/1‑150 (MQ=255) gagaTCTTAACTTCCTCAGGATAACGCTGGCACATAGACAACTTGGTCACGTCCACCATGTTTGGCCGCATAGAGACCGCCGTCAGCGCGCTTTATTGCATCATCAATCTGCTCACCCTGGGTGAGCAGGGCGCACCCAATGCTCATGGt > 2:341208/1‑150 (MQ=255) agaTCTTAACTTCCTCAGGATAACGCTGGCACATAGACAACTTGGTCACGTCCACCATGTTTGGCCGCATAGAGACCGCCGTCAGCGCGCTTTATTGCATCATCAATCTGCTCACCCTGGGTGAGCAGGGCGCACCCAATGCTCATGGTc < 1:419924/150‑1 (MQ=255) ttCCTCAGGATAACGCTGGCACATAGACAACTTGGTCACGTCCACCATGTTTGGCCGCATAGAGACCGCCGTCAGCGCGCTTTATTGCATCATCAATCTGCTCACCCTGGGTGAGCAGGGCGCACCCAATGCTCATGGTCATCTTCCGAg > 1:1017726/1‑150 (MQ=255) cTCAGGATAACGCTGGCACATAGACAACTTGGTCACGTCCACCATGTTTGGCCGCATAGAGACCGCCGTCAGCGCGCTTTATTGCATCATCAATCTGCTCACCCTGGGTGAGCAGGGCGCACCCAATGCTCATGGTCATCTTCCGAGGgc > 1:129930/1‑150 (MQ=255) tgaCGCTGGCACATAGACAACTTGGTCACGTCCACCATGTTTGGCCGCATAGAGACCGCCGTCAGCGCGCTTTATGGCATCATCAATCTGCTCACCCTGGGTGAGCAGGGCGCACCCAATGCTCATGGTCATCTTCCGAGGGCCGCCggg < 1:447914/148‑1 (MQ=255) cTGGCACATAGACAACTTGGTCACGTCCACCATGTTTGGCCGCATAGAGACCGCCGTCAGCGCGCTTTATTGCATCATCAATCTGCTCACCCTGGGTGAGCAGGGCGCACCCAATGCTCATGGTCATCTTCCGAGGGCCGCCGGGTAACa < 1:792659/150‑1 (MQ=255) cTTGGTCACGTCCACCATGTTTGGCCGCATAGAGACCGCCGTCAGCGCGCTTTATTGCATCATCAATCTGCTCACCCTgg < 1:1111625/80‑1 (MQ=255) cTTGGTCACGTCCACCATGTTTGGCCGCATAGAGACCGCCGTCAGCGCGCTTTATTGCATCATCAATCTGCTCACCCTgg > 2:1111625/1‑80 (MQ=255) ggTCACGTCCACCATGTTTGGCCGCATAGAGACCGCCGTCAGCGCGCTTTATTGCATCATCAATCTGCTCACCCTGGGTGAGCAGGGCGCACCCAATGCTCATGGTCATCTTCCGAGGGCCGCCGGGTAACAGAAACGGCGCGCTGCCGa < 2:1017726/150‑1 (MQ=255) caccaTGTTTGGCCGCATAGAGACCGCCGTCAGCGCGCTTTATTGCATCATCAACCTGCTCACCCAGGGTGAGCAga > 1:119676/1‑76 (MQ=255) caccaTGTTTGGCCGCATAGAGACCGCCGTCAGCGCGCTTTATTGCATCATCAACCTGCTCACCCAGGGTGAGCAga < 2:119676/77‑2 (MQ=39) cGCATAGAGACCGCCGTCAGCGCGCTTTATTGCATCATCAATCTGCTCACCCTGGGTGAGCAGGGCGCACCCAATGCTCATGGTCATCTTCCGAGGGCCGCCGGGTAACAGAAACGGCGCGCTGCCGATTGACTGACGCAAGCGCTCCGc > 2:75235/1‑150 (MQ=255) cGCATAGAGACCGCCGTCAGCGCGCTTTATTGCATCATCAATCTGCTCACCCTGGGTGAGCAGGGCGCACCCAATGCTAATGGTCATCTTCCGAGGGCCGCCGGGTAACAGAAAc > 1:129608/1‑115 (MQ=255) cGCATAGAGACCGCCGTCAGCGCGCTTTATTGCATCATCAATCTGCTCACCCTGGGTGAGCAGGGCGCACCCAATGCTAATGGTCATCTTCCGAGGGCCGCCGGGTAACAGAAAc < 2:129608/115‑1 (MQ=255) aTAGAGACCGCCGTCAGCGCGCTTTATTGCATCATCAATCTGCTCACCCTGGGTGAGCAGGGCGCACCCAATGCTCATGGTCATCTTCCGAGGGCCGCCGGGTAACAGAAACGGCGCGCTGCCGATTGACTGACGCAAGCGCTCCGCAAt < 2:109338/150‑1 (MQ=255) aTAGAGACCGCCGTCAGCGCGCTTTATTGCATCATCAATCTGCTCACCCAGGGTGAGCAga > 1:217222/1‑60 (MQ=255) aTAGAGACCGCCGTCAGCGCGCTTTATTGCATCATCAATCTGCTCACCCAGGGTGAGCAga < 2:217222/61‑2 (MQ=255) gcgcgcTTTATTGCATCATCAATCTGCTCACCCTGGGTGAGCAGGGCGCACCCAATGCTCATGGTCATCTTCCGAGGGCCGCCGGGTAACAGAAACGGCGCGCTGCCGATTGACTGACGCAAGCGCTCCGCAATCTGCGCCAATTCTTgg < 1:1120034/150‑1 (MQ=255) gCATCATCAATCTGCTCACCCTGGGTGAGCAGGGCGCACCCAATGCTCATGGTCATCTTCCGAGGGCCGCCGGGTAACAGAAACGGCGCGCTGCCGATTGACTGACGCAAGCGCTCCGCAATCTGCGCCAATTCTTGGCTATCTTTAAcc > 2:543026/1‑150 (MQ=255) gCATCATCAATCTGCTCACCCTGGGTGAGCAGGGCGCACCCAATGCTCATGGTCATCTTCCGAGGGCCGCCGGGTAACAGAAACGGCGCGCTGCCGATTGACTGACGCAAGCGCTCCGCAATCTGCGCCAATTCTTGGCTATCTTTAAcc < 1:341208/150‑1 (MQ=255) catcaATCTGCTCACCCTGGGTGAGCAGGGCGCACCCAATGCTCATGGTCATCTTCCGAGGGCCGCCGGGTAACAGAAACGGCGCGCTGCCGATTGACTGACGCAAGCGCTCCGCAATCTGCGCCAATTCTTGGCTATCTTTAACCCGTa > 2:631963/1‑150 (MQ=255) atcaATCTGCTCACCCTGGGTGAGCAGGGCGCACCCAATGCTCATGGTCATCTTCCGAGGGCCGCCGGGt < 1:847853/70‑1 (MQ=255) atcaATCTGCTCACCCTGGGTGAGCAGGGCGCACCCAATGCTCATGGTCATCTTCCGAGGGCCGCCGGGt > 2:847853/1‑70 (MQ=255) atcaATCTGCTCACCCTGGGTGAGCAGGGCGCACCCAATGCTCATGGTCATCTTCCGAGGGCCGCCGGGTAACAGAAACGGCGCGCTGCCGATTGACTGACGCAAGCGCTCCGCAATCTGCGCCAATTCTTGGCTATCTTTAACCCGTAg > 1:68426/1‑150 (MQ=255) ccTGCTCACCCAGGGTGAGCAGGGCGCACCCAATGCTCATGGTCATCTTCCGAGGGCCGCCGGGTAACAGAAACGGCGCGCTGCCGATTGACTGACGCAAGCGCTCCGCAATCTgc > 1:994649/2‑116 (MQ=255) ccTGCTCACCCAGGGTGAGCAGGGCGCACCCAATGCTCATGGTCATCTTCCGAGGGCCGCCGGGTAACAGAAACGGCGCGCTGCCGATTGACTGACGCAAGCGCTCCGCAATCTgc < 2:994649/115‑1 (MQ=255) gCTCACCCTGGGTGAGCAGGGCGCACCCAATGCTCATGGTCATCTTCCGAGGGCCGCCGGGTAACAGAAACGGCGCGCTGCCGATTGACTGACGCAAg < 1:115620/98‑1 (MQ=255) gCTCACCCTGGGTGAGCAGGGCGCACCCAATGCTCATGGTCATCTTCCGAGGGCCGCCGGGTAACAGAAACGGCGCGCTGCCGATTGACTGACGCAAg > 2:115620/1‑98 (MQ=255) gCTCACCCTGGGTGAGCAGGGCGCACCCAATGCTCATGGTCATCTTCCGAGGGCCGCCGGGTAACAGAAACGGCGCGCTGCCGATTGACTGACGCAAGCGCTCCGCAATCTGCGCCAATTCTTGGCTATCTTTAACCCGTAGCAAGGCTa > 2:188261/1‑150 (MQ=255) aCCCTGGGTGAGCAGGGCGCACCCAATGCTCATGGTCATCTTCCGAGGGCCGCCGGGTAACAGAAACGGCGCGCTGCCGATTGACTGACGCAAGCGCTCCGCAATCTGCGCCAATTCTTGGCTATCTTTAACCCGTAGCAAGGCTACGaa > 1:95601/1‑150 (MQ=255) ccTGGGTGAGCAGGGCGCACCCAATGCTCATGGTCATCTTCCGAGGGCCGCCGGGTAACAGAAACGGCGCGCTGCCGATTGACTGACGCAAGCGCTCCGCAATCTGCGCCAATTCTTGGCTATCTTTAACCCGTAGCAAGGCTACGAAtt > 1:995404/1‑150 (MQ=255) cTGGGTGAGCAGGGCGCACCCAATGCTCATGGTCATCTTCCGAGGGCCGCCGGGTAACAGAAACGGCGCGCTGCCGATTGACTGACGCAAGCGCTCCGCAATCTGCGCCAATTCTTGGCTATCTTTAACCCGTAGCAAGGCTACGAATtc < 2:343835/150‑1 (MQ=255) cTGGGTGAGCAGGGCGCACCCAATGCTCATGGTCATCTTCCGAGGGCCGCCGGGTAACAGAAACGGCGCGCTGCCGATTGACTGACGCAAGCGCTCCGCAATCTGCGCCAATTCTTGGCTATCTTTAACCCGTAGCAAGGCTACGAATtc < 2:995404/150‑1 (MQ=255) | TACCCTCCCGGATCCGGAGATCTTAACTTCCTCAGGATAACGCTGGCACATAGACAACTTGGTCACGTCCACCATGTTTGGCCGCATAGAGACCGCCGTCAGCGCGCTTTATTGCATCATCAATCTGCTCACCCTGGGTGAGCAGGGCGCACCCAATGCTCATGGTCATCTTCCGAGGGCCGCCGGGTAACAGAAACGGCGCGCTGCCGATTGACTGACGCAAGCGCTCCGCAATCTGCGCCAATTCTTGGCTATCTTTAACCCGTAGCAAGGCTACGAATTC > NC_002947/2852306‑2852588 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |