Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,491,878	G→A	intergenic (+28/+44)	ydcA → / ← hokB	putative periplasmic protein/toxic polypeptide, small

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,491,878	0	G	A	100.0%	112.6 / NA	36	intergenic (+28/+44)	ydcA/hokB	putative periplasmic protein/toxic polypeptide, small
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (18/18);  total (18/18)

AAACCCTGTTCTGGTTCGAAAGGGGGGATCTCACACTGTACGGCAGGCGGCAAATTTGTCTGTAATGATGGTTCTATTAGTGCATCGAAAAAAACATGCACTAACTGAAGTGTAAAAGGGGTGCCATGAGAAAATGAATCTGCTAGTCAAATGCGCGGGGAAAATCCCCGCGCTTGCCCTTACCTGGACGTGCAGGCCATGAGCGCAGCAACCTCCTTATCACCGTCCCGGAACCGCAGTTCGTAGAGCGTTTGTCGGGT  >  NC_000913/1491743‑1492002                                                                                                                                        |                                                                                                                             aaaCCCTGTTCTGGTTCGAAAGGGGGGATCTCACACTGTACGGCAGGCGGCAAATTTGTCTGTAATGATGGTTCTATTAGTGCATCGAAAAAAACATGCACTAACTGAAGTGTAAAAGGGGTGCCATGAGAaaataaat                                                                                                                           >  2:462506/1‑139 (MQ=255)                  gAAAGGGGGGATCTCACACTGTACGGCAGGCGGCAAATTTGTCTGTAATGATGGTTCTATTAGTGCATCGAAAAAAACATGCACTAACTGAAGTGTAAAAGGGGTGCCATGAGAAAATAAATCTGCTAGTCAAATgcgc                                                                                                          <  2:560027/139‑1 (MQ=255)                   aaaGGGGGGATCTCACACTGTACGGCAGGCGGCAAATTTGTCTGTAATGATGGTTCTATTAGTGCATCGAAAAAAACATGCACTAACTGAAGTGTAAAAGGGGTGCCATGAGAAAATAAATCTGCTAGTCAAATGcgcg                                                                                                         <  1:554406/139‑1 (MQ=255)                                cacaCTGTACGGCAGGCGGCAAATTTGTCTGTAATGATGGTTCTATTAGTGCATCGAAAAAAACATGCACTAACTGAAGTGTAAAAGGGGTGCCATGAGAAAATAAATCTGCTAGTCAAATGCGCGGGGAAAATCCCcg                                                                                            <  1:462506/139‑1 (MQ=255)                                    cTGTACGGCAGGCGGCAAATTTGTCTGTAATGATGGTTCTATTAGTGCATCGAAAAAAACATGCACTAACTGAAGTGTAAAAGGGGTGCCATGAGAAAATAAATCTGCTAGTCAAATGCGCGGGGAAAATCCCCGCGCt                                                                                        <  2:559290/139‑1 (MQ=255)                                       tACGGCAGGCGGCAAATTTGTCTGTAATGATGGTTCTATTAGTGCATCGAAAAAAACATGCACTAACTGAAGTGTAAAAGGGGTGCCATGAGAAAATAAATCTGCTAGTCAAATGCGCGGGGAAAATCCCCGCGCTTGc                                                                                     >  2:359651/1‑139 (MQ=255)                                       tACGGCAGGCGGCAAATTTGTCTGTAATGATGGTTCTATTAGTGCATCGAAAAAAACATGCACTAACTGAAGTGTAAAAGGGGTGCCATGAGAAAATAAATCTGCTAGTCAAATGCGCGGGGAAAATCCCCGCGCTTGc                                                                                     >  1:380729/1‑139 (MQ=255)                                       tACGGCAGGCGGCAAATTTGTCTGTAATGATGGTTCTATTAGTGCATCGAAAAAAACATGCACTAACTGAAGTGTAAAAGGGGTGCCATGAGAAAATAAATCTGCTAGTCAAATGCGCGGGGAAAATCCCCGCGCTTGc                                                                                     >  1:487686/1‑139 (MQ=255)                                       tACGGCAGGCGGCAAATTTGTCTGTAATGATGGTTCTATTAGTGCATCGAAAAAAACATGCACTAACTGAAGTGTAAAAGGGGTGCCATGAGAAAATAAATCTGCTAGTCAAATGCGCGGGGAAAATCCCCGCGCTTGc                                                                                     >  1:770468/1‑139 (MQ=255)                                                 ggcAAATTTGTCTGTAATGATGGTTCTATTAGTGCATCGAAAAAAACATGCACTAACTGAAGTGTAAAAGGGGTGCCATGAGAAAATAAATCTGCTAGTCAAATGCGCGGGGAAAATCCCCGCGCTTGCCCTTACCTgg                                                                           >  2:333221/1‑139 (MQ=255)                                                 ggcAAATTTGTCTGTAATGATGGTTCTATTAGTGCATCGAAAAAAACATGCACTAACTGAAGTGTAAAAGGGGTGCCATGAGAAAATAAATCTGCTAGTCAAATGCGCGGGGAAAATCCCCGCGCTTGCCCTTACCTgg                                                                           <  2:23695/139‑1 (MQ=255)                                                                    atgGTTCTATTAGTGCATCGAAAAAAACATGCACTAACTGAAGTGTAAAAGGGGTGCCATGAGAAAATAAATCTGCTAGTCAAATGCGCGGGGAAAATCCCCGCGCTTGCCCTTACCTGGACg                                                                        <  2:742992/123‑1 (MQ=255)                                                                    atgGTTCTATTAGTGCATCGAAAAAAACATGCACTAACTGAAGTGTAAAAGGGGTGCCATGAGAAAATAAATCTGCTAGTCAAATGCGCGGGGAAAATCCCCGCGCTTGCCCTTACCTGGACg                                                                        >  1:742992/1‑123 (MQ=255)                                                                        ttCTATTAGTGCATCGAAAAAAACATGCACTAACTGAAGTGTAAAAGGGGTGCCATGAGAAAATAAATCTGCTAGTCAAATGCGCGGGGAAAATCCCCGCGCTTGCCCTTACCTGGACGTGCAGGCCATGAGCgcagca                                                    <  1:490888/139‑1 (MQ=255)                                                                         tCTATTAGTGCATCGAAAAAAACATGCACTAACTGAAGTGTAAAAGGGGTGCCATGAGAAAATAAATCTGCTAGTCAAATGCGCGGGGAAAATCCCCGCGCTTGCCCTTACCTGGACGTGCAGGCCATGAGCGCAGCaa                                                   >  2:84214/1‑139 (MQ=255)                                                                          cTATTAGTGCATCGAAAAAAACATGCACTAACTGAAGTGTAAAAGGGGTGCCATGAGAAAATAAATCTGCTAGTCAAATGCGCGGGGAAAATCCCCGCGCTTGCCCTTACCTGGACGTGCAGGCCATGAGCGCAGCAAc                                                  <  2:418075/139‑1 (MQ=255)                                                                           tATTAGTGCATCGAAAAAAACATGCACTAACTGAAGTGTAAAAGGGGTGCCATGAGAAAATAAATCTGCTAGTCAAATGCGCGGGGAAAATCCCCGCGCTTGCCCTTACCTGGACGTGCAGGCCATGAGCGCAGCAAcc                                                 <  2:487686/139‑1 (MQ=255)                                                                            aTTAGTGCATCGAAAAAAACATGCACTAACTGAAGTGTAAAAGGGGTGCCATGAGAAAATAAATCTGCTAGTCAAATGCGCGGGGAAAATCCCCGCGCTTGCCCTTACCTGGACGTGCAGGCCATGAGCGCAGCAAcct                                                >  2:108188/1‑139 (MQ=255)                                                                                  gCATCGATAAAAACATGCACTAACAGAAGTGTAAAAGGGGCGCCATGAGAAAATAAATCTGCTAGTCAAATGCGCGGGGAAAATCCCCGCGCTTGCCCTTACCTGGACGTGCAGGCCATGAGCGCAGCAACCTCCTTAt                                          >  2:490888/1‑139 (MQ=255)                                                                                  gCATCGAAAAAAACATGCACTAACTGAAGTGTAAAAGGGGTGCCATGAGAAAATAAATCTGCTAGTCAAATGCGCGGGGAAAATCCCCGCGCTTGCCCTTACCTGGACGTGCAGGCCATGACCGCAGCAAc                                                  <  1:585635/131‑1 (MQ=255)                                                                                  gCATCGAAAAAAACATGCACTAACTGAAGTGTAAAAGGGGTGCCATGAGAAAATAAATCTGCTAGTCAAATGCGCGGGGAAAATCCCCGCGCTTGCCCTTACCTGGACGTGCAGGCCATGACCGCAGCAAc                                                  >  2:585635/1‑131 (MQ=255)                                                                                   cATCGAAAAAAACATGCACTAACTGAAGTGTAAAAGGGGTGCCATGAGAAAATAAATCTGCTAGTCAAATGCGCGGGGAAAATCCCCGCGCTTGCCCTTACCTGGACGTGCAGGCCATGAGCGCAGCAACCTCCTTATc                                         >  1:205310/1‑139 (MQ=255)                                                                                    aTCGAAAAAAACATGCACTAACTGAAGTGTAAAAGGGGTGCCATGAGAAAATAAATCTGCTAGTCAAATGCGCGGGGAAAATCCCCGCGCTTGCCCTTACCTGGACGTGCAGGCCATGAGCGCAGCAACCTCCTTATCa                                        >  2:662937/1‑139 (MQ=255)                                                                                        aaaaaaaCATGCACTAACTGAAGTGTAAAAGGGGTGCCATGAGAAAATAAATCTGCTAGTCAAATGCGCGGGGAAAATCCCCGCGCTTGCCCTTACCTgg                                                                           >  1:603197/1‑100 (MQ=255)                                                                                        aaaaaaaCATGCACTAACTGAAGTGTAAAAGGGGTGCCATGAGAAAATAAATCTGCTAGTCAAATGCGCGGGGAAAATCCCCGCGCTTGCCCTTACCTgg                                                                           <  2:603197/100‑1 (MQ=255)                                                                                                             aGTGTAAAAGGGGTGCCATGAGAAAATAAATCTGCTAGTCAAATGCGCGGGGAAAATCCCCGCGCTTGCCCTTACCTGGACGTGCAGGCCATGAGCGCAGCAACCTCCTTATCACCGTCCCGGAACCGCAGTTCGTaga               >  1:791334/1‑139 (MQ=255)                                                                                                                gtAAAAGGGGTGCCATGAGAAAATAAATCTGCTAGTCAAATGCGCGGGGAAAATCCCCGCGCTTGCCCTTACCTGGACGTGCAGGCCATGAGCGCAGCAACCTCCTTATCACCGTCCCGGAACCGCAGTTCGTAGAGCg            >  1:770487/1‑139 (MQ=255)                                                                                                                  aaaaGGGGTGCCATGAGAAAATAAATCTGCTAGTCAAATGCGCGGGGAAAATCCCCGCGCTTGCCCTTACCTGGACGTGCAGGCCATGAGCGCAGCAACCTCCTTATCACCGTCCCGGAACCGCAGTTCGTAGAGCGtt          >  1:355714/1‑139 (MQ=255)                                                                                                                     aGGGGTGCCATGAGAAAATAAATCTGCTAGTCAAATGCGCGGGGAAAATCCCCGCGCTTGCCCTTACCTGGACGTGCAGGCCATGAGCGCAGCAACCTCCTTATCACCGTCCCGGAACCGCAGTTCGTAGAGCGTTTGt       <  2:355714/139‑1 (MQ=255)                                                                                                                     aGGGGTGCCATGAGAAAATAAATCTGCTAGTCAAATGCGCGGGGAAAATCCCCGCGCTTGCCCTTACCTGGACGTGCAGGCCATGAGCGCAGCAACCTCCTTATCACCGTCCCGGAACCGCAGTTCGTAGAGCGTTTGt       <  2:770468/139‑1 (MQ=255)                                                                                                                     aGGGGTGCCATGAGAAAATAAATCTGCTAGTCAAATGCGCGGGGAAAATCCCCGCGCTTGCCCTTACCTGGACGTGCAGGCCATGAGCGCAGCAACCTCCTTATCACCGTCCCGGAACCGCAGTTCGTAGAGCGTTTGt       <  2:791334/139‑1 (MQ=255)                                                                                                                     aGGGGTGCCATGAGAAAATAAATCTGCTAGTCAAATGCGCGGGGAAAATCCCCGCGCTTGCCCTTACCTGGACGTGCAGGCCATGAGCGCAGCAACCTCCTTATCACCGTCCCGGAACCGCAGTTCGTAGAGCGTTTGt       <  1:35039/139‑1 (MQ=255)                                                                                                                      ggggTGCCATGAGAAAATAAATCTGCTAGTCAAATGCGCGGGGAAAATCCCCGCGCTTGCCCTTACCTGGACGTGCAGGCCATGAGCGCAGCAACCTCCTTATCACCGTCCCGGAACCGCAGTTCGTAGAGCGTTTGTc      >  1:320323/1‑139 (MQ=255)                                                                                                                          tGCCATGAGAAAATAAATCTGCTAGTCAAATGCGCGGGGAAAATCCCCGCGCTTGCCCTTACCTGGACGTGCAGGCCATGAGCGCAGCAACCTCCTTATCACCGTCCCGGAACCGCAGTTCGTAGAGCGTTTGTCGGGt  <  2:205310/139‑1 (MQ=255)                                                                                                                          tGCCATGAGAAAATAAATCTGCTAGTCAAATGCGCGGGGAAAATCCCCGCGCTTGCCCTTACCTGGACGTGCAGGCCATGAGCGCAGCAACCTCCTTATCACCGTCCCGGAACCGCAGTTCGTAGAGCGTTTGTCGGGt  <  1:662937/139‑1 (MQ=255)                                                                                                                           gCCATGAGAAAATAAATCTGCTAGTCAAATGCGCGGGGAAAATCCCCGCGCTTGCCCTTACCTGGACGTg                                                                      >  2:659510/1‑70 (MQ=255)                                                                                                                           gCCATGAGAAAATAAATCTGCTAGTCAAATGCGCGGGGAAAATCCCCGCGCTTGCCCTTACCTGGACGTg                                                                      <  1:659510/70‑1 (MQ=255)                                                                                                                                        |                                                                                                                             AAACCCTGTTCTGGTTCGAAAGGGGGGATCTCACACTGTACGGCAGGCGGCAAATTTGTCTGTAATGATGGTTCTATTAGTGCATCGAAAAAAACATGCACTAACTGAAGTGTAAAAGGGGTGCCATGAGAAAATGAATCTGCTAGTCAAATGCGCGGGGAAAATCCCCGCGCTTGCCCTTACCTGGACGTGCAGGCCATGAGCGCAGCAACCTCCTTATCACCGTCCCGGAACCGCAGTTCGTAGAGCGTTTGTCGGGT  >  NC_000913/1491743‑1492002

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