Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NZ_CP009273	3,223,174	A→T	I310F (ATC→TTC)	ygjK →	alpha‑glucosidase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NZ_CP009273	3,223,174	0	A	T	93.5%	222.0 / 5.8	62	I310F (ATC→TTC)	ygjK	alpha‑glucosidase
Reads supporting (aligned to +/- strand):  ref base A (2/2);  new base T (30/28);  total (32/30)
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.35e-01

CAGATCCGCGATATTCTGGCGCGTCCGGCGTTTTATCTCACCGCCTCGCAGCAACGCTGGGAAGAATACCTGAAGAAAGGGTTAACCAATCCGGATGCGACGCCGGAACAGACGCGCGTCGCGGTGAAAGCGATCGAAACGCTCAACGGTAACTGGCGCTCACCGGGCGGTGCGGTGAAATTTAACACCGTTACACCGTCGGTGACCGGGCGCTGGTTCTCCGGCAATCAGACCTGGCCGTGGGATACCTGGAAGCAGGCGT  >  NZ_CP009273/3223042‑3223303                                                                                                                                     |                                                                                                                                  cAGATCCGCGATATTCTGGCGCGTCCGGCGTTTTATCTCACCGCCTCGCAGCAACGCTGGGAAGAATACCTGAAGAAAGGGTTAACCAATCCGGATGCGACGCCGGAACAGACGCGCGTCGCGGTGAAAGCGTTCGaa                                                                                                                              <  2:1007259/138‑1 (MQ=255)       cgcgATATTCTGGCGCGTCCGGCGTTTTATCTCACCGCCTCGCAGCAACGCTGGGAAGAATACCTGAAGAAAGGGTTAACCAATCCGGATGCGACGCCGGAACAGACGCGCGTCGCGGTGAAAGCGTTCGAAACGCTCa                                                                                                                       <  2:1093603/139‑1 (MQ=255)                   gcgcGTCCGGCGTTTTATCTCACCGCCTCGCAGCAACGCTGGGAAGAATACCTGAAGAAAGGGTTAACCAATCCGGATGCGACGCCGGAACAGACGCGCGTCGCGGTGAAAGCGTTCGAAACGCTCAACGGTAACTGgc                                                                                                           <  1:1283416/139‑1 (MQ=255)                      cgTCCGGCGTTTTATCTCACCGCCTCGCAGCAACGCTGGGAAGAATACCTGAAGAAAGGGTTAACCAATCCGGATGCGACGCCGGAACAGACGCGCGTCGCGGTGAAAGCGTTCGAAACGCTCAACGGTAACTGGCGCt                                                                                                        >  2:379317/1‑139 (MQ=255)                        tCCGGCGTTTTATCTCACCGCCTCGCAGCAACGCTGGGAAGAATACCTGAAGAAAGGGTTAACCAATCCGGATGCGACGCCGGAACAGACGCGCGTCGCGGTGAAAGCGTTCGAAACGCTCAACGGTAACTGGCGCTCa                                                                                                      >  1:1218023/1‑139 (MQ=255)                        tCCGGCGTTTTATCTCACCGCCTCGCAGCAACGCTGGGAAGAATACCTGAAGAAAGGGTTAACCAATCCGGATGCGACGCCGGAACAGACGCGCGTCGCGGTGAAAGCGTTCGAAACGCTCAACGGTAACTGGCGCTCa                                                                                                      >  1:1219010/1‑139 (MQ=255)                        tCCGGCGTTTTATCTCACCGCCTCGCAGCAACGCTGGGAAGAATACCTGAAGAAAGGGCTAACCAATCCGGATGCGACGCCGGAACAGACGCGCGTCGCGGTGAAAGCGTTCGAAACGCTCAACGGTAACTGGCGCTCa                                                                                                      >  2:207240/1‑139 (MQ=255)                             cGTTTTATCTCACCGCCTCGCAGCAACGCTGGGAAGAATACCTGAAGAAAGGGTTAACCAATCCGGATGCGACGCCGGAACAGACGCGCGTCGCGGTGAAAGCGTTCGAAACGCTCAACGGTAACTGGCGCTCACCggg                                                                                                 <  1:1071879/139‑1 (MQ=255)                             cGTTTTATCTCACCGCCTCGCAGCAACGCTGGGAAGAATACCTGAAGAAAGGGTTAACCAATCCGGATGCGACGCCGGAACAGACGCGCGTCGCGGTGAAAGCGTTCGAAACGCTCAACGGTAACTGGCGCTCACCggg                                                                                                 <  1:1231531/139‑1 (MQ=255)                                       cACCGCCTCGCAGCAACGCTGGGAAGAATACCTGAAGAAAGGGTTAACCAATCCGGATGCGACGCCGGAACAGACGCGCGTCGCGGTGAAAGCGTTCGAAACGCTCAACGGTAACTGGCGCTCACCGGgcggtgcggt                                                                                        >  1:925585/1‑138 (MQ=255)                                        aCCGCCTCGCAGCAACGCTGGGAAGAATACCTGAAGAAAGGGTTAACCAATCCGGATGCGACGCCGGAACAGACGCGCGTCGCGGTGAAAGCGTTCGAAACGCTCAACGGTAACTGGCGCTCACCGGGcggtgcggtg                                                                                       <  2:760102/138‑1 (MQ=255)                                        aCCGCCTCGCAGCAACGCTGGGAAGAATACCTGAAGAAAGGGTTAACCAATCCGGATGCGACGCCGGAACAGACGCGCGTCGCGGTGAAAGCGTTCGAAACGCTCAACGGTAACTGGCGCTCACCGGGcggtgcggtg                                                                                       <  2:657186/138‑1 (MQ=255)                                            ccTCGCAGCAACGCTGGGAAGAATACCTGAAGAAAGGGTTAACCAATCCGGATGCGACGCCGGAACAGACGCGCGTCGCGGTGAAAGCGTTCGAAACGCTCAACGGTAACTGGCGCTCACCGGGCGGTGCGGTGAAAtt                                                                                  >  2:519800/1‑139 (MQ=255)                                                gcagcaACGCTGGGAAGAATACCTGAAGAAAGGGTTAACCAATCCGGATGCGACGCCGGAACAGACGCGCGTCGCGGTGAAAGCGTTCGAAACGCTCAACGGTAACTGGCGCTCACCGGGCGGTGCGGTGAAATTTaa                                                                               <  2:928572/138‑1 (MQ=255)                                                gcagcaACGCTGGGAAGAATACCTGAAGAAAGGGTTAACCAATCCGGATGCGACGCCGGAACAGACGCGCGTCGCGGTGAAAGCGTTCGAAACGCTCAACGGTAACTGGCGCTCACCGGGCGGTGCGGTGAAATTTaa                                                                               <  1:90371/138‑1 (MQ=255)                                                gcagcaACGCTGGGAAGAATACCTGAAGAAAGGGTTAACCAATCCGGATGCGACGCCGGAACAGACGCGCGTCGCGGTGAAAGCGTTCGAAACGCTCAACGGTAACTGGCGCTCACCGGGCGGTGCGGTGAAATTTaa                                                                               <  1:519800/138‑1 (MQ=255)                                                gcagcaACGCTGGGAAGAATACCTGAAGAAAGGGTTAACCAATCCGGATGCGACGCCGGAACAGACGCGCGTCGCGGTGAAAGCGTTCGAAACGCTCAACGGTAACTGGCGCTCACCGGGCGGTGCGGTGAAATTTaa                                                                               <  1:521280/138‑1 (MQ=255)                                                gcagcaACGCTGGGAAGAATACCTGAAGAAAGGGTTAACCAATCCGGATGCGACGCCGGAACAGACGCGCGTCGCGGTGAAAGCGTTCGAAACGCTCAACGGTAACTGGCGCTCACCGGGCGGTGCGGTGAAATTTaa                                                                               >  2:675870/1‑138 (MQ=255)                                                gcagcaACGCTGGGAAGAATACCTGAAGAAAGGGTTAACCAATCCGGATGCGACGCCGGAACAGACGCGCGTCGCGGTGAAAGCGTTCGAAACGCTCAACGGTAACTGGCGCTCACCGGGCGGTGCGGTGAAATTTaa                                                                               <  2:870634/138‑1 (MQ=255)                                                gcagcaACGCTGGGAAGAATACCTGAAGAAAGGGTTAACCAATCCGGATGCGACGCCGGAACAGACGCGCGTCGCGGTGAAAGCGTTCGAAACGCTCAACGGTAACTGGCGCTCACCGGGCGGTGCGGTGAAATTTaa                                                                               <  2:978010/138‑1 (MQ=255)                                                         cTGGGAAGAATACCTGAAGAAAGGGTTAACCAATCCGGATGCGACGCCGGAACAGACGCGCGTCGCGGTGAAAGCGTTCGAAACGCTCAACGGTAACTGGCGCTCACCGGGCGGTGCGGTGAAATTTAACACCGTTaca                                                                     <  2:925585/139‑1 (MQ=255)                                                         cTGGGAAGAATACCTGAAGAAAGGGTTAACCAATCCGGATGCGACGCCGGAACAGACGCGCGTCGCGGTGAAAGCGTTCGAAACGCTCAACGGTAACTGGCGCTCACCGGGCGGTGCGGTGAAATTTAACACCGTTaca                                                                     <  2:1222443/139‑1 (MQ=255)                                                            ggAAGAATACCTGAAGAAAGGGTTAACCAATCCGGATGCGACGCCGGAACAGACGCGCGTCGCGGTGAAAGCGTTCGAAACGCTCAACGGTAACTGGCGCTCACCGGGCGGTGCGGTGAAATTTAACACCGTTACACCg                                                                  >  2:206622/1‑139 (MQ=255)                                                                   tACCTGAAGAAAGGGTTAACCAATCCGGATGCGACGCCGGAACAGACGCGCGTCGCGGTGAAAGCGTTCGAAACGCTCAACGGTAACTGGCGCTCACCGGGCGGTGCGGTGAAATTTAACACCGTTACACCGTCGGTGa                                                           >  2:470651/1‑139 (MQ=255)                                                                     ccTGAAGAAAGGGTTAACCAATCCGGATGCGACGCCGGAACAGACGCGCGTCGCGGTGAAAGCGTTCGAAACGCTCAACGGTAACTGGCGCTCACCGGGCGGTGCGGTGAAATTTAACACCGTTACACCGTCGGTGAc                                                          <  1:207240/138‑1 (MQ=255)                                                                      cTGAAGAAAGGGTTAACCAATCCGGATGCGACGCCGGAACAGACGCGCGTCGCGGTGAAAGCGTTCGAAACGCTCAACGGTAACTGGCGCTCACCGGGCGGTGCGGTGAAATTTaa                                                                               <  1:1044395/116‑1 (MQ=255)                                                                      cTGAAGAAAGGGTTAACCAATCCGGATGCGACGCCGGAACAGACGCGCGTCGCGGTGAAAGCGTTCGAAACGCTCAACGGTAACTGGCGCTCACCGGGCGGTGCGGTGAAATTTaa                                                                               >  2:1044395/1‑116 (MQ=255)                                                                      cTGAAGAAAGGGTTAACCAATCCGGATGCGACGCCGGAACAGACGCGCGTCGCGGTGAAAGCGTTCGAAACGCTCAACGGTAACTGGCGCTCACCGGGCGGTGCGGTGAAATTTAACACCGTTACACCGTCGGTGACCg                                                        >  2:169780/1‑139 (MQ=255)                                                                      cTGAAGAAAGGGTTAACCAATCCGGATGCGACGCCGGAACAGACGCGCGTCGCGGTGAAAGCGATCGAAACGCTCAACGGTAACTGGCGCTCACCGGGCGGTGCGGTGAAATTTaa                                                                               <  2:688827/116‑1 (MQ=255)                                                                      cTGAAGAAAGGGTTAACCAATCCGGATGCGACGCCGGAACAGACGCGCGTCGCGGTGAAAGCGATCGAAACGCTCAACGGTAACTGGCGCTCACCGGGCGGTGCGGTGAAATTTaa                                                                               >  1:688827/1‑116 (MQ=255)                                                                        gaagaaAGGGTTAACCAATCCGGATGCGACGCCGGAACAGACGCGCGTCGCGGTGAAAGCGTTCGAAACGCTCAACGGTAACTGGCGCTCACCGGGCGGTGCGGTGAAATTTAACACCGTTACACCGTCGGTGACCggg                                                      >  1:1179247/1‑139 (MQ=255)                                                                          agaaAGGGTTAACCAATCCGGATGCGACGCCGGAACAGACGCGCGTCGCGGTGAAAGCGTTCGAAACGCTCAACGGTAACTGGCGCTCACCGGGCGGTGCGGTGAAATTTAACACCGTTACACCGTCGGTGACCGGgcg                                                    <  1:588572/139‑1 (MQ=255)                                                                           gaaAGGGTTAACCAATCCGGATGCGACGCCGGAACAGACGCGCGTCGCGGTGAAAGCGTTCGAAACGCTCAACGGTAACTGGCGCTCACCGGGCGGTGCGGTGAAATTTAACACCGTTACACCGTCGGTGACCGGgcgc                                                   >  1:1131898/1‑139 (MQ=255)                                                                              aGGGTTAACCAATCCGGATGCGACGCCGGAACAGACGCGCGTCGCGGTGAAAGCGTTCGAAACGCTCAACGGTAACTGGCGCTCACCGGGCGGTGCGGTGAAATTTAACACCGTTACACCGTCGGTGACCGGGCGCTgg                                                >  1:279651/1‑139 (MQ=255)                                                                              aGGGTTAACCAATCCGGATGCGACGCCGGAACAGACGCGCGTCGCGGTGAAAGCGTTCGAAACGCTCAACGGTAACTGGCGCTCACCGGGCGGTGCGGTGAAATTTAACACCGTTACACCGTCGGTGACCGGGCGCTgg                                                >  2:667399/1‑139 (MQ=255)                                                                                     aCCAATCCGGATGCGACGCCGGAACAGACGCGCGTCGCGGTGAAAGCGTTCGAAACGCTCAACGGTAACTGGCGCTCACCGGGCGGTGCGGTGAAATTTAACACCGTTACACCGTCGGTGACCGGGCGCTGGTTCTCCg                                         >  2:1016838/1‑139 (MQ=255)                                                                                          tCCGGATGCGACGCCGGAACAGACGCGCGTCGCGGTGAAAGCGTTCGAAACGCTCAACGGTAACTGGCGCTCACCGGGCGGTGCGGTGAAATTTAACACCGTTaca                                                                     >  2:1071070/1‑106 (MQ=255)                                                                                          tCCGGATGCGACGCCGGAACAGACGCGCGTCGCGGTGAAAGCGTTCGAAACGCTCAACGGTAACTGGCGCTCACCGGGCGGTGCGGTGAAATTTAACACCGTTaca                                                                     <  1:1071070/106‑1 (MQ=255)                                                                                          tCCGGATGCGACGCCGGAACAGACGCGCGTCGCGGTGAAAGCGTTCGAAACGCTCAACGGTAACTGGCGCTCACCGGGCGGTGCGGTGAAATTTAACACCGTTACACCGTCGGTGACCGGGCGCTGGTTCTCCGGCAAt                                    <  2:1179247/139‑1 (MQ=255)                                                                                          tCCGGATGCGACGCCGGAACAGACGCGCGTCGCGGTGAAAGCGTTCGAAACGCTCAACGGTAACTGGCGCTCACCGGGCGGTGCGGTGAAATTTAACACCGTTACACCGTCGGTGACCGGGCGCTGGTTCTCCGGCAAt                                    <  1:269123/139‑1 (MQ=255)                                                                                          tCCGGATGCGACGCCGGAACAGACGCGCGTCGCGGTGAAAGCGTTCGAAACGCTCAACGGTAACTGGCGCTCACCGGGCGGTGCGGTGAAATTTAACACCGTTACACCGTCGGTGACCGGGCGCTGGTTCTCCGGCAAt                                    >  1:635881/1‑139 (MQ=255)                                                                                          tCCGGATGCGACGCCGGAACAGACGCGCGTCGCGGTGAAAGCGTTCGAAACGCTCAACGGTAACTGGCGCTCACCGGGCGGTGCGGTGAAATTTAACACCGTTACACCGTCGGTGACCGGGCGCTGGTTCTCCGGCAAt                                    >  1:1193020/1‑139 (MQ=255)                                                                                          tCCGGATGCGACGCCGGAACAGACGCGCGTCGCGGTGAAAGCGTTCGAAACGCTCAACGGTAACTGGCGCTCACCGGGCGGTGCGGTGAAATTTAACACCGTTACACCGTCGGTGACCGGGCGCTGGTTCTCCGGCAAt                                    >  2:269123/1‑139 (MQ=255)                                                                                          tCCGGATGCGACGCCGGAACAGACGCGCGTCGCGGTGAAAGCGTTCGAAACGCTCAACGGTAACTGGCGCTCACCGGGCGGTGCGGCGAAATTTAACACCGTTACACCGTCGGTGACCGGGCGCTGGTTCTCCGGCAAt                                    >  1:829518/1‑139 (MQ=255)                                                                                          tCCGGATGCGACGCCGGAACAGACGCGCGTCGCGGTGAAAGCGATCGAAACGCt                                                                                                                         <  2:1044544/54‑1 (MQ=255)                                                                                          tCCGGATGCGACGCCGGAACAGACGCGCGTCGCGGTGAAAGCGATCGAAACGCt                                                                                                                         >  1:1044544/1‑54 (MQ=255)                                                                                               aTGCGACGCCGGAACAGACGCGCGTCGCGGTGAAAGCGTTCGAAACGCTCAACGGTAACTGGCGCTCACCGGGCGGTGCGGTGAAATTTAACACCGTTACACCGTCGGTGACCGGGCGCTGGTTCTCCGGCAATCAGAc                               >  1:315170/1‑139 (MQ=255)                                                                                               aTGCGACGCCGGAACAGACGCGCGTCGCGGTGAAAGCGTTCGAAACGCTCAACGGTAACTGGCGCTCACCGGGCGGTGCGGTGAAATTTAACACCGTTACACCGTCGGTGACCGGGCGCTGGTTCTCCGGCAATCAGAc                               <  1:832230/139‑1 (MQ=255)                                                                                                   gACGCCGGAACAGACGCGCGTCGCGGTGAAAGCGTTCGAAACGCTCAACGGTAACTGGCGCTCACCGGGCGGTGCGGTGAAATTTAACACCGTTACACCGTCGGTGACCGGGCGCTGGTTCTCCGGCAATCAGACCTgg                           >  2:291435/1‑139 (MQ=255)                                                                                                   gACGCCGGAACAGACGCGCGTCGCGGTGAAAGCGTTCGAAACGCTCAACGGTAACTGGCGCTCACCGGGCGGTGCGGTGAAATTTAACACCGTTACACCGTCGGTGACCGGGCGCTGGTTCTCCGGCAATCAGACCTgg                           <  1:1016838/139‑1 (MQ=255)                                                                                                    aCGCCGGAACAGACGCGCGTCGCGGTGAAAGCGTTCGAAACGCTCAACGGTAACTGGCGCTCACCGGGcggtgcggtg                                                                                       >  1:868194/1‑78 (MQ=255)                                                                                                    aCGCCGGAACAGACGCGCGTCGCGGTGAAAGCGTTCGAAACGCTCAACGGTAACTGGCGCTCACCGGGcggtgcggtg                                                                                       <  2:868194/78‑1 (MQ=255)                                                                                                    aCGCCGGAACAGACGCGCGTCGCGGTGAAAGCGTTCGAAACGCTCAACGGTAACTGGCGCTCACCGGGCGGTGCGGTGAAATTTAACACCGTTACACCGTCGGTGACCGGGCGCTGGTTCTCCGGCAATCAGACCTGGc                          >  1:191091/1‑139 (MQ=255)                                                                                                    aCGCCGGAACAGACGCGCGTCGCGGTGAAAGCGTTCGAAACGCTCAACGGTAACTGGCCCTCACCGGGCGGAGCGGAGAAATTTAACACCGTTACACCGTCGGTGACCGGGCGCTGGTTCTCCTGCACTCAGACATGGc                          >  2:317281/1‑139 (MQ=255)                                                                                                          gAACAGACGCGCGTCGCGGTGAAAGCGTTCGAAACGCTCAACGGTAACTGGCGCTCACCGGGCGGTGCGGTGAAATTTAACACCGTTACACCGTCGGTGACCGGGCGCTGGTTCTCCGGCAATCAGACCTGGCCGTggg                    <  2:635881/139‑1 (MQ=255)                                                                                                           aaCAGACGCGCGTCGCGGTGAAAGCGTTCGAAACGCTCAACGGTAACTGGCGCTCACCGGGCGGTGCGGTGAAATTTAACACCGTTACACCGTCGGTGACCGGGCGCTGGTTCTCCGGCAATCAGACCTGGCCGTGGGa                   <  2:1219010/139‑1 (MQ=255)                                                                                                               gACGCGCGTCGCGGTGAAAGCGTTCGAAACGCTCAACGGTAACTGGCGCTCACCGGGCGGTGCGGTGAAATTTAACACCGTTACACCGTCGGTGACCGGGCGCTGGTTCTCCGGCAATCAGACCTGGCCGTGGGATAcc               <  1:161627/139‑1 (MQ=255)                                                                                                                aCGCGCGTCGCGGTGAAAGCGTTCGAAACGCTCAACGGTAACTGGCGCTCACCGGGCGGTGCGGTGAAATTTAACACCGTTACACCGTCGGTGACCGGGCGCTGGTTCTCCGGCAATCAGACCTGGCCGTGGGATACCt              >  2:596386/1‑139 (MQ=255)                                                                                                                     cgTCGCGGTGAAAGCGTTCGAAACGCTCAACGGTAACTGGCGCTCACCGGGCGGTGCGGTGAAATTTAACACCGTTACACCGTCGGTGACCGGGCGCTGGTTCTCCGGCAATCAGACCTGGCCGTGGGATACCTGGAAg         >  1:942180/1‑139 (MQ=255)                                                                                                                     cgTCGCGGTGAAAGCGTTCGAAACGCTCAACGGTAACTGGCGCTCACCGGGCGGTGCGGTGAAATTTAACACCGTTACACCGTCGGTGACCGGGCGCTGGTTCTCCGGCAATCAGACCTGGCCGTGGGATACCTGGAAg         >  1:879401/1‑139 (MQ=255)                                                                                                                          cgGTGAAAGCGTTCGAAACGCTCAACGGTAACTGGCGCTCACCGGGCGGTGCGGTGAAATTTAACACCGTTACACCGTCGGTGACCGGGCGCTGGTTCTCCGGCAATCAGACCTGGCCGTGGGATACCTGGAAGCAGGc    >  1:1037409/1‑139 (MQ=255)                                                                                                                            gTGAAAGCGTTCGAAACGCTCAACGGTAGCTGGCGCTCACCGGGCGGTGCGGTGAAATTTAACACCGTTACACCGTCGGTGACCGGGCGCTGGTTCTCCGGCAATCAGACCTGGCCGTGGGATACCTGGAAGCAGGCGt  <  1:317281/139‑1 (MQ=255)                                                                                                                                     |                                                                                                                                  CAGATCCGCGATATTCTGGCGCGTCCGGCGTTTTATCTCACCGCCTCGCAGCAACGCTGGGAAGAATACCTGAAGAAAGGGTTAACCAATCCGGATGCGACGCCGGAACAGACGCGCGTCGCGGTGAAAGCGATCGAAACGCTCAACGGTAACTGGCGCTCACCGGGCGGTGCGGTGAAATTTAACACCGTTACACCGTCGGTGACCGGGCGCTGGTTCTCCGGCAATCAGACCTGGCCGTGGGATACCTGGAAGCAGGCGT  >  NZ_CP009273/3223042‑3223303

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑