Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,952,661	G→A	*132* (TAG→TAA)	yecN →	MAPEG family inner membrane protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,952,661	0	G	A	100.0%	169.7 / NA	53	*132* (TAG→TAA)	yecN	MAPEG family inner membrane protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (28/25);  total (28/25)

TGCATTATTACGGTTTTCATCACCGTCTGTTCCGCTGGCGACGTTCCGGCATGAGCGCCACCTGGTGTGCGCTGTTGCTGATGGTGCTGGCGAATCTTTGGTATATGCCCTGGGAGTTGGTTTTCTCCCTGCGTTAGCGCACAATACGCCACTTTCTTTTTCCCGGATTTTTACGTTATGTCTCACCGCGACACGCTATTTTCTGCCCCTATCGCCAGACTGGGCGACTGGACCTTTGATGAACGGGTAGCTGAAGTCTTCCCGGATATG  >  NC_000913/1952525‑1952794                                                                                                                                         |                                                                                                                                      tGCATTATTACGGTTTTCATCACCGTCTGTTCCGCTGGCGACGTTCCGGCATGAGCGCCACCTGGTGTGCGCTGTTGCTGATGGTGCTGGCGAATCTTTGGTATATGCCCTGGGAGTTGGTTTTCTCCCTGCGTTAacg                                                                                                                                     <  1:951847/139‑1 (MQ=255)      tattaCGGTTTTCATCACCGTCTGTTCCGCTGGCGACGTTCCGGCATGAGCGCCACCTGGTGTGCGCTGTTGCTGATGGTGCTGGCGAATCTTTGGTATATGCCCTGGGAGTTGGTTTTCTCCCTGCGTTAACGCACaa                                                                                                                                >  1:1132817/1‑139 (MQ=255)                ttCATCACCGTCTGTTCCGCTGGCGACGTTCCGGCATGAGCGCCACCTGGTGTGCGCTGTTGCTGATGGTGCTGGCGAATCTTTGGTATATGCCCTGGGAGTTGGTTTTCTCCCTGCGTTAACGCACAATACGCCActt                                                                                                                      >  1:1204927/1‑139 (MQ=255)                ttCATCACCGTCTGTTCCGCTGGCGACGTTCCGGCATGAGCGCCACCTGGTGTGCGCTGTTGCTGATGGTGCTGGCGAATCTTTGGTATATGCCCTGGGAGTTGGTTTTCTCCCTGCGTTAACGCACAATACGCCActt                                                                                                                      >  2:862016/1‑139 (MQ=255)                ttCATCACCGTCTGTTCCGCTGGCGACGTTCCGGCATGAGCGCCACCTGGTGTGCGCTGTTGCTGATGGTGCTGGCGAATCTTTGGTATATGCCCTGGGAGTTGGTTTTCTCCCTGCGTTAACGCACAATACGCCActt                                                                                                                      >  1:1432802/1‑139 (MQ=255)                       ccGTCTGTTCCGCTGGCGACGTTCCGGCATGAGCGCCACCTGGTGTGCGCTGTTGCTGATGGTGCTGGCGAATCTTTGGTATATGCCCTGGGAGTTGGTTTTCTCCCTGCGTTAACGCACAATACGCCACTTTCttttt                                                                                                               >  2:1139105/1‑139 (MQ=255)                          tCTGTTCCGCTGGCGACGTTCCGGCATGAGCGCCACCTGGTGTGCGCTGTTGCTGATGGTGCTGGCGAATCTTTGGTATATGCCCTGGGAGTTGGTTTTCTCCCTGCGTTAacg                                                                                                                                     <  2:56918/114‑1 (MQ=255)                          tCTGTTCCGCTGGCGACGTTCCGGCATGAGCGCCACCTGGTGTGCGCTGTTGCTGATGGTGCTGGCGAATCTTTGGTATATGCCCTGGGAGTTGGTTTTCTCCCTGCGTTAacg                                                                                                                                     >  1:56918/1‑114 (MQ=255)                          tCTGTTCCGCTGGCGACGTTCCGGCATGAGCGCCACCTGGTGTGCGCTGTTGCTGATGGTGCTGGCGAATCTTTGGTATATGCCCTGGGAGTTGGTTTTCTCCCTGCGTTAACGCACAATACGCCa                                                                                                                         <  1:72551/126‑1 (MQ=255)                          tCTGTTCCGCTGGCGACGTTCCGGCATGAGCGCCACCTGGTGTGCGCTGTTGCTGATGGTGCTGGCGAATCTTTGGTATATGCCCTGGGAGTTGGTTTTCTCCCTGCGTTAACGCACAATACGCCa                                                                                                                         >  2:72551/1‑126 (MQ=255)                          tCTGTTCCGCTGGCGACGTTCCGGCATGAGCGCCACCTGGTGTGCGCTGTTGCTGATGGTGCTGGCGAATCTTTGGTATATGCCCTGGGAGTTGGTTTTCTCCCTGCGTTAACGCACAATACGCCACTTTCTTTTTccc                                                                                                            >  1:1342411/1‑139 (MQ=255)                                   cTGGCGACGTTCCGGCATGAGCGCCATCTGGTGTGCGCTGTTGCTGATGGTGCTGGCGAATCTTTGGTATATGCCCTGGGAGTTGGTTTTCTCCCTGCGTTAACGCACAATACGCCActttct                                                                                                                   >  2:1263822/1‑123 (MQ=255)                                   cTGGCGACGTTCCGGCATGAGCGCCATCTGGTGTGCGCTGTTGCTGATGGTGCTGGCGAATCTTTGGTATATGCCCTGGGAGTTGGTTTTCTCCCTGCGTTAACGCACAATACGCCActttct                                                                                                                   <  1:1263822/123‑1 (MQ=255)                                   cTGGCGACGTTCCGGCATGAGCGCCACCTGGTGTGCGCTGTTGCTGATGGTGCTGGCGAATCTTTGGTATATGCCCTGGGAGTTGGTTTTCTCCCTGCGTTAACGCACAATACGCCACTTTCTTTTTCCCGGATTTTTa                                                                                                   >  1:1144304/1‑139 (MQ=255)                                      gCGACGTTCCGGCATGAGCGCCACCTGGTGTGCGCTGTTGCTGATGGTGCTGGCGAATCTTTGGTATATGCCCTGGGAGTTGGTTTTCTCCCTGCGTTAACGCACAATACGCCACTTTCTTTTTCCCGGAtt                                                                                                       <  2:601255/132‑1 (MQ=255)                                      gCGACGTTCCGGCATGAGCGCCACCTGGTGTGCGCTGTTGCTGATGGTGCTGGCGAATCTTTGGTATATGCCCTGGGAGTTGGTTTTCTCCCTGCGTTAACGCACAATACGCCACTTTCTTTTTCCCGGAtt                                                                                                       >  1:601255/1‑132 (MQ=255)                                       cGACGTTCCGGCATGAGCGCCACCTGGTGTGCGCTGTTGCTGATGGTGCTGGCGAATCTTTGGTATATGCCCTGGGAGTTGGTTTTCTCCCTGCGTTAACGCACAATACGCCACTTTCTTTTTCCCGGATTTTTACGtt                                                                                               <  2:1204927/139‑1 (MQ=255)                                        gACGTTCCGGCATGAGCGCCACCTGGTGTGCGCTGTTGCTGATGGTGCTGGCGAATCTTTGGTATATGCCCTGGGAGTTGGTTTTCTCCCTGCGTTAACGCACAATACGCCACTTTCTTTTTCCCGGATTTTTACGTTa                                                                                              >  2:717477/1‑139 (MQ=255)                                             tCCGGCATGAGCGCCACCTGGTGTGCGCTGTTGCTGATGGTGCTGGCGAATCTTTGGTATATGCCCTGGGAGTTGGTTTTCTCCCTGCGTTAACGCACAATACGCCACTTTCTTTTTCCCGGATTTTTACGTTATGtct                                                                                         <  2:1342411/139‑1 (MQ=255)                                                      aGCGCCACCTGGTGTGCGCTGTTGCTGATGGTGCTGGCGAATCTTTGGTATATGCCCTGGGAGTTGGTTTTCTCCCTGCGTTAACGCACAATACGCCACTTTCTTTTTCCCGGATTTTTACGTTATGTCTCACCGCGac                                                                                >  2:562527/1‑139 (MQ=255)                                                         gcCACCTGGTGTGCGCTGTTGCTGATGGTGCTGGCGAATCTTTGGTATATGCCCTGGGAGTTGGTTTTCTCCCTGCGTTAACGCACAATACGCCACTTTCTTTTTCCCGGATTTTTACGTTATGTCTCACCGCGACACg                                                                             <  1:562527/139‑1 (MQ=255)                                                          ccACCTGGTGTGCGCTGTTGCTGATGGTGCTGGCGAATCTTTGGTATATGCCCTGGGAGTTGGTTTTCTCCCTGCGTTAACGCACAATACGCCACTTTCTTTTTCCCGGATTTTTACGTTATGTCTCACCGCGACACGc                                                                            <  2:702908/139‑1 (MQ=255)                                                                       gcTGTTGCTGATGGTGCTGGCGAATCTTTGGTATATGCCCTGGGAGTTGGTTTTCTCCCTGCGTTAACGCACAATACGCCACTTTCTTTTTCCCGGATTTTTACGTTATGTCTCACCGCGACACGCTATTTTCTGcccc                                                               <  2:696987/139‑1 (MQ=255)                                                                       gcTGTTGCTGATGGTGCTGGCGAATCTTTGGTATATGCCCTGGGAGTTGGTTTTCTCCCTGCGTTAACGCACAATACGCCACTTTCTTTTTCCCGGATTTTTACGTTATGTCTCACCGCGACACGCTATTTTCTGcccc                                                               <  2:1144304/139‑1 (MQ=255)                                                                                     tGCTGGCGAATCTTTGGTATATGCCCTGGGAGTTGGTTTTCTCCCTGCGTTAACGCACAATACGCCACTTTCTTTTTCCCGGATTTTTACGTTATGTCTCACCGCGa                                                                                 >  1:813652/1‑107 (MQ=255)                                                                                     tGCTGGCGAATCTTTGGTATATGCCCTGGGAGTTGGTTTTCTCCCTGCGTTAACGCACAATACGCCACTTTCTTTTTCCCGGATTTTTACGTTATGTCTCACCGCGa                                                                                 <  2:813652/107‑1 (MQ=255)                                                                                     tGCTGGCGAATCTTTGGTATATGCCCTGGGAGTTGGTTTTCTCCCTGCGTTAACGCACAATACGCCACTTTCTTTTTCCCGGATTTTTACGTTATGTCTCACCGCGACACGc                                                                            >  1:1423680/1‑112 (MQ=255)                                                                                     tGCTGGCGAATCTTTGGTATATGCCCTGGGAGTTGGTTTTCTCCCTGCGTTAACGCACAATACGCCACTTTCTTTTTCCCGGATTTTTACGTTATGTCTCACCGCGACACGc                                                                            <  2:1423680/112‑1 (MQ=255)                                                                                       cTGGCGAATCTTTGGTATATGCCCTGGGAGTTGGTTTTCTCCCTGCGTTAACGCACAATACGCCACTTTCTTTTTCCCGGATTTTTACGTTATGTCTCACCGCGa                                                                                 >  2:1514792/1‑105 (MQ=255)                                                                                       cTGGCGAATCTTTGGTATATGCCCTGGGAGTTGGTTTTCTCCCTGCGTTAACGCACAATACGCCACTTTCTTTTTCCCGGATTTTTACGTTATGTCTCACCGCGa                                                                                 <  1:1514792/105‑1 (MQ=255)                                                                                        tGGCGAATCTTTGGTATATGCCCTGGGAGTTGGTTTTCTCCCTGCGTTAACGCACAATACGCCACTTTCTTTTTCCCGGAttt                                                                                                      >  1:1380695/1‑83 (MQ=255)                                                                                        tGGCGAATCTTTGGTATATGCCCTGGGAGTTGGTTTTCTCCCTGCGTTAACGCACAATACGCCACTTTCTTTTTCCCGGAttt                                                                                                      <  2:1380695/83‑1 (MQ=255)                                                                                        tGGCGAATCTTTGGTATATGCCCTGGGAGTTGGTTTTCTCCCTGCGTTAACGCACAATACGCCACTTTCTTTTTCCCGGATTTTTACGTTATGTCTCACCGCGACACGCTATTTTCTGCCCCTATCGCCAGACTGGGCg                                              >  1:851151/1‑139 (MQ=255)                                                                                                         atGCCCTGGGAGTTGGTTTTCTCCCTGCGTTAACGCACAATACGCCCCTTTCTTTTTCCCGGATTTTTACGTTATGTCTCACCGCGACACGCTATTTTCTGCCCCTATCGCCTGACTGGGCGACTGGACCTTtgatgag                             <  2:1132817/139‑2 (MQ=255)                                                                                                          tGCCCTGGGAGTTGGTTTTCTCCCTGCGTTAACGCACAATACGCCACTTTCTTTTTCCCGGATTTTTACGTTATGTCTCACCGCGACACGCTATTTTCTGCCCCTATCGCCAGACTGGGCGACTg                                          >  2:788836/1‑125 (MQ=255)                                                                                                          tGCCCTGGGAGTTGGTTTTCTCCCTGCGTTAACGCACAATACGCCACTTTCTTTTTCCCGGATTTTTACGTTATGTCTCACCGCGACACGCTATTTTCTGCCCCTATCGCCAGACTGGGCGACTg                                          <  1:788836/125‑1 (MQ=255)                                                                                                          tGCCCTGGGAGTTGGTTTTCTCCCTGCGTTAACGCACAATACGCCACTTTCTTTTTCCCGGATTTTTACGTTATGTCTCACCGCGACACGCTATTTTCTGCCCCTATCGCCAGACTGGGCGACTGGACCTTTGATGAAc                            >  1:1991/1‑139 (MQ=255)                                                                                                           gCCCTGGGAGTTGGTTTTCTCCCTGCGTTAACGCACAATACGCCACTTTCTTTTTCCCGGATTTTTACGTTATGTCTCACCGCGACACGCTATTTTCTGCCCCTATCGCCAGACTGGGCGACTGGACCTTTGATGAACg                           <  1:585854/139‑1 (MQ=255)                                                                                                           gCCCTGGGAGTTGGTTTTCTCCCTGCGTTAACGCACAATACGCCACTTTCTTTTTCCCGGATTTTTACGTTATGTCTCACCGCGACACGCTATTTTCTGCCCCTATCGCCAGACTGGGCGACTGGACCTTTGATGAACg                           >  1:143032/1‑139 (MQ=255)                                                                                                             ccTGGGAGTTGGTTTTCTCCCTGCGTTAACGCACAATACGCCACTTTCTTTTTCCCGGATTTTTACGTTATGTCTCACCGCGACACGCt                                                                           >  2:1234067/1‑89 (MQ=255)                                                                                                             ccTGGGAGTTGGTTTTCTCCCTGCGTTAACGCACAATACGCCACTTTCTTTTTCCCGGATTTTTACGTTATGTCTCACCGCGACACGCt                                                                           <  1:1234067/89‑1 (MQ=255)                                                                                                               tGGGAGTTGGTTTTCTCCCTGCGTTAACGCACAATACGCCACTTTCTTTTTCCCGGATTTTTACGTTATGTCTCACCGCGACACGCTATTTTCTGcccc                                                               <  2:702133/99‑1 (MQ=255)                                                                                                               tGGGAGTTGGTTTTCTCCCTGCGTTAACGCACAATACGCCACTTTCTTTTTCCCGGATTTTTACGTTATGTCTCACCGCGACACGCTATTTTCTGcccc                                                               >  1:702133/1‑99 (MQ=255)                                                                                                                gggAGTTGGTTTTCTCCCTGCGTTAACGCACAATACGCCACTTTCTTTTTCCCGGATTTTTACg                                                                                                 <  1:540958/64‑1 (MQ=255)                                                                                                                gggAGTTGGTTTTCTCCCTGCGTTAACGCACAATACGCCACTTTCTTTTTCCCGGATTTTTACg                                                                                                 >  2:540958/1‑64 (MQ=255)                                                                                                                     ttGGTTTTCTCCCTGCGTTAACGCACAATACGCCACTTTCTTTTTCCCGGATTTTTACGTTATGTCTCACCGCGac                                                                                >  2:1052581/1‑76 (MQ=255)                                                                                                                     ttGGTTTTCTCCCTGCGTTAACGCACAATACGCCACTTTCTTTTTCCCGGATTTTTACGTTATGTCTCACCGCGac                                                                                <  1:1052581/76‑1 (MQ=255)                                                                                                                      tGGTTTTCTCCCTGCGTTAACGCACAATACGCCACTTGCTTTTTCACGGATTTTTACGTTATGTCTCACCGCGACACGCTATTTTCTGCCCCTATCGCCAGACTGGGCGACTGGACCTTTGATGAACGGGTAGCTGAAg                <  1:1443623/139‑1 (MQ=255)                                                                                                                       ggTTTTCTCCCTGCGTTAACGCACAATACGCCACTTTCTTTTTCCCGGATTTTTACGTTATGTCTCACCGCGACACGCTATTTTCTGCCCCTATCGCCAGACTGGGCGACTGGACCTTTGATGAACGGGTAGCTGAAGt               >  2:793861/1‑139 (MQ=255)                                                                                                                           ttCTCCCTGCGTTAACGCACAATACGCCACTTTCTTTTTCCCGGATTTTTACGTTATGCTTCACCGCGACACGCTATTTTCTGCCCCTATCGCCAGACTGGGCGACTGGACCTTTGATGAACGGGTAGCTGAAGTCTTc           <  2:928457/139‑1 (MQ=255)                                                                                                                             ctcCCTGCGTTAACGCACAATACGCCACTTTCTTTTTCCCGGATTTTTACGTTATGTCTCACCGCGac                                                                                <  2:1375970/68‑1 (MQ=255)                                                                                                                             ctcCCTGCGTTAACGCACAATACGCCACTTTCTTTTTCCCGGATTTTTACGTTATGTCTCACCGCGac                                                                                >  1:1375970/1‑68 (MQ=255)                                                                                                                             ctcCCTGCGTTAACGCACAATACGCCACTTTCTTTTTCCCGGATTTTTACGTTATGTCTCACCGCGACACGCTATTTTCTGCCCCTATCGCCAGACTGGGCGACTGGACCTTTGATGAACGGGTAGCTGAAGTCTTccc         >  1:1287817/1‑139 (MQ=255)                                                                                                                                   gCGTTAACGCACAATACGCCACTTTCTTTTTCCCGGATTTTTACg                                                                                                 >  2:199491/1‑45 (MQ=255)                                                                                                                                   gCGTTAACGCACAATACGCCACTTTCTTTTTCCCGGATTTTTACg                                                                                                 <  1:199491/45‑1 (MQ=255)                                                                                                                                    cGTTAACGCACAATACGCCACTTTCTTTTTCCCGGATTTTTACGTTATGTCTCACCGCGACACGCTATTTTCTGCCCCTATCGCCAGACTGGGCGACTGGACCTTTGATGAACGGGTAGCTGAAGTCTTCCCGGATATg  <  2:696947/139‑1 (MQ=255)                                                                                                                                         |                                                                                                                                      TGCATTATTACGGTTTTCATCACCGTCTGTTCCGCTGGCGACGTTCCGGCATGAGCGCCACCTGGTGTGCGCTGTTGCTGATGGTGCTGGCGAATCTTTGGTATATGCCCTGGGAGTTGGTTTTCTCCCTGCGTTAGCGCACAATACGCCACTTTCTTTTTCCCGGATTTTTACGTTATGTCTCACCGCGACACGCTATTTTCTGCCCCTATCGCCAGACTGGGCGACTGGACCTTTGATGAACGGGTAGCTGAAGTCTTCCCGGATATG  >  NC_000913/1952525‑1952794

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: ‑