Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,281,472	A→G	E498G (GAG→GGG)	yjcE →	putative cation/proton antiporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,281,472	0	A	G	100.0%	213.0 / NA	56	E498G (GAG→GGG)	yjcE	putative cation/proton antiporter
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (27/29);  total (27/29)

GGAGGAGCGTCTGGCGGCGGATACCGAAGAGAATATCGATAACCAGCTGCTTACGGAGGTCAGTTCTCGCGTCATTGGTAACCTGCGTCGTCGCGCCGATGGACGTAACGACGTTGAAAGTTCCGTGCAGGAAGAGAACCTTGAGCGTCGCTTCCGTCTGGCGGCATTGCGTTCTGAACGTGCTGAACTTTACCACCTGCGCGCCACACGGGAGATCAGCAACGAAACGCTGCAAAAATTACTGCACGATCTCGATTTGCTTGAAGCGT  >  NC_000913/4281338‑4281606                                                                                                                                       |                                                                                                                                       ggaggaGCGTCTGGCGGCGGATACCGAAGAGAATATCGATAACCAGCTGCTTACGGAGGTCAGTTCTCGCGTCATTGGTAACCTGCGTCGTCGCGCCGATGGACGTAACGACGTTGAAAGTTCCGTGCAGGAAGGGAAc                                                                                                                                    >  1:953187/1‑139 (MQ=255)   aggagCGTCTGGCGGCGGATACCGAAGAGAATATCGATAACCAGCTGCTTACGGAGGTCAGTTCTCGCGTCATTGGTAACCTGCGTCGTCGCGCCGATGGACGTAACGACGTTGAAAGTTCCGTGCAGGAAGGGAACCt                                                                                                                                  <  2:676277/139‑1 (MQ=255)      agCGTCTGGCGGCGGATACCGAAGAGAATATCGATAACCAGCTGCTTACGGAGGTCAGTTCTCGCGTCATTGGTAACCTGCGTCGTCGCGCCGATGGACGTAACGACGTTGAAAGTTCCGTGCAGGAAGGGAACCTTGa                                                                                                                               >  2:859042/1‑139 (MQ=255)                    gATACCGAAGAGAATATCGATAACCAGCTGCTTACGGAGGTCAGTTCTCGCGTCATTGGTAACCTGCGTCGTCGCGCCGATGGACGTAACGACGTTGAAAGTTCCGTGCAGGAAGGGAACCTTGAGCGTCGCTTCCGTc                                                                                                                 >  2:496910/1‑139 (MQ=255)                       aCCGAAGAGAATATCGATAACCAGCTGCTTACGGAGGTCAGTTCTCGCGTCATTGGTAACCTGCGTCGTCGCGCCGATGGACGTAACGACGTTGAAAGTTCCGTGCAGGAAGGGAACCTTGAGCGTCGCTTCCGTCTgg                                                                                                              <  2:440413/139‑1 (MQ=255)                                  tatCGATAACCAGCTGCTTACGGAGGTCAGTTCTCGCGTCATTGGTAACCTGCGTCGTCGCGCCGATGGACGTAACGACGTTGAAAGTTCCGTGCAGGAAGGGAACCTTGAGCGTCGCTTCCGTCTggcg                                                                                                            <  1:664093/130‑1 (MQ=255)                                  tatCGATAACCAGCTGCTTACGGAGGTCAGTTCTCGCGTCATTGGTAACCTGCGTCGTCGCGCCGATGGACGTAACGACGTTGAAAGTTCCGTGCAGGAAGGGAACCTTGAGCGTCGCTTCCGTCTggcg                                                                                                            >  2:664093/1‑130 (MQ=255)                                        tAACCAGCTGCTTACGGAGGTCAGTTCTCGCGTCATTGGTAACCTGCATCGTCGCGCCGATGGACGTAACGACGTTGAAAGTTCCGTGCAGGAAGGGAACCTTGAGCGTCGCTTCCGTCTGGCGGCATTGCGTTCTGaa                                                                                             <  2:704110/139‑1 (MQ=255)                                          aCCAGCTGCTTACGGAGGTCAGTTCTCGCGTCATTGGTAACCTGCGTCGTCGCGCCGATGGACGTAACGACGTTGAAAGTTCCGTGCAGGAAGGGAACCTTGAGCGTCGCTTCCGTCTGGCGGCATTGCGTTCTGAACg                                                                                           >  1:1012598/1‑139 (MQ=255)                                            cAGCTGCTTACGGAGGTCAGTTCTCGCGTCATTGGTAACCTGCGTCGTCGCGCCGATGGACGTAACGACGTTGAAAGTTCCGTGCAGGAAGGGAACCTTGAGCGTCGCTTCCGTCTGGCGGCATTGCGTTCTGAACGTg                                                                                         <  1:1078962/139‑1 (MQ=255)                                             aGCTGCTTACGGAGGTCAGTTCTCGCGTCATTGGTAACCTGCGTCGTCGCGCCGATGGACGTAACGAGGTTGAAAGTTCCGTGCAGGAAGGGAACCTTGAGCGTCGCTTCCGTCTGGCGGCATTGCGTTCTGAACGTGc                                                                                        >  2:734945/1‑139 (MQ=255)                                              gctgctTACGGAGGTCAGTTCTCGCGTCATTGGTAACCTGCGTCGTCGCGCCGATGGACGTAACGACGTTGAAAGTTCCGTGCAGGAAGGGAACCTTGAGCGTCGCTTCCGTCTGGCGGCATTGCGTTCTGAACGTGCt                                                                                       >  1:1173387/1‑139 (MQ=255)                                              gctgctTACGGAGGTCAGTTCTCGCGTCATTGGTAACCTGCGTCGTCGCGCCGATGGACGTAACGACGTTGAAAGTTCCGTGCAGGAAGGGAACCTTGAGCGTCGCTTCCGTCTGGCGGCATTGCGTTCTGAACGTGCt                                                                                       >  1:1137552/1‑139 (MQ=255)                                               ctgctTACGGAGGTCAGTTCTCGCGTCATTGGTAACCTGCGTCGTCGCGCCGATGGACGTAACGACGTTGAAAGTTCCGTGCAGGAAGGGAACCTTGAGCGTCGCTTCCGTCTGGCGGCATTGCGTTCTGAACGTGCTg                                                                                      >  2:694476/1‑139 (MQ=255)                                                tgctTACGGAGGTCAGTTCTCGCGTCATTGGTAACCTGCGTCGTCGCGCCGATGGACGTAACGACGTTGAAAGTTCCGTGCAGGAAGGGAACCTTGAGCGTCGCTTCCGTCTGGCGGCATTGCGTTCTGAACGTGCTGa                                                                                     <  2:32363/139‑1 (MQ=255)                                                tgctTACGGAGGTCAGTTCTCGCGTCATTGGTAACCTGCGTCGTCGCGCCGATGGACGTAACGACGTTGAAAGTTCCGTGCAGGAAGGGAACCTTGAGCGTCGCTTCCGTCTGGCGGCATTGCGTTCTGAACGTGCTGa                                                                                     <  2:526427/139‑1 (MQ=255)                                                tgctTACGGAGGTCAGTTCTCGCGTCATTGGTAACCTGCGTCGTCGCGCCGATGGACGTAACGACGTTGAAAGTTCCGTGCAGGAAGGGAACCTTGAGCGTCGCTTCCGTCTGGCGGCATTGCGTTCTGAACGTGCTGa                                                                                     <  1:694476/139‑1 (MQ=255)                                                            tCAGTTCTCGCGTCATTGGTAACCTGCGTCGTCGCGCCGATGGACGTAACGACGTTGAAAGTTCCGTGCAGGAAGGGAACCTTGAGCGTCGCTTCCGTCTGGCGGCa                                                                                                         >  2:682133/1‑107 (MQ=255)                                                            tCAGTTCTCGCGTCATTGGTAACCTGCGTCGTCGCGCCGATGGACGTAACGACGTTGAAAGTTCCGTGCAGGAAGGGAACCTTGAGCGTCGCTTCCGTCTGGCGGCa                                                                                                         <  1:682133/107‑1 (MQ=255)                                                            tCAGTTCTCGCGTCATTGGTAACCTGCGTCGTCGCGCCGATGGACGTAACGACGTTGAAAGTTCCGTGCAGGAAGGGAACCTTGAGCGTCGCTTCCGTCTGGCGGCATTGCGTTCTGAACGTGCt                                                                                       >  2:641813/1‑125 (MQ=255)                                                            tCAGTTCTCGCGTCATTGGTAACCTGCGTCGTCGCGCCGATGGACGTAACGACGTTGAAAGTTCCGTGCAGGAAGGGAACCTTGAGCGTCGCTTCCGTCTGGCGGCATTGCGTTCTGAACGTGCt                                                                                       <  1:641813/125‑1 (MQ=255)                                                            tCAGTTCTCGCGTCATTGGTAACCTGCGTCGTCGCGCCGATGGACGTAACGACGTTGAAAGTTCCGTGCAGGAAGGGAACCTTGAGCGTCGCTTCCGTCTGGCGGCATTGCGTTCTGAACGTGCTGAACTTTACCACCt                                                                         >  1:384028/1‑139 (MQ=255)                                                                  ctcGCGTCATTGGTAACCTGCGTCGTCGCGCCGATGGACGTAACGACCTTGAAAGTTCCGTGCAGGAAGGGAACCTTGAGCGTCGCTTCCGTCTGGCGGCATTGCGTTCTGAACGTGCTGAACTTTACCACCTgcgcgc                                                                   <  2:703603/139‑1 (MQ=255)                                                                    cgcgTCATTGGTAACCTGCGTCGTCGCGCCGATGGACGTAACGACGTTGAAAGTTCCGTGCAGGAAGGGAACCTTGAGCGTCGCTTCCGTCTGGCGGCATTGCGTTCTGAACGTGCTGAACTTTACCACCTGCGCGCca                                                                 <  1:631473/139‑1 (MQ=255)                                                                      cgTCATTGGTAACCTGCGTCGTCGCGCCGATGGACGTAACGACGTTGAAAGTTCCGTGCAGGAAGGGAACCTTGAGCGTCGCTTCCGTCTGGCGGCATTGCGTTCTGAACGTGCTGAACTTTACCACCTGCGCGCcaca                                                               <  2:1173387/139‑1 (MQ=255)                                                                       gTCATTGGTAACCTGCGTCGTCGCGCCGATGGACGTAACGACGTTGAAAGTTCCGTGCAGGAAGGGAACCTTGAGCGTCGCTTCCGTCTGGCGGCATTGCGTTCTGAACGTGCTGAACTTTACCACCTGCGCGCCacac                                                              <  1:734945/139‑1 (MQ=255)                                                                       gTCATTGGTAACCTGCGTCGTCGCGCCGATGGACGTAACGACGTTGAAAGTTCCGTGCAGGAAGGGAACCTTGAGCGTCGCTTCCGTCTGGCGGCATTGCGTTCTGAACGTGCTGAACTTTACCACCTGCGCGCCacac                                                              <  2:1137552/139‑1 (MQ=255)                                                                         cATTGGTAACCTGCGTCGTCGCGCCGATGGACGTAACGACGTTGAAAGTTCCGTGCAGGAAGGGAACCTTGAGCGTCGCTTCCGTCTGGCGGCATTGCGTTCTGAACGTGCTGAACTTTACCACCTGCGCGCCACACgg                                                            <  2:953187/139‑1 (MQ=255)                                                                               tAACCTGCGTCGTCGCGCCGATGGACGTAACGACGTTGAAAGTTCCGTGCAGGAAGGGAACCTTGAGCGTCGCTTCCGTCTGGCGGCATTGCGTTCTGAACGTGCTGAACTTTACCACCTGCGCGCCACACGGGAGATc                                                      <  2:384028/139‑1 (MQ=255)                                                                                      cgtcgtCGCGCCGATGGACGTAACGACGTTGAAAGTTCCGTGCAGGAAGGGAACCTTGAGCGTCGCTTCCGTCTGGCGGCATTGCGTTCTGa                                                                                              <  2:1097003/92‑1 (MQ=255)                                                                                      cgtcgtCGCGCCGATGGACGTAACGACGTTGAAAGTTCCGTGCAGGAAGGGAACCTTGAGCGTCGCTTCCGTCTGGCGGCATTGCGTTCTGa                                                                                              >  1:1097003/1‑92 (MQ=255)                                                                                      cgtcgtCGCGCCGATGGACGTAACGACGTTGAAAGTTCCGTGCAGGAAGGGAACCTTGAGCGTCGCTTCCGTCTGGCGGCATTGCGTTCTGAACGTGCTGAACTTTACCACCTgcgcg                                                                    >  2:81313/1‑118 (MQ=255)                                                                                      cgtcgtCGCGCCGATGGACGTAACGACGTTGAAAGTTCCGTGCAGGAAGGGAACCTTGAGCGTCGCTTCCGTCTGGCGGCATTGCGTTCTGAACGTGCTGAACTTTACCACCTgcgcg                                                                    <  1:81313/118‑1 (MQ=255)                                                                                      cgtcgtCGCGCCGATGGACGTAACGACGTTGAAAGTTCCGTGCAGGAAGGGAACCTTGAGCGTCGCTTCCGTCTGGCGGCATTGCGTTCTGAACGTGCTGAACTTTACCACCTGCGCGCca                                                                 <  1:1007865/121‑1 (MQ=255)                                                                                      cgtcgtCGCGCCGATGGACGTAACGACGTTGAAAGTTCCGTGCAGGAAGGGAACCTTGAGCGTCGCTTCCGTCTGGCGGCATTGCGTTCTGAACGTGCTGAACTTTACCACCTGCGCGCca                                                                 >  2:1007865/1‑121 (MQ=255)                                                                                      cgtcgtCGCGCCGATGGACGTAACGACGTTGAAAGTTCCGTGCAGGAAGGGAACCTTGAGCGTCGCTTCCGTCTGGCGGCATTGCGTTCTGAACGTGCTGAACTTTACCACCTGCGCGCCacac                                                              >  1:288840/1‑124 (MQ=255)                                                                                      cgtcgtCGCGCCGATGGACGTAACGACGTTGAAAGTTCCGTGCAGGAAGGGAACCTTGAGCGTCGCTTCCGTCTGGCGGCATTGCGTTCTGAACGTGCTGAACTTTACCACCTGCGCGCCacac                                                              >  2:525367/1‑124 (MQ=255)                                                                                      cgtcgtCGCGCCGATGGACGTAACGACGTTGAAAGTTCCGTGCAGGAAGGGAACCTTGAGCGTCGCTTCCGTCTGGCGGCATTGCGTTCTGAACGTGCTGAACTTTACCACCTGCGCGCCacac                                                              <  1:525367/124‑1 (MQ=255)                                                                                      cgtcgtCGCGCCGATGGACGTAACGACGTTGAAAGTTCCGTGCAGGAAGGGAACCTTGAGCGTCGCTTCCGTCTGGCGGCATTGCGTTCTGAACGTGCTGAACTTTACCACCTGCGCGCCacac                                                              <  2:288840/124‑1 (MQ=255)                                                                                      cgtcgtCGCGCCGATGGACGTAACGACGTTGAAAGTTCCGTGCAGGAAGGGAACCTTGAGCGTCGCTTCCGTCTGGCGGCATTGCGTTCTGAACGTGCTGAACTTTACCACCTGCGCGCCACACGGGAGATCAGCAACg                                               >  2:336611/1‑139 (MQ=255)                                                                                      cgtcgtCGCGCCGATGGACGTAACGACGTTGAAAGTTCCGTGCAGGAAGGGAACCTTGAGCGTCGCTTCCGTCTGGCGGCATTGCGTTCTGAACGTGCTGAACTTTACCACCTGCGCGCCACACGGGAGATCAGCAACg                                               >  2:69300/1‑139 (MQ=255)                                                                                                ccGATGGACGTAACGACGTTGAAAGTTCCGTGCAGGAAGGGAACCTTGAGCGTCGCTTCCGTCTGGCGGCATTGCGTTCTGAACGTGCTGAACTTTACCACCTGCGCGCCACACGGGAGATCAGCAACGAAACGCTGCa                                     >  1:648131/1‑139 (MQ=255)                                                                                                  gATGGACGTAACGACGTTGAAAGTTCCGTGCAGGAAGGGAACCTTGAGCGTCGCTTCCGTCTGGCGGCATTGCGTTCTGAACGTGCTGAACTTTACCACCTGCGCGCCACACGGGAGATCAGCAACGAAACGCTGCaaa                                   <  2:648131/139‑1 (MQ=255)                                                                                                     ggACGTAACGACGTTGAAAGTTCCGTGCAGGAAGGGAACCTTGAGCGTCGCTTCCGTCTGGCGGCATTGCGTTCTGAACGTGCTGAACTTTACCACCTGCGCGCCACACGGGAGATCAGCAACGAAACGCTGCAAAAAt                                <  1:1152168/139‑1 (MQ=255)                                                                                                            acgacgTTGAAAGTTCCGTGCAGGAAGGGAACCTTGAGCGTCGCTTCCGTCTGGCGGCATTGCGTTCTGAACGTGCTGAACTTTACCACCTGCGCGCCACACGGGAGATCAGCAACGAAACGCTGCAAAAATTACTGCa                         <  1:336611/139‑1 (MQ=255)                                                                                                                cgTTGAAAGTTCCGTGCAGGAAGGGAACCTTGAGCGTCGCTTCCGTCTGGCTGCATTGCGTTCTGAACGTGCTTAACTTTACCACCTGCGCGCCACACGGGAGATCAGCAACGAAACGCTGCAAGAATTACTGCACGAt                     >  2:721650/1‑139 (MQ=255)                                                                                                                   tGAAAGTTCCGCGCAGGAAGGGAACCTTGAGCGTCGCTTCCGTCTGGCGGCATTGCGTTCTGAACGTGCTGAACTTTACCACCTGCGCGCCACACGGGAGATCAGCAACGAAACGCTGCAAAAATTACTGCACGAtctc                  >  1:243457/1‑139 (MQ=255)                                                                                                                     aaaGTTCCGTGCAGGAAGGGAACCTTGAGCGTCGCTTCCGTCTGGCGGCATTGCGTTCTGAACGTGCTGa                                                                                     >  1:255579/1‑70 (MQ=255)                                                                                                                     aaaGTTCCGTGCAGGAAGGGAACCTTGAGCGTCGCTTCCGTCTGGCGGCATTGCGTTCTGAACGTGCTGa                                                                                     <  2:255579/70‑1 (MQ=255)                                                                                                                     aaaGTTCCGTGCAGGAAGGGAACCTTGAGCGTCGCTTCCGTCTGGCGGCATTGCGTTCTGAACGTGCTGa                                                                                     >  1:1083678/1‑70 (MQ=255)                                                                                                                     aaaGTTCCGTGCAGGAAGGGAACCTTGAGCGTCGCTTCCGTCTGGCGGCATTGCGTTCTGAACGTGCTGa                                                                                     <  2:1083678/70‑1 (MQ=255)                                                                                                                            cGTGCAGGAAGGGAACCTTGAGCGTCGCTTCCGTCTGGCGGCATTGCGTTCTGAACGTGCTGAACTTTACCACCTGCGCGCCACACGGGAGATCAGCAACGAAACGCTGCAAAAATTACTGCACGATCTCGATTTGCtt         >  1:182863/1‑139 (MQ=255)                                                                                                                             gTGCAGGAAGGGAACCTTGAGCGTCGCTTCCGTCTGGCGGCATTGCGTTCTGAACGTGCTGAACTTTACCACCTGCGCGCCACACGGGAGATCAGCAACGAAACGCTGCAAAAATTACTGCACGATCTCGATTTGCTTg        <  1:946052/139‑1 (MQ=255)                                                                                                                                   gAAGGGAACCTTGAGCGTCGCTTCCGTCTGGCGGCATTGCGTTCTGAACGTGc                                                                                        <  2:1065254/53‑1 (MQ=255)                                                                                                                                   gAAGGGAACCTTGAGCGTCGCTTCCGTCTGGCGGCATTGCGTTCTGAACGTGc                                                                                        >  1:1065254/1‑53 (MQ=255)                                                                                                                                   gAAGGGAACCTTGAGCGTCGCTTCCGTCTGGCGGCATTGCGTTCTGAACGTGCTGAACTTTACCACCTGCGCGCCACACGGGAGATCAGCAACGAAACGCTGCAAAAATTACTGCACGATCTCGATTTGCTTGAAGCGt  >  1:523676/1‑139 (MQ=255)                                                                                                                                       |                                                                                                                                       GGAGGAGCGTCTGGCGGCGGATACCGAAGAGAATATCGATAACCAGCTGCTTACGGAGGTCAGTTCTCGCGTCATTGGTAACCTGCGTCGTCGCGCCGATGGACGTAACGACGTTGAAAGTTCCGTGCAGGAAGAGAACCTTGAGCGTCGCTTCCGTCTGGCGGCATTGCGTTCTGAACGTGCTGAACTTTACCACCTGCGCGCCACACGGGAGATCAGCAACGAAACGCTGCAAAAATTACTGCACGATCTCGATTTGCTTGAAGCGT  >  NC_000913/4281338‑4281606

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