Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,495,050	C→T	*81* (TAG→TAA)	yfdY ←	DUF2545 family putative inner membrane protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,495,050	0	C	T	100.0%	123.3 / NA	39	*81* (TAG→TAA)	yfdY	DUF2545 family putative inner membrane protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (20/19);  total (20/19)

CCCCCTTTTCCGACAAAGGCGAATTATCTTATTCATCCTGTGCGATAGCGTCTTTCACATGGCAGGTGAACAATCATGTCTGATGCTGATGTGAACACCTGCCCAAATGAGCCATCCTGACCCACTACCGCCTGCCATACCCCACCATAAAGCGAATCATCAAGCCACAAGCCGGGGCAAGAAACATAAATATCCCCGTCACTAGTTCATTACCGCCAGTTAACCACATGGTGAAGTAATAAATGAGTGCA  >  NC_000913/2494926‑2495176                                                                                                                             |                                                                                                                               tccccGTTTCCGACGAAGGCGATTCTTCTTATTCATCCTGTGCGATAACGTCTTTCACATGGCAGGTGAACAATCATGTCTGATGCTGATGTGAACACCTGCCCAAAGGAGCCATCCTGACCCATTACCGCCTGCCATa                                                                                                                  <  2:306463/138‑1 (MQ=255)              cAAAGGCGAATCATCTTATTCATCCTGTGCGATAACGTCTTTCACATGGCAGGTGAACAATCATGTCTGATGCTGATGTGAACACCTGCCCAAATGAGCCATCCTGACCCATTACCGCCTGCCATACCCCACCATAAAg                                                                                                     >  1:123516/1‑139 (MQ=255)                        tcaTCTTATTCATCCTGTGCGATAACGTCTTTCACATGGCAGGTGAACAATCATGTCTGATGCTGATGTGAACACCTGCCCAAATGAGCCATCCTGACCCATTACCGCCTGCCATACCCCACCATAAAGCGAAtcatca                                                                                           <  1:228056/137‑1 (MQ=255)                        tcaTCTTATTCATCCTGTGCGATAACGTCTTTCACATGGCAGGTGAACAATCATGTCTGATGCTGATGTGAACACCTGCCCAAATGAGCCATCCTGACCCATTACCGCCTGCCATACCCCACCATAAAGCGAAtcatca                                                                                           <  1:51896/137‑1 (MQ=255)                                     ccTGTGCGATAACGTCTTTCACATGGCAGGTGAACAATCATGTCTGATGCTGATGTGAACACCTGCCCAAATGAGCCATCCTGACCCATTACCGCCTGCCATACCCCACCATAAAGCGAATCATCAAGCCACAAGCCgg                                                                              >  1:265286/1‑139 (MQ=255)                                     ccTGTGCGATAACGTCTTTCACATGGCAGGTGAACAATCATGTCTGATGCTGATGTGAACACCTGCCCAAATGAGCCATCCTGACCCATTACCGCCTGCCATACCCCACCATAAAGCGAATCATCAAGCCACAAGCCgg                                                                              >  1:646231/1‑139 (MQ=255)                                       tgtgCGATAACGTCTTTCACATGGCAGGTGAACAATCATGTCTGATGCTGATGTGAACACCTGCCCAAATGAGCCATCCTGACCCATTACCGCCTGCCATACCCCACCATAAAGCGAATCATCAAGCCACAAGCCgggg                                                                            >  2:67177/1‑139 (MQ=255)                                                 cGTCTTTCTCATGGCAGGTGAACAATCATGTCTGATGCTGATGTGAACACCTGCCCAAATGAGCCATCCTGACCCATTACCGCCTGCCATACCCCACCATAAAGCGAATCATCAAGCCACAAGCCGGGGCAAGAAACAt                                                                  <  2:433057/139‑1 (MQ=255)                                                  gTCTTTCACATGGCAGGTGAACAATCATGTCTGATGCTGATGTGAACACCTGCCCAAATGAGCCATCCTGACCCATTACCGCCTGCCATAcc                                                                                                                <  2:253646/92‑1 (MQ=255)                                                  gTCTTTCACATGGCAGGTGAACAATCATGTCTGATGCTGATGTGAACACCTGCCCAAATGAGCCATCCTGACCCATTACCGCCTGCCATAcc                                                                                                                >  1:253646/1‑92 (MQ=255)                                                          caTGGCAGGTGAACAATCATGTCTGATGCTGATGTGAACACCTGCCCAAATGAGCCATCCTGACCCATTACCGCCTGCCATACCCCACCATAAAGCGAATCATCAAGCCACAAGCCGGGGCAAGAAACATAAATATccc                                                         <  2:123516/139‑1 (MQ=255)                                                           aTGGCAGGTGAACAATCATGTCTGATGCTGATGTGAACACCTGCCCAAATGAGCCATCCTGACCCATTACCGCCTGCCATACCCCACCATAAAGCGAATCATCAAGCCACAAGCCGGGGCAAGAAACATAAATATcccc                                                        <  2:523216/139‑1 (MQ=255)                                                                aGGTGAACAATCATGTCTGATGCTGATGTGAACACCTGCCCAAATGAGCCATCCTGACCCATTACCGCCTGCCATACCCCACCATAAAGCGAATCATCAAGCCACAAGCCGGGGCAAGAAACATAAATATCCCCGTCAc                                                   <  2:357769/139‑1 (MQ=255)                                                                             tGTCTGATGCTGATGTGAACACCTGCCCAAATGAGCCATCCTGACCCATTACCGCCTGCCATACCCCACCATAAAGCGAATCATCAAGCcac                                                                                     >  1:396469/1‑92 (MQ=255)                                                                             tGTCTGATGCTGATGTGAACACCTGCCCAAATGAGCCATCCTGACCCATTACCGCCTGCCATACCCCACCATAAAGCGAATCATCAAGCcac                                                                                     <  2:396469/92‑1 (MQ=255)                                                                               tCTGATGCTGATGTGAACACCTGCCCAAATGAGCCATCCTGACCCATTACCGCCTGCCATACCCCACCATAAAGCGAATCATCAAGCCACAAGCCGGGGCAAGAAACATAAATATCCCCGTCACTAGTTCATTACCGcc                                    <  1:67177/139‑1 (MQ=255)                                                                               tCTGATGCTGATGTGAACACCTGCCCAAATGAGCCATCCTGACCCATTACCGCCTGCCATACCCCACCATAAAGCGAATCATCAAGCCACAAGCCGGGGCAAGAAACATAAATATCCCCGTCACTAGTTCATTACCGcc                                    <  1:710853/139‑1 (MQ=255)                                                                               tCTGATGCTGATGTGAACACCTGCCCAAATGAGCCATCCTGACCCATTACCGCCTGCCATACCCCACCATAAAGCGAATCATCAAGCCACAAGCCGGGGCAAGAAACATAAATATCCCCGTCACTAGTTCATTACCGcc                                    >  2:710853/1‑139 (MQ=255)                                                                                    tgctgatgTGAACACCTGCCCAAATGAGCCATCCTGACCCATTACCGCCTGCCATACCCCACCATAAAGCGAATCATCAAGCCACAAGCCGGGGCAAGAAACATAAATATCCCCGTCACTAGTTCATTAcc                                       >  1:330181/1‑131 (MQ=255)                                                                                    tgctgatgTGAACACCTGCCCAAATGAGCCATCCTGACCCATTACCGCCTGCCATACCCCACCATAAAGCGAATCATCAAGCCACAAGCCGGGGCAAGAAACATAAATATCCCCGTCACTAGTTCATTAcc                                       <  2:330181/131‑1 (MQ=255)                                                                                    tgctgatgTGAACACCTGCCCAAATGAGCCATCCTGACCCATTACCGCCTGCCATACCCCACCATAAAGCGAATCATCAAGCCACAAGCCGGGGCAAGAAACATAAATATCCCCGTCACTAGTTCATTACCGCCAGTTa                               >  1:508870/1‑139 (MQ=255)                                                                                        gatgTGAACACCTGCCCAAATGAGCCATCCTGACCCATTACCGCCTGCCATACCCCACCATAAAGCGAATCATCAAGCCACAAGCCGGGGCAAGAAACATAAATATCCCCGTCACTAGTTCATTACCGCCAGTTaa                              >  2:170893/1‑136 (MQ=255)                                                                                        gatgTGAACACCTGCCCAAATGAGCCATCCTGACCCATTACCGCCTGCCATACCCCACCATAAAGCGAATCATCAAGCCACAAGCCGGGGCAAGAAACATAAATATCCCCGTCACTAGTTCATTACCGCCAGTTaa                              <  1:170893/136‑1 (MQ=255)                                                                                        gatgTGAACACCTGCCCAAATGAGCCATCCTGACCCATTACCGCCTGCCATACCCCACCATAAAGCGAATCATCAAGCCACAAGCCGGGGCAAGAAACATAAATACCCCCGTCACTAGTTCATTACCGCCAGTTAACca                           >  2:616037/1‑139 (MQ=255)                                                                                              aaCACCTGCCCAAATGAGCCATCCTGACCCATTACCGCCTGCCATACCCCACCATAAAGCGAATCATCAAGCCACAAGCCGGGGCAAGAAACATAAATATCCCCGTCACTAGTTCATTACCGCCAGTTAACCACATGGt                     >  2:390432/1‑139 (MQ=255)                                                                                                   cTGCCCAAATGAGCCATCCTGACCCATTACCGCCTGCCATACCCCACCATAAAGCGAATCATCAAGCcac                                                                                     >  1:317802/1‑70 (MQ=255)                                                                                                   cTGCCCAAATGAGCCATCCTGACCCATTACCGCCTGCCATACCCCACCATAAAGCGAATCATCAAGCcac                                                                                     <  2:317802/70‑1 (MQ=255)                                                                                                   cTGCCCAAATGAGCCATCCTGACCCATTACCGCCTGCCATACCCCACCATAAAGCGAATCATCAAGCCACAAGCCGGGGCa                                                                          >  1:147687/1‑81 (MQ=255)                                                                                                   cTGCCCAAATGAGCCATCCTGACCCATTACCGCCTGCCATACCCCACCATAAAGCGAATCATCAAGCCACAAGCCGGGGCa                                                                          <  2:157508/81‑1 (MQ=255)                                                                                                   cTGCCCAAATGAGCCATCCTGACCCATTACCGCCTGCCATACCCCACCATAAAGCGAATCATCAAGCCACAAGCCGGGGCa                                                                          <  2:147687/81‑1 (MQ=255)                                                                                                   cTGCCCAAATGAGCCATCCTGACCCATTACCGCCTGCCATACCCCACCATAAAGCGAATCATCAAGCCACAAGCCGGGGCa                                                                          >  1:157508/1‑81 (MQ=255)                                                                                                   cTGCCCAAATGAGCCATCCTGACCCATTACCGCCTGCCATACCCCACCATAAAGCGAATCATCAAGCCACAAGCCGGGGCAAGAAACATAAATATCCCCGTCACTAGTTCATTACCGCCAGTTAACCACATGGTGAAGt                >  2:494573/1‑139 (MQ=255)                                                                                                     gCCCAAATGAGCCATCCTGACCCATTACCGCCTGCCATACCCCACCATAAAGCGAATCATCAAGCCACAAGCCGGGGCAAGAAACATAAATATCCCCGTCACTAGTTCATTACCGCCAGTTAACCACATGGTGAAGtaa              >  2:766000/1‑139 (MQ=255)                                                                                                        cAAATGAGCCATCCTGACCCATTACCGCCTGCCATACCCCACCATAAAGCGAATCATCAAGCCACAAGCCGGGGCAAGAAACATAAATATCCCCGTCACTAGTTCATTACCGCCAGTTAACCACATGGTg                    <  2:365622/130‑1 (MQ=255)                                                                                                        cAAATGAGCCATCCTGACCCATTACCGCCTGCCATACCCCACCATAAAGCGAATCATCAAGCCACAAGCCGGGGCAAGAAACATAAATATCCCCGTCACTAGTTCATTACCGCCAGTTAACCACATGGTg                    >  1:365622/1‑130 (MQ=255)                                                                                                                ccATCCTGACCCATTACCGCCTGCCATACCCCACCATAAAGCGAATCATCAAGCCACAAGCCGGGGCAAGAAACATAAATATCCCCGTCACTAGTTCATTACCGCCAGTTAACCACATGGTGAAGTAATAAATGAGTGc   >  1:603591/1‑139 (MQ=255)                                                                                                                ccATCCTGACCCATTACCGCCTGCCACACCCCACCATAAAGCGAATCATCAAGCCACAAGCCGGGGCAAGAAACATAAATATCCCCGTCACTAGTTCATTACCGCCAGTTAACCACATGGTGa                   <  2:547305/123‑1 (MQ=255)                                                                                                                ccATCCTGACCCATTACCGCCTGCCACACCCCACCATAAAGCGAATCATCAAGCCACAAGCCGGGGCAAGAAACATAAATATCCCCGTCACTAGTTCATTACCGCCAGTTAACCACATGGTGa                   >  1:547305/1‑123 (MQ=255)                                                                                                                 cATCCTGACCCATTACCGCCTGCCATACCCCACCATAAAGCGAATCATCAAGCCACAAGCCGGGGCAAGAAACATAAATATCCCCGTCACTAGTTCATTACCGCCAGTTAACCACATGGTGAAGTAATAAATGAGTGCa  <  2:581133/139‑1 (MQ=255)                                                                                                                             |                                                                                                                               CCCCCTTTTCCGACAAAGGCGAATTATCTTATTCATCCTGTGCGATAGCGTCTTTCACATGGCAGGTGAACAATCATGTCTGATGCTGATGTGAACACCTGCCCAAATGAGCCATCCTGACCCACTACCGCCTGCCATACCCCACCATAAAGCGAATCATCAAGCCACAAGCCGGGGCAAGAAACATAAATATCCCCGTCACTAGTTCATTACCGCCAGTTAACCACATGGTGAAGTAATAAATGAGTGCA  >  NC_000913/2494926‑2495176

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