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

ATTACGGTTTTCATCACCGTCTGTTCCGCTGGCGACGTTCCGGCATGAGCGCCACCTGGTGTGCGCTGTTGCTGATGGTGCTGGCGAATCTTT‑GGTATATGCCCTGGGAGTTGGTTTTCTCCCTGCGTTAGCGCACAATACGCCACTTTCTTTTTCCCGGATTTTTACGTTATGTCTCACCGCGACACGCTATTTTCTGCCCCTATCGCCAGACTGGGCGACTGGACCTTTGATGAACGGGTAGCTGAAGTCTTCCCGGATAT  >  NC_000913/1952531‑1952793                                                                                                                                    |                                                                                                                                     attaCGGTTTTCATCACCGTCTGTTCCGCTGGCGACGTTCCGGCATGAGCGCCACCTGGTGTGCGCTGTTGCTGATGGTGCTGGCGAATCTTT‑GGTATATGCCCTGGGAGTTGGTTTTCTCCCTGCGTTAACGCACAAt                                                                                                                              >  2:285042/1‑139 (MQ=255)          ttCATCACCGTCTGTTCCGCTGGCGACGTTCCGGCATGAGCGCCACCTGGTGTGCGCTGTTGCTGATGGTGCTGGCGAATCTTT‑GGTATATGCCCTGGGAGTTGGTTTTCTCCCTGCGTTAACGCACAATACGCCActt                                                                                                                     >  1:1122098/1‑139 (MQ=255)          ttCATCACCGTCTGTTCCGCTGGCGACGTTCCGGCATGAGCGCCACCTGGTGTGCGCTGTTGCTGATGGTGCTGGCGAATCTTT‑GGTATATGCCCTGGGAGTTGGTTTTCTCCCTGCGTTAACGCACAATACGCCActt                                                                                                                     >  1:1294310/1‑139 (MQ=255)                aCCGTCTGTTCCGCTGGCGACGTTCCGGCATGAGCGCCACCTGGTGTGCGCTGTTGCTGATGGTGCTGGCGAATCTTT‑GGTATATGCCCTGGGAGTTGGTTTTCTCCCTGCGTTAACGCACAATACGCCACTTTCtttt                                                                                                               >  1:1337062/1‑139 (MQ=255)                    tCTGTTCCGCTGGCGACGTTCCGGCATGAGCGCCACCTGGTGTGCGCTGTTGCTGATGGTGCTGGCGAATCTTT‑GGTATATGCCCTGGGAGTTGGTTTTCTCCCTGCGTTAACGCACAATACGCCACTTTCTTTTTccc                                                                                                           >  2:322426/1‑139 (MQ=255)                    tCTGTTCCGCTGGCGACGTTCCGGCATGAGCGCCACCTGGTGTGCGCTGTTGCTGATGGTGCTGGCGAATCTTT‑GGTATATGCCCTGGGAGTTGGTTTTCTCCCTGCGTTAACGCACAATACGCCACTTTCTTTTTccc                                                                                                           >  2:1027374/1‑139 (MQ=255)                     cTGTTCCGCTGGCGACGTTCCGGCATGAGCGCCACCTGGTGTGCGCTGTTGCTGATGGTGCTGGCGAATCTTT‑GGTATATGCCCTGGGAGTTGGTTTTCTCCCTGCGTTAACGCACAATACGCCACTTTCTTTTTCCCg                                                                                                          <  1:322426/139‑1 (MQ=255)                             cTGGCGACGTTCCGGCATGAGCGCCACCTGGTGTGCGCTGTTGCTGATGGTGCTGGCGAATCTTT‑GGTATATGCCCTGGGAGTTGGTTTTCTCCCTGCGTTAACGCACAATACGCCACTTTCTTTTTccc                                                                                                           <  1:1024426/130‑1 (MQ=255)                             cTGGCGACGTTCCGGCATGAGCGCCACCTGGTGTGCGCTGTTGCTGATGGTGCTGGCGAATCTTT‑GGTATATGCCCTGGGAGTTGGTTTTCTCCCTGCGTTAACGCACAATACGCCACTTTCTTTTTccc                                                                                                           >  2:1024426/1‑130 (MQ=255)                             cTGGCGACGTTCCGGCATGAGCGCCACCTGGTGTGCGCTGTTGCTGATGGTGCTGGCGAATCTTT‑GGTATATGCCCTGGGAGTTGGTTTTCTCCCTGCGTTAACGCACAATACGCCACTTTCTTTTTCCCGGAtttttg                                                                                                  <  2:1294310/139‑2 (MQ=255)                             cTGGCGACGTTCCGGCATGAGCGCCACCTGGTGTGCGCTGTTGCTGATGGTGCTGGCGAATCTTT‑GGTATATGCCCTGGGAGTTGGTTTTCTCCCTGCGTTAACGCACAATACGCCACTTTCTTTTTCCCGGATTTTTa                                                                                                  >  2:1738163/1‑139 (MQ=255)                                              tGAGCGCCACCTGGTGTGCGCTGTTGCTGATGGTGCTGGCGAATCTTT‑GGTATATGCCCTGGGAGTTGGTTTTCTCCCTGCGTTAACGCACAATACGCCACTTTCTTTTTCCCGGATTTTTACGTTATGTCTCACcgcg                                                                                 <  2:1337062/139‑1 (MQ=255)                                                aGCGCCACCTGGTGTGCGCTGTTGCTGATGGTGCTGGCGAATCTTT‑GGTATATGCCCTGGGAGTTGGTTTTCTCCCTGCGTTAACGCACAATACGCCACTTTCTTTTTCCCGGATTTTTACGTTATGTCTCACCGCGac                                                                               <  1:1027374/139‑1 (MQ=255)                                                 gcgcCACCTGGTGTGCGCTGTTGCTGATGGTGCTGGCGAATCTTT‑GGTATATGCCCTGGGAGTTGGTTTTCTCCCTGCGTTAACGCACAATACGCCACTTTCTTTTTCCCGGATTTTTACGTTATGTCTCACCGCGaca                                                                              >  2:313986/1‑139 (MQ=255)                                                  cgcCACCTGGTGTGCGCTGTTGCTGATGGTGCTGGCGAATCTTT‑GGTATATGCCCTGGGAGTTGGTTTTCTCCCTGCGTTAACGCACAATACGCCACTTTCTTTTTCCCGGATTTTTACGTTATGTCTCACCGCGacac                                                                             <  1:59551/139‑1 (MQ=255)                                                    ccACCTGGTGTGCGCTGTTGCTGATGGTGCTGGCGAATCTTT‑GGTATATGCCCTGGGAGTTGGTTTTCTCCCTGCGTTAACGCACAATACGCCACTTTCTTTTTCCCGGATTTTTACGTTATGTCTCACCGCGACACGc                                                                           <  1:285042/139‑1 (MQ=255)                                                            tgtgCGCTGTTGCTGATGGTGCTGGCGAATCTTT‑GGTATATGCCCTGGGAGTTGGTTTTCTCCCTGCGTTAACGCACAATACGCCACTTTCTTTTTccc                                                                                                           >  2:1939905/1‑99 (MQ=255)                                                            tgtgCGCTGTTGCTGATGGTGCTGGCGAATCTTT‑GGTATATGCCCTGGGAGTTGGTTTTCTCCCTGCGTTAACGCACAATACGCCACTTTCTTTTTccc                                                                                                           <  1:1939905/99‑1 (MQ=255)                                                                cgcTGTTGCTGATGGTGATGGCGAATCTTT‑GGTATATGCCCTGGGAGTTGGTTTTCTCCCTGCGTTAACGCACAATACGCCACTTTCTTTTTCCCGGATTTTTACGTTATGTCTCACCGCGACACTCTATTTTCTGccc                                                               >  2:1285596/1‑139 (MQ=255)                                                                 gcTGTTGCTGATGGTGCTGGCGAATCTTT‑GGTATATGCCCTGGGAGTTGGTTTTCTCCCTGCGTTAACGCACAATACGCCACTTTCTTTTTCCCGGATTTTTACGTTATGTCTCAcc                                                                                    >  2:966621/1‑117 (MQ=255)                                                                 gcTGTTGCTGATGGTGCTGGCGAATCTTT‑GGTATATGCCCTGGGAGTTGGTTTTCTCCCTGCGTTAACGCACAATACGCCACTTTCTTTTTCCCGGATTTTTACGTTATGTCTCAcc                                                                                    <  1:966621/117‑1 (MQ=255)                                                                 gcTGTTGCTGATGGTGCTGGCGAATCTTT‑GGTATATGCCCTGGGAGTTGGTTTTCTCCCTGCGTTAACGCACAATACGCCACTTTCTTTTTCCCGGATTTTTACGTTATGTCTCACCGCGACACGCTATTTTCTGcccc                                                              <  2:1677492/139‑1 (MQ=255)                                                                   tGTTGCTGATGGTGCTGGCGAATCTTT‑GGTATATGCCCTGGGAGTTGGTTTTCTCCCTGCGTTAacg                                                                                                                                    >  1:244356/1‑67 (MQ=255)                                                                     ttGCTGATGGTGCTGGCGAATCTTTGGGTATATGCCCTGGGAGTTGGTTTTCTCCCTGCGTTAACGCACAATACGCCACTTTCTTTTTCCCGGATTTTTACGTTATGTCTCACCGCGACACGCTATTTTCTGCCCCTAt                                                           >  2:1434297/1‑139 (MQ=255)                                                                     ttGCTGATGGTGCTGGCGAATCTTT‑GGTATATGCCCTGGGAGTTGGTTTTCTCCCTGCGTTAACGCACAATACGCCACTTTCTTTTTCCCGGATTTTTACGTTATGTCTCACCGCGACACGCTATTTTCTGCCCCTATc                                                          >  2:96348/1‑139 (MQ=255)                                                                     ttGCTGATGGTGCTGGCGAATCTTT‑GGTATATGCCCTGGGAGTTGGTTTTCTCCCTGCGTTAACGCACAATACGCCACTTTCTTTTTCCCGGATTTTTACGTTATGTCTCACCGCGACACGCTATTTTCTGCCCCTATc                                                          >  1:1148842/1‑139 (MQ=255)                                                                      tGCTGATGGTGCTGGCGAATCTTT‑GGTATATGCCCTGGGAGTTGGTTTTCTCCCTGCGTTAACGCACAATACGCCACTTTCTTTTTCCCGGATTTTTACGTTATg                                                                                           <  2:1676621/105‑1 (MQ=255)                                                                      tGCTGATGGTGCTGGCGAATCTTT‑GGTATATGCCCTGGGAGTTGGTTTTCTCCCTGCGTTAACGCACAATACGCCACTTTCTTTTTCCCGGATTTTTACGTTATg                                                                                           >  1:1676621/1‑105 (MQ=255)                                                                       gCTGATGGTGCTGGCGAATCTTT‑GGTATATGCCCTGGGAGTTGGTTTTCTCCCTGCGTTAACGCACAATACGCCACTTTCTTTTTCCCGGATTTTTACGTTATGTCTCACCGCGACACGCTATTTTCTGCCCCTATCGc                                                        <  2:1122098/139‑1 (MQ=255)                                                                             ggTGCTGGCGAATCTTT‑GGTATATGCCCTGGGAGTTGGTTTTCTCCCCGCGTTAACGCACAATACGCCACTTTCTTTTTCCCGGATTTTTACGTTATGTCTCACCGCGACACGCTATTTTCTGCCCCTATCGCCAGACt                                                  >  2:336891/1‑139 (MQ=255)                                                                                gCTGGCGAATCTTT‑GGTATATGCCCTGGGAGTTGGTTTTCTCCCTGCGTTAACGCACAATACGCCACTTTCTTTTTCCCGGATTTTTACGTTATGTCTCACcgcg                                                                                 >  2:92736/1‑105 (MQ=255)                                                                                gCTGGCGAATCTTT‑GGTATATGCCCTGGGAGTTGGTTTTCTCCCTGCGTTAACGCACAATACGCCACTTTCTTTTTCCCGGATTTTTACGTTATGTCTCACcgcg                                                                                 <  1:92736/105‑1 (MQ=255)                                                                                 cTGGCGAATCTTT‑GGTATATGCCCTGGGAGTTGGTTTTCTCCCTGCGTTAACGCACAATACGCCACTTTCTTTTTCCCGGATTTTTACGTTATGTCTCACCGCGACACGCTATTTTCTGCCCCTATCGCCAGACTGGGc                                              >  2:1204324/1‑139 (MQ=255)                                                                                  tGGCGAATCTTT‑GGTATATGCCCTGGGAGTTGGTTTTCTCCCTGCGTTAACGca                                                                                                                                  >  1:49997/1‑54 (MQ=255)                                                                                  tGGCGAATCTTT‑GGTATATGCCCTGGGAGTTGGTTTTCTCCCTGCGTTAACGca                                                                                                                                  <  2:49997/54‑1 (MQ=255)                                                                                    gCGAATCTTT‑GGTATATGCCCTGGGAGTTGGTTTTCTCCCTGCGTTAACGCACAATACGCCACTTTCTTTTTCCCGGATTTTTACGTTATGTCTCACCGCGACACGCTATTTTCTGCCCCTATCGCCAGACTGGGCGAc                                           <  1:1285596/139‑1 (MQ=255)                                                                                                    atGCCCTGGGAGTTGGTTTTCTCCCTGCGTTAACGCACAATACGCCACTTTCTTTTTCCCGGATTTTTACGTTATGTCTCACCGCGACACGCTATTTTCTGCCCCTATCGcc                                                       <  2:306667/112‑1 (MQ=255)                                                                                                    atGCCCTGGGAGTTGGTTTTCTCCCTGCGTTAACGCACAATACGCCACTTTCTTTTTCCCGGATTTTTACGTTATGTCTCACCGCGACACGCTATTTTCTGCCCCTATCGcc                                                       >  1:306667/1‑112 (MQ=255)                                                                                                     tGCCCTGGGAGTTGGTTTTCTCCCTGCGTTAACGCACAATACGCCACTTTCTTTTTCCCGGATTTTTACGTTATGTCTCACCGCGaca                                                                              >  2:96570/1‑88 (MQ=255)                                                                                                     tGCCCTGGGAGTTGGTTTTCTCCCTGCGTTAACGCACAATACGCCACTTTCTTTTTCCCGGATTTTTACGTTATGTCTCACCGCGaca                                                                              <  1:96570/88‑1 (MQ=255)                                                                                                     tGCCCTGGGAGTTGGTTTTCTCCCTGCGTTAACGCACAATACGCCACTTTCTTTTTCCCGGATTTTTACGTTATGTCTCACCGCGACACGCt                                                                          >  1:1635988/1‑92 (MQ=255)                                                                                                     tGCCCTGGGAGTTGGTTTTCTCCCTGCGTTAACGCACAATACGCCACTTTCTTTTTCCCGGATTTTTACGTTATGTCTCACCGCGACACGCt                                                                          <  2:1635988/92‑1 (MQ=255)                                                                                                        ccTGGGAGTTGGTTTTCTCCCTGCGTTAACGCACAATACGCCACTTTCTTTTTCCCGGATTTTTACGTTATGTCTCACCGCGACACGCTATTTTCTGCCCCTATCGcc                                                       >  1:1330565/1‑108 (MQ=255)                                                                                                        ccTGGGAGTTGGTTTTCTCCCTGCGTTAACGCACAATACGCCACTTTCTTTTTCCCGGATTTTTACGTTATGTCTCACCGCGACACGCTATTTTCTGCCCCTATCGcc                                                       <  2:1330565/108‑1 (MQ=255)                                                                                                          tGGGAGTTGGTTTTCTCCCTGCGTTAACGCACAATACGCCACTTTCTTTTTCCCGGATTTTTACGTTATGTCTCACCGCGACACGCTATTTTCTGCCCCTATCGCCAGACTGGGCGACTGGACCTTTGATGAACGGGTa                      >  2:1059132/1‑139 (MQ=255)                                                                                                           gggAGTTGGTTTTCTCCCTGCGTTAACGCACAATACGCCACTTTCTTTTTCCCGGAtttt                                                                                                    <  1:1665556/60‑1 (MQ=255)                                                                                                           gggAGTTGGTTTTCTCCCTGCGTTAACGCACAATACGCCACTTTCTTTTTCCCGGAtttt                                                                                                    >  2:1665556/1‑60 (MQ=255)                                                                                                           gggAGTTGGTTTTCTCCCTGCGTTAACGCACAATACGCCACTTTCTTTTTCCCGGATTTTTACGTTATGtctc                                                                                       >  2:659401/1‑73 (MQ=255)                                                                                                           gggAGTTGGTTTTCTCCCTGCGTTAACGCACAATACGCCACTTTCTTTTTCCCGGATTTTTACGTTATGtctc                                                                                       <  1:659401/73‑1 (MQ=255)                                                                                                           gggAGTTGGTTTTCTCCCTGCGTTAACGCACAATACGCCACTTTCTTTTTCCCGGATTTTTACGTTATGTCTCACCGCGac                                                                               <  1:724328/81‑1 (MQ=255)                                                                                                           gggAGTTGGTTTTCTCCCTGCGTTAACGCACAATACGCCACTTTCTTTTTCCCGGATTTTTACGTTATGTCTCACCGCGac                                                                               >  2:724328/1‑81 (MQ=255)                                                                                                           gggAGTTGGTTTTCTCCCTGCGTTAACGCACAATACGCCACTTTCTTTTTCCCGGATTTTTACGTTATGTCTCACCGCGACACGCTATTTTCTGCCCCTATCGc                                                        >  2:1344104/1‑104 (MQ=255)                                                                                                           gggAGTTGGTTTTCTCCCTGCGTTAACGCACAATACGCCACTTTCTTTTTCCCGGATTTTTACGTTATGTCTCACCGCGACACGCTATTTTCTGCCCCTATCGc                                                        <  1:1344104/104‑1 (MQ=255)                                                                                                               gTTGGTTTTCTCCCTGCGTTAACGCACAATACGCCACTTTCTTTTTCCCGGAtttt                                                                                                    <  2:602288/56‑1 (MQ=255)                                                                                                               gTTGGTTTTCTCCCTGCGTTAACGCACAATACGCCACTTTCTTTTTCCCGGAtttt                                                                                                    >  1:602288/1‑56 (MQ=255)                                                                                                                ttGGTTTTCTCCCTGCGTTAACGCACAATACGCCACTTTCTTTTTCCCGGATTTTTACGTTATGTCTCACCGCGACACGCTATTTTCTGTCCCTATCGCCAGACTGGGCGACTGGACCTTTGATGAACGGGTAGCTGaa                >  1:349448/1‑139 (MQ=255)                                                                                                                  ggTTTTCTCCCTGCGTTAACGCACAATACGCCACTTTCTTTTTCCCGGAtttt                                                                                                    <  2:618468/53‑1 (MQ=255)                                                                                                                  ggTTTTCTCCCTGCGTTAACGCACAATACGCCACTTTCTTTTTCCCGGAtttt                                                                                                    >  1:618468/1‑53 (MQ=255)                                                                                                                  ggTTTTCTCCCTGCGTTAACGCACAATACGCCACTTTCTTTTTCCCGGATTTTTa                                                                                                  <  2:1248222/55‑1 (MQ=255)                                                                                                                  ggTTTTCTCCCTGCGTTAACGCACAATACGCCACTTTCTTTTTCCCGGATTTTTa                                                                                                  >  1:1248222/1‑55 (MQ=255)                                                                                                                  ggTTTTCTCCCTGCGTTAACGCACAATACGCCACTTTCTTTTTCCCGGATTTTTACGTTATGTCTCACCGCGACACGCTATTTTCTGCCCCTATCGCCAGACTGGGCGACTGGACCTTTGATGAACGGGTAGCTGAAGt              >  1:1111642/1‑139 (MQ=255)                                                                                                                  ggTTTTCTCCCTGCGTTAACGCACAATACACCACTTTCTTTTTCCCg                                                                                                          >  1:1968287/1‑47 (MQ=255)                                                                                                                  ggTTTTCTCCCTGCGTTAACGCACAATACACCACTTTCTTTTTCCCg                                                                                                          <  2:1968287/47‑1 (MQ=255)                                                                                                                     tttCTCCCTGCGTTAACGCACAATACGCCACTTTCTTTTTCCCGGATTTTTACGTTATGTCTCACCGCGACACGCTATTTTCTGCCCCTATCGCCAGACTGGGCGACTGGACCTTTGATGAACGGGTAGCTGAAGTCtt           <  2:1111642/139‑1 (MQ=255)                                                                                                                      ttCTCCCTGCGTTAACGCAGAATACGCCACTTTTTTTTTCCCGGATTTTTACGTTATGTCTCACCGCGACACGCTATTTTCTGCCCCTATCGCCAGACTGGGCGACTGGACCTTTGATGAACGGGTAGCTGAAGTCTTc          <  2:311723/139‑1 (MQ=255)                                                                                                                      ttCTCCCTGCGTTAACGCACAATACGCCATTTTCTTTTTCCCGGATTTTTACGTTATGTCTCACCGCGACACGCTATATTCTGCCCCTACCGCCAGACTGGGCGACTGGATCTTTGATGAACGGGTAGCTGAAGTCTTc          <  1:1329291/139‑1 (MQ=255)                                                                                                                       tctcCCTGCGTTAACGCACAATACGCCACTTTCTTTTTCCCGGATTTTTACGTTATGTCTCACCGCGACACGCTATTTTCTGCCCCTATCGCCAGACTGGGCGACTGGACCTTTGATGAACGGGTAGCTGAAGTCTTcc         <  1:1434297/139‑1 (MQ=255)                                                                                                                       tctcCCTGCGTTAACGCACAATACGCCACTTTCTTTTTCCCGGATTTTTACGTTATGTCTCACCGCGACACGCTATTTTCTGCCCCTATCGCCAGACTGGGCGACTGGACCTTTGATGAACGGGTAGCTGAAGTCTTcc         <  1:1204324/139‑1 (MQ=255)                                                                                                                              gCGTTAACGCACAATACGCCACTTTTTTTTTCCCGGATTTTTACGTTATGTCTCACCGCGACACGCTATTTTCTGCCCCTATCGCCAGACTGGGCGACTGGACTTTTGATGACTGGGTATCCGAAGTATTCCCGGatat  >  1:311723/1‑139 (MQ=255)                                                                                                                                    |                                                                                                                                     ATTACGGTTTTCATCACCGTCTGTTCCGCTGGCGACGTTCCGGCATGAGCGCCACCTGGTGTGCGCTGTTGCTGATGGTGCTGGCGAATCTTT‑GGTATATGCCCTGGGAGTTGGTTTTCTCCCTGCGTTAGCGCACAATACGCCACTTTCTTTTTCCCGGATTTTTACGTTATGTCTCACCGCGACACGCTATTTTCTGCCCCTATCGCCAGACTGGGCGACTGGACCTTTGATGAACGGGTAGCTGAAGTCTTCCCGGATAT  >  NC_000913/1952531‑1952793

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