Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,655,800	C→T	7.3%	G649D (GGC→GAC)	tktA ←	transketolase A

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,655,800	0	C	T	7.3%	158.8 / 2.5	55	G649D (GGC→GAC)	tktA	transketolase A
Reads supporting (aligned to +/- strand):  ref base C (31/20);  new base T (2/2);  total (33/22)
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.89e-01

CGGGCAAGCCCGCTCCCACAGGGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTTCTCCAGGGTGAAGCCAAACTCTTCGAACAATGCCGAAGCCGGCGCCGACTCACCGTAGGTGGTCATGCCGATGATACGACCTTCCAGGCCGACGTACTTGTACCAGAAGTCGGCGTGGGCAGCTTCGATGGCGATGCGCGCACCGACTTCCAGCGGCAGCACG  >  NC_002947/5655657‑5655949                                                                                                                                                |                                                                                                                                                      cgGGCAAGCCCGCTCCCACAGGGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTTCTCCAGGGTGAAGCCAAAct                                                                                                                                                 <  2:637995/150‑1 (MQ=255)          ccGCTCCCACAGGGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACGGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTTCTCCAGGGTGAAGCCAAACTCTTCGAACa                                                                                                                                        >  1:863180/1‑150 (MQ=255)          ccGCTCCCACAGGGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTTCTCCAGGGTGAAGCCAAACTCTTCGAACa                                                                                                                                        >  1:48285/1‑150 (MQ=255)          ccGCTCACACAGGGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTTCTCCAGGGTGAAGCCAAACTCTTCGAACa                                                                                                                                        >  1:751357/1‑150 (MQ=255)          cagCTCCCACAGAGAATCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTTCTCCAGGGTGAAGCCAAACTCTTCGAACa                                                                                                                                        >  2:68395/3‑150 (MQ=255)                ccACAGGGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTTCTCCAGGGTGAAGCCAAACTCtt                                                                                                                                              <  1:987436/138‑1 (MQ=255)                ccACAGGGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTTCTCCAGGGTGAAGCCAAACTCtt                                                                                                                                              >  2:987436/1‑138 (MQ=255)                cacagaggTTGCATGAAACC‑CTGTGGGAGCGGGCTTGCCCGCGAAT‑‑‑‑GCCACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTTCTCCAGGGTGAAGCCAAACTCTTCGAACAATGCCGAAGcc                                                                                                                             >  1:817628/7‑150 (MQ=25)                  acaGGGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTTCTCCAGGGTGAAGCCAAACTCTTCGAACAATGCCGaa                                                                                                                                >  1:573511/1‑150 (MQ=255)                              gCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATGTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTTCTCCAGGGTGAAGCCAAACTCTTCGAACAATGCCGAAGCCGGCGCCGAc                                                                                                                    >  2:597310/1‑150 (MQ=255)                              gCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTTCTCCAGGGTGAAGCCAAACTCTTCGAACAATGCCGAAGCCGGCGCCGAc                                                                                                                    >  2:751448/1‑150 (MQ=255)                                  ccTCTGTGGGAGCGGGCTTGCCCTCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTTCTCCAGGGTGAAGTCAAACTCTTCGAACAATGCCGAAGCCGGCGCCGACTCAc                                                                                                                >  2:551557/1‑150 (MQ=255)                                  ccTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTTCTCCAGGGTGAAGTCAAACTCTTCGAACAATGCCGAAGCCGGCGCCGACTCAc                                                                                                                >  2:551656/1‑150 (MQ=255)                                    tctGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTTCTCCAGGGTGAAGCCAAACTCTTCGAACAATGCCGAAGCCGGCGCCGACTCACCg                                                                                                              >  1:784167/1‑150 (MQ=255)                                    cctGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTTCTCCAGGGTGAAGCCAAACTCTTCGAACAATGCCGAAGCCGGCGCCGACTCACCg                                                                                                              >  1:443530/2‑150 (MQ=255)                                         gggAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTTCTCCAGGGTGAAGCCAAACTCTTCGAACAATGCCGAAGCCGGCGCCGACTCACCGTAggt                                                                                                         <  2:751357/150‑1 (MQ=255)                                                ggCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTTCTCCAGGGTGAAGCCAAACTCTTCGAACAATGCCGAAGCCGGCGCCGACTCACCGTAGGTGGTCATg                                                                                                  >  1:1010255/1‑150 (MQ=255)                                                       ccGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTTCTCCAGGGTGAAGCCAAACTCTTCGAACAATGCCGAAGCCGGCGCCGACTCACCGTAggt                                                                                                         <  2:906750/136‑1 (MQ=255)                                                       ccGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTTCTCCAGGGTGAAGCCAAACTCTTCGAACAATGCCGAAGCCGGCGCCGACTCACCGTAggt                                                                                                         >  1:906750/1‑136 (MQ=255)                                                         gcgAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTTCTCCAGGGTGAAGCCAAACTCTTCGAACAATGCCGAAGCCGGCGCCGACTCACCGTAGGTGGTCATGCCGATGATa                                                                                         <  2:573511/150‑1 (MQ=255)                                                         gcgAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTTCTCCAGGGTGAAGCCAAACTCTTCGAACAATACCGAAGCCGGCGCCGACTCACCGTAGGTGGTCATGCCGATGATa                                                                                         >  2:805731/1‑150 (MQ=255)                                                            aaTGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCTGCAGTACCGAGGATGTTCTCCAGGGTGAAGCCAAACTCTTCGAACAATGCCGAAGCCGGCGCCGACTCACCGTAGGTGGTCATGCCGATGATACGa                                                                                      <  1:68395/150‑1 (MQ=255)                                                                 aGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTTCTCCAGGGTGAAGCCAAACTCTTCGAACAATGCCGAAGCCGGCGCCGACTCACCGTAGGTGGTCATGCCGATGATACGACCTTc                                                                                 >  1:384330/1‑150 (MQ=255)                                                                   gCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTTCTCCAGGGTGAAGCCAAACTCTTCGAACAATGCCGAAGCCGGCGCCGACTCACCGTAGGTGGTCATGCCGATGATACGACCTTCCa                                                                               <  1:751448/150‑1 (MQ=255)                                                                        acCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTTCTCCAGGGTGAAGCCAAACTCTTCGAACAATGCCGAAGCCGGCGCCGACTCACCGTAGGTGGTCATGCCGATGATACGACCTTCCAGGCCg                                                                          >  1:894194/1‑150 (MQ=255)                                                                        acCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTTCTCCAGGGTGAAGCCAAACTCTTCGAACAATGCCGAAGCCGGCGCCGACTCACCGTAGGTGGTCATGCCGATGATACGACCTTCCAGGCCg                                                                          >  1:184600/1‑150 (MQ=255)                                                                           gcccgcAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTTCTCCAGGGTGAAGCCAAACTCTTCGAACAATGCCGAAGCCGGCGCCGACTCACCGTAg                                                                                                           >  2:48430/1‑114 (MQ=255)                                                                           gcccgcAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTTCTCCAGGGTGAAGCCAAACTCTTCGAACAATGCCGAAGCCGGCGCCGACTCACCGTAg                                                                                                           <  1:48430/114‑1 (MQ=255)                                                                                cAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTTCTCCAGGGTGAAGCCAAACTCTTCGAACAATGCCGAAGCCGGCGCCGACTCACCGTAGGTGGTCATGCCGATGATACGACCTTCCAGGCCGACGTACtt                                                                  <  2:184600/150‑1 (MQ=255)                                                                                 aGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTTCTCCAGGGTGAAGCCAAACTCTTCGAACAATGCCGAAGCCGGCGCCGACTCACCGTAGGTGGTCATGCCGATGATACGACCTTCCAGGCCGACGTACTTg                                                                 >  2:364765/1‑150 (MQ=255)                                                                                    tCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTTCTCCAGGGTGAAGCCAAACTCTTCGAACAATGCCGAAGCCGGCGCCGACTCACCGTAGGTGGTCATGCCGATGATACGACCTTCCAGGCCGACGTACTTGTAc                                                              <  2:48285/150‑1 (MQ=255)                                                                                            tCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTTCTCCAGGGTGAAGCCAAACTCTTCGAACAATGCCGAAGCCGGCGCCGACTCACCGTAGGTGGTCATGCCGATGATACGACCTTCCAGGCCGACGTACTTGTACCAGAAGTc                                                      >  2:860180/1‑150 (MQ=255)                                                                                             cAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTTCTCCAGGGTGAAGCCAAACTCTTCGAACAATGCCGAAGCCGGCGCCGACTCACCGTAGGTGGTCATGCCGATGATACGACCTTCCAGGCCGACGTACTTGTACCAGAAGTCg                                                     <  1:239683/150‑1 (MQ=255)                                                                                                tCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTTCTCCAGGGTGAAGCCAAACTCTTCGAACAATGCCGAAGCCGGCGCCGACTCACCGTAGGTGGTCATGCCGCTGATACGACCTTCCGGGCCGACGTACTGGTACGAGAAGTCGGCg                                                  >  2:796964/1‑150 (MQ=255)                                                                                                 cTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTTCTCCAGGGTGAAGCCAAACTCTTCGAACAATGCCGAAGCCGGCGCCGACTCACCGTAGGTGGTCATGCCGATGATACGACCTTCCAGGCCGACGTACTTGTACCAGAAGTCGGCGt                                                 <  1:796964/150‑1 (MQ=255)                                                                                                  ttCCAGCAGCTCTTCGGCAGTACCGAGGATGTTCTCCAGGGTGAAGCCAAACTCTTCGAACAATGCCGAAGCCGGCGCCGACTCACCGTAGGTGGTCATGCCGATGATACGACCTTCCAGGCCGACGTACTTGTACCAGAAGTCGGCGTg                                                >  1:768016/1‑150 (MQ=255)                                                                                                  ttCCAGCAGCTCTTCGGCAGTACCGAGGATGTTCTCCAGGGTGAAGCCAAACTCTTCGAACAATGCCGAAGCCGGCGCCGACTCACCGTAGGTGGTCATGCCGATGATACGACCTTCCAGGCCGACGTACTTGTAACAGAAGTCGGCGTg                                                >  2:64728/1‑150 (MQ=255)                                                                                                   tCCAGCAGCTCTTCGGCAGTACCGAGGATGTTCTCCAGGGTGAAGTCAAACTCTTCGAACAATGCCGAAGCCGGCGCCGACTCACCGTAGGTGGTCATGCCGATGATACGACCTTCCAGGCCGACGTACTTGTACCAGAAGTCGGCGTgg                                               <  1:551557/150‑1 (MQ=255)                                                                                                   tCCAGCAGCTCTTCGGCAGTACCGAGGATGTTCTCCAGGGTGAAGTCAAACTCTTCGAACAATGCCGAAGCCGGCGCCGACTCACCGTAGGTGGTCATGCCGATGATACGACCTTCCAGGCCGACGTACTTGTACCAGAAGTCGGCGTgg                                               <  1:551656/150‑1 (MQ=255)                                                                                                   tCCAGCAGCTCTTCGGCAGTACCGAGGATGTTCTCCAGGGTGAAGCCAAACTCTTCGAACAATGCCGAAGCCGGCGCCGACTCACCGTAGGTGGTCATGCCGATGATACGACCTTCCAGGCCGACGTACTTGTACCAGAAGTCGGCGTgg                                               <  2:443530/150‑1 (MQ=255)                                                                                                       gcagcTCTTCGGCAGTACCGAGGATGTTCTCCAGGGTGAAGCCAAACTCTTCGAACAATGCCGAAGCCGGCGCCGACTCACCGTAGGTGGTCATGCCGATGATACGACCTTCCAGGCCGAc                                                                        >  1:1069166/1‑121 (MQ=255)                                                                                                       gcagcTCTTCGGCAGTACCGAGGATGTTCTCCAGGGTGAAGCCAAACTCTTCGAACAATGCCGAAGCCGGCGCCGACTCACCGTAGGTGGTCATGCCGATGATACGACCTTCCAGGCCGAc                                                                        <  2:1069166/121‑1 (MQ=255)                                                                                                         agcTCTTCGGCAGTACCGAGGATGTTCTCCAGGGTGAAGCCAAACTCTTCGAACAATGCCGAAGCCGGCGCCGACTCACCGTAGGTGGTCATGCCGATGATACGACCTTCCAGGCCGACGTACTTGTACCAGAAGTCGGCGTGGGCAGCt                                         <  1:597310/150‑1 (MQ=255)                                                                                                              ttCGGCAGTACCGAGGATGTTCTCCAGGGTGAAGCCAAACTCTTCGAACAATGCCGAAGCCGGCGCCGACTCACCGTAGGTGGTCATGCCGATGATACGACCTTCCAGGCCGACGTACTTGTACCAGAAGTCGGCGTGGGCAGCTTCGAt                                    >  1:266950/1‑150 (MQ=255)                                                                                                                cGGCAGTACCGAGGATGTTCTCCAGGGTGAAGCCAAACTCTTCGAACAATGCCGAAGCCGGCGCCGACTCACCGTAGGTGGTCATGCCGATGATACGACCTTCCAGGCCGACGTACTTGTACCAGAAGTCGGCGTGGGCAGCTTCGATgg                                  <  2:1010255/150‑1 (MQ=255)                                                                                                                        ccGAGGATGTTCTCCAGGGTGAAGCCAAACTCTTCGAACAATGCCGAAGCCGGCGCCGACTCACCGTAGGTGGTCATGCCGATGATACGACCTTCCAGGCCGACGTACTTGTACCAGAAGTCGGCGTGGGCAGCTTCGATGGCGATgcgc                          >  1:1080464/1‑150 (MQ=255)                                                                                                                          gAGGATGTTCTCCAGGGTGAAGCCAAACTCTTCGAACAATGCCGAAGCCGGCGCCGACTCACCGTAGGTGGTCATGCCGATGATACGACCTTCCAGGCCGACGTACTTGTACCAGAAGTCGGCGTGGGCAGCTTCGATGGCGATgcgcgc                        >  2:39244/1‑150 (MQ=255)                                                                                                                             gATGTTCTCCAGGGTGAAGCCAAACTCTTCGAACAATGCCGAAGCCGGCGCCGACTCACCGTAGGTGGTCATGCCGATGATACGACCTTCCAGGCCGACGTACTTGTACCAGAAGTCGGCGTGGGCAGCTTCGATGGCGATGCGCGCAcc                     <  1:364765/150‑1 (MQ=255)                                                                                                                                     ccAGGGTGAAGCCAAACTCTTCGAACAATGCCGAAGCCGGCGCCGACTCACCGTAGGTGGTCATGCCGATGATACGACCTTCCAGGCCGACGTACTTGTACCAGAAGTCGGCGTGGGCAGCTTCGATGGCGATGCGCGCACCGACTTCCa             >  1:70140/1‑150 (MQ=255)                                                                                                                                       aGGGTGAAGCCAAACTCTTCGAACAATGCCGAAGCCGGCGCCGACTCACCGTAGGTGGTCATGCCGATGATACGACCTTCCAGGCCGACGTACTTGTAACAGAAGTCGGCGTGGGCAGCTTCGATGGCGATGCGCGCACCGACTTCCAGc           <  1:64728/150‑1 (MQ=255)                                                                                                                                           tGAAGCCAAACTCTTCGAACAATGCCGAAGCCGGCGCCGACTCACCGTAGGTGGTCATGCCGATGATACGACCTTCCAGGCCGACGTACTTGTACCAGAAGTCGGCGTGGGCAGCTTCGATGGCGATGCGCGCACCGACTTCCAGCGgca       >  2:804050/1‑150 (MQ=255)                                                                                                                                            gAAGCCAAACTCTTCGAACAATGCCGAAGCCGGCGCCGACTCACCGTAGGTGGTCATGCCGATGATACGACCTTCCAGGCCGACGTACTTGTACCAGAAGTCGGCGTGGGCAGCTTCGATGGCGATGCGCGCACCGACTTCCAGCGgcag      <  1:860180/150‑1 (MQ=255)                                                                                                                                            gAAGCCAAACTCTTCGAACAATGCCGAAGCCGGCGCCGACTCACCGTAGGTGGTCATGCCGATGATACGACCTTCCAGGCCGACGTACTTGTACCAGAAGTCGGCGTGGGCAGCTTCGATGGCGATGCGCGCACCGACTTCCAGCGgcag      <  2:784167/150‑1 (MQ=255)                                                                                                                                               gCCAAACTCTTCGAACAATGCCGAAGCCGGCGCCGACTCACCGTAGGTGGTCATGCCGATGATACGACCTTCCAGGCCGACGTACTTGTACCAGAAGTCGGCGTGGGCAGCTTCGATGGCGATGCGCGCACCGACTTCCAGCGGCAGCAc   >  1:993179/1‑150 (MQ=255)                                                                                                                                               gCCAAACTCTTCGAACAATGCCGAAGCCGGCGCCGACTCACCGTAGGTGGTCATGCCGATGATACGACCTTCCAGGCCGACGTACTTGTACCAGAAGTCGGCGTGGGCAGCTTCGATGGCGATGCGCGCACCGACTTCCAGCGGCAGCAc   >  2:960458/1‑150 (MQ=255)                                                                                                                                                ccAAACTCTTCGAACAATGCCGAAGCCGGCGCCGACTCACCGTAGGTGGTCATGCCGATGATACGACCTTCCAGGCCGACGTACTTGTACCAGAAGTCGGCGTGGGCAGCTTCGATGGCGATGCGCGCACCGACTTCCAGCGGCAGCACg  >  2:283632/1‑150 (MQ=255)                                                                                                                                                |                                                                                                                                                      CGGGCAAGCCCGCTCCCACAGGGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTTCTCCAGGGTGAAGCCAAACTCTTCGAACAATGCCGAAGCCGGCGCCGACTCACCGTAGGTGGTCATGCCGATGATACGACCTTCCAGGCCGACGTACTTGTACCAGAAGTCGGCGTGGGCAGCTTCGATGGCGATGCGCGCACCGACTTCCAGCGGCAGCACG  >  NC_002947/5655657‑5655949

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