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%	109.1 / NA	35	*132* (TAG→TAA)	yecN	MAPEG family inner membrane protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (18/17);  total (18/17)

GTTTTCATCACCGTCTGTTCCGCTGGCGACGTTCCGGCATGAGCGCCACCTGGTGTGCGCTGTTGCTGATGGTGCTGGCGAATCTTTGGTATATGCCCTGGGAGTTGGTTTTCTCCCTGCGTTAGCGCACAATACGCCACTTTCTTTTTCCCGGATTTTTACGTTATGTCTCACCGCGACACGCTATTTTCTGCCCCTATCGCCAGACTGGGCGACTGGACCTTTGATGAACGGGTAGCTGAAGTCTTCCCGG  >  NC_000913/1952537‑1952789                                                                                                                             |                                                                                                                                 gTTTTCATCACCGTCTGTTCCGCTGGCGACGTTCCGGCATGAGCGCCACCTGGTGTGCGCTGTTGCTGATGGTGCTGGCGAATCTTTGGTATATGCCCTGGGAGTTGGTTTTCTCCCTGCGTTAACGCACAATACGCCa                                                                                                                    >  2:915721/1‑139 (MQ=255)  ttttCATCACCGTCTGTTCCGCTGGCGACGTTCCGGCATGAGCGCCACCTGGTGTGCGCTGTTGCTGATGGTGCTGGCGAATCTTTGGTATATGCCCTGGGAGTTGGTTTTCTCCCTGCGTTAACGCACAATACGCCAc                                                                                                                   >  1:113984/1‑139 (MQ=255)      catcaCCGTCTGTTCCGCTGGCGACGTTCCGGCATGAGCGCCACCTGGTGTGCGCTGTTGCTGATGGTGCTGGCGAATCTTTGGTATATGCCCTGGGAGTTGGTTTTCTCCCTGCGTTAACGCACAATACGCCActttc                                                                                                               >  2:823106/1‑139 (MQ=255)              tCTGTTCCGCTGGCGACGTTCCGGCATGAGCGCCACCTGGTGTGCGCTGTTGCTGATGGTGCTGGCGAATCTTTGGTATATGCCCTGGGAGTTGGTTTTCTCCCTGCGTTAacg                                                                                                                                >  1:352588/1‑114 (MQ=255)              tCTGTTCCGCTGGCGACGTTCCGGCATGAGCGCCACCTGGTGTGCGCTGTTGCTGATGGTGCTGGCGAATCTTTGGTATATGCCCTGGGAGTTGGTTTTCTCCCTGCGTTAacg                                                                                                                                <  2:398244/114‑1 (MQ=255)              tCTGTTCCGCTGGCGACGTTCCGGCATGAGCGCCACCCGGTGTGCGCTGTTGCTGATGGTGCTGGCGAATCTTTGGTATATGCCCTGGGAGTTGGTTTTCTCCCTGCGTTAACGCACAATACGCCACTTTCTTTTTccc                                                                                                       >  2:492126/1‑139 (MQ=255)                   tCCGCTGGCGACGTTCCGGCATGAGCGCCACCCGGTGTGCGCTGTTGCTGATGGTGCTGGCGAATCTTTGGTATATGCCCTGGGAGTTGGTTTTCTCCCTGCGTTAACGCACAATACGCCACTTTCTTTTTCCCGGAtt                                                                                                  <  1:492126/139‑1 (MQ=255)                                            cgcCACCTGGTGTGCGCTGTTGCTGATGGTGCTGGCGAATCTTTGGTATATGCCCTGGGAGTTGGTTTTCTCCCTGCGTTAACGCACAATACGCCACTTTCTTTTTCCCGGATTTTTACg                                                                                            <  2:568049/120‑1 (MQ=255)                                            cgcCACCTGGTGTGCGCTGTTGCTGATGGTGCTGGCGAATCTTTGGTATATGCCCTGGGAGTTGGTTTTCTCCCTGCGTTAACGCACAATACGCCACTTTCTTTTTCCCGGATTTTTACg                                                                                            >  1:568049/1‑120 (MQ=255)                                            cgcCACCTGGTGTGCGCTGTTGCTGATGGTGCTGGCGAATCTTTGGTATATGCCCTGGGAGTTGGTTTTCTCCCTGCGTTAACGCACAATACGCCACTTTCTTTTTCCCGGATTTTTACGTTATGTCTCACCGCGacac                                                                         >  1:992282/1‑139 (MQ=255)                                                 ccTGGTGTGCGCTGTTGCTGATGGTGCTGGCGAATCTTTGGTATATGCCCTGGGAGTTGGTTTTCTCCCTGCGTTAACGCACAATACGCCACTTTCTTTTTCCCGGATTTTTACGTTATGTCTCACCGCGACACGCTAt                                                                    >  1:5119/1‑139 (MQ=255)                                                           gcTGTTGCTGATGGTGCTGGCGAATCTTTGGTATATGCCCTGGGAGTTGGTTTTCTCCCTGCGTTAACGCACAATACGCCACTTTCTTTTTCCCGGATTTTTACGTTATGTCTCACCGCGACACGCTATTTTCTGcccc                                                          <  1:823106/139‑1 (MQ=255)                                                            cTGTTGCTGATGGTGCTGGCGAATCTTTGGTATATGCCCTGGGAGTTGGTTTTCTCCCTGCGTTAACGCACAATACGCCACTTTGTTTTTCCCGGATTTTTACGTTATGTCTCACCGCGACACGCTATTTTCTGCCCCt                                                         >  1:664360/1‑139 (MQ=255)                                                            cTGTTGCTGATGGTGCTGGCGAATCTTTGGTATATGCCCTGGGAGTTGGTTTTCTCCCTGCGTTAACGCACAATACGCCACTTTCTTTTTCCCGGATTTTTACGTTATGTCTCACCGCGACACGCTATTTTCTGCCCCt                                                         >  2:168253/1‑139 (MQ=255)                                                               ttGCTGATGGTGCTGGCGAATCTTTGGTATATGCCCTGGGAGTTGGTTTTCTCCCTGCGTTAACGCACAATACGCCACTTTCTTTTTCCCGGATTTTTACGTTATGTCTCACCGCGACACGCTATTTTCTGCCCCTATc                                                      >  2:246559/1‑139 (MQ=255)                                                                       ggTGCTGGCGAATCTTTGGTATATGCCCTGGGAGTTGGTTTTCTCCCTGCGTTAACGCACAATACGCCACTTTCTTTTTCCCGGATTTTTACGTTATGTCTCACCGCGACACGCTATTTTCTg                                                              >  1:914796/1‑123 (MQ=255)                                                                       ggTGCTGGCGAATCTTTGGTATATGCCCTGGGAGTTGGTTTTCTCCCTGCGTTAACGCACAATACGCCACTTTCTTTTTCCCGGATTTTTACGTTATGTCTCACCGCGACACGCTATTTTCTg                                                              <  2:914796/123‑1 (MQ=255)                                                                        gTGCTGGCGAATCTTTGGTATATGCCCTGGGAGTTGGTTTTCTCCCTGCGTTAACGCACAATACGCCACTTTCTTTTTCCCGGATTTTTACGTTATGTCTCACCGCGACACGCTATTTTCTGCCCCTATCGCCAGACTg                                             <  2:5119/139‑1 (MQ=255)                                                                         tGCTGGCGAATCTTTGGTATATGCCCTGGGAGTTGGTTTTCTCCCTGCGTTAACGCACAATACGCCACTTTCTTTTTCCCGGATTTTTACGTTATGTCTCACCGCGACACGCTATTTTCTGCCCCTATCGCCAGACTgg                                            >  1:207795/1‑139 (MQ=255)                                                                              gCGAATCTTTGGTATATGCCCTGGGAGTTGGTTTTCTCCCTGCGTTAACGCACAATACGCCACTTTCTTTTTCCCGGATTTTTACGTTATGTCTCACCGCGACACGCTATTTTCTGCCCCTATCGCCAGACTGGGCGAc                                       <  2:113984/139‑1 (MQ=255)                                                                                gAATCTTTGGTATATGCCCTGGGAGTTGGTTTTCTCCCTGCGTTAACGCACAATACGCCACTTTCTTTTTCCCGGATTTTTACGTTATGTCTCACCGCGACACGCTATTTTCTGCCCCTAACGCCAGACTGGGCGACTg                                     <  2:664360/139‑1 (MQ=255)                                                                                        ggTATATGCCCTGGGAGTTGGTTTTCTCCCTGCGTTAATGCACAATACGCCACTTTCTTTTTCCCGGATTTTTACGTTa                                                                                         <  2:614963/79‑1 (MQ=255)                                                                                        ggTATATGCCCTGGGAGTTGGTTTTCTCCCTGCGTTAACGCACAATACGCCACTTTCTTTTTCCCGGATTTTTACGTTa                                                                                         >  1:614963/1‑79 (MQ=255)                                                                                            tatGCCCTGGGAGTTGGTTTTCTCCCTGCGTTAACGCACAATACGCCACTTTCTTTTTCCCGGATTTTTACGTTATGTCTCACCGCGACACGCTATTTTCTGCCCCTATCGCCAGACTGGGCGACTGGACCTTtgatga                         <  1:168253/139‑1 (MQ=255)                                                                                             atGCCCTGGGAGTTGGTTTTCTCCCTGCGTTAACGCACAATACGCCACTTTCTTTTTCCCGGATTTTTACGTTATGTCTCACCGCGACACGCTATTTTCTGCCCCTATCGCCAGACTGGGCGACTGGACCTTTGATGaa                        >  1:161067/1‑139 (MQ=255)                                                                                              tGCCCTGGGAGTTGGTTTTCTCCCTGCGTTAACGCACAATACGCCACTTTCTTTTTCCCGGATTTTTACGTTATGTCTCACCGCGACACGCTATTTTCTGCCCCTATCGCCAGACTGGGCGACTGGACCTTTGATGaa                        <  1:528186/138‑1 (MQ=255)                                                                                              tGCCCTGGGAGTTGGTTTTCTCCCTGCGTTAACGCACAATACGCCACTTTCTTTTTCCCGGATTTTTACGTTATGTCTCACCGCGACACGCTATTTTCTGCCCCTATCGCCAGACTGGGCGACTGGACCTTTGATGaa                        >  2:528186/1‑138 (MQ=255)                                                                                                    gggAGTTGGTTTTCTCCCTGCGTTAACGCACAATACGCCACTTTCTTTTTCCCGGATTTTTACGTTATGTCTCACCGCGACACGc                                                                       <  2:43396/85‑1 (MQ=255)                                                                                                    gggAGTTGGTTTTCTCCCTGCGTTAACGCACAATACGCCACTTTCTTTTTCCCGGATTTTTACGTTATGTCTCACCGCGACACGc                                                                       >  1:43396/1‑85 (MQ=255)                                                                                                         ttGGTTTTCTCCCTGCGTTAACGCACAATACGCCACTTTCTTTTTCCCGGATTTTTACGTTATGTCTCACCGCGACACGCTATTTTCTGCCCCTATCGCCAGACTGGGCGGCTGGACCTTTGATGAACGGGTAGCTGaa            <  2:146249/139‑1 (MQ=255)                                                                                                         ttGGTTTTCTCCCTGCGTTAACGCACAATACGCCACTTTCTTTTTCCCGGATTTTTACGTTATGTCTCACCGCGACACGCTATTTTCTGCCCCTATCGCCAGACTGGGCGACTGGACCTTTGATGAACGGGTAGCTGaa            <  2:358393/139‑1 (MQ=255)                                                                                                         ttGGTTTTCTCCCTGCGTTAACGCACAATACGCCACTTTCTTTTTCCCGGATTTTTACGTTATGTCTCACCGCGACACGCTATTTTCTGCCCCTATCGCCAGACTGGGCGACTGGACCTTTGATGAACGGGTAGCTGaa            <  2:207795/139‑1 (MQ=255)                                                                                                           ggTTTTCTCCCTGCGTTAACGCACAATACGCCACTTTCTTTTTCCCGGATTTTTACGTTATGTCTCACCGCGACACGCTATTTTCTGCCCCTATCGCCAGACTGGGCGACTGGACCTTTGATGAACGGGTAGCTGAAGt          >  1:996306/1‑139 (MQ=255)                                                                                                            gTTTTCTCCCTGCGTTAACGCACAATACGCCACTTTCTTTTTCCCGGATTTTTACGTTATGTCTCACCGCGACACGCTATTTTCTGCCCCTATCGCCAGACTGGGCGACTGGACCTTTGATGAACGGGTAGCTGAAGTc         <  1:915721/139‑1 (MQ=255)                                                                                                               ttCTCCCTGCGTTAACGCACAATACGCCACTTTCTTTTTCCCGGATTTTTACGTTATGTCTCACCGCGACACGCTATTTTCTGCCCCTATCGCCAGACTGGGCGACTGGACCTTTGATGAACGGGTAGCTGAAGTCTTc      <  2:992282/139‑1 (MQ=255)                                                                                                                   cccTGCGTTAACGCACAATACGCCACTTTCTTTTTCCCGGATTTTTACGTTATGTCTCACCGCGACACGCTATTTTCTGCCCCTATCGCCAGACTGGGCGACTGGACCTTTGATGAACGGGTAGCTGAAGTCTTCCCgg  <  2:701800/139‑1 (MQ=255)                                                                                                                   cccTGCGTTAACGCACAATACGCCACTTTCTTTTTCCCGGATTTTTACGTTATGTCTCACCGCGACACGCTATTTTCTGCCCCTATCGCCAGACTGGGCGACTGGACCTTTGATGAACGGGTAGCTGAAGTCTTCCCgg  >  2:912311/1‑139 (MQ=255)                                                                                                                             |                                                                                                                                 GTTTTCATCACCGTCTGTTCCGCTGGCGACGTTCCGGCATGAGCGCCACCTGGTGTGCGCTGTTGCTGATGGTGCTGGCGAATCTTTGGTATATGCCCTGGGAGTTGGTTTTCTCCCTGCGTTAGCGCACAATACGCCACTTTCTTTTTCCCGGATTTTTACGTTATGTCTCACCGCGACACGCTATTTTCTGCCCCTATCGCCAGACTGGGCGACTGGACCTTTGATGAACGGGTAGCTGAAGTCTTCCCGG  >  NC_000913/1952537‑1952789

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