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%	196.4 / NA	60	*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 (33/27);  total (33/27)

TATGATGATATCATTGAAGCTGCGGGTGAATTAGGTGATAAAACGCTACTTCCTGTTTTAGATACTATGTTGTACAAGTTTGATGACAATGAAATTATAACTTCCGCTATTGATAAGCTGAAGCGTTCATAGCGTGAGTTGCCTATGCACAGTGGGGGATTCCCGCCGGCACGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGACCGTG  >  NC_000913/3767048‑3767313                                                                                                                                    |                                                                                                                                       tATGATGATATCATTGAAGCTGCGGGTGAATTAGGTGATAAAACGCTACTTCCTGTTTTAGATACTATGTTGTACAAGTTTGATGACAATGAAATTATAACTTCCGCTATTGATAAGCTGAAGCGTTCATAACGTGAGt                                                                                                                                 <  2:501977/139‑1 (MQ=255)    gatgatATCATTGAAGCTGCGGGTGAATTAGGTGATAAAACGCTACTTCCTGTTTTAGATACTATGTTGTACAAGTTTGATGACAATGAAATTATAACTTCCGCTATTGATAAGCTGAAGCGTTCATAACGTGAGTTGc                                                                                                                              <  2:753773/139‑1 (MQ=255)    gatgatATCATTGAAGCTGCGGGTGAATTAGGTGATAAAACGCTACTTCCTGTTTTAGATACTATGTTGTACAAGTTTGATGACAATGAAATTATAACTTCCGCTATTGATAAGCTGAAGCGTTCATAACGTGAGTTGc                                                                                                                              <  1:1239299/139‑1 (MQ=255)                       cGGGTGAATTAGGTGATAAAACGCTACTTCCTGTTTTAGATACTATGTTGTACAAGTTTGATGACAATGAAATTATAACTTCCGCTATTGATAAGCTGAAGCGTTCATAACGTGAGt                                                                                                                                 <  1:103019/117‑1 (MQ=255)                       cGGGTGAATTAGGTGATAAAACGCTACTTCCTGTTTTAGATACTATGTTGTACAAGTTTGATGACAATGAAATTATAACTTCCGCTATTGATAAGCTGAAGCGTTCATAACGTGAGt                                                                                                                                 >  2:103019/1‑117 (MQ=255)                                    tGATAAAACGCTACTTCCTGTTTTAGATACTATGTTGTACAAGTTTGATGACAATGAAATTATAACTTCCGCTATTGATAAGCTGAAGCGTTCATAACGTGAGTTGCCTATGCACAGTGGGGGATTCCCGCCGGCACgg                                                                                              >  1:133136/1‑139 (MQ=255)                                       tAAAACGCTACTTCCTGTTTTAGATACTATGTTGTACAAGTTTGATGACAATGAAATTATAACTTCCGCTATTGATAAGCTGAAGCGTTCATAACGTGAGTTGCCTATGCACAGTGGGGGATTCCCGCCGGCACGGTGc                                                                                           >  1:729464/1‑139 (MQ=255)                                                        ttttAGATACTATGTTGTACAAGTTTGATGACAATGAAATTATAACTTCCGCTATTGATAAGCTGAAGCGTTCATAACGTGAGTTGCCTATGCACAGTGGGGGATTCCCGCCGGCACGGTGCAGATGCTCGACCGGCTg                                                                          >  2:1247394/1‑139 (MQ=255)                                                                aCTATGTTGTACAAGTTTGATGACAATGAAATTATAACTTCCGCTATTGATAAGCTGAAGCGTTCATAACGTGAGTTGCCTATGCACAGTGGGGGATTCCCGCCGGCACGGTGCAGATGCTCGACCGGCTGGAAAGTGa                                                                  >  1:540335/1‑139 (MQ=255)                                                                    tgttgtACAAGTTTGATGACAATGAAATTATAACTTCCGCTATTGATAAGCTGAAGCGTTCATAACGTGAGTTGCCTATGCACAGt                                                                                                                   >  1:1246961/1‑86 (MQ=255)                                                                    tgttgtACAAGTTTGATGACAATGAAATTATAACTTCCGCTATTGATAAGCTGAAGCGTTCATAACGTGAGTTGCCTATGCACAGt                                                                                                                   <  2:1246961/86‑1 (MQ=255)                                                                    tgttgtACAAGTTTGATGACAATGAAATTATAACTTCCGCTATTGATAAGCTGAAGCGTTCATAACGTGAGTTGCCTATGCACAGTGGGGGATTCCCGCCGGCACGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATc                                                              <  2:540335/139‑1 (MQ=255)                                                                       tgtACAAGTTTGATGACAATGAAATTATAACTTCCGCTATTGATAAGCTGAAGCGTTCATAACGTGAGTTGCCTATGCACAGTGGGGGATTCCCGCCGGCACGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTg                                                           >  2:36385/1‑139 (MQ=255)                                                                        gtACAAGTTTGATGACAATGAAATTATAACTTCCGCTATTGATAAGCTGAAGCGTTCATAACGTGAGTTGCCTATGCACAGTGGGGGATTCCCGCCGGCACGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTgg                                                          <  2:133136/139‑1 (MQ=255)                                                                         tACAAGTTTGATGACAATGAAATTATAACTTCCGCTATTGATAAGCTGAAGCGTTCATAACGTGAGTTGCCTATGCACAGTGGGGGATTCCCGCCGGCACGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTGGc                                                         >  2:1231204/1‑139 (MQ=255)                                                                             aGTTTGATGACAATGAAATTATAACTTCCGCTATTGATAAGCTGAAGCGTTCATAACGTGAGTTGCCTATGCACAGTGGGGGATTCCCGCCGGCACGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGAc                                                     <  1:36385/139‑1 (MQ=255)                                                                              gTTTGATGACAATGAAATTATAACTTCCGCTATTGATAAGCTGAAGCGTTCATAACGTGAGTTGCCTATGCACAGTGGGGGATTCCCGCCGGCACGGTGCAGATGCTCGACCGGCTGGAAAGt                                                                    <  2:1275637/123‑1 (MQ=255)                                                                              gTTTGATGACAATGAAATTATAACTTCCGCTATTGATAAGCTGAAGCGTTCATAACGTGAGTTGCCTATGCACAGTGGGGGATTCCCGCCGGCACGGTGCAGATGCTCGACCGGCTGGAAAGt                                                                    >  1:1275637/1‑123 (MQ=255)                                                                              gTTTGATGACAATGAAATTATAACTTCCGCTATTGATAAGCTGAAGCGTTCATAACGTGAGTTGCCTATGCACAGTGGGGGATTCCCGCCGGCACGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGAcc                                                    >  2:793695/1‑139 (MQ=255)                                                                                ttGATGACAATGAAATTATAACTTCCGCTATTGATAAGCTGAAGCGTTCATAACGTGAGTTGCCTATGCACAGTGGGGGATTCCCGCCGGCACGGTGCAGATGCTCGACCgg                                                                             <  1:1181086/112‑1 (MQ=255)                                                                                ttGATGACAATGAAATTATAACTTCCGCTATTGATAAGCTGAAGCGTTCATAACGTGAGTTGCCTATGCACAGTGGGGGATTCCCGCCGGCACGGTGCAGATGCTCGACCgg                                                                             >  2:1181086/1‑112 (MQ=255)                                                                                 tgatgaCAATGAAATTATAACTTCCGCTATTGATAAGCTGAAGCGTTCATAACGTGAGTTGCCTATGCACAGTGGGGGATTCCCGCCGGCACGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCggg                                                 <  2:729464/139‑1 (MQ=255)                                                                                  gatgaCAATGAAATTATAACTTCCGCTATTGATAAGCTGAAGCGTTCATAACGTGAGTTGCCTATGCACAGTGGGGGATTCCCGCCGGCACGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGgt                                                <  1:1247394/139‑1 (MQ=255)                                                                                   atgaCAATGAAATTATAACTTCCGCTATTGATAAGCTGAAGCGTTCATAACGTGAGTTGCCTATGCACAGTGGGGGATTCCCGCCGGCACGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTgg                                                          <  1:1503836/128‑1 (MQ=255)                                                                                   atgaCAATGAAATTATAACTTCCGCTATTGATAAGCTGAAGCGTTCATAACGTGAGTTGCCTATGCACAGTGGGGGATTCCCGCCGGCACGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTgg                                                          >  2:1503836/1‑128 (MQ=255)                                                                                   atgaCAATGAAATTATAACTTCCGCTATTGATAAGCTGAAGCGTTCATAACGTGAGTTGCCTATGCACAGTGGGGGATTCCCGCCGGCACGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGgtg                                               >  1:1182050/1‑139 (MQ=255)                                                                                   atgaCAATGAAATTATAACTTCCGCTATTGATAAGCTGAAGCGTTCATAACGTGAGTTGCCTATGCACAGTGGGGGATTCCCGCCGGCACGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGgtg                                               >  2:1358076/1‑139 (MQ=255)                                                                                   atgaCAATGAAATTATAACTTCCGCTATTGATAAGCTGAAGCGTTCATAACGTGAGTTGCCTATGCACAGTGGGGGATTCCCGCCGGCACGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGgtg                                               >  1:1239067/1‑139 (MQ=255)                                                                                          tGAAATTATAACTTCCGCTATTGATAAGCTGAAGCGTTCATAACGTGAGTTGCCTATGcac                                                                                                                      <  2:367762/61‑1 (MQ=255)                                                                                          tGAAATTATAACTTCCGCTATTGATAAGCTGAAGCGTTCATAACGTGAGTTGCCTATGcac                                                                                                                      >  1:367762/1‑61 (MQ=255)                                                                                                    aCTTCCGCTATTGATAAGCTGAAGCGTTCATAACGTGAGTTGCCTATGCACAGTGGGGGATTCCCGCCGGCACGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGccgccg                              <  1:793695/139‑1 (MQ=255)                                                                                                       tCCGCTATTGATAAGCTGAAGCGTTCATAACGTGAGTTGCCTATGCACAGTGGGGGATTCCCGCCGGCa                                                                                                 <  1:1470282/69‑1 (MQ=255)                                                                                                       tCCGCTATTGATAAGCTGAAGCGTTCATAACGTGAGTTGCCTATGCACAGTGGGGGATTCCCGCCGGCa                                                                                                 >  2:1470282/1‑69 (MQ=255)                                                                                                        ccGCTATTGATAAGCTGAAGCGTTCATAACGTGAGTTGCCTATGCACAGTGGGGGATTCCCGCCGGCACGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAg                                    <  2:733644/129‑1 (MQ=255)                                                                                                        ccGCTATTGATAAGCTGAAGCGTTCATAACGTGAGTTGCCTATGCACAGTGGGGGATTCCCGCCGGCACGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAg                                    >  1:733644/1‑129 (MQ=255)                                                                                                                gATAAGCTGAAGCGTTCATAACGTGAGTTGCCTATGCACAGTGGGGGATTCccgccg                                                                                                    <  2:1429396/57‑1 (MQ=255)                                                                                                                gATAAGCTGAAGCGTTCATAACGTGAGTTGCCTATGCACAGTGGGGGATTCccgccg                                                                                                    >  1:1429396/1‑57 (MQ=255)                                                                                                                      cTGAAGCGTTCATAACGTGAGTTGCCTATGCACAGTGGGGGATTCCCGCCGGCACGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGgt                                                >  1:109013/1‑103 (MQ=255)                                                                                                                      cTGAAGCGTTCATAACGTGAGTTGCCTATGCACAGTGGGGGATTCCCGCCGGCACGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGgt                                                <  2:109013/103‑1 (MQ=255)                                                                                                                        gAAGCGTTCATAACGTGAGTTGCCTATGCACAGTGGGGGATTCCCGCCGGCACGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCgg                                                  >  2:946733/1‑99 (MQ=255)                                                                                                                        gAAGCGTTCATAACGTGAGTTGCCTATGCACAGTGGGGGATTCCCGCCGGCACGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCgg                                                  <  1:946733/99‑1 (MQ=255)                                                                                                                        gAAGCGTTCATAACGTGAGTTGCCTATGCACAGTGGGGGATTCCCGCCGGCACGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGgt                                                >  2:218822/1‑101 (MQ=255)                                                                                                                        gAAGCGTTCATAACGTGAGTTGCCTATGCACAGTGGGGGATTCCCGCCGGCACGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGgt                                                <  1:218822/101‑1 (MQ=255)                                                                                                                        gAAGCGTTCATAACGTGAGTTGCCTATGCACAGTGGGGGATTCCCGCCGGCACGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCGCtggctggc                     >  1:716088/1‑128 (MQ=255)                                                                                                                        gAAGCGTTCATAACGTGAGTTGCCTATGCACAGTGGGGGATTCCCGCCGGCACGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCGCtggctggc                     <  2:716088/128‑1 (MQ=255)                                                                                                                        gAAGCGTTCATAACGTGAGTTGCCTATGCACAGTGGGGGATTCCCGCCGGCACGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGcc          >  1:850432/1‑139 (MQ=255)                                                                                                                         aaGCGTTCATAACGTGAGTTGCCTATGCACAGTGGGGGATTCCCGCCGGCACGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCt         <  2:850432/139‑1 (MQ=255)                                                                                                                         aaGCGTTCATAACGTGAGTTGCCTATGCACAGTGGGGGATTCCCGCCGGCACGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCt         >  2:439671/1‑139 (MQ=255)                                                                                                                          aGCGTTCATAACGTGAGTTGCCTATGCACAGTGGGGGATTCCCGCCGGCACGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTg        >  2:778600/1‑139 (MQ=255)                                                                                                                           gCGTTCATAACGTGAGTTGCCTATGCACAGTGGGGGATTCCCGCCGGCa                                                                                                 <  2:1426570/49‑1 (MQ=255)                                                                                                                           gCGTTCATAACGTGAGTTGCCTATGCACAGTGGGGGATTCCCGCCGGCa                                                                                                 >  1:1426570/1‑49 (MQ=255)                                                                                                                           gCGTTCATAACGTGAGTTGCCTATGCACAGTGGGGGATTCCCGCCGGCACGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGa       >  1:723156/1‑139 (MQ=255)                                                                                                                              ttCATAACGTGAGTTGCCTATGCACAGTGGGGGATTCCCGCCGGCACGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTgg                                                          <  2:1162840/85‑1 (MQ=255)                                                                                                                              ttCATAACGTGAGTTGCCTATGCACAGTGGGGGATTCCCGCCGGCACGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTgg                                                          >  1:1162840/1‑85 (MQ=255)                                                                                                                              ttCATAACGTGAGTTGCCTATGCACAGTGGGGGATTCCCGCCGGCACGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGACCg    >  1:162030/1‑139 (MQ=255)                                                                                                                               tCATAACGTGAGTTGCCTATGCACAGTGGGGGATTCCCGCCGGCACGGTGc                                                                                           >  2:857778/1‑51 (MQ=255)                                                                                                                               tCATAACGTGAGTTGCCTATGCACAGTGGGGGATTCCCGCCGGCACGGTGc                                                                                           <  1:857778/51‑1 (MQ=255)                                                                                                                                cATAACGTGAGTTGCCTATGCACAGTGGGGGATTCCCGCCGGCa                                                                                                 <  1:211165/44‑1 (MQ=255)                                                                                                                                cATAACGTGAGTTGCCTATGCACAGTGGGGGATTCCCGCCGGCa                                                                                                 >  2:211165/1‑44 (MQ=255)                                                                                                                                cATAACGTGAGTTGCCTATGCACAGTGGGGGATTCCCGCCGGCACGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGACCGTg  >  1:49715/1‑139 (MQ=255)                                                                                                                                    |                                                                                                                                       TATGATGATATCATTGAAGCTGCGGGTGAATTAGGTGATAAAACGCTACTTCCTGTTTTAGATACTATGTTGTACAAGTTTGATGACAATGAAATTATAACTTCCGCTATTGATAAGCTGAAGCGTTCATAGCGTGAGTTGCCTATGCACAGTGGGGGATTCCCGCCGGCACGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGACCGTG  >  NC_000913/3767048‑3767313

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: ‑