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

GAAATGCCGCCCCCTTTTCCGACAAAGGCGAATTATCTTATTCATCCTGTGCGATAGCGTCTTTCACATGGCAGGTGAACAATCATGTCTGATGCTGATGTGAACACCTGCCCAAATGAGCCATCCTGACCCACTACCGCCTGCCATACCCCACCATAAAGCGAATCATCAAGCCACAAGCCGGGGCAAGAAACATAAATATCCCCGTCACTAGTTCATTACCGCCAGTTAACCACATGGTGAAGTAATAAATGAGTGCAG  >  NC_000913/2494917‑2495177                                                                                                                                      |                                                                                                                                gaaATGCCGCCCCCTTTTCCGACAAAGGCGAATCATCTTATTCATCCTGTGCGATAACGTCTTTCACATGGCAGGTGAACAATCATGTCTGATGCTGATGTGAACACCTGCCCAAATGAGCCATCCTGACCCATTACCg                                                                                                                            <  1:498111/139‑1 (MQ=255)                tttCCGACAAAGGCGAATCATCTTATTCATCCTGTGCGATAACGTCTTTCACATGGCAGGTGAACAATCATGTCTGATGCTGATGTGAACACCTGCCCAAATGAGCCATCCTGACCCATTACCGCCTGCCATACCccac                                                                                                             <  1:375642/139‑1 (MQ=255)                    cGACAAAGGCGAATCATCTTATTCATCCTGTGCGATAACGTCTTTCACATGGCAGGTGAACAATCATGTCTGATGCTGATGTGAACACCTGCCCAAATGAGCCATCCTGACCCATTACCGCCTGCCATACCCCACCATa                                                                                                         <  2:363785/139‑1 (MQ=255)                                atcaTCTTATTCATCCTGTGCGATAACGTCTTTCACATGGCAGGTGAACAATCATGTCTGATGCTGATGTGAACACCTGCCCAAATGAGCCATCCTGACCCATTACCGCCTGCCATACCCCACCATAAAGCGa                                                                                                   >  2:184036/4‑133 (MQ=255)                                atcaTCTTATTCATCCTGTGCGATAACGTCTTTCACATGGCAGGTGAACAATCATGTCTGATGCTGATGTGAACACCTGCCCAAATGAGCCATCCTGACCCATTACCGCCTGCCATACCCCACCATAAAGCGa                                                                                                   <  1:184036/130‑1 (MQ=255)                                              ccTGTGCGATAACGTCTTTCACATGGCAGGTGAACAATCATGTCTGATGCTGATGTGAACACCTGTCCAAATGAGCCATCCTGACCCATTACCGCCTGCCATACCCCACCATAAAGCGAATCATCAAGCCACAAGCCgg                                                                               >  2:599859/1‑139 (MQ=255)                                               cTGTGCGATAACGTCTTTCACATGGCAGGTGAACAATCATGTCTGATGCTGATGTGAACACCTGCCCAAATGAGCCATCCTGACCCATTACCGCCTGCCATACCCCACCATAAAGCGAATCATCAAGCCACAAGCCggg                                                                              >  1:302635/1‑139 (MQ=255)                                                  tgCGATAACGTCTTTCACATGGCAGGTGAACAATCATGTCTGATGCTGATGTGAACACCTGCCCAAATGAGCCATCCTGACCCATTACCGCCTGCCATACCCCACCATAAAGCGAATCATCAAGCCACAAGCCGGGGCa                                                                           >  1:354540/1‑139 (MQ=255)                                                  tgCGATAACGTCTTTCACATGGCAGGTGAACAATCATGTCTGATGCTGATGTGAACACCTGCCCAAATGAGCCATCCTGACCCATTACCGCCTGCCATACCCCACCATAAAGCGAATCATCAAGCCACAAGCCGGGGCa                                                                           >  1:530521/1‑139 (MQ=255)                                                   gCGATAACGTCTTTCACATGGCAGGTGAACAATCATGTCTGATGCTGATGTGAACACCTGCCCAAATGAGCCATCCTGACCCATTACCGCCTGCCATACCCCACCATAAAGCGAATCATCAAGCCACAAGCCGGGGCaa                                                                          >  1:1167602/1‑139 (MQ=255)                                                    cGATAACGTCTTTCACATGGCAGGTGAACAATCATGTCTGATGCTGATGTGAACACCTGCCCAAATGAGCCATCCTGACCCATTACCGCCTGCCATACCCCACCATAAAGCGAATCATCAAGCCACAAGCCGGGGCAAg                                                                         <  2:386566/139‑1 (MQ=255)                                                    cGATAACGTCTTTCACATGGCAGGTGAACAATCATGTCTGATGCTGATGTGAACACCTGCCCAAATGAGCCATCCTGACCCATTACCGCCTGCCATACCCCACCATAAAGCGAATCATCAAGCCACAAGCCGGGGCAAg                                                                         <  2:1167602/139‑1 (MQ=255)                                                     gATAACGTCTTTCACATGGCAGGTGAACAATCATGTCTGATGCTGATGTGAACACCTGCCCAAATGAGCCATCCTGACCCATTACCGCCTGCCATACCCCACCATAAAGCGAATCATCAAGCCACAAGCCGGGGCAAGa                                                                        <  1:750593/139‑1 (MQ=255)                                                       taacGTCTTTCACATGGCAGGTGAACAATCATGTCTGATGCTGATGTGAACACCTGCCCAAATGAGCCATCCTGACCCATTACCGCCTGCCATACCCCACCATAAAGCGAATCATCAAGCCACAAGCCGGGGCAAGaaa                                                                      <  1:142170/136‑1 (MQ=255)                                                          cGTCTTTCACATGGCAGGTGAACAATCATGTCTGATGCTGATGTGAACACCTGCCCAAATGAGCCATCCTGACCCATTACCGCCTGCCATACCCCACCATAAAGCGAATCATCAAGCCACAAGCCGGGGCAAGAAACAt                                                                   <  2:1079012/139‑1 (MQ=255)                                                                    aTGGCAGGTGAACAATCATGTCTGATGCTGATGTGAACACCTGCCCAAATGAGCCATCCTGACCCATTAcc                                                                                                                             <  1:1216982/71‑1 (MQ=255)                                                                    aTGGCAGGTGAACAATCATGTCTGATGCTGATGTGAACACCTGCCCAAATGAGCCATCCTGACCCATTAcc                                                                                                                             >  2:1216982/1‑71 (MQ=255)                                                                    aTGGCAGGTGAACAATCATGTCTGATGCTGATGTGAACACCTGCCCAAATGAGCCATCCTGACCCATTACCGCCTGCCATACCCCACCATAAAGCGAATCATCAAGCCACAAGCCGGGGCAAGAAACATAAATATcccc                                                         <  1:287356/139‑1 (MQ=255)                                                                      ggCAGGTGAACAATCATGTCTGATGCTGATGTGAACACCTGCCCAAATGAGCCATCCTGACCCATTACCGCCTGCCATACCCCACCATAAAGCGAATCATCAAGCCACAAGCCGGGGCAAGAAACATAAATATCCCCGt                                                       >  1:1217992/1‑139 (MQ=255)                                                                           gTGAACAATCATGTCTGATGCTGATGTGAACACCTGCCCAAATGAGCCATCCTGACCCATTACCGCCTGCCa                                                                                                                     >  1:919637/1‑72 (MQ=255)                                                                           gTGAACAATCATGTCTGATGCTGATGTGAACACCTGCCCAAATGAGCCATCCTGACCCATTACCGCCTGCCa                                                                                                                     <  2:919637/72‑1 (MQ=255)                                                                              aaCAATCATGTCTGATGCTGATGTGAACACCTGCCCAAATGAGCCATCCTGACCCATTACCGCCTGCCATACCCCACCATAAAGCGAATCATCAAGCCACAAGCCGGGGCAAGAAACATAAATATCCCCGTCACTAGtt                                               >  2:348653/1‑139 (MQ=255)                                                                                  aTCATGTCTGATGCTGATGTGAACACCTGCCCAAATGAGCCATCCTGACCCATTACCGCCTGCCATACCCCACCATAAAGCGAATCATCAAGCCACAAGCCGGGGCAAGAAACATAAATATCCCCGTCAc                                                    >  1:506039/1‑130 (MQ=255)                                                                                  aTCATGTCTGATGCTGATGTGAACACCTGCCCAAATGAGCCATCCTGACCCATTACCGCCTGCCATACCCCACCATAAAGCGAATCATCAAGCCACAAGCCGGGGCAAGAAACATAAATATCCCCGTCAc                                                    <  2:506039/130‑1 (MQ=255)                                                                                  aTCATGTCTGATGCTGATGTGAACACCTGCCCAAATGAGCCATCCTGACCCATTACCGCCTGCCATACCCCACCATAAAGCGAATCATCAAGCCACAAGCCGGGGCAAGAAACATAAATATCCCCGTCACTAGTTCAtt                                           <  1:499904/139‑1 (MQ=255)                                                                                     aTGTCTGATGCTGATGTGAACACCTGCCCAAATGAGCCATCCTGACCCATTACCGCCTGCCATACCTCACCATAAAGCGAATCATCAAGCCACAAGCCGGGGCAAGAAACATAAATATCCCCGTCACTAGTTCATTAcc                                        <  2:530521/139‑1 (MQ=255)                                                                                     aTGTCTGATGCTGATGTGAACACCTGCCCAAATGAGCCATCCTGACCCATTACCGCCTGCCATACCCCACCATAAAGCGAATCATCAAGCCACAAGCCGGGGCa                                                                           <  1:145932/104‑1 (MQ=255)                                                                                     aTGTCTGATGCTGATGTGAACACCTGCCCAAATGAGCCATCCTGACCCATTACCGCCTGCCATACCCCACCATAAAGCGAATCATCAAGCCACAAGCCGGGGCa                                                                           >  2:145932/1‑104 (MQ=255)                                                                                        tCTGATGCTGATGTGAACACCTGCCCAAATGAGCCATCCTGACCCATTACCGCCTGCCATACCCCACCATAAAGCGAATCATCAAGCCACAAGCCGGGGCAAGAAACATAAATATCCCCGTCACTAGTTCATTACCGcc                                     >  2:790053/1‑139 (MQ=255)                                                                                        tCTGATGCTGATGTGAACACCTGCCCAAATGAGCCATCCTGACCCATTACCGCCTGCCATACCCCACCATAAAGCGAATCATCAAGCCACAAGCCGGGGCAAGAAACATAAATATCCCCGTCACTAGTTCATTACCGcc                                     >  1:296683/1‑139 (MQ=255)                                                                                           gatgctgatgTGAACACCTGCCCAAATGAGCCATCCTGACCCATTACCGCCTGCCATACCCCACCATAAAGCGAATCATCAAGCCACAAGCCGGGGCAAGAAATATAAATATCCCCGTCACTAGTTCATTACCGCCAGt                                  >  1:1154985/1‑139 (MQ=255)                                                                                             tgctgatgTGAACACCTGCCCAAATGAGCCATCCTGACCCATTACCGCCTGCCATACCCCACCATAAAGCGAATCATCAAGCCACAAGCCGGGGCAAGAAACATAAATATCCCCGTCACTAGTTCATTACCGCCAGTTa                                >  1:937154/1‑139 (MQ=255)                                                                                                  atgTGAACACCTGCCCAAATGAGCCATCCTGACCCATTACCGCCTGCCATACCCCACCATAAAGCGAATCATCAAGCCACAAGCCGGGGCAAGAAACATAAATATcccc                                                         <  1:276660/109‑1 (MQ=255)                                                                                                  atgTGAACACCTGCCCAAATGAGCCATCCTGACCCATTACCGCCTGCCATACCCCACCATAAAGCGAATCATCAAGCCACAAGCCGGGGCAAGAAACATAAATATcccc                                                         >  2:276660/1‑109 (MQ=255)                                                                                                  atgTGAACACCTGCCCAAATGAGCCATCCTGACCCATTACCGCCTGCCATACCCCACCATAAAGCGAATCATCAAGCCACAAGCCGGGGCAAGAAACATAAATAT‑CCCGTCACTAGTTCATTACCGCCAGTTAACcaca                          <  2:354540/139‑1 (MQ=255)                                                                                                   tgtgAACACCTGCCCAAATGAGCCATCCTGACCCATTACCGc                                                                                                                           >  2:711479/1‑42 (MQ=255)                                                                                                   tgtgAACACCTGCCCAAATGAGCCATCCTGACCCATTACCGc                                                                                                                           <  1:711479/42‑1 (MQ=255)                                                                                                    gtgAACACCTGCCCAAATGAGCCATCCTGACCCATTACCGCCTGCCATACCCCACCATAAAGCGAATCATCAAGCCACAAGCCGGGGCAAGAAACATAAATATCCCCGTCACTAGTTCATTACCGCCAGTTAACCACAt                         <  2:1217992/139‑1 (MQ=255)                                                                                                       aaCACCTGCCCAAATGAGCCATCCTGACCCATTACCGCCTGCCATACCCCACCATAAAGCGAATCATCAAGCCACAAGCCGGGGCAAGAAACATAAATATCCCCGTCACTAGTTCATTACCGc                                      >  1:195890/1‑123 (MQ=255)                                                                                                       aaCACCTGCCCAAATGAGCCATCCTGACCCATTACCGCCTGCCATACCCCACCATAAAGCGAATCATCAAGCCACAAGCCGGGGCAAGAAACATAAATATCCCCGTCACTAGTTCATTACCGc                                      <  2:195890/123‑1 (MQ=255)                                                                                                       aaCACCTGCCCAAATGAGCCATCCTGACCCATTACCGCCTGCCATACCCCACCATAAAGCGAATCATCAAGCCACAAGCCGGGGCAAGAAACATAAATATCCCCGTCACTAGTTCATTACCGCCAGTTAACCACATGGt                      <  2:302635/139‑1 (MQ=255)                                                                                                         cacCTGCCCAAATGAGCCATCCTGACCCATTACCGCCTGCCATACCCCACCATAAAGCGAATCATCAAGCCACAAGCCGGGGCAAGAAATATAAATATCCCCGTCACTAGTTCATTACCGCCAGTTAACCACATGGTGa                    <  2:1154985/139‑1 (MQ=255)                                                                                                            cTGCCCAAATGAGCCATCCTGACCCATTACCGCCTGCCATACCCCACCATAAAGCGAATCATCAAGCCACAAGc                                                                                  >  1:821249/1‑74 (MQ=255)                                                                                                            cTGCCCAAATGAGCCATCCTGACCCATTACCGCCTGCCATACCCCACCATAAAGCGAATCATCAAGCCACAAGc                                                                                  <  2:821249/74‑1 (MQ=255)                                                                                                              gCCCAAATGAGCCATCCTGACCCATTACCGCCTGCCATACCCCACCATAAAGCGAATCATCAAGCCACAAGCCGGGGCAAGAAACATAAATATCCCCGTCACTAGTTCATTACCGCCAGTTAACCACATGGTGa                    >  2:900532/1‑134 (MQ=255)                                                                                                              gCCCAAATGAGCCATCCTGACCCATTACCGCCTGCCATACCCCACCATAAAGCGAATCATCAAGCCACAAGCCGGGGCAAGAAACATAAATATCCCCGTCACTAGTTCATTACCGCCAGTTAACCACATGGTGa                    <  1:900532/134‑1 (MQ=255)                                                                                                               cccAAATGAGCCATCCTGACCCATTACCGCCTGCCATACCCCACCATAAAGCGAATCATCAGGCCACAAGCCGGGGCa                                                                           >  2:524909/1‑78 (MQ=255)                                                                                                               cccAAATGAGCCATCCTGACCCATTACCGCCTGCCATACCCCACCATAAAGCGAATCATCAGGCCACAAGCCGGGGCa                                                                           <  1:524909/78‑1 (MQ=255)                                                                                                               cccAAATGAGCCATCCTGACCCATTACCGCCTGCCATACCCCACCATAAAGCGAATCATCAAGCCACAAGCCGGGGCa                                                                           <  1:202957/78‑1 (MQ=255)                                                                                                               cccAAATGAGCCATCCTGACCCATTACCGCCTGCCATACCCCACCATAAAGCGAATCATCAAGCCACAAGCCGGGGCa                                                                           >  2:202957/1‑78 (MQ=255)                                                                                                               cccAAATGAGCCATCCTGACCCATTACCGCCTGCCATACCCCACCATAAAGCGAATCATCAAGCCACAAGCCGGGGCAAGAAACATAAATATCCCCGTCACTAGTTCATTACCGCCAGTTAACCACATGGTGa                    <  2:444204/133‑1 (MQ=255)                                                                                                               cccAAATGAGCCATCCTGACCCATTACCGCCTGCCATACCCCACCATAAAGCGAATCATCAAGCCACAAGCCGGGGCAAGAAACATAAATATCCCCGTCACTAGTTCATTACCGCCAGTTAACCACATGGTGa                    >  1:444204/1‑133 (MQ=255)                                                                                                                         ccATCCTGACCCATTACCGCCTGCCATACCCCACCATAAAGCGAATCATCAAGCCACAAGCCGGGGCAAGAAACATAAATATCCCCGTCACTAGTTCATTAc                                         >  1:1296536/1‑102 (MQ=255)                                                                                                                         ccATCCTGACCCATTACCGCCTGCCATACCCCACCATAAAGCGAATCATCAAGCCACAAGCCGGGGCAAGAAACATAAATATCCCCGTCACTAGTTCATTAc                                         <  2:1296536/102‑1 (MQ=255)                                                                                                                         ccATCCTGACCCATTACCGCCTGCCATACCCCACCATAAAGCGAATCATCAAGCCACAAGCCGGGGCAAGAAACATAAATATCCCCGTCACTAGTTCATTACCGCCAGTTAACCACATGGTGa                    <  2:1250973/123‑1 (MQ=255)                                                                                                                         ccATCCTGACCCATTACCGCCTGCCATACCCCACCATAAAGCGAATCATCAAGCCACAAGCCGGGGCAAGAAACATAAATATCCCCGTCACTAGTTCATTACCGCCAGTTAACCACATGGTGa                    >  1:1250973/1‑123 (MQ=255)                                                                                                                          cATCCTGACCCATTACCGCCTGCCATACCCCACCATAAAGCGAATCATCAAGCCACAAGCCGGGGCAAGAAACATAAATATCCCCGTCACTAGTTCATTACCGCCAGTTAACcaca                          <  2:1142140/116‑1 (MQ=255)                                                                                                                          cATCCTGACCCATTACCGCCTGCCATACCCCACCATAAAGCGAATCATCAAGCCACAAGCCGGGGCAAGAAACATAAATATCCCCGTCACTAGTTCATTACCGCCAGTTAACcaca                          >  1:1142140/1‑116 (MQ=255)                                                                                                                           aTCCTGACCCATTACCGCCTGCCATACCCCACCATAAAGCGAATCATCAAGCCACAAGCCGGGGCAAGAAACATAAATATCCCCGTCACTAGTTCATTACCGCCAGTTAACCACATGGTGAAGTAATAAATGAGTGCAg  <  1:348653/139‑1 (MQ=255)                                                                                                                            tCCTGACCCATTACCGCCTTCCATACCCCACCATAAAGCGAATCATCAAGCCACAAGCCGGGGCa                                                                           >  1:862501/1‑65 (MQ=255)                                                                                                                            tCCTGACCCATTACCGCCTGCCATACCCCACCATAAAGCGAATCATCAAGCCACAAGCCGGGGCa                                                                           <  2:862501/65‑1 (MQ=255)                                                                                                                                 aCCCATTACCGCCTGCCATACCCCACCATAAAGCGAATCATCAAGCCACAAGCCGGGGCAAGaa                                                                       <  1:5924/64‑1 (MQ=255)                                                                                                                                 aCCCATTACCGCCTGCCATACCCCACCATAAAGCGAATCATCAAGCCACAAGCCGGGGCAAGaa                                                                       >  2:5924/1‑64 (MQ=255)                                                                                                                                      |                                                                                                                                GAAATGCCGCCCCCTTTTCCGACAAAGGCGAATTATCTTATTCATCCTGTGCGATAGCGTCTTTCACATGGCAGGTGAACAATCATGTCTGATGCTGATGTGAACACCTGCCCAAATGAGCCATCCTGACCCACTACCGCCTGCCATACCCCACCATAAAGCGAATCATCAAGCCACAAGCCGGGGCAAGAAACATAAATATCCCCGTCACTAGTTCATTACCGCCAGTTAACCACATGGTGAAGTAATAAATGAGTGCAG  >  NC_000913/2494917‑2495177

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