Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | CP009273 | 374,648 | Δ1 bp | 25.5% | coding (380/1110 nt) | frmA ← | alcohol dehydrogenase class III/glutathione‑dependent formaldehyde dehydrogenase |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | CP009273 | 374,648 | 0 | T | . | 25.5% | 143.6 / 47.2 | 51 | coding (380/1110 nt) | frmA | alcohol dehydrogenase class III/glutathione‑dependent formaldehyde dehydrogenase |
Reads supporting (aligned to +/- strand): ref base T (18/20); new base . (6/7); total (24/27) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.63e-01 |
AATACCGGTGGTCACGCCACAGCCCAGCAGGCAGACGTGTTCATGGTTTGCTTCTGGATTAATTTTTGCCAGAGACACTTCCGCGACTACGGTGTATTCACTGAATGTAGAGCACCCCATGTAGTGATAAAGCGGCTGCCCGTTGTAAGAAAAACGGGTGGTGCCGTCTGGCATCAGGCCTTTACCCTGGGTTTCGCGAACCGCAACACAGAGGTTAGTTTTGCCAGAACGACAGAACTCACACTCGCCGCACTCCGCGGTGTAAAGCGGGATCACATGGTCGCCAGGTTTGAC > CP009273/374503‑374796 | ccggtggtcacgccacagcccagcaggcagacgtgttcatggtttgcttctggattaatttttgccagagacacttccgcgactacggtgtattcactgaatgtagagcaccccatgtagtgataaagcggctgcccGCCCGTTg‑aagaa < 1:13526‑M2/13‑1 (MQ=255) cgccacagcccagcaggcagacgtgttcatggtttgcttccggattaatttttgccagagacacttccgcgactacggtgtattcactgaatgtagagcaccccatgtagtgataaagcggctgcccGCCCGTTGTAAGAAAAACGggtg < 1:65120‑M2/23‑1 (MQ=255) gccacagcccagcaggcagacgtgttcatggtttgcttctggattaatttttgccagagacacttccgcgactacggtgtattcactgaatgtagagcaccccatgtagtgataaagcggctgcccGCCCGTTGTAAGAAAAACGggtgg > 2:870375‑M2/127‑150 (MQ=255) ccacagcccagcaggcagacgtgttcatggtttgcttctggattaatttttgccagagacacttccgcgactacggtgtattcactgaatgtagagcaccccatgtagtgataaagcggctgcccGCCCGTTG‑AAGAAAAACGGgtggtg < 1:1058272‑M2/25‑1 (MQ=255) acagcccagcaggcagacgtgttcatggtttgcttctggattaatttttgccagagacacttccgcgactacggtgtattcactgaatgtagagcaccccatgtagtgataaagcggctgcccGCCCGTTGTAAGAAAAACGGGTGGTGc > 2:573906‑M2/124‑150 (MQ=255) cagcccagcaggcagacgtgttcatggtttgcttctggattaatttttgccagagacacttccgcgactacggtgtattcactgaatgtagagcaccccatgtagtgataaagcggctgcccGCCCGTTGTAAGAAAAACGGGTGGTGcc < 1:421109‑M2/28‑1 (MQ=255) gacgtgttcatggtttgcttctggattaatttttgccagagacacttccgcgactacggtgtattcactgaatgtagagcaccccatgtagtgataaagcggctgcccGCCCGTTGTAAGAAAAACGGGTGGTGCCGTCTGGCATCAGGc < 2:291337‑M2/42‑1 (MQ=255) gacgtgttcatggtttgcttctggattaatttttgccagagacacttccgcgactacggtgtattcactgaatgtagagcaccccatgtagtgataaagcggctgcccGCCCGTTG‑AAGAAAAACGGGTGGTGCCGTCTGGCATCAGGcc > 2:1035975‑M2/109‑150 (MQ=255) atggtttgcttctggattaatttttgccagagacacttccgcgactacggtgtattcactgaatgtagagcaccccatgtagtgataaagcggctgcccGCCCGTTGTAAGAAAAACGGGTGGTGCCGTCTGGCATCAGGCCTTTACCCt > 1:843158‑M2/100‑150 (MQ=255) gtttgcttctggattaatttttgccagagacacttccgcgactacggtgtattcactgaatgtagagcaccccatgtagtgataaagcggctgcccGCCCGTTGTAAGAAAAACGGGTGGTGCCGTCTGGCATCAGGCCTTTACCCTggg > 2:97011‑M2/97‑150 (MQ=255) gtttgcttctggattaatttttgccagagacacttccgcgactacggtgtattcactgaatgtagagcaccccatgtagtgataaagcggctgcccGCCCGTTGTAAGAAAAACGGGTGGTGCCGTCTGGCATCAGGCCTTTACCCTggg < 1:85183‑M2/54‑1 (MQ=255) gattaatttttgccagagacacttccgcgactacggtgtattcactgaatgtagagcaccccatgtagtgataaagcggctgcccGCCCGTTG‑AAGAAAAACGGGTGGTGCCGTCTGGCATCAGGCCTTTACCCTGGGTTTCGCGAAcc < 1:647257‑M2/64‑1 (MQ=255) aatttttgccagagacacttccgcgactacggtgtattcactgaatgtagagcaccccatgtagtgataaagcggctgcccGCCCGTTGTAAGAAAAACGGGTGGTGCCGTCTGGCATCAGGCCTTTACCCTGGGTTTCGCGAACCGCa < 1:488471‑M2/68‑1 (MQ=255) atttttgccagagacacttccgcgactacggtgtattcactgaatgtagagcaccccatgtagtgataaagcggctgcccGCCCGTTGTAAGAAAAACGGGTGGTGCCGTCTGGCATCAGGCCTTTACCCTGGGTTTCGCGAACCGCaa < 1:639783‑M2/69‑1 (MQ=255) tttgccagagacacttccgcgactacggtgtattcactgaatgtagagcaccccatgtagtgataaagcggctgcccGCCCGTTGTAAGAAAAACGGGTGGTGCCGTCTGGCATCAGGCCTTTACCCTGGGTTTCGCGAACCGCAaca > 1:307692‑M2/78‑148 (MQ=255) tgccagagacacttccgcgactacggtgtattcactgaatgtagagcaccccatgtagtgataaagcggctgcccGCCCGTTGTAAGAAAAACGGGTGGTGCCGTCTGGCATCAGGCCTTTACCCTGGGTTTCGCGAACCGCAACACag < 1:398153‑M2/74‑1 (MQ=255) tgccagagacacttccgcgactacggtgtattcactgaatgtagagcaccccatgtagtgataaagcggctgcccGCCCGTTGTAAGAAAAACGGGTGGTGCCGTCTGGCATCAGGCCTTTACCCTGGGTTTCGCGAACCGCAACACag < 1:870375‑M2/74‑1 (MQ=255) ccagagacacttccgcgactacggtgtattcactgaatgtagagcaccccatgtagtgataaggcggctgcccGCCCGTTGTAAGAAAAACGGGTGGTGCCGTCTGGCATCAGGCCTTTACCCTGGGTTTCGCGAACCGCAACACagag < 1:97011‑M2/76‑1 (MQ=255) cagagacacttccgcgactacggtgtattcactgaatgtagagcaccccatgtagtgataaagcggctgcccGCCCGTTGTAAGAAAAACGGGTGGTGCCGTCTGGCATCAGGCCTTTACCCTGGGTTTCGCGAACCGCAACACAGAGGt < 2:678982‑M2/78‑1 (MQ=255) agagactcttccgcgactacggtgtattcactgaatgtagagcaccccatgtagtgataaagcggctgcccGCCCGTTGTAAGAAAAACGGGTGGTGCCGTCTGGCATCAGGCCTTTACCCTGGGTTTCGCGAAGCGGAACACAGAGGtt < 2:674496‑M2/79‑1 (MQ=255) agagacacttccgcgactacggtgtattcactgaatgtagagcaccccatgtagtgataaagcggctgcccGCCCGTTG‑AAGAAAAACGGGTGGTGCCGTCTGGCATCAGGCCTTTACCCTGGGTTTCGCGAACCGCAACACAGAGGt > 1:139326‑M2/72‑148 (MQ=255) gagacacttccgcgactacggtgtattcactgaatgtagagcaccccatgtagtgataaagcggctgcccGCCCGTTGTAAGAAAAACGGGTGGTGCCGTCTGGCATCAGGCCTTTACCCTGGGTTTCGCGAACCGCAACACAGAGGt > 2:798058‑M2/71‑148 (MQ=255) gacacttccgcgactacggtgtattcactgaatgtagagcaccccatgtagtgataaagcggctgcccGCCCGTTG‑AAGAAAAACGGGTGGTGCCGTCTGGCATCAGGCCTTTACCCTGGGTTTCGCGAACCGCAACACAGAGGTTAGtt > 2:975567‑M2/69‑150 (MQ=255) acacttccgcgactacggtgtattcactgaatgtagagcaccccatgtagtgataaagcggctgcccGCCCGTTGTAAGAAAAACGGGTGGTGCCGTCTGGCATCAGGCCTTTACCCTGGGTTTCGCGAACCGCAACACAGAGGTTAGtt > 2:812891‑M2/68‑150 (MQ=255) cttccgcgactacggtgtattcactgaatgtagagcaccccatgtagtgataaagcggctgcccGCCCGTTGTAAGAAAAACGGGTGGTGCCGTCTGGCATCAGGCCTTTACCCTGGGTTTCGCGAACCGCAACACAGAGGTTAGTTTTg < 1:873745‑M2/86‑1 (MQ=255) actacggtgtattcactgaatgtagagcaccccatgtagtgataaagcggctgcccGCCCGTTGTAAGAAAAACGGGTGGTGCCGTCTGGCATCAGGCCTTTACCCTGGGTTTCGCGAACCGCAACACAGAGGTTAGTTTTGCCAGAACg < 1:909678‑M2/94‑1 (MQ=255) ctacggtgtattcactgaatgtagagcaccccatgtagtgataaagcggctgcccGCCCGTTGTAAGAAAAACGGGTGGTGCCGTCTGGCATCAGGCCTTTACCCTGGGTTTCGCGAACCGCAACACAGAGGTTAGTTTTGCCAGAACGa > 2:802714‑M2/56‑150 (MQ=255) cggtgtattcactgaatgtagagcaccccatgtagtgataaagcggctgcccGCCCGTTGTAAGAAAAACGGGTGGTGCCGTCTGGCATCAGGCCTTTACCCTGGGTTTCGCGAACCGCAACACAGAGGTTAGTTTTGCCAGAACGAc < 1:573906‑M2/96‑1 (MQ=255) ttcactgaatgtagagcaccccatgtagtgataaagcggctgcccGCCCGTTGTAAGAAAAACGGGTGGTGCCGTCTGGCATCAGGCCTTTACCCTGGGTTTCGCGAACCGCGACACAGAGGTTAGTTTTGCCAGAACGACAGAACTca > 1:643050‑M2/46‑149 (MQ=255) actgaatgtagagcaccccatgtagtgataaagcggctgcccGCCCGTTGTAAGAAAAACGGGTGGTGCCGTCTGGCATCAGGCCTTTACCCTGGGTTTCGCGAACCGCAACACAGAGGTTAGTTTTGCCAGAACGACAGAACTCACACt > 1:803809‑M2/43‑150 (MQ=255) tgtagagcaccccatgtagtgataaagcggctgcccGCCCGTTG‑AAGAAAAACGGGTGGTGCCGTCTGGCATCAGGCCTTTACCCTGGGTTTCGCGAACCGCAACACAGAGGTTAGTTTTGCCAGAACGACAGAACTCACACTCGCCGCa < 2:201552‑M2/114‑1 (MQ=255) agagcaccccatgtagtgataaagcggctgcccGCCCGTTGTAAGAAAAACGGGTGGTGCCGTCTGGCATCAGGCCTTTACCCTGGGTTTCGCGAACCGCAACACAGAGGTTAGTTTTGCCAGAACGACAGAACTCACACTCGCCGCa > 2:857800‑M2/34‑148 (MQ=255) gagcaccccatgtagtgataaagcggctgcccGCCCGTTGTAAGAAAAACGGGTGGTGCCGTCTGGCATCAGGCCTTTACCCTGGGTTTCGCGAACCGCAACACAGAGGTTAGTTTTGCCAGAACGACAGAACTCACACTCGCCGCAc < 1:857800‑M2/116‑1 (MQ=255) gcaccccatgtagtgataaagcggctgcccGCCCGTTGTAAGAAAAACGGGTGGTGCCGTCTGGCATCAGGCCTTTACCCTGGGTTTCGCGAACCGCAACACAGAGGTTAGTTTTGCCAGAACGACAGAACTCACACTCGCCGCACTCcg < 2:560285‑M2/120‑1 (MQ=255) taccccatgtagtgataaagcggctgcccGCCCGTTGTAAGAAAAACGGGTGGTGCCGTCTGGCATCAGGCCTTTACCCTGGGTTTCGCGAACCGCAACACAGAGGTTAGTTTTGCCAGAACGACAGAACTCACACTCGCCGCACTCcgc < 2:413824‑M2/121‑1 (MQ=255) accccatgtagtgataaagcggctgcccGCCCGTTGTAAGAAAAACGGGTGGTGCCGTCTGGCATCAGGCCTTTACCCTGGGTTTCGCGAACCGCAACACAGAGGTTAGTTTTGCCAGAACGACAGAACTCACACTCGCCGCACTCcgcg > 2:103579‑M2/29‑150 (MQ=255) ccccatgtagtgataaagcggctgcccGCCCGTTGTAAGAAAAACGGGTGGTGCCGTCTGGCATCAGGCCTTTACCCTGGGTTTCGCGAACCGCAACACAGAGGTTAGTTTTGCCAGAACGACAGAACTCACACTCGCCGCACTCCGCgg > 2:726775‑M2/28‑150 (MQ=255) ccatgtagtgataaagcggctgcccGCCCGTTGTAAGAAAAACGGGTGGTGCCGTCTGGCATCAGGCCTTTACCCTGGGTTTCGCGAACCGCAACACAGAGGTTAGTTTTGCCAGAACGACAGAACTCACACTCGCCGCACTCCGCGgtg < 2:843158‑M2/125‑1 (MQ=255) ccatgtagtgataaagcggctgcccGCCCGTTG‑AAGAAAAACGGGTGGTGCCGTCTGGCATCAGGCCTTTACCCTGGGTTTCGCGAACCGCAACACAGAGGTTAGTTTTGCCAGAACGACAGAACTCACACTCGCCGCACTCCGCGgtg > 2:922192‑M2/26‑149 (MQ=255) atgtagtgataaagcggctgcccGCCCGTTGTAAGAAAAACGGGTGGTGCCGTCTGGCATCAGGCCTTTACCCTGGGTTTCGCGAACCGCAACACAGAGGTTAGTTTTGCCAGAACGACAGAACTCACACTCGCCGCACTCCGCGgtgt > 1:907951‑M2/24‑149 (MQ=255) tagtgataaagcggctgcccGCCCGTTGTAAGAAAAACGGGTGGTGCCGTCTGGCATCAGGCCTTTACCCTGGGTTTCGCGAACCGCAACACAGAGGTTAGTTTTGCCAGAACGACAGAACTCACACTCGCCGCACTCCGCGGTGTAAAg < 2:1024869‑M2/130‑1 (MQ=255) tagtgataaagcggctgcccGCCCGTTG‑AAGAAAAACGGGTGGTGCCGTCTGGCATCAGGCCTTTACCCTGGGTTTCGCGAACCGCAACACAGAGGTTAGTTTTGCCAGAACGACAGAACTCACACTCGCCGCACTCCGCGGTGTAAAGc < 2:480290‑M2/130‑1 (MQ=255) gataaagcggctgcccGCCCGTTG‑AAGAAAAACGGGTGGTGCCGTCTGGCATCAGGCCTTTACCCTGGGTTTCGCGAACCGCAACACAGAGGTTAGTTTTGCCAGAACGACAGAACTCACACTCGCCGCACTCCGCGGTGTAAAGCGGGa > 1:619163‑M2/17‑150 (MQ=255) ataaagcggctgcccGCCCGTTG‑AAGAAAAACGGGTGGTGCCGTCTGGCATCAGGCCTTTACCCTGGGTTTCGCGAACCGCAACACAGAGGTTAGTTTTGCCAGAACGACAGAACTCACACTCGCCGCACTCCGCGGTGTAAAGCGGGa < 1:113029‑M2/134‑1 (MQ=255) ataaagcggctgcccGCCCGTTG‑AAGAAAAACGGGTGGTGCCGTCTGGCATCAGGCCTTTACCCTGGGTTTCGCGAACCGCAACACAGAGGTTAGTTTTGCCAGAACGACAGAACTCACACTCGCCGCACTCCGCGGTGTAAAGCGGGAt < 2:560106‑M2/135‑1 (MQ=255) aaagcggctgcccGCCCGTTGTAAGAAAAACGGGTGGTGCCGTCTGGCATCAGGCCTTTACCCTGGGTTTCGCGAACCGCAACACAGAGGTTAGTTTTGCCAGAACGACAGAACTCACACTCGCCGCACTCCGCGGTGTAAAGCGGGAt > 1:993473‑M2/14‑149 (MQ=255) tgcccGCCCGTTG‑AAGAAAAACGGGTGGTGCCGTCTGGCATCAGGCCTTTACCCTGGGTTTCGCGAACCGCAACACAGAGGTTAGTTTTGCCAGTACGACAGAACTCACACTCGCCGCACTCCGCGGTGTAAAGCGGGATCACATGGTCg > 2:867632‑M2/6‑150 (MQ=255) gcccGCCCGTTGTAAGAAAAACGGGTGGTGCCGTCTGGCATCAGGCCTTTACCCTGGGTTTCGCGAACCGCAACACAGAGGTTAGTTTTGCCAGAACGACAGAACTCACACTCGCCGCACTCCGCGGTGTAAAGCGGGATCACATGGTCg > 1:281704‑M2/5‑150 (MQ=255) cccGCCCGTTGTAAGAAAAACGGGTGGTGCCGTCTGGCATCAGGCCTTTACCCTGGGTTTCGCGAACCGCAACACAGAGGTTAGTTTTGCCAGAACGACAGAACTCACACTCGCCGCACTCCGCGGTGTAAAGCGGGATCACATGGTCg < 2:993473‑M2/146‑1 (MQ=255) cGCCCGTTG‑AAGAAAAACGGGTGGTGCCGTCTGGCATCAGGCCTTTACCCTGGGTTTCGCGAACCGCAACACAGAGGTTAGTTTTGCCAGAACGACAGAACTCACACTCGCCGCACTCCGCGGTGTAAAGCGGGATCACATGGTCGCCAg < 1:975567‑M2/149‑1 (MQ=255) cccGTTGTAAGAAAAACGGGTGGTGCCGTCTGGCATCAGGCCTTTACCCTGGGTTTCGCGAACCGCAACACAGAGGTTAGTTTTGCCAGAACGACAGAACTCACACTCGCCGCACTCCGCGGTGTAAAGCGGGATCACATGGTCGCCAgg > 1:906633/1‑150 (MQ=255) ttGTAAGAAAAACGGGTGGTGCCGTCTGGCATCAGGCCTTTACCCTGGGTTTCGCGAACCGCAACACAGAGGTTAGTTTTGCCAGAACGACAGAACTCACACTCGCCGCACTCCGCGGTGTAAAGCGGGATCACATGGTCGCCAGGttt > 2:391121/1‑149 (MQ=255) tAAGAAAAACGGGTGGTGCCGTCTGGCATCAGGCCTTTACCCTGGGTTTCGCGAACCGCAACACAGAGGTTAGTTTTGCCAGAACGACAGAACTCACACTCGCCGCACTCCGCGGTGTAAAGCGGGATCACATGGTCGCCAGGTTTGAc < 1:798058/149‑1 (MQ=255) | AATACCGGTGGTCACGCCACAGCCCAGCAGGCAGACGTGTTCATGGTTTGCTTCTGGATTAATTTTTGCCAGAGACACTTCCGCGACTACGGTGTATTCACTGAATGTAGAGCACCCCATGTAGTGATAAAGCGGCTGCCCGTTGTAAGAAAAACGGGTGGTGCCGTCTGGCATCAGGCCTTTACCCTGGGTTTCGCGAACCGCAACACAGAGGTTAGTTTTGCCAGAACGACAGAACTCACACTCGCCGCACTCCGCGGTGTAAAGCGGGATCACATGGTCGCCAGGTTTGAC > CP009273/374503‑374796 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 17 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |