Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,604,202	G→A	*242* (TAG→TAA)  E15K (GAA→AAA)	yjjQ → bglJ →	putative transcriptional regulator bgl operon transcriptional activator

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,604,202	0	G	A	100.0%	240.9 / NA	74	*242* (TAG→TAA)  E15K (GAA→AAA)	yjjQ bglJ	putative transcriptional regulator bgl operon transcriptional activator
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (39/35);  total (39/35)

GCAACGCTGAAATGTATGCATGGATAAATAGCGCGCAGGGTGCAAGAGAACTTAACTTGCCTTCTGTTTATGGAGATGCCGCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAGAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGT  >  NC_000913/4604075‑4604335                                                                                                                                |                                                                                                                                      gCAACGCTGAAATGTATGCATGGATAAATAGCGCGCAGGGTGCAAGAGAACTTAACTTGCCTTCTGTTTATGGAGATGCCGCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAAAAAAAt                                                                                                                                 >  2:930328/1‑134 (MQ=255) gCAACGCTGAAATGTATGCATGGATAAATAGCGCGCAGGGTGCAAGAGAACTTAACTTGCCTTCTGTTTATGGAGATGCCGCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAAAAAAAt                                                                                                                                 <  1:930328/134‑1 (MQ=255)     cGCTGAAATGTATGCATGGATAAATAGCGCGCAGGGTGCAAGAGAACTTAACTTGCCTTCTGTTTATGGAGATGCCGCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAAAAAAATGCGTCATGa                                                                                                                        >  1:1468237/1‑139 (MQ=255)             tGTATGCATGGATAAATAGCGCGCAGGGTGCAAGAGAACTTAACTTGCCTTCTGTTTATGGAGATGCCGCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAAAAAAATGCGTCATGAGTAGTATc                                                                                                                <  2:1468237/139‑1 (MQ=255)               tATGCATGGATAAATAGCGCGCAGGGTGCAAGAGAACTTAACTTGCCTTCTGTTTATGGAGATGCCGCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAAAAAAATGCGTCATGAgtagta                                                                                                                  <  1:1061131/135‑1 (MQ=255)               tATGCATGGATAAATAGCGCGCAGGGTGCAAGAGAACTTAACTTGCCTTCTGTTTATGGAGATGCCGCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAAAAAAATGCGTCATGAgtagta                                                                                                                  >  2:1061131/1‑135 (MQ=255)                      ggATAAATAGCGCGCAGGGTGCAAGAGAACTTAACTTGCCTTCTGTTTATGGAGATGCCGCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAAAAAAATGCGTCATGAGTAGTATCGGTATTg                                                                                                         >  1:1446401/1‑137 (MQ=255)                      ggATAAATAGCGCGCAGGGTGCAAGAGAACTTAACTTGCCTTCTGTTTATGGAGATGCCGCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAAAAAAATGCGTCATGAGTAGTATCGGTATTg                                                                                                         <  2:1446401/137‑1 (MQ=255)                           aaTAGCGCGCAGGGTGCAAGAGAACTTAACTTGCCTTCTGTTTATGGAGATGCCGCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAAAAAAATGCGTAATGAGTAGTATCGGTATTGAGAGttt                                                                                                  >  1:1652510/1‑139 (MQ=255)                                cgcgcAGGGTGCAAGAGAACTTAACTTGCCTTCTGTTTATGGAGATGCCGCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAAAAAAATGCGTCATGAGTAGTATCg                                                                                                               >  1:555909/1‑121 (MQ=255)                                cgcgcAGGGTGCAAGAGAACTTAACTTGCCTTCTGTTTATGGAGATGCCGCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAAAAAAATGCGTCATGAGTAGTATCg                                                                                                               <  2:555909/121‑1 (MQ=255)                                cgcgcAGGGTGCAAGAGAACTTAACTTGCCTTCTGTTTATGGAGATGCCGCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAAAAAAATGCGTCATGAGTAGTATCGGTAtt                                                                                                          <  2:1580219/126‑1 (MQ=255)                                cgcgcAGGGTGCAAGAGAACTTAACTTGCCTTCTGTTTATGGAGATGCCGCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAAAAAAATGCGTCATGAGTAGTATCGGTAtt                                                                                                          >  1:1580219/1‑126 (MQ=255)                                  cgcAGGGTGCAAGAGAACTTAACTTGCCTTCTGTTTATGGAGATGCCGCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAAAAAAATGCGTCATGag                                                                                                                       >  2:1599212/1‑111 (MQ=255)                                  cgcAGGGTGCAAGAGAACTTAACTTGCCTTCTGTTTATGGAGATGCCGCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAAAAAAATGCGTCATGag                                                                                                                       <  1:1599212/111‑1 (MQ=255)                                  cgcAGGGTGCAAGAGAACTTAACTTGCCTTCTGTTTATGGAGATGCCGCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAAAAAAATGCGTCATGAgtagta                                                                                                                  >  2:181008/1‑116 (MQ=255)                                  cgcAGGGTGCAAGAGAACTTAACTTGCCTTCTGTTTATGGAGATGCCGCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAAAAAAATGCGTCATGAgtagta                                                                                                                  <  1:181008/116‑1 (MQ=255)                                  cgcAGGGTGCAAGAGAACTTAACTTGCCTTCTGTTTATGGAGATGCCGCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAAAAAAATGCGTCATGAGTAGTAt                                                                                                                 >  1:1045552/1‑117 (MQ=255)                                  cgcAGGGTGCAAGAGAACTTAACTTGCCTTCTGTTTATGGAGATGCCGCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAAAAAAATGCGTCATGAGTAGTAt                                                                                                                 <  2:1045552/117‑1 (MQ=255)                                       ggTGCAAGAGAACTTAACTTGCCTTCTGTTTATGGAGATGCCGCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAAAAAAATGCGTCATGAGTAGTATCg                                                                                                               >  2:359222/1‑114 (MQ=255)                                       ggTGCAAGAGAACTTAACTTGCCTTCTGTTTATGGAGATGCCGCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAAAAAAATGCGTCATGAGTAGTATCg                                                                                                               <  1:359222/114‑1 (MQ=255)                                              gagaACTTAACTTGCCTTCTGTTTATGGAGATGCCGCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAAAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGt                                                                               <  2:1652510/139‑1 (MQ=255)                                                   cttaacttGCCTTCTGTTTATGGAGATGCCGCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAAAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGAAAAcca                                                                          >  1:383404/1‑138 (MQ=255)                                                             cTTCTGTTTATGGAGATGCCGCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAAAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTc                                                                <  2:383404/139‑1 (MQ=255)                                                                  gTTTATGGAGATGCCGCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAAAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATAc                                                           <  1:557801/139‑1 (MQ=255)                                                                   tttATGGAGATGCCGCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAAAAAAATGCGTCATGag                                                                                                                       >  1:573786/1‑78 (MQ=255)                                                                   tttATGGAGATGCCGCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAAAAAAATGCGTCATGag                                                                                                                       <  2:573786/78‑1 (MQ=255)                                                                           gaTGCCGCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAACAAAATGCGTCATGag                                                                                                                       <  1:1115144/70‑1 (MQ=255)                                                                           gaTGCCGCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAAAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTc                                                                >  1:1083440/1‑125 (MQ=255)                                                                           gaTGCCGCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAAAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTc                                                                <  2:1083440/125‑1 (MQ=255)                                                                           gaTGCCGCAGAATGGAACAAAGCCGAATTCAGAAGAGAAATGTCGCACTCATAAAAAAATGCGTCATGag                                                                                                                       >  2:1115144/1‑70 (MQ=255)                                                                            aTGCCGCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAAAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGt                                                 >  1:1152913/1‑139 (MQ=255)                                                                             tGCCGCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAAAAAAATGCGTCATGagtag                                                                                                                    >  1:1264961/1‑71 (MQ=255)                                                                             tGCCGCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAAAAAAATGCGTCATGagtag                                                                                                                    <  2:1264961/71‑1 (MQ=255)                                                                             tGCCGCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAAAAAAATGCGTCATGag                                                                                                                       <  2:1594028/68‑1 (MQ=255)                                                                             tGCCGCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAAAAAAATGCGTCATGag                                                                                                                       >  1:1594028/1‑68 (MQ=255)                                                                             tGCCGCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAAAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGtt                                                                                      >  2:443575/1‑101 (MQ=255)                                                                             tGCCGCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAAAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGtt                                                                                      <  1:443575/101‑1 (MQ=255)                                                                             tGCCGCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAAAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTc                                                >  1:980107/1‑139 (MQ=255)                                                                               ccGCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAAAAAAATGCGTCATGag                                                                                                                       >  1:245078/1‑66 (MQ=255)                                                                               ccGCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAAAAAAATGCGTCATGag                                                                                                                       <  2:245078/66‑1 (MQ=255)                                                                                cGCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAAAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCggg                                                                                >  1:62401/1‑104 (MQ=255)                                                                                cGCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAAAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCggg                                                                                <  2:62401/104‑1 (MQ=255)                                                                                          aaCACAGCCGAATTAAGAAGAGGAATGTCGCACTCATAAAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGa                                            <  1:1134851/130‑1 (MQ=255)                                                                                          aaCACAGCCGAATTAAGAAGAGGAATGTCGCACTCATAAAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGa                                            >  2:1134851/1‑130 (MQ=255)                                                                                          aaCACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAAAAAAATGCGTCATTAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCggg                                                                                <  2:1026806/94‑1 (MQ=255)                                                                                          aaCACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAAAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCggg                                                                                >  1:1026806/1‑94 (MQ=255)                                                                                            cacaGCCGAATTAAGAAGAGAAATGTCGCACTCATAAAAAAATGCGTc                                                                                                                            >  2:415216/1‑48 (MQ=255)                                                                                            cacaGCCGAATTAAGAAGAGAAATGTCGCACTCATAAAAAAATGCGTc                                                                                                                            <  1:415216/48‑1 (MQ=255)                                                                                               aGCCGAATTAAGAAGAGAAATGTCGCACTCATAAAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCggg                                                                                >  2:423741/1‑89 (MQ=255)                                                                                               aGCCGAATTAAGAAGAGAAATGTCGCACTCATAAAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCggg                                                                                <  1:423741/89‑1 (MQ=255)                                                                                               aGCCGAATTAAGAAGAGAAATGTCGCACTCATAAAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTcc                                                               <  1:1653142/106‑1 (MQ=255)                                                                                               aGCCGAATTAAGAAGAGAAATGTCGCACTCATAAAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTcc                                                               >  2:1653142/1‑106 (MQ=255)                                                                                               aGCCGAATTAAGAAGAGAAATGTCGCACTCATAAAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATg                              >  2:970265/1‑139 (MQ=255)                                                                                                 ccGAATTAAGAAGAGAAATGTCGCACTCATAAAAAAATGCg                                                                                                                              <  2:1302187/41‑1 (MQ=255)                                                                                                 ccGAATTAAGAAGAGAAATGTCGCACTCATAAAAAAATGCg                                                                                                                              >  1:1302187/1‑41 (MQ=255)                                                                                                 ccGAATTAAGAAGAGAAATGTCGCACTCATAAAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTc                                                                <  2:512914/103‑1 (MQ=255)                                                                                                 ccGAATTAAGAAGAGAAATGTCGCACTCATAAAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTc                                                                >  1:512914/1‑103 (MQ=255)                                                                                                  cGAATTAAGAAGAGAAATGTCGCACTCATAAAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCa                           <  2:980107/139‑1 (MQ=255)                                                                                                    aaTTAAGAAGAGAAATGTCGCACTCATAAAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTc                                                                >  2:1525182/1‑100 (MQ=255)                                                                                                    aaTTAAGAAGAGAAATGTCGCACTCATAAAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTc                                                                <  1:1525182/100‑1 (MQ=255)                                                                                                     aTTAAGAAGAGAAATGTCGCACTCATAAAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGt                        <  1:970265/139‑1 (MQ=255)                                                                                                      ttAAGAAGAGAAATGTCGCACTCATAAAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTc                       >  1:1523702/1‑139 (MQ=255)                                                                                                              agaAATGTCGCACTCATAAAAAAATGCGTCATGAGTAGTATCGGTATTg                                                                                                         >  1:824838/1‑49 (MQ=255)                                                                                                              agaAATGTCGCACTCATAAAAAAATGCGTCATGAGTAGTATCGGTATTg                                                                                                         <  2:824838/49‑1 (MQ=255)                                                                                                                  aTGTCGCACTCATAAAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGtt           >  1:670030/1‑139 (MQ=255)                                                                                                                      cGCACTCATAAAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTgcg       <  2:424864/139‑1 (MQ=255)                                                                                                                        cACTCATAAAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTgcggc     >  1:656551/1‑139 (MQ=255)                                                                                                                          cTCATAAAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGtt                                                                                      <  1:686417/56‑1 (MQ=255)                                                                                                                          cTCATAAAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGtt                                                                                      >  2:686417/1‑56 (MQ=255)                                                                                                                          cTCATAAAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGcggcgg   >  2:742959/1‑139 (MQ=255)                                                                                                                           tCATAAAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGGTTATTCAGAAAGTTTGCGGGTAACCCTTCTAAGCTCCATCCCTTTACCAGTCAGGAGTCATGTCAGGATGCCATGTCGAGGAACTCGTTTGCGGCGGt  >  1:424864/1‑139 (MQ=255)                                                                                                                             aTAAAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGa                                            <  2:1449691/95‑1 (MQ=255)                                                                                                                             aTAAAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGa                                            >  1:1449691/1‑95 (MQ=255)                                                                                                                                |                                                                                                                                      GCAACGCTGAAATGTATGCATGGATAAATAGCGCGCAGGGTGCAAGAGAACTTAACTTGCCTTCTGTTTATGGAGATGCCGCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAGAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGT  >  NC_000913/4604075‑4604335

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