Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,798,752	Δ1 bp	7.8%	coding (191/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,752	0	T	.	7.8%	157.6 / 7.9	51	coding (191/831 nt)	PP_3358	hydroxycinnamoyl‑CoA hydratase‑lyase
Reads supporting (aligned to +/- strand):  ref base T (20/26);  new base . (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 = 6.06e-01

GACCATGGCGATGGTCGGTTTGGCATACAGACGCAGCAACTTCCATTGCCATTGCGAGGCTTCGCGACGAATCTTTTCCTGGAGGATTTCCGGGCCGGCGTCCACCTCGCGGAAGTACTCCTTCAGGTCCATGCCGGCG‑GTCCAGGACTCGCCGGCACCGGTCAATACCAGCACGCCAGCGTCAGCGTCCTGCTCAAGGGTTTCCAGCACGTCGACCATTTCCCGGTTCAGGGTGGGGCTCATGGCATTGCGTTTTTCCGGGCGATTGAGGGTCACCCAGGCGATGCCCGCTTCCAGT  >  NC_002947/3798604‑3798901                                                                                                                                                      |                                                                                                                                                      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)                                   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                                                                 <  1:357368/150‑1 (MQ=255)                                                                                      aTTTCCGGGCCGGCGTCCACCTCGCGGAAGTACTCCTTCAGGTCCATGCCGGCG‑GTCCAGGACTCGCCGGCACCGGTCAATACCAGCACGCCAGCGTCAGCGTCCTGCTCAAGGGTTTCCAGCACGTCGACCATTTCCCGGTTCAGGGTg                                                                 <  2:1473021/150‑1 (MQ=255)                                                                                                          cTCGCGGAAGTACTCCTTCAGGTCCATGCCGGCG‑GTCCAGGACTCGCCGGCACCGGTCAATACCAGCACGCCAGCGTCAGCGTCCTGCTCAAGGGTTTCCAGCACGTCGACCATTTCCCGGTTCAGGGTGGGGCTCATGGCATTGCGttt                                             >  1:891730/1‑150 (MQ=255)                                                                                                          cTCGCGGAAGTACTCCTTCAGGTCCATGCCGGCG‑GTCCAGGACTCGCCGGCACCGGTCAATACCAGCACGCCAGCGTCAGCGTCCTGCTCAAGGGTTTCCAGCACGTCGACCATTTCCCGGTTCAGGGTGGGGCTCATGGCATTGCGttt                                             >  2:264684/1‑150 (MQ=255)                                                                                                          cTCGCGGAAGTACTCCTTCAGGTCCATGCCGGCG‑GTCCAGGACTCGCCGGCACCGGTCAATACCAGCACGCCAGCGTCAGCGTCCTGCTCAAGGGTTTCCAGCACGTCGACCATTTCCCGGTTCAGGGTGGGGCTCATGGCATTGCGttt                                             >  1:123104/1‑150 (MQ=255)                                                                                                               ggAAGTACTCCTTCAGGTCCATGCCGGCG‑GTCCAGGACTCGCCGGCACCGGTCAATACCAGCACGCca                                                                                                                          >  1:437936/1‑68 (MQ=255)                                                                                                               ggAAGTACTCCTTCAGGTCCATGCCGGCG‑GTCCAGGACTCGCCGGCACCGGTCAATACCAGCACGCca                                                                                                                          <  2:437936/68‑1 (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                                                                                    >  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)                                                                                                                                                cAGGACTCGCCGGCACCGGTCAATACCAGCACGCCAGCGTCAGCGTCCTGCTCAAGGGTTTCCAGCACGTCGACCATTTCCCgg                                                                          <  2:355522/84‑1 (MQ=255)                                                                                                                                                cAGGACTCGCCGGCACCGGTCAATACCAGCACGCCAGCGTCAGCGTCCTGCTCAAGGGTTTCCAGCACGTCGACCATTTCCCgg                                                                          >  1:355522/1‑84 (MQ=255)                                                                                                                                                    aCTCGCCGGCACCGGTCAATACCAGCACGCCAGCGTCAGCGTCCTGCTCAAGGGTTTCCAGCACGTCGACCATTTCCCGGTTCAGGGTGGGGCTCATGGCATTGCGTTTTTCCGGGCGATTGAgg                             >  2:1226854/1‑125 (MQ=255)                                                                                                                                                    aCTCGCCGGCACCGGTCAATACCAGCACGCCAGCGTCAGCGTCCTGCTCAAGGGTTTCCAGCACGTCGACCATTTCCCGGTTCAGGGTGGGGCTCATGGCATTGCGTTTTTCCGGGCGATTGAgg                             <  1:1226854/125‑1 (MQ=255)                                                                                                                                                    aCTCGCCGGCACCGGTCAATACCAGCACGCCAGCGTCAGCGTCCTGCTCAAGGGTTTCCAGCACGTCGACCATTTCCCGGTTCAGGGTGGGGCTCATGGCATTGCGTTTTTCCGGGCGATTGAGGGTCACCCAGGCGATGCCCGCTTCCa    <  2:1542911/150‑1 (MQ=255)                                                                                                                                                      tCGCCGGCACCGGTCAATACCAGCACGCCAGCGTCAGCGTCCTGCTCAAGGGTTTCCAGCACGTCGACCATTTCCCGGTTCAGGGTGGGGCTCATGGCATTGCGTTTTTCCGGGCGATTGAGGGTCACCCAGGCGATGCCCGCTTCCAGt  >  2:1344475/1‑150 (MQ=255)                                                                                                                                                      tCGCCGGCACCGGTCAATACCAGCACGCCAGCGTCAGCGTCCTGCTCAAGGGTTTCCAGCACGTCGACCATTTCCCGGTTCAGGGTGGGGCTCATGGCATTGCGTTTTTCCGGGCGATTGAGGGTCACCCAGGCGATGCCCGCTTCCAGt  >  1:1578504/1‑150 (MQ=255)                                                                                                                                                      |                                                                                                                                                      GACCATGGCGATGGTCGGTTTGGCATACAGACGCAGCAACTTCCATTGCCATTGCGAGGCTTCGCGACGAATCTTTTCCTGGAGGATTTCCGGGCCGGCGTCCACCTCGCGGAAGTACTCCTTCAGGTCCATGCCGGCG‑GTCCAGGACTCGCCGGCACCGGTCAATACCAGCACGCCAGCGTCAGCGTCCTGCTCAAGGGTTTCCAGCACGTCGACCATTTCCCGGTTCAGGGTGGGGCTCATGGCATTGCGTTTTTCCGGGCGATTGAGGGTCACCCAGGCGATGCCCGCTTCCAGT  >  NC_002947/3798604‑3798901

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: ‑