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%	213.7 / NA	66	*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 (40/26);  total (40/26)

TATGATGATATCATTGAAGCTGCGGGTGAATTAGGTGATAAAACGCTACTTCCTGTTTTAGATACTATGTTGTACAAGTTTGATGACAATGAAATTATAACTTCCGCTATTGATAAGCTGAAGCGTTCATAGCGTGAGTTGCCTATGCACAGTGGGGGATTCCCGCCGGCACGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGACCGTG  >  NC_000913/3767048‑3767313                                                                                                                                    |                                                                                                                                       tATGATGATATCATTGAAGCTGCGGGTGAATTAGGTGATAAAACGCTACTTCCTGTTTTAGATACTATGTTGTACAAGTTTGATGACAATGAAATTATAACTTCCGCTATTGATAAGCTGAAGCGTTCATAACGTGAGt                                                                                                                                 <  2:772608/139‑1 (MQ=255) tATGATGATATCATTGAAGCTGCGGGTGAATTAGGTGATAAAACGCTACTTCCTGTTTTAGATACTATGTTGTACAAGTTTGATGACAATGAAATTATAACTTCCGCTATTGATAAGCTGAAGCGTTCATAACGTGAGt                                                                                                                                 <  1:219728/139‑1 (MQ=255)         tatCATTGAAGCTGCGGGTGAATTAGGTGATAAAACGCTACTTCCTGTTTTAGATACTATGTTGTACAAGTTTGATGACAATGAAATTATAACTTCCGCTATTGATAAGCTGAAGCGTTCATAACGTGAGt                                                                                                                                 <  1:443739/131‑1 (MQ=255)         tatCATTGAAGCTGCGGGTGAATTAGGTGATAAAACGCTACTTCCTGTTTTAGATACTATGTTGTACAAGTTTGATGACAATGAAATTATAACTTCCGCTATTGATAAGCTGAAGCGTTCATAACGTGAGt                                                                                                                                 >  2:443739/1‑131 (MQ=255)          atCATTGAAGCTGCGGGTGAATTAGGTGATAAAACGCTACTTCCTGTTTTAGATACTATGTTGTACAAGTTTGATGACAATGAAATTATAACTTCCGCTATTGATAAGCTGAAGCGTTCATAACGTGAGTTGCCTATGc                                                                                                                        >  1:791410/1‑139 (MQ=255)                 aaGCTGCGGGTGAATTAGGTGATAAAACGCTACTTCCTGTTTTAGATACTATGTTGTACAAGTTTGATGACAATGAAATTATAACTTCCGCTATTGATAAGCTGAAGCGTTCATAACGTGAGTTGCCTATGCACAGTgg                                                                                                                 >  2:32362/1‑139 (MQ=255)                  aGCTGCGGGTGAATTAGGTGATAAAACGCTACTTCCTGTTTTAGATACTATGTTGTACAAGTTTGATGACAATGAAATTATAACTTCCGCTATTGATAAGCTGAGGCGTTCATAACGTGAGTTGCCTATGca                                                                                                                       >  2:127483/1‑132 (MQ=255)                  aGCTGCGGGTGAATTAGGTGATAAAACGCTACTTCCTGTTTTAGATACTATGTTGTACAAGTTTGATGACAATGAAATTATAACTTCCGCTATTGATAAGCTGAAGCGTTCATAACGTGAGTTGCCTATGca                                                                                                                       <  1:127483/132‑1 (MQ=255)                  aGCTGCGGGTGAATTAGGTGATAAAACGCTACTTCCTGTTTTAGATACTATGTTGTACAAGTTTGATGACAATGAAATTATAACTTCCGCTATTGATAAGCTGAAGCGTTCATAACGTGAGTTGCCTATGCACAGTggg                                                                                                                >  2:115843/1‑139 (MQ=255)                     tGCGGGTGAATTAGGTGATAAAACGCTACTTCCTGTTTTAGATACTATGTTGTACAAGTTTGATGACAATGAAATTATAACTTCCGCTATTGATAAGCTGAAGCGTTCATAACGTGAGTTGCCTATGCACAGTGGGGGa                                                                                                             >  1:657404/1‑139 (MQ=255)                       cGGGTGAATTAGGTGATAAAACGCTACTTCCTGTTTTAGATACTATGTTGTACAAGTTTGATGACAATGAAATTATAACTTCCGCTATTGATAAGCTGAAGCGTTCATAACGTGAGTTGCCTATGCACAGTGGGGGAtt                                                                                                           >  1:558083/1‑139 (MQ=255)                               ttAGGTGATAAAACGCTACTTCCTGTTTTAGATACTATGTTGTACAAGTTTGATGACAATGAAATTATAACTTCCGCTATTGATAAGCTGAAGCGTTCATAACGTGAGTTGCCTATGCACAGTGGGGGATTCCCGCCgg                                                                                                   <  1:364852/139‑1 (MQ=255)                                 aGGTGATAAAACGCTACTTCCTGTTTTAGATACTATGTTGTACAAGTTTGATGACAATGAAATTATAACTTCCGCTATTGATAAGCTGAAGCGTTCATAACGTGAGTTGCCTATGCACAGTGGGGGATTCCCGCCGGCa                                                                                                 <  2:325328/139‑1 (MQ=255)                                             gCTACTTCCTGTTTTAGATACTATGTTGTACAAGTTTGATGACAATGAAATTATAACTTCCGCTATTGATAAGCTGAAGCGTTCATAACGTGAGTTgt                                                                                                                              >  1:274440/1‑97 (MQ=255)                                             gCTACTTCCTGTTTTAGATACTATGTTGTACAAGTTTGATGACAATGAAATTATAACTTCCGCTATTGATAAGCTGAAGCGTTCATAACGTGAGTTGc                                                                                                                              <  2:274440/98‑1 (MQ=255)                                             gCTACTTCCTGTTTTAGATACTATGTTGTACAAGTTTGATGACAATGAAATTATAACTTCCGCTATTGATAAGCTGAAGCGTTCATAACGTGAGTTGCCTATGCACAGTGGGGGATTCCCGCCGGCACGGTGCAGATGc                                                                                     >  1:224858/1‑139 (MQ=255)                                                aCTTCCTGTTTTAGATACTATGTTGTACAAGTTTGATGACAATGAAATTATAACTTCCGCTATTGATAAGCTGAAGCGTTCATAACGTGAGTTGCCTATGCACAGTGGGGGATTCCCGCCGGCACGGTGCAGATGCTCg                                                                                  <  1:115843/139‑1 (MQ=255)                                                 cTTCCTGTTTTAGATACTATGTTGTACAAGTTTGATGACAATGAAATTATAACTTCCGCTATTGATAAGCTGAAGCGTTCATAACGTGAGTTGCCTATGCACAGTGGGGGATTCCCGCCGGCACGGTGCAGATGCTCGa                                                                                 >  2:23250/1‑139 (MQ=255)                                                       gTTTTAGATACTATGTTGTACAAGTTTGATGACAATGAAATTATAACTTCCGCTATTGATAAGCTGAAGCGTTCATAACGTGAGTTGCCTATGCACAGTGGGGGATTCCCGCCGGCACGGTGCAGATGCTCGACCGGCt                                                                           >  2:472458/1‑139 (MQ=255)                                                        ttttAGATACTATGTTGTACAAGTTTGATGACAATGAAATTATAACTTCCGCTATTGATAAGCTGAAGCGTTCATAACGTGAGTTGCCTATGCACAGTGGGGGATTCCCGCCGGCACGGTGCAGATGCTCGACCGGCTg                                                                          >  2:750001/1‑139 (MQ=255)                                                         tttAGATACTATGTTGTACAAGTTTGATGACAATGAAATTATAACTTCCGCTATTGATAAGCTGAAGCGTTCATAACGTGAGTTGCCTATGCACAGTGGGGGATTCCCGCCGGCACGGTGCAGATGCTCGCCCGGCTgg                                                                         >  1:108681/1‑139 (MQ=255)                                                              aTACTATGTTGTACAAGTTTGATGACAATGAAATTATAACTTCCGCTATTGATAAGCTGAAGCGTTCATAACGTGAGTTGCCTATGCACAGTGGGGGATTCCCGCCGGCACGGTGCAGATGCTCGACCGGCTGGAAAGt                                                                    <  2:791410/139‑1 (MQ=255)                                                                  tATGTTGTACAAGTTTGATGACAATGAAATTATAACTTCCGCTATTGATAAGCTGAAGCGTTCATAACGTGAGTTGCCTATGCACAGTGGGGGATTCCCGCCGGCACGGTGCAGATGCTCGACCGGCTGGAAAGTGaaa                                                                >  1:190108/1‑139 (MQ=255)                                                                  tATGTTGTACAAGTTTGATGACAATGAAATTATAACTTCCGCTATTGATAAGCTGAAGCGTTCATAACGTGAGTTGCCTATGCACAGTGGGGGATTCCCGCCGGCACGGTGCAGATGCTCGACCGGCTGGAAAGTGaaa                                                                >  2:463706/1‑139 (MQ=255)                                                                    tgttgtACAAGTTTGATGACAATGAAATTATAACTTCCGCTATTGATAAGCTGAAGCGTTCATAACGTGAGTTg                                                                                                                               <  2:630564/74‑1 (MQ=255)                                                                    tgttgtACAAGTTTGATGACAATGAAATTATAACTTCCGCTATTGATAAGCTGAAGCGTTCATAACGTGAGTTg                                                                                                                               >  1:630564/1‑74 (MQ=255)                                                                    tgttgtACAAGTTTGATGACAATGAAATTATAACTTCCGCTATTGATAAGCTGAAGCGTTCATAACGTGAGTTGCCTATGCACAGTGGGGGATTCCCGCCGGCACGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATc                                                              >  1:433708/1‑139 (MQ=255)                                                                         tACAAGTTTGATGACAATGAAATTATAACTTCCGCTATTGATAAGCTGAAGCGTTCATAACGTGAGTTGCCTATGCACAGTGGGGGATTCCCGCCGGCACGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTGGc                                                         >  1:516401/1‑139 (MQ=255)                                                                            aaGTTTGATGACAATGAAATTATAACTTCCGCTATTGATAAGCTGAAGCGTTCATAACGTGAGTTGCCTATGCACAGTGGGGGATTCCCGCCGGCACGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGa                                                      >  1:100634/1‑139 (MQ=255)                                                                              gTTTGATGACAATGAAATTATAACTTCCGCTATTGATAAGCTGAAGCGTTCATAACGTGAGTTGCCTATGCACAGTGGGGGATTCCCGCCgg                                                                                                   <  1:95932/92‑1 (MQ=255)                                                                              gTTTGATGACAATGAAATTATAACTTCCGCTATTGATAAGCTGAAGCGTTCATAACGTGAGTTGCCTATGCACAGTGGGGGATTCCCGCCgg                                                                                                   >  2:95932/1‑92 (MQ=255)                                                                              gTTTGATGACAATGAAATTATAACTTCCGCTATTGATAAGCTGAAGCGTTCATAACGTGAGTTGCCTATGCACAGTGGGGGATTCCCGCCGGCa                                                                                                 <  2:763141/94‑1 (MQ=255)                                                                              gTTTGATGACAATGAAATTATAACTTCCGCTATTGATAAGCTGAAGCGTTCATAACGTGAGTTGCCTATGCACAGTGGGGGATTCCCGCCGGCa                                                                                                 >  1:763141/1‑94 (MQ=255)                                                                              gTTTGATGACAATGAAATTATAACTTCCGCTATTGATAAGCTGAAGCGTTCATAACGTGAGTTGCCTATGCACAGTGGGGGATTCCCGCCGGCACGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGAcc                                                    >  2:481849/1‑139 (MQ=255)                                                                                ttGATGACAATGAAATTATAACTTCCGCTATTGATAAGCTGAAGCGTTCATAACGTGAGTTGc                                                                                                                              >  2:725110/1‑63 (MQ=255)                                                                                ttGATGACAATGAAATTATAACTTCCGCTATTGATAAGCTGAAGCGTTCATAACGTGAGTTGc                                                                                                                              <  1:725110/63‑1 (MQ=255)                                                                                ttGATGACAATGAAATTATAACTTCCGCTATTGATAAGCTGAAGCGTTCATAACGTGAGTTGCCTATGCACAGTGGGGGATTCCCGCCGGCACGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTGGCt                                                        <  1:267140/133‑1 (MQ=255)                                                                                ttGATGACAATGAAATTATAACTTCCGCTATTGATAAGCTGAAGCGTTCATAACGTGAGTTGCCTATGCACAGTGGGGGATTCCCGCCGGCACGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTGGCt                                                        >  2:267140/1‑133 (MQ=255)                                                                                  gatgaCAATGAAATTATAACTTCCGCTATTGATAAGCTGAAGCGTTCATAACGTGAGTTGCCTATGCACAGTGGGGGATTCCCGCCGGCACGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGgt                                                >  2:611581/1‑139 (MQ=255)                                                                                  gatgaCAATGAAATTATAACTTCCGCTATTGATAAGCTGAAGCGTTCATAACGTGAGTTGCCTATGCACAGTGGGGGATTCCCGCCGGCACGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGgt                                                <  2:153847/139‑1 (MQ=255)                                                                                   atgaCAATGAAATTATAACTTCCGCTATTGATAAGCTGAAGCGTTCATAACGTGAGTTGCCTATGCACAGTGGGGGATTCCCGCCGGCACGGTGCAGATGCTCGACCGGCTGGAAAGTGaaa                                                                <  1:186652/122‑1 (MQ=255)                                                                                   atgaCAATGAAATTATAACTTCCGCTATTGATAAGCTGAAGCGTTCATAACGTGAGTTGCCTATGCACAGTGGGGGATTCCCGCCGGCACGGTGCAGATGCTCGACCGGCTGGAAAGTGaaa                                                                >  2:186652/1‑122 (MQ=255)                                                                                        aaTGAAATTATAACTTCCGCTATTGATAAGCTGAAGCGTTCATAACGTGAGTTGCCTATGCACAGTGGGGGATTCCCGCCGGCACGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGgtgggtga                                          >  1:183510/1‑139 (MQ=255)                                                                                        aaTGAAATTATAACTTCCGCTATTGATAAGCTGAAGCGTTCATAACGTGAGTTGCCTATGCACAGTGGGGGATTCCCGCCGGCACGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGgtgagtga                                          >  1:351994/1‑139 (MQ=255)                                                                                        aaTGAAATTATAACTTCCGCTATTGATAAGCTGAAGCGTTCATAACGTGAGTTGCCTATGCACAGTGGGGGATTCCCGCCGGCACGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGgtgagtga                                          >  1:801853/1‑139 (MQ=255)                                                                                                  taACTTCCGCTATTGATAAGCTGAAGCGTTCATAACGTGAGTTGCCTATGCACAGTGGGGGATTCCCGCCGGCACGGTGCAGATGCTCGACCGGCTGGa                                                                        <  1:716819/99‑1 (MQ=255)                                                                                                  taACTTCCGCTATTGATAAGCTGAAGCGTTCATAACGTGAGTTGCCTATGCACAGTGGGGGATTCCCGCCGGCACGGTGCAGATGCTCGACCGGCTGGa                                                                        >  2:716819/1‑99 (MQ=255)                                                                                                  taACTTCCGCTATTGATAAGCTGAAGCGTTCATAACGTGAGTTGCCTATGCACAGTGGGGGATTCCCGCCGGCACGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAgccgc                                >  2:125218/1‑139 (MQ=255)                                                                                                    aCTTCCGCTATTGATAAGCTGAAGCGTTCATAACGTGAGTTGCCTATGCACAGTGGGGGATTCCCGCCGGCACGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGccgccg                              <  1:297947/139‑1 (MQ=255)                                                                                                    aCTTCCGCTATTGATAAGCTGAAGCGTTCATAACGTGAGTTGCCTATGCACAGTGGGGGATTCCCGCCGGCACGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGccgccg                              <  1:750001/139‑1 (MQ=255)                                                                                                        ccGCTATTGATAAGCTGAAGCGTTCATAACGTGAGTTGCCTATGCACAGTGGGGGATTCCCGCCGGCACGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCgctgg                          >  1:113916/1‑139 (MQ=255)                                                                                                         cGCTATTGATAAGCTGAAGCGTTCATAACGTGAGTTGCCTATGca                                                                                                                       >  1:333125/1‑45 (MQ=255)                                                                                                         cGCTATTGATAAGCTGAAGCGTTCATAACGTGAGTTGCCTATGca                                                                                                                       <  2:333125/45‑1 (MQ=255)                                                                                                         cGCTATTGATAAGCTGAAGCGTTCATAACGTGAGTTGCCTATGCACAGTGGGGGATTCCCGCCGGCACGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCgctggc                         <  2:679164/139‑1 (MQ=255)                                                                                                                     gCTGAAGCGTTCATAACGTGAGTTGCCTATGCACAGTGGGGGATTCCCGCCGGCACGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTGCg             <  2:635179/139‑1 (MQ=255)                                                                                                                      cTGAAGCGTTCATAACGTGAGTTGCCTATGCACAGTGGGGGATTCCCGCCGGCACGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTGGCTCACCGGGTGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTGCgg            >  1:519412/1‑139 (MQ=255)                                                                                                                        gAAGCGTTCATAACGTGAGTTGCCTATGCACAGTGGGGGATTCCCGCCGGCACGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGa                                                      <  1:302189/95‑1 (MQ=255)                                                                                                                        gAAGCGTTCATAACGTGAGTTGCCTATGCACAGTGGGGGATTCCCGCCGGCACGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGa                                                      >  2:302189/1‑95 (MQ=255)                                                                                                                        gAAGCGTTCATAACGTGAGTTGCCTATGCACAGTGGGGGATTCCCGCCGGCACGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGcc          >  1:773955/1‑139 (MQ=255)                                                                                                                        gAAGCGTTCATAACGTGAGTTGCCTATGCACAGTGGGGGATTCCCGCCGGCACGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGcc          >  1:86566/1‑139 (MQ=255)                                                                                                                        gAAGCGTTCATAACGTGAGTTGCCTATGCACAGTGGGGGATTCCCGCCGGCACGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGcc          <  1:32362/139‑1 (MQ=255)                                                                                                                         aaGCGTTCATAACGTGAGTTGCCTATGCACAGTGGGGGATTCCCGCCGGCACGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCt         >  1:184430/1‑139 (MQ=255)                                                                                                                           gCGTTCATAACGTGAGTTGCCTATGCACAGTGGGGGATTCCCGCCGGCACGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGa       <  1:257377/139‑1 (MQ=255)                                                                                                                              ttCATAACGTGAGTTGCCTATGCACAGTGGGGGATTCCCGCCGGCACGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGACCg    >  2:442219/1‑139 (MQ=255)                                                                                                                                cATAACGTGAGTTGCCTATGCACAGTGGGGGATTCCCGCCGGCACGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGACCGTg  >  1:143594/1‑139 (MQ=255)                                                                                                                                cATAACGTGAGTTGCCTATGCACAGTGGGGGATTCCCGCCGGCACGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGACCGTg  <  2:224858/139‑1 (MQ=255)                                                                                                                                    |                                                                                                                                       TATGATGATATCATTGAAGCTGCGGGTGAATTAGGTGATAAAACGCTACTTCCTGTTTTAGATACTATGTTGTACAAGTTTGATGACAATGAAATTATAACTTCCGCTATTGATAAGCTGAAGCGTTCATAGCGTGAGTTGCCTATGCACAGTGGGGGATTCCCGCCGGCACGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGACCGTG  >  NC_000913/3767048‑3767313

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