Predicted mutation
evidence seq id position mutation freq annotation gene description
RA KT2440 86,554 G→A 6.9% Y188Y (TAC→TAT limB ← Limonene 1,2‑monooxygenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*KT244086,5540GA 6.9% 179.3 / 3.6 58Y188Y (TAC→TATlimBLimonene 1,2‑monooxygenase
Reads supporting (aligned to +/- strand):  ref base G (32/22);  new base A (2/2);  total (34/24)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.63e-01

CTGCGCCACCACCATGGGAATGCCCAGCATCACATACGGCTTGTCCAGCGTGGTCGAGGGCTTGAAATGGTTGCGGTACACGCGAATGGCCTCGTGCATGTAGCGTGGCGCGAAATGCGAGGCGAAGGCGTAGGGCATGCCGCGCATGCCCGCCAGCTGGGCACTGAACAGGCTGGATCCCAGCAACCACATCGGCACGTCGGTGTCGTACCCCGGTACGGCGATCACTTTTTGATCATCAGTGCGTGGCCCGAGGTAGCGC  >  KT2440/86425‑86686
                                                                                                                                 |                                                                                                                                    
cTGCGCCACCACCATGGGAATGCCCAGCATCACATACGGCTTGTCCAGCGTGGTCGAGGGCTTGAAATGGTTGCGGTACACGCGAATGGCCTCGTGCATGTAGCGTGGCGCGAAATGCGAGGCGAAGGCGTAgg                                                                                                                                  >  3:109880/1‑134 (MQ=255)
    gcCACCACCATGGGAATGCCCAGCATCACATACGGCTTGTCCAGCGTGGTCGAGGGCTTGAAATGGTTGCGGTACACGCGAATGGCCTCGTGCATGTAGCGTGGCGCGAAATGCGAGGCGAAGGCGTAGGGc                                                                                                                                <  5:219115/132‑1 (MQ=255)
                ggAATGCCCAGCATCACATACGGCTTGTCCAGCGTGGTCGAGGGCTTGAAATGGTTGCGGTACACGCGAATGGCCTCGTGCATGTAGCGTGGCGCGAAATGCGAGGCGAAGGCGTAGGGCATGCCGCGCATGccc                                                                                                                 >  6:141339/1‑135 (MQ=255)
                 gAATGCCCAGCATCACATAGGGCTTGTCCAGCGTGGTCGAGGGCTTGAAATGGTTGCGGTACACGCGAATGGCCTCGTGCATGTAGCGCGGCGCGAAATGCTAGGATAAGGCGTAGGGCTTGCCGCGCATGcc                                                                                                                  <  8:220484/133‑1 (MQ=255)
                    tGCCCCGCATCACATACGGCTGGTCCAGCGTGGTCGCGGGCTTGAAATGGTTGCGGTACCCGCGAATGGCCTCGTGCATGTAGCGTGGCGCGAAATGCGCGGCGAAGGCGTAGGGCATGCCGCGCATGCCCGCCa                                                                                                             <  8:97979/135‑1 (MQ=255)
                    tGCCCAGCATCACATACGGCTTGGCCAGCGTGGTCGAGGGCTTGAAATGGTTGCGGTACACGCGAATGGCCTCGTGCATGTAGCGTGGCGCGAAATGCGAGGCGAAGGCGTAGGGCATGCCGCGCATGCCCGCCa                                                                                                             <  3:142424/135‑1 (MQ=255)
                      cccAGCATCACATACGGCTTGTCCAGCGTGGTCGAGGGCTTGAAATGGTTGCGGTACACGCGAATGGCCTCGTGCATGTAGCGTGGCGCGAAATGCGAGGCGAAGGCGTAGGGCATGCCGCGCATGCCCGCCAGc                                                                                                           >  4:157564/1‑135 (MQ=255)
                         aGCATCACATACGGCTTGTCCAGCGTGGTCGAGGGCTTGAAATGGTTGCGGTACACGCGAATGGCCTCGTGCATGTAGCGTGGCGCGAAATGCGAGGCGAAGGCGTAGGGCATGCCGCGCATGCCCGCCAGCTgg                                                                                                        >  1:232231/1‑135 (MQ=255)
                         aGCATCACATACGGCTTGTCCAGCGTGGTCGAGGGCTTGAAATGGTTGCGGTACACGCGAATGGCCTCGTGCATGTAGCGTGGCGCGAAATGCGAGGCGAAGGCGTAGGGCATGCCGCGCATGCCCGCCAGCTgg                                                                                                        >  1:283909/1‑135 (MQ=255)
                         aGCATCACATACGGCTTGTCCAGCGTGGTCGAGGGCTTGAAATGGTTGCGGTACACGCGAATGGCCTCGTGCATGTAGCGTGGCGCGAAATGCGAGGCGAAGGCGTAGGGCATGCCGCGCATGCCCGCCAGCTgg                                                                                                        >  1:294126/1‑135 (MQ=255)
                         aGCATCACATACGGCTTGTCCAGCGTGGTCGAGGGCTTGAAATGGTTGCGGTACACGCGAATGGCCTCGTGCATGTAGCGTGGCGCGAAATGCGAGGCGAAGGCGTAGGGCATGCCGCGCATGCCCGCCAGCTgg                                                                                                        <  5:153789/135‑1 (MQ=255)
                         aGCATCACATACGGCTTGTCCAGCGTGGTCGAGGGCTTGAAATGGTTGCGGTACACGCGAATGGCCTCGTGCATGTAGCGTGGCGCGAAATGCGAGGCGAAGGCGTAGGGCATGCCGCGCATGCCCGCCAGCTgg                                                                                                        >  6:301837/1‑135 (MQ=255)
                            aTCACATACGGCTTGTCCAGCGTGGTCGAGGGCTTGAAATGGTTGCGGTACACGCGAATGGCCTCGTGCATGTAGCGTGGCGCGAAATGCGAGGCGAAGGCGTAGGGCATGCCGCGCATGCCCGCCAGCTGGGCa                                                                                                     >  3:87532/1‑135 (MQ=255)
                              cacaTACGGCTTGTCCCGCGTGGTCGCGGGCTTGAAATGGTTGCGGTACCCGCGAATGGCCTCGTGCATGTAGCGTGGCGCGAAATGCGAGGCGAAGGCGTAGGGCATGCCGCGCATGCCCGCCAGCTGGGCACt                                                                                                   <  6:301663/135‑1 (MQ=255)
                              cacaTACGGCTTGTCCAGCGTGGTCGAGGGCTTGAAATGGTTGCGGTACACGCGATTGGCCTCGTGCATGTAGGGTGGCGCGAAATGCGAGGCGAAGGCGTAGGGAATGCCGCGCATGCCCGCCAGCTGGGCACt                                                                                                   <  3:142756/135‑1 (MQ=255)
                              caaaTACGGCTGGTCCAGCGTGGTCGAGGGCTTTAAATGGTTGCGGAACACGCGAATGGCCTCTGGCATGTAGCGTGGCGCTAAATGCGCGGCTAAGGCGTAGGGCAGGCCGCGCATGCCCGCCAGCTGGGCACt                                                                                                   <  6:322994/132‑1 (MQ=255)
                                        ttGTCCAGCGTGGTCGAGTGCTTGAAATGGTTGCGGTACACGCGAATGGCCTCGTGCATGTAGCGTGGCGCGAAATGCGAGGCGAAGGCGTAGTGCATGCCGCGCATGCCCGACAGCTGGGCACTGAACAGGCt                                                                                          >  5:16210/1‑134 (MQ=255)
                                        ttGTCCAGCGTGGTCGAGGGCTTGAAATGGTTGCGGTACACGCGAATGGCCTCGTGCATGTAGCGTGGCGCGAAATGCGAGGCGAAGGCGTAGGGCATGCCGCGCATGCCCGCCAGCTGGGCACTGAACAGGCt                                                                                          >  3:206973/1‑134 (MQ=255)
                                          gTCCAGCGTGGTCGAGGGCTTGAAATGGTTGCGGTACACGCGAATGGCCTCGTGCATGTAGCGTGGCGCGAAATGCGAGGCGAAGGCGTAGGGCATGCCGCGCATGCCCGCCAGCTGGGCACTGAACAGGCTGGa                                                                                       >  2:238937/1‑135 (MQ=255)
                                          gTCCAGCGTGGTCGAGGGCTTGAAATGGTTGCGGTACACGCGAATGGCCTCGTGCATGTAGCGTGGCGCGAAATGCGAGGCGAAGGCGTAGGGCATGCCGCGCATGCCCGCCAGCTGGGCACTGAACAGGCTGGa                                                                                       >  8:78011/1‑135 (MQ=255)
                                             cAGCGTGGTCGAGGGCTTGAAATGGTTGCGGTACACGCGAATGGCCTCGTGCATGTAGCGTGGCGCGAAATGCGAGGCGAAGGCGTAGGGCATGCCGCGCATGCCCGCCAGCTGGGCACTGAACAGGCTGGATcc                                                                                    >  5:324593/1‑135 (MQ=255)
                                             cAGCGTGGTCGAGGGCTTGAAATGGTTGCGGTACACGCGAATGGCCTCGTGCATGTAGCGTGGCGCGAAATGCGAGGCGAAGGCGTAGGGCATGCCGCGCATGCCCGCCAGCTGGGCACTGAACAGGCTGGATc                                                                                     >  1:42554/1‑134 (MQ=255)
                                              aGCGTGGTCGAGGGCTTGAAATGGTTGCGGTACACGCGAATGGCCTCGTGCATGTAGCGTGGCGCGAAATGCGAGGCGAAGGCGTAGGGCATGCCGCGCATGCCCGCCAGCTGGGCACTGAACAGGCTGGATc                                                                                     >  1:161312/1‑133 (MQ=255)
                                                        aGGGCTTGAAATGGTTGCGGTACACGCGAATGGCCTCGTGCATGTAGCGTGGCGCGAAATGCGAGGCGAAGGCGTAGGGCATGCCGCGCATGCCCGCCAGCTGGGCACTGAACAGGCTGGa                                                                                       <  4:39977/121‑1 (MQ=255)
                                                        aGGGCTTGAAATGGTTGCGGTACACGCGAATGGCCTCGTGCATGTAGCGTGGCGCGAAATGCGAGGCGAAGGCGTAGGGCATGCCGCGCATGCCCGCCAGCTGGGCACTGAACAGGCTGGa                                                                                       >  3:39977/1‑121 (MQ=255)
                                                        aGGGCTTGAAATGGTTGCGGTACACGCGAATGGCCTCGTGCATGTAGCGTGGCGCGAAATGCGAGGCGAAGGCGTAGGGCATGCCGCGCATGCCCGCCAGCTGGGCACTGAACAGGCTGGATCCCAGCAACcaca                                                                         >  4:288880/1‑135 (MQ=255)
                                                        aGGGCTTGAAATGGTTGCGGTACACGCGAATGGCCTCGTGCATGTAGCGTGGCGCGAAATGCGAGGCGAAGGCGTAGGGCATGCCGCGCATGCCCGCCAGCTGGGCACTGAACAGGCTGGATCCCAGCAACcaca                                                                         >  1:187317/1‑135 (MQ=255)
                                                          ggCTTGAAATGGTTGCGGTACACGCGAATGGCCTCGTGCATGTAGCGTGGCGCGAAATGCGAGGCGAAGGCGTAGGGCATGCCGCGCATGCCCGCCAGCTGGGCACTGAACAGGCTGGATCCCAGCAACCACATc                                                                       >  4:159108/1‑135 (MQ=255)
                                                                 aaTGGTTGCGGTACACGCGAATGGCCTCGTGCATGTAGCGTGGCGCGAAATGCGAGGCGAAGGCGTAGGGCATGCCGCGCATGCCCGCCAGCTGGGCACTGAACAGGCTGGATCCCAGCAACCACATCGGCACg                                                                 <  6:324593/134‑1 (MQ=255)
                                                                    ggTTGCGGTACACGCGAATGGCCTCGTGCATGTAGCGTGGCGCGAAATGCGAGGCGAAGGCGTAGGGCATGCCGCGCATGCCCGCCAGCTGGGCACTGAACAGGCTGGATCCCAGCAACCACATCGGCACGTCgg                                                             <  2:42554/135‑1 (MQ=255)
                                                                        gCGGTACACGCGAATGGCCTCGTGCATGTAGCGTGGCGCGAAATGCGAGGCGAAGGCGTAGGGCATGCCGCGCATGCCCGCCAGCTGGGCACTGAACAGGCTGGATCCCAGCAACCACATCGGCACGTCGGTGTc                                                         >  4:15598/1‑135 (MQ=255)
                                                                        gCGGTACACGCGAATGGCCTCGTGCATGTAGCGTGGCGCGAAATGCGAGGCGAAGGCGTAGGGCATGCCGCGCATGCCCGCCAGCTGGGCACTGAACAGGCTGGATCCCAGCAACCACATCGGCACGTCGGTGTc                                                         >  3:80705/1‑135 (MQ=255)
                                                                        gCGGTACACGCGAATGGCCTCGTGCATGTAGCGTGGCGCGAAATGCGAGGCGAAGGCATAGGGCATGCCGCGCATGcc                                                                                                                  <  4:325836/78‑1 (MQ=255)
                                                                        gCGGTACACGCGAATGGCCTCGTGCATGTAGCGTGGCGCGAAATGCGAGGCGAAGGCATAGGGCATGCCGCGCATGcc                                                                                                                  <  4:330973/78‑1 (MQ=255)
                                                                        gCGGTACACGCGAATGGCCTCGTGCATGTAGCGTGGCGCGAAATGCGAGGCGAAGGCATAGGGCATGCCGCGCATGcc                                                                                                                  >  3:330973/1‑78 (MQ=255)
                                                                        gCGGTACACGCGAATGGCCTCGTGCATGTAGCGTGGCGCGAAATGCGAGGCGAAGGCATAGGGCATGCCGCGCATGcc                                                                                                                  >  3:325836/1‑78 (MQ=255)
                                                                         cGGTACACGCGAATGGCCTCGTGCATGTAGCGTGGCGCGAAATGCGAGGCGAAGGCGTAGGGCATGCCGCGCATGCCCGCCAGCTGGGCACTGAACAGGCTGGATCCCAGCAACCACATCGGCACGTCGGTGTc                                                         >  8:57239/1‑134 (MQ=255)
                                                                         cGGTACACGCGAATGGCCTCGTGCATGTAGCGTGGCGCGAAATGCGAGGCGAAGGCGTAGGGCAGGCCGCGCATGCCCGCCAGCTGGGCACTGAACAGGCTGGATCCCAGCAACCACATCGGCACGTCTGTGTc                                                         <  5:301837/134‑1 (MQ=255)
                                                                                       ggCCTCGTGCATGTAGCGTGGCGCGAAATGCGAGGCGAAGGCGTAGGGCATGCCGCGCATGCCCGCCAGCTGGGCACTGAACAGGCTGGATCCCAGCAACCACATCGGCACGTCGGTGTCGTACCCCGGTACgg                                           >  6:298336/1‑134 (MQ=255)
                                                                                        gCCTCGGGCATGTCGCGGGGCGCTAAATGCGTGGCGAAGGCGTAGGGCATGCCGCGCAGCCCCGCCAGCTGGGCACTGAACAGGCTGGATCCCAGCAACCACATCGGCACGTCGGTGTCGTACCCCGGTACGGc                                          <  4:223337/134‑1 (MQ=255)
                                                                                         ccTCGTGCATGTAGCGTGGCGCGAAATGCGAGGCGAAGGCGTAGGGCATGCCGCGCATGCCCGCCAGCTGGGCACTGAACAGGCTGGATCCCAGCAACCACATCGGCACGTCGGTGTCGTACCCCGGTACGGCGa                                        >  3:333423/1‑135 (MQ=255)
                                                                                         ccTCGTGCATGTAGCGTGGCGCGAAATGCGAGGCGAAGGCGTAGGGCATGCCGCGCATGCCCGCCAGCTGGGCACTGAACAGGCTGGATCCCAGCAACCACATCGGCACGTCGGTGTCGTACCCCGGTACGGCGa                                        >  8:100756/1‑135 (MQ=255)
                                                                                                cATGTAGCGTGGCGCGAAATGCGAGGCGAAGGCGTAGGGCATGCCGCGCATGCCCGCCAGCTGGGCACTGAACAGGCTGGATCCCAGCAACCACATCGGCACGTCGGTGTCGTACCCCGGTACGGCGATCACttt                                 >  4:203458/1‑135 (MQ=255)
                                                                                                      gCGTGGCGCGAAATGCGAGGCGAAGGCGTAGGGCATGCCGCGCATGCCCGCCAGCTGGGCACTGAACAGGCTGGATCCCAGCAACCACATCGGCACGTCGGTGTCGTACCCCGGTACGGCGATCACTTTTTGatc                           >  1:268725/1‑135 (MQ=255)
                                                                                                          ggCGCGAAATGCGAGGCGAAGGCGTAGGGCATGCCGCGCATGCCCGCCAGCTGGGCACTGAACAGGCTGGATCCCAGCAACCACATCGGCACg                                                                 <  1:70451/93‑1 (MQ=255)
                                                                                                          ggCGCGAAATGCGAGGCGAAGGCGTAGGGCATGCCGCGCATGCCCGCCAGCTGGGCACTGAACAGGCTGGATCCCAGCAACCACATCGGCACg                                                                 >  2:70451/1‑93 (MQ=255)
                                                                                                          ggCGCGAAATGCGAGGCGAAGGCGTAGGGCATGCCGCGCATGCCCGCCAGCTGGGCACTGAACAGGCTGGATCCCAGCAACCACATCGGCACGTCGGTGTCGTACCCCGGTACTGCGATCACTTTTTGAtcatca                       >  8:159848/1‑135 (MQ=255)
                                                                                                          ggCGCGAAATGCGAGGCGAAGGCGTAGGGCATGCCGCGCATGCCCGCCAGCTGGGCACTGAACAGGCTGGATCCCAGCAACCACATCGGCACGTCGGTGTCGTACCCCGGTACGGCGATCACTTTTTGAtcatca                       >  7:212489/1‑135 (MQ=255)
                                                                                                           gcggGATATGCGCGGGGAAGGCGTAGGGCATGCCGCGCTTGCCCGCGAGCTGGGCACTGAACAGGCTGGATCCCAGCAACCACATCGGCACGTCGGTGTCGTACCCGGGTACGGCGATCTCTTTTTGATCATCAg                      <  3:83729/135‑1 (MQ=255)
                                                                                                            cgctAATTGCGAGGCGAAGGCGTAGGGCATGCCGGGCTTGCCCGCCAGCTGGGCACTGAACAGGCTGGATCCCAGCAACCACATCGGCACGTCGGTGTCGTACCCCGGTACGGCGATCACTTTTTGATCATCAGt                     <  4:206973/135‑1 (MQ=255)
                                                                                                              cgAAATGCGAGGCGAAGGCGTAGGGCATGCCGCGCATGCCCGCCAGCTGGGCACTGAACAGGCTGGATCCCAGCAACCACATCGGCACGTCGGTGTCGTACCCCGGTACGGCGATCACTTTTTGATCATCAGTGc                   >  7:291356/1‑135 (MQ=255)
                                                                                                               gAAATGCGAGGCGAAGGCGTAGGGCATGCCGCGCATGCCCGCCAGCTGGGCACTGAACAGGCTGGATCCCAGc                                                                                <  2:68485/73‑1 (MQ=255)
                                                                                                               gAAATGCGAGGCGAAGGCGTAGGGCATGCCGCGCATGCCCGCCAGCTGGGCACTGAACAGGCTGGATCCCAGc                                                                                >  1:68485/1‑73 (MQ=255)
                                                                                                                       cggggAAGGCGTAGGGCATGCCGCGCATGCCCGCCAGCTGGGCACTGAACAGGCTGGATCCCAGCAACCACATCGGCACGTCGGTGTCGTACCCCGGTACGGCGATCACTTTTTGATCATCAGTGCGTGGCCCGa          <  8:229924/131‑1 (MQ=255)
                                                                                                                       aGGCGAAGGCGTAGGGCATGCCGCGCATGCCCGCCAGCTGGGCACTGAACAGGCTGGATCCCAGCAACCACATCGGCACGTCGGTGTCGTACCCCGGTACGGCGATCACTTTTTGATCATCAGTGCGTGGCCCGa          <  3:15598/135‑1 (MQ=255)
                                                                                                                       aGGCGAAGGCGTAGGGCATGCCGCGCATGCCCGCCAGCTGGGCACTGAACAGGCTGGATCCCAGCAACCACATCGGCACGTCGGTGTCGTACCCCGGTACGGCGATCACTTTTTGATCATCAGTGCGTGGCCCGa          <  3:157564/135‑1 (MQ=255)
                                                                                                                       aGGCGAAGGCGTAGGGCATGCCGCGCATGCCCGCCAGCTGGGCACTGAACAGGCTGGATCCCAGCAACCACATCGGCACGTCGGTGTCGTACCCCGGTACGGCGATCACTTTTTGATCATCAGTGCGTGGCCCGa          <  4:87532/135‑1 (MQ=255)
                                                                                                                               gCGTAGGGCATGCCGCGCATGCCCGCCAGCTGGGCACTGAACAGGCTGGATCCCAGCAACCACATCGGCACGTCGGTGTCGTACCCCGGTACGGCGATCACTTTTTGATCATCAGTGCGTGGCCCGAGGTAgcgc  <  3:159108/135‑1 (MQ=255)
                                                                                                                                 |                                                                                                                                    
CTGCGCCACCACCATGGGAATGCCCAGCATCACATACGGCTTGTCCAGCGTGGTCGAGGGCTTGAAATGGTTGCGGTACACGCGAATGGCCTCGTGCATGTAGCGTGGCGCGAAATGCGAGGCGAAGGCGTAGGGCATGCCGCGCATGCCCGCCAGCTGGGCACTGAACAGGCTGGATCCCAGCAACCACATCGGCACGTCGGTGTCGTACCCCGGTACGGCGATCACTTTTTGATCATCAGTGCGTGGCCCGAGGTAGCGC  >  KT2440/86425‑86686

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