Predicted mutation
evidence seq id position mutation freq annotation gene description
RA ACB122 4,135,612 T→G 33.4% W34G (TGG→GGG)  tgnC → (Z)‑2‑((N‑methylformamido)methylene)‑5‑hydroxybu tyrolactone dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*ACB1224,135,6120TG33.4% 30.5 / 4.1 44W34G (TGG→GGG) tgnC(Z)‑2‑((N‑methylformamido)methylene)‑5‑hydroxybu tyrolactone dehydrogenase
Reads supporting (aligned to +/- strand):  ref base T (10/16);  new base G (6/10);  total (18/26)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGCACCCTCACCGAATCCGCCAAGGGGATTGTTGCATGACCCTCGTTCGTTTCCAGATGTGCATCGACGGCCAATGGCGCGATGCCCAGAGCGGCAAGACCTTCGACAGCCTCAACCCGGCCACCGCCCAGGCCTGGGCGCAACTGCCCGACGCCGATGAAGCCGATGTGGAGCTGGCCGTTCAGGCCGCCCAGCGTGCCTTCGACAGCAAGGCGTGGCGCAGCATCACCGCTACCGCCCGCGGCAAGCTGCTGCGTCGTCTGG  >  ACB122/4135478‑4135741
                                                                                                                                      |                                                                                                                                 
cGCACCCTCACCGAATCCGCCAAGGGGATTGTTGCATGACCCTCGTTCGTTTCCAGATGTGCATCGACGGCCAATGGCGCGATGCCCAGAGCGGCAAGACCTTCGACAGCCTCAACCCGGCCACCGCCCAGGCCt                                                                                                                                   >  2:345488/1‑135 (MQ=255)
     ccTCACCGAATCCGCCAAGGGGATTGTTGCATGACCCTCGTTCGTTTCCAGATGTGCATCGACGGCCAATGGCGCGATGCCCAGAGCGGCAAGACCTTCGACAGCCTCAACCCGGCCACCGCCCAGGCCTGGgc                                                                                                                               >  4:109068/1‑134 (MQ=255)
        cACCGAATCCGCCAAGGGGATTGTTGCATGACCCTCGTTCGTTTCCAGATGTGCATCGACGGCCAATGGCGCGATGCCCAGAGCGGCAAGACCTTCGACAGCCTCAACCCGGCCACCGCCCACGCCGGGgcgcc                                                                                                                            >  8:112346/1‑133 (MQ=255)
          ccGAATCCGCCAAGGGGATTGTTGCATGACCCTCGTTCGTTTCCAGATGTGCATCGACGGCCAATGGCGCGATGCCCAGAGCGGCAAGACCTTCGACAGCCTCAACCCGGCCACCGCCCACGCCGGGGCGCCAcg                                                                                                                         >  5:165529/1‑134 (MQ=255)
            gAATCCGCCAAGGGGATTGTTGCATGACCCTCGTTCGTTTCCAGATGTGCATCGACGGCCAATGGCGCGATGCCCAGAGCGGCAAGACCTTCGACAGCCTCAACCCGGCCACCGCCCCCGCCGGGGCGCaccagc                                                                                                                       >  3:328606/1‑130 (MQ=255)
             aaTCGGCCAAGGGGATTGTTGCATGACCGTGGTTGTTTTCCAGATGTGAAACGACGGCCAATGGCGCGATGCCCAGAGCGGCAAGACCTTCGACAGCCTCAACCCGGCCACCGCCCCGGCGGGGGCGCAACTGcc                                                                                                                      <  8:246271/135‑1 (MQ=255)
               tCCGCCACGGGGTTTTTTGAATGACCGTCGTGCGTTTCCAGTTGTGCAGCGACGGCCAATGGGGCGATCCCCAGAGCGCCAAGACCTTCGACAGCCACACCCCGCCCCCCGCCCAGGCCGGGGCGCAACTGCCCg                                                                                                                    <  4:229113/135‑1 (MQ=255)
               gccGCCCAGGGGATTGTTGCATGACCCTCGTGCGTTGCCAGATGGGCATCGACGGCCAATGGCGCGAGGCCCAGAGCGCCAAGACCTTCGACAGCCTCACCCCGGCCACCGCCCAGGCCTGGGCGCAACTGCCCg                                                                                                                    <  3:324778/134‑1 (MQ=255)
                   ccAAGGGGATTGTTGCATGACCCTCGTTCGTTTCCAGATGTGCAGCGACGGCCAATGGCGCGATGCCCAGAGCGGCAAGACCTTCGACAGCCTCAACCCGGCCACCGCCCAGGCCTGGGCGCAACTGCCCGACGc                                                                                                                <  4:36011/135‑1 (MQ=255)
                    cAAGGGGATTGTTGCATGACCCTCGTTCGTTTCCAGATGTGCATCGACGGCCAATGGCGCGATGCCCAGAGCGGCAAGACCTTCGACAGCCTCAACCCGGCCACCGCCCCGGCCTGGGCGCAACTGCCCgccgcc                                                                                                               >  3:310725/1‑135 (MQ=255)
                      cggggATTGTTGCTGGCCCCTCGTTCGTTTCCAGATGTGCATCGCCGCCCACTGGCGCGATGCCCAGAGCGGCAAGACCTTCGACAGCCTCAACCCGGCCACCGCCCAGGCCTGGGCGCAACTGCCCGACGCCGa                                                                                                             <  8:327020/134‑1 (MQ=255)
                           aTTGTTGCATGACCCTCGTTCTTTTCCTCATGTGCATCTACGGCCTAGGGGGCGATGCCCAGAGCGGCAAGACCTTCGACAGCCTAAACCCGGCCACCGCCCAGGCCTGGGCGCAACTGCCCGACGCCGATGAAg                                                                                                        <  4:283858/135‑1 (MQ=255)
                              gttgCATGACCCTCGTTCGTTTCCAGATGTGCATCGACGGCCAATGGCGCGATGCCCAGAGCGGCAAGACCTTCGACAGCCTCAACCCGGCCACCGCCCCGGCCGGGGCGCCACTGCCCGCCGCCGCTGAAGCCg                                                                                                     >  5:152978/1‑135 (MQ=255)
                                  cATGACCCTCGTTCGTTTCCAGATGTGCATCGACGGCCAATGGCGCGATGCCCAGAGCGGCAAGACCTTCGACAGCCTCAACCCGGCCACCGCCCAGGCCTGGGCGCAACTGCCCGACGCCGATGAAGCCGAtgt                                                                                                 >  4:40792/1‑135 (MQ=255)
                                     gACCCTCGTTCGTTTCCAGATGTGCATCGACGGCCAATGGCGCGATGCCCAGAGCGGCAAGACCTTCGACAGCCTCACCCCGGCCACCGCCCCCGCCAGGGCGCAACAGCCCGACGCCGAGGAACCCGCAGTGGa                                                                                              >  8:290341/1‑135 (MQ=255)
                                       cccTCGTTCGTTTCCAGATGTGCATCGACGGCCAATGGCGCGATGCCCAGAGCGGCAAGACCTTCGACAGCCTCAACCCGGCCACCGCCCAGGCCTGGGCGCAACTGCCCGACGCCGATGAAGCCGATGTGGAg                                                                                             <  8:450377/134‑1 (MQ=255)
                                                  tgcccgATGTGCAGCGCCGGCCACTGGCGCGAGGCCCAGAGCGGCAAGACCTTCGACAGCCTCAACCCGGCCACCGCCCCGGCCGGGGCGCAACTGCCCGACGCCGATGAAGCCGATGTGGAGCTGGCCGTTCAg                                                                                 <  8:159193/130‑1 (MQ=255)
                                                    ccctttgtcCATCGCCGGCCAATGGCGCGAGGCCAAGAGCGCCAAGACCTTCGACCGCCCCACCCCGCCCACCGCCCCGGCCGGGGCGCAACTGCCCGACGCCGATGAAGCCGATGTGGAGCTGGCCGTTCAgg                                                                                <  3:109068/129‑1 (MQ=255)
                                                     cAGATGTGCATCGACGGCCAATGGCGCGATGCCCAGAGCGGCAAGACCTTCGACAGCCTCAACCCGGCCACCCCCCAGGCCAGGGCGCAACTGCCCCCCGCCGAATAAGACGAATTGGAGCTGGCCGGTAAGgcc                                                                              >  7:408462/1‑135 (MQ=255)
                                                             cATCGACGGCCAATGGCGCGATGCCCAGAGCGGCAAGACCTTCGACAGCCTCAACCCGGCCACCGCCCCGGCCGGGGCGGCCCTGCCCGCCGCCGCGGGAGCCGCGGTGGGGCTGGCCGGTGAGGCCGccccccg                                                                      >  7:271907/1‑131 (MQ=255)
                                                              aTCGACGGCCAATGGCGCGATGCCCAGAGCGGCAAGACCTTCGACAGCCTCAACCCGGCCACCGCCCCGGCCTGGGCGCCACGGCCCGCCGCCGAGGCAGCCGCTGTGGAGCTGGCCGGTCAGGCCGCCCCGCGt                                                                     >  1:86476/1‑135 (MQ=255)
                                                              aTCGACGGCCAATGGCGCGATGCCCAGAGCGGCAAGACCTTCGACAGCCTCAACCCGGCCACCGCCCAGGCCTGGGCGCAACTGCCCGCCGCCGCTGAAGCCGATGTGGAGCTGGCCGTTCACGCCGCCCAGCGt                                                                     >  7:155546/1‑135 (MQ=255)
                                                                cgccgcccACTGGCGCGATGCCCCGAGCGGCACCCCCTCCGACAGCCTCACCCCGGCCACCGCCCAGGCCTGGGCGCAACTGCCCGCCGCCGATGAAGCCGATGTGGAGCTGGCCGTTCAGGCCGCCCAGCGTGc                                                                   <  5:238837/129‑1 (MQ=255)
                                                                  aCGCCCACTGGCGCGAGGCCCCGAGCGGAACGCCCTTCGACCGCCTCACCCCGCCCCCCGCCCCGGCCTGGGCGCAACTGCCCGACGCCGATGAAGCCGATGTGGAGCTGGCCGTTCAGGCCGCCCAGCGTGc                                                                   <  2:172302/133‑1 (MQ=255)
                                                                   cGGCCAATGGCGCGATGCCCAGAGCGGCAAGACCTTCGACAGCCTCAACCCGGCCCCCGCCCCGGCCTGGGCGCCACTGCCCGCCGCCGCTGAAGCCGATGTGGAGCTGGCCGGTCAGGCCGCCCACCGGGccgt                                                                >  2:203105/1‑133 (MQ=255)
                                                                        aaTGGCGCGAGGCCCAGAGCGGCAAGACCTGCGACAGCCTCAACCCGGCCCCCGCCCAGGCCTGGGCGCAACTGCCCGACGCCGATGAAGCCGATGTGGAGCTGGCCGTTCAGGCCGCCCAGCGTGCCTTCGACa                                                           <  6:464864/135‑1 (MQ=255)
                                                                             cgcgATGCCCAGAGCGGCAAGACCTTCGACAGCCTCCACCCGGCCACCGCCCCGGCCGGGGCGCAACTGCCCGACGCCGATGAAGCCGATGTGGAGCTGGCCGTTCAGGCCGCCCAGCGTGCCTTCGACAGCAAg                                                      <  7:328316/135‑1 (MQ=255)
                                                                                  ggCCCAGCGCGGCAAGACCTTCGACAGCCGCACCCCGGCCCCCCCCCAGGCCTGGGCGCAACTGCCCGACGCCGATGAAGCCGATGTGGAGCTGGCCGTTCAGGCCGCCCAGCGTGCCTTCGACAGCAAGGCGTg                                                 <  6:192886/134‑1 (MQ=255)
                                                                                      cAGCGCCGCAACACCTTCGACAGCCTCACCCCGCCCCCCGCCCCGGCGGGGGCGCAACTGCCCGACGCCGATGAAGCCGATGTGGAGCTGGCCGTTCAGGCCGCCCAGCGTGCCTTCGACAGCAAGGCGTGgcgc                                             <  8:60432/135‑1 (MQ=255)
                                                                                      cAGAGCGGCAAGACCTGCGACAGCCTCACCCCGGCCACCGCCCCGGCCGGGGCGCAACTGCCCGACGCCGATGAAGCCGATGTGGAGCTGGCCGTTCAGGCCGCCCAGCGGGCCTTCGACAGCAAGGCGTGgcgc                                             <  6:170843/135‑1 (MQ=255)
                                                                                      cAGAGCGGCAAGACCTCCGACAGCCTCACCCCGGCCACCGCCCTGGCCGGGGCGCAACTGCCCGACGCCGATGAAGCCGATGTGGAGCTGGCCGTTCAGGCCGCCCAGCGTGCCTTCGACAGCAAGGCGTGgcgc                                             <  6:163055/135‑1 (MQ=255)
                                                                                           cGGCAAGACCTGCGACAGCCTCAACCCGGCCACCGCCCAGGCCTGGGCGCAACTGCCCGACGCCGATGAAGCCGATGTGGAGCTGGCCGTTCAGGCCGCCCAGCGTGCCTTCGACAGCAAGGCGTGGCGCAGCAt                                        <  4:328606/135‑1 (MQ=255)
                                                                                                agcccTTCGACAGCCCCACCCCGGCCCCCGCCCAGGCCTGGGCGCAACTGCCCGACGCCGATGAAGCCGATGTGGAGCTGGCCGTTCAGGCCGCCCAGCGTGCCTTCGACAGCAAGGCGTGGCGCAGCATCACCg                                   <  2:151982/132‑1 (MQ=255)
                                                                                                     ttCGACCGCCTCACCCCGCCCACCCCCCAGGCCTGGGCGCAACTGCCCGACGCCGATGAAGCCGATGTGGAGCTGGCCGTTCAGGCCGCCCAGCGTGCCTTCGACAGCAAGGCGTGGCGCAGCATCACCGCTAcc                              <  1:82660/135‑1 (MQ=255)
                                                                                                     ttCGACAGCCGCACCCCGGCCACCGCCCCGGCCGGGGCGCAACTGCCCGACGCCGATGAAGCCGATGTGGAGCTGGCCGTTCAGGCCGCCCAGCGTGCCTTCGACAGCAAGGCGTGGCGCAGCATCACCGCTAcc                              <  2:40996/135‑1 (MQ=255)
                                                                                                            gCCTCACCCCGGCCCCCGCCCCGGCCTGGGCGCAACTGCCCGACGCCGATGAAGCCGATGTGGAGCTGGCCGTTCAGGCCGCCCAGCGTGCCTTCGACAGCAAGGCGTGGCGCAGCATCACCGCTACCGCCcgcg                       <  5:287772/135‑1 (MQ=255)
                                                                                                            gCCTCACCCCGGCCCCCGCCCAGGCCTGGGCGCAACTGCCCGACGCCGATGAAGCCGATGTGGAGCTGGCCGTTCAGGCCGCCCAGCGTGCCTTCGACAGCAAGGCGTGGCGCAGCATCACCGCTACCGCCcgcg                       <  1:148946/135‑1 (MQ=255)
                                                                                                            cccccAACCCGGCCACCGCCCAGGCCTGGGCGCAACTGCCCGACGCCGATGAAGCCGATGTGGAGCTGGCCGTTCAGGCCGCCCAGCGTGCCTTCGACAGCAAGGCGTGGCGCAGCATCACCGCTACCGCCcgcg                       <  2:150318/131‑1 (MQ=255)
                                                                                                                cAACCCGGCCACCGCCCCGGCCTGGGCGCAACTGCCCGACGCCGATGAAGCCGCTGTGGAGCTGGCCGTTGAGGCCCCCCCGCGTGCCGTCGCCCGCAAGGCGTGGCGCAGCATCCCCGCGACCGCCCCCGGCaa                   >  8:189562/1‑135 (MQ=255)
                                                                                                                cAACCCGGCCACCGCCCAGGCCTGGGCGCAACTGCCCGCCGCCGATGAAGCCGATGTGGAGCTGGCCGTTCAGGCCGCCCAGCGTGCCTTCGACAGCAAGGCGTGGCGCAGCATCACCGCTACCGCCCGCGGCaa                   >  4:321367/1‑135 (MQ=255)
                                                                                                                          aCCGCCCAGGCCGGGGCGCAACTGCCCGACGCCGATGAAGCCGATGTGGAGCTGGCCGTTCAGGCCGCCCAGCGTGCCTTCGCCAGCAAGGCGTGGCGCAGCATCCCCGCTACCGCCCGCCGCCAGCTgcggcgg         >  4:41711/1‑134 (MQ=255)
                                                                                                                                cAGGCCTGGGCGCAACTGCCCGACGCCGATGAAGCCGATGTGGAGCTGGCCGTTCAGGCCGCCCAGCGTGCCTTCGACAGCAAGGCGTGGCGCAGCATCACCGCTACCGCCCGCGGCAAGCTGCTGCGTCGTCTg   <  7:343202/135‑1 (MQ=255)
                                                                                                                                cAGGCCGGGGCGCAACTGCCCGACGCCGATGAAGCCGATGTGGAGCTGGCCGTTCAGGCCGCCCAGCGTGCCTTCGACAGCAAGGCGTGGCGCAGCATCACCGCTACCGCCCGCGGCAAGCTGCTGCGTCGTCTg   <  5:70470/135‑1 (MQ=255)
                                                                                                                                 cggCCTGGGCGCAACTGCCCGCCGCCGCTGAAGCCGATGTGGAGCTGGCCGTTCAGGCCGCCCAGCGTGCCTTCGACAGCAAGGCGTGGCGCAGCATCACCGCTACCGCCCGCGGCAAGCTGCTGCGTCGTCTgg  >  5:374470/2‑135 (MQ=255)
                                                                                                                                     cTGGGCGCAACTGCCCGACGCCGATGAAGCCGATGTGGAGCTGGCCGTTCAGGCCGCCCAGCGTGCCTTCGACAGCAAGGCGTGGCGCAGCATCACCGCTACCGCCCGCGGCAAgctgct             >  8:392604/1‑120 (MQ=255)
                                                                                                                                      |                                                                                                                                 
CGCACCCTCACCGAATCCGCCAAGGGGATTGTTGCATGACCCTCGTTCGTTTCCAGATGTGCATCGACGGCCAATGGCGCGATGCCCAGAGCGGCAAGACCTTCGACAGCCTCAACCCGGCCACCGCCCAGGCCTGGGCGCAACTGCCCGACGCCGATGAAGCCGATGTGGAGCTGGCCGTTCAGGCCGCCCAGCGTGCCTTCGACAGCAAGGCGTGGCGCAGCATCACCGCTACCGCCCGCGGCAAGCTGCTGCGTCGTCTGG  >  ACB122/4135478‑4135741

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑