Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,798,742:1	+A	8.0%	coding (201/831 nt)	PP_3358 ←	hydroxycinnamoyl‑CoA hydratase‑lyase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,798,742	1	.	A	8.0%	191.2 / 11.6	50	coding (201/831 nt)	PP_3358	hydroxycinnamoyl‑CoA hydratase‑lyase
Reads supporting (aligned to +/- strand):  ref base . (20/26);  new base A (2/2);  total (22/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.06e-01

ACCAGCCGTTGACCATGGCGATGGTCGGTTTGGCATACAGACGCAGCAACTTCCATTGCCATTGCGAGGCTTCGCGACGAATCTTTTCCTGGAGGATTTCCGGGCCGGCGTCCACCTCGCGGAAGTACTCCTTCAGGTCCATGCCGGCG‑GTCCAGGACTCGCCGGCACCGGTCAATACCAGCACGCCAGCGTCAGCGTCCTGCTCAAGGGTTTCCAGCACGTCGACCATTTCCCGGTTCAGGGTGGGGCTCATGGCATTGCGTTTTTCCGGGCGATTGAGGGTCACCCAGGCGATGC  >  NC_002947/3798594‑3798890                                                                                                                                                      |                                                                                                                                                     aCCAGCCGTTGACCATGGCGATGGTCGGTTTGGCATACAGNCGCAGCNACTTCCATTGCCATTGCGAGGCTTCGCGACGAACCGTTTCCTGGAGGATTTCCGGGCCGGCGTCCACCTCGCGGAAGTACTCCTTCAGGTCCATGCcggcg‑g                                                                                                                                                     >  1:1336237/1‑150 (MQ=255)  ccAGCCGTTGACCATGGCGATGGTCGGTTTGGCATACAGACGCAGCAACTTCCATTGCCATTGCGAGGCTTCGCGACGAATCTTTTCCTGGAGGATTTCCGGGCCGGCGTCCACCTCGCGGAAGTACTCCTTCAGGTCCATACCGGCG‑Gt                                                                                                                                                    >  2:529447/1‑150 (MQ=255)     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                                                                                                                                      >  1:1416/1‑146 (MQ=255)                    gATGGTCGGTTTGGCATACAGACGCAGCAACTTCCATTGCCATTGCGAGGCTTCGCGACGAATCTTTTCCTGGAGGATTTCCGGGCCGGCGTCCACCTCGCGGAAGTACTCCTTCAGGTCCGTGCCGGCGAGTCCTGGAC‑CGCCgg                                                                                                                                      <  2:1416/146‑1 (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                                                                                                                                              <  2:1361145/118‑1 (MQ=255)                                        gACGCAGCAACTTCCATTGCCATTGCGAGGCTTCGCGACGAATCTTTTCCTGGAGGATTTCCGGGCCGGCGTCCACCTCGCGGAAGTACTCCTTCAGGTCCATGCCGGCG‑GTCCAGGa                                                                                                                                              >  1:1361145/1‑118 (MQ=255)                                             agcaACTTCCATTGCCATTGCGAGGCTTCGCGACGAATCTTTTCCTGGAGGATTTCCGGGCCGGCGTCCACCTCGCGGAAGTACTCCTTCAGGTCCATGCCGGCG‑GTCCAGGACTCGCCGGCACCGGTCAATACCAGCACGCcagcgtca                                                                                                         <  2:1018211/150‑1 (MQ=255)                                                  cTTCCATTGCCATTGCGAGGCTTCGCGACGAATCTTTTCCTGGAGGATTTCCGGGCCGGCGTCCACCTCGCGGAAGTACTCCTTCAGGTCCATGCCGGCG‑GTCCAGGACTCGCCGGCACCGGTCaa                                                                                                                            >  2:175740/1‑126 (MQ=255)                                                  cTTCCATTGCCATTGCGAGGCTTCGCGACGAATCTTTTCCTGGAGGATTTCCGGGCCGGCGTCCACCTCGCGGAAGTACTCCTTCAGGTCCATGCCGGCG‑GTCCAGGACTCGCCGGCACCGGTCaa                                                                                                                            <  1:175740/126‑1 (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                                                                                                                                   <  2:598304/79‑2 (MQ=37)                                                                                          tGGAGGATTTCCGGGCCGGCGTCCACCTCGCGGAAGTACTCCTTCAGGTCCATGCCGGCGAGTCCTGGAC‑CGCCGGCat                                                                                                                                   >  1:598304/1‑78 (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:1473021/150‑1 (MQ=255)                                                                                                aTTTCCGGGCCGGCGTCCACCTCGCGGAAGTACTCCTTCAGGTCCATGCCGGCG‑GTCCAGGACTCGCCGGCACCGGTCAATACCAGCACGCCAGCGTCAGCGTCCTGCTCAAGGGTTTCCAGCACGTCGACCATTTCCCGGTTCAGGGTg                                                      <  2:364655/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                                                                                                               >  1:437936/1‑68 (MQ=255)                                                                                                                         ggAAGTACTCCTTCAGGTCCATGCCGGCG‑GTCCAGGACTCGCCGGCACCGGTCAATACCAGCACGCca                                                                                                               <  2:437936/68‑1 (MQ=255)                                                                                                                              tACTCCTTCAGGTCCATGCCGGCG‑GTCCAGGACTCGCCGGCACCGGTCAATACCAGCACGCCAGCGTCAGCGTCCTGCTCAAGGGTTTCCAGCACGTCGACCATTTCCCgg                                                               <  2:1566097/111‑1 (MQ=255)                                                                                                                              tACTCCTTCAGGTCCATGCCGGCG‑GTCCAGGACTCGCCGGCACCGGTCAATACCAGCACGCCAGCGTCAGCGTCCTGCTCAAGGGTTTCCAGCACGTCGACCATTTCCCgg                                                               >  1:1566097/1‑111 (MQ=255)                                                                                                                                  ccTTCAGGTCCATGCCGGCG‑GTCCAGGACTCGCCGGCACCGGTCAATACCAGCACGCCAGCGTCAGCGTCCTGCTCAAGGGTTTCCAGCACGTCGACCATTTCCCGGTTCAGGGTGGGGCTCATGGCATTGCGTTTTTCCGGGCGATTGa                    <  2:891730/150‑1 (MQ=255)                                                                                                                                   cTTCAGGTCCATGCCGGCG‑GTCCAGGACTCGCCGGCACCGGTCAATACCAGCACGCCAGCGTCAGCGTCCTGCTCAAGGGTTTCCAGCACGTCGAc                                                                         >  1:525449/1‑96 (MQ=255)                                                                                                                                   cTTCAGGTCCATGCCGGCG‑GTCCAGGACTCGCCGGCACCGGTCAATACCAGCACGCCAGCGTCAGCGTCCTGCTCAAGGGTTTCCAGCACGTCGAc                                                                         <  2:525449/96‑1 (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)                                                                                                                                                      |                                                                                                                                                     ACCAGCCGTTGACCATGGCGATGGTCGGTTTGGCATACAGACGCAGCAACTTCCATTGCCATTGCGAGGCTTCGCGACGAATCTTTTCCTGGAGGATTTCCGGGCCGGCGTCCACCTCGCGGAAGTACTCCTTCAGGTCCATGCCGGCG‑GTCCAGGACTCGCCGGCACCGGTCAATACCAGCACGCCAGCGTCAGCGTCCTGCTCAAGGGTTTCCAGCACGTCGACCATTTCCCGGTTCAGGGTGGGGCTCATGGCATTGCGTTTTTCCGGGCGATTGAGGGTCACCCAGGCGATGC  >  NC_002947/3798594‑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: ‑