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%	204.6 / NA	54	E498G (GAG→GGG)	yjcE	putative cation/proton antiporter
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (31/23);  total (31/23)

AGGAGCGTCTGGCGGCGGATACCGAAGAGAATATCGATAACCAGCTGCTTACGGAGGTCAGTTCTCGCGTCATTGGTAACCTGCGTCGTCGCGCCGATGGACGTAACGACGTTGAAAGTTCCGTGCAGGAAGAGAACCTTGAGCGTCGCTTCCGTCTGGCGGCATTGCGTTCTGAACGTGCTGAACTTTACCACCTGCGCGCCACACGGGAGATCAGCAACGAAACGCTGCAAAAATTACTGCACGATCTCGATTTGCT  >  NC_000913/4281340‑4281598                                                                                                                                     |                                                                                                                               aggagCGTCTGGCGGCGGATACCGAAGAGAATATCGATAACCAGCTGCTTACGGAGGTCAGTTCTCGCGTCATTGGTAACCTGCGTCGTCGCGCCGATGGACGTAACGACGTTGAAAGTTCCGTGCAGGAAGGGAACCt                                                                                                                          >  1:672455/1‑139 (MQ=255)     gCGTCTGGCGGCGGATACCGAAGAGAATATCGATAACCAGCTGCTTACGGAGGTCAGTTCTCGCGTCATTGGTAACCTGCGTCGTCGCGCCGATGGACGTAACGACGTTGAAAGTTCCGTGCAGGAAGGGAACCTTGAg                                                                                                                      <  2:1511055/139‑1 (MQ=255)      cGTCTGGCGGCGGATACCGAAGAGAATATCGATAACCAGCTGCTTACGGCGGTCAGTTCTCGCGTCATTGGTAACCTGCGTCGTCGCGCCGATGGACGTAACGACGTTGAAAGTTCCGTGCAGGAAGGGAACCTTGAGc                                                                                                                     >  1:1411846/1‑139 (MQ=255)      cGTCTGGCGGCGGATACCGAAGAGAATATCGATAACCAGCTGCTTACGGAGGTCAGTTCTCGCGTCATTGGTAACCTGCGTCGTCGCGCCGATGGACGTAACGACGTTGAAAGTTCCGTGCAGGAAGGGAACCTTGAGc                                                                                                                     >  2:705284/1‑139 (MQ=255)                          agagAATATCGATAACCAGCTGCTTACGGAGGTCAGTTCTCGCGTCATTGGTAACCTGCGTCGTCGCGCCGATGGACGTAACGACGTTGAAAGTTCCGTGCAGGAAGGGAACCTTGAGCGTCGCTTCCGTCTGGCGGCa                                                                                                 <  1:705284/139‑1 (MQ=255)                             gaATATCGATAACCAGCTGCTTACGGAGGTCAGTTCTCGCGTCATTGGTAACCTGCGTCGTCGCGCCGATGGACGTAACGACGTTGAAAGTTCCGTGCAGGAAGGGAACCTTGAGCGTCGCTTCCGTCTGGCGGCATTg                                                                                              <  2:150270/139‑1 (MQ=255)                                 atCGATAACCAGCTGCTTACGGAGGTCAGTTCTCGCGTCATTGGTAACCTGCGTCGTCGCGCCGATGGACGTAACGACGTTGAAAGTTCCGTGCAGGAAGGGAACCTTGAGCGTCGCTTCCGTCTGGCGGCATTGCGtt                                                                                          >  1:1490560/1‑139 (MQ=255)                                      tAACCAGCTGCTTACGGAGGTCAGTTCTCGCGTCATTGGTAACCTGCGTCGTCGCGCCGATGGACGTAACGACGTTGAAAGTTCCGTGCAGGAAGGGAACCTTGAGCGTCGCTTCCGTCTGGCGGCAtt                                                                                               <  2:490102/129‑1 (MQ=255)                                      tAACCAGCTGCTTACGGAGGTCAGTTCTCGCGTCATTGGTAACCTGCGTCGTCGCGCCGATGGACGTAACGACGTTGAAAGTTCCGTGCAGGAAGGGAACCTTGAGCGTCGCTTCCGTCTGGCGGCAtt                                                                                               >  1:490102/1‑129 (MQ=255)                                        aCCAGCTGCTTACGGAGGTCAGTTCTCGCGTCATTGGTAACCTGCGTCGTCGCGCCGATGGACGTAACGACGTTGAAAGTTCCGTGCAGGAAGGGAACCTTGAGCGTCGCTTCCGTCTGGCGGCATTGCGTTCTGAACg                                                                                   >  1:1362916/1‑139 (MQ=255)                                            gctgctTACGGAGGTCAGTTCTCGCGTCATTGGTAACCTGCGTCGTCGCGCCGATGGACGTAACGACGTTGAAAGTTCCGTGCAGGAAGGGAACCTTGAGCGTCGCTTCCGTCTGGCGGCAt                                                                                                <  1:1201880/122‑1 (MQ=255)                                            gctgctTACGGAGGTCAGTTCTCGCGTCATTGGTAACCTGCGTCGTCGCGCCGATGGACGTAACGACGTTGAAAGTTCCGTGCAGGAAGGGAACCTTGAGCGTCGCTTCCGTCTGGCGGCAt                                                                                                >  2:1201880/1‑122 (MQ=255)                                             ctgctTACGGAGGTCAGTTCTCGCGTCATTGGTAACCTGCGTCGTCGCGCCGATGGACGTAACGACGTTGAAAGTTCCGTGCAGGAAGGGAACCTTGAGCGTCGCTTCCGTCTGGCGGCATTGCGTTCTGAACGTGCTg                                                                              <  2:1490560/139‑1 (MQ=255)                                             ctgctTACGCAGGTCAGTTCTCGCGTCATTGGTAACCTGCGTCGTCGCGCCGATGGACGTAACGACGTTGAAAGTTCCGTGCAGGAAGGGAACCTTGAGCGTCGCTTCCGTCTGGCGGCATTGCGTTCTGAACGTGCTg                                                                              >  2:1106295/1‑139 (MQ=255)                                              tgctTACGGAGGTCAGTTCTCGCGTCATTGGTAACCTGCGTCGTCGCGCCGATGGACGTAACGACGTTGAAAGTTCCGTGCAGGAAGGGAACCTTGAGCGTCGCTTCCGTCTGGCGGCATTGCGTTCTGAACGTGCTGa                                                                             <  1:1031154/139‑1 (MQ=255)                                              tgctTACGGAGGTCAGTTCTCGCGTCATTGGTAACCTGCGTCGTCGCGCCGATGGACGTAACGACGTTGAAAGTTCCGTGCAGGAAGGGAACCTTGAGCGTCGCTTCCGTCTGGCGGCATTGCGTTCTGAACGTGCTGa                                                                             <  1:1297008/139‑1 (MQ=255)                                              tgctTACGGAGGTCAGTTCTCGCGTCATTGGTAACCTGCGTCGTCGCGCCGATGGACGTAACGACGTTGAAAGTTCCGTGCAGGAAGGGAACCTTGAGCGTCGCTTCCGTCTGGCGGCATTGCGTTCTGAACGTGCTGa                                                                             <  2:672455/139‑1 (MQ=255)                                                         gTCAGTTCTCGCGTCATTGGTAACCTGCGTCGTCGCGCCGATGGACGTAACGACGTTGAAAGTTCCGTGCAGGAAGGGAACCTTGAGCGTCGCTTCCGTCTGGCGGCATTGCGTTCTGAACGTGCTGAACTTTaccacc                                                                  >  1:462611/1‑139 (MQ=255)                                                          tCAGTTCTCGCGTCATTGGTAACCTGCGTCGTCGCGCCGATGGACGTAACGACGTTGAAAGTTCCGTGCAGGAAGGGAACCTTGAGCGTCGc                                                                                                                <  2:604759/92‑1 (MQ=255)                                                          tCAGTTCTCGCGTCATTGGTAACCTGCGTCGTCGCGCCGATGGACGTAACGACGTTGAAAGTTCCGTGCAGGAAGGGAACCTTGAGCGTCGc                                                                                                                >  1:604759/1‑92 (MQ=255)                                                          tCAGTTCTCGCGTCATTGGTAACCTGCGTCGTCGCGCCGATGGACGTAACGACGTTGAAAGTTCCGTGCAGGAAGGGAACCTTGAGCGTCGCTTCCGTCTGGCGGCATTGCGTTCTGAACGTGCTGAACTTTACCACCt                                                                 >  1:366312/1‑139 (MQ=255)                                                          tCAGTTCTCGCGTCATTGGTAACCTGCGTCGTCGCGCCGATGGACGTAACGACGTTGAAAGTTCCGTGCAGGAAGGGAACCTTGAGCGTCGCTTCCGTCTGGCGGCATTGCGTTCTGAACGTGCTGAACTTTACCACCt                                                                 >  2:1602859/1‑139 (MQ=255)                                                              ttCTCGCGTCATTGGTAACCTGCGTCGTCGCGCCGATGGACGTAACGACGTTGAAAGTTCCGTGCAGGAAGGGAACCTTGAGCGTCGCTTCCGTCTGGCGGCATTGCGTTCTGaa                                                                                     >  1:393551/1‑115 (MQ=255)                                                              ttCTCGCGTCATTGGTAACCTGCGTCGTCGCGCCGATGGACGTAACGACGTTGAAAGTTCCGTGCAGGAAGGGAACCTTGAGCGTCGCTTCCGTCTGGCGGCATTGCGTTCTGaa                                                                                     <  2:393551/115‑1 (MQ=255)                                                                   gcgTCATTGGTAACCTGCGTCGTCGCGCCGATGGACGTAACGACGTTGAAAGTTCCGTGCAGGAAGGGAACCTTGAGCGTCGCTTCCGTCTGGCGGCATTGCGTTCTGAACGTGCTGAACTTTACCACCTgcgcgc                                                           <  2:853499/136‑1 (MQ=255)                                                                   gcgTCATTGGTAACCTGCGTCGTCGCGCCGATGGACGTAACGACGTTGAAAGTTCCGTGCAGGAAGGGAACCTTGAGCGTCGCTTCCGTCTGGCGGCATTGCGTTCTGAACGTGCTGAACTTTACCACCTgcgcgc                                                           >  1:853499/1‑136 (MQ=255)                                                                     gTCATTGGTAACCTGCGTCGTCGCGCCGATGGACGTAACGACGTTGAAAGTTCCGTGCAGGAAGGGAACCTTGAGCGTCGCTTCCGTCTGGCGGCATTGCGTTCTGAACGTGCTGAACTTTACCACCTGCGCGCCacac                                                      <  1:1106295/139‑1 (MQ=255)                                                                             tAACCTGCGTCGTCGCGCCGATGGACGTAACGACGTTGAAAGTTCCGTGCAGGAAGGGAACCTTGAGCGTCGCTTCCGTCTGGCGGCAt                                                                                                >  1:1389896/1‑89 (MQ=255)                                                                             tAACCTGCGTCGTCGCGCCGATGGACGTAACGACGTTGAAAGTTCCGTGCAGGAAGGGAACCTTGAGCGTCGCTTCCGTCTGGCGGCAt                                                                                                <  2:1389896/89‑1 (MQ=255)                                                                               aCCTGCGTCGTCGCGCCGATGGACGTAACGACGTTGAAAGTTCCGTGCAGGAAGGGAACCTTGAGCGTCGCTTCCGTCTGGCGGCATTGCGTTCTGAACGTGCTGAACTTTACCACCTGCGCGCCACACGGGAGATCAg                                            <  1:1122292/139‑1 (MQ=255)                                                                                   gCGTCGTCGCGCCGATGGACGTAACGACGTTGAAAGTTCCGTGCAGGAAGGGAACCTTGAGCGTCGCTTCCGTCTGGCGGCATTGCGTTCTGAACGTGCTGAACTTTACCACCTGCGCGCCACACGGGAGATCAGCAAc                                        >  1:1411268/1‑139 (MQ=255)                                                                                    cgtcgtCGCGCCGATGGACGTAACGACGTTGAAAGTTCCGTGCAGGAAGGGAACCTTGAGCGTCGCTTCCGTCTGGCGGCATTGCGTTCTGAACGTGCTGAAGTTTACCACCTGCGCGCCACACGGGAGATCAGCAACg                                       >  2:1036045/1‑139 (MQ=255)                                                                                    cgtcgtCGCGCCGATGGACGTAACGACGTTGAAAGTTCCGTGCAGGAAGGGAACCTTGAGCGTCGCTTCCGTCTGGCGGCATTGCGTTCTGAACGTGCTGAACTTTACCACCTgcgcg                                                            <  2:521461/118‑1 (MQ=255)                                                                                    cgtcgtCGCGCCGATGGACGTAACGACGTTGAAAGTTCCGTGCAGGAAGGGAACCTTGAGCGTCGCTTCCGTCTGGCGGCATTGCGTTCTGAACGTGCTGAACTTTACCACCTgcgcg                                                            >  1:521461/1‑118 (MQ=255)                                                                                    cgtcgtCGCGCCGATGGACGTAACGACGTTGAAAGTTCCGTGCAGGAAGGGAACCTTGAGCGTCGCTTCCGTCTGGCGGCATTGCGTTCTGAACGTGCTGAACTTTACCACCTGCGCGCCACACGGGAGATCAGCAACg                                       >  2:1018114/1‑139 (MQ=255)                                                                                    cgtcgtCGCGCCGATGGACGTAACGACGTTGAAAGTTCCGTGCAGGAAGGGAACCTTGAGCGTCGCTTCCGTCTGGCGGCATTGCGTTCTGAACGTGCTGAACTTTACCACCTGCGCGCCACACGGGAGATCAGCAACg                                       >  2:1546939/1‑139 (MQ=255)                                                                                    cgtcgtCGCGCCGATGGACGTAACGACGTTGAAAGTTCCGTGCAGGAAGGGAACCTTGAGCGTCGCTTCCGTCTGGCGGCATTGCGTTCTGAACGTGCTGAACTTTACCACCTGCGCGCCACACGGGAGATCAGCAACg                                       >  1:54605/1‑139 (MQ=255)                                                                                    cgtcgtCGCGCCGATGGACGTAACGACGTTGAAAGTTCCGTGCAGGAAGGGAACCTTGAGCGTCGCTTCCGTCTGGCGGCATTGCGTTCTGAACGTGCTGAACTTTACCACCTGCGCGCCACACGGGAGATCAGCAACg                                       >  1:671882/1‑139 (MQ=255)                                                                                                gATGGACGTAACGACGTTGAAAGTTCCGTGCAGGAAGGGAACCTTGAGCGTCGCTTCCGTCTGGCGGCATTGCGTTCTGAACGTGCTGAACTTTACCACCTGCGCGCCACACGGGAGATCAGCAACGAAACGCTGCaaa                           <  2:1315680/139‑1 (MQ=255)                                                                                                      cGTAACGACGTTGAAAGTTCCGTGCAGGAAGGGAACCTTGAGCGTCGCTTCCGTCTGGCGGCATTGCGTTCTGAACGTGCTGAACTTTACCACCTGCGCGCCACACGGGAGATCAGCAACGAAACGCTGCAAAAATTAc                     >  1:930822/1‑139 (MQ=255)                                                                                                         aaCGACGTTGAAAGTTCCGTGCAGGAAGGGAAGCTTGAGCGt                                                                                                                   >  2:848953/1‑42 (MQ=255)                                                                                                         aaCGACGTTGAAAGTTCCGTGCAGGAAGGGAACCTTGAGCGt                                                                                                                   <  1:848953/42‑1 (MQ=255)                                                                                                           cgacgTTGAAAGTTCCGTGCAGGAAGGGAACCTTGAGCGTCGCTTCCGTCTGGCGGCATTGCGTTCTGAACGTGCTGAACTTTACCACCTGCGCGCCACACGGGAGATCAGCAACGAAACGCTGCAAAAATTACTGCAc                >  2:1258833/1‑139 (MQ=255)                                                                                                              cgTTGAAAGTTCCGTGCAGGAAGGGAACCTTGAGCGTCGCTTCCGTCTGGCGGCATTGCGTTCTGAACGTGCTGAACTTTACCACCTGCGCGCca                                                         >  1:114575/1‑95 (MQ=255)                                                                                                              cgTTGAAAGTTCCGTGCAGGAAGGGAACCTTGAGCGTCGCTTCCGTCTGGCGGCATTGCGTTCTGAACGTGCTGAACTTTACCACCTGCGCGCca                                                         <  2:114575/95‑1 (MQ=255)                                                                                                              cgTTGAAAGTTCCGTGCAGGAAGGGAACCTTGAGCGTCGCTTCCGTCTGGCGGCATTGCGTTCTGAACGTGCTGAACTTTACCACCTGCGCGCCACACGGGAGAt                                               >  2:141361/1‑105 (MQ=255)                                                                                                              cgTTGAAAGTTCCGTGCAGGAAGGGAACCTTGAGCGTCGCTTCCGTCTGGCGGCATTGCGTTCTGAACGTGCTGAACTTTACCACCTGCGCGCCACACGGGAGAt                                               <  1:141361/105‑1 (MQ=255)                                                                                                                 tGAAAGTTCCGTGCAGGAAGGGAACCTTGAGCGTCGCTTCCGTCTGGCGGCATTGCGTTCTGAACGTGCTGAACTTTACCACCTGCGCGCCACACGGGAGATCAGCAACGAAACGCTGCAAAAATTACTGCACGAtctc          >  1:111513/1‑139 (MQ=255)                                                                                                                   aaaGTTCCGTGCAGGAAGGGAACCTTGAGCGTCGCTTCCGTCTGGCGGCAtt                                                                                               <  2:280164/52‑1 (MQ=255)                                                                                                                   aaaGTTCCGTGCAGGAAGGGAACCTTGAGCGTCGCTTCCGTCTGGCGGCAtt                                                                                               >  1:280164/1‑52 (MQ=255)                                                                                                                   aaaGTTCCGTGCAGGAAGGGAACCTTGAGCGTCGCTTCCGTCTGGCGGCATTGCGTTCTGAACGTGCTGa                                                                             <  1:1333080/70‑1 (MQ=255)                                                                                                                   aaaGTTCCGTGCAGGAAGGGAACCTTGAGCGTCGCTTCCGTCTGGCGGCATTGCGTTCTGAACGTGCTGa                                                                             >  2:1333080/1‑70 (MQ=255)                                                                                                                   aaaGTTCCGTGCAGGAAGGGAACCTTGAGCGTCGCTTCCGTCTGGCGGCATTGCGTTCTGAACGTGCTGAACTTTACCACCTGCGCGCCACACGGGAGATCAGCAACGAAACGCTGCAAAAATTACTGCACGATCTCGa        <  2:769983/139‑1 (MQ=255)                                                                                                                         ccGTGCAGGAAGGGAACCTTGAGCGTCGCTTCCGTCTGGCGGCATTGCGTTCTGAACGTGCTGAACTTTACCACCTGCGCGCCACACGGGAGATCAGCAACGAAACGCTGCAAAAATTACTGCACGATCTCGATTTGCt  >  2:1563690/1‑139 (MQ=255)                                                                                                                                     |                                                                                                                               AGGAGCGTCTGGCGGCGGATACCGAAGAGAATATCGATAACCAGCTGCTTACGGAGGTCAGTTCTCGCGTCATTGGTAACCTGCGTCGTCGCGCCGATGGACGTAACGACGTTGAAAGTTCCGTGCAGGAAGAGAACCTTGAGCGTCGCTTCCGTCTGGCGGCATTGCGTTCTGAACGTGCTGAACTTTACCACCTGCGCGCCACACGGGAGATCAGCAACGAAACGCTGCAAAAATTACTGCACGATCTCGATTTGCT  >  NC_000913/4281340‑4281598

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