Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,530,092	A→G	V51A (GTC→GCC)	sgcC ←	putative PTS enzyme IIC component SgcC

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,530,092	0	A	G	100.0%	217.4 / NA	62	V51A (GTC→GCC)	sgcC	putative PTS enzyme IIC component SgcC
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (30/32);  total (30/32)

GATAATCATGGCGCCAATGGCTGTCGCATAACCAACGCCGGAGGCAGGGCCTGCGCCAACGTCAAAGACATGCAGGGCTAGGCCGAAACGCTCAATCATCACTTTGATAGGGGGACTGAGGCTGTCGATCGCCATGACGATCACCAGTCCCATACCGACAAAACCAATCCCAACCGTTACCCCTGCCTTAATCGCCTGCAGCCAGGGGATGCGGAAGATAAGGCCAATCACGATCATGATGATAGGGACGAAAACGGTACCGCCGAGAGAC  >  NC_000913/4529956‑4530226                                                                                                                                         |                                                                                                                                       gATAATCATGGCGCCAATGGCTGTCGCATAACCAACGCCGGAGGCAGGGCCTGCGCCAACGTCAAAGACATGCAGGGCTAGGCCGAAACGCTCAATCATCACTTTGATAGGGGGACTGAGGCTGTCGATCGCCATGGCGa                                                                                                                                     >  1:590479/1‑140 (MQ=255)     aTCATGGCGCCAATGGCTGTCGCATAACCAACGCCGGAGGCAGGGCCTGCGCCAACGTCAAAGACATGCAGGGCTAGGCCGAAACGCTCAATCATCACTTTGATAGGGGGACTGAGGCTGTCGATCGCCATGGCGATCAc                                                                                                                                 <  2:380548/140‑1 (MQ=255)     aTCATGGCGCCAATGGCTGTCGCATAACCAACGCCGGAGGCAGGGCCTGCGCCAACGTCAAAGACATGCAGGGCTAGGCCGAAACGCTCAATCATCACTTTGATAGGGGGACTGAGGCTGTCGATCGCCATGGCGATCAc                                                                                                                                 <  1:575415/140‑1 (MQ=255)     aTCATGGCGCCAATGGCTGTCGCATAACCAACGCCGGAGGCAGGGCCTGCGCCAACGTCAAAGACATGCAGGGCTAGGCCGAAACGCTCAATCATCACTTTGATAGGGGGACTGAGGCTGTCGATCGCCATGGCGATCAc                                                                                                                                 <  1:739989/140‑1 (MQ=255)     aTCATGGCGCCAATGGCTGTCGCATAACCAACGCCGGAGGCAGGGCCTGCGCCAACGTCAAAGACATGCAGGGCTAGGCCGAAACGCTCAATCATCACTTTGATAGGGGGACTGAGGCTGTCGATCGCCATGGCGATCAc                                                                                                                                 <  1:437411/140‑1 (MQ=255)          ggCGCCAATGGCTGTCGCATAACCAACGCCGGAGGCAGGGCCTGCGCCAACGTCAAAGACATGCAGGGCTAGGCCGAAACGCTCAATCATCACTTTGATAGGGGGACTGAGGCTGTCGATCGCCATGGCGATCACCAGTc                                                                                                                            >  1:573319/1‑140 (MQ=255)           gcgcCAATGGCTGTCGCATAACCAACGCCGGAGGCAGGGCCTGCGCCAACGTCAAAGACATGCAGGGCTAGGCCGAAACGCTCAATCATCACTTTGATAGGGGGACTGAGGCTGTCGATCGCCATGGCGATCAc                                                                                                                                 <  1:387831/134‑1 (MQ=255)           gcgcCAATGGCTGTCGCATAACCAACGCCGGAGGCAGGGCCTGCGCCAACGTCAAAGACATGCAGGGCTAGGCCGAAACGCTCAATCATCACTTTGATAGGGGGACTGAGGCTGTCGATCGCCATGGCGATCAc                                                                                                                                 >  2:387831/1‑134 (MQ=255)           gcgcCAATGGCTGTCGCATAACCAACGCCGGAGGCAGGGCCTGCGCCAACGTCAAAGACATGCAGGGCTAGGCCGAAACGCTCAATCATCACTTTGATAGGGGGACTGAGGCTGTCGATCGCCATGGCGATCACCAGTcc                                                                                                                           >  1:644990/1‑140 (MQ=255)                    gCTGTCGCATAACCAACGCCGGAGGCAGGGCCTGCGCCAACGTCAAAGACATGCAGGGCTAGGCCGAAACGCTCAATCATCACTTTGATAGGGGGACTGAGGCTGTCGATCGCCATGGCGATCACCAGTCCCATACCGAc                                                                                                                  <  2:668093/140‑1 (MQ=255)                       gTCGCATAACCAACGCCGGAGGCAGGGCCTGCGCCAACGTCAAAGACATGCAGGGCTAGGCCGAAACGCTCAATCATCACTTTGATAGGGGGACTGAGGCTGTCGATCGCCATGGCGATCACCAGTCCCATACCGACaaa                                                                                                               >  2:633096/1‑140 (MQ=255)                          gCATAACCAACGCCGGAGGCAGGGCCTGCGCCAACGTCAAAGACATGCAGGGCTAGGCCGAAACGCTCAATCATCACTTTGATAGGGGGACTGAGGCTGTCGATCGCCATGGCGATCACCAGTCCCATACCGACAAAAcc                                                                                                            <  1:714467/140‑1 (MQ=255)                             tAACCAACGCCGGAGGCAGGGCCTGCGCCAACGTCAAAGACATGCAGGGCTAGGCCGAAACGCTCAATCATCACTTTGATAGGGGGACTGAGGCTGTCGATCGCCATGGCGATCACCAGTCCCATACCGACAAAACCAAt                                                                                                         >  2:865097/1‑140 (MQ=255)                               aCCAACGCCGGAGGCAGGGCCTGCGCCAACGTCAAAGACATGCAGGGCTAGGCCGAAACGCTCAATCATCACTTTGATAGGGGGACTGAGGCTGTCGATCGCCATGGCGATCACCAGTCCCATACCGACAAAACCAATcc                                                                                                       <  2:433281/140‑1 (MQ=255)                                ccAACGCCGGAGGCAGGGCCTGCGCCAACGTCAAAGACATGCAGGGCTAGGCCGAAACGCTCAATCATCACTTTGATAGGGGGACTGAGGCTGTCGATCGCCATGGCGATCACCAGTCCCATACCGACAAAACCAATccc                                                                                                      <  2:342684/140‑1 (MQ=255)                                 cAACGCCGGAGGCAGGGCCTGCGCCAATGTCAAAGACATGCAGGGCTAGGCCGAAACGCTCAATCATCACTTTGATAGGGGGACTGAGGCTGTCGATCGCCATGGCGATCACCAGTCCCATACCGACAAAACCAATCCCa                                                                                                     <  2:862330/140‑1 (MQ=255)                                 cAACGCCGGAGGCAGGGCCTGCGCCAACGTCAAAGACATGCAGGGCTAGGCCGAAACGCTCAATCATCACTTTGATAGGGGGACTGAGGCTGTCGATCGCCATGGCGATCACCAGTCCCATACCGACAAAACCAATCCCa                                                                                                     <  2:590479/140‑1 (MQ=255)                                 cAACGCCGGAGGCAGGGCCTGCGCCAACGTCAAAGACATGCAGGGCTAGGCCGAAACGCTCAATCATCACTTTGATAGGGGGACTGAGGCTGTCGATCGCCATGGCGATCACCAGTCCCATACCGACAAAACCAATCCCa                                                                                                     <  2:148995/140‑1 (MQ=255)                                      ccGGAGGCAGGGCCTGCGCCAACGTCAAAGACATGCAGGGCTAGGCCGAAACGCTCAATCATCACTTTGATAGGGGGACTGAGGCTGTCGATCGCCATGGCGATCACCAGTCCCATACCGACAAAACCAATCCCAACCGt                                                                                                <  2:846552/140‑1 (MQ=255)                                         gAGGCAGGGCCTGCGCCAACGTCAAAGACATGCAGGGCTAGGCCGAAACGCTCAATCATCACTTTGATAGGGGGACTGAGGCTGTCGATCGCCATGGCGATCACCAGTCCCATACCGACAAAACCAATCCCAACCGTTAc                                                                                             >  2:87339/1‑140 (MQ=255)                                                 gCCTGCGCCAACGTCAAAGACATGCAGGGCTAGGCCGAAACGCTCAATCATCACTTTGATAGGGGGACTGAGGCTGTCGATCGCCATGGCGATCACCAGTCCCATACCGACAAAACCAATCCCAACCGTTACCCCTGCCt                                                                                     <  1:87339/140‑1 (MQ=255)                                                     gcgcCAACGTCAAAGACATGCAGGGCTAGGCCGAAACGCTCAATCATCACTTTGATAGGGGGACTGAGGCTGTCGATCGCCATGGCGATCACCAGTCCCATACCGACAAAACCAATCCCAACCGTTACCCCTGCCTTAAt                                                                                 >  1:539422/1‑140 (MQ=255)                                                      cgcCAACGTCAAAGACATGCAGGGCTAGGCCGAAACGCTCAATCATCACTTTGATAGGGGGACTGAGGCTGTCGATCGCCATGGCGATCACCAGTCCCATACCGACAAAACCAATCCCAACCGTTACCCCTGCCTTAATc                                                                                >  1:110806/1‑140 (MQ=255)                                                        ccAACGTCAAAGACATGCAGGGCTAGGCCGAAACGCTCAATCATCACTTTGATAGGGGGACTGAGGCTGTCGATCGCCATGGCGATCACCAGTCCCATACCGACAAAACCAATCCCAACCGTTACCCCTGCCTTAATCGc                                                                              <  1:379820/140‑1 (MQ=255)                                                         cAACGTCAAAGACATGCAGGGCTAGGCCGAAACGCTCAATCATCACTTTGATAGGGGGACTGAGGCTGTCGATCGCCATGGCGATCACCAGTCCCATACCGACAAAACCAATCCCAACCGTTACCCCTGCCTTAATCGcc                                                                             <  2:539422/140‑1 (MQ=255)                                                          aaCGTCAAAGACATGCAGGGCTAGGCCGAAACGCTCAATCATCACTTTGATAGGGGGACTGAGGCTGTCGATCGCCATGGCGATCACCAGTCCCATACCGACAAAACCAATCCCAACCGTTACCCCTGCCTTAATCGCCt                                                                            <  1:658911/140‑1 (MQ=255)                                                           aCGTCAAAGACATGCAGGGCTAGGCCGAAACGCTCAATCATCACTTTGATAGGGGGACTGAGGCTGTCGATCGCCATGGCGATCACCAGTCCCATACCGACAAAACCAATCCCAACCGTTACCCCTGCCTTAATCGcc                                                                             >  1:822338/1‑138 (MQ=255)                                                           aCGTCAAAGACATGCAGGGCTAGGCCGAAACGCTCAATCATCACTTTGATAGGGGGACTGAGGCTGTCGATCGCCATGGCGATCACCAGTCCCATACCGACAAAACCAATCCCAACCGTTACCCCTGCCTTAATCGcc                                                                             <  2:822338/138‑1 (MQ=255)                                                           aCGTCAAAGACATGCAGGGCTAGGCCGAAACGCTCAATCATCACTTTGATAGGGGGACTGAGGCTGTCGATCGCCATGGCGATCACCAGTCCCATACCGACAAAACCAATCCCAACCGTTACCCCTGCCTTAATCGCCTg                                                                           >  2:581340/1‑140 (MQ=255)                                                            cGTCAAAGACATGCAGGGCTAGGCCGAAACGCTCAATCATCACTTTGATAGGGGGACTGAGGCTGTCGATCGCCATGGCGATCACCAGTCCCATACCGACAAAACCAATCCCAACCGTTACCCCTGCCTTAATCGCCTGc                                                                          >  1:117301/1‑140 (MQ=255)                                                                aaaGACATGCAGGGCTAGGCCGAAACGCTCAATCATCACTTTGATAGGGGGACTGAGGCTGTCGATCGCCATGGCGATCACCAGTCCCATACCGACAAAACCAATCCCAACCGTTACCCCTGCCTTAATCGCCTGCAGcc                                                                      >  1:461929/1‑140 (MQ=255)                                                                       tGCAGGGCTAGGCCGAAACGCTCAATCATCACTTTGATAGGGGGACTGAGGCTGTCGATCGCCATGGCGATCACCAGTCCCATACCGACAAAACCAATCCCAACCGTTACCCCTGCCTTAATCGCCTGCAGCCAGGGGAt                                                               >  1:210531/1‑140 (MQ=255)                                                                          aGGGCTAGGCCGAAACGCTCAATCATCACTTTGATAGGGGGACTGAGGCTGTCGATCGCCATGGCGATCACCAGTCCCATACCGACAAAACCAATCCCAACCGTTACCCCTGCCTTAATCGCCTGCAGCCAGGGGATGCg                                                            <  1:470760/140‑1 (MQ=255)                                                                           gggCTAGGCCGAAACGCTCAATCATCACTTTGATAGGGGGACTGAGGCTGTCGATCGCCATGGCGATCACCAGTCCCATACCGACAAAACCAATCCCAACCGTTACCCCTGCCTTAATCGCCTGCAGCCAGGGGATGCgg                                                           <  1:468976/140‑1 (MQ=255)                                                                           gggCTAGGCCGAAACGCTCAATCATCACTTTGATAGGGGGACTGAGGCTGTCGATCGCCATGGCGATCACCAGTCCCATACCGACAAAACCAATCCCAACCGTTACCCCTGCCTTAATCGCCTGCAGCCAGGGGATGCgg                                                           <  1:148329/140‑1 (MQ=255)                                                                             gCTAGGCCGAAACGCTCAATCATCACTTTGATAGGGGGACTGAGGCTGTCGATCGCCATGGCGATCACCAGTCCCATACCGACAAAACCAATCCCAACCGTTACCCCTGCCTTAATCGCCTGCAGCCAGGGGATGCGGaa                                                         >  1:649670/1‑140 (MQ=255)                                                                             gCTAGGCCGAAACGCTCAATCATCACTTTGATAGGGGGACTGAGGCTGTCGATCGCCATGGCGATCACCAGTCCCATACCGACAAAACCAATCCCAACCGTTACCCCTGCCTTAATCGCCTGCAGCCAGGGGATGCGGaa                                                         >  1:250643/1‑140 (MQ=255)                                                                              cTAGGCCGAAACGCTCAATCATCACTTTGATAGGGGGACTGAGGCTGTCGATCGCCATGGCGATCACCAGTCCCATACCGACAAAACCAATCCCAACCGTTACCCCTGCCTTAATCGCCTGCAGCCAGGGGATGCGGAAg                                                        >  1:488183/1‑140 (MQ=255)                                                                                aGGCCGAAACGCTCAATCATCACTTTGATAGGGGGACTGAGGCTGTCGATCGCCATGGCGATCACCAGTCCCATACCGACAAAACCAATCCCAACCGTTACCCCTGCCTTAATCGCCTGCAGCCAGGGGATGCGGAAGAt                                                      >  2:388711/1‑140 (MQ=255)                                                                                   ccGAAACGCTCAATCATCACTTTGATAGGGGGACTGAGGCTGTCGATCGCCATGGCGATCACCAGTCCCATACCGACAAAACCAATCCCAACCGTTACCCCTGCCTTAATCGCCTGCAGCCAGGGGATGCGGAAGATAAg                                                   >  2:466283/1‑140 (MQ=255)                                                                                    cGAAACGCTCAATCATCACTTTGATAGGGGGACTGAGGCTGTCGATCGCCATGGCGATCACCAGTCCCATACCGACAAAACCAATCCCAACCGTTACCCCTGCCTTAATCGCCTGCAGCCAGGGGATGCGGAAGATAAgg                                                  <  2:649670/140‑1 (MQ=255)                                                                                     gAAACGCTCAATCATCACTTTGATAGGGGGACTGAGGCTGTCGATCGCCATGGCGATCACCAGTCCCATACCGACAAAACCAATCCCAACCGTTACCCCTGCCTTAATCGCCTGCAGCCAGGGGATGCGGAAGATAAGGc                                                 <  2:488183/140‑1 (MQ=255)                                                                                         cGCTCAATCATCACTTTGATAGGGGGACTGAGGCTGTCGATCGCCATGGCGATCACCAGTCCCATACCGACAAAACCAATCCCAACCGTTACCCCTGCCTTAATCGCCTGCAGCCAGGGGATGCGGAAGATAAGGCCAAt                                             <  2:461929/140‑1 (MQ=255)                                                                                          gCTCAATCATCACTTTGATAGGGGGACTGAGGCTGTCGATCGCCATGGCGATCACCAGTCCCATACCGACAAAACCAATCCCAACCGTTACCCCTGCCTTAATCGCCTGCAGCCAGGGGATGCGGAAGATAAGGCCAATc                                            <  1:215368/140‑1 (MQ=255)                                                                                           cTCAATCATCACTTTGATAGGGGGACTGAGGCTGTCGATCGCCATGGCGATCACCAGTCCCATACCGACAAAACCAATCCCAACCGTTACCCCTGCCTTAATCGCCTGCAGCCAGGGGATGCg                                                            <  2:242709/123‑1 (MQ=255)                                                                                           cTCAATCATCACTTTGATAGGGGGACTGAGGCTGTCGATCGCCATGGCGATCACCAGTCCCATACCGACAAAACCAATCCCAACCGTTACCCCTGCCTTAATCGCCTGCAGCCAGGGGATGCg                                                            >  1:242709/1‑123 (MQ=255)                                                                                           cTCAATCATCACTTTGATAGGGGGACTGAGGCTGTCGATCGCCATGGCGATCACCAGTCCCATACCGACAAAACCAATCCCAACCGTTACCCCTGCCTTAATCGCCCGCAGCCAGGGGATGCGGAAGATAAGGCCAATCa                                           >  2:472838/1‑140 (MQ=255)                                                                                             cAATCATCACTTTGATAGGGGGACTGAGGCTGTCGATCGCCATGGCGATCACCAGTCCCATACCGACAAAACCAATCCCAACCGTTACCCCTGCCTTAATCGCCTGCAGCCAGGGGATGCGGAAGATAAGGCCAATCACg                                         <  2:69773/140‑1 (MQ=255)                                                                                                  atcaCTTTGATAGGGGGACTGAGGCTGTCGATCGCCATGGCGATCACCAGTCCCATACCGACAAAACCAATCCCAACCGTTACCCCTGCCTTAATCGCCTGCAGCCAGGGGATGCGGAAGATAAGGCCAATCACGATCat                                    <  1:388711/140‑1 (MQ=255)                                                                                                        ttGATAGGGGGACTGAGGCTGTCGATCGCCATGGCGATCACCAGTCCCATACCGACAAAACCAATCCCAACCGTTACCCCTGCCTTAATCGCCTGCAGCCAGGGGATGCGGAAGATAAGGCCAATCACGATCATgatgat                              <  1:865097/140‑1 (MQ=255)                                                                                                        ttGATAGGGGGACTGAGGCTGTCGATCGCCATGGCGATCACCAGTCCCATACCGACAAAACCAATCCCAACCGTTACCCCTGCCTTAATCGCCTGCAGCCAGGGGATGCGGAAGATAAGGCCAATCACGATCATgatgat                              >  2:319772/1‑140 (MQ=255)                                                                                                          gATAGGGGGACTGAGGCTGTCGATCGCCATGGCGATCACCAGTCCCATACCGACAAAACCAATCCCAACCGTTACCCCTGCCTTAATCGCCTGCAGCCAGGGGATGCGGAAGATAAGGCCAATCACg                                         <  2:406727/127‑1 (MQ=255)                                                                                                          gATAGGGGGACTGAGGCTGTCGATCGCCATGGCGATCACCAGTCCCATACCGACAAAACCAATCCCAACCGTTACCCCTGCCTTAATCGCCTGCAGCCAGGGGATGCGGAAGATAAGGCCAATCACg                                         >  1:406727/1‑127 (MQ=255)                                                                                                            tAGGGGGACTGAGGCTGTCGATCGCCATGGCGATCACCAGTCCCATACCGACAAAACCAATCCCAACCGTTACCCCTGCCTTAATCGCCTGCAGCCAGGGGATGCGGAAGATAAGGCCAATc                                            >  2:678181/1‑122 (MQ=255)                                                                                                            tAGGGGGACTGAGGCTGTCGATCGCCATGGCGATCACCAGTCCCATACCGACAAAACCAATCCCAACCGTTACCCCTGCCTTAATCGCCTGCAGCCAGGGGATGCGGAAGATAAGGCCAATc                                            <  1:678181/122‑1 (MQ=255)                                                                                                               ggggACTGAGGCTGTCGATCGCCATGGCGATCACCAGTCCCATACCGACAAAACCAATCCCAACCGTTACCCCTGCCTTAATCGCCTGCAGCCAGGGGATGCGGAAGATAAGGCCAATCACGATCATGATGATAGGGACg                       >  2:656934/1‑140 (MQ=255)                                                                                                                   aCTGAGGCTGTCGATCGCCATGGCGATCACCAGTCCCATACCGACAAAACCAATCCCAACCGTTACCCCTGCCTTAATCGCCTGCAGCCAGGGGATGCGGAAGATAAGGCCAATCACGATCATGATGATAGGGACGaaaa                   >  1:414720/1‑140 (MQ=255)                                                                                                                     tGAGGCTGTCGATCGCCATGGCGATCACCAGTCCCATACCGACAAAACCAATCCCAACCGTTACCCCTGCCTTAATCGCCTGCAGCCAGGGGATGCGGAAGATAAGGCCAATCACGATCATGATGATAGGGACGAAAACg                 >  2:510245/1‑140 (MQ=255)                                                                                                                         gCTGTCGATCGCCATGGCGATCACCAGTCCCATACCGACAAAACCAATCCCAACCGTTACCCCTGCCTTAATCGCCTGCAGCCAGGGGATGCGGAAGATAAGGCCAATCACGATCATGATGATAGGGACGAAAACGGTAc             >  2:705376/1‑140 (MQ=255)                                                                                                                             tCGATCGCCATGGCGATCACCAGTCCCATACCGACAAAACCAATCCCAACCGTTACCCCTGCCTTAATCGCCTGCAGCCAGGGGATGCGGAAGATAAGGCCAATCACGATCATGATGATAGGGACGAAAACGGTAccgcc         >  2:693662/1‑140 (MQ=255)                                                                                                                             tCGATCGCCATGGCGATCACCAGTCCCATACCGACAAAACCAATCCCAACCGTTACCCCTGCCTTAATCGCCTGCAGCCAGGGGATGCGGAAGATAAGGCCAATCACGATCATGATGATAGGGACGAAAACGGTAccgcc         >  2:241260/1‑140 (MQ=255)                                                                                                                               gATCGCCATGGCGATCACCAGCCCCATCCCGACAAAACCAACCCCAACCGTTCCCCCTGCCTTAATCGCCTGCAGCCAGGGGATGCGGAAGATAGGGCAAATCACGATCATGATGATAGGGAGGAAAACGGTCCCGCCga       <  2:573319/140‑1 (MQ=255)                                                                                                                                    ccATGGCGATCACCAGTCCCATACCGACAAAACCAATCCCAACCGTTACCCCTGCCTTAATCGCCTGCAGCCAGGGGATGCGGAAGATAAGGCCAATCACGATCATGATGATAGGGACGAAAACGGTACCGCCGAGAGAc  <  1:705376/140‑1 (MQ=255)                                                                                                                                         |                                                                                                                                       GATAATCATGGCGCCAATGGCTGTCGCATAACCAACGCCGGAGGCAGGGCCTGCGCCAACGTCAAAGACATGCAGGGCTAGGCCGAAACGCTCAATCATCACTTTGATAGGGGGACTGAGGCTGTCGATCGCCATGACGATCACCAGTCCCATACCGACAAAACCAATCCCAACCGTTACCCCTGCCTTAATCGCCTGCAGCCAGGGGATGCGGAAGATAAGGCCAATCACGATCATGATGATAGGGACGAAAACGGTACCGCCGAGAGAC  >  NC_000913/4529956‑4530226

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