Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 5,705,255 | G→A | 12.1% | C243Y (TGC→TAC) | phaA → | poly(3‑hydroxyalkanoate) polymerase |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 5,705,255 | 0 | G | A | 12.1% | 89.1 / 4.0 | 33 | C243Y (TGC→TAC) | phaA | poly(3‑hydroxyalkanoate) polymerase |
Reads supporting (aligned to +/- strand): ref base G (18/11); new base A (2/2); total (20/13) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.64e-01 |
CGATGTGCTGGAACTGATCCAGTACAGCCCCATCACCGAGCAGGTGCATGCCCGTCCGCTGCTGGTGGTGCCACCGCAGATCAACAAGTTCTACGTGTTCGACCTCAGCCCGGAAAAGAGCCTGGCGCGCTTCTGCCTGCGCTCGCAGCAGCAGACCTTCATCATCAGCTGGCGCAACCCGACCAAGGCCCAGCGTGAATGGGGCCTGTCCACCTACATCGATGCGCTGAAAGAAGCCGTCGACGCGGTGCTGTCGATTACCGGCA > NC_002947/5705121‑5705386 | cGATGTGCTGGAACTGATCCAGTACAGCCCCATCACCGAGCAGGTGCATGCCCGTCCGCTGCTGGTGGTGCCACCGCAGATCAACAAGTTCTACGTGTTCGACCTCAGCCCGGAAAAGAGCCTGGCGCGCTTctgc > 1:385410/1‑136 (MQ=255) tGGAACTGATCCAGTACAGCCCCATCACCGAGCAGGTGCATGCCCGTCCGCTGCTGGTGGTGCCACCGCAGATCAACAAGTTCTACGAGTTCGACCTCAGCCCGGAAAAGAGCCTGGCGCGCTTCTGCCTGCGCTc > 2:260186/1‑136 (MQ=255) tGGAACTGATCCAGTACAGCCCCATCACCGAGCAGGTGCATGCCCGTCCGCTGCTGGTGGTGCCACCGCAGATCAACAAGATCTACGTGATCGACCTCAGCCCGGAAAAGAGCCTGGCGCGCTTCTGCCTGCGCTc > 2:110320/1‑136 (MQ=255) ccAGTACAGCCCCATCACCGAGCAGGTGCATGCCCGTCCGCTGCTGGTGGTGCCACCGCAGATCAACAAGTTCTACGTGTTCGACCTCAGCCCGGAAAAGAGCCTGGCGCGCTTCTGCCTGCGCTCgcagcagcac > 1:51826/1‑135 (MQ=255) aGTACAGCCCCATCACCGAGCAGGTGCATGCCCGTCCGCTGCTGGTGGTGCCACCGCAGATCAACAAGTTCTACGTGTTCGACCTCAGCCCGGAAAAGAGCCTGGCGCGCTTCTACCTGCGCTCGCAGCAGCAGAc > 2:112025/1‑136 (MQ=255) aGTACAGCCCCATCACCGAGCAGGTGCATGCCCGTCCGCTGCTGGTGGTGCCACCGCAGATCAACAAGTTCTACGTGTTCGACCTCAGCCCGGAAAAGAGCCTGGCGCGCTTCTACCTGCGCTCGCAGCAGCAGAc > 2:109957/1‑136 (MQ=255) cAGCCCCATCACCGAGCAGGTGCATGCCCGTCCGCTGCTGGTGGTGCCACCGCAGATCAACAAGTTCTACGTGTTCGACCTCAGCCCGGAAAAGAGCCTGGCGCGCTTCTGCCTGCGCTCGCAGCAGCAGACCTtc > 2:75972/1‑136 (MQ=255) cAGCCCCATCACCGAGCAGGTGCATGCCCGTCCGCTGCTGGTGGTGCCACCGCAGATCAACAAGTTCTACGTGTTCGACCTCAGCCCGGAAAAGAGCCTGGCGCGCTTCTGCCAGCGCTCGCAGCAGCAGACCTtc > 1:293926/1‑136 (MQ=255) cATCACCGAGCAGGTGCATGCCCGTCCGCTGCTGGTGGTGCCACCGCAGATCAACAAGTTCTACGTGTTCGACCTCAGCCCGGAAAAGAGCCTGGCGCGCTTCTGCCTGCGCTCGCAGCAGCAGACCTtcatcatc > 2:364470/1‑136 (MQ=255) aCCGAGCAGGTGCATGCCCGTCCGCTGCTGGTGGTGCCACCGCAGATCAACAAGTTCTACGTGTTCGACCTCAGCCCGGAAAAGAGCCTGGCGCGCTTCTGCCTGCGCTCGCAGCAGCAGACCTTCATCATCAGCt > 2:83769/1‑136 (MQ=255) ccGAGCAGGTGCATGCCCGTCCGCTGCTGGTGGTGCCACCGCAGATCAACAAGTTCTACGTGTTCGACCTCAGCCCGGAAAAGAGCCTGGCGCGCTTCTGCCTGCGCTCGCAGCAGCAGACCTTCATCATCAGCTg > 2:420106/1‑136 (MQ=255) gCAGGTGCATGCCCGTCCGCTGCTGGTGGTGCCACCGCAGATCAACAAGTTCTACGTGTTCGACCTCAGCCCGGAAAAGAGCCTGGCGCGCTTCTGCCTGCGCTCGCAGCAGCAGACCTTCATCATCAGCTGgtgc < 1:162774/136‑1 (MQ=255) gCAGGTGCATGCCCGTCCGCTGCTGGTGGTGCCACCGCAGATCAACAAGTTCTACGTGTTCGACCTCAGCCCGGAAAAGAGCCTGGCGCGCTTCTGCCTGCGCTCGCAGCAGCAGACCTTCATCATCAGCTGgcgc < 1:420106/136‑1 (MQ=255) cATGCCCGTCCGCTGCTGGTGGTGCCACCGCAGATCAACAAGTTCTACGTGTTCGACCTCAGCCCGGAAAAGAGCCTGGCGCGCTTCTGCCTGCGCTCGCAGCAGCAGACCTTCATCATCAGCTGGCGCAACCCGa < 2:81614/136‑1 (MQ=255) cATGCCCGTCCGCTGCTGGTGGTGCCACCGCAGATCAACAAGTTCTACGTGTTCGACCTCAGCCCGGAAAAGAGCCTGGCGCGCTTCTGCCTGCGCTCGCAGCAGCAGACCTTCATCATCAGCTGGCGCAACCCGa > 1:215876/1‑136 (MQ=255) cATGCCCGTCCGCTGCTGGTGGTGCCACCGCAGATCAACAAGTTCTACGTGTTCGACCTCAGCCCGGAAAAGAGCCTGGCGCGCTTCTGCCTGCGCTCGCAGCAGCAGACCTTCATCATCAGCTGGCGCAACCCGa < 2:385410/136‑1 (MQ=255) gctgctGGTGGTGCCACCGCAGATCAACAAGTTCTACGTGTTCGACCTCAGCCCGGAAAAGAGCCTGGCGCGCTTCTGCCTGCGCTCGCAGCAGCAGACCTTCATCATCAGCTGGCGCAACCCGACCAAGGCCCAg > 1:46574/1‑136 (MQ=255) gctgctGGTGGTGCCACCGCAGATCAACAAGTTCTACGTGTTCGACCTCAGCCCGGAAAAGAGCCTGGCGCGCTTCTGCCTGCGCTCGCAGCAGCAGACCTTCATCATCAGCTGGCGCAACCCGACCAAGGCCCAg > 1:307861/1‑136 (MQ=255) gctgctGGTGGTGCCACCGCAGATCAACAAGTTCTACGTGTTCGACCTCAGCCCGGAAAAGAGCCTGGCGCGCTTCTGCCTGCGCTCGCAGCAGCAGACCTTCATCATCAGCTGGCGCAACCCGACCAAGGCCCAg > 1:13062/1‑136 (MQ=255) ctgctgGTGGTGCCACCGCAGATCAACAAGTTCTACGTGTTCGACCTCAGCCCGGAAAAGAGCCTGGCGCGCTTCTGCCTGCGCTCGCAGCAGCAGACCTTCATCATCAGCTGGCGCAACCCGACCAAGGCCCAGc < 1:92453/136‑1 (MQ=255) gctgGTGGTGCCACCGCAGATCAACAAGTTCTACGTGTTCGACCTCAGCCCGGAAAAGAGCCTGGCGCGCTTCTGCCTGCGCTCGCAGCAGCAGACCTTCATCATCAGCTGGCGCAACCCGACCAAGGCCCAGCGt < 2:291045/136‑1 (MQ=255) ctgGTGGTGCCACCGCAGATCAACAAGTTCTACGTGTTCGACCTCAGCCCGGAAAAGAGCCTGGCGCGCTTCTACCTGCGCTCGCAGCAGCAGACCTTCATCATCAGCTGGCGCAACCCGACCAAGGCCCAGCGTg < 1:109957/136‑1 (MQ=255) ctgGTGGTGCCACCGCAGATCAACAAGTTCTACGTGTTCGACCTCAGCCCGGAAAAGAGCCTGGCGCGCTTCTACCTGCGCTCGCAGCAGCAGACCTTCATCATCAGCTGGCGCAACCCGACCAAGGCCCAGCGTg < 1:112025/136‑1 (MQ=255) ggtgCCACCGCAGATCAACAAGTTCTACGTGTTCGACCTCAGCCCGGAAAAGAGCCTGGCGCGCTTCTGCCTGCGCTCGCAGCAGCAGACCTTCATCATCAGCTGGCGCAACCCGACCAAGGCCCAGCGTGAATgg < 2:215876/136‑1 (MQ=255) cACCGCAGATCAACAAGTTCTACGTGTTCGACCTCAGCCCGGAAAAGAGCCTGGCGCGCTTCTGCCTGCGCTCGCAGCAGCAGACCTTCATCATCAGCTGGCGCAACCCGACCAAGGCCCAGCGTGAATGGGGCCt < 1:390063/136‑1 (MQ=255) cACCGCAGATCAACAAGTTCTACGTGTTCGACCTCAGCCCGGAAAAGAGCCTGGCGCGCTTCTGCCTGCGCTCGCAGCAGCAGACCTTCATCATCAGCTGGCGCAACCCGACCAAGGCCCAGCGTGAATGGGGCCt < 1:297545/136‑1 (MQ=255) aCCGCAGATCAACAAGTTCTACGTGTTCGACCTCAGCCCGGAAAAGAGCCTGGCGCGCTTCTGCCTGCGCTCGCAGCAGCAGACCTTCATCATCAGCTGGCGCAACCCGACCAAGGCCCAGCGTGAATGGGGCCTg > 2:111562/1‑136 (MQ=255) gTTCTACGTGTTCGACCTCAGCCCGGAAAAGAGCCTGGCGCGCTTCTGCCTGCGCTCGCAGCAGCAGACCTTCATCATCAGCTGGCGCAACCCGACCAAGGCCCAGCGTGAATGGGGCCTGTCCACCTACATCGAt > 2:219502/1‑136 (MQ=255) gtTCGACCTCAGCCCGGAAAAGAGCCTGGCGCGCTTCTGCCTGCGCTCGCAGCAGCAGACCTTCATCATCAGCTGGCGCAACCCGACCAAGGCCCAGCGTGAATGGGGCCTGTCCACCTACATCGATGCGCTGaaa > 2:307732/1‑136 (MQ=255) aCCTCAGCCCGGAAAAGAGCCTGGCGCGCTTCTGCCTGCGCTCGCAGCAGCAGACCTTCATCATCAGCTGGCGCAACCCGACCAAGGCCCAGCGTGAATGGGGCCTGTCCACCTACATCGATGCGCTGAAAGAAGc > 2:384999/1‑136 (MQ=255) aCCTCAGCCCGGAAAAGAGCCTGGCGCGCTTCTGCCTGCGCTCGCAGCAGCAGACCTTCATCATCAGCTGGCGCAACCCGACCAAGGCCCAGCGTGAATGGGGCCTGTCCACCTACATCGATGCGCTGAAAGAAGc > 1:32893/1‑136 (MQ=255) aaGAGCCTGGCGCGCTTCTGCCTGCGCTCGCAGCAGCAGACCTTCATCATCAGCTGGCGCAACCCGACCAAGGCCCAGCGTGAATGGGGCCTGTCCACCTACATCGATGCGCTGAAAGAAGCCGTCGACGCGGTGc > 2:234580/1‑136 (MQ=255) cgcgcTTCTGCCTGCGCTCGCAGCAGCAGACCTTCATCATCAGCTGGCGCAACCCGACCAAGGCCCAGCGTGAATGGGGCCTGTCCACCTACATCGATGCGCTGAAAGAAGCCGTCGACGCGGTGCTGTCGATTAc < 2:149948/136‑1 (MQ=255) cTTCTGCCTGCGCTCGCAGCAGCAGACCTTCATCATCAGCTGGCGCAACCCGACCAAGG‑CCAGCGTGAATGGGGCCTGTCCACCTACATCGATGCGCTGAAAGAAGCCGTCGACGCGGTGCTGTCGATTACCGgca < 1:364470/136‑1 (MQ=255) | CGATGTGCTGGAACTGATCCAGTACAGCCCCATCACCGAGCAGGTGCATGCCCGTCCGCTGCTGGTGGTGCCACCGCAGATCAACAAGTTCTACGTGTTCGACCTCAGCCCGGAAAAGAGCCTGGCGCGCTTCTGCCTGCGCTCGCAGCAGCAGACCTTCATCATCAGCTGGCGCAACCCGACCAAGGCCCAGCGTGAATGGGGCCTGTCCACCTACATCGATGCGCTGAAAGAAGCCGTCGACGCGGTGCTGTCGATTACCGGCA > NC_002947/5705121‑5705386 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |