Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,568,907	C→T	G141D (GGC→GAC)	glgC ←	glucose‑1‑phosphate adenylyltransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,568,907	0	C	T	100.0%	182.8 / NA	47	G141D (GGC→GAC)	glgC	glucose‑1‑phosphate adenylyltransferase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (23/24);  total (23/24)

GATAATTTTATCGTTCTCATCAACCGCCATAACGCCAAATGCGGAGGCTTCTTCAATCGGTACTGGCATACAAGCAACGGTGCAACGTGCGCCTTTTTCGACGTGATCGATAAGCATACGCGAGTAGTCTTGCTTGTAGATATGGTCGCCCGCCAGGATCACCACGTATTCCGCTTTATAACGGCGGATAATGTCGAGGTTTTGGGTGACCGCATCTGCGGTGCCGCGATACCAGTTTTCCCCTTTCATTCTCTGCTGTGCTGGCAGCAGATCGACAAACTCGTTCATTTCTTCA  >  NC_000913/3568759‑3569053                                                                                                                                                     |                                                                                                                                                   gATAATTTTATCGTTCTCATCAACCGCCATAACGCCAAATGCGGAGGCTTCTTCAATCGGTACTGGCATACAAGCAACGGTGCAACGTGCGCCTTTTTCGACGTGATCGATAAGCATACGCGAGTAGTCTTGCTTGTAGATATGGTcgtcc                                                                                                                                                  >  1:704385/1‑151 (MQ=255) gATAATTTTATCGTTCTCATCAACCGCCATAACGCCAAATGCGGAGGCTTCTTCAATCGGTACTGGCATACAAGCAACGGTGCAACGTGCGCCTTTTTCGACGTGATCGATAAGCATACGCGAGTAGTCTTGCTTGTAGATATGGTcgtcc                                                                                                                                                  >  1:475430/1‑151 (MQ=255)   tAATTTTATCGTTCTCATCAACCGCCATAACGCCAAATGCGGAGGCTTCTTCAATCGGTACTGGCATACAAGCAACGGTGCAACGTGCGCCTTTTTCGACGTGATCGATAAGCATACGCGAGTAGTCTTGCTTGTAGATATGGTcgtccgc                                                                                                                                                >  2:605210/1‑151 (MQ=255)   tAATTTTATCGTCCTCATCAACCGCCATAACGCCAAATGCGGAGGCTTCTTCAATCGGTACTGGCATACAAGCAACGGTGCAACGTGCGCCTTTTTCGACGTGATCGATAAGCATACGCGAGTAGTCTTGCTTGTAGATATGGTcgtccg                                                                                                                                                 <  2:669677/150‑1 (MQ=255)      ttttATCTTTCTAAACAGCCGCCATAACGCCAAATGCGGAGGCTTCTTCAATCGGTACTGGCATACAAGCAACGGTGCAACGTGCGCCTTTTTCGACGTGATCGATAAGCATACGCGAGTAGTCTTGCTTGTAGATATGGTcgtccgc                                                                                                                                                <  1:279729/148‑1 (MQ=255)            cGTTCTCATCAACCGCCATAACGCCAAATGCGGAGGCTTCTTCAATCGGTACTGGCATACAAGCAACGGTGCAACGTGCGCCTTTTTCGACGTGATCGATAAGCATACGCGAGTAGTCTTGCTTGTAGATATGGTCGTCCGCCAGGATcac                                                                                                                                       >  2:589131/1‑151 (MQ=255)             gTTCTCATCAACCGCCATAACGCCAAATGCGGAGGCTTCTTCAATCGGTACTGGCATACAAGCAACGGTGCAACGTGCGCCTTTTTCGACGTGATCGATAAGCATACGCGAGTAGTCTTGCTTGTAGATATGGTCGTCCGCCAGGATcacc                                                                                                                                      >  2:40660/1‑151 (MQ=255)              ttCTCATCAACCGCCATAACGCCAAATGCGGAGGCTTCTTCAATCGGTACTGGCATACAAGCAACGGTGCAACGTGCGCCTTTTTCGACGTGATCGATAAGCATACGCGAGTAGTCTTGCTTGTAGATATGGTCGTCCGCCAGGATcacc                                                                                                                                      >  2:581396/1‑150 (MQ=255)                   atcaCCCGCCATAACGCCAAATGCGGAGGCTTCTTCAATCGGTACTGGCATACAAGCAACGGTGCAACGTGCGCCTTTTTCGACGTGATCGATAAGCATACGCGAGTAGTCTTGCTTGTAGATATGGTCGTCCGCCAGGATCACCACGTa                                                                                                                                 <  2:70806/150‑1 (MQ=255)                    tcaACCGCCATAACGCCAAATGCGGAGGCTTCTTCAATCGGTACTGGCATACAAGCAACGGTGCAACGTGCGCCTTTTTCGACGTGATCGATAAGCATACGCGAGTAGTCTTGCTTGTAGATATGGTCGTCCGCCAGGATCACCACGTAtt                                                                                                                               >  1:13676/1‑151 (MQ=255)                     caACAGCCATAACGCCAAATGCGGAGGCTTCTTCAATCGGTACTGGCATACAAGCAACGGTGCAACGTGCGCGTTTCTCGACGTGATCGATAAGCATACGCGAGTAGTCTTGCTTGTAGATATGGTCGTCCGCCAGGATCACCACGTATTc                                                                                                                              <  2:587545/151‑1 (MQ=255)                      aaCCGCCATAACGCCAAATGCGGAGGCTTCTTCAATCGGTACTGGCATACAAGCAACGGTGCAACGTGCGCCTTTTTCGACGTGATCGATAAGCATACGCGAGTAGTCTTGCTTGTAGATATGGTCGTCCGCCAGGATCACCACGTATTcc                                                                                                                             >  2:190548/1‑151 (MQ=255)                         cGCCATAACGCCAAATGCGGAGGCTTCTTCAATCGGTACTGGCATACAAGCAACGGTGCAACGTGCGCCTTTTTCGACGTGATCGATAAGCATACGCGAGTAGTCTTGCTTGTAGATATGGTCGTCCGCCAGGATCACCACGTATTCCGCt                                                                                                                          >  1:70701/1‑151 (MQ=255)                          gCCATAACGCCAAATGCGGAGGCTTCTTCAATCGGTACTGGCATACAAGCAACGGTGCAACGTGCGCCTTTTTCGACGTGATCGATAAGCATACGCGAGTAGTCTTGCTTGTAGATATGGTCGTCCGCCAGGATCACCACGTATTCCGCtt                                                                                                                         <  1:717778/151‑1 (MQ=255)                              tAACGCCAAATGCGGAGGCTTCTTCAATCGGTACTGGCATACAAGCAACGGTGCAACGTGCGCCTTTTTCGACGTGATCGATAAGCATACGCGAGTAGTCTTGCTTGTAGATATGGTCGTCCGCCAGGATCACCACGTATTCCGCTTtat                                                                                                                      >  2:456604/1‑150 (MQ=255)                                      aaTGCGGAGGCTTCTTCAATCGGTACTGGCATACAAGCAACGGTGCAACGTGCGCCTTTTTCGACGTGATCGATAAGCATACGCGAGTAGTCTTGCTTGTAGATATGGTCGTCCGCCAGGATCACCACGTATTCCGCTTTATAACGGCGGa                                                                                                             <  1:295316/151‑1 (MQ=255)                                          cGGAGGCTTCTTCAATCGGTACTGGCATACAAGCAACGGTGCAACGTGCGCCTTTTTCGACGTGATCGATAAGCATACGCGAGTAGTCTTGCTTGTAGATATGGTCGTCCGCCAGGATCACCACGTATTCCGCTTTATAAcggcgg                                                                                                              >  2:343180/1‑146 (MQ=255)                                                  tcttcACTCGGCACTGGCATACAAGCAACGGTGCAACGTGCGCCTTTTTCGACGGGATCGATAAGCATACGCGAGTAGTCTTGCTTGCAGATATGGTCGTCCGCCAGGCTCACCACGTATTCCGCCTTATAACGGCGGTTAAt                                                                                                         >  1:57284/1‑143 (MQ=255)                                                      cAATCGGTACTGGCATACAAGCAACGGTGCAACGTGCGCCTTTTTCGACGTGATCGATAAGCATACGCGAGTAGTCTTGCTTGTAGATATGGTCGTCCGCCAGGATCACCACGTATTCCGCTTTATAACGGCGGATAATGTCGAGGtttt                                                                                              <  1:463145/150‑1 (MQ=255)                                                            gTACTGGCATACAAGCAACGGTGCAACGTGCGCCTTTTTCGACGTGATCGATAAGCATACGCGAGTAGTCTTGCTTGTAGATATGGTCGTCCGCCAGGATCACCACGTATTCCGCTTTATAACGGCGGATAATGTCGAGGTTTTGGGTg                                                                                         >  1:671014/1‑149 (MQ=255)                                                             tACTGGCATACAAGCAACGGTGCAACGTGCGCCTTTTTCGACGTGATCGATAAGCATACGCGAGTAGTCTTGCTTGTAGATATGGTCGTCCGCCAGGATCACCACGTATTCCGCTTTATAACGGCGGATAATGTCGAGGTTTTGGGTGAcc                                                                                      >  1:699880/1‑151 (MQ=255)                                                                    aTACAAGCAACGGTGCAACGTGCGCCTTTTTCGACGTGATCGATAAGCATACGCGAGTAGTCTTGCTTGTAGATATGGTCGTCCGCCAGGATCACCACGTATTCCGCTTTATAACGGCGGATAATGTCGAGGTTTTGGGTGACCGCATCt                                                                                >  1:554509/1‑150 (MQ=255)                                                                              cGGTGCAACGTGCGCCTTTTTCGACGTGATCGATAAGCATACGCGAGTAGTCTTGCTTGTAGATATGGTCGTCCGCCAGGATCACCACGTATTCCGCTTTATAACGGCGGATAATGTCGAGGTTTTGGGTGACCGCATCTGCGGTGCcgcg                                                                     <  1:772971/151‑1 (MQ=255)                                                                                  gCAACGTGCGCCTTTTTCGACGTGATCGATAAGCATACGCGAGTAGTCTTGCTTGTAGATATGGTCGTCCGCCAGGATCACCACGTATTCCGCTTTATAACGGCGGATAATGTCGAGGTTTTGGGTGACCGCATCTGCGGTGCCGCGATa                                                                  <  2:344307/150‑1 (MQ=255)                                                                                   cAACGTGCGCCTTTTTCGACGTGATCGATAAGCATACGCGAGTAGTCTTGCTTGTAGATATGGTCGTCCGCCAGGATCACCACGTATTCCGCTTTATAACGGCGGATAATGTCGAGGTTTTGGGTGACCGCATCTGCGGTGCCGCGATAcc                                                                <  2:461385/151‑1 (MQ=255)                                                                                       gTGCGCCTTTTTCGACGTGATCGATAAGCATACGCGAGTAGTCTTGCTTGTAGATATGGTCGTCCGCCAGGATCACCACGTATTCCGCTTTATAACGGCGGATAATGTCGAGGTTTTGGGTGACCGCATCTGCGGTGCCGCGATACCAGtt                                                            >  2:94769/1‑151 (MQ=255)                                                                                                    gACGTGATCGATAAGCATACGCGAGTAGTCTTGCTTGTAGATATGGTCGTCCGCCAGGATCACCACGTATTCCGCTTTATAACGGCGGATAATGTCGAGGTTTTGGGTGACCGCATCTGCGGTGCCGCGATACCAGTTTTCCCCTTTCAtt                                               <  2:267539/151‑1 (MQ=255)                                                                                                    gACGTGATCGATAAGCATACGCGAGTAGTCTTGCTTGTAGATATGGTCGTCCGCCAGGATCACCACGTATTCCGCTTTATAACGGCGGATAATGTCGAGGTTTTGGGTGACCGCATCTGCGGTGCCGCGATACCAGTTTTCCCCTTTCAtt                                               >  1:360320/1‑151 (MQ=255)                                                                                                         gATCGATAAGCATACGCGAGTAGTCTTGCTTGTAGATATGGTCGTCCGCCAGGATCACCACGTATTCCGCTTTATAACGGCGGATAATGTCGAGGTTTTGGGTGACCGCATCTGCGGTGCCGCGATACCAGTTTTCCCCTTTCATTctct                                           <  2:747377/150‑1 (MQ=255)                                                                                                            cGATAAGCATACGCGAGTAGTCTTGCTTGTAGATATGGTCGTCCGCCAGGATCACCACGTATTCCGCTTTATAACGGCGGATAATGTCGAGGTTTTGGGTGACCGCATCTGCGGTGCCGCGATACCAGTTTTCCCCTTTCATTCTctgct                                        >  2:480338/1‑150 (MQ=255)                                                                                                            cGATAAGCATACGCGAGTAGTCTTGCTTGTAGATATGGTCGTCCGCCAGGATCACCACGTATTCCGCTTTATAACGGCGGATAATGTCGAGGTTTTGGGTGACCGCATCTGCGGGGCCGCGATACCAGTTTTCCCCTTTCATTCTctgctg                                       >  2:421632/1‑151 (MQ=255)                                                                                                                      aCGCGAGTAGTCTTGCTTGTAGATATGGTCGTCCGCCAGGATCACCACGTATTCCGCTTTATAACGGCGGATAATGTCGAGGTTTTGGGTGACCGCATCTGCGGTGCCGCGATACCAGTTTTCCCCTTTCATTCTCTGCTGTGCTGgcag                              <  1:132194/150‑1 (MQ=255)                                                                                                                      aCGCGAGTAGTCTTGCTTGTAGATATGGTCGTCCGCCAGGATCACCACGTATTCCGCTTTATAACGGCGGATAATGTCGAGGTTTTGGGTGACCGCATCTGCGGTGCCGCGATACCAGTTTTCCCCTTTCATTCTCTGCTGTGCTGgcag                              <  2:519188/150‑1 (MQ=255)                                                                                                                           agtagtCTTGCTTGTAGATATGGTCGTCCGCCAGGATCACCACGTATTCCGCTTTATAACGGCGGATAATGTCGAGGTTTTGGGTGACCGCATCTGCGGTGCCGCGATACCAGTTTTCCCCTTTCATTCTCTGCTGTGCTGGCAGCAGATc                        <  2:616287/151‑1 (MQ=255)                                                                                                                           agtagtCTTGCTTGTAGATATGGTCGTCCGCCAGGATCACCACGTATTCCGCTTTATAACGGCGGATAATGTCGAGGGTTTGGGTGACCGCATCTGCGGGGCCGCGATACCAGGTTTCCCCTTTCATTCTCTGCTGTGCTGgcagc                             >  1:299577/1‑146 (MQ=255)                                                                                                                            gtagtCTTGCTTGTAGATATGGTCGTCCGCCAGGATCACCACGTATTCCGCTTTATAACGGCGGATAATGTCGAGGTTTTGGGTGACCGCATCTGCGGTGCCGCGATACCAGTTTTCCCCTTTCATTCTCTGCTGTGCTGGCAGCAGATCg                       <  1:52641/151‑1 (MQ=255)                                                                                                                            gtagtCTTGCTTGTAGATATGGTCGTCCGCCAGGATCACCACGTATTCCGCTTTATAACGGCGGATAATGTCGAGGTTTTGGGTGACCGCATCTGCGGTGCCGCGATACCAGTTTTCCCCTTTCATTCTCTGCTGTGCTGGCAGCAGATCg                       <  1:419607/151‑1 (MQ=255)                                                                                                                              agtCTTGCTTGTAGATATGGTCGTCCGCCAGGATCACCACGTATTCCGCTTTATAACGGCGGATAATGTCGAGGTTTTGGGTGACCGCATCTGCGGTGCCGCGATACCAGTTTTCCCCTTTCATTCTCTGCTGTGCTGGCAGCAGATCGa                      <  1:190515/150‑1 (MQ=255)                                                                                                                              agtCTTGCTTGTAGATATGGTCGTCCGCCAGGATCACCACGTATTCCGCTTTATAACGGCGGATAATGTCGAGGTTTTGGGTGACCGCATCTGCGGTGCCGCGATACCAGTTTTCCCCTTTCATTCTCTGCTGTGCTGGCAGCAGATCGAc                     >  2:542878/1‑151 (MQ=255)                                                                                                                               gtCTTGCTTGTAGATATGGTCGTCCGCCAGGATCACCACGTATTCCGCTTTATAACGGCGGATAATGTCGAGGTTTTGGGTGACCGCATCTGCGGTGCCGCGATACCAGTTTTCCCCTTTCATTCTCTGCTGTGCTGGCAGCAGATCGACa                    <  1:287828/151‑1 (MQ=255)                                                                                                                                   tgcttgTAGATATGGTCGTCCGCCAGGATCACCACGTATTCCGCTTTATAACGGCGGATAATGTCGAGGTTTTGGGTGACCGCATCTGCGGTGCCGCGATACCAGTTTTCCCCTTTCATTCTCTGCTGTGCTGGCAGCAGATCGACaaa                  <  2:464723/149‑1 (MQ=255)                                                                                                                                   tgcttgTAGATATGGTCGTCCGCCAGGATCACCACGTATTCCGCTTTATAACGGCGGATAATGTCGAGGTTTTGGGTGACCGCATCTGCGGTGCCGCGATACCAGTTTTCCCCTTTCATTCTCTGCTGTGCTGGCAGCAGATCGACaa                   <  1:94738/148‑1 (MQ=255)                                                                                                                                   tgcttgTAGATATGGTCGTCCGCCAGGATCACCACGTATTCCGCTTTATAACGGCGGATAATGTCGAGGTTTTGGGTGACCGCATCTGCGGTGCCGCGATACCAGTTTTCCCCTTTCATTCTCTGCTGTGCTGGCAGCAGATCGACAAAc                 <  2:137146/150‑1 (MQ=255)                                                                                                                                   tgcttgTAGATATGGTCGTCCGCCAGGATCACCACGTATTCCGCTTTATAACGGCGGATAATGTCGAGGTTTTGGGTGACCGCATCTGCGGTGCCGCGATACCAGTTTTCCCCTTTCATTCTCTGCTGTGCTGGCAGCAGATCGACAAAc                 <  2:751684/150‑1 (MQ=255)                                                                                                                                   tgcttgTAGATATGGTCGTCCGCCAGGATCACCACGTATTCCGCTTTATAACGGCGGATAATGTCGAGGGTTTGGGTGACCGCATCTGCGGTGCCGCGATACCAGTTTTCCCCTTTCATTCTCTGCTGTGCTGGCAGCAGATCGACAAAc                 >  1:548257/1‑150 (MQ=255)                                                                                                                                    gcttgTAGATATGGTCGTCCGCCAGGATCACCACGTATTCCGCTTTATAACGGCGGATAATGTCGAGGTTTTGGGTGACCGCATCTGCGGTGCCGCGATACCAGTTTTCCCCTTTCATTCTCTGCTGTGCTGGCAGCAGATCGACAAACTc               >  1:162389/1‑151 (MQ=255)                                                                                                                                     cttgTAGATATGGTCGTCCGCCAGGATCACCACGTATTCCGCTTTATAACGGCGGATAATGTCGAGGTTTTGGGTGACCGCATCTGCGGTGCCGCGATACCAGTTTTCCCCTTTCATTCTCTGCTGTGCTGGCAGCAGATCGACAAACTc               >  1:749501/1‑150 (MQ=255)                                                                                                                                        gTAGATATGGTCGTCCGCCAGGATCACCACGTATTCCGCTTTATAACGGCGGATAATGTCGAGGTTTTGGGTGACCGCATCTGCGGTGCCGCGATACCAGTTTTCCCCTTTCATTCTCTGCTGTGCTGGCAGCAGATCGACAAACTCGtt            <  2:332073/150‑1 (MQ=255)                                                                                                                                        gTAGATATGGTCGTCCGCCAGGATCACCACGTATTCCGCTTTATAACGGCGGATAATGTCGAGGTTTTGGGTGACCGCATCTGCGGTGCCGCGATACCAGTTTTCCCCTTTCATTCTCTGCTGTGCTGGCAGCAGATCGACAAACTCGTTc           <  1:16177/151‑1 (MQ=255)                                                                                                                                              atGGTCGTCCGCCAGGATCACCACGTATTCCGCTTTATAACGGCGGATAATGTCGAGGTTTTGGGTGACCGCATCTGCGGTGCCGCGATACCAGTTTTCCCCTTTCATTCTCTGCTGTGCTGGCAGCAGATCGACAAACTCGTTCATttct     >  1:334926/1‑151 (MQ=255)                                                                                                                                              atGGTCGTCCGCCAGGATCACCACGTATTCCGCTTTATAACGGCGGATAATGTCGAGGTTTTGGGTGACCGCATCTGCGGTGCCGCGATACCAGTTTTCCCCTTTCATTCTCTGCTGTGCTGGCAGCAGATCGACAAACTCGTTCATttc      <  1:178168/150‑1 (MQ=255)                                                                                                                                                  tCGTCCGCCAGGATCACCACGTATTCCGCTTTATAACGGCGGATAATGTCGAGGTTTTGGGTGACCGCATCTGCGGTGCCGCGATACCAGTTTTCCCCTTTCATTCTCTGCTGTGCTGGCAGCAGATCGACAAACTCGTTCATTTCTTCa  >  2:45309/1‑150 (MQ=255)                                                                                                                                                     |                                                                                                                                                   GATAATTTTATCGTTCTCATCAACCGCCATAACGCCAAATGCGGAGGCTTCTTCAATCGGTACTGGCATACAAGCAACGGTGCAACGTGCGCCTTTTTCGACGTGATCGATAAGCATACGCGAGTAGTCTTGCTTGTAGATATGGTCGCCCGCCAGGATCACCACGTATTCCGCTTTATAACGGCGGATAATGTCGAGGTTTTGGGTGACCGCATCTGCGGTGCCGCGATACCAGTTTTCCCCTTTCATTCTCTGCTGTGCTGGCAGCAGATCGACAAACTCGTTCATTTCTTCA  >  NC_000913/3568759‑3569053

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