| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | CP009273 | 3,178,277 | G→T | 21.1% | E27* (GAA→TAA) | yqiC → | BMFP family putative fusogenic protein |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | CP009273 | 3,178,277 | 0 | G | T | 21.1% | 59.5 / 7.8 | 33 | E27* (GAA→TAA) | yqiC | BMFP family putative fusogenic protein |
| Reads supporting (aligned to +/- strand): ref base G (17/9); new base T (4/3); total (21/12) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 6.86e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
GTTTATGACTTACCCTTATCGCACTACAATGGCACTCAACACCTATCAGTACAGGGAAACTACAATGATTGACCCGAAAAAAATTGAGCAAATCGCTCGCCAGGTTCACGAATCAATGCCTAAAGGAATCAGGGAGTTCGGGGAAGATGTGGAGAAAAAAATCCGCCAAACCCTACAAGCGCAGCTGACGCGCCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGTCCTGTTACGGACTCGTGAAAAACTGGCGCTGCTTGAACAAAGAATCAGCGAAC > CP009273/3178135‑3178424 | aagaaacggtctaaataggctgattcaaggcatttacgggagaaaaaatcggctcaaacatgaagaaatgaaatgactgagtcagccgagaagaatttccccacttattcgcaccttcaCTAAAGGAATCCGGGAGTTCgggg > 2:152522‑M2/120‑143 (MQ=255)aagaaacggtctaaataggctgattcaaggcatttacgggagaaaaaatcggctcaaacatgaagaaatgaaatgactgagtcagccgagaagaatttacccgcttattcgcaccttccCTAAAGGAATCAGGGAGTTCgggg > 2:109893‑M2/120‑143 (MQ=255) tATGACTTACCCTTATCGCACTACAATGGCACTCAACACCTATCAGTACAGGGAAACTACAATGATTGACCCGAAAAAAATTGAGCAAATCGCTCGCCAGGTTCACGAATCAATGCCTAAAGGAATCAGGGAGTTCGGGTAAGATGTgg > 1:525737/1‑149 (MQ=255) aacggtctaaataggctgattcaaggcatttacgggagaaaaaatcggctcaaacatgaagaaatgaaatgactgagtcagccgagaagaatttccccgcttattcgcaccttccCTAACGGAATCAGGGCGCTCGGGGTAGa > 2:139805‑M2/116‑143 (MQ=255) tctaaataggctgattcaaggcatttacgggagaaaaaatcggctcaaacatgaagaaatgaaatgactgagtcagccgagaagaatttccccgcttattcgcaccttccCTAAAGGAATCAGGGAGTTCGGGGAAGATGTGGAGaaaa > 1:1491‑M2/111‑149 (MQ=255) gctgattcaaggcatttacgggagaaaaaatcggctcaaacatgaagaaatgaaatgactgagtcagccgagaagaatttccccgcttattcgcaccttccCTAAAGGAATCAGGGAGTTCGGGGAAGATGTGGAGAAAAAAATCCGcc > 1:469071‑M2/102‑149 (MQ=255) gcatttacgggagaaaaaatcggctcaaacatgaagaaatgaaatgactgagtcagccgagaagaatttccccgcttattcgcaccttccCTAAAGGACTCAGGGAGTTCGGGGAAGATGTGGAGAAAAAAATCAGCCAAAcc > 2:414873‑M2/91‑143 (MQ=255) tttacgggagaaaaaatcggctcaaacatgaagaaatgaaatgactgagtcagccgagaagaatttccccgcttattcgcaccttccCTAAAGGAATCAGGGAGTTCGGGGAAGATGTGGAGAAAAAAATCCGCCAAACCCTACAAgcg > 1:347266‑M2/88‑149 (MQ=255) tttacgggagaaaaaatcggctcaaacatgaagaaatgaaatgactgagtcagccgagaagaatttccccgcttattcgcaccttccCTAAAGGAATCAGGGAGTTCGGGGAAGATGTGGAGAAAAAAATCCGCCAAACCCTACAAgcg > 1:497669‑M2/88‑149 (MQ=255) gggagaaaaaatcggctcaaacatgaagaaatgaaatgactgagtcagccgcgaagaatttccccgcttattcgcaccttccCTAAAGGAATCAGGGAGTTCGGGGAAGATGTGGAGAAAAAAATCCGCCAAACCCTACAAGCGCAGCt > 1:216548‑M2/83‑149 (MQ=255) acCTACCAGTACAGGGAAACTACAATGATTGACCCGAAAAAAATTGAGCAAATCGCTCGCCAGGTTCACGAATCAATGCCTAAAGGAATCAGGGAGTTCGGGTAAGATGTGGAGAAAAAAATCCGCCAAACCCTACAAGCGCa < 2:33877/143‑1 (MQ=255) gaaaaaatcggctcaaacatgaagaaatgaaatgactgagtcagccgagaagaatttccccgcttattcgcaccttccCTAAAGGAATCAGGGAGTTCGGGGAAGATGTGGAGAAAAAAATCCGCCAAACCCTACAAGCGCAGCTGAcg < 1:425871‑M2/71‑1 (MQ=255) aaaaaatcggctcaaacatgaagaaatgaaatgactgagtcagccgagaagaatttccccgcttattcgcaccttccCTAAAGGAATCAGGGAGTTCGGGGAAGATGTGGAGAAAAAAATCCGCCAAACCCTACAAGCGCAGCTGAcgc > 1:175707‑M2/78‑149 (MQ=255) aaaaaatcggctcaaacatgaagaaatgaaatgactgagtcagccgagaagaatttccccgcttattcgcaccttccCTAAAGGAATCAGGGAGTGCGGGGAAGATGTGGAGAAAAAAATCCGCCAAACCCTACAAGCGCAGCTGAcgc > 1:175694‑M2/78‑149 (MQ=255) aaaaatcggctcaaacatgaagaaatgaaatgactgagtcagccgagaagaatttccccgcttcttcgcaccttccCTAAAGGAATCAGGGAGTTCGGGGAAGATGTGGAGAAAAAAGTCCGCCAAACCCTACAAGCGCAGCTGAcgcg > 1:353715‑M2/77‑149 (MQ=255) aaaatcggctcaaacatgaagaaatgaaatgactgagtcagccgagaagaatttccccgcttattcgcaccttccCTAAAGGAATCAGGGAGTTCGGGGAAGATGTGGAGAAAAAAATCCGCCAAACCCTACAAGCGCAGCTg > 2:96324‑M2/76‑143 (MQ=255) ggctcaaacatgaagaaatgaaatgactgagtcagccgagaagaatttccccgcttattcgcaccttccCTAAAGGAATCAGGGAGTTCGGGGAAGATGTGGAGAAAAAAATCCGCCAAACCCTGCAAGCGCAGCTGACGCGCCCCGAc > 1:21933‑M2/70‑149 (MQ=255) aGGGAAACTACAATGATTGACCCGAAAAAAATTGAGCAAATCGCTCGCCAGGTTCACGAATCAATGCCTAAAGGAATCAGGGAGTTCGGGTAAGCTGTGGAGAAAACAATCCGCCAAACCCTACAAGACCACCTGCCGCGCCt > 2:14389/1‑143 (MQ=255) tGATTGACCCGAAAAAAATTGAGCAAATCGCTCGCCAGGTTCACGAATCAATGCCTAAAGGAATCAGGGAGTTCGGGTAAGATGTGGAGAAAAAAATCCGCCAAACCCTACAAGCGCAGCTGACGCGCCTCGACCTTGTAAGc > 2:4797/1‑143 (MQ=255) ttGACCCAAAAAAAACTGAGCAACTCGCTCGCCAGTTCCCCGAATCAATGCCTAAAGGAATCAGGGAGTTCGGGTAAGATGTGGAGAAAAAAATCCGCCACCCCCTACGAGCGCAGCTGACGCGCCTCGACCTTGTAAGCCGt < 2:376108/143‑1 (MQ=255) aaatgactgagtcagccgagaagaatttccccgcttattcgcaccttccCTAAAGGAATCAGGGAGTTCGGGGAAGATGTGGAGAAAAAAATCCGCCAAACCCTACAAGCGCAGCTGACGCGCCTCGACCTTGTAAGCCGTGa > 2:180913‑M2/50‑143 (MQ=255) aatgactgagtcagccgagaagaatttccccgcttattcgcaccttccCTAAAGGAATCAGGGAGTTCGGGGAAGATGTGGAGAAAAAAATCCGCCAAACCCTACAAGCGCAGCTGACGCGCCTCGACCTTGTAAGCCGTGAAGAGTTc < 1:339639‑M2/101‑1 (MQ=255) actgagtcagccgagaagaatttccccgcttattcgcaccttccCTAAAGGAATCAGGGAGTTCGGGGACGATGAGGAGAAAAAAATCCGCCAAACCCTACAAGCGCAGCTGACGCGCCTCGACCTTGTAAGCCGTGAACAGTTCCACg > 1:197489‑M2/45‑149 (MQ=255) ctgagtcagccgagaagaatttccccgcttattcgcaccttccCTAAAGGAATCAGGGAGTTCGGGGAAGATGTGGAGAAAAAAATCCGCCAAACCCTACAAGCGCAGCTGACGCGCCTCGACCTTGTAAGCCGTGAAGAGtt < 2:314747‑M2/100‑1 (MQ=255) aacaaaTTGAGCAACTCCGTGGCCAGGTTCACGACTCAATGCCTAAAGGAATCAGGGAGTGCGGGTAAGATGTGGAGAAAAAAATCCGCCAAACCCTACAAGCGCAGCTGACGCGCCTCTACCTTGTAAGCCGTGAAGAGTTc < 2:9237/140‑1 (MQ=255) gagtcagccgagaagaatttccccgcttattcgcaccttccCTAAAGGAATCAGGGAGTTCGGGGAAGATGTGGAGAAAAAAATCCACCAAACCCTACAAGCGCAGCTGACGCGCCTCGCCCTGCTAAGCCGTGAAGAGTTCGACGGcc > 1:357344‑M2/42‑149 (MQ=255) agtcagccgagaagaatttccccgcttattcgcaccttccCTAAAGGAATCAGGGAGTTCGGGGAAGATGTGGAGAAAAAAATCCGCCAAACCCTACAAGCGCAGCTGACGCGCCTCGACCTTGTAAGCCGTGAAGAGTTCGa > 2:425224‑M2/41‑143 (MQ=255) tcagccgagaagaatttccccgcttattcgcaccttccCTAAAGGAATCAGGGAGTTCGGGGAAGATGTGGAGAAAAAAATCCGCCAAACCCTACAAGCGCAGCTGACGCGCCTCGACCTTGTAAGCCGTGAAGAGTTCGACg < 2:53298‑M2/105‑1 (MQ=255) cagccgagaagaattcccccgcttattcgcaccttccCTAAAGGAATCAGGGAGTTCGGGGAAGATGTGGAGAAACAAATCCGCCAAACCCTACAAGCGCAGCTGACGCGCCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAAc < 1:338884‑M2/112‑1 (MQ=255) agaatttccccgcttattcgcacgttccCTAGAGGAATCAGGGCGTCCGTGGAATATGTGGAGAAAAAAATCCGCCAAACCCTACAAGCGCAGCTGACGCGCCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGTCCt < 1:516423‑M2/121‑1 (MQ=255) aTCGCTCGCCAGGTTCACGAATCAATGCCTAAAGGAATCAGGGAGTTCGGGTAAGATGTGGAGAAAAAAATCCGCCAAACCCTACAAGCGCAGCTGACGCGCCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCa > 2:182799/1‑143 (MQ=255) gcttattcgcaccttccCTCACGTAATCAGGGAGTTCGGGGAAGATGTGGAGAAAAAAATCCGCCAAACCCTACAAGCGCAGCTGACGCGCCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGTCCTGTTAc < 2:191208‑M2/126‑1 (MQ=255) tattcgcaccttccCTAAAGGAATCAGGGAGTTCGGGGAAGATGTGGAGAAAAAAATCCGCCAAACCCTACAAGCGCAGCTGACGCGCCTCGACCTTGTAAGCCGTGAAGAGTGCGACGGCCAAACGCAAGTCCTGTTACGGACTCGTg > 1:492479‑M2/15‑149 (MQ=255) cCTAAAGGAATCAGGGAGTTCGGGGAAGATGTGGAGAAAAAAATCCGCCAAGCCCTACAAGCGCAGCTGACGCGCCTCGACCTTGTAAGCCGTGAAGAGTTCGGCGTCCAAACGCAAGTCCTGTTACGGACTCGTGAAAAACt < 2:436647‑M2/142‑1 (MQ=255) gggAGTTCGGGGAAGATGTGGAGAAAAAAATCCGCCAAACCCTACAAGCGCAGCTGACGCGCCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGTCCTGTTACGGACTCGTGAAAAACTGGCGCTGCTTGaa < 2:110652/143‑1 (MQ=255) ggAAGATGTGGAGAAAAAAATCCGCCAAACCCTACAAGCGCAGCTGACGCGCCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGTCCTGTTACGGGCTCGTGAAAAACTGGCGCTGCTTGAACAAAGAATCAGCGAAc < 1:393139/149‑1 (MQ=255) | GTTTATGACTTACCCTTATCGCACTACAATGGCACTCAACACCTATCAGTACAGGGAAACTACAATGATTGACCCGAAAAAAATTGAGCAAATCGCTCGCCAGGTTCACGAATCAATGCCTAAAGGAATCAGGGAGTTCGGGGAAGATGTGGAGAAAAAAATCCGCCAAACCCTACAAGCGCAGCTGACGCGCCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGTCCTGTTACGGACTCGTGAAAAACTGGCGCTGCTTGAACAAAGAATCAGCGAAC > CP009273/3178135‑3178424 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |