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
*CP009273374,6480T.25.5% 143.6 / 47.2 51coding (380/1110 nt)frmAalcohol 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
                                                                                                                                                 |                                                                                                                                                    
ccggtggtcacgccacagcccagcaggcagacgtgttcatggtttgcttctggattaatttttgccagagacacttccgcgactacggtgtattcactgaatgtagagcaccccatgtagtgataaagcggctgcccGCCCGTTgaagaa                                                                                                                                                 <  1:13526‑M2/13‑1 (MQ=255)
          cgccacagcccagcaggcagacgtgttcatggtttgcttccggattaatttttgccagagacacttccgcgactacggtgtattcactgaatgtagagcaccccatgtagtgataaagcggctgcccGCCCGTTGTAAGAAAAACGggtg                                                                                                                                        <  1:65120‑M2/23‑1 (MQ=255)
           gccacagcccagcaggcagacgtgttcatggtttgcttctggattaatttttgccagagacacttccgcgactacggtgtattcactgaatgtagagcaccccatgtagtgataaagcggctgcccGCCCGTTGTAAGAAAAACGggtgg                                                                                                                                       >  2:870375‑M2/127‑150 (MQ=255)
            ccacagcccagcaggcagacgtgttcatggtttgcttctggattaatttttgccagagacacttccgcgactacggtgtattcactgaatgtagagcaccccatgtagtgataaagcggctgcccGCCCGTTGAAGAAAAACGGgtggtg                                                                                                                                     <  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)
                             gacgtgttcatggtttgcttctggattaatttttgccagagacacttccgcgactacggtgtattcactgaatgtagagcaccccatgtagtgataaagcggctgcccGCCCGTTGAAGAAAAACGGGTGGTGCCGTCTGGCATCAGGcc                                                                                                                    >  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)
                                                    gattaatttttgccagagacacttccgcgactacggtgtattcactgaatgtagagcaccccatgtagtgataaagcggctgcccGCCCGTTGAAGAAAAACGGGTGGTGCCGTCTGGCATCAGGCCTTTACCCTGGGTTTCGCGAAcc                                                                                              <  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)
                                                                  agagacacttccgcgactacggtgtattcactgaatgtagagcaccccatgtagtgataaagcggctgcccGCCCGTTGAAGAAAAACGGGTGGTGCCGTCTGGCATCAGGCCTTTACCCTGGGTTTCGCGAACCGCAACACAGAGGt                                                                                 >  1:139326‑M2/72‑148 (MQ=255)
                                                                   gagacacttccgcgactacggtgtattcactgaatgtagagcaccccatgtagtgataaagcggctgcccGCCCGTTGTAAGAAAAACGGGTGGTGCCGTCTGGCATCAGGCCTTTACCCTGGGTTTCGCGAACCGCAACACAGAGGt                                                                                 >  2:798058‑M2/71‑148 (MQ=255)
                                                                     gacacttccgcgactacggtgtattcactgaatgtagagcaccccatgtagtgataaagcggctgcccGCCCGTTGAAGAAAAACGGGTGGTGCCGTCTGGCATCAGGCCTTTACCCTGGGTTTCGCGAACCGCAACACAGAGGTTAGtt                                                                            >  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)
                                                                                                     tgtagagcaccccatgtagtgataaagcggctgcccGCCCGTTGAAGAAAAACGGGTGGTGCCGTCTGGCATCAGGCCTTTACCCTGGGTTTCGCGAACCGCAACACAGAGGTTAGTTTTGCCAGAACGACAGAACTCACACTCGCCGCa                                            <  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)
                                                                                                                ccatgtagtgataaagcggctgcccGCCCGTTGAAGAAAAACGGGTGGTGCCGTCTGGCATCAGGCCTTTACCCTGGGTTTCGCGAACCGCAACACAGAGGTTAGTTTTGCCAGAACGACAGAACTCACACTCGCCGCACTCCGCGgtg                                  >  2:922192‑M2/26‑149 (MQ=255)
                                                                                                                  atgtagtgataaagcggctgcccGCCCGTTGTAAGAAAAACGGGTGGTGCCGTCTGGCATCAGGCCTTTACCCTGGGTTTCGCGAACCGCAACACAGAGGTTAGTTTTGCCAGAACGACAGAACTCACACTCGCCGCACTCCGCGgtgt                                 >  1:907951‑M2/24‑149 (MQ=255)
                                                                                                                     tagtgataaagcggctgcccGCCCGTTGTAAGAAAAACGGGTGGTGCCGTCTGGCATCAGGCCTTTACCCTGGGTTTCGCGAACCGCAACACAGAGGTTAGTTTTGCCAGAACGACAGAACTCACACTCGCCGCACTCCGCGGTGTAAAg                             <  2:1024869‑M2/130‑1 (MQ=255)
                                                                                                                     tagtgataaagcggctgcccGCCCGTTGAAGAAAAACGGGTGGTGCCGTCTGGCATCAGGCCTTTACCCTGGGTTTCGCGAACCGCAACACAGAGGTTAGTTTTGCCAGAACGACAGAACTCACACTCGCCGCACTCCGCGGTGTAAAGc                            <  2:480290‑M2/130‑1 (MQ=255)
                                                                                                                         gataaagcggctgcccGCCCGTTGAAGAAAAACGGGTGGTGCCGTCTGGCATCAGGCCTTTACCCTGGGTTTCGCGAACCGCAACACAGAGGTTAGTTTTGCCAGAACGACAGAACTCACACTCGCCGCACTCCGCGGTGTAAAGCGGGa                        >  1:619163‑M2/17‑150 (MQ=255)
                                                                                                                          ataaagcggctgcccGCCCGTTGAAGAAAAACGGGTGGTGCCGTCTGGCATCAGGCCTTTACCCTGGGTTTCGCGAACCGCAACACAGAGGTTAGTTTTGCCAGAACGACAGAACTCACACTCGCCGCACTCCGCGGTGTAAAGCGGGa                        <  1:113029‑M2/134‑1 (MQ=255)
                                                                                                                          ataaagcggctgcccGCCCGTTGAAGAAAAACGGGTGGTGCCGTCTGGCATCAGGCCTTTACCCTGGGTTTCGCGAACCGCAACACAGAGGTTAGTTTTGCCAGAACGACAGAACTCACACTCGCCGCACTCCGCGGTGTAAAGCGGGAt                       <  2:560106‑M2/135‑1 (MQ=255)
                                                                                                                            aaagcggctgcccGCCCGTTGTAAGAAAAACGGGTGGTGCCGTCTGGCATCAGGCCTTTACCCTGGGTTTCGCGAACCGCAACACAGAGGTTAGTTTTGCCAGAACGACAGAACTCACACTCGCCGCACTCCGCGGTGTAAAGCGGGAt                       >  1:993473‑M2/14‑149 (MQ=255)
                                                                                                                                    tgcccGCCCGTTGAAGAAAAACGGGTGGTGCCGTCTGGCATCAGGCCTTTACCCTGGGTTTCGCGAACCGCAACACAGAGGTTAGTTTTGCCAGTACGACAGAACTCACACTCGCCGCACTCCGCGGTGTAAAGCGGGATCACATGGTCg             >  2:867632‑M2/6‑150 (MQ=255)
                                                                                                                                     gcccGCCCGTTGTAAGAAAAACGGGTGGTGCCGTCTGGCATCAGGCCTTTACCCTGGGTTTCGCGAACCGCAACACAGAGGTTAGTTTTGCCAGAACGACAGAACTCACACTCGCCGCACTCCGCGGTGTAAAGCGGGATCACATGGTCg             >  1:281704‑M2/5‑150 (MQ=255)
                                                                                                                                      cccGCCCGTTGTAAGAAAAACGGGTGGTGCCGTCTGGCATCAGGCCTTTACCCTGGGTTTCGCGAACCGCAACACAGAGGTTAGTTTTGCCAGAACGACAGAACTCACACTCGCCGCACTCCGCGGTGTAAAGCGGGATCACATGGTCg             <  2:993473‑M2/146‑1 (MQ=255)
                                                                                                                                        cGCCCGTTGAAGAAAAACGGGTGGTGCCGTCTGGCATCAGGCCTTTACCCTGGGTTTCGCGAACCGCAACACAGAGGTTAGTTTTGCCAGAACGACAGAACTCACACTCGCCGCACTCCGCGGTGTAAAGCGGGATCACATGGTCGCCAg         <  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: