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%	169.0 / NA	53	*81* (TAG→TAA)	yfdY	DUF2545 family putative inner membrane protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (27/26);  total (27/26)

AATGCCGCCCCCTTTTCCGACAAAGGCGAATTATCTTATTCATCCTGTGCGATAGCGTCTTTCACATGGCAGGTGAACAATCATGTCTGATGCTGATGTGAACACCTGCCCAAATGAGCCATCCTGACCCACTACCGCCTGCCATACCCCACCATAAAGCGAATCATCAAGCCACAAGCCGGGGCAAGAAACATAAATATCCCCGTCACTAGTTCATTACCGCCAGTTAACCACATGGTGAAGTAATAAATGAGTG  >  NC_000913/2494919‑2495174                                                                                                                                    |                                                                                                                             aaTGCCGCCCCCTTTTCCGACAAAGGCGAATCATCTTATTCATCCTGTGCGATAACGTCTTTCACATGGCAGGTGAACAATCATGTCTGATGCTGATGTGAACACCTGCCCAAATGAGCCATCCTGACCCATTACCGcc                                                                                                                       <  2:780296/139‑1 (MQ=255)   tGCCGCCCCCTTTTCCGACAAAGGCGAATCATCTTATTCATCCTGTGCGATAACGTCTTTCACATGGCAGGTGAACAATCATGTCTGATGCTGATGTGAACACCTGCCCAAATGAGCCATCCTGACCCATTACCGCCTg                                                                                                                     >  2:554499/1‑139 (MQ=255)         ccccTTTTCCGACAAAGGCGAATCATCTTATTCATCCTGTGCGATAACGTCTTTCACATGGCAGGTGAACAATCATGTCTGATGCTGATGTGAACACCTGCCCAAATGAGCCATCCTGACCCATTACCGCCTGCCATAc                                                                                                               >  1:1042750/1‑139 (MQ=255)               ttCCGACAAAGGCGAATCATCTTATTCATCCTGTGCGATAACGTCTTTCACATGGCAGGTGAACAATCATGTCTGATGCTGATGTGAACACCTGCCCAAATGAGCCATCCTGACCCATTACCGCCTGCCATACCccacc                                                                                                         <  1:128721/139‑1 (MQ=255)                              atcaTCTTATTCATCCTGTGCGATAACGTCTTTCACATGGCAGGTGAACAATCATGTCTGATGCTGATGTGAACACCTGCCCAAATGAGCCATCCTGACCCATTACCGCCTGCCATACCCCACCATAAAGCGAatcatc                                                                                          <  2:808583/136‑1 (MQ=255)                                    ttATTCATCCTGTGCGATTACGTCTTTCACATGGCAGGTGAACAATCATGTCTGATGCTGATGTGAACACCTGCCCAAATGAGCCATCCTGACCCATTACCGCCTGCCATACCCCACCATAAAGCGAATCATCAAGCca                                                                                    <  1:708797/139‑1 (MQ=255)                                            ccTGTGCGATAACGTCTTTCACATGGCAGGTGAACAATCATGTCTGATGCTGATGTGAACACCTGCCCAAATGAGCCATCCTGACCCATTACCGCCTGCCATACCCCACCATAAAGCGAATCATCAAGCCACAAGc                                                                               >  1:14817/1‑136 (MQ=255)                                            ccTGTGCGATAACGTCTTTCACATGGCAGGTGAACAATCATGTCTGATGCTGATGTGAACACCTGCCCAAATGAGCCATCCTGACCCATTACCGCCTGCCATACCCCACCATAAAGCGAATCATCAAGCCACAAGc                                                                               <  2:14817/136‑1 (MQ=255)                                            ccTGTGCGATAACGTCTTTCACATGGCAGGTGAACAATCATGTCTGATGCTGATGTGAACACCTGCCCAAATGAGCCATCCTGACCCATTACCGCCTGCCATACCCCACCATAAAGCGAATCATCAAGCCACAAGCCgg                                                                            >  2:173796/1‑139 (MQ=255)                                            ccTGTGCGATAACGTCTTTCACATGGCAGGTGAACAATCATGTCTGATGCTGATGTGAACACCTGCCCAAATGAGCCATCCTGACCCATTACCGCCTGCCATACCCCACCATAAAGCGAATCATCAAGCCACAAGCCgg                                                                            >  1:972211/1‑139 (MQ=255)                                             cTGTGCGATAACGTCTTTCACATGGCAGGTGAACAATCATGTCTGATGCTGATGTGAACACCTGCCCAAATGAGCCATCCTGACCCATTACCGCCTGCCATACCCCACCATAAAGCGAATCATCAAGCCACAAGCCggg                                                                           <  1:183734/139‑1 (MQ=255)                                               gtgCGATAACGTCTTTCACATGGCAGGTGAACAATCATGTCTGATGCTGATGTGAACACCTGCCCAAATGAGCCATCCTGACCCATTACCGCCTGCCATACCCCACCATAAAGCGAATCATCAAGCCACAAGCCGGGGc                                                                         >  2:671546/1‑139 (MQ=255)                                                tgCGATAACGTCTTTCACATGGCAGGTGAACAATCATGTCTGATGCTGATGTGAACACCTGCCCAAATGAGCCATCCTGACCCATTACCGCCTGCCATACCCCACCATAAAGCGAATCATCAAGCCACAAGCCGGGGCa                                                                        >  1:262154/1‑139 (MQ=255)                                                  cGATAACGTCTTTCACATGGCAGGTGAACAATCATGTCTGATGCTGATGTGAACACCTGCCCAAATGAGCCATCCTGACCCATTACCGCCTGCCATACCCCACCATAAAGCGAATCATCAAGCCACAAGCCGGGGCAAg                                                                      <  2:875044/139‑1 (MQ=255)                                                  cGATAACGTCTTTCACATGGCAGGTGAACAATCATGTCTGATGCTGATGTGAACACCTGCCCAAATGAGCCATCCTGACCCATTACCGCCTGCCATACCCCACCATAAAGCGAATCATCAAGCCACAAGCCGGGGCAAg                                                                      >  1:856659/1‑139 (MQ=255)                                                     taacGTCTTTCACATGGCAGGTGAACAATCATGTCTGATGCTGATGTGAACACCTGCCCAAATGAGCCATCCTGACCCATTACCg                                                                                                                         >  2:528900/4‑85 (MQ=255)                                                     taacGTCTTTCACATGGCAGGTGAACAATCATGTCTGATGCTGATGTGAACACCTGCCCAAATGAGCCATCCTGACCCATTACCg                                                                                                                         <  1:528900/82‑1 (MQ=255)                                                          tCTTTCACATGGCAGGTGAACAATCATGTCTGATGCTGATGTGAACACCTGCCCAAATGAGCCATCCTGACCCATTACCGCCTGCCATACCCCACCATAAAGCGAATCATCAAGCCACAAGCCGGGGCAAGAAACATaa                                                              >  1:504932/1‑139 (MQ=255)                                                             ttCACATGGCAGGTGAACAATCATGTCTGATGCTGATGTGAACACCTGCCCAAATGAGCCATCCTGACCCATTACCGCCTGCCATACCCCACCATAAAGCGAATCATCAATCCACAAGCCGGGGCAAGAAACATAAata                                                           >  1:633948/1‑139 (MQ=255)                                                                  aTGGCAGGTGAACAATCATGTCTGATGCTGATGTGAACACCTGCCCAAATGAGCCATCCTGACCCATTACCGCCTGCCATACCCCACCATAAAGCGAATCATCAAGCCACAAGCCGGGGCAAGAAACATAAATATcccc                                                      <  2:951490/139‑1 (MQ=255)                                                                       aGGTGAACAATCATGTCTGATGCTGATGTGAACACCTGCCCAAATGAGCCATCCTGACCCATTACCGCCTGCCATACCCCACCATAAAGCGAATCATCAAGCCACAAGCCGGGGCAAGAAACATAAATATCCCCGTCAc                                                 <  1:91209/139‑1 (MQ=255)                                                                        tgTGAACAATCATGTCTGATGCTGATGTGAACACCTGCGCAAATGAGCCATCCTGACCCATTACCGCCTGCCATACCCCACCATAAAGCGAATCATCAAGCCACAAGCCGGGGTAAGTAACATAAATATCCCCGTCACt                                                >  2:253564/2‑139 (MQ=255)                                                                              cAATCATGTCTGATGCTGATGTGAACACCTGCCCAAATGAGCCATCCTGACCCATTACCGCCTGCCATACCCCACCATAAAGCGAATCATCAAGCCACAAGCCGGGGCAAGAAACATAAATATCCCCGTCACTAGTTCa                                          >  1:348175/1‑139 (MQ=255)                                                                                aTCATGTCTGATGCTGATGTGAACACCTGCCCAAATGAGCCATCCTGACCCATTACCGCCTGCCATACCCCACCATAAAGCGAATCATCAAGCCACAAGCCGGGGCAAGAAACATAAATATCCCCGTCACTAGTTCAtt                                        >  1:191476/1‑139 (MQ=255)                                                                                   aTGTCTGATGCTGATGTGAACACCTGCCCAAATGAGCCATCCTGACCCATTACCGCCTGCCATACCCCACCATAAAGCGa                                                                                                <  1:177498/80‑1 (MQ=255)                                                                                   aTGTCTGATGCTGATGTGAACACCTGCCCAAATGAGCCATCCTGACCCATTACCGCCTGCCATACCCCACCATAAAGCGa                                                                                                >  2:177498/1‑80 (MQ=255)                                                                                     gTCTGATGCTGATGTGAACACCTGCCCAAATGAGCCATCCTGACCCATTACCGCCTGCCATACCCCACCATAAAGCGAATCATCAAGCCACAAGCCGGGGCAAGAAACATAAATATCCCCGTCACTAGTTCATTACCGc                                   <  2:882800/139‑1 (MQ=255)                                                                                         gatgctgatgTGAACACCTGCCCAAATGAGCCCTCCTGACCCATTACCGCCTGCCATACCCCAACATTAAGCGAATCATCAAGCCACAAGCCGGGGCAAGAAACATAAATATCCCCGTCACTAGTTCATTACCGCCAGt                               <  2:183898/139‑1 (MQ=255)                                                                                         gatgctgatgTGAACACCTGCCCAAATGAGCCATCCTGACCCATTACCGCCTGCCATACCCCACCATAAAGCGAATCATCAAGCCACAAGCCGGGGCAAGAAACATAAATATCCCCGTCACTAGTTCATTACCGCCAGt                               <  2:504932/139‑1 (MQ=255)                                                                                         gatgctgatgTGAACACCTGCCCAAATGAGCCATCCTGACCCATTACCGCCTGCCATACCCCACCATAAAGCGAATCATCAAGCCACAAGCCGGGGCAAGAAACATAAATATCCCCGTCACTAGTTCATTACCGCCAGt                               <  1:554499/139‑1 (MQ=255)                                                                                           tgctgatgTGAACACCTGCCCAAATGAGCCATCCTGACCCATTACCGCCTGCCATACCCCACCATAAAGCGAATCATCAAGCCACAAGCCGGGGCAAGAAACATAAATATCCCCGTCACTAGTTCATTACCGCCAGTTa                             >  1:749272/1‑139 (MQ=255)                                                                                             ctgatgTGAACACCTGCCCAAATGAGCCATCCTGACCCATTACCGCCTGCCATACCCCACCATAAAGCGAATCATCAAGCCACAAGCCGGGGCAAGAAACATAAATATCCCCGTCACTAGTTCATTACCGCCAGTTAAc                           >  1:765692/1‑139 (MQ=255)                                                                                              tgatgTGAACACCTGCCCAAATGAGCCATCCTGACCCATTACCGCCTGCCATACCCCACCATAAAGCGAATCATCAAGCCACAAGCCGGGGCAAGAAACATAAATATCCCCGTCACTAGTTCATTACCGCCAGTTAAcc                          <  1:31730/139‑1 (MQ=255)                                                                                               gatgTGAACACCTGCCCAAATGAGCCATCCTGACCCATTACCGCCTGCCATACCCCACCATAAAGCGAATCATCAAGCCACAAGCCGGGGCAAGAAACATAAATATCCCCGTCACTAGTTCATTACCGCCAGTTAACca                         <  1:197044/139‑1 (MQ=255)                                                                                                atgTGAACACCTGCCCAAATGAGCCATCCTGACCCATTACCGCCTGCCATACCCCACCATAAAGCGAATCATCAAGCCACAAGCCGGGGCAAGAAACATAAATATCCCCGTCACTAGTTCATTACCGCCAGtt                              >  1:1335/1‑133 (MQ=255)                                                                                                atgTGAACACCTGCCCAAATGAGCCATCCTGACCCATTACCGCCTGCCATACCCCACCATAAAGCGAATCATCAAGCCACAAGCCGGGGCAAGAAACATAAATATCCCCGTCACTAGTTCATTACCGCCAGtt                              <  2:1335/133‑1 (MQ=255)                                                                                                  gtgAACACCTGCCCAAATGAGCCATCCTGACCCATTACCGCCTGCCATACCCCACCATAAAGCGAATCATCAAGCCACAAGCCGGGGCAAGAAACATAAATATCCCCGTCACTAGTTCATTACCGCCAGTTAACCACAt                      <  2:749749/139‑1 (MQ=255)                                                                                                     aaCACCTGCCCAAATGAGCCATCCTGACCCATTACCGCCTGCCATACCCCACCATAAAGCGAATCATCAAGCCACAAGCCGGGGCAAGAAACATAAATATCCCCGTCACTAGTTCATTACCGCCAGTTAACCACATGGt                   <  1:173796/139‑1 (MQ=255)                                                                                                       caccgtactaaaTGAGTCAGCCTGACCCATTACCGCCTTCCATACCCCCGCACGAAGCGAATCATCAAGCCAAAAGCCGGGTCAAGAACCTCAAAGATCACCGTGGCTATTTCAATACCGCCAGTTAAACACATGTTGa                 <  2:79823/130‑1 (MQ=255)                                                                                                          cTGCCCAAATGAGCCATCCTGACCCATTACCGCCTGCCATACCCCACCATAAAGCGAATCATCAAGCCACAAGCCGGGGCAAGAAACATAAATATCCCCGTCACTAGTTCATTACCGCCAGTTAACCACATGGTGAAGt              <  2:348175/139‑1 (MQ=255)                                                                                                            gCCCAAATGAGCCATCCTGACCCATTACCGCCTGCCATACCCCACCATaa                                                                                                     <  2:902133/50‑1 (MQ=255)                                                                                                            gCCCAAATGAGCCATCCTGACCCATTACCGCCTGCCATACCCCACCATaa                                                                                                     >  1:902133/1‑50 (MQ=255)                                                                                                            gCCCAAATGAGCCATCCTGACCCATTACCGCCTGCCATACCCCACCATAAAGCGAATCATCAAGCCACAAGCCGGGGCAAGAAACATAAATATCCCCGTCACTAGTTCa                                          <  2:722196/109‑1 (MQ=255)                                                                                                            gCCCAAATGAGCCATCCTGACCCATTACCGCCTGCCATACCCCACCATAAAGCGAATCATCAAGCCACAAGCCGGGGCAAGAAACATAAATATCCCCGTCACTAGTTCa                                          >  1:722196/1‑109 (MQ=255)                                                                                                            gCCCAAATGAGCCATCCTGACCCATTACCGCCTGCCATACCCCACCATAAAGCGAATCATCAAGCCACAAGCCGGGGCAAGAAACATAAATATCCCCGTCACTAGTTCATTACCGCCAGTTAACCACATGGTGAAGtaa            >  2:1366/1‑139 (MQ=255)                                                                                                            gCCCAAATGAGCCATCCTGACCCATTACCGCCTGCCATACCCCACCATAAAGCGAATCATCAAGCCACAAGCCGGGGCAAGAAACATAAATATCCCCGTCACTAGTTCATTACCGACAGTTAACCACATGGTGAAGtaa            >  2:557007/1‑139 (MQ=255)                                                                                                             cccAAATGAGCCATCCTGACCCATTACCGCCTGCCATACCCCACCATAAAGCGAATCATCAAGCCACAAGCCGGGGCAAGAAACATAAATATcccc                                                      >  1:47007/1‑96 (MQ=255)                                                                                                             cccAAATGAGCCATCCTGACCCATTACCGCCTGCCATACCCCACCATAAAGCGAATCATCAAGCCACAAGCCGGGGCAAGAAACATAAATATcccc                                                      <  2:47007/96‑1 (MQ=255)                                                                                                             cccAAATGAGCCATCCTGACCCATTACCGCCTGCCATACCCCACCATAAAGCGAATCATCAAGCCACAAGCCGGGGCAAGAAACATAAATATCCCCGTCACTAGTTCATTACCGCCAGTTAACCACATGGTGAAGtaat           >  1:957353/1‑139 (MQ=255)                                                                                                               cAAATGAGCCATCCTGACCCATTACCGCCTGCCATACCCCACCATAGAGCGAATCATCAAGCCACAAGCCGGGGCAAGAAACATAAATATCCCCGTCACTAGTTCATTACCGCCAGTTAACCACATGGTGAAGtaataa         >  2:738215/1‑139 (MQ=255)                                                                                                                aaaTGAGCCATCCTGACCCATTACCGCCTGCCATACCCCACCATAGAGCGAATCATCAAGCCACAAGCCGGGGCAAGAAACATAAATATCCCCGTCACTAGTTCATTACCGCCAGTTAACCACATGGTGAAGTAATaaa        <  1:738215/139‑1 (MQ=255)                                                                                                                  aTGAGCCATCCTGACCCATTACCGCCTGCCATACCCCACCATAAAGCGAATCATCAAGTCACAAGCCGGGGCAAGAAACATAAATATCCCCGTCACTAGTTCATTACCGCCAGTTAACCACATGGTGAAGTAATAAATg      >  2:1026032/1‑139 (MQ=255)                                                                                                                    gAGCCATCCTGACCCATTACCGCCTGCCATACCCCACCATAAAGCGAATCATCAAGCCACAAGCCGGGGCAAGAAACATAAATATCCCCGTCACTAGTTCATTACCGCCAGTTAACCACATGGTGAAGTAATAAATGAg    >  2:572188/1‑139 (MQ=255)                                                                                                                      gCCATCCTGACCCATTACCGCCTGCCATACCCCACCATAAAGCGAATCATCAAGCCACAAGCCGGGGCAAGAAACATAAATATCCCCGTCACTAGTTCATTACCGCCAGTTAACCACATGGTGAAGTAATAAATGAGTg  <  2:749272/139‑1 (MQ=255)                                                                                                                                    |                                                                                                                             AATGCCGCCCCCTTTTCCGACAAAGGCGAATTATCTTATTCATCCTGTGCGATAGCGTCTTTCACATGGCAGGTGAACAATCATGTCTGATGCTGATGTGAACACCTGCCCAAATGAGCCATCCTGACCCACTACCGCCTGCCATACCCCACCATAAAGCGAATCATCAAGCCACAAGCCGGGGCAAGAAACATAAATATCCCCGTCACTAGTTCATTACCGCCAGTTAACCACATGGTGAAGTAATAAATGAGTG  >  NC_000913/2494919‑2495174

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