Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,061,157	G→A	*237* (TAG→TAA)	yihL →	putative DNA‑binding transcriptional regulator

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,061,157	0	G	A	100.0%	230.0 / NA	70	*237* (TAG→TAA)	yihL	putative DNA‑binding transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (36/34);  total (36/34)

TATCTTGGCGGTTCTGCGGGTATGCCAGCGATCTTAATTGAAAAGCATCGCGCCGACCAGCAAGGCAATATTGTCCAAATAGATATTGAATATTGGCGATTTGAGGCCGTAGACCTCATCATTAATCTGTAGGTGTTTTATGGTGACAATAAATAACGCAAGAAAGATTCTACAACGTGTCGACACTCTTCCTCTTTATTTACATGCTTATGCCTTTCATTTAAATATGCGGCTGGAAAGAGTGTTGCCTGCTGATTTACTTG  >  NC_000913/4061026‑4061288                                                                                                                                    |                                                                                                                                    tATCTTGGCGGTTCTGCGGGTATGCCAGCGATCTTAATTGAAAAGCATCGCGCCGACCAGCAAGGCAATATTGTCCAAATAGATATTGAATATTGGCGATTTGAGGCCGTAGACCTCATCATTAATCTGTAAGTGtttt                                                                                                                              >  2:1160364/1‑139 (MQ=255)     ttGGCGGTTCTGCGGGTATGCCAGCGATCTTAATTGAAAAGCATCGCGCCGACCAGCAAGGCAATATTGTCCAAATAGATATTGAATATTGGCGATTTGAGGCCGTAGACCTCATCATTAATCTGTAAGTGTTTTATgg                                                                                                                          <  1:1949546/139‑1 (MQ=255)          ggTTCTGCGGGTATGCCAGCGATCTTAATTGAAAAGCATCGCGCCGACCAGCAAGGCAATATTGTCCAAATAGATATTGAATATTGGCGATTTGAGGCCGTAGACCTCATCATTAATCTGTAAGTGTTTTATGGTGACa                                                                                                                     >  1:409431/1‑139 (MQ=255)           gTTCTGCGGGTATGCCAGCGATCTTAATTGAAAAGCATCGCGCCGACCAGCAAGGCAATATTGTCCAAATAGATATTGAATATTGGCGATTTGAGGCCGTAGACCTCATCATTAATCTGTAAGTGTTTTATGGTGACaa                                                                                                                    >  1:572652/1‑139 (MQ=255)            ttCTGCGGGTATGCCAGCGATCTTAATTGAAAAGCATCGCGCCGACCAGCAAGGCAATATTGTCCAAATAGATATTGAATATTGGCGATTTGAGGCCGTAGACCTCATCATTAATCTGTAAGTGTTTTATGGTg                                                                                                                        >  1:1268375/1‑134 (MQ=255)            ttCTGCGGGTATGCCAGCGATCTTAATTGAAAAGCATCGCGCCGACCAGCAAGGCAATATTGTCCAAATAGATATTGAATATTGGCGATTTGAGGCCGTAGACCTCATCATTAATCTGTAAGTGTTTTATGGTg                                                                                                                        <  2:1268375/134‑1 (MQ=255)                      aTGCCAGCGATCTTAATTGAAAAGCATCGCGCCGACCAGCAAGGCAATATTGTCCAAATAGATATTGAATATTGGCGATTTGAGGCCGTAGACCTCATCATTAATCTGTAAGTGTTTTATGGTGACAATAAATAACGCa                                                                                                         >  1:1714352/1‑139 (MQ=255)                          cAGCGATCTTAATTGAAAAGCATCGCGCCGACCAGCAAGGCAATATTGTCCAAATAGATATTGAATATTGGCGATTTGAGGCCGTAGACCTCATCATTAATCTGTAAGTGTTTTATGGTGACAATAAATAACGCaagaa                                                                                                     <  2:632021/139‑1 (MQ=255)                             cGATCTTAATTGAAAAGCATCGCGCCGACCAGCAAGGCAATATTGTCCAAATAGATATTGAATATTGGCGATTTGAGGCCGTAGACCTCATCATTAATCTGTAAGTGTTTTATGGTGACAATAAATAACGCaagaaaga                                                                                                  >  2:1543281/1‑139 (MQ=255)                                 cTTAATTGAAAAGCATCGCGCCGACCAGCAAGGCAATATTGTCCAAATAGATATTGAATATTGGCGATTTGAGGCCGTAGACCTCATCATTAATCTGTAAGTGTTTTATGGTGACAATAAATAACGCAAGAAAGATTCt                                                                                              >  2:197793/1‑139 (MQ=255)                                      ttGAAAAGCATCGCGCCGACCAGCAAGGCAATATTGTCCAAATAGATATTGAATATTGGCGATTTGAGGCCGTAGACCTCATCATTAATCTGTAAGTGTTTTATGGTGACAATAAATAACGCAAGAAAGATTCTACAAc                                                                                         >  2:1622057/1‑139 (MQ=255)                                           aaGCATCGCGCCGACCAGCAAGGCAATATTGTCCAAGTAGATATTGAATATTGGCGATTTGAGGCCGTAGACCTCATCATTAATCTGTAAGTGTTTTATgg                                                                                                                          >  1:641678/1‑101 (MQ=255)                                           aaGCATCGCGCCGACCAGCAAGGCAATATTGTCCAAGTAGATATTGAATATTGGCGATTTGAGGCCGTAGACCTCATCATTAATCTGTAAGTGTTTTATgg                                                                                                                          <  2:641678/101‑1 (MQ=255)                                           aaGCATCGCGCCGACCAGCAAGGCAATATTGTCCAAATAGATATTGAATATTGGCGATTTGAGGCCGTAGACCTCATCATTAATCTGTAAGTGTTTTATGGTGACAATAAATAACGCaagaaag                                                                                                   <  2:535199/124‑1 (MQ=255)                                           aaGCATCGCGCCGACCAGCAAGGCAATATTGTCCAAATAGATATTGAATATTGGCGATTTGAGGCCGTAGACCTCATCATTAATCTGTAAGTGTTTTATGGTGACAATAAATAACGCaagaaag                                                                                                   >  1:535199/1‑124 (MQ=255)                                           aaGCATCGCGCCGACCAGCAAGGCAATATTGTCCAAATAGATATTGAATATTGGCGATTTGAGGCCGTAGACCTCATCATTAATCTGTAAGTGTTTTATGGTGACAATAAATAACGCaagaa                                                                                                     >  2:240585/1‑122 (MQ=255)                                           aaGCATCGCGCCGACCAGCAAGGCAATATTGTCCAAATAGATATTGAATATTGGCGATTTGAGGCCGTAGACCTCATCATTAATCTGTAAGTGTTTTATGGTGACAATAAATAACGCaagaa                                                                                                     <  1:240585/122‑1 (MQ=255)                                             gCATCGCGCCGACCAGCAAGGCAATATTGTCCAAATAGATATTGAATATTGGCGATTTGAGGCCGTAGACCTCATCATTAATCTGTAAGTGTTTTATGGTGACAATAAATAACGCAAGAAAGATTCTACAACGTGTCGa                                                                                  <  1:1160364/139‑1 (MQ=255)                                             gCATCGCGCCGACCAGCAAGGCAATATTGTCCAAATAGATATTGAATATTGGCGATTTGAGGCCGTAGACCTCATCATTAATCTGTAAGTGTTTTATGGTGACAATAAATAACGCAAGAAAGATTCTACAACGTGTCGa                                                                                  <  2:572652/139‑1 (MQ=255)                                              cATCGCGCCGACCAGCAAGGCAATATTGTCCAAATAGATATTGAATATTGGCGATTTGAGGCCGTAGACCTCATCATTAATCTGTAAGTGTTTTATgg                                                                                                                          >  1:1507373/1‑98 (MQ=255)                                              cATCGCGCCGACCAGCAAGGCAATATTGTCCAAATAGATATTGAATATTGGCGATTTGAGGCCGTAGACCTCATCATTAATCTGTAAGTGTTTTATgg                                                                                                                          <  2:1507373/98‑1 (MQ=255)                                              cATCGCGCCGACCAGCAAGGCAATATTGTCCAAATAGATATTGAATATTGGCGATTTGAGGCCGTAGACCTCATCATTAATCTGTAAGTGTTTTATGGTGACAATAAATAACGCAAGAAAGATTCTACAACGTGTCGac                                                                                 <  2:1714352/139‑1 (MQ=255)                                                 cgcgCCGACCAGCAAGGCAATATTGTCCAAATAGATATTGAATATTGGCGATTTGAGGCCGTAGACCTCATCATTAATCTGTAAGTGTTTTATGGTGACAATAAATAACGCAAGAAAGATTCTACAACGTGTCGACAct                                                                              >  1:998792/1‑139 (MQ=255)                                                     ccGACCAGCAAGGCAATATTGTCCAAATAGATATTGAATATTGGCGATTTGAGGCCGTAGACCTCATCATTAATCTGTAAGTGTTTTATGGTGACAATAAATAACGCAAGAAAGATTCTACAACGTGTCGACACTCTTc                                                                          <  1:197793/139‑1 (MQ=255)                                                          cAGCAAGGCAATATTGTCCAAATAGATATTGAATATTGGCGATTTGAGGCCGTAGACCTCATCATTAATCTGTAAGTGTTTTATGGTGACAATAAATAACGCAAGAAAGATTCTACAACGTGTCGACACTCTTCCTCtt                                                                     <  1:1543281/139‑1 (MQ=255)                                                             cAAGGCAATATTGTCCAAATAGATATTGAATATTGGCGATTTGAGGCCGTAGACCTCATCATTAATCTGTAAGTGTTTTATGGTGACAATAAATAACGCAAGAAAGATTCTACAACGTGTCGACACTCTTCCTCtttat                                                                  <  2:998792/139‑1 (MQ=255)                                                                 gCAATATTGTCCAAATAGATATTGAATATTGGCGATTTGAGGCCGTAGACCTCATCATTAATCTGTAAGTGTTTTATGGTGACAATAAATAACGCAAGAAAGATTCTACAACGTGTCGACACTCTTCCTCTTTATTTAc                                                              <  1:1622057/139‑1 (MQ=255)                                                                      atTGTCCAAATAGATATTGAATATTGGCGATTTGAGGCCGTAGACCTCATCATTAATCTGTAAGTGTTTTATgg                                                                                                                          >  2:432822/1‑74 (MQ=255)                                                                      atTGTCCAAATAGATATTGAATATTGGCGATTTGAGGCCGTAGACCTCATCATTAATCTGTAAGTGTTTTATgg                                                                                                                          <  1:432822/74‑1 (MQ=255)                                                                        tGTCCAAATAGATATTGAATATTGGCGATTTGAGGCCGTAGACCTCATCATTAATCTGTAAGTGTTTTATGGt                                                                                                                         <  1:1720003/73‑1 (MQ=255)                                                                        tGTCCAAATAGATATTGAATATTGGCGATTTGAGGCCGTAGACCTCATCATTAATCTGTAAGTGTTTTATGGt                                                                                                                         >  2:1720003/1‑73 (MQ=255)                                                                          tCCAAATAGATATTGAATATTGGCGATTTGAGGCCGTAGA‑CTCATCATTAATCTGTAAGTGTTTTATGGTGACAATAAATAACGCAAGAAAGATTCTACAACGTGTCGACACTCTTCCTCtttatttg                                                               >  2:1342638/1‑127 (MQ=255)                                                                          tCCAAATAGATATTGAATATTGGCGATTTGAGGCCGTAGA‑CTCATCATTAATCTGTAAGTGTTTTATGGTGACAATAAATAACGCAAGAAAGATTCTACAACGTGTCGACACTCTTCCTCtttattta                                                               <  1:1342638/128‑1 (MQ=255)                                                                           ccAAATAGATATTGAATATTGGCGATTTGAGGCCGTAGACCTCATCATTAATCTGTAAGTGTTTTATgg                                                                                                                          <  1:1730666/69‑1 (MQ=255)                                                                           ccAAATAGATATTGAATATTGGCGATTTGAGGCCGTAGACCTCATCATTAATCTGTAAGTGTTTTATgg                                                                                                                          >  2:1730666/1‑69 (MQ=255)                                                                           ccAAATAGATATTGAATATTGGCAATTTGAGGCCGTAGACCTCATCATTAATCTGTAAGTGTTTTATgg                                                                                                                          >  2:867333/1‑69 (MQ=255)                                                                           ccAAATAGATATTGAATATTGGCAATTTGAGGCCGTAGACCTCATCATTAATCTGTAAGTGTTTTATgg                                                                                                                          <  1:867333/69‑1 (MQ=255)                                                                            cAAATAGATATTGAATATTGGCGATTTGAGGCCGTAGACCTCATCATTAATCTGTAAGTGTTTTATGGTGACAATAAATAACGCAAGAAAGATTCTACAACGTGTCGACACTCTTCCTCTTTATTTACATGCTTATGcc                                                   >  2:1783129/1‑139 (MQ=255)                                                                               aTAGATATTGAATATTGGCGATTTGAGGCCGTAGACCTCATCATTAATCTGTAAGTGTTTTATGGTGACAATAAATAACGCaagaa                                                                                                     >  1:24838/1‑86 (MQ=255)                                                                               aTAGATATTGAATATTGGCGATTTGAGGCCGTAGACCTCATCATTAATCTGTAAGTGTTTTATGGTGACAATAAATAACGCaagaa                                                                                                     <  2:24838/86‑1 (MQ=255)                                                                                    tatTGAATATTGGCGATTTGAGGCCGTAGACCTCATCATTAATCTGTAAGTGTTTTATGGTGACAATAAATAACGCAAGAAAGATTCTACAACGTGTCGACACTCTTCCTCTTTATTTACATGCTTATGCCTTTCAttt                                           >  1:155923/1‑139 (MQ=255)                                                                                         aaTATTGGCGATTTGAGGCCGTAGACCTCATCATTAATCTGTAAGTGTTTTATGGTGACAATAAATAACGCAAGAAAGATTCTACAACGTGTCGACACTCTTCCTCTTTATTTACATGCTTATGCCTTTCATTTAAata                                      >  2:1945302/1‑139 (MQ=255)                                                                                           tatTGGCGATTTGAGGCCGTAGACCTCATCATTAATCTGTAAGTGTTTTATGGTGACAATAAATAACGCAAGAAAGATTCTACAACGTGTCGACACTCTTCCTCTTTATTTACATGCTTATGCCTTTCATTTAAATATg                                    <  2:409431/139‑1 (MQ=255)                                                                                                gCGATTTGAGGCCGTAGACCTCATCATTAATCTGTAAGTGTTTTATGGTGACAATAAATAACGCAAGAAAGATTCTACAACGTGTCGa                                                                                  >  2:912622/1‑88 (MQ=255)                                                                                                gCGATTTGAGGCCGTAGACCTCATCATTAATCTGTAAGTGTTTTATGGTGACAATAAATAACGCAAGAAAGATTCTACAACGTGTCGa                                                                                  >  1:1124728/1‑88 (MQ=255)                                                                                                gCGATTTGAGGCCGTAGACCTCATCATTAATCTGTAAGTGTTTTATGGTGACAATAAATAACGCAAGAAAGATTCTACAACGTGTCGa                                                                                  <  1:912622/88‑1 (MQ=255)                                                                                                gCGATTTGAGGCCGTAGACCTCATCATTAATCTGTAAGTGTTTTATGGTGACAATAAATAACGCAAGAAAGATTCTACAACGTGTCGa                                                                                  <  2:1124728/88‑1 (MQ=255)                                                                                                   aTTTGAGGCCGTAGACCTCATCATTAATCTGTAAGTGTTTTATgg                                                                                                                          <  1:1175623/45‑1 (MQ=255)                                                                                                   aTTTGAGGCCGTAGACCTCATCATTAATCTGTAAGTGTTTTATgg                                                                                                                          >  2:1175623/1‑45 (MQ=255)                                                                                                            cGTAGACCTCATCATTAATCTGTAAGTGTTTTATGGTGACAATAAATAACGCAAGAAAGATTCTACAACGTGTCGACACTCTTCCTCTTTATTTACATGCTTATGCCTTTCATTTAAATATGCGGCTGGAAAGAGTGtt                   <  2:155923/139‑1 (MQ=255)                                                                                                              tAGACCTCATCATTAATCTGTAAGTGTTTTATGGTGACAATAAATAACGCAAGAAAGATTCTACaa                                                                                          <  1:1745570/66‑1 (MQ=255)                                                                                                              tAGACCTCATCATTAATCTGTAAGTGTTTTATGGTGACAATAAATAACGCAAGAAAGATTCTACaa                                                                                          >  2:1745570/1‑66 (MQ=255)                                                                                                              tAGACATCATCATTAATCTGTAAGTGTTTTATGGTGACAATAAATAACGCaagaa                                                                                                     >  1:1650024/1‑55 (MQ=255)                                                                                                              tAGACATCATCATTAATCTGTAAGTGTTTTATGGTGACAATAAATAACGCaagaa                                                                                                     <  2:1650024/55‑1 (MQ=255)                                                                                                                  ccTCATCATTAATCTGTAAGTGTTTTATGGTGACAATAAATAACGCAAGAAAGATTCTACaa                                                                                          <  2:1944465/62‑1 (MQ=255)                                                                                                                  ccTCATCATTAATCTGTAAGTGTTTTATGGTGACAATAAATAACGCAAGAAAGATTCTACaa                                                                                          >  1:1944465/1‑62 (MQ=255)                                                                                                                     catcatTAATCTGTAAGTGTTTTATGGTGACAATAAATAACGCAAGAAAGATTCTACAACGTGTCGa                                                                                  >  2:1271472/1‑67 (MQ=255)                                                                                                                     catcatTAATCTGTAAGTGTTTTATGGTGACAATAAATAACGCAAGAAAGATTCTACAACGTGTCGa                                                                                  <  1:1271472/67‑1 (MQ=255)                                                                                                                     catcatTAATCTGTAAGTGTTTTATGGTGACAATAAATAACGCAAGAAAGATTCTACAACGTGTCGACACTCTTCCTCTTTATTTACATGCTTATGCCTTTCATTTAAATATGCGGCTGGAAAGAGTGTTGCCTGCTGa          >  2:1913973/1‑139 (MQ=255)                                                                                                                     catcatTAATCTGTAAGTGTTTTATGGTGACAATAAATAACGCAAGAAAGATTCTACAACGTGTCGACACTCTTCCTCTTTATTTACATGCTTATGCCTTTCATTTAAATATGCGGCTGGAAAGAGTGTTGCCTGCTGa          >  1:414923/1‑139 (MQ=255)                                                                                                                      atcatTAATCTGTAAGTGTTTTATGGTGACAATAAATAACGCAAGAAAGATTCTAc                                                                                            >  1:1703410/1‑56 (MQ=255)                                                                                                                      atcatTAATCTGTAAGTGTTTTATGGTGACAATAAATAACGCAAGAAAGATTCTAc                                                                                            <  2:1703410/56‑1 (MQ=255)                                                                                                                      atcatTAATCTGTAAGTGTTTTATGGTGACAATAAATAACGCAAGAAAGATTCTACAACGTGTCGaca                                                                                <  2:1195332/68‑1 (MQ=255)                                                                                                                      atcatTAATCTGTAAGTGTTTTATGGTGACAATAAATAACGCAAGAAAGATTCTACAACGTGTCGaca                                                                                >  1:1195332/1‑68 (MQ=255)                                                                                                                      atcatTAATCTGTAAGTGTTTTATGGTGACAATAAATAACGCAAGAAAGATTCTACAACGTGTCGa                                                                                  >  1:454008/1‑66 (MQ=255)                                                                                                                      atcatTAATCTGTAAGTGTTTTATGGTGACAATAAATAACGCAAGAAAGATTCTACAACGTGTCGa                                                                                  <  2:454008/66‑1 (MQ=255)                                                                                                                      atcatTAATCTGTAAGTGTTTTATGGTGACAATAAATAACGCAAGAAAGATTCTACAACGTGTCGACACTCTTCCTCTTTATTTACATGCTTATGCCTTTCATTTAAATATGCGGCTGGAAAGAGTGTTGCCTGCTGAt         <  2:414923/139‑1 (MQ=255)                                                                                                                             aTCTGTAAGTGTTTTATGATGACAATAAATAACGCAAGAAAGATTCTACAACGTGTCGACACTCTTCCTCTTTATTTACATGCTTATGCCTTTCATTTAAATATGCGGCTGGAAAGAGTGTTGCCTGCTGATTTACTTg  >  1:965975/1‑139 (MQ=255)                                                                                                                                tGTAAGTGTTTTATGGTGACAATAAATAACGCAAGAAAGATTCTAc                                                                                            <  1:1767384/46‑1 (MQ=255)                                                                                                                                tGTAAGTGTTTTATGGTGACAATAAATAACGCAAGAAAGATTCTAc                                                                                            >  2:1767384/1‑46 (MQ=255)                                                                                                                                    |                                                                                                                                    TATCTTGGCGGTTCTGCGGGTATGCCAGCGATCTTAATTGAAAAGCATCGCGCCGACCAGCAAGGCAATATTGTCCAAATAGATATTGAATATTGGCGATTTGAGGCCGTAGACCTCATCATTAATCTGTAGGTGTTTTATGGTGACAATAAATAACGCAAGAAAGATTCTACAACGTGTCGACACTCTTCCTCTTTATTTACATGCTTATGCCTTTCATTTAAATATGCGGCTGGAAAGAGTGTTGCCTGCTGATTTACTTG  >  NC_000913/4061026‑4061288

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