Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,767,179	G→A	*281* (TAG→TAA)	yibA →	putative immunity protein for polymorphic toxin RhsA, HEAT‑domain protein, lethality reduction protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,767,179	0	G	A	100.0%	184.6 / NA	58	*281* (TAG→TAA)	yibA	putative immunity protein for polymorphic toxin RhsA, HEAT‑domain protein, lethality reduction protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (34/24);  total (34/24)

TATGATGATATCATTGAAGCTGCGGGTGAATTAGGTGATAAAACGCTACTTCCTGTTTTAGATACTATGTTGTACAAGTTTGATGACAATGAAATTATAACTTCCGCTATTGATAAGCTGAAGCGTTCATAGCGTGAGTTGCCTATGCACAGTGGGGG‑ATTCCCGCCGGCACGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGACCG  >  NC_000913/3767048‑3767311                                                                                                                                    |                                                                                                                                      tATGATGATATCATTGAAGCTGCGGGTGAATTAGGTGATAAAACGCTACTTCCTGTTTTAGATACTATGTTGTACAAGTTTGATGACAATGAAATTATAACTTCCGCTATTGATAAGCTGAAGCGTTCATAACGTGAGt                                                                                                                                <  1:984434/139‑1 (MQ=255) tATGATGATATCATTGAAGCTGCGGGTGAATTAGGTGATAAAACGCTACTTCCTGTTTTAGATACTATGTTGTACAAGTTTGATGACAATGAAATTATAACTTCCGCTATTGATAAGCTGAAGCGTTCATAACGTGAGt                                                                                                                                <  1:1052891/139‑1 (MQ=255)   tgttgatATCATTGAAGCTGTGTGTGAATTAGGTGATAAAACGCTACTTCCTGTTGTAGGTACTATGTTGTACAAGTTTGATGACAATGAAATTATAACTTCCGCTATTGATAAGCTGAAGCGTTCATAACGTGAGTTg                                                                                                                              <  2:976497/136‑1 (MQ=255)    gatgatATCATTGAAGCTGCGGGTGAATTAGGTGATAAAACGCTACTTCCTGTTTTAGATACTATGTTGTACAAGTTTGATGACAATGAAATTATAACTTCCGCTATTGATAAGCTGAAGCGTTCATAACGTGAGTTGc                                                                                                                             <  1:809273/139‑1 (MQ=255)         tatCATTGAAGCTGCGGGTGAATTAGGTGATAAAACGCTACTTCCTGTTTTAGATACTATGTTGTACAAGTTTGATGACAATGAAATTATAACTTCCGCTATTGATAAGCTGAAGCGTTCATAACGTGAGTTGCCTATg                                                                                                                        >  1:26924/1‑139 (MQ=255)         tatCATTGAAGCTGCGGGTGAATTAGGTGATAAAACGCTACTTCCTGTTTTAGATACTATGTTGTACAAGTTTGATGACAATGAAATTATAACTTCCGCTATTGATAAGCTGAAGCGTTCATAACGTGAGTTGCCTATg                                                                                                                        >  2:123584/1‑139 (MQ=255)         tatCATTGAAGCTGCGGGTGAATTAGGTGATAAAACGCTACTTCCTGTTTTAGATACTATGTTGTACAAGTTTGATGACAATGAAATTATAACTTCCGCTATTGATAAGCTGAAGCGTTCATAACGTGAGTTGCCTATg                                                                                                                        >  1:443150/1‑139 (MQ=255)           tCATTGAAGCTGCGGGTGAATTAGGTGATAAAACGCTACTTCCTGTTTTAGATACTATGTTGTACAAGTTTGATGACAATGAAATTATAACTTCCGCTATTGATAAGCTGAAGCGTTCATAACGTGAGTTGCCTATGca                                                                                                                      <  1:203809/139‑1 (MQ=255)           tCATTGAAGCTGCGGGTGAATTAGGTGATAAAACGCTACTTCCTGTTTTAGATACTATGTTGTACAAGTTTGATGACAATGAAATTATAACTTCCGCTATTGATAAGCTGAAGCGTTCATAACGTGAGTTGCCTATGca                                                                                                                      <  1:919813/139‑1 (MQ=255)            cATTGAAGCTGCGGGTGAATTAGGTGATAAAACGCTACTTCCTGTTTTAGATACTATGTTGTACAAGTTTGATGACAATGAAATTATAACTTCCGCTATTGATAAGCTGAAGCGTTCATAACGTGAGTTGCCTATGcac                                                                                                                     <  2:1027505/139‑1 (MQ=255)                 aaGCTGCGGGTGAATTAGGTGATAAAACGCTACTTCCTGTTTTAGATACTATGTTGTACAAGTTTGATGACAATGAAATTATAACTTCCGCTATTGATAAGCTGAAGCGTTCATAACGTGAGTTGCCTATGCACAGTgg                                                                                                                >  2:693059/1‑139 (MQ=255)                   gCTGCGGGTGAATTAGGTGATAAAACGCTACTTCCTGTTTTAGATACTATGTTGTACAAGTTTGATGACAATGAAATTATAACTTCCGCTATTGATAAGCTGAAGCGTTCATAACGTGAGTTGCCTATGCACAGTgggg                                                                                                              >  2:144900/1‑139 (MQ=255)                   gCTGCGGGTGAATTAGGTGATAAAACGCTACTTCCTGTTTTAGATACTATGTTGTACAAGTTTGATGACAATGAAATTATAACTTCCGCTATTGATAAGCTGAAGCGTTCATAACGTGAGTTGCCTATGCACAGTgggg                                                                                                              >  2:776160/1‑139 (MQ=255)                    cTGCGGGTGAATTAGGTGATAAAACGCTACTTCCTGTTTTAGATACTATGTTGTACAAGTTTGATGACAATGAAATTATAACTTCCGCTATTGATAAGCTGAAGCGTTCATAACGTGAGTTGCCTATGCACAGTggggg                                                                                                             >  2:863757/1‑139 (MQ=255)                    cTGCGGGTGAATTAGGTGATAAAACGCTACTTCCTGTTTTAGATACTATGTTGTACAAGTTTGATGACAATGAAATTATAACTTCCGCTATTGATAAGCTGAAGCGTTCATAACGTGAGTTGCCTATGCACAGTggggg                                                                                                             >  1:659178/1‑139 (MQ=255)                       cGGGTGAATTAGGTGATAAAACGCTACTTCCTGTTTTAGATACTATGTTGTACAAGTTTGATGACAATGAAATTATAACTTCCGCTATTGATAAGCTGAAGCGTTCATAACGTGAGTTGc                                                                                                                             <  1:489114/120‑1 (MQ=255)                       cGGGTGAATTAGGTGATAAAACGCTACTTCCTGTTTTAGATACTATGTTGTACAAGTTTGATGACAATGAAATTATAACTTCCGCTATTGATAAGCTGAAGCGTTCATAACGTGAGTTGc                                                                                                                             >  2:1045658/1‑120 (MQ=255)                       cGGGTGAATTAGGTGATAAAACGCTACTTCCTGTTTTAGATACTATGTTGTACAAGTTTGATGACAATGAAATTATAACTTCCGCTATTGATAAGCTGAAGCGTTCATAACGTGAGTTGc                                                                                                                             <  1:1045658/120‑1 (MQ=255)                       cGGGTGAATTAGGTGATAAAACGCTACTTCCTGTTTTAGATACTATGTTGTACAAGTTTGATGACAATGAAATTATAACTTCCGCTATTGATAAGCTGAAGCGTTCATAACGTGAGTTGc                                                                                                                             >  2:489114/1‑120 (MQ=255)                            gAATTAGGTGATAAAACGCTACTTCCTGTTTTAGATACTATGTTGTACAAGTTTGATGACAATGAAATTATAACTTCCGCTATTGATAAGCTGAAGCGTTCATAACGTGAGTTGCCTATGCACAGTGGGGG‑ATTCccgc                                                                                                    <  1:693059/139‑1 (MQ=255)                               ttAGGTGATAAAACGCTACTTCCTGTTTTAGATACTATGTTGTACAAGTTTGATGACAATGAAATTATAACTTCCGCTATTGATAAGCTGAAGCGTTCATAACGTGAGTTGCCTATGCACAGTGGGGG‑ATTCCCGCCgg                                                                                                 <  2:808617/139‑1 (MQ=255)                                  ggTGATAAAACGCTACTTCCTGTTTTAGATACTATGTTGTACAAGTTTGATGACAATGAAATTATAACTTCCGCTATTGATAAGCTGAAGCGTTCATAACGTGAGTTGCCTATGCACAGTGGGGG‑ATTCCCGCCGGCAc                                                                                              <  2:443150/139‑1 (MQ=255)                                     gATAAAACGCTACTTCCTGTTTTAGATACTATGTTGTACAAGTTTGATGACAATGAAATTATAACTTCCGCTATTGATAAGCTGAAGCGTTCATAACGTGAGTTGCCTATGCACAGTGGGGG‑ATTCCCGCCGGCACGGt                                                                                           >  2:335508/1‑139 (MQ=255)                                                aCTTCCTGTTTTAGATACTATGTTGTACAAGTTTGATGACAATGAAATTATAACTTCCGCTATTGATAAGCTGAAGCGTTCATAACGTGAGTTGCCTATGCACAGTGGGGG‑ATTCCCGCCGGCACGGTGCAGATGCTCg                                                                                <  1:1040980/139‑1 (MQ=255)                                                 cTTCCTGTTTTAGATACTATGTTGTACAAGTTTGATGACAATGAAATTATAACTTCCGCTATTGATAAGCTGAAGCGTTCATAACGTGAGTTGCCTATGCACAGTGGGGG‑ATTCCCGCCGGCACGGTGCAGATGCTCGa                                                                               >  1:168836/1‑139 (MQ=255)                                                               tACTATGTTGTACAAGTTTGATGACAATGAAATTATAACTTCCGCTATTGATAAGCTGAAGCGTTCATAACGTGAGTTGCCTATGCACAGTGGGGG‑ATTCCCGCCGGCACGGTGCAGATGCTCGACCGGCTGGAAAGTg                                                                 >  2:847915/1‑139 (MQ=255)                                                                    tgttgtACAAGTTTGATGACAATGAAATTATAACTTCCGCTATTGATAAGCTGAAGCGTTCATAACGTGAGt                                                                                                                                <  1:711459/72‑1 (MQ=255)                                                                    tgttgtACAAGTTTGATGACAATGAAATTATAACTTCCGCTATTGATAAGCTGAAGCGTTCATAACGTGAGt                                                                                                                                >  2:711459/1‑72 (MQ=255)                                                                        gttcAAATTTGATGACTTTGAAATTATAACTTCCGTTATTGATAAGCTGAAGCGTTCATAACGCCAGATGCATATGCACAGTGTTAG‑ACTCCCCCCGGCACGCCGCAGATGCTCGACCGGTTGGAAAGTGAAATCCTgg                                                        >  2:423637/4‑139 (MQ=255)                                                                         tACAAGTTTGATGACAATGAAATTATAACTTCCGCTATTGATAAGCTGAAGCGTTCATAACGTGAGTTGCCTATGCACAGTGGGGG‑ATTCCCGCCGGCACGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTGGc                                                       >  2:781557/1‑139 (MQ=255)                                                                         tACAAGTTTGATGACAATGAAATTATAACTTCCGCTATTGATAAGCTGAAGCGTTCATAACGTGAGTTGCCTATGCACAGTGGGGG‑ATTCCCGCCGGCACGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTGGc                                                       >  2:798909/1‑139 (MQ=255)                                                                         tACAAGTTTGATGACAATGAAATTATAACTTCCGCTATTGATAAGCTGAAGCGTTCATAACGTGAGTTGCCTATGCACAGTGGGGG‑ATTCCCGCCGGCACGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTGGc                                                       >  1:100098/1‑139 (MQ=255)                                                                          aCAAGTTTGATGACAATGAAATTATAACTTCCGCTATTGATAAGCTGAAGCGTTCATAACGTGAGTTGCCTATGCACAGTGGGGG‑ATTCCCGCCGGCACGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTGGCt                                                      <  2:100098/139‑1 (MQ=255)                                                                          aCAAGTTTGATGACAATGAAATTATAACTTCCGCTATTGATAAGCTGAAGCGTTCATAACGTGAGTTGCCTATGCACAGTGGGGG‑ATTCCCGCCGGCACGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTGGCt                                                      <  1:776160/139‑1 (MQ=255)                                                                            aaGTTTGATGACAATGAAATTATAACTTCCGCTATTGATAAGCTGAAGCGTTCATAACGTGAGt                                                                                                                                <  1:632391/64‑1 (MQ=255)                                                                            aaGTTTGATGACAATGAAATTATAACTTCCGCTATTGATAAGCTGAAGCGTTCATAACGTGAGt                                                                                                                                >  2:632391/1‑64 (MQ=255)                                                                            aaGTTTGATGACAATGAAATTATAACTTCCGCTATTGATAAGCTGAAGCGTTCATAACGTGAGTTGCCTATGCACAGTGGGGG‑ATTCCCGCCGGCACGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGa                                                    >  1:213346/1‑139 (MQ=255)                                                                              gTTTGATGACAATGAAATTATAACTTCCGCTATTGATAAGCTGAAGCGTTCATAACGTGAGTTGCCTATGCACAGTGGGGG‑ATTCCCGCCGGCACGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGAcc                                                  >  2:317476/1‑139 (MQ=255)                                                                                ttGATGACAATGAAATTATAACTTCCGCTGTTGATAAGCTGAAGCGTTCATAACGTGAGTTGCCTATGCACAGTGGGGG‑ATTCCCGCCGGCACGGTGCGGATGCTTGACCGGCTGGAAAGTGAAATCCTGGCTGACAgg                                                >  1:585741/1‑139 (MQ=38)                                                                                 tgatgaCAATGAAATTATAACTTCCGCTATTGATAAGCTGAAGCGTTCATAACGTGAGTTGCCTATGCACAGTGGGGG‑ATTCCCGCCGGCACGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCggg                                               >  2:373075/1‑139 (MQ=255)                                                                                 tgatgaCAATGAAATTATAACTTCCGCTATTGATAAGCTGAAGCGTTCATAACGTGAGTTGCCTATGCACAGTGGGGG‑ATTCCCGCCGGCACGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCggg                                               >  1:970536/1‑139 (MQ=255)                                                                                   atgaCAATGAAATTATAACTTCCGCTATTGATAAGCTGAAGCGTTCATAACGTGAGTTGc                                                                                                                             >  2:412725/1‑60 (MQ=255)                                                                                   atgaCAATGAAATTATAACTTCCGCTATTGATAAGCTGAAGCGTTCATAACGTGAGTTGc                                                                                                                             <  1:412725/60‑1 (MQ=255)                                                                                   atgaCAATGAAATTATAACTTCCGCTATTGATAAGCTGAAGCGTTCATAACGTGAGTTGCCTATGCACAGTGGGGG‑ATTCCCGCCGGCACGGTGCAGATGCTCGACCGGCTGGAAGGTGAAATCCTGGCTGACCGGgtg                                             >  2:239852/1‑139 (MQ=255)                                                                                   atgaCAATGAAATTATAACTTCCGCTATTGATAAGCTGAAGCGTTCATAACGTGAGTTGCCTATGCACAGTGGGGG‑ATTCCCGCCGGCACGGTGCAGATGCTCGACCGGCGGGAAAGTGAAATCCTGGCTGACCGGgtg                                             >  1:294795/1‑139 (MQ=255)                                                                                        aaTGAAATTATAACTTCCGCTATTGATAAGCTGAAGCGTTCATAACGTGAGTTGCCTATGCACAGTGGGGG‑ATTCCCGCCGGCACGGTGCAGATGCTCGACCGGCTgg                                                                       >  1:872170/1‑108 (MQ=255)                                                                                        aaTGAAATTATAACTTCCGCTATTGATAAGCTGAAGCGTTCATAACGTGAGTTGCCTATGCACAGTGGGGG‑ATTCCCGCCGGCACGGTGCAGATGCTCGACCGGCTgg                                                                       <  2:872170/108‑1 (MQ=255)                                                                                              aTTATAACTTCCGCTATTGATAAGCTGAAGCGTTCATAACGTGAGTTGCCTATGCACAGTGGGGG‑ATTCCCGCCGGCACGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAg                                  <  1:123584/139‑1 (MQ=255)                                                                                                    aCTTCCGCTATTGATAAGCTGAAGCGTTCATAACGTGAGTTGCCTATGCACAGTGGGGG‑ATTCCCGCCGGCACGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGccgccg                            <  2:659178/139‑1 (MQ=255)                                                                                                       tCCGCTATTGATAAGCTGAAGCGTTCATAACGTGAGTTGCCTATGCACAGTGGGGG‑ATTCCCGCCGGCACGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGGTGAGTGAGGaa                                    >  1:457762/1‑128 (MQ=255)                                                                                                       tCCGCTATTGATAAGCTGAAGCGTTCATAACGTGAGTTGCCTATGCACAGTGGGGG‑ATTCCCGCCGGCACGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGGTGAGTGAGGaa                                    <  2:457762/128‑1 (MQ=255)                                                                                                                 aTAAGCTGAAGCGTTCATAACGTGAGTTGCCTATGCACAGTGGGGG‑ATTCCCGCCGGCACGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCGCTGGCTGGCATCg               <  1:373075/139‑1 (MQ=255)                                                                                                                      cTGAAGCGTTCATAACGTGAGTTGCCTATGCACAGTGGGGG‑ATTCCCGCCGGCACGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTGCgg          >  2:549983/1‑139 (MQ=255)                                                                                                                        gAAGCGTTCATAACGTGAGTTGCCTATGCACAGTGGGGGAATTCCCGCCGGCACGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGc         >  1:351765/1‑139 (MQ=255)                                                                                                                        gAAGCGTTCATAACGTGAGTTGCCTATGCACAGTGGGGG‑ATTCCCGCCGGCACGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGcc        >  1:195454/1‑139 (MQ=255)                                                                                                                        gAAGCGTTCATAACGTGAGTTGCCTATGCACAGTGGGGG‑ATTCCCGCCGGCACGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGcc        >  2:40338/1‑139 (MQ=255)                                                                                                                           gCGTTCATAACGTGAGTTGCCTATGCACAGTGGGGG‑ATTCCCGCCGGCACGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGa     >  1:328656/1‑139 (MQ=255)                                                                                                                              ttCATAACGTGAGTTGCCTATGCACAGTGGGGG‑ATTCCCGCCGGCACGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGACCg  <  1:423637/139‑1 (MQ=255)                                                                                                                                    |                                                                                                                                      TATGATGATATCATTGAAGCTGCGGGTGAATTAGGTGATAAAACGCTACTTCCTGTTTTAGATACTATGTTGTACAAGTTTGATGACAATGAAATTATAACTTCCGCTATTGATAAGCTGAAGCGTTCATAGCGTGAGTTGCCTATGCACAGTGGGGG‑ATTCCCGCCGGCACGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGACCG  >  NC_000913/3767048‑3767311

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