Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,798,747 A→T 7.8% W66R (TGG→AGG)  PP_3358 ← hydroxycinnamoyl‑CoA hydratase‑lyase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,798,7470AT 7.8% 156.6 / 4.2 51W66R (TGG→AGG) PP_3358hydroxycinnamoyl‑CoA hydratase‑lyase
Reads supporting (aligned to +/- strand):  ref base A (20/27);  new base T (2/2);  total (22/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.87e-01

GCCGTTGACCATGGCGATGGTCGGTTTGGCATACAGACGCAGCAACTTCCATTGCCATTGCGAGGCTTCGCGACGAATCTTTTCCTGGAGGATTTCCGGGCCGGCGTCCACCTCGCGGAAGTACTCCTTCAGGTCCATGCCGGCG‑GTCCAGGACTCGCCGGCACCGGTCAATACCAGCACGCCAGCGTCAGCGTCCTGCTCAAGGGTTTCCAGCACGTCGACCATTTCCCGGTTCAGGGTGGGGCTCATGGCATTGCGTTTTTCCGGGCGATTGAGGGTCACCCAGGCGATGC  >  NC_002947/3798598‑3798890
                                                                                                                                                      |                                                                                                                                               
gCCGTTGACCATGGCGATGGTCGGTTTGGCATACAGACGCAGCAACTTCCATTGCCATTGCGAGGCTTCGCGACGAATCTTTTCCTGGAGGATTTCCGGGCCGGCGTCCACCTCGCGGAAGTACTCCTTCAGGTCCATGCCGGCG‑GTCCa                                                                                                                                                 <  1:1494645/150‑1 (MQ=255)
   gTTGACCATGGCGCTGGTCGGTTTGGCATACAGACGCAGCAACTTCCATTGCCATTGCGAGGCTTCGCGACGAATCTTTTCCTGGAGGATTTCCGGGCCGGCGTCCACCTCGCGGAAGTGCTCCTTCAGGTCCATGCCGGCG‑GTCCAGGa                                                                                                                                              <  2:1611918/150‑1 (MQ=255)
      gACCATGGCGATGGTCGGTTTGGCATACAGACGCAGCAACTTCCATTGCCATTGCGAGGCTTCGCGACGAATCTTTTCCTGGAGGATTTCCGGGCCGGCGTCCACCTCGCGGAAGTACTCCTTCAGGTCCATGCCGGCG‑GTCCAGGACTc                                                                                                                                           <  2:1401619/150‑1 (MQ=255)
       aCCATGGCGATGGTCGGTTTGGCATACAGACGCAGCAACTTCCATTGCCATTGCGAGGCTTCGCGACGAATCTTTTACTTGAGGATTTCCGGGCCGGCGTCCACCTCGCGGAAGTACTCCTTCAGGTCCATGCCGGCG‑GTCCAGGACTCg                                                                                                                                          >  1:364655/1‑150 (MQ=255)
         cATGGCGATGGTCGGTTTGGCATACAGACGCAGCAACTTCCATTGCCATTGCGAGGCTTCGCGACGAATCTTTTCCTGGAGGATTTCCGGGCCGGCGTCCACCTCGCGGAAGTACTCCTTCAGGTCCATGCCGGCG‑GTCCAGGACTCGcc                                                                                                                                        >  1:1049398/1‑150 (MQ=255)
          aTGGCGATGGTCGGTTTGGCATACAGACGCAGCAACTTCCATTGCCATTGCGAGGCTTCGCGACGAATCTTTTCCTGGAGGATTTCCGGGCCGGCGTCCACCTCGCGGAAGTACTCCTTCAGGTCCATGCCGGCG‑GTCCAGGACTCGCCg                                                                                                                                       >  1:1473021/1‑150 (MQ=255)
               gATGGTCGGTTTGGCATACAGACGCAGCAACTTCCATTGCCATTGCGAGGCTTCGCGACGAATCTTTTCCTGGAGGATTTCCGGGCCGGCGTCCACCTCGCGGAAGTACTCCTTCAGGTCCGTGCCGGCGAGTCCTGGAC‑CGCCgg                                                                                                                                      <  2:1416/146‑1 (MQ=255)
               gATGGTCGGTTTGGCATACAGACGCAGCAACTTCCATTGCCATTGCGAGGCTTCGCGACGAATCTTTTCCTGGAGGATTTCCGGGCCGGCGTCCACCTCGCGGAAGTACTCCTTCAGGTCCGTGCCGGCGAGTCCTGGAC‑CGCCgg                                                                                                                                      >  1:1416/1‑146 (MQ=255)
                           ggCATACAGACGCAGCAACTTCCATTGCCATTGCGAGGCTTCGCGACGAATCTTTTCCTGGAGGATTTCCGGGCCGGCGTCCACCTCGCGGAAGTACTCCTTCAGGTCCATGCCGGCG‑GTCCAGGACTCGCCGGCACCGGTCAATACCAg                                                                                                                      <  1:1095667/150‑1 (MQ=255)
                           ggCATACAGACGCAGCAACTTCCATTGCCATTGCGAGGCTTCGCGACGAATCTTTTCCTGGAGGATTTCCGGGCCGGCGTCCACCTCGCGGAAGTACTCCTTCAGGTCCATGCCGGCG‑GTCCAGGACTCGCCGGCACCGGTCAATACCAg                                                                                                                      <  1:1587299/150‑1 (MQ=255)
                            gCATACAGACGCAGCAACTTCCATTGCCATTGCGAGGCTTCGCGACGAATCTTTTCCTGGAGGATTTCCGGGCCGGCGTCCACCTCGCGGAAGTACTCCTTCAGGTCCATGCCGGCG‑GTCCAGGACTCGc                                                                                                                                         <  1:662800/130‑1 (MQ=255)
                            gCATACAGACGCAGCAACTTCCATTGCCATTGCGAGGCTTCGCGACGAATCTTTTCCTGGAGGATTTCCGGGCCGGCGTCCACCTCGCGGAAGTACTCCTTCAGGTCCATGCCGGCG‑GTCCAGGACTCGc                                                                                                                                         >  2:662800/1‑130 (MQ=255)
                              aTACAGACGCAGCAACTTCCATTGCCATTGCGAGGCTTCGCGACGAATCTTTTCCTGGAGGATTTCCGGGCCGGCGTCCACCTCGCGGAAGTACTCCTTCAGGTCCATGCCGGCG‑GTCCAGGACTCGCCGGCACCGGTCAATACCAGCAc                                                                                                                   >  2:623395/1‑150 (MQ=255)
                                   gACGCATCAACTTCCATTGCCATTGCGAGGCTTCGCGACGAATCTTTTCCTGGAGGATTTCCGGGCCGGCGTCCACCTCGCGGAAGTACTCCTTCAGGTCCATGCCGGCG‑GTCCAGGACTCGCCGGCACCGGTCAATACCAGCACGCcag                                                                                                              >  2:1002252/1‑150 (MQ=255)
                                   gACGCAGCAACTTCCATTGCCATTGCGAGGCTTCGCGACGAATCTTTTCCTGGAGGATTTCCGGGCCGGCGTCCACCTCGCGGAAGTACTCCTTCAGGTCCATGCCGGCG‑GTCCAGGa                                                                                                                                              >  1:1361145/1‑118 (MQ=255)
                                   gACGCAGCAACTTCCATTGCCATTGCGAGGCTTCGCGACGAATCTTTTCCTGGAGGATTTCCGGGCCGGCGTCCACCTCGCGGAAGTACTCCTTCAGGTCCATGCCGGCG‑GTCCAGGa                                                                                                                                              <  2:1361145/118‑1 (MQ=255)
                                        agcaACTTCCATTGCCATTGCGAGGCTTCGCGACGAATCTTTTCCTGGAGGATTTCCGGGCCGGCGTCCACCTCGCGGAAGTACTCCTTCAGGTCCATGCCGGCG‑GTCCAGGACTCGCCGGCACCGGTCAATACCAGCACGCcagcgtca                                                                                                         <  2:1018211/150‑1 (MQ=255)
                                             cTTCCATTGCCATTGCGAGGCTTCGCGACGAATCTTTTCCTGGAGGATTTCCGGGCCGGCGTCCACCTCGCGGAAGTACTCCTTCAGGTCCATGCCGGCG‑GTCCAGGACTCGCCGGCACCGGTCaa                                                                                                                            <  1:175740/126‑1 (MQ=255)
                                             cTTCCATTGCCATTGCGAGGCTTCGCGACGAATCTTTTCCTGGAGGATTTCCGGGCCGGCGTCCACCTCGCGGAAGTACTCCTTCAGGTCCATGCCGGCG‑GTCCAGGACTCGCCGGCACCGGTCaa                                                                                                                            >  2:175740/1‑126 (MQ=255)
                                                       cattgcGAGGCTTCGCGACGAATCTTTTCCTGGAGGATTTCCGGGCCGGCGTCCACCTCGCGGAAGTACTCCTTCAGGTCCATGCCGGCG‑GTCCAGGACTCGCCGGCACCGGTCAATACCAGCACGCCAGCGTCAGCGTCCTGCTCAAgg                                                                                          <  1:656230/150‑1 (MQ=255)
                                                            cGAGGCTTCGCGACTAATCTTTTCCTGGAGGATTTCCGGGCCGGCGTCCACCTCGCGGAAGTACTCCTTCAGGTCCATGCCGGCG‑GTCCAGGACTCGCCGGCACCGGTCAATACCAGCACGCCAGCGTCAGCGTCCTGCTCAAGGGTTTc                                                                                     <  2:860525/150‑1 (MQ=255)
                                                                 cTTCGCGACGAATCTTTTCCTGGAGGATTTCCGGGCCGGCGTCCACCTCGCGGAAGTACTCCTTCAGGTCCATGCCGGCG‑GTCCAGGACTCGCCGGCACCGGTCAATACCAGCACGCCAGCGTCAGCGTCCTGCTCAAGGGTTTCCAGCa                                                                                >  1:1421836/1‑150 (MQ=255)
                                                                        acgaATCTTTTCCTGGAGGATTTCCGGGCCGGCGTCCACCTCGCGGAAGTACTCCTTCAGGTCCATGCCGGCG‑GTCCAGGACTCGCCGGCACCGGTCAATACCAGCACGCCAGCGTCAGCGTCCTGCTCAAGGGTTTCCAGCACGTCGAc                                                                         <  1:623395/150‑1 (MQ=255)
                                                                        acgaATCTTCTCCTTGAGGATTTCCGGGCCGGCGTCCACCTCGCGGAAGTACTCAATTTGCTGGCGAACGGCT‑GTCCAGGCCCCGCCGGCACCGGTCAATACCAGCACGCCAGCGTCAGCGTCCTGCTCAAGGGTTTCCAGCACGTCGAc                                                                         <  1:529447/150‑1 (MQ=255)
                                                                               ttttCCTGGAGGATTTCCGGGCCGGCGTCCACCTCGCGGAAGTACTCCTTCAGGTCCATGCCGGCG‑GTCCAGGACTCACCGGCACCGGTCAATACCAGCACGCCAGCGTCAGCGTCCTGCTCAAGGGTTTCCAGCACGTCGACCATTTcc                                                                  >  2:266701/1‑150 (MQ=255)
                                                                                tttCCTTTAGGATGTCCTGGTCGGCGTCCACCTCGCGGAAGTACTCCTTCAGGTCCATGCCGGCG‑GTCCAGGACTCGCCGGCACCGGTCAATACCAGCACGCCAGCGTCAGCGTCCTGCTCAAGGGTTTCCAGCACGTCGACCATTTccc                                                                 <  2:721445/150‑1 (MQ=255)
                                                                                     tGGAGGATTTCCGGGCCGGCGTCCACCTCGCGGAAGTACTCCTTCAGGTCCATGCCGGCGAGTCCTGGAC‑CGCCGGCat                                                                                                                                   >  1:598304/1‑78 (MQ=37)
                                                                                     tGGAGGATTTCCGGGCCGGCGTCCACCTCGCGGAAGTACTCCTTCAGGTCCATGCCGGCGAGTCCTGGAC‑CGCCGGCat                                                                                                                                   <  2:598304/79‑2 (MQ=37)
                                                                                       gaggaTTTCCGGGCCGGCGTCCACCTCGCGGAAGTACTCCTTCAGGTCCATGCCGGCG‑GTCCAGGACTCGCCGGCACCGGTCAATACCAGCACGCCAGCGTCAGCGTCCTGCTCAAGGGTTTCCAGCACGTCGACCATTTCCCGGTTCAg                                                          >  1:1369311/1‑150 (MQ=255)
                                                                                       gaggaTTTCCGGGCCGGCGTCCACCTCGCGGAAGTACTCCTTCAGGTCCATGCCGGCG‑GTCCAGGACTCACCGGCACCGGTCAATACCAGCACGCCAGCGTCAGCGTCCTGCTCAAGGGTTTCCAGCACGTCGACCATTTCCCGGTTCAg                                                          <  1:266701/150‑1 (MQ=255)
                                                                                        aggaTTTCCGGGCCGGCGTCCACCTCGCGGAAGTACTCCTTCAGGTCCATGCCGGCG‑GTCCAGGACTCGCCGGCACCGGTCAATACCAGCACGCCAGCGTCAGCGTCCTGCTCAAGGGTTTCCAGCACGTCGACCATTTCCCGGTTCAgg                                                         <  2:1336237/150‑1 (MQ=255)
                                                                                         ggaTTTCCGGGCCGGCGTCCACCTCGCGGAAGTACTCCTTCAGGTCCATGCCGGCG‑GTCCAGGACTCGCCGGCACCGGTCAATACCAGCACGCCAGCGTCAGCGTCCTGCTCAAGGGTTTCCAGCACGTCGACCATTTCCCGGTTCAggg                                                        <  2:1049398/150‑1 (MQ=255)
                                                                                           aTTTCCGGGCCGGCGTCCACCTCGCGGAAGTACTCCTTCAGGTCCATGCCGGCG‑GTCCAGGACTCGCCGGCACCGGTCAATACCAGCACGCCAGCGTCAGCGTCCTGCTCAAGGGTTTCCAGCACGTCGACCATTTCCCGGTTCAGGGTg                                                      <  2:364655/150‑1 (MQ=255)
                                                                                           aTTTCCGGGCCGGCGTCCACCTCGCGGAAGTACTCCTTCAGGTCCATGCCGGCG‑GTCCAGGACTCGCCGGCACCGGTCAATACCAGCACGCCAGCGTCAGCGTCCTGCTCAAGGGTTTCCAGCACGTCGACCATTTCCCGGTTCAGGGTg                                                      <  2:1473021/150‑1 (MQ=255)
                                                                                           aTTTCCGGGCCGGCGTCCACCTCGCGGAAGTACTCCTTCAGGTCCATGCCGGCG‑GTCCAGGACTCGCCGGCACCGGTCAATACCAGCACGCCAGCGTCAGCGTCCTGCTCAAGGGTTTCCAGCACGTCGACCATTTCCCGGTTCAGGGTg                                                      <  1:357368/150‑1 (MQ=255)
                                                                                                               cTCGCGGAAGTACTCCTTCAGGTCCATGCCGGCG‑GTCCAGGACTCGCCGGCACCGGTCAATACCAGCACGCCAGCGTCAGCGTCCTGCTCAAGGGTTTCCAGCACGTCGACCATTTCCCGGTTCAGGGTGGGGCTCATGGCATTGCGttt                                  >  1:123104/1‑150 (MQ=255)
                                                                                                               cTCGCGGAAGTACTCCTTCAGGTCCATGCCGGCG‑GTCCAGGACTCGCCGGCACCGGTCAATACCAGCACGCCAGCGTCAGCGTCCTGCTCAAGGGTTTCCAGCACGTCGACCATTTCCCGGTTCAGGGTGGGGCTCATGGCATTGCGttt                                  >  1:891730/1‑150 (MQ=255)
                                                                                                               cTCGCGGAAGTACTCCTTCAGGTCCATGCCGGCG‑GTCCAGGACTCGCCGGCACCGGTCAATACCAGCACGCCAGCGTCAGCGTCCTGCTCAAGGGTTTCCAGCACGTCGACCATTTCCCGGTTCAGGGTGGGGCTCATGGCATTGCGttt                                  >  2:264684/1‑150 (MQ=255)
                                                                                                                    ggAAGTACTCCTTCAGGTCCATGCCGGCG‑GTCCAGGACTCGCCGGCACCGGTCAATACCAGCACGCca                                                                                                               <  2:437936/68‑1 (MQ=255)
                                                                                                                    ggAAGTACTCCTTCAGGTCCATGCCGGCG‑GTCCAGGACTCGCCGGCACCGGTCAATACCAGCACGCca                                                                                                               >  1:437936/1‑68 (MQ=255)
                                                                                                                         tACTCCTTCAGGTCCATGCCGGCG‑GTCCAGGACTCGCCGGCACCGGTCAATACCAGCACGCCAGCGTCAGCGTCCTGCTCAAGGGTTTCCAGCACGTCGACCATTTCCCgg                                                               >  1:1566097/1‑111 (MQ=255)
                                                                                                                         tACTCCTTCAGGTCCATGCCGGCG‑GTCCAGGACTCGCCGGCACCGGTCAATACCAGCACGCCAGCGTCAGCGTCCTGCTCAAGGGTTTCCAGCACGTCGACCATTTCCCgg                                                               <  2:1566097/111‑1 (MQ=255)
                                                                                                                             ccTTCAGGTCCATGCCGGCG‑GTCCAGGACTCGCCGGCACCGGTCAATACCAGCACGCCAGCGTCAGCGTCCTGCTCAAGGGTTTCCAGCACGTCGACCATTTCCCGGTTCAGGGTGGGGCTCATGGCATTGCGTTTTTCCGGGCGATTGa                    <  2:891730/150‑1 (MQ=255)
                                                                                                                              cTTCAGGTCCATGCCGGCG‑GTCCAGGACTCGCCGGCACCGGTCAATACCAGCACGCCAGCGTCAGCGTCCTGCTCAAGGGTTTCCAGCACGTCGAc                                                                         <  2:525449/96‑1 (MQ=255)
                                                                                                                              cTTCAGGTCCATGCCGGCG‑GTCCAGGACTCGCCGGCACCGGTCAATACCAGCACGCCAGCGTCAGCGTCCTGCTCAAGGGTTTCCAGCACGTCGAc                                                                         >  1:525449/1‑96 (MQ=255)
                                                                                                                               ttCAGGTCCATGCCGGCG‑GTCCAGGACTCGCCGGCACCGGTCAATACCAGCACGCCAGCGTCAGCGTCCTGCTCAAGGGTTTCCAGCACGTCGACCATTTCCCGGTTCAGGGTGGGGCTCATGGCATTGCGTTTTTCCGGGCGATTGAgg                  <  2:1369311/150‑1 (MQ=255)
                                                                                                                                tCAGGTCCATGCCGGCG‑GTCCAGGACTCGCCGGCACCGGTCAATACCAGCACGCCAGCGTCAGCGTCCTGCTCAAGGGTTTCCAGCACGTCGACCATTTCCCGGTTCAGGGTGGGGCTCATGGCATTGCGTTTTTCCGGGCGATTGAggg                 >  1:627753/1‑150 (MQ=255)
                                                                                                                                    gTCCATGCCGGCG‑GGCCAGGACTCGCCGGCACCGGTCAATACCAGCACGCCAGCGTCAGCGTCCTGCTCAAGGGTTTCCAGCACGTCGACCATTTCCCGGTTCAGGGTGGGGCTCATGGCATTGCGTTTTTCCGGGCGATTGAGGGTCAc             >  1:1542911/1‑150 (MQ=255)
                                                                                                                                          gCCGGCG‑GTCCAGGACTCGCCGGCACCGGTCAATACCAGCACGCCAGCGTCAGCGTCCTGCTCAAGGGTTTCCAGCACGTCGACCATTTCCCGGTTCAGGGTGGGGCTCATGGCATTGCGTTTTTCCGGGCGATTGAGGGTCACCCAGGc       <  1:1002252/150‑1 (MQ=255)
                                                                                                                                               cg‑gTCCAGGACTCGCCGGCACCGGTCAATACCAGCACGCCAGCGTCAGCGTCCTGCTCAAGGGTTTCCAGCACGTCGACCATTTCCCGGTTCAGGGTGGGGCTCATGGCATTGCGTTTTTCCGGGCGATTGAGGGTCACCCAGGCGATGc  <  1:264684/150‑1 (MQ=255)
                                                                                                                                                     cAGGACTCGCCGGCACCGGTCAATACCAGCACGCCAGCGTCAGCGTCCTGCTCAAGGGTTTCCAGCACGTCGACCATTTCCCgg                                                               <  2:355522/84‑1 (MQ=255)
                                                                                                                                                     cAGGACTCGCCGGCACCGGTCAATACCAGCACGCCAGCGTCAGCGTCCTGCTCAAGGGTTTCCAGCACGTCGACCATTTCCCgg                                                               >  1:355522/1‑84 (MQ=255)
                                                                                                                                                      |                                                                                                                                               
GCCGTTGACCATGGCGATGGTCGGTTTGGCATACAGACGCAGCAACTTCCATTGCCATTGCGAGGCTTCGCGACGAATCTTTTCCTGGAGGATTTCCGGGCCGGCGTCCACCTCGCGGAAGTACTCCTTCAGGTCCATGCCGGCG‑GTCCAGGACTCGCCGGCACCGGTCAATACCAGCACGCCAGCGTCAGCGTCCTGCTCAAGGGTTTCCAGCACGTCGACCATTTCCCGGTTCAGGGTGGGGCTCATGGCATTGCGTTTTTCCGGGCGATTGAGGGTCACCCAGGCGATGC  >  NC_002947/3798598‑3798890

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