Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | AE015451 | 2,414,802 | C→G | 22.8% | intergenic (+344/‑1) | PP_2114 → / → PP_2115 | transposase/conserved protein of unknown function |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | AE015451 | 2,414,802 | 0 | C | G | 22.8% | 52.1 / 13.4 | 22 | intergenic (+344/‑1) | PP_2114/PP_2115 | transposase/conserved protein of unknown function |
Reads supporting (aligned to +/- strand): ref base C (8/9); new base G (2/3); total (10/12) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.30e-03 |
TCAGCGAACAGTCCATTCGCGATGGCGACACTGACTCACCTTTTCGCCCTTACGGCGAGTCACTTTTTGTCAAACGCGACAAAAAGTAACCAAAAAACGCTGCGCTCCCATCATCCGGCCCCTGCGCTGCGCTCCGGGGTCCCCTCACTCCGGCCTTGCTCCCGGCAGGACCGCGCCGAAGGCCCCATCCTGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTTCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCGAGATCAAGATCAAGATCGAGATCAAGAGCGCTGGGTTTGAAAGCCGTGCTTGGACCGCGTCGCG > AE015451/2414586‑2414984 | tCAGCGAACAGTCCATTCGCGATGGCGACACTGACTCACCTTTTCGCCCTTACGGCGAGTCACTTTTTGTCAAACGCGACAAAAAGTAACCAAAAAACGCTGCGCTCCCATCATCCGGCCCCTGCGCTGCGCTCCGGGGTCCCCTCACTCCGGCCTTGCTCCCGGCAGGACCGCGCCGAAGGCCCCATCCTGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTTc < 1:116290/236‑1 (MQ=38) aaCAGTCCATTCGCGATGGCGACACTGACTCACCTTTTCGCCCTTACGGCGAGTCACTTTTTGTCAAACGCGACAAAAAGTAACCAAAAAACGCTGCGCTCCCATCATCCGGCCCCTGCGCTGCGCTCCGGGGTCCCCTCACTCCGGCCTTGCTCCCGGCAGGACCGCGCCGAAGGCCCCATCCTGGGGCCTCAGCGCTTGACGGGCATCCATg > 1:429356/1‑214 (MQ=39) aaCAGTCCATTCGCGATGGCGACACTGACTCACCTTTTCGCCCTTACGGCGAGTCACTTTTTGTCAAACGCGACAAAAAGTAACCAAAAAACGCTGCGCTCCCATCATCCGGCCCCTGCGCTGCGCTCCGGGGTCCCCTCACTCCGGCCTTGCTCCCGGCAGGACCGCGCCGAAGGCCCCATCCTGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTTCGCAAgg > 1:391473/1‑236 (MQ=40) tCCATTCGCGATGGCGACACTGACTCACCTTTTCGCCCTTACGGCGAGTCACTTTTTGTCAAACGCGACAAAAAGTAACCAAAAAACGCTGCGCTCCCATCATCCGGCCCCTGCGCTGCGCTCCGGGGTCCCCTCACTCCGGCCTTGCTCCCGGCAGGACCGCGCCGAAGGCCCCATCCTGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTTCGCAAGGCCTGc < 1:470236/236‑1 (MQ=38) cacTGACTCACCTTTTCGCCCTTACGGCGAGTCACTTTTTGTCAAACGCGACAAAAAGTAACCAAAAAACGCTGCGCTCCCATCATCCGGCCCCTGCGCTGCGCTCCGGGGTCCCCTCACTCCGGCCTTGCTCCCGGCAGGACCGCGCCGAAGGCCCCATCCTGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTTCGCAAGGCCTGCGtt < 1:550008/222‑1 (MQ=34) aCCTTTTCGCCCTTACGGCGAGTCACTTTTTGTCAAACGCGACAAAAAGTAACCAAAAAACGCTGCGCTCCCATCATCCGGCCCCTGCGCTGCGCTCCGGGGTCCCCTCACTCCGGCCTTGCTCCCGGCAGGACCGCGCCGAAGGCCCCATCCTGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTTc < 1:56070/199‑1 (MQ=32) aGTCACTTTTTGTCAAACGCGACAAAAAGTAACCAAAAAACGCTGCGCTCCCATCATCCGGCCCCTGCGCTGCGCTCCGGGGTCCCCTCACTCCGGCCTTGCTCCCGGCAGGACCGCGCCGAAGGCCCCATCCTGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTTCGCAAg < 2:189818/184‑1 (MQ=31) cGCTGCGCTCCCATCATCCGGCCCCTGCGCTGCGCTCCGGGGTCCCCTCACTCCGGCCTTGCTCCCGGCAGGACCGCGCCGAAGGCCCCATCCTGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTTCGCAAGGCCTGCGTTCGGcctcct < 1:19967/162‑1 (MQ=31) cGGCCCCTGCGCTGCGCTCCGGGGTCCCCTCACTCCGGCCTTGCTCCCGGCAGGACCGCGCCGAAGGCCCCATCCTGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTTCGCAAGGCCTGCGCTCGGCCTCCTGAAGCCGCGAAGGTCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAagatcaagatcaagatcgagatcaagatcaagat > 1:175238/1‑236 (MQ=255) ggCCCCTGCGCTGCGCTCCGGGGTCCCCTCACTCCGGCCTTGCTCCCGGCAGGACCGCGCCGAAGGCCCCATCCTGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTTCGCAAGGCCTGCGTTCGGCCTCCTGAAGTc < 1:250480/149‑1 (MQ=32) ccccTGCGCTGCGCTCCGGGGTCCCCTCACTCCGGCCTTGCTCCC‑‑‑‑GGACCGCGCCGAAGGCCCCATCCTGGGGCCTCAGCGCTTGACGGGCATGGATGCCCGTCACCTCCCTTCGCAAGCCCTACGCTCGGCCTCCTGAAGTcgcaatct < 1:358405/150‑6 (MQ=11) ccccTGCGCTGCGCTCCGGGGTCCCCTCACTCCGGCCTTGCTCCC‑‑‑‑GGACCGCGCCGAAGGCCCCATCCTGGGGCCTCAGCGCTTGACGGGCATGGATGCCCGTCACCTCCCTTCGCAAGCCCTACGCTCGGCCTCCTGAAGTcgcaatct > 2:358404/1‑145 (MQ=11) ccTGCGCTGCGCTCCGGGGTCCCCTCACTCCGGCCTTGCTCCCGGCAGGACCGCGCCGAAGGCCCCATCCTGGGGCCTCAGCGCTTGACGGGCATGGATGCCCGTCACCTGCCTTc > 1:486087/1‑116 (MQ=9) ccTGCGCTGCGCTCCGGGGTCCCCTCACTCCGGCCTTGCTCCCGGCAGGACCGCGCCGAAGGCCCCATCCTGGGGCCTCAGCGCTTGACGGGCATGGATGCCCGTCACCTGCCTTc < 2:486086/116‑1 (MQ=9) ggTCCCCTCACTCCGGCCTTGCTCCCGGCAGGACCGCGCCGAAGGCCCCATCCTGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTTCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGaagatcaa > 2:260436/1‑139 (MQ=33) ccTCACTCCGGCCTTGCTCCCGGCAGGACCGCGCCGAAGGCCCCATCCTGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTTCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCGAGATCAAGATCAAGATCGAGATCAAGAGCGCTGGGATTGa > 1:71540/1‑233 (MQ=255) cAGGACCGCGCCGAAGGCCCCATCCTGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTTCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCGAGATCAAGATCAAGATCGAGATCAAGAGCGCTGGGTTTGaa < 2:688259/211‑1 (MQ=255) cGAAGGCCCCATCCTGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTTCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCGAGATCAAGATCAAGATCGAGATCAAGAGCGCTGGGTTTGAAAGCCGTGCTTGGACCGAGTcgag > 1:61621/1‑221 (MQ=255) cATCCTGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTTCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCGAGATCAAGATCAAGAGCGAGATCAAGAGCGCTGGGTTTGaa < 1:260437/191‑1 (MQ=255) gcTTGACGGGCATCCATGCCCGTCACCTGCCTTCGCAAGGGCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAGATGAAGATCAAGATCAAGATCAAGATCAagatcaagatcaagatcgagatga > 2:681538/1‑137 (MQ=255) aggtGACGGGCATCCATGCCCGTCACCTGCCTTCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCGAGATCAAGATCAAGATCGAGATCAAGAGCGCTGGGATTGAAAGCCggg > 2:100536/4‑179 (MQ=255) gggCATGGATGCCCGTCACCTGCCTTCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGaagatcaagatcaagatcaagatcaagatcaagatcaag < 1:606557/98‑1 (MQ=255) | TCAGCGAACAGTCCATTCGCGATGGCGACACTGACTCACCTTTTCGCCCTTACGGCGAGTCACTTTTTGTCAAACGCGACAAAAAGTAACCAAAAAACGCTGCGCTCCCATCATCCGGCCCCTGCGCTGCGCTCCGGGGTCCCCTCACTCCGGCCTTGCTCCCGGCAGGACCGCGCCGAAGGCCCCATCCTGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTTCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCGAGATCAAGATCAAGATCGAGATCAAGAGCGCTGGGTTTGAAAGCCGTGCTTGGACCGCGTCGCG > AE015451/2414586‑2414984 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |