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

GAGGAGCGTCTGGCGGCGGATACCGAAGAGAATATCGATAACCAGCTGCTTACGGAGGTCAGTTCTCGCGTCATTGGTAACCTGCGTCGTCGCGCCGATGGACGTAACGACGTTGAAAGTTCCGTGCAGGAAGAGAACCTTGAGCGTCGCTTCCGTCTGGCGGCATTGCGTTCTGAACGTGCTGAACTTTACCACCTGCGCGCCACACGGGAGATCAGCAACGAAACGCTGCAAAAATTACTGCACGATCTCGATTTGCTTGAA  >  NC_000913/4281339‑4281602                                                                                                                                      |                                                                                                                                   gaggagCGTCTGGCGGCGGATACCGAAGAGAATATCGATAACCAGCTGCTTACGGAGGTCAGTTCTCGCGTCATTGGTAACCTGCGTCGTCGCGCCGATGGACGTAACGACGTTGAAAGTTCCGTGCAGGAAGGGAAcc                                                                                                                               >  1:208409/1‑139 (MQ=255)     agCGTCTGGCGGCGGATACCGAAGAGAATATCGATAACCAGCTGCTTACGGAGGTCAGTTCTCGCGTCATTGGTAACCTGCGTCGTCGCGCCGATGGACGTAACGACGTTGAAAGTTCCGTGCAGGAAGGGAACCTTGa                                                                                                                           >  1:112955/1‑139 (MQ=255)                   gATACCGAAGAGAATATCGATAACCAGCTGCTTACGGAGGTCAGTTCTCGCGTCATTGGTAACCTGCGTCGTCGCGCCGATGGACGTAACGACGTTGAAAGTTCCGTGCAGGAAGGGAACCTTGAGCGTCGCTTCCGTc                                                                                                             >  2:938322/1‑139 (MQ=255)                      aCCGAAGAGAATATCGATAACCAGCTGCTTACGGAGGTCAGTTCTCGCGTCATTGGTAACCTGCGTCGTCGCGCCGATGGACGTAACGACGTTGAAAGTTCCGTGCAGGAAGGGAACCTTGAGCGTCGCTTCCGTCTgg                                                                                                          <  1:318535/139‑1 (MQ=255)                          aaGAGAATATCGATAACCAGCTGCTTACGGAGGTCAGTTCTCGCGTCATTGGTAACCTGCGTCGTCGCGCCGATGGACGTAACGACGTTGAAAGTTCCGTGCAGGAAGGGAACCTTGAGCGTCGCTTCCGTCTggcggc                                                                                                      <  1:1348535/139‑1 (MQ=255)                             agaATATCGATAACCAGCTGCTTACGGAGGTCAGTTCTCGCGTCATTGGTAACCTGCGTCGTCGCGCCGATGGACGTAACGACGTTGAAAGTTCCGTGCAGGAAGGGAACCTTGAGCGTCGCTTCCGTCTGGCGGCAtt                                                                                                   <  2:112955/139‑1 (MQ=255)                                         aCCAGCTGCTTACGGAGGTCAGTTCTCGCGTCATTGGTAACCTGCGTCGTCGCGCCGATGGACGTAACGACGTTGAAAGTTCCGTGCAGGAAGGGAACCTTGAGCGTCGCTTCCGTCTGGCGGCATTGCGTTCTGAACg                                                                                       >  1:461044/1‑139 (MQ=255)                                           cAGCTGCTTACGGAGGTCAGTTCTCGCGTCATTGGTAACCTGCGTCGTCGCGCCGATGGACGTAACGACGTTGAAAGTTCCGTGCAGGAAGGGAACCTTGAGCGTCGCTTCCGTCTGGCGGCATTGCGTTCTGAACGTg                                                                                     <  1:1078750/139‑1 (MQ=255)                                             gctgctTACGGAGGTCAGTTCTCGCGTCATTGGTAACCTGCGTCGTCGCGCCGATGGACGTAACGACGTTGAAAGTTCCGTGCAGGAAGGGAACCTTg                                                                                                                            >  2:1261259/1‑98 (MQ=255)                                             gctgctTACGGAGGTCAGTTCTCGCGTCATTGGTAACCTGCGTCGTCGCGCCGATGGACGTAACGACGTTGAAAGTTCCGTGCAGGAAGGGAACCTTg                                                                                                                            <  1:1261259/98‑1 (MQ=255)                                             gctgctTACGGAGGTCAGTTCTCGCGTCATTGGTAACCTGCGTCGTCGCGCCGATGGACGTAACGACGTTGAAAGTTCCGTGCAGGAAGGGAACCTTGAGCGTCGCTTCCGTCTGGCGGCATTGCGTTCTGAACGTGCt                                                                                   >  2:992304/1‑139 (MQ=255)                                             gctgctTACGGAGGTCAGTTCTCGCGTCATTGGTAACCTGCGTCGTCGCGCCGATGGACGTAACGACGTTGAAAGTTCCGTGCAGGAAGGGAACCTTGAGCGTCGCTTCCGTCTGGCGGCATTGCGTTCTGAACGTGCt                                                                                   >  1:860659/1‑139 (MQ=255)                                              ctgctTACGGAGGTCAGTTCTCGCGTCATTGGTAACCTGCGTCGTCGCGCCGATGGACGTAACGACGTTGAAAGTTCCGTGCAGGAAGGGAACCTTGAGCGTCGCTTCCGTCTGGCGGCATTGCGTTCTGAACGTGCTg                                                                                  >  1:1259535/1‑139 (MQ=255)                                              ctgctTACGGAGGTCAGTTCTCGCGTCATTGGTAACCTGCGTCGTCGCGCCGATGGACGTAACGACGTTGAAAGTTCCGTGCAGGAAGGGAACCTTGAGCGTCGCTTCCGTCTGGCGGCATTGCGTTCTGAACGTGCTg                                                                                  >  2:149335/1‑139 (MQ=255)                                               tgctTACGGAGGTCAGTTCTCGCGTCATTGGTAACCTGCGTCGTCGCGCCGATGGACGTAACGACGTTGAAAGTTCCGTGCAGGAAGGGAACCTTGAGCGTCGCTTCCGTCTGGCGGCATTGCGTTCTGAACGTGCTGa                                                                                 <  1:992304/139‑1 (MQ=255)                                                             aGTTCTCGCGTCATTGGTAACCTGCGTCGTCGCGCCGATGGACGTAACGACGTTGAAAGTTCCGTGCAGGAAGGGAACCTTGAGCGTCGCTTCCGTCTGGCGGCATTGCGTTCTGAACGTGCTGAACTTTACCACCTgc                                                                   >  1:666212/1‑139 (MQ=255)                                                               ttCTCGCGTCATTGGTAACCTGCGTCGTCGCGCCGATGGACGTAACGACGTTGAAAGTTCCGTGCAGGAAGGGAACCTTGAGCGTCGCTTCCg                                                                                                               <  1:1367284/93‑1 (MQ=255)                                                               ttCTCGCGTCATTGGTAACCTGCGTCGTCGCGCCGATGGACGTAACGACGTTGAAAGTTCCGTGCAGGAAGGGAACCTTGAGCGTCGCTTCCg                                                                                                               >  2:1367284/1‑93 (MQ=255)                                                                tctcGCGTCATTGGTAACCTGCGTCGTCGCGCCGATGGACGTAACGACGTTGAAAGTTCCGTGCAGGAAGGGAACCTTGAGCGTCGCTTCCGTCTGGCGGCATTGCGTTCTGAACGTGCTGAACTTTACCACCTgcgcg                                                                <  1:1172459/139‑1 (MQ=255)                                                                tctcGCGTCATTGGTAACCTGCGTCGTCGCGCCGATGGACGTAACGACGTTGAAAGTTCCGTGCAGGAAGGGAACCTTGAGCGTCGCTTCCGTCTGGCGGCATTGCGTTCTGAACGTGCTGAACTTTACCACCTgcgcg                                                                <  2:461044/139‑1 (MQ=255)                                                                   cgcgTCATTGGTAACCTGCGTCGTCGCGCCGATGGACGTAACGACGTTGAAAGTTCCGTGCAGGAAGGGAACCTTGAGCGTCGCTTCCGTCTGGCGGCATTGCGTTCTGAACGTGCTGAACTTTACCACCTGCGCGCca                                                             <  2:554708/139‑1 (MQ=255)                                                                    gcgTCATTGGTAACCTGCGTCGTCGCGCCGATGGACGTAACGACGTTGAAAGTTCCGTGCAGGAAGGGAACCTTGAGCGTCGCTTCCGTCTGGCGGCATTGCGTTCTGAACGTGCTGAACTTTACCACCTGCGCGCcac                                                            >  2:950012/1‑139 (MQ=255)                                                                    gcgTCATTGGTAACCTGCGTCGTCGCGCCGATGGACGTAACGACGTTGAAAGTTCCGTGCAGGAAGGGAACCTTGAGCGTCGCTTCCGTCTGGCGGCATTGCGTTCTGAACGTGCTGAACTTTACCACCTGCGCGCcac                                                            >  2:1372122/1‑139 (MQ=255)                                                                    gcgTCATTGGTAACCTGCGTCGTCGCGCCGATGGACGTAACGACGTTGAAAGTTCCGTGCAGGAAGGGAACCTTGAGCGTCGCTTCCGTCTGGCGGCATTGCGTTCTGAACGTGCTGAACTTTACCACCTGCGCGCcac                                                            >  2:116867/1‑139 (MQ=255)                                                                      gTCATTGGTAACCTGCGTCGTCGCGCCGATGGACGTAACGACGTTGAAAGTTCCGTGCAGGAAGGGAACCTTGAGCGTCGCTTCCGTCTGGCGGCATTGCGTTCTGAACGTGCTGAACTTTACCACCTGCGCGCCacac                                                          <  1:1083683/139‑1 (MQ=255)                                                                         aTTGGTAACCTGCGTCGTCGCGCCGATGGACGTAACGACGTTGAAAGTTCCGTGCAGGAAGGGAACCTTGAGCGTCGCTTCCg                                                                                                               >  2:548321/1‑83 (MQ=255)                                                                         aTTGGTAACCTGCGTCGTCGCGCCGATGGACGTAACGACGTTGAAAGTTCCGTGCAGGAAGGGAACCTTGAGCGTCGCTTCCg                                                                                                               <  1:548321/83‑1 (MQ=255)                                                                                     cgtcgtCGCGCCGATGGACGTAACGACGTTGAAAGTTCCGTGCAGGAAGGGATCCTTGAGCGTCGCTTCCGTCTGGCGGCATTGCGTTCTGAACGTGCTg                                                                                  >  1:1017063/1‑100 (MQ=255)                                                                                     cgtcgtCGCGCCGATGGACGTAACGACGTTGAAAGTTCCGTGCAGGAAGGGAACCTTGAGCGTCGCTTCCGTCTGGCGGCAt                                                                                                    <  1:14407/82‑1 (MQ=255)                                                                                     cgtcgtCGCGCCGATGGACGTAACGACGTTGAAAGTTCCGTGCAGGAAGGGAACCTTGAGCGTCGCTTCCGTCTGGCGGCAt                                                                                                    >  2:14407/1‑82 (MQ=255)                                                                                     cgtcgtCGCGCCGATGGACGTAACGACGTTGAAAGTTCCGTGCAGGAAGGGAACCTTGAGCGTCGCTTCCGTCTGGCGGCATTGCGTTCTGAACGTGCTg                                                                                  <  2:1017063/100‑1 (MQ=255)                                                                                     cgtcgtCGCGCCGATGGACGTAACGACGTTGAAAGTTCCGTGCAGGAAGGGAACCTTGAGCGTCGCTTCCGTCTGGCGGCATTGCGTTCTGAACGTGCTGAACTTTACCACCTgcgcg                                                                <  1:535483/118‑1 (MQ=255)                                                                                     cgtcgtCGCGCCGATGGACGTAACGACGTTGAAAGTTCCGTGCAGGAAGGGAACCTTGAGCGTCGCTTCCGTCTGGCGGCATTGCGTTCTGAACGTGCTGAACTTTACCACCTgcgcg                                                                >  2:535483/1‑118 (MQ=255)                                                                                     cgtcgtCGCGCCGATGGACGTAACGACGTTGAAAGTTCCGTGCAGGAAGGGAACCTTGAGCGTCGCTTCCGTCTGGCGGCATTGCGTTCTGAACGTGCTGAACTTTACCACCTGCGCGCCACACGGGAGATCAGCAACg                                           >  1:776863/1‑139 (MQ=255)                                                                                     cgtcgtCGCGCCGATGGACGTAACGACGTTGAAAGTTCCGTGCAGGAAGGGAACCTTGAGCGTCGCTTCCGTCTGGCGGCATTGCGTTCTGAACGTGCTGAACTTTACCACCTGCGCGCCACACGGGAGATCAGCAACg                                           >  2:178671/1‑139 (MQ=255)                                                                                                       cGTAACGACGTTGAAAGTTCCGTGCAGGAAGGGAACCTTGAGCGTCGCTTCCGTCTGGCGGCATTGCGTTCTGAACGTGCTGAACTTTACCACCTGCGCGCca                                                             <  1:1140032/103‑1 (MQ=255)                                                                                                       cGTAACGACGTTGAAAGTTCCGTGCAGGAAGGGAACCTTGAGCGTCGCTTCCGTCTGGCGGCATTGCGTTCTGAACGTGCTGAACTTTACCACCTGCGCGCca                                                             >  2:1140032/1‑103 (MQ=255)                                                                                                                gTTGAAAGTTCCGTGCAGGAAGGGAACCTTGAGCGTCGCTTCCGTCTGGCGGCATTGCGTTCTGAACGTGCTGAACTTTACCACCTGCGCGCCACACGGGAGATCAGCAACGAAACGCTGCAAAAATTACTGCACGAtc                >  1:680772/1‑139 (MQ=255)                                                                                                                 ttGAAAGTTCCGTGCAGGAAGGGAACCTTGAGCGTCGCTTCCGTCTGGCGGCATTGCGTTCTGAACGTGCTGAACTTTACCACCTGCGCGCCACATGGGAGATCAGCAACGAAACGCTGCAAAAATTACTGCACGAtct               >  2:215687/1‑139 (MQ=255)                                                                                                                    aaaGTTCCGTGCAGGAAGGGAACCTTGAGCGTCGCTTCCGTCTGGCGGCATTGCGTTCTGAACGTGCTGAACTTTACCACCTGCGCGCCACACGGGAGATCAGCAACGAAACGCTGCAAAAATTACTGCACGATCTCGa            <  2:666212/139‑1 (MQ=255)                                                                                                                      aGTTCCGTGCAGGAAGGGAACCTTGAGCGTCGCTTCCGTCTGGCGGCATTGCGTTCTGAACGTGCTGAACTTTACCACCTGCGCGCCACACGGGAGATCAGCAACGAAACGCTGCAAAAATTACTGCACGATCTCGAtt          <  1:116867/139‑1 (MQ=255)                                                                                                                         tCCGTGCAGGAAGGGAACCTTGAGCGTCGCTTCCGTCTGGCGGCATTGCGTTCTGAACGTGCTGAACTTTACCACCTGCGCGCCACACGGGAGATCAGCAACGAAACGCTGCAAAAATTACTGCACGATCTCGATTTGc       >  2:1505422/1‑139 (MQ=255)                                                                                                                         tCCGTGCAGGAAGGGAACCTTGAGCGTCGCTTCCGTCTGGCGGCATTGCGTTCTGAACGTGCTGAACTTTACCACCTGCGCGCCACACGGGAGATCAGCAACGAAACGCTGCAAAAATTACTGCACGATCTCGATTTGc       >  1:108718/1‑139 (MQ=255)                                                                                                                             tGCAGGAAGGGAACCTTGAGCGTCGCTTCCGTCTg                                                                                                           <  2:319047/35‑1 (MQ=255)                                                                                                                             tGCAGGAAGGGAACCTTGAGCGTCGCTTCCGTCTg                                                                                                           >  1:319047/1‑35 (MQ=255)                                                                                                                             tGCAGGAAGGGAACCTTGAGCGTCGCTTCCGTCTGGCGGCATTGCGTTCTGAACGTGCTGAACTTTACCACCTGCGCGCCACACGGGAGATCAGCAACGAAACGCTGCAAAAATTACTGCACGATCTCGATTTGCTTGa   >  1:383623/1‑139 (MQ=255)                                                                                                                              gCAGGAAGGGAACCTTGAGCGTCGCTTCCGTCTGGCGGCATTGCGTTCTGAACGTGCTGAACTTTACCACCTGCGCGCCACACGGGAGAt                                                   >  2:276986/1‑90 (MQ=255)                                                                                                                              gCAGGAAGGGAACCTTGAGCGTCGCTTCCGTCTGGCGGCATTGCGTTCTGAACGTGCTGAACTTTACCACCTGCGCGCCACACGGGAGAt                                                   <  1:276986/90‑1 (MQ=255)                                                                                                                              gCAGGAAGGGAACCTTGAGCGTCGCTTCCGTCTGGCGGCATTGCGTTCTGAACGTGCTGAACTTTACCACCTGCGCGCCACACGGGAGATCAGCAACGAAACGCTGCAAAAATTACTGCACGAGCTCGATTTGCTTGaa  <  2:680772/139‑1 (MQ=255)                                                                                                                                      |                                                                                                                                   GAGGAGCGTCTGGCGGCGGATACCGAAGAGAATATCGATAACCAGCTGCTTACGGAGGTCAGTTCTCGCGTCATTGGTAACCTGCGTCGTCGCGCCGATGGACGTAACGACGTTGAAAGTTCCGTGCAGGAAGAGAACCTTGAGCGTCGCTTCCGTCTGGCGGCATTGCGTTCTGAACGTGCTGAACTTTACCACCTGCGCGCCACACGGGAGATCAGCAACGAAACGCTGCAAAAATTACTGCACGATCTCGATTTGCTTGAA  >  NC_000913/4281339‑4281602

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