Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,796,177 G→T 9.3% P107Q (CCG→CAG)  fcs ← feruloyl‑CoA‑synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,796,1770GT 9.3% 151.1 / 3.4 51P107Q (CCG→CAG) fcsferuloyl‑CoA‑synthetase
Reads supporting (aligned to +/- strand):  ref base G (24/22);  new base T (3/2);  total (27/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 = 9.96e-01

TCGCAGACATGGCGCAACTTGGCGAAGTCTTGCGACAACAGCGCGTAGGCCGGCGACACCGGGCAATAGGCAATACCGGCATACATGGCGCCGAGGGCGATTTGCAGGTGTTCGATGTCGTTGCCGGAAAGCAGCGCCAGCGGGCGCTCGGCACTGAGGCCCAGTCCTAGCAAGTTGGCGGCGATGGTGCGCACATCGGCGAGCATCTGCACGTAGCTGATCGAACGCCAGGCACCGTCTGCCTGGCGTGCCGCGATGAAAGTGGTGTCCGGGCGCACCTGGGCCCAATGC  >  NC_002947/3796036‑3796326
                                                                                                                                             |                                                                                                                                                     
tCGCAGACATGGCGCAACTTGGCGAAGTCTTGCGACAACAGCGCGTAGGCCGGCGACACCGGGCAATAGGCAATACCGGCATACATGGCGCCGAGGGCGATTTGCAGGTGTTCGATGTCGTTGCCTGAAAGCAGCGCCAGCGGGCGCTCg                                                                                                                                               >  1:391214/1‑150 (MQ=255)
      aCATGGCCAAAATTGGAAAAGTTTTGCGATAACATCGCGTAGGCCGGCGACACCGGGCAATAGGCAATACCGGCATACATGGCGCCGAGGGCGATTTGCAGGTGTTCGATGTCGTTGCCGGAAAGCAGCGCCAGCGGGCGCTCGGCACTg                                                                                                                                         >  2:728783/1‑150 (MQ=255)
        aTGGCGCAACTTGGCGAGGTCTTGCGACAACAGCGCGTAGGCCGGCGACACCGGGCAATAGGCAATACCGGCATACATGGCGCCGAGGGCGATTTGCAGGTGTTCGATGTCGTTGCCGGAAAGCAGCGCCAGCGGGCGCTCGGCACTGAg                                                                                                                                       >  1:956048/1‑150 (MQ=255)
         tGGCGCAACCTGGCGAAGTCTTGCGACAACAGCGCGTAGGCCGGCGACACCGGGCAATCGGTATTACCGGCATACCTAGCGCCGAGGGCGATTTGCAGGTGTTCGATGTCGTTGCCGGAAAGCAGCGCCAGCGGGCGCTCGGCACTGAgg                                                                                                                                      <  1:460953/150‑1 (MQ=255)
           gcgcAACTTGGCGAAGTCTTGCGACAACAGCGCGTAGGCCGGCGACACCGGGCAATAGGCAATACCGGCATACATGGCGCCGAGGGCGATTTGCAGGTGTTCGATGTCGTTGCCGGAAAGCAGCGCCAGCGGGCGCTCGGCACTGAGGcc                                                                                                                                    >  2:1297798/1‑150 (MQ=255)
             gcAACTTGGCGAAGTCTTGCGACAACAGCGCGTAGGCCGGCGACACCGGGCAATAGGCAATACCGGCATACATGGCGCCGAGGGCGATTTGCAGGTGTTCGATGTCGTTGCCGGAAAGCAGCGCCAGCGGGCGCTCGGCACTGAGGCCCa                                                                                                                                  <  1:1020372/150‑1 (MQ=255)
                   tGGCGAAGTCTTGCGACAACAGCGCGTAGGCCGGCGACACCGGGCAATAGGCAATACCGGCATACATGGCGCCGAGGGCGATTTGCAGGTGTTCGATGTCGTTGCCGGAAAGCAGCGCCAGCGGGCGCTCGGCACTGAGGCCCAGTCCTa                                                                                                                            <  1:655542/150‑1 (MQ=255)
                       gAAGTCTTGCGACAACAGCGCGTAGGCCGGCGACACCGGGCAATAGGCAATACCGGCATACATGGCGCCGAGGGCGATTTGCAGGTGTTCGATGTCGTTGCCGGAAAGCAGCGCCAGCGGGCGCTCGGCACTGAGGCCCAGTCCTAGCaa                                                                                                                        >  2:796838/1‑150 (MQ=255)
                          gTCTTGCGACAACAGCGCGTAGGCCGGCGACACCGGGCAATAGGCAATACCGGCATACATGGCGCCGAGGGCGATTTGCAGGTGTTCGATGTCGTTGCCGGAAAGCAGCGCCAGCGGGCGCTCGGCACTGAGGCCCAGTCCTAGCAAGtt                                                                                                                     >  2:720904/1‑150 (MQ=255)
                                 gACAACAGCGCGTAGGCCGGCGACAGCGGGCAATAGGCAATACCGGCATACATGGCGCCGAGGGCGATTTGCAGGTGTTCGATGTCGTTGCCGGAAAGCAGCGCCAGCGGACGCTCGGCACTGAGGCCCAGTCCTAGCAAGTTGgcggcg                                                                                                              >  1:1087283/1‑150 (MQ=255)
                                  acaacaGCGCGTAGGCCGGCGACACCGGGCAATAGGCAATACCGGCATACATGGCGCCGAGGGCGATTTGCAGGTGTTCGATGTCGTTGCCGGAAAGCAGCGCCAGCGGGCGCTCGGCACTGAGGCCCAGTCCTAGCAAGTTGGCGGCGa                                                                                                             <  1:686622/150‑1 (MQ=255)
                                     acaGCGCGTAGGCCGGCGACACCGGGCAATAGTCAATACCGGCATACATGGCGCCGAGGGCGATTTGCAGGTATTCGATGTCGTTGCCGGAAAGCAGCGCCAGCTGGCGCTCGGCACCGAGGCCCAGTCCTAGCAAGTTGGCGGCGATgg                                                                                                          >  2:290081/1‑150 (MQ=255)
                                      caGCGCGTAGGCCGGCGACACCGGGCAATAGGCAATACCGGCATACATGGCGCCGAGGGCGATTTGCAGGTGTTCGATGTCGTTGCCGGAAAGCAGCGCCAGCGGGCGCTCGGCACTGAGGCCCAGTCCTAGCAAGTTGGCGGCGATGGt                                                                                                         <  1:1297798/150‑1 (MQ=255)
                                       aGCGCGTAGGCCGGCGACACCGGGCACTAGGCAATCCCGGCATACATGGCGCCGAGGGCGCTTTGCAGGTGTTCGATGTCGTTGCCGGAAAGCACCGCCCGCGGGCGCTCGGCCCTGAGGCCCAGTCCTAGCACGCTGGCGGCGgtggtg                                                                                                        >  1:616689/1‑150 (MQ=255)
                                                  cGGCGACACCGGGCAATAGGCAATACCGGCATACATGGCGCCGAGGGCGATTTGCAGTTGTTCGATGTCGTTGCCGGAAAGCAGCGCCAGCGGGCGCTCGGCACTGAGGCCCAGTCCTAGCAAGTTGGCGGCGATGGTGCGCACATCGGc                                                                                             >  2:843886/1‑150 (MQ=255)
                                                  cGGCGACACCGGGCAATAGGCAATACCGGCATACATGGCGCCGAGGGCGATTTGCAGGTGTTCGATGTCGTTGCCGGAAAGCAGCGCCAGCGGGCGCTCGGCACTGAGGCCCAGTCCTAGCAAGTTGGCGGCGATGGTGCGCACATCGGc                                                                                             >  1:1401949/1‑150 (MQ=255)
                                                  cGGCGACACCGGGCAATAGGCAATACCGGCATACATGGCGCCGAGGGCGATTTGCAGGTGTTCAATGTCGTTGCCGGAAAGCAGCGCCAGCGGGCGCTCGGCACTGAGGCCCAGTCCTAGCAAGTTGGCGGCGATGGTGCGCACATCGGc                                                                                             <  2:64075/150‑1 (MQ=255)
                                                          ccGGGCAATAGGCAATACCGGCATACATGGCGCCGAGGGCGATTTGCAGGTGTTCGATGTCGTTGCCGGAAAGCAGCGCCAGCGGGCGCTCGGCACTGAGGCCCAGTCCTAGCAAGTTGGCGGCGATGGTGCGCACATCGGCGAGCATCt                                                                                     >  1:509525/1‑150 (MQ=255)
                                                           cGGGCAATAGGCAATACCGGCATACATGGCGCCGAGGGCGATTTGCAGGTGTTCGATGTCGTTGCCGGAAAGCAGCGCCAGCGGGCGCTCGGCACTGAGGCCCAGTCCTAGCAAGTTGGCGGCGATGGTGCGCACATCGGCGAGCATCTg                                                                                    >  1:1137445/1‑150 (MQ=255)
                                                             ggcaataggcaataCCGGCATACATGGCGCCGAGGGCGATTTGCAGGTGTTCGATGTCGTTGCCGGAAAGCAGCGCCAGCGGGCGCTCGGCACTGAGGCCCAGTCCTAGCAAGTTGGCGGCGATGGTGCGCACATCGGCGAGCATCTGCa                                                                                  >  1:444630/1‑150 (MQ=255)
                                                                aataggcaataCCGGCATACATGGCGCCGAGGGCGATTTGCAGGTGTTCGATGTCGTTGCCGGAAAGCAGCGCCAGCGGGCGCTCGGCACTGAGGCCCAGTCCTAGCAAGTTGGCGGCGATGGTGCGCACATCGGCGAGCATCTGCACGt                                                                               <  1:728783/150‑1 (MQ=255)
                                                                aataggcaataCCGGCATACATGGCGCCGAGGGCGATTTGCAGGTGTTCGATGTCGTTGCCGGAAAGCAGCGCCAGCGGGCGCTCGGCACTGAGGCCCAGTCCTAGCAAGTTGGCGGCGATGGTGCGCACATCGGCGAGCATCTGCACGt                                                                               <  1:1172901/150‑1 (MQ=255)
                                                                  taggcaataCCGGCATACATGGCGCCGAGGGCGATTTGCAGGTGTTCGATGTCGTTGCCGGAAAGCAGCGCCAGCGGGCGCTCGGCACTGAGGCCCAGTCCTAGCAAGTTGGCGGCGATGGTGCGCACATCGGCGAGCATCTGCACGTAg                                                                             <  2:1401949/150‑1 (MQ=255)
                                                                    ggcaataCCGGCATACATGGCGCCGAGGGCGATTTGCAGGTGTTCGATGTCGTTGCCGGAAAGCAGCGCCAGCGGGCGCTCGGCACTGAGGCCCAGTCCTAGCAAGTTGGCGGCGATGGTGCGCACATCGGCGAGCATCTGCACGTAGCt                                                                           <  1:720904/150‑1 (MQ=255)
                                                                     gcaataCCGGCATACATGGCGCCGAGGGCGATTTGCAGGTGTTCGATGTCGTTGCCGGAAAGCAGCGCCAGCggg                                                                                                                                                     <  2:854831/75‑1 (MQ=255)
                                                                     gcaataCCGGCATACATGGCGCCGAGGGCGATTTGCAGGTGTTCGATGTCGTTGCCGGAAAGCAGCGCCAGCggg                                                                                                                                                     >  1:854831/1‑75 (MQ=255)
                                                                     gcaataCCGGCATACATGGCGCCGAGGGCGATTTGCAGGTGTTCGATGTCGTTGCCGGAAAGCAGCGCCAGCGGGCGCTCGGCACTGAGGCCCAGTCCTAGCAAGTTGGCGGCGATGGTGCGCACATCGGCGAGCATCTGCACGTAGCTg                                                                          >  1:1278837/1‑150 (MQ=255)
                                                                     gcaataCCGGCATACATGGCGCCGAGGGCGATTTGCAGGTGTTCGATGTCGTTGCCGGAAAGCAGCGCCAGCGGGCGCTCGGCACTGAGGCCCAGTCCTAGCAAGTTGGCGGCGATGGTGCGCACATCGGCGAGCATCTGCACGTAGCTg                                                                          >  1:160383/1‑150 (MQ=255)
                                                                            cGGCATACATGGCGCCGAGGGCGATTTGCAGGTGTTCGATGTCGTTGCCGGAAAGCAGCGCCAGCGGGCGCTCGGCACTGAGGCCCAGTCCTAGCAAGTTGGCGGCGATGGTGCGCACATCGGCGAGCATCTGCACGTAGCTGATCGAAc                                                                   >  2:345095/1‑150 (MQ=255)
                                                                            cGGCATACATGGCGCCGAGGGCGATTTGCAGGTGTTCGATGTCGTTGCCGGAAAGCAGCGCCAGCGGGCGCTCGGCACTGAGGCCCAGTCCTAGCAAGTTGGCGGCGATGGTGCGCACATCGGCGAGCATCTGCACGTAGCTGATCGAAc                                                                   <  2:525474/150‑1 (MQ=255)
                                                                              gCATACATGGCGCCGAGGGCGATTTGCAGGTGTTCGATGTCGTTGCCGGAAAGCAGCGCCAGCGGGCGCTCGGCACTGAGGCCCAGTCCTAGCAAGTTGGCGGCGATGGTGCGCACATCGGCGAGCATCTGCACGTAGCTGATCGAACGc                                                                 >  2:883456/1‑150 (MQ=255)
                                                                                          ccGAGGGCGATTTGCAGGTGTTCGATGTCGTTGCGCGAAAGCAGCGCCAGCGGGCGCTCGGCACTGAGGCCCAGTCCTAGCAAGTTGGCGGCGATGGTGCGCACATCGGCGAGCATCTGCACGTAGCTGATCGAACGCCAGGCACCGTCt                                                     <  2:616689/150‑1 (MQ=255)
                                                                                          ccGAGGGCGATTTGCAGGTGTTCGATGTCGTTGCCGGAAAGCAGCGCCAGCGGGCGCTCGGCACTGAGGCCCAGTCCTAGCAAGTTGGCGGCGATGGTGCGCACATCGGCGAGCATCTGCACGTAGCTGATCGAACGCCAGGCACCGTCt                                                     <  2:444630/150‑1 (MQ=255)
                                                                                             aGGGCGATTTGCAGGTGTTCGATGTCGTTGCCGGAAAGCAGCGCCAGCGGGCGCTCGGCACTGAGGCCCAGTCCTAGCaa                                                                                                                        <  1:345585/80‑1 (MQ=255)
                                                                                             aGGGCGATTTGCAGGTGTTCGATGTCGTTGCCGGAAAGCAGCGCCAGCGGGCGCTCGGCACTGAGGCCCAGTCCTAGCaa                                                                                                                        >  2:345585/1‑80 (MQ=255)
                                                                                               ggCGATTTGCAGGTGTTCGATGTCGTTGCCGGAAAGCAGCGCCAGCGGGCGCTCGGCACTGAGGCCCAGTCCTAGCAAGTTGGCGGCGATGGTGCGCACATCGGCGAGCATCTGCACGTAGCTGATCGAACGCCAGGCACCGTCTGCCTg                                                <  2:956048/150‑1 (MQ=255)
                                                                                                gCGATTTGCAGGTGTTCGATGTCGTTGCCGGAAAGCAGCGCCAGCGGGCGCTCGGCACTGAGGCCCAGTCCTAGCAAGTTGGCGGCGATGGTGCGCACATCGGCGAGCATCTGCACGTAGCTGATCGAACGCCAGGCACCGTCTGCCTgg                                               >  2:256893/1‑150 (MQ=255)
                                                                                                           gtgtTCGATGTCGTTGCCGGAAAGCAGCGCCAGCGGGCGCTCGGCACTGAGGCCCAGTCCTAGCAAGTTGGCGGCGATGGTgc                                                                                                       <  1:1392760/83‑1 (MQ=255)
                                                                                                           gtgtTCGATGTCGTTGCCGGAAAGCAGCGCCAGCGGGCGCTCGGCACTGAGGCCCAGTCCTAGCAAGTTGGCGGCGATGGTgc                                                                                                       >  2:1392760/1‑83 (MQ=255)
                                                                                                                 gATGTCGTTGCCGGAAAGCAGCGCCAGCGGGCGCTCGGCACTGAGGCCCAGTCCTAGCAAGTTGGCGGCGATGGTGCGCACATCGGCGAGCATCTGCACGTAGCTGATCGAACGCCAGGCACCGTCTGCCTGGCGTGCCGCGATGAAAgt                              >  2:41227/1‑150 (MQ=255)
                                                                                                                 gATGTCGTTGCCGGAAAGCAGCGCCAGCGGGCGCTCGGCACTGAGGCCCAGTCCTAGCAAGTTGGCGGCGATGGTGCGCACATCGGCGAGCATCTGCACGTAGCTGATCGAACGCCAGGCACCGTCTGCCTGGCGTGCCGCGATGAAAgt                              <  1:843886/150‑1 (MQ=255)
                                                                                                                  aTGTCGTTGCCGGAAAGCAGCGCCAGCGGGCGCTCGGCACTGAGGCCCAGTCCTAGCAAGTTGGCGGCGATGGTGCGCACATCGGCGAGCATCTGCACGTAGCTGATCGAACGCCAGGCACCGTCTGCCTGGCGTGCCGCGATGAAAgtg                             <  1:41227/150‑1 (MQ=255)
                                                                                                                   tGTCGTTGCCGGAAAGCAGCGCCAGCGGGCGCTCGGCACTGAGGCCCAGTCCTAGCAAGTTGGCGGCGATGGTGCGCACATCGGCGAGCATCTGCACGTAGCTGATCGAACGCCAGGCACCGTCTGCCTGGCGTGCCGCGATGAAAgtgg                            >  1:668437/1‑150 (MQ=255)
                                                                                                                    gTCGTTGCCGGAAAGCAGCGCCAGCGGGCGCTCGGCACTGAGGCCCAGTCCTAGCAAGTTGGCGGCGATGGTGCGCACATCGGCGAGCATCTGCACGTAGCTGATCGAACGCCAGGCACCGTCTGCCTGGCGTGCCGCGATGAAAgtggt                           <  2:1137445/150‑1 (MQ=255)
                                                                                                                           ccGGAAAGCAGCGCCAGCGGGCGCTCGGCACTGAGGCCCAGTCCTAGCACGTTGGCGGCGATGGTGCGCACATCGGCGAGCATCTGCACGTAGCTGATCGAACGCCAGGCACCGTCTGCCTGGCGTGCCGCGATGGAAGTGGTGTCCggg                    <  2:1087283/150‑1 (MQ=255)
                                                                                                                             ggAAAGCAGCGCCAGCGGGCGCTCGGCACTGAGGCCCAGTCCTAGCAAGTTGGCGGCGATGGTGCGCACATCGGCGAGCATCTGCACGTAGCTGATCGAACGCCAGGCACCGTCTGCCTGGCGTGCCGCGATGAAAGTGGTGTCCGGgcg                  >  2:932317/1‑150 (MQ=255)
                                                                                                                                  cgagcGCCCGCTGGCGCTCGGCACTGAGGCCCAGTCCTAGCAAGTTGGCGGCGATGGTGCGCACATCGGCGAGCATCTGCACGTAGCTGATCGAACGCCAGGCACCGTCTGCCTGGCGTGCCGCGATGAAAGTGGTGTCCGGGCGCACCt             >  2:759430/3‑150 (MQ=255)
                                                                                                                                   gagcGCCCGCTGGCGCTCGGCACTGAGGCCCAGTCCTAGCAAGTTGGCGGCGa                                                                                                             <  2:362240/52‑1 (MQ=255)
                                                                                                                                   gagcGCCCGCTGGCGCTCGGCACTGAGGCCCAGTCCTAGCAAGTTGGCGGCGa                                                                                                             >  1:362240/2‑53 (MQ=255)
                                                                                                                                   cagcGCCAGCGGGCGCTCGGCACTGAGGCCCAGTCCTAGCAACTTGGCGGCGATGGTGCGCACATCGGCGAGCATCTGCACGTAGCTGATCGAACGCCAGGCACCGTCTGCCTGGCGTGCCGCGATGAAAGTGGTGTCCGGGCGCACCTg            <  1:80903/150‑1 (MQ=255)
                                                                                                                                    agcGCCCGCTGGCGCTCGGCACTGAGGCCCAGTCCTAGCAAGTTGGCGGCGATGGTGCGCACATCGGCGAGCATCTGCACGTAGCTGATCGAACGCCAGGCACCGTCTGCCTGGCGTGCCGCGATGAAAGTGGTGTCCGGGCGCACCTgg           <  1:759430/150‑1 (MQ=255)
                                                                                                                                        ccAGCGGGCGCTCGGCACTGAGGCCCAGTCCTAGCAAGTTGGCGGCGATGGTGCGCACATCGGCGAGCATCTGCACGTAGCTGATCGAACGCCAGGCACCGTCTGCCTGGCGTGCCGCGATGAAAGTGGTGTCCGGGCGCACCTGGGccc       <  1:796838/150‑1 (MQ=255)
                                                                                                                                            cGGGCGCTCGGCACTGAGGCCCAGTCCTAGCAAGTTGGCGGCGATGGTGCGCACATCGGCGAGCATCTGCACGTAGCTGATCGAACG‑CAGGCACCGTCTGCCTGGCGTGCCGCGATGAAAGTGGTGTCCGGGCGCACCTGGGCCCAATGc  >  2:438420/1‑150 (MQ=255)
                                                                                                                                             |                                                                                                                                                     
TCGCAGACATGGCGCAACTTGGCGAAGTCTTGCGACAACAGCGCGTAGGCCGGCGACACCGGGCAATAGGCAATACCGGCATACATGGCGCCGAGGGCGATTTGCAGGTGTTCGATGTCGTTGCCGGAAAGCAGCGCCAGCGGGCGCTCGGCACTGAGGCCCAGTCCTAGCAAGTTGGCGGCGATGGTGCGCACATCGGCGAGCATCTGCACGTAGCTGATCGAACGCCAGGCACCGTCTGCCTGGCGTGCCGCGATGAAAGTGGTGTCCGGGCGCACCTGGGCCCAATGC  >  NC_002947/3796036‑3796326

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