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%	153.1 / NA	48	*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 (25/23);  total (25/23)

GATGATATCATTGAAGCTGCGGGTGAATTAGGTGATAAAACGCTACTTCCTGTTTTAGATACTATGTTGTACAAGTTTGATGACAATGAAATTATAACTTCCGCTATTGATAAGCTGAAGCGTTCATAGCGTGAGTTGCCTATGCACAGTGGGGGATTCCCGCCGGCACGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGACCGTG  >  NC_000913/3767051‑3767313                                                                                                                                 |                                                                                                                                       gatgatATCATTGAAGCTGCGGGTGAATTAGGTGATAAAACGCTACTTCCTGTTTTAGATACTATGTTGTACAAGTTTGATGACAATGAAATTATAACTTCCGCTATTGATAAGCTGAAGCGTTCATAACGTGAGTTGc                                                                                                                              <  2:150743/139‑1 (MQ=255)    gatATCATTGAAGCTGCGGGTGAATTAGGTGATAAAACGCTACTTCCTGTTTTAGATACTATGTTGTACAAGTTTGATGACAATGAAATTATAACTTCCGCTATTGATAAGCTGAAGCGTTCATAACGTGAGTTGCCTa                                                                                                                           >  1:713913/1‑139 (MQ=255)      tatCATTGAAGCTGCGGGTGAATTAGGTGATAAAACGCTACTTCCTGTTTTAGATACTATGTTGTACAAGTTTGATGACAATGAAATTATAACTTCCGCTATTGATAAGCTGAAGCGTTCATAACGTGa                                                                                                                                   >  2:441918/1‑129 (MQ=255)      tatCATTGAAGCTGCGGGTGAATTAGGTGATAAAACGCTACTTCCTGTTTTAGATACTATGTTGTACAAGTTTGATGACAATGAAATTATAACTTCCGCTATTGATAAGCTGAAGCGTTCATAACGTGa                                                                                                                                   <  1:441918/129‑1 (MQ=255)      tatCATTGAAGCTGCGGGTGAATTAGGTGATAAAACGCTACTTCCTGTTTTAGATACTATGTTGTACAAGTTTGATGACAATGAAATTATAACTTCCGCTATTGATAAGCTGAAGCGTTCATAACGTGAGt                                                                                                                                 <  2:691694/131‑1 (MQ=255)      tatCATTGAAGCTGCGGGTGAATTAGGTGATAAAACGCTACTTCCTGTTTTAGATACTATGTTGTACAAGTTTGATGACAATGAAATTATAACTTCCGCTATTGATAAGCTGAAGCGTTCATAACGTGAGt                                                                                                                                 >  1:691694/1‑131 (MQ=255)      tatCATTGAAGCTGCGGGTGAATTAGGTGATAAAACGCTACTTCCTGTTTTAGATACTATGTTGTACAAGTTTGATGACAATGAAATTATAACTTCCGCTATTGATAAGCTGAAGCGTTCATAACGTGAGTTGCCTATg                                                                                                                         >  1:303761/1‑139 (MQ=255)        tCATTGAAACTGCGGGTGAATTAGGTGATAAAACGCTACTTCCTGTTTTAGATACTATGTTGTACAAGTTTGATGACAATGAAATTATAACTTCCGCTATTGATAAGCTGAAGCGTTCATAACGTGAGTTGCCTATGca                                                                                                                       <  1:1033413/139‑1 (MQ=255)            tGAAGCTGCGGGTGAATTAGGTGATAAAACGCTACTTCCTGTTTTAGATACTATGTTGTACAAGTTTGATGACAATGAAATTATAACTTCCGCTATTGATAAGCTGAAGCGTTCATAACGTGAGTTGCCTATg                                                                                                                         <  2:533790/133‑1 (MQ=255)            tGAAGCTGCGGGTGAATTAGGTGATAAAACGCTACTTCCTGTTTTAGATACTATGTTGTACAAGTTTGATGACAATGAAATTATAACTTCCGCTATTGATAAGCTGAAGCGTTCATAACGTGAGTTGCCTATg                                                                                                                         >  1:533790/1‑133 (MQ=255)               aGCTGCGGGTGAATTAGGTGATAAAACGCTACTTCCTGTTTTAGATACTATGTTGTACAAGTTTGATGACAATGAAATTATAACTTCCGCTATTGATAAGCTGAAGCGTTCATAACGTGAGt                                                                                                                                 >  1:1028342/1‑122 (MQ=255)               aGCTGCGGGTGAATTAGGTGATAAAACGCTACTTCCTGTTTTAGATACTATGTTGTACAAGTTTGATGACAATGAAATTATAACTTCCGCTATTGATAAGCTGAAGCGTTCATAACGTGAGt                                                                                                                                 <  2:1028342/122‑1 (MQ=255)                    cGGGTGAATTAGGTGATAAAACGCTACTTCCTGTTTTAGATACTATGTTGTACAAGTTTGATGACAATGAAATTATAACTTCCGCTATTGATAAGCTGAAGCGTTCATAACGTGAGTTGCCTATGCACAGTGGGGGAtt                                                                                                           >  2:133482/1‑139 (MQ=255)                              aGGTGATAAAACGCTACTTCCTGTTTTAGATACTATGTTGTACAAGTTTGATGACAATGAAATTATAACTTCCGCTATTGATAAGCTGAAGCGTTCATAACGTGAGTTGCCTATGCACAGTGGGGGATTCCCGCCGGCa                                                                                                 <  2:303761/139‑1 (MQ=255)                                gTGAT‑AAACGCTACTTCCTGTTTTAGATACTATGTTGTACAAGTTTGATGACAATGAAATTATAACTTCCGCTATTGATAAGCTGAAGCGTTCATAACGTGAGTTGCCTATGCACAGTGGGGGATTCCCGCCGGCACgg                                                                                              >  2:804555/1‑139 (MQ=255)                                     aaaaCGCTACTTCCTGTTTTAGATACTATGTTGTACAAGTTTGATGACAATGAAATTATAACTTCCGCTATTGATAAGCTGAAGCGTTCATAACGTGAGTTGCCTATGCACAGTGGGGGATTCCCGCCGGCACGGTGCa                                                                                          <  1:882834/139‑1 (MQ=255)                                         cGCTACTTCCTGTTTTAGATACTATGTTGTACAAGTTTGATGACAATGAAATTATAACTTCCGCTATTGATAAGCTGAAGCGTTCATAACGTGAGTTGCCTATGCACAGTGGGGGATTCCCGCCGGCACGGTGCAGATg                                                                                      >  1:315299/1‑139 (MQ=255)                                                    gTTTTAGATACTATGTTGTACAAGTTTGATGACAATGAAATTATAACTTCCGCTATTGATAAGCTGAAGCGTTCATAACGTGAGTTGCCTATGca                                                                                                                       >  2:612988/1‑95 (MQ=255)                                                    gTTTTAGATACTATGTTGTACAAGTTTGATGACAATGAAATTATAACTTCCGCTATTGATAAGCTGAAGCGTTCATAACGTGAGTTGCCTATGca                                                                                                                       <  1:612988/95‑1 (MQ=255)                                                               tATGTTGTACAAGTTTGATGACAATGAAATTATAACTTCCGCTATTGATAAGCTGAAGCGTTCATAACGTGAGTTGCCTATGCACAGTGGGGGATTCCCGCCGGCACGGTGCAGATGCTCGACCGGCTGGAAAGTGaaa                                                                >  1:915813/1‑139 (MQ=255)                                                                    tgtACAAGTTTGATGACAATGAAATTATAACTTCCGCTATTGATAAGCTGAAGCGTTCATAACGTGAGTTGCCTATGCACAGTGGGGGATTCCCGCCGGCACGGTGCAGATGCTCGACCGGCTGGAAAGTg                                                                   >  1:773476/1‑131 (MQ=255)                                                                    tgtACAAGTTTGATGACAATGAAATTATAACTTCCGCTATTGATAAGCTGAAGCGTTCATAACGTGAGTTGCCTATGCACAGTGGGGGATTCCCGCCGGCACGGTGCAGATGCTCGACCGGCTGGAAAGTg                                                                   <  2:773476/131‑1 (MQ=255)                                                                     gtACAAGTTTGATGACAATGAAATTATAACTTCCGCTATTGATAAGCTGAAGCGTTCATAACGTGAGTTGCCTATGCACAGTGGGGGATTCCCGCCGGCACGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTgg                                                          >  1:257915/1‑139 (MQ=255)                                                                      tACAAGTTTGATGACAATGAAATTATAACTTCCGCTATTGATAAGCTGAAGCGTTCATAACGTGAGTTGCCTATGCACAGTGGGGGATTCCCGCCGGCACGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTGGc                                                         >  1:982014/1‑139 (MQ=255)                                                                        cAAGTTTGATGACAATGAAATTATAACTTCCGCTATTGATAAGCTGAAGCGTTCATAACGTGAGTTGCCTATGCACAGAGGGGGATTCCCGCCGGCACGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTGGCTg                                                       >  1:857035/1‑139 (MQ=255)                                                                         aaGTTTGATGACAATGAAATTATAACTTCCGCTATTGATAAGCTGAAGCGTTCATAACGTGAGTTGCCTATGCACAGTGGGGGATTCCCGCCGGCACGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGa                                                      <  2:915813/139‑1 (MQ=255)                                                                          aGTTTGATGACAATGAAATTATAACTTCCGCTATTGATAAGCTGAAGCGTTCATAACGTGAGTTGCCTATGCACAGTGGGGGATTCCCGCCGGCACGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGAc                                                     <  2:315299/139‑1 (MQ=255)                                                                           gTTTGATGACAATGAAATTATAACTTCCGCTATTGATAAGCTGAAGCGTTCATAACGTGAGTTGCCTATGCACAGTGGGGGATTCCCGCCGGCACGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGAcc                                                    >  2:416255/1‑139 (MQ=255)                                                                             ttGATGACAATGAAATTATAACTTCCGCTATTGATAAGCTGAAGCGTTCATAACGTGAGTTGCCTa                                                                                                                           <  2:607293/66‑1 (MQ=255)                                                                             ttGATGACAATGAAATTATAACTTCCGCTATTGATAAGCTGAAGCGTTCATAACGTGAGTTGCCTa                                                                                                                           >  1:607293/1‑66 (MQ=255)                                                                              tgatgaCAATGAAATTATAACTTCCGCTATTGATAAGCTGAAGCGTTCATAACGTGAGTTGCCTATGCACAGTGGGGGATTCCCGCCGGCACGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCggg                                                 <  1:416255/139‑1 (MQ=255)                                                                                atgaCAATGAAATTATAACTTCCGCTATTGATAAGCTGAAGCGTTCATAACGTGAGTTGCCTATGCACAGTGGGGGATTCCCGCCGGCACGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCgg                                                  <  2:444454/136‑1 (MQ=255)                                                                                atgaCAATGAAATTATAACTTCCGCTATTGATAAGCTGAAGCGTTCATAACGTGAGTTGCCTATGCACAGTGGGGGATTCCCGCCGGCACGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCgg                                                  >  1:444454/1‑136 (MQ=255)                                                                                     aaTGAAATTATAACTTCCGCTATTGATAAGCTGAAGCGTTCATAACGTGAGTTGCCTATGCACAGTGGGGGATTCCCGCCGGCACGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGgtgagtga                                          >  1:403240/1‑139 (MQ=255)                                                                                     aaTGAAATTATAACTTCCGCTATTG‑‑AAGCTGAAGCGTTCATAACGTGAGTTGCCTATGCACAGTGGGGGATTCCCGCCGGc                                                                                                  >  1:464320/1‑81 (MQ=255)                                                                                     aaTGAAATTATAACTTCCGCTATTG‑‑AAGCTGAAGCGTTCATAACGTGAGTTGCCTATGCACAGTGGGGGATTCCCGCCGGc                                                                                                  <  2:464320/81‑1 (MQ=255)                                                                                                 aCTTCCGCTATTGATAAGCTGAAGCGTTCATAACGTGAGTTGCCTATGCACAGTGGGGGATTCCCGCCGGCACGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGccgccg                              <  2:982014/139‑1 (MQ=255)                                                                                                 aCTTCCGCTATTGATAAGCTGAAGCGTTCATAACGTGAGTTGCCTATGCACAGTGGGGGATTCCCGCCGGCACGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGccgccg                              <  1:804555/139‑1 (MQ=255)                                                                                                 aCTTCCGCTATTGATAAGCTGAAGCGTTCATAACGTGAGTTGCCTATGCACAGTGGGGGATTCCCGCCGGCACGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGccgccg                              <  1:133482/139‑1 (MQ=255)                                                                                                     ccGCTATTGATAAGCTGAAGCGTTCATAACGTGAGTTGCCTATGCACAGTGGGGGATTCCCGCCGGCACGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGGTGAGTGAGGa                                       <  2:103495/126‑1 (MQ=255)                                                                                                     ccGCTATTGATAAGCTGAAGCGTTCATAACGTGAGTTGCCTATGCACAGTGGGGGATTCCCGCCGGCACGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGGTGAGTGAGGa                                       >  1:103495/1‑126 (MQ=255)                                                                                                      cGCTATTGATAAGCTGAAGCGTTCATAACGTGAGTTGCCTATGCACAGTGGGGGATTCCCGCCGGCACGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCgctggc                         >  2:628362/1‑139 (MQ=255)                                                                                                                aaGCTGAAGCGTTCATAACGTGAGTTGCCTATGCACAGTGGGGGATTCCCGCCGGCa                                                                                                 >  1:774116/1‑57 (MQ=255)                                                                                                                aaGCTGAAGAGTTCATAACGCGAGTTGCCTATGCACAGTGGGGGATTCCCGCCGGCa                                                                                                 <  2:774116/57‑1 (MQ=255)                                                                                                                      aaGCGTTCATAACGTGAGTTGCCTATGCACAGTGGGGGATTCCCGCCGGCACGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCt         <  2:295829/139‑1 (MQ=255)                                                                                                                           ttCATAACGTGAGTTGCCTATGCACAGTGGGGGATTCCCGCCGGCACGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGgtga                                              <  1:981486/97‑1 (MQ=255)                                                                                                                           ttCATAACGTGAGTTGCCTATGCACAGTGGGGGATTCCCGCCGGCACGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGgtga                                              >  2:981486/1‑97 (MQ=255)                                                                                                                             cATAACGTGAGTTGCCTATGCACAGTGGGGGATTCCCGCCGGCACGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGACCGTg  >  1:809544/1‑139 (MQ=255)                                                                                                                                 |                                                                                                                                       GATGATATCATTGAAGCTGCGGGTGAATTAGGTGATAAAACGCTACTTCCTGTTTTAGATACTATGTTGTACAAGTTTGATGACAATGAAATTATAACTTCCGCTATTGATAAGCTGAAGCGTTCATAGCGTGAGTTGCCTATGCACAGTGGGGGATTCCCGCCGGCACGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGACCGTG  >  NC_000913/3767051‑3767313

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