Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,076,627	G→C	S161T (AGT→ACT)	yihY →	BrkB family putative transporter, inner membrane protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,076,627	0	G	C	100.0%	207.0 / NA	58	S161T (AGT→ACT)	yihY	BrkB family putative transporter, inner membrane protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base C (25/33);  total (25/33)

GATGTACTCCATCGATAGCGCGTTGAATACCATCTGGCGCAGTAAACGAGCGCGACCAAAAATTTACTCGTTCGCCGTGTACTGGATGATTTTAACGCTGGGGCCGCTGCTGGCAGGGGCCAGTCTGGCGATCAGTTCCTATTTGCTCTCCCTGCGCTGGGCGAGCGATCTCAATACTGTCATCGACAACGTGCTGCGTATTTTTCCGTTGCTGTTGTCGTGGATCTCCTTCTGGTTGCTGTACAGCATTGTTCCTACC  >  NC_000913/4076493‑4076751                                                                                                                                       |                                                                                                                             gATGTACTCCATCGATAGCGCGTTGAATACCATCTGGCGCAGTAAACGAGCGCGACCAAAAATTTACTCGTTCGCCGTGTACTGGATGATTTTAACGCTGGGGCCGCTGCTGGCAGGGGCCAGTCTGGCGATCACTTcc                                                                                                                          >  2:438784/1‑139 (MQ=255)    gTACTCCATCGATAGCGCGTTGAATACCATCTGGCGCAGTAAACGAGCGCGACCAAAAATTTACTCGTTCGCCGTGTACTGGATGATTTTAACGCTGGGGCCGCTGCTGGCAGGGGTCAGTCTGGCGATCACTTCCTAt                                                                                                                       >  1:1113549/1‑139 (MQ=255)        tCCATCGATAGCGCGTTGAATACCATCTGGCGCAGTAAACGAGCGCGACCAAAAATTTACTCGTTCGCCGTGTACTGGATGATTTTAACGCTGGGGCCGCTGCTGGCAGGGGCCAGTCTGGCGATCACTTCCTATTTGc                                                                                                                   <  2:1383589/139‑1 (MQ=255)         ccATCGATAGCGCGTTGAATACCATCTGGCGCAGTAAACGAGCGCGACCAAAAATTTACTCGTTCGCCGTGTACTGGATGATTTTAACGCTGGGGCCGCTGCTGGCAGGGGCCAGTCTGGCGATCACTTCCTATTTGct                                                                                                                  <  2:1658396/139‑1 (MQ=255)         ccATCGATAGCGCGTTGAATACCATCTGGCGCAGTAAACGAGCGCGACCAAAAATTTACTCGTTCGCCGTGTACTGGATGATTTTAACGCTGGGGCCGCTGCTGGCAGGGGCCAGTCTGGCGATCACTTCCTATTTGct                                                                                                                  >  2:1439309/1‑139 (MQ=255)         ccATCGATAGCGCGTTGAATACCATCTGGCGCAGTAAACGAGCGCGACCAAAAATTTACTCGTTCGCCGTGTACTGGATGATTTTAACGCTGGGGCCGCTGCTGGCAGGGGCCAGTCTGGCGATCACTTCCTATTTGct                                                                                                                  <  1:1274153/139‑1 (MQ=255)         ccATCGATAGCGCGTTGAATACCATCTGGCGCAGTAAACGAGCGCGACCAAAAATTTACTCGTTCGCCGTGTACTGGATGATTTTAACGCTGGGGCCGCTGCTGGCAGGGGCCAGTCTGGCGATCACTTCCTATTTGct                                                                                                                  <  1:438784/139‑1 (MQ=255)         ccATCGATAGCGCGTTGAATACCATCTGGCGCAGTAAACGAGCGCGACCAAAAATTTACTCGTTCGCCGTGTACTGGATGATTGTAACGCTGGGGCCGCTGCTGGCAGGGGCCAGTCTGGCGATCACTTCCTATTTGct                                                                                                                  >  2:977111/1‑139 (MQ=255)            tCGATAGCGCGTTGAATAACATCTGGCGCAGTAAACGAGCGCGACCAAAAATTTACTCGTTCGCCGTGTACTGGATGATTTTAACGCTGGGGCCGCTGCTGGCAGGGGCCAGTCTGGCGATCACTTCCTATTTGCtctc                                                                                                               >  1:1019677/1‑139 (MQ=255)                  gcgcGTTGAATACCATCTGGCGCAGTAAACGAGCGCGACCAAAAATTTACTCGTTCGCCGTGTACTGGATGATTTTAACGCTGGGGCCGCTGCTGGCAGGGGCCAGTCTGGCGATCACTTCCTATTTGCTCTCCCTgcg                                                                                                         >  2:1033334/1‑139 (MQ=255)                   cgcgTTGAATACCATCTGGCGCAGTAAACGAGCGCGACCAAAAATTTACTCGTTCGCCGTGTACTGGATGATTTTAACGCTGGGGCCGCTGCTGGCAGGGGCCAGTCTGGCGATCACTTCCTATTTGCTCTCCCTgcgc                                                                                                        >  2:969014/1‑139 (MQ=255)                       ttGAATACCATCTGGCGCAGTAAACGAGCGCGACCAAAAATTTACTCGTTCGCCGTGTACTGGATGTTTTTAACGCTGGGGCCGCTGCTGGCAGGGGCCAGTCTGGCGATCACTTCCTATTTGCTCTCCCTGCGCTggg                                                                                                    <  1:1512515/139‑1 (MQ=255)                          aaTACCATCTGGCGCAGTAAACGAGCGCGACCAAAAATTTACTCGTTCGCCGTGTACTGGATGATTTTAACGCTGGGGCCGCTGCTGGCAGGGGCCAGTCTGGCGATCACTTCCTATTTGCTCTCCCTGCGCTGGgcga                                                                                                 <  1:1033334/139‑1 (MQ=255)                          aaTACCATCTGGCGCAGTAAACGAGCGCGACCAAAAATTTACTCGTTCGCCGTGTACTGGATGATTTTAACGCTGGGGCCGCTGCTGGCAGGGGCCAGTCTGGCGATCACTTCCTATTTGCTCTCCCTGCGCTGGgcga                                                                                                 >  2:497905/1‑139 (MQ=255)                           aTACCATCTGGCGCAGTAAACGAGCGCGACCAAAAATTTACTCGTTCGCCGTGTACTGGATGATTTTAACGCTGGGGCCGCTGCTGGCAGGGGCCAGTCTGGCGATCACTTCCTATTTGCTCTCCCTGCGCTGGgcgag                                                                                                <  1:497905/139‑1 (MQ=255)                            tACCATCTGGCGCAGTAAACGAGCGCGACCAAAAATTTACTCGTTCGCCGTGTACTGGATGATTTTAACGCTGGGGCCGCTGCTGGCAGGGGCCAGTCTGGCGATCACTTCCTATTTGCTCTCCCTGCGCTGGgcgagc                                                                                               <  2:1019677/139‑1 (MQ=255)                               cATCTGGCGCAGTAAACGAGCGCGACCAAAAATTTACTCGTTCGCCGTGTACTGGATGATTTTAACGCTGGGGCCGCTGCTGGCAGGGGCCAGTCTGGCGATCACTTCCTATTTGCTCTCCCTGCGCTGGGCGAGCGAt                                                                                            >  2:1104505/1‑139 (MQ=255)                                    ggCGCAGTAAACGAGCGCGACCAAAAATTTACTCGTTCGCCGTGTACTGGATGATTTTAACGCTGGGGCCGCTGCTGGCAGGGGCCAGTCTGGCGATCACTTCCTATTTGCTCTCCCTGCGCTGGgcgagc                                                                                               <  1:163057/131‑1 (MQ=255)                                    ggCGCAGTAAACGAGCGCGACCAAAAATTTACTCGTTCGCCGTGTACTGGATGATTTTAACGCTGGGGCCGCTGCTGGCAGGGGCCAGTCTGGCGATCACTTCCTATTTGCTCTCCCTGCGCTGGgcgagc                                                                                               >  2:163057/1‑131 (MQ=255)                                    ggCGCAGTAAACGAGCGCGACCAAAAATTTACTCGTTCGCCGTGTACTGGATGATTTTAACGCTGGGGCCGCTGCTGGCAGGGGCCAGTCTGGCGATCACTTCCTATTTGCTCTCCCTGCGCTGGGCGAGCGATCTCaa                                                                                       >  2:1254507/1‑139 (MQ=255)                                      cgcAGTAAACGAGCGCGACCAAAAATTTACTCGTTCGCCGTGTACTGGATGATTTTAACGCTGGGGCCGCTGCTGGCAGGGGCCAGTCTGGCGATCACTTCCTATTTGCTCTCCCTGCGCTGGGCGAGCGATCTCAATa                                                                                     >  2:59436/1‑139 (MQ=255)                                      cgcAGTAAACGAGCGCGACCAAAAATTTACTCGTTCGCCGTGTACTGGATGATTTTAACGCTGGGGCCGCTGCTGGCAGGGGCCAGTCTGGCGATCACTTCCTATTTGCTCTCCCTGCGCTGGGCGAGCGATCTCAATa                                                                                     <  2:1464973/139‑1 (MQ=255)                                                         cAAAAATTTACTCGTTCGCCGTGTACTGGATGATTTTAACGCTGGGGCCGCTGCTGGCAGGGGCCAGTCTGGCGATCACTTCCTATTTGCTCTCCCTGCGCTGGGCGAGCGATCTCAATACTGTCATCGACAACGtgct                                                                  <  1:977111/139‑1 (MQ=255)                                                            aaaTTTACTCGTTCGCCGTGTACTGGATGATTTTAACGCTGGGGCCGCTGCTGGCAGGGGCCAGTCTGGCGATCACTTCCTa                                                                                                                        >  1:1033932/1‑82 (MQ=255)                                                            aaaTTTACTCGTTCGCCGTGTACTGGATGATTTTAACGCTGGGGCCGCTGCTGGCAGGGGCCAGTCTGGCGATCACTTCCTa                                                                                                                        <  2:1033932/82‑1 (MQ=255)                                                                  aCTCGTTCGCCGTGTACTGGATGATTTTAACGCTGGGGCCGCTGCTGGCAGGGGCCAGTCTGGCGATCACTTCCTATTTGCTCTCCCTGCGCTGGGCGAGCGATCTCAATACTGTCATCGACAACGTGCTGCGTAtttt                                                         >  1:172398/1‑139 (MQ=255)                                                                       ttCGCCGTGTACTGGATGATTTTAACGCTGGGGCCGCTGCTGGCAGGGGCCAGTCTGGCGATCACTTCCTATTTGCTCTCCCTGCGCTGGGCGAGCGATCTCAATACTGTCATCGACAACGTGCTGCGTAttt                                                          <  1:1545558/133‑1 (MQ=255)                                                                       ttCGCCGTGTACTGGATGATTTTAACGCTGGGGCCGCTGCTGGCAGGGGCCAGTCTGGCGATCACTTCCTATTTGCTCTCCCTGCGCTGGGCGAGCGATCTCAATACTGTCATCGACAACGTGCTGCGTAttt                                                          >  2:1545558/1‑133 (MQ=255)                                                                            cGTGTACTGGATGATTTTAACGCTGGGGCCGCTGCTGGCAGGGGCCAGTCTGGCGATCACTTCCTATTTGCTCTCCCTGCGCTGGGCGAGCGATCTCAATACTGTCATCGACAACGTGCTGCGTATTTTTCCGTTGCtg                                               <  1:1439309/139‑1 (MQ=255)                                                                                        gatTTTAACGCTGGGGCCGCTGCTGGCAGGGGCCAGTCTGGCGATCACTTCCTATTTGCTCTCCCTGCGCTGGGCGAGCGATCTCAATACTGTCATCGACAACGTGCTGCGTATTTTTCCGTTGCTGTTGTCGTGGAtc                                   >  1:1608010/1‑139 (MQ=255)                                                                                          ttttAACGCTGGGGCCGCTGCTGGCAGGGGCCAGTCTGGCGATCACTTCCTATTTGCTCTCCCTGCGCTGGGCGAGCGATCTCAATACTGTCATCGACa                                                                         >  2:562189/1‑99 (MQ=255)                                                                                          ttttAACGCTGGGGCCGCTGCTGGCAGGGGCCAGTCTGGCGATCACTTCCTATTTGCTCTCCCTGCGCTGGGCGAGCGATCTCAATACTGTCATCGACa                                                                         <  1:562189/99‑1 (MQ=255)                                                                                                  cTGGGGCCGCTGCTGGCAGGGGCCAGTCTGGCGATCACTTCCTATTTGCTCTCCCTGCGCTGGGCGAGCGATCTCAATACTGTCATCGACAACGTGCTGCGTATTTTTCCGTTGCTGTTGTCGTGGATCTCCTTCTGGt                         <  1:990447/139‑1 (MQ=255)                                                                                                    ggggCCGCTGCTGGCAGGGGCCAGTCTGGCGATCACTTCCTATTTGCTCTCCCTGCGCTGGGCGAGCGATCTCAATACTGTCATCGACAACGtgct                                                                  >  2:96482/1‑96 (MQ=255)                                                                                                    ggggCCGCTGCTGGCAGGGGCCAGTCTGGCGATCACTTCCTATTTGCTCTCCCTGCGCTGGGCGAGCGATCTCAATACTGTCATCGACAACGtgct                                                                  <  1:96482/96‑1 (MQ=255)                                                                                                    ggggCCGCTGCTGGCAGGGGCCAGTCTGGCGATCACTTCCTATTTGCTCTCCCTGCGCTGGGCGAGCGATCTCAATACTGTCATCGACAACGTGCTGCGTATTTTTCCGTTGCTGTTGTCGTGGa                                     <  1:36673/125‑1 (MQ=255)                                                                                                    ggggCCGCTGCTGGCAGGGGCCAGTCTGGCGATCACTTCCTATTTGCTCTCCCTGCGCTGGGCGAGCGATCTCAATACTGTCATCGACAACGTGCTGCGTATTTTTCCGTTGCTGTTGTCGTGGa                                     >  2:36673/1‑125 (MQ=255)                                                                                                    ggggCCGCTGCTGGCAGGGGCCAGTCTGGCGATCACTTCCTATTTGCTCTCCCTGCGCTGGGCGAGCGATCTCAATACTGTCATCGACAACGTGCTGCGTATTTTTCCGTTGCTGTTGTCGTGGATCTCCTTCTGGTTg                       <  1:1104505/139‑1 (MQ=255)                                                                                                     gggCCGCTGCTGGCAGGGGCCAGTCTGGCGATCACTTCCTATTTGCTCTCCCTGCGCTGGGCGAGCGATCTCAATa                                                                                     >  1:1236773/1‑76 (MQ=255)                                                                                                     gggCCGCTGCTGGCAGGGGCCAGTCTGGCGATCACTTCCTATTTGCTCTCCCTGCGCTGGGCGAGCGATCTCAATa                                                                                     <  2:1236773/76‑1 (MQ=255)                                                                                                     gggCCGCTGCTGGCAGGGGCCAGTCTGGCGATCACTTCCTATTTGCTCTCCCTGCGCTGGGCGAGCGATCTCAATACTGTCATCGACAACGTGCTGCGTATTTTTCCGTTGCt                                                <  1:1243652/113‑1 (MQ=255)                                                                                                     gggCCGCTGCTGGCAGGGGCAAGTCTGGCGATCACTTCCTATTTGCTCTCCCTGCGCTGGGCGAGCGATCTCAATACTGTCATCGACAACGTGCTGCGTATTTTTCCGTTGCt                                                >  2:1243652/1‑113 (MQ=255)                                                                                                             gctgGCAGGGGCCAGTCTGGCGATCACTTCCTATTTGCTCTCCCTGCGCTGGGCGAGCGATCTCAATACTGTCATCGACAACGTGCTGCGTATTTTTCCGTTGCTGTTGTCGTGGATCTCCTTCTGGTTGCTGTACAGc              <  2:763622/139‑1 (MQ=255)                                                                                                              ctgGCAGGGGTCAGTCTGGCGATCACTTCCTATTTGCTCTCCCTGCGCTGGGCGAGCGATCTCAATACTGTCATCGACAACGTGCTGCGTATTTTTCCGTTGCTGTTGTCGTGGATCTCCTTCTGGTTGCTGTACAGCa             <  2:1113549/139‑1 (MQ=255)                                                                                                              ctgGCAGGGGCCAGTCTGGCGATCACTTCCTATTTGCTCTCCCTGCGCTGGGCGAGCGATCTCAATACTGTCATCGACAACGTGCTGCGTATTTTTCCGTTGCTGTTGTCGTGGATCTCCTTCTGGTTGCTGTACAGCa             <  2:1608010/139‑1 (MQ=255)                                                                                                              ctgGCAGGGGCCAGTCTGGCGATCACTTCCTATTTGCTCTCCCTGCGCTGGGCGAGCGATCTCAATACTGTCATCGACAACGTGCTGCGTATTTTTCCGTTGCTGTTGTCGTGGATCTCCTTCTGGTTGCTGTACAGCa             <  2:172398/139‑1 (MQ=255)                                                                                                              ctgGCAGGGGCCAGTCTGGCGATCACTTCCTATTTGCTCTCCCTGCGCTGGGCGAGCGATCTCAATACTGTCATCGACAACGTGCTGCGTATTTTTCCGTTGCTGTTGTCGTGGATCTCCTTCTGGTTGCTGTACAGCa             <  1:969014/139‑1 (MQ=255)                                                                                                                 gCAGGGGCCAGTCTGGCGATCACTTCCTATTTGCTCTCCCTGCGCTGGGCGAGCGATCTCAATACTGTCATCGACAACGTGCTGCGTATTTTTCCGTTGCTGTTGTCGTGGATCTCCTTCTGGTTGCTGTACAGCATTg          <  2:263018/139‑1 (MQ=255)                                                                                                                     gggCCAGTCTGGCGATCACTTCCTATTTGCTCTCCCTGCGCTGGGCGAGCGATCTCAATACTGTCATCGACAACGTGCTGCGTATTTTTCCGTTGCTGTTGTCGTGGATCTCCTTCTGGTTGCTGTACAGCATTGTTcc      <  1:1721075/139‑1 (MQ=255)                                                                                                                        ccAGTCTGGCGATCACTTCCTATTTGCTCTCCCTGCGCTGGgcgag                                                                                                <  1:1760618/46‑1 (MQ=255)                                                                                                                        ccAGTCTGGCGATCACTTCCTATTTGCTCTCCCTGCGCTGGgcgag                                                                                                >  2:1760618/1‑46 (MQ=255)                                                                                                                        ccAGTCTGGCGATCACTTCCTATTTGCTCTCCCTGCGCTGGGCGAGCGATCTCAATACTGTCATCGACAACGTGCTGCGTATTTTTCCGTTGCTGTTGTCGTGGATCTCCTTCTGGTTGCTGTACAGCATTGTTCCTAc   >  2:1614947/1‑139 (MQ=255)                                                                                                                         cAGTCTGGCGATCACTTCCTATTTGCTCTCCCTGCGCTGGGCGAGCGATCTCAATACTGTCATCGACAACGTGCTGCGTATTTTTCCGTTGCTGTTGTCGTGGATCTCCTTCTGGTTGCTGTACAGCATTGTTCCTAcc  <  1:1614947/139‑1 (MQ=255)                                                                                                                         cAGTCTGGCGATCACTTCCTATTTGCTCTCCCTGCGCTGGGCGAGCGATCTCAATACTGTCATCGACAACGTGCTGCGTATTTTTCCGTTGCTGTTGTCGTGGATCTCCTTCTGGTTGCTGTACAGCATTGTTCCTAcc  <  1:1206906/139‑1 (MQ=255)                                                                                                                                gCGATCACTTCCTATTTGCTCTCCCTGCGCTGGGCGAGCGATCTCAATACTGTCATCg                                                                            <  1:202139/58‑1 (MQ=255)                                                                                                                                gCGATCACTTCCTATTTGCTCTCCCTGCGCTGGGCGAGCGATCTCAATACTGTCATCg                                                                            >  2:202139/1‑58 (MQ=255)                                                                                                                                 cGATCACTTCCTATTTGCTCTCCCTGCGCTGGGCGAGCGATCTCAATACTGTCATCGACAACGTGCTGCGTATTTTTCCGt                                                    >  1:1409784/1‑81 (MQ=255)                                                                                                                                 cGATCACTTCCTATTTGCTCTCCCTGCGCTGGGCGAGCGATCTCAATACTGTCATCGACAACGTGCTGCGTATTTTTCCGt                                                    <  2:1409784/81‑1 (MQ=255)                                                                                                                                       |                                                                                                                             GATGTACTCCATCGATAGCGCGTTGAATACCATCTGGCGCAGTAAACGAGCGCGACCAAAAATTTACTCGTTCGCCGTGTACTGGATGATTTTAACGCTGGGGCCGCTGCTGGCAGGGGCCAGTCTGGCGATCAGTTCCTATTTGCTCTCCCTGCGCTGGGCGAGCGATCTCAATACTGTCATCGACAACGTGCTGCGTATTTTTCCGTTGCTGTTGTCGTGGATCTCCTTCTGGTTGCTGTACAGCATTGTTCCTACC  >  NC_000913/4076493‑4076751

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