Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 374,648 Δ1 bp 54.4% 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.54.4% 99.9 / 142.1 57coding (380/1110 nt)frmAalcohol dehydrogenase class III/glutathione‑dependent formaldehyde dehydrogenase
Reads supporting (aligned to +/- strand):  ref base T (17/9);  new base . (14/17);  total (31/26)
Fisher's exact test for biased strand distribution p-value = 1.83e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCCAATACCGGTGGTCACGCCACAGCCCAGCAGGCAGACGTGTTCATGGTTTGCTTCTGGATTAATTTTTGCCAGAGACACTTCCGCGACTACGGTGTATTCACTGAATGTAGAGCACCCCATGTAGTGATAAAGCGGCTGCCCGTTGTAAGAAAAACGGGTGGTGCCGTCTGGCATCAGGCCTTTACCCTGGGTTTCGCGAACCGCAACACAGAGGTTAGTTTTGCCAGAACGACAGAACTCACACTCGCCGCACTCCGCGGTGTAAAGCGGGATCACATGGTCGCCAGGTTTGAC  >  CP009273/374500‑374796
                                                                                                                                                    |                                                                                                                                                    
ataccggtggtcacgccacagcccagcaggcagacgtgttcatggtttgcttctggattaatttttgccagagacacttccgcgactacggtgtattcactgaatgtagagcaccccatgtagtgataaagcggctgcccGCCCGTTGt                                                                                                                                                      >  2:965888‑M2/141‑149 (MQ=255)
            acgccacagcccagcaggcagacgtgttcatggtttgcttctggattaatttttgccagagacacttccgcgactacggtgtattcactgaatgtagagcaccccatgtagtgataaagcggctgcccGCCCGTTGTAAGAAAAACGggt                                                                                                                                         <  1:479518‑M2/22‑1 (MQ=255)
            acgccacagcccagcaggcagacgtgttcatggtttgcttctggattaatttttgccagagacacttccgcgactacggtgtattcactgaatgtagagcaccccatgtagtgataaagcggctgcccGCCCGTTGAAGAAAAACGggtg                                                                                                                                        <  2:880208‑M2/22‑1 (MQ=255)
                cacagcccagcaggcagacgtgttcatggtttgcttctggattaatttttgccagagacacttccgcgactacggtgtattcactgaatgtagagcaccccatgtagtgataaagcggctgcccGCCCGTTGAAGAAAAACGGGTGGTGc                                                                                                                                    <  2:952949‑M2/26‑1 (MQ=255)
                   agcccagcaggcagacgtgttcatggtttgcttctggattaatttttgccagagacacttccgcgactacggtgtattcactgaatgtagagcaccccatgtagtgataaagcggctgcccGCCCGTTGTAAGAAAAACGGGTGGTGCCg                                                                                                                                  >  2:718842‑M2/122‑150 (MQ=255)
                     cccagcaggcagacgtgttcatggtttgcttctggattaatttttgccagagacacttccgcgactacggtgtattcactgaatgtagagcaccccatgtagtgataaagcggctgcccGCCCGTTGAAGAAAAACGGGTGGTGCCGTCt                                                                                                                               <  2:718781‑M2/31‑1 (MQ=255)
                       cagcaggcagacgtgttcatggtttgcttctggattaatttttgccagagacacttccgcgactacggtgtattcactgaatgtagagcaccccatgtagtgataaagcggctgcccGCCCGTTGTAAGAAAAACGGGTGGTGCCGTCTg                                                                                                                              >  2:231057‑M2/118‑150 (MQ=255)
                           aggcagacgtgttcatggtttgcttctggattaatttttgccagagacacttccgcgactacggtgtattcactgaatgtagagcaccccatgtagtgataaagcggctgcccGCCCGTTGAAGAAAAACGGGTGGTGCCGTCTGGCATc                                                                                                                         <  1:401458‑M2/37‑1 (MQ=255)
                                      ttcatggtttgcttctggattaatttttgccagagacacttccgcgactacggtgtattcactgaatgtagagcaccccatgtagtgataaagcggctgcccGCCCGTTGAAGAAAAACGGGTGGTGCCGTCTGGCATCAGGCCTTTAcc                                                                                                              >  1:391837‑M2/103‑150 (MQ=255)
                                         atggtttgcttctggattaatttttgccagagacacttccgcgactacggtgtattcactgaatgtagagcaccccatgtagtgataaagcggctgcccGCCCGTTGTAAGAAAAACGGGTGGTGCCGTCTGGCATCAGGCCTTTAccc                                                                                                             <  2:221591‑M2/50‑1 (MQ=255)
                                          tggtttgcttctggattaatttttgccagagacacttccgcgactacggtgtattcactgaatgtagagcaccccatgtagtgataaagcggctgcccGCCCGTTGAAGAAAAACGGGTGGTGCCGTCTGGCATCAGGCCTTTACCCTgg                                                                                                          <  1:405400‑M2/52‑1 (MQ=255)
                                                  ttctggattaatttttgccagagtcacttccgcgactacggtgtattcactgaatgtagagcaccccatgtagtgatatagcggctgcccGCCCGTTGTAAGAAAAACGGGTGGTGCCGTCTGGCATCAGGCCTTTACCCTGGGTTTcgc                                                                                                   >  2:126656‑M2/91‑150 (MQ=255)
                                                  ttctggattaatttttgccagagacacttccgcgactacggtgtattcactgaatgtagagcaccccatgtagtgataaagcggctgcccGCCCGTTGAAGAAAAACGGGTGGTGCCGTCTGGCATCAGGCCTTTACCCTGGGTTTcgc                                                                                                   >  2:805110‑M2/91‑149 (MQ=255)
                                                  ttctggattaatttttgccagagacacttccgcgactacggtgtattcactgaatgtagagcaccccatgtagtgataaagcggctgcccGCCCGTTGAAGAAAAACGGGTGGTGCCGTCTGGCATCAGGCCTTTACCCTGGGTTTcgc                                                                                                   >  2:147378‑M2/91‑149 (MQ=255)
                                                       gattaatttttgccagagacacttccgcgactacggtgtattcactgaatgtagagcaccccatgtagtgataaagcggctgcccGCCCGTTGTAAGAAAAACGGGTGGTGCCGTCTGGCATCAGGCCTTTACCCTGGGTTTCGCGAAcc                                                                                              >  2:654997‑M2/86‑150 (MQ=255)
                                                           aatttttgccagagacacttccgcgactacggtgtattcactgaatgtagagcaccccatgtagtgataaagcggctgcccGCCCGTTGTAAGAAAAACGGGTGGTGCCGTCTGGCATCAGGCCTTTACCCTGGGTTTCGCGAACCGCaa                                                                                          >  2:701094‑M2/82‑150 (MQ=255)
                                                           aatttttgccagagacacttccgcgactacggtgtattcactgaatgtagagcaccccatgtagtgataaagcggctgcccGCCCGTTGAAGAAAAACGGGTGGTGCCGTCTGGCATCAGGCCTTTACCCTGGGTTTCGCGAACCGCaa                                                                                          >  2:755617‑M2/82‑149 (MQ=255)
                                                           aatttttgccagagacacttccgcgactacggtgtattcactgaatgtagagcaccccatgtagtgataaagcggctgcccGCCCGTTGAAGAAAAACGGGTGGTGCCGTCTGGCATCAGGCCTTTACCCTGGGTTTCGCGAACCGCaa                                                                                          >  2:415382‑M2/82‑149 (MQ=255)
                                                                ttgccagagacacttccgcgactacggtgtattcactgaatgtagagcaccccatgtagtgataaagcggctgcccGCCCGTTGAAGAAAAACGGGTGGTGCCGTCTGGCATCAGGCCTTTACCCTGGGTTTCGCGAACCGCAACACag                                                                                     <  2:102075‑M2/73‑1 (MQ=255)
                                                                 tgccagagacacttccgcgactacggtgtattcactgaatgtagagcaccccatgtagtgataaagcggctgcccGCCCGTTGAAGAAAAACGGGTGGTGCCGTCTGGCATCAGGCCTTTACCCTGGGTTTCGCGAACCGCAACACagag                                                                                   >  1:959947‑M2/76‑150 (MQ=255)
                                                                      gagacacttccgcgactacggtgtattcactgaatgtagagcaccccatgtagtgataaagcggctgcccGCCCGTTGTAAGAAAAACGGGTGGTGCCGTCTGGCATCAGGCCTTTACCCTGGGTTTCGCGAACCGCAACACAGAGGTTa                                                                               <  2:490826‑M2/80‑1 (MQ=255)
                                                                       agacacttccgcgactacggtgtattcactgaatgtagagcaccccatgtagtgataaagcggctgcccGCCCGTTGAAGAAAAACGGGTGGTGCCGTCTGGCATCAGGCCTTTACCCTGGGTTTCGCGAACCGCAACACAGAGGTTAGt                                                                             >  2:759229‑M2/70‑150 (MQ=255)
                                                                        gacacttccgcgactacggtgtattcactgaatgtagagcaccccatgtagtgataaagcggctgcccGCCCGTTGTAAGAAAAACGGGTGGTGCCGTCTGGCATCAGGCCTTTACCCTGGGTTTCGCGAACCGCAACACAGAGGTTAGt                                                                             >  1:766987‑M2/69‑150 (MQ=255)
                                                                        gacacttccgcgactacggtgtattcactgaatgtagagcaccccatgtagtgataaagcggctgcccGCCCGTTGTAAGAAAAACGGGTGGTGCCGTCTGGCATCAGGCCTTTACCCTGGGTTTCGCGAACCGCAACACAGAGGTTAGt                                                                             >  2:653803‑M2/69‑150 (MQ=255)
                                                                          cacttccgcgactacggtgtattcactgaatgtagagcaccccatgtagtgataaagcggctgcccGCCCGTTGAAGAAAAACGGGTGGTGCCGTCTGGCATCAGGCCTTTACCCTGGGTTTCGCGAACCGCAACACAGAGGTTAGtttt                                                                          <  1:805110‑M2/84‑1 (MQ=255)
                                                                           acttccgcgactacggtgtattcactgaatgtagagcaccccatgtagtgataaagcggctgcccGCCCGTTGTAAGAAAAACGGGTGGTGCCGTCTGGCATCAGGCCTTTACCCTGGGTTTCGCGAACCGCAACACAGAGGTTAGtttt                                                                          <  2:246956‑M2/85‑1 (MQ=255)
                                                                           acttccgcgactacggtgtattcactgaatgtagagcaccccatgtagtgataaagcggctgcccGCCCGTTGAAGAAAAACGGGTGGTGCCGTCTGGCATCAGGCCTTTACCCTGGGTTTCGCGAACCGCAACACAGAGGTTAGTTTTg                                                                         >  1:976522‑M2/66‑150 (MQ=255)
                                                                             ttccgcgactacggtgtattcactgaatgtagagcaccccatgtagtgataaagcggctgcccGCCCGTTGAAGAAAAACGGGTGGTGCCGTCTGGCATCAGGCCTTTACCCTGGGTTTCGCGAACCGCAACACAGAGGTTAGTTTTGc                                                                        >  2:983699‑M2/64‑149 (MQ=255)
                                                                               ccgcgactacggtgtattcactgaatgtagagcaccccatgtagtgataaagcggctgcccGCCCGTTGAAGAAAAACGGGTGGTGCCGTCTGGCATCAGGCCTTTACCCTGGGTTTCGCGAACCGCAACACAGAGGTTAGTTTTGCCa                                                                      <  2:97764‑M2/88‑1 (MQ=255)
                                                                                cgcgactacggtgtattcactgaatgtagagcaccccatgtagtgataaagcggctgcccGCCCGTTGAAGAAAAACGGGTGGTGCCGTCTGGCATCAGGCCTTTACCCTGGGTTTCGCGAACCGCAACACAGAGGTTAGTTTTGCCAg                                                                     <  1:714435‑M2/89‑1 (MQ=255)
                                                                                   gactacggtgtattcactgaatgtagagcaccccatgtagtgataaagcggctgcccGCCCGTTGTAAGAAAAACGGGTGGTGCCGTCTGGCATCAGGCCTTTACCCTGGGTTTCGCGAACCGCAACACAGAGGTTAGTTTTGCCAGAAc                                                                  <  1:890655‑M2/93‑1 (MQ=255)
                                                                                         ggtgtattcactgaatgtagagcaccccatgtagtgataaagcggctgcccGCCCGTTGTAAGAAAAACGGGTGGTGCCGTCTGGCATCAGGCCTTTACCCTGGGTTTCGCGAACCGCAACACAGAGGTTAGTTTTGCCAGAACGACAg                                                             >  1:576550‑M2/52‑149 (MQ=255)
                                                                                         ggtgtattcactgaatgtagagcaccccatgtagtgataaagcggctgcccGCCCGTTGAAGAAAAACGGGTGGTGCCGTCTGGCATCAGGCCTTTACCCTGGGTTTCGCGAACCGCAACACAGAGGTTAGTTTTGCCAGAACGACAGaa                                                           >  2:671983‑M2/52‑150 (MQ=255)
                                                                                               ttcactgaatgtagagcaccccatgtagtgataaagcggctgcccGCCCGTTGTAAGAAAAACGGGTGGTGCCGTCTGGCATCAGGCCTTTACCCTGGGTTTCGCGAACCGCAACACAGAGGTTAGTTTTGCCAGAACGACAGAACTca                                                       >  1:101797‑M2/46‑149 (MQ=255)
                                                                                               ttcactgaatgtagagcaccccatgtagtgataaagcggctgcccGCCCGTTGAAGAAAAACGGGTGGTGCCGTCTGGCATCAGGCCTTTACCCTGGGTTTCGCGAACCGCAACACAGAGGTTAGTTTTGCCAGAACGACAGAACTcac                                                      <  2:208830‑M2/104‑1 (MQ=255)
                                                                                                 cactgaatgtagagcaccccatgtagtgataaagcggctgcccGCCCGTTGTAAGAAAAACGGGTGGTGCCGTCTGGCATCAGGCCTTTACCCTGGGTTTCGCGAACCGCAACACAGAGGTTAGTTTTGCCAGAACGACAGAACTCacac                                                    >  2:129864‑M2/44‑150 (MQ=255)
                                                                                                 cactgaatgtagagcaccccatgtagtgataaagcggctgcccGCCCGTTGTAAGAAAAACGGGTGGTGCCGTCTGGCATCAGGCCTTTACCCTGGGTTTCGCGAACCGCAACACAGAGGTTAGTTTTGCCAGAACGACAGAACTCacac                                                    >  2:750085‑M2/44‑150 (MQ=255)
                                                                                                      aatgtagagcaccccatgtagtgataaagcggctgcccGCCCGTTGAAGAAAAACGGGTGGTGCCGTCTGGCATCAGGCCTTTACCCTGGGTTTCGCGAACCGCAACACAGAGGTTAGTTTTGCCAGAACGACAGAACTCACACTcgc                                                <  2:866699‑M2/110‑1 (MQ=255)
                                                                                                          tagagcaccccatgtagtgataaagcggctgcccGCCCGTTGTAAGAAAAACGGGTGGTGCCGTCTGGCATCAGGCCTTTACCCTGGGTTTCGCGAACCGCAACACAGAGGTTAGTTTTGCCAGAACGACAGAACTCAAACTCGCCGCAc                                           <  1:965888‑M2/116‑1 (MQ=255)
                                                                                                            gagcaccccatgtagtgataaagcggctgcccGCCCGTTGAAGAAAAACGGGTGGTGCCGTCTGGCATCAGGCCTTTACCCTGGGTTTCGCGAACCGCAACACAGAGGTTAGTTTTGCCAGAACGACAGAACTCACACTCGCCGCACTcc                                        <  1:415814‑M2/118‑1 (MQ=255)
                                                                                                              gcaccccatgtagtgataaagcggctgcccGCCCGTTGAAGAAAAACGGGTGGTGCCGTCTGGCATCAGGCCTTTACCCTGGGTTTCGCGAACCGCAACACAGAGGTTAGTTTTGCCAGAACGACAGAACTCACACTCGCCGCACTCcgc                                      >  2:520462‑M2/31‑150 (MQ=255)
                                                                                                               caccccatgtagtgataaagcggctgcccGCCCGTTGTAAGAAAAACGGGTGGTGCCGTCTGGCATCAGGCCTTTACCCTGGGTTTCGCGAACCGCAACACAGAGGTTAGTTTTGCCAGAACGACAGAACTCACACTCGCCGCACTCcgc                                      >  2:214183‑M2/30‑150 (MQ=255)
                                                                                                                accccatgtagtgataaagcggctgcccGCCCGTTGTAAGAAAAACGGGTGGTGCCGTCTGGCATCAGGCCTTTACCCTGGGTTTCGCGAACCGCAACACAGAGGTTAGTTTTGCCAGAACGACAGAACTCACACTCGCCGCACTCcgcg                                     <  2:904957‑M2/122‑1 (MQ=255)
                                                                                                                accccatgtagtgataaagcggctgcccGCCCGTTGTAAGAAAAACGGGTGGTGCCGTCTGGCATCAGGCCTTTACCCTGGGTTTCGCGAACCGCAACACAGAGGTTAGTTTTGCCAGAACGACAGAACTCACACTCGCCGCACTCcgcg                                     >  1:676200‑M2/29‑150 (MQ=255)
                                                                                                                 ccccatgtagtgataaagcggctgcccGCCCGTTGAAGAAAAACGGGTGGTGCCGTCTGGCATCAGGCCTTTACCCTGGGTTTCGCGAACCGCAACACAGAGGTTAGTTTTGCCAGAACGACAGAACTCACACTCGCCGCACTCCGCGgt                                   <  1:643596‑M2/123‑1 (MQ=255)
                                                                                                                          gtgataaagcggctgcccGCCCGTTGTAAGAAAAACGGGTGGTGCCGTCTGGCATCAGGCCTTTACCCTGGGTTTCGCGAACCGCAACACAGAGGTTAGTTTTGCCAGAACGACAGAACTCACACTCGCCGCACTCCGCGGTGTAAAGCg                           >  1:122623‑M2/19‑150 (MQ=255)
                                                                                                                          gtgataaagcggctgcccGCCCGTTGTAAGAAAAACGGGTGGTGCCGTCTGGCATCAGGCCTTTACCCTGGGTTTCGCGAACCGCAACACAGAGGTTAGTTTTGCCAGAACGACAGAACTCACACTCGCCGCACTCCGCGGTGTAAAGCg                           >  2:812397‑M2/19‑150 (MQ=255)
                                                                                                                            gataaagcggctgcccGCCCGTTGAAGAAAAACGGGTGGTGCCGTCTGGCATCAGGCCTTTACCCTGGGTTTCGCGAACCGCAACACAGAGGTTAGTTTTGCCAGAACGACAGAACTCACACTCGCCGCACTCCGCGGTGTAAAGCGGGa                        >  1:813000‑M2/17‑150 (MQ=255)
                                                                                                                              taaagcggctgcccGCCCGTTGAAGAAAAACGGGTGGTGCCGTCTGGCATCAGGCCTTTACCCTGGGTTTCGCGAACCGCAACACAGAGGTTAGTTTTGCCAGAACGACAGAACTCACACTCGCCGCACTCCGCGGTGTAAAGCGGGATc                      <  1:415382‑M2/136‑1 (MQ=255)
                                                                                                                                aagcggctgcccGCCCGTTGAAGAAAAACGGGTGGTGCCGTCTGGCATCAGGCCTTTACCCTGGGTTTCGCGAACCGCAACACAGAGGTTAGTTTTGCCAGAACGACAGAACTCACACTCGCCGCACTCCGCGGTGTAAAGCGGGATcac                    <  1:48583‑M2/138‑1 (MQ=255)
                                                                                                                                aagcggctgcccGCCCGTTGAAGAAAAACGGGTGGTGCCGTCTGGCATCAGGCCTTTACCCTGGGTTTCGCGAACCGCAACACAGAGGTTAGTTTTGCCAGAACGACAGAACTCACACTCGCCGCACTCCGCGGTGTAAAGCGGGATcac                    <  2:976522‑M2/138‑1 (MQ=255)
                                                                                                                                aagcggctgcccGCCCGTTGAAGAAAAACGGGTGGTGCCGTCTGGCATCAGGCCTTTACCCTGGGTTTCGCGAACCGCAACACAGAGGTTAGTTTTGCCAGAACGACAGAACTCACACTCGCCGCACTCCGCGGTGTAAAGCGGGATca                     >  1:930013‑M2/13‑149 (MQ=255)
                                                                                                                                  gcggctgcccGCCCGTTGTAAGAAAAACGGGTGGTGCCGTCTGGCATCAGGCCTTTACCCTGGGTTTCGCGAACCGCAACACAGAGGTTAGTTTTGCCAGAACGACAGAACTCACACTCGCCGCACTCCGCGGTGTAAAGCGGGATcaca                   <  2:743749‑M2/140‑1 (MQ=255)
                                                                                                                                  gcggctgcccGCCCGTTGTAAGAAAAACGGGTGGTGCCGTCTGGCATCAGGCCTTTACCCTGGGTTTCGCGAACCGCAACACAGAGGTTAGTTTTGCCAGAACGACAGAACTCACACTCGCCGCACTCCGCGGTGTAAAGCGGGATcaca                   <  2:101797‑M2/140‑1 (MQ=255)
                                                                                                                                   cggctgcccGCCCGTTGAAGAAAAACGGGTGGTGCCGTCTGGCATCAGGCCTTTACCCTGGGTTTCGCGAACCGCAACACAGAGGTTAGTTTTGCCAGAACGACAGAACTCACACTCGCCGCACTCCGCGGTGTAAAGCGGGATCACATg                 >  1:259242‑M2/10‑150 (MQ=255)
                                                                                                                                       tgcccGCCCGTTGTAAGAAAAACGGGTGGTGCCGTCTGGCATCAGGCCTTTACCCTGGGTTTCGCGAACCGCAACACAGAGGTTAGTTTTGCCAGAACGACAGAACTCACACTCGCCGCACTCCGCGGTGTAAAGCGGGATCACATGGTc              >  1:583459‑M2/6‑150 (MQ=255)
                                                                                                                                         cccGCCCGTTGTAAGAAAAACGGGTGGTGCCGTCTGGCATCAGGCCTTTACCCTGGGTTTCGCGAACCGCAACACAGAGGTTAGTTTTGCCAGAACGACAGAACTCACACTCGCCGCACTCCGCGGTGTAAAGCGGGATCACATGGTCGc            >  1:113801‑M2/4‑150 (MQ=255)
                                                                                                                                              ccGTTGAAGAAAAACGGGTGGTGCCGTCTGGCATCAGGCCTTTACCCTGGGTTTCGCGAACCGCAACACAGAGGTTAGTTTTGCCAGAACGACAGAACTCACACTCGCCGCACTCCGCGGTGTAAAGCGGGATCACATGGTCGCCAGGtt      <  1:983699/150‑1 (MQ=255)
                                                                                                                                                    tAAGAAAAACGGGTGGTGCCGTCTGGCATCAGGCCTTTACCCTGGGTTTCGCGAACCGCAACACAGAGGTTAGTTTTGCCAGAACGACAGAACTCACACTCGCCGCACTCCGCGGTGTAAAGCGGGATCACATGGTCGCCAGGTTTGAc  <  1:129864/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
GCCAATACCGGTGGTCACGCCACAGCCCAGCAGGCAGACGTGTTCATGGTTTGCTTCTGGATTAATTTTTGCCAGAGACACTTCCGCGACTACGGTGTATTCACTGAATGTAGAGCACCCCATGTAGTGATAAAGCGGCTGCCCGTTGTAAGAAAAACGGGTGGTGCCGTCTGGCATCAGGCCTTTACCCTGGGTTTCGCGAACCGCAACACAGAGGTTAGTTTTGCCAGAACGACAGAACTCACACTCGCCGCACTCCGCGGTGTAAAGCGGGATCACATGGTCGCCAGGTTTGAC  >  CP009273/374500‑374796

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 29 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: