Predicted mutation
evidence seq id position mutation freq annotation gene description
RA KT2440 4,135,612 T→G 46.2% 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
*KT24404,135,6120TG46.2% 3.5 / 9.1 38W34G (TGG→GGG) tgnC(Z)‑2‑((N‑methylformamido)methylene)‑5‑hydroxybu tyrolactone dehydrogenase
Reads supporting (aligned to +/- strand):  ref base T (9/7);  new base G (7/9);  total (22/16)
Fisher's exact test for biased strand distribution p-value = 7.24e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCTCACCGAATCCGCCAAGGGGATTGTTGCATGACCCTCGTTCGTTTCCAGATGTGCATCGACGGCCAATGGCGCGATGCCCAGAGCGGCAAGACCTTCGACAGCCTCAACCCGGCCACCGCCCAGGCCTGGGCGCAACTGCCCGACGCCGATGAAGCCGATGTGGAGCTGGCCGTTCAGGCCGCCCAGCGTGCCTTCGACAGCAAGGCGTGGCGCAGCATCACCGCTACCGCCCGCGGCAAGCTGCTGCGTCGTC  >  KT2440/4135483‑4135738
                                                                                                                                 |                                                                                                                              
ccTCACCGAATCCGCCAAGGGGATTGTTGCATGACCCTCGTTCGTTTCCAGATGTGCATCGACGGCCAATGGCGCGATGCCCAGAGCGGCAAGACCTTCGACAGCCTCAACCCGGCCCCCGCCCACCCCGGGgccc                                                                                                                          >  8:307441/1‑134 (MQ=255)
 cTCGCCGAATCCGCCAAGGGGATTGTTGCATGACCCTCGATCGTTTCCAGATGTGCATCGACGGCCAAAGGCGCGATGCCCAGCGCGTCTAGACCTTCGACCGCCGCCACCCGGCCACCCCCCCGGCCAggggcc                                                                                                                          >  2:180975/1‑132 (MQ=255)
     ccGAATCCGCCAAGGGGATTGTTGCATGACCCTCGTTCGTTTCCAGATGTGCATCGACGGCCAATGGCGCGATGCCCAGAGCGGCAAGACCTTCGACAGCCTCAACCCGGCCACCGCCCAGGCCTGGGCGCAACt                                                                                                                      >  4:95317/1‑135 (MQ=255)
           ccGCCAAGGGGATTGTTGCATGACCCTCGTTCGTTTCCAGATGTGCATCGACGGCCAATGGCGCGATGCCCAGAGCGGCAAGACCTTCGACAGCCTCCACCCGGCCACCGCCCCCCCCCGGGCGCCCCAGCCCgc                                                                                                                >  1:170712/1‑134 (MQ=255)
                 aGGGGATTGTTGCATGACCCTCGTTCGTTTCCAGATGTGCATCGACGGCCAATGGCGCGATGCCCAGAGCGGCAAGACCTTCGACAGCCTCAACCCGGCCACCGCCCACGCCCGGGCGCACCTGCCCGCCGCCGAt                                                                                                         >  7:208184/1‑136 (MQ=255)
                    ggATTGTTGCATGACCCTCGTTCGTTTCCAGATGTGCATCGACGGCCAATGGCGCGATGCCCAGAGCGGCAAGACCTTCGACAGCCTCAACCCGGCCACCGCCCAGGCCTGGGCGCAACTGCCCGACGCCGATGa                                                                                                       <  2:292369/135‑1 (MQ=255)
                    ggATGGTTGCTTGACCCTCGTGCGTTCCCATATGTGCAACGACGGCCAATGGCGCGTGGCCCAGAGCGGCAAGACCTTCGACAGCCTCAACCCGGCCACCGCCCAGGCCGGGGCGCAACTGCCCGACGCCGATGa                                                                                                       <  6:1943/135‑1 (MQ=255)
                         gttgCATGACCCTCGTTCGTTTCCAGATGTGCATCGACGGCCAATGGCGCGATGCCCAGAGCGGCAAGACCTTCGACAGCCTCAACCCGGCCACCGCCCACGCCAGGGCGCCACTGCCCGCCCCCGATGCAGcc                                                                                                   >  4:438101/1‑134 (MQ=255)
                           tgCATGACCCTCGTTCGTTTCCAGATGTGCATCGACGGCCAATGGCGCGATGCCCAGAGCGGCAAGACCTTCGACAGCCTCAACCCGGCCACCGCCCAGGCCTGGGCGCAACAGCCCGCCGCCGATGAATcccaag                                                                                               >  7:147755/1‑132 (MQ=255)
                                gACCCTCGTTCGTTTCCAGATGTGCATCGACGGCCAATGGCGCGATGCCCAGAGCGGCAAGACCTTCGACAGCCTCAACCCGGCCACCGCCCCCCCCGGGGCGCCCCTCCCCCCCGCCCACCAAGCCCAAGTGGAg                                                                                          >  7:80024/1‑136 (MQ=255)
                                       gttcgttTCCAGATGTGCATCGACGGCCAATGGCGCGATGCCCAGAGCGGCAAGACCTTCGCCAGCCTCAACCCGGCCACCGCCCCGGCCGGGGCGCCACGGCCCGCCGCCGAGGAAGCCGCTGTGGGGCTGGcc                                                                                    >  1:91552/1‑135 (MQ=255)
                                               ccAGATGTGCATCGACGGCCAATGGCGCGATGCCCAGAGCGGCAAGACCTTCGACAGCCTCACCCCGGCCACCGCCCAGGCCTGGGCGCCACTGCCCGCCGCCGCTGCAGCCGATGTGGCGCTGGCCGGTGAGgcc                                                                           >  8:446469/1‑136 (MQ=255)
                                               ccAGATGTGCATCGACGGCCAATGGCGCGATGCCCAGAGCGGCAAGACCTTCGACAGCCTCAACCCGGCCACCGCCCCGGCCTGGGCGCCACAGCCCGCCGCCCCTGAAGCCCCTGTGGAGCTGGCCGTTCAGgc                                                                            >  2:6971/1‑135 (MQ=255)
                                               ccAGATGTGCATCGACGGCCAATGGCGCGATGCCCAGAGCGGCAAGACCTTCGACAGCCTCAACCCGGCCACCGCCCAGGCCGGGGCGCAACAGCCCGACGCCGATGAAGCCGCTGTGGAGCTGGCCGTTCAGgc                                                                            >  6:339277/1‑135 (MQ=255)
                                               ccAGATGTGCATCGACGGCCAATGGCGCGATGCCCAGAGCGGCAAGACCTTCGACAGCCTCAACCCGGCCACCGCCCACGCCAGGGCGCCACTGCCCGCCCCCCAAGGAGCCGATGTGGAGCTGGCCGATAAGgcc                                                                           >  8:295165/1‑136 (MQ=255)
                                               ccAGATGTGCATCGACGGCCAATGGCGCGATGCCAGAGCGGCAAGACCTTCGACAGCCTCAACCCGGCCACCGCCCCCGCCGGGGCGCACCTGCCCCACGCCGATGCAGCCGCTGTGGAGCTGGCCGTTCAGgcc                                                                           >  4:425510/1‑135 (MQ=255)
                                                   aggcgCATCCCCGGCCACTGGCGCGCGGCCCAGAGCGGCAAGCCCTACGACCCCCTCACCCCGGCCCCCGCCCCGGCCGGGCCGCAACTGCCCGACGCCGATGAAGCCGATGTGGAGCTGGCCGTTCAGGCCGccc                                                                       <  7:107368/132‑1 (MQ=255)
                                                          tCGACGGCCAATGGCGCGATGCCCAGAGCGGCAAGACCTTCGACAGCCTCACCCCGGCCACCGCCCCCGCCACGGCGCCCCTGCCCCCCGCCGCTGAAGCCGATGTGGGGCGGGCGGTTCTGGCCGCCCCGCGAGc                                                                >  8:21796/1‑136 (MQ=255)
                                                                   aaTGGCGCGATGCCCAGAGCGGCAAGACCTTCGACAGCCTCAACCCGGCCACCGCCCCCGCCGGGGGGCCCCTGCCCGCCGCCGCTGAAGCCGCTGTGGAGCTGGGCGGTGAGGCCGCCCCGCGGGCCGTCGCCa                                                        >  1:247214/1‑135 (MQ=255)
                                                                       ccgcgAGGCCCCGCGCGGCAAGACCTTCGCCAGCCTCACCCCGCCCCCCGCCCAGGCCGGGGCGCAACTGCCCGACGCCGATGAAGCCGATGTGGAGCTGGCCGTTCAGGCCGCCCAGCGTGCCTTCGACAGCAAg                                                   <  8:233975/135‑1 (MQ=255)
                                                                        cgcgCTGCCCACCGCGGCACGACCAGCGACCGCCTCACCCCGGCCACCGCCCAGGCCGGGGCGCAACTGCCCGCCGCCGATGAAGCCGAGGTGGAGCTGGCCGTTCAGGCCGCCCAGCGTGCCTTCGACAGCAAg                                                   <  6:275076/135‑1 (MQ=255)
                                                                           gATGCCCCGAGCGGCAAGCCCTTCGACCGCCGCCACCCGCCCCCCGCCCAGGCCGGGGCGCAACTGCCCGACGCCGATGAAGCCGATGTGGAGCTGGCCGTTCAGGCCGCCCAGCGTGCCTTCGACAGCAAGGCg                                                <  5:354613/135‑1 (MQ=255)
                                                                             tGCCCAGAGCGGCAAGACCTTCGACAGCCTCAACCCGGCCACCGCCCCGGCCTGGGCGCCACTGCCCGCCGCCGATGAAGCCGCTGTGGCGCTGGCCGTTCAGGCCGCCCCGCGGGCCGTCGCCAGCAAGGCGTg                                              >  5:223691/1‑135 (MQ=255)
                                                                             tGCCCAGAGCGGCAAGACCTTCGACAGCCTCAACCCGGCCACCGCCCCGGCCTGGGCGCAACTGCCCGCCGCCGATGACTCCGATGTGGAGCTGGCCGTTCAGGCCGCCCAGCGTGCCGTGGCCAGCAAGGCGTgg                                             >  7:322178/1‑136 (MQ=255)
                                                                                 cAGAGCGCCAAGACCTTCGACAGCCTCACCCCCGCCACCCCCCAGGCCGGCGCGCAACTGCCCGACGCCGATGAAGCCGATGTGGAGCTGGCCGTTCAGGCCGCCCAGCGTGCCTTCGACAGCAAGGCGTGgcgc                                          <  1:157485/135‑1 (MQ=255)
                                                                                   gcgCGGCAAGACCTTCGACAGCCTCAACCCGGCCACCGCCCCGGCCTGGGCGCAACTGCCCGACGCCGATGAAGCCGATGTGGAGCTGGCCGTTCAGGCCGCCCAGCGTGCCTTCGACAGCAAGGCGTGGCgcag                                        <  2:247214/133‑1 (MQ=255)
                                                                                          aaGCCCTTCGACAGCCTCAACCCGGCCACACCCCCGGCCGGGGCCCAACAGCCCGCCGCCGAGGACGGCGCTGTGGAGCGGGCCGTTCAGGCCCCCCAGCGGTCCTGCGCCCGCCAGGCGTGGCGCGCCGTCAcc                                 >  6:112932/1‑135 (MQ=255)
                                                                                                    aCAGCCTCAACCCGGCCACCGCCCAGGCCTGGGCGCAACTGCCCGACGCCGATGAAGCCGATGTGGAGCTGGCCGTTCAGGCCGCCCAGCGTGCCTTCGACAGCAAGGCGTGGCGCAGCATCACCGCTAccgccc                       <  4:278921/135‑1 (MQ=255)
                                                                                                      cgCCTCACCCCGCCCCCCGCCCAGCCCGGGGCGCAACTGCCCGACGCCGATGCAGCCGAGGTGGAGCTGGCCGGTCAGGCCGCCCAGCGTGCCTTCGACAGCAAGGCGTGGCGCAGCATCACCGCTAccgcccgc                     <  5:131529/134‑1 (MQ=255)
                                                                                                       gCCTCACCCCGGCCACCGCCCAGGCCTGGGCGCAACTGCCCGACGCCGATGAAGCCGATGTGGAGCTGGCCGTTCAGGCCGCCCAGCGTGCCTTCGACAGCAAGGCGTGGCGCAGCATCACCGCTACCGCCcgcg                    <  2:250779/135‑1 (MQ=255)
                                                                                                           cAACCCGGCCACCGCCCCGGCCTGGGCGCAACTGCCCGCCGCCGATGAAGCCGATGTGGGGCTGGCCGTTCAGGCCGCCCAGCGTGGCGTCGACAGCAAGGCGGGGCGCAGCATAACCCCTACCCCCCGCGGCCAg               >  8:230744/1‑136 (MQ=255)
                                                                                                               gcGGCCCCCGCCCAGGCGGGGCCGCAACTGCCCGACGCCGATGAAGCCGATGTGGAGCTGGCCGGTCAGGCCGCCCAGCGTGCCTTCGACAGCAAGGCGTGGCGCAGCATCACCGCTACCGCCCGCGGCAAgctg            <  5:304998/134‑1 (MQ=255)
                                                                                                                 tgCCACCGCCCCGGCCTGGGCGCAACTGCCCGACGCCGATGAAGCCGATGTGGAGCTGGCCGTTCAGGCCGCCCAGCGTGCGTTCGACAGCAAGGCGTGGCGCAGCATCACCGCTACCGCCCGCGGCAAgctgct          <  2:170712/134‑1 (MQ=255)
                                                                                                                 ggCCACCGCCCAGGCCTGGGCGCAACTGCCCGACGCCGATGAAGCCGATGTGGAGCTGGCCGTTCAGGCCGCCCAGCGTGCCTTCGACAGCAAGGCGTGGCGCAGCATCACCGCTACCGCCCGCGGCAAgctgct          <  6:340284/135‑1 (MQ=255)
                                                                                                                    ccccGCCCCGGCGGGGGCGCAACTGCCCGACGCCGATGAAGCCGATGTGGAGCTGGCCGTTCAGGCCGCCCAGCGTGCCTTCGACAGCAAGGCGTGGCGCAGCATCACCGCTACCGCCCGCGGCAAGCTGCTGcgt      <  8:80024/134‑1 (MQ=255)
                                                                                                                     aCCGCCCCGGCCTGGGCGCCACTGCCCGCCGCCGATGAAGCCGATGTGGAGCTGGCCGTTCAGGCCGCCCCGCGTGCCGTCGCCAGCAAGGCGTGGCGCCGCATCACCGCGACCGCCCGCGGCCAGCGGCTGcgt      >  2:59710/1‑135 (MQ=255)
                                                                                                                     aCCGCCCCGGCCTGGGCGCAACTGCCCCCCGCCGCGGAAGCCGATGTGGAGCTGGCCGTTCAGGCCGCCCAGCGTGGCGGGGACAGCAAGGCGTGGCGCCGCCTCACCGCTACCGCCCCCGGCCAgcggcggcgg      >  2:15706/1‑134 (MQ=255)
                                                                                                                         cccTGGCCTGGGCGCAACTGCCCGACGCCGATGAATCCGTTGTGGAGCTGGCCGTTCAGGCCGCCCAGCGTGCCTTCGACAGCAAGGCGTGGCGCAGCATCACCTCTACCGCCCGCGGCAAGCTGCTGCgtcgtc  <  4:382722/135‑1 (MQ=255)
                                                                                                                                 |                                                                                                                              
CCTCACCGAATCCGCCAAGGGGATTGTTGCATGACCCTCGTTCGTTTCCAGATGTGCATCGACGGCCAATGGCGCGATGCCCAGAGCGGCAAGACCTTCGACAGCCTCAACCCGGCCACCGCCCAGGCCTGGGCGCAACTGCCCGACGCCGATGAAGCCGATGTGGAGCTGGCCGTTCAGGCCGCCCAGCGTGCCTTCGACAGCAAGGCGTGGCGCAGCATCACCGCTACCGCCCGCGGCAAGCTGCTGCGTCGTC  >  KT2440/4135483‑4135738

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