Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,282,649	A→C	K225T (AAA→ACA)	agaS →	tagatose‑6‑phosphate ketose/aldose isomerase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,282,649	0	A	C	100.0%	191.5 / NA	52	K225T (AAA→ACA)	agaS	tagatose‑6‑phosphate ketose/aldose isomerase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base C (22/30);  total (22/30)

TACCACCATGATGGCCAGCTGCCTCGCGGTTTTCGCACCTGAGACGATCAACAGCCAAACCTTCCGCGACGTGGCGGATCGTTGCCAGGCGATCCTGACCTCACTGGGCGATTTCAGCGAAGGTGTGTTTGGTTACGCACCGTGGAAACGGATCGTTTATCTCGGTAGCGGTGGCTTACAGGGCGCAGCACGCGAGTCGGCGCTGAAAGTGCTGGAACTGACGGC‑GGGTAAACTGGCGGCCTTTTATGATTCTCCAACCGGATTCCGTCATGGACCAAAATCGCTGGTCG  >  NC_000913/3282503‑3282792                                                                                                                                                   |                                                                                                                                                 tACCACCATGATGGCCAGCTGCCTCGCGGTTTTCGCACCTGAGACGATCAACAGCCAAACCTTCCGCGACGTGGCGGATCGTTGCCAGGCGATCCTGACCTCACTGGGCGATTTCAGCGAAGGTGTGTTTGGTTACGCACCGTGGACACg                                                                                                                                               <  2:432008/150‑1 (MQ=255)        atgatgGCCAGCTGCCTCGCGGTTTTCGCACCTGAGTCGATCAACAGCCAAACCTTCCGCGACGTGGCGGATCGTTGCCAGGCGATCCTGACCTCACTGGGCGATTTCAGCGAAGGTGTGTTTGGTTACGCACCGTGGACACGGATCGttt                                                                                                                                       <  2:403572/151‑1 (MQ=255)                      ccTCGCGGTTTTCGCACCTGAGACGATCAACAGCCAAACCTTCCGCGACGTGGCGGATCGTTGCCAGGCGATCCTGACCTCACTGGGCGATTTCAGCGAAGGTGTGTTTGGTTACGCACCGTGGACACGGATCGTTTATCTCGGTAGCgg                                                                                                                          >  1:744408/1‑150 (MQ=255)                         cgcgGTTTTCGCACCTGAGACGATCAACAGCCAAACCTTCCGCGACGTGGCGGATCGTTGCCAGGCGATCCTGACCTCACTGGGCGATTTCAGCGAAGGTGTGTTTGGTTACGCACCGTGGACACGGATCGTTTATCTCGGTAGCGGTGGc                                                                                                                      <  2:191522/151‑1 (MQ=255)                            ggTTTTCGCACCTGAGACGATCAACAGCCAAACCTTCCGCGACGTGGCGGATCGTTGCCAGGCGATCCTGACCTCACTGGGCGATTTCAGCGAAGGTGTGTTTGGTTACGCACCGTGGACACGGATCGTTTATCTCGGTAGCGGTGGCtt                                                                                                                    >  2:492208/1‑150 (MQ=255)                                ttCGCACCTGAGACGATCAACAGCCAAACCTTCCGCGACGTGGCGGATCGTTGCCAGGCGATCCTGACCTCACTGGGCGATTTCAGCGAAGGTGTGTTTGGTTACGCACCGTGGACACGGATCGTTTATCTCGGTAGCGGTGGCTTACAg                                                                                                                <  2:500834/150‑1 (MQ=255)                                  cGCACCTGAGACGATCAACAGCCAAACCTTCCGCGACGTGGCGGATCGTTGCCAGGCGATCCTGACCTCACTGGGCGATTTCAGCGAAGGTGTGTTTGGTTACGCACCGTGGACACGGATCGTTTATCTCGGTAGCGGTGGCTTACAggg                                                                                                              <  1:601894/150‑1 (MQ=255)                                         gagaCGATCAACAGCCAAACCTTCCGCGACGTGGCGGATCGTTGCCAGGCGATCCTGACCTCACTGGGCGATTTCAGCGAAGGTGTGTTTGGTTACGCACCGTGGACACGGATCGTTTATCTCGGTAGCGGTGGCTTACAGGGCGCAGCAc                                                                                                      <  2:727503/151‑1 (MQ=255)                                         gagaCGATCAACAGCCAAACCTTCCGCGACGTGGCGGATCGTTGCCAGGCGATCCTGACCTCACTGGGCGATTTCAGCGAAGGTGTGTTTGGTTACGCACCGTGGACACGGATCGTTTATCTCGGTAGCGGTGGCTTACAGGGCGCAGCAc                                                                                                      >  2:293091/1‑151 (MQ=255)                                                 cAACAGCCAAACCTTCCGCGACGTGGCGGATCGTTGCCAGGCGATCCTGACCTCACTGGGCGATTTCAGCGAAGGTGTGTTTGGTTACGCACCGTGGACACGGATCGTTTATCTCGGTAGCGGTGGCTTACAGGGCGCAGCACGCGAGTc                                                                                               <  1:844643/150‑1 (MQ=255)                                                 cAACAGCCAAACCTTCCGCGACGTGGCGGATCGTTGCCAGGCGATCCTGACCTCACTGGGCGATTTCAGCGAAGGTGTGTTTGGTTACGCACCGTGGACACGGATCGTTTATCTCGGTAGCGGTGGCTTACAGGGCGCAGCACGCGAGTc                                                                                               <  1:676259/150‑1 (MQ=255)                                                  aaCAGCCAAACCTTCCGCGACGTGGCGGATCGTTGCCAGGCGATCCTGACCTCACTGGGCGATTTCAGCGAAGGTGTGTTTGGTTACGCACCGTGGACACGGATCGTTTATCTCGGTAGCGGTGGCTTACAGGGCGCAGCACGCGAGTCgg                                                                                             <  2:570979/151‑1 (MQ=255)                                                   aCAGCCAAACCTTCCGCGACGTGGCGGATCGTTGCCAGGCGATCCTGACCTCACTGGGCGATTTCAGCGAAGGTGTGTTTGGTTACGCACCGTGGACACGGATCGTTTATCTCGGTAGCGGTGGCTTACAGGGCGCAGCACGCGAGTCGgc                                                                                            <  2:140912/151‑1 (MQ=255)                                                    cAGCCAAACCTTCCGCGACGTGGCGGATCGTTGCCAGGCGATCCTGACCTCACTGGGCGATTTCAGCGAAGGTGTGTTTGGTTACGCACCGTGGACACGGATCGTTTATCTCGGTAGCGGTGGCTTACAGGGCGCAGCACGCGAGTCGgcg                                                                                           >  1:142533/1‑151 (MQ=255)                                                          aaCCTTCCGCGACGTGGCGGATCGTTGCCAGGCGATCCTGACCTCACTGGGCGATTTCAGCGAAGGTGTGTTTGGTTACGCACCGTGGACACGGATCGTTTATCTCGGTAGCGGTGGCTTACAGGGCGCAGCACGCGAGTCGGCGCTGaa                                                                                      >  1:82569/1‑150 (MQ=255)                                                                 cgcgACGTGGCGGATCGTTGCCAGGCGATCCTGACCTCACTGGGCGATTTCAGCGAAGGTGTGTTTGGTTACGCACCGTGGACACGGATCGTTTATCTCGGTAGCGGTGGCTTACAGGGCGCAGCACGCGAGTCGGCGCTGAAAGTGCTgg                                                                              >  1:501702/1‑151 (MQ=255)                                                                 cgcgACGTGGCGGATCGTTGCCAGGCGATCCTGACCTCACTGGGCGATTTCAGCGAAGGTGTGTTTGGTTACGCACCGTGGACACGGATCGTTTATCTCGGTAGCGGTGGCTTACAGGGCGCAGCACGCGAGTCGGCGCTGAAAGTGCTgg                                                                              >  2:175939/1‑151 (MQ=255)                                                                       gTGGCGGATCGTTGCCAGGCGATCCTGACCTCACTGGGCGATTTCAGCGAAGGTGTGTTTGGTTACGCACCGTGGACACGGATCGTTTATCTCGGTAGCGGTGGCTTACAGGGCGCAGCACGCGAGTCGGCGCTGAAAGTGCTGGAACTg                                                                         >  1:131682/1‑150 (MQ=255)                                                                       gTGGCGGATCGTTGCCAGGCGATCCTGACCTCACTGGGCGATTTCAGCGAAGGTGTGTTTGGTTACGCACCGTGGACACGGATCGTTTATCTCGGTAGCGGTGGCTTACAGGGCGCAGCACGCGAGTCGGCGCTGAAAGTGCTGGAACTg                                                                         <  1:759683/150‑1 (MQ=255)                                                                        tGGCGGATCGTTGCCAGGCGATCCTGACCTCACTGGGCGATTTCAGCGAAGGTGTGTTTGGTTACGCACCGTGGACACGGATCGTTTATCTCGGTAGCGGTGGCTTACAGGGCGCAGCACGCGAGTCGGCGCTGAAAGTGCTGGAACTGAc                                                                       >  1:755965/1‑151 (MQ=255)                                                                           cGGATCGTTGCCAGGCGATCCTGACCTCACTGGGCGATTTCAGCGAAGGTGTGTTTGGTTACGCACCGTGGACACGGATCGTTTATCTCGGTAGCGGTGGCTTACAGGGCGCAGCACGCGAGTCGGCGCTGAAAGTGCTGGAACTGAcggc                                                                    >  2:880562/1‑151 (MQ=255)                                                                              aTCGTTGCCAGGCGATCCTGACCTCACTGGGCGATTTCAGCGAAGGTGTGTTTGGTTACGCACCGTGGACACGGATCGTTTATCTCGGTAGCGGTGGCTTACAGGGCGCAGCACGCGAGTCGGCGCTGAAAGTGCTGGAACTGACGGC‑ggg                                                                <  2:727810/151‑1 (MQ=255)                                                                                 gTTGCCAGGCGATCCTGACCTCACTGGGCGATTTCAGCGAAGGTGTGTTTGGTTACGCACCGTGGACACGGATCGTTTATCTCGGTAGCGGTGGCTTACAGGGCGCAGCACGCGAGTCGGCGCTGAAAGTGCTGGAACTGACGGC‑GGGTaa                                                             <  2:176583/151‑1 (MQ=255)                                                                                 gTTGCCAGGCGATCCTGACCTCACTGGGCGATTTCAGCGAAGGTGTGTTTGGTTACGCACCGTGGACACGGATCGTTTATCTCGGTAGCGGTGGCTTACAGGGCGCAGCACGCGAGTCGGCGCTGAAAGTGCTGGAACTGACGGC‑GGGTa                                                              >  1:592844/1‑150 (MQ=255)                                                                                    gCCAGGCGGTCCTGACCTCACTGGGCGATTTCAGCGAAGGTGTGTTTGGTTACGCACCGTGGACACGGATCGTTTATCTCGGTAGCGGTGGCTTACAGGGCGCAGCACGCGAGTCGGCGCTGAAAGTGCTGGAACTGACGGC‑GGGTAAAc                                                           <  2:256292/150‑1 (MQ=255)                                                                                     ccAGGCGATCCTGACCTCACTGGGCGATTTCAGCGAAGGTGTGTTTGGTTACGCACCGTGGACACGGATCGTTTATCTCGGTAGCGGTGGCTTACAGGGCGCAGCACGCGAGTCGGCGCTGAAAGTGCTGGAACTGACGGC‑GGGTAAACTg                                                         <  2:489834/151‑1 (MQ=255)                                                                                                 gACCTCACTGGGCGATTTCAGCGAAGGTGTGTTTGGTTACGCACCGTGGACACGGATCGTTTATCTCGGTAGCGGTGGCTTACAGGGCGCAGCACGCGAGTCGGCGCTGAAAGTGCTGGAACTGACGGC‑GGGTAAACTGGCGGCCtttt                                               <  1:91469/149‑1 (MQ=255)                                                                                                 gACCTCACTGGGCGATTTCAGCGAAGGTGTGTTTGGTTACGCACCGTGGACACGGATCGTTTATCTCGGTAGCGGTGGCTTACAGGGCGCAGCACGCGAGTCGGCGCTGAAAGTGCTGGAACTGACGGC‑GGGTAAACTGGCGGCCtttt                                               <  1:43709/149‑1 (MQ=255)                                                                                                 gACCTCACTGGGCGATTTCAGCGAAGGTGTGTTTGGTTACGCACCGTGGACACGGATCGTTTATCTCGGTAGCGGTGGCTTACAGGGCGCAGCACGCGAGTCGGCGCTGAAAGTGCTGGAACTGACGGC‑GGGTAAACTGGCGGCCTTTTAt                                             <  2:595489/151‑1 (MQ=255)                                                                                                    cTCACTGGGCGATTTCAGCGAAGGTGTGTTTGGTTACGCACCGTGGACACGGATCGTTTATCTCGGTAGCGGTGGCTTACAGGGCGCAGCACGCGAGTCGGCGCTGAAAGTGCTGGAACTGACGGC‑GGGTAAACTGGCGGCCTTTTATGAt                                          >  1:140451/1‑151 (MQ=255)                                                                                                              gATTTCAGCGAAGGTGTGTTTGGTTACGCACCGTGGACACGGATCGTTTATCTCGGTAGCGGTGGCTTACAGGGCGCAGCACGCGAGTCGGCGCTGAAAGTGCTGGAACTGACGGC‑GGGTAAACTGGCGGCCTTTTATGATTCTCCAACCg                                <  1:280922/151‑1 (MQ=255)                                                                                                                 ttCAGCGAAGGTGTGTTTGGTTACGCACCGTGGACACGGATCGTTTATCTCGGTAGCGGTGGCTTACAGGGCGCAGCACGCGAGTCGGCGCTGAAAGTGCTGGAACTGACGGC‑GGGTAAACTGGCGGCCTTTTATGATTCTCCAACCGGAt                             >  1:865651/1‑151 (MQ=255)                                                                                                                 ttCAGCGAAGGTGTGTTTGGTTACGCACCGTGGACACGGATCGTTTATCTCGGTAGCGGTGGCTTACAGGGCGCAGCACGCGAGTCGGCGCTGAAAGTGCTGGAACTGACGGC‑GGGGAAACTGGCGGCCTTTTATGATTCTCCAACCGGAt                             >  1:696556/1‑151 (MQ=255)                                                                                                                  tCAGCGAAGGTGTGTTTGGTTACGCACCGTGGACACGGATCGTTTATCTCGGTAGCGGTGGCTTACAGGGCGCAGCACGCGAGTCGGCGCTGAAAGTGCTGGAACTGACGGC‑GGGTAAACTGGCGGCCTTTTTTGATTCTCCAACCGGAtt                            >  2:318501/1‑151 (MQ=255)                                                                                                                    aGCGAAGGTGTGTTTGGTTACGCACCGTGGACACGGATCGTTTATCTCGGTAGCGGTGGCTTACAGGGCGCAGCACGCGAGTCGGCGCTGAAAGTGCTGGAACTGACGGC‑GGGTAAACTGGCGGCCTTTTATGATTCTCCAACCGGATTcc                          <  2:554385/151‑1 (MQ=255)                                                                                                                      cGAAGGTGTGTTTGGTTACGCACCGTGGACACGGATCGTTTATCTCGGTAGCGGTGGCTTACAGGGCGCAGCACGCGAGTCGGCGCTGAAAGTGCTGGAACTGACGGC‑GGGTAAACTGGCGGCCTTTTATGATTCTCCAACCGGATTcc                          >  1:416106/1‑149 (MQ=255)                                                                                                                      cGAAGGTGTGTTTGGTTACGCACCGTGGACACGGATCGTTTATCTCGGTAGCGGTGGCTTACAGGGCGCAGCACGCGAGCCGGCGCTGAAAGTGCTGGAACTGACGGC‑GGGGATCAACGAGGGCTTTTATGATTTTCCaaaccggt                             >  2:151132/1‑140 (MQ=255)                                                                                                                       gAAGGTGTGTTTGGTTACGCACCGTGGACACGGATCGTTTATCTCGGTAGCGGTGGCTTACAGGGCGCAGCACGCGAGTCGGCGCTGAAAGTGCTGGAACTGACGGC‑GGGTAAACTGGCGGCCTTTTATGATTCTCCAACCGGATTCCGTc                       <  1:719730/151‑1 (MQ=255)                                                                                                                       gAAGGTGTGTTTGGTTACGCACCGTGGACACGGATCGTTTATCTCGGTAGCGGTGGCTTACAGGGCGCAGCACGCGAGTCGGCGCTGAAAGTGCTGGAACTGACGGC‑GGGGAAACTGGCGGCCTTTTATGATTCTCCAACCGGATTcc                          >  1:522918/1‑148 (MQ=255)                                                                                                                          ggTGTGTTTGGTTACGCACCGTGGACACGGATCGTTTATCTCGGTAGCGGTGGCTTACAGGGCGCAGCACGCGAGTCGGCGCTGAAAGTGCTGGAACTGACGGC‑GGGTAAACTGGCGGCCTTTTATGATTCTCCAACCGGATTCCGTCATg                    <  2:307050/151‑1 (MQ=255)                                                                                                                           gtgtgtTTGGTTACGCACCGTGGACACGGATCGTTTATCTCGGTAGCGGTGGCTTACAGGGCGCAGCACGCGAGTCGGCGCTGAAAGTGCTGGAACTGACGGC‑GGGTAAACTGGGGGGCTTTTTTGATTTTCCAAACGGATTCCGTCATgg                   >  2:551854/1‑151 (MQ=255)                                                                                                                           gtgtgtTTGGTTACGCACCGTGGACACGGATCGTTTATCTCGGTAGCGGTGGCTTACAGGGCGCAGCACGCGAGTCGGCGCTGAAAGTGCTGGAACTGACGGC‑GGGTAAACTGGCGGCCTTTTATGATTCTCCAACCGGATTCCGTCATgg                   <  2:281389/151‑1 (MQ=255)                                                                                                                                 ttGGTTACGCACCGTGGACACGGATCGTTTATCTCGGTAGCGGTGGCTTACAGGGCGCAGCACGCGAGTCGGCGCTGAAAGTGCTGGAACTGACGGC‑GGGTAAACTGGCGGCCTTTTATGATTCTCCAACCGGATTCCGTCATGGACCaa              <  1:41016/150‑1 (MQ=255)                                                                                                                                   ggTTACGCACCGTGGACACGGATCGTTTATCTCGGTAGCGGTGGCTTACAGGGCGCAGCACGCGAGTCGGCGCTGGAAGTGCTGGAACTGACGGC‑GGGTAAACAGGCGGCCTTTTATGATTCTCCAACCGGATTCCGTCATGGACCAAAAt           <  2:532744/151‑1 (MQ=255)                                                                                                                                      tACGCACCGTGGACACGGATCGTTTATCTCGGTAGCGGTGGCTTACAGGGCGCAGCACGCGAGTCGGCGCTGAAAGTGCTGGAACTGACGGC‑GGGTAAACTGGCGGCCTTTTATGATTCTCCAACCGGATTCCGTCATGGACCAAAATCg         <  2:198577/150‑1 (MQ=255)                                                                                                                                      tACGCACCGTGGACACGGATCGTTTATCTCGGTAGCGGTGGCTTACAGGGCGCAGCACGCGAGTCGGCGCTGAAAGTGCTGGAACTGACGGC‑GGGGAAACTGGCGGCCTTTTATGATTCTCCAACCGGATTCCGTCATGGACCaaaa            >  2:578252/1‑147 (MQ=255)                                                                                                                                           aCCGTGGACACGGATCGTTTATCTCGGTAGCGGTGGCTTACAGGGCGCAGCACGCGAGTCGGCGCTGAAAGTGCTGGAACTGACGGC‑GGGTAAACTGGCGGCCTTTTATGATTCTCCAACCGGATTCCGTCATGGACCAAAATCGCTGGTc   <  2:521359/151‑1 (MQ=255)                                                                                                                                           aCCGTGGACACGGATCGTTTATCTCGGTAGCGGTGGCTTACAGGGCGCAGCACGCGAGTCGGCGCTGAAAGTGCTGGAACTGACGGC‑GGGTAAACTGGCGGCCTTTTATGATTCTCCAACCGGATTCCGTCATGGACCAAAATCGCTGGTc   <  2:501689/151‑1 (MQ=255)                                                                                                                                           aCCGTGGACACGGATCGTTTATCTCGGTAGCGGTGGCTTACAGGGCGCAGCACGCGAGTCGGCGCTGAAAGTGCTGGAACTGACGGC‑GGGTAAACTGGCGGCCTTTTATGATTCTCCAACCGGATTCCGTCATGGACCAAAATCGCTGGTc   <  1:371053/151‑1 (MQ=255)                                                                                                                                            ccGTGGACACGGATCGTTTATCTCGGTAGCGGTGGCTTACAGGGCGCAGCACGCGAGTCGGCGCTGAAAGTGCTGGAACTGACGGC‑GGGTAAACTGGCGGCCTTTTATGATTCTCCAACCGGATTCCGTCATGGACCAAAATCGCTGGTCg  <  1:51317/151‑1 (MQ=255)                                                                                                                                             cGTGGACACGGATCGTTTATCTCGGTAGCGGTGGCTTACAGGGCGCAGCACGCGAGTCGGCGCTGAAAGTGCTGGAACTGACGGC‑GGGGAAACTGGCGGCCTTTTATGATTCTCCAACCGGGTTCCGGCATGGACCAAAAACGCTGGTCg  >  2:119687/1‑150 (MQ=255)                                                                                                                                              gTGGACACGGATCGTTTATCTCGGTAGCGGTGGCTTACAGGGCGCAGCACGCGAGTCGGCGCTGGAAGTGCTGGAACTGACGGCGGGGAAAACTGGCGCCCTTTT‑TGATTCTTCAACCCGGTTTCCTGCTGGACCAAAATCACCGGt    >  2:495246/1‑147 (MQ=255)                                                                                                                                                   |                                                                                                                                                 TACCACCATGATGGCCAGCTGCCTCGCGGTTTTCGCACCTGAGACGATCAACAGCCAAACCTTCCGCGACGTGGCGGATCGTTGCCAGGCGATCCTGACCTCACTGGGCGATTTCAGCGAAGGTGTGTTTGGTTACGCACCGTGGAAACGGATCGTTTATCTCGGTAGCGGTGGCTTACAGGGCGCAGCACGCGAGTCGGCGCTGAAAGTGCTGGAACTGACGGC‑GGGTAAACTGGCGGCCTTTTATGATTCTCCAACCGGATTCCGTCATGGACCAAAATCGCTGGTCG  >  NC_000913/3282503‑3282792

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