Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,075,454	G→A	*262* (TAG→TAA)	yihW →	putative transcriptional regulator for sulphoquinovose utilization

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,075,454	0	G	A	100.0%	145.9 / NA	46	*262* (TAG→TAA)	yihW	putative transcriptional regulator for sulphoquinovose utilization
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (18/28);  total (18/28)

ATGCTTCGGCGGCGGTTGAAATTGGTAACGTGGCACAGGTCACTGCGCTCTTTACCGACGAGCTGCCGCCCGCTGCGCTAAAATCACGCTTACAAGACAGCCAAATTGAAATCATTCTTCCCCAGGAAGACGCGTAGATTTTTGTGACCGTTAACCTGGCTTCATACCTTGCCACATAGCCAAACCCATCCTTTCCCGCTACAGTTAATTTCTTGTGGCGCGAAAGGAGGCAAAAATGCTCTATATCTTTGATTTAGGTAATGT  >  NC_000913/4075318‑4075581                                                                                                                                         |                                                                                                                                aTGCTTCGGCGGCGGTTGAAATTGGTAACGTGGCACAGGTCACTGCGCTCTTTACCGACGAGCTGCCGCCCGCTGCGCTAAAATCACGCTTACAAGACAGCCAAATTGAAATCATTCTTCCCCAGGAAGACGCGTAAAt                                                                                                                               >  2:191854/1‑139 (MQ=255) aTGCTTCGGCGGCGGTTGAAATTGGTAACGTGGCACAGGTCACTGCGCTCTTTACCGACGAGCTGCCGCCCGCTGCGCTAAAATCACGCTTACAAGACAGCCAAATTGAAATCATTCTTCCCCAGGAAGACGCGTAAAt                                                                                                                               >  1:628498/1‑139 (MQ=255) aTGCTTCGGCGGCGGTTGAAATTGGTAACGTGGCACAGGTCACTGCGCTCTTTACCGACGAGCTGCCGCCCGCTGCGCTAAAATCACGCTTACAAGACAGCCAAATTGAAATCATTCTTCCCCAGGAAGACGCGTAAAt                                                                                                                               >  1:395074/1‑139 (MQ=255)       cggcggcggTTGACATTGTTAACGTGCCACAGGAGACTGCGCTCTTTACCGACGAGCTGCCGCCCGCTGCGCTAAAATAACGCTTAGAAGACAGCCAAATTGAAATCATTCTTCCCCAGTAAGAAGCGTAAATTTTagt                                                                                                                         >  1:324772/1‑139 (MQ=255)       cggcggcggTTGAAATTGGTAACGTGGCACAGGTCACTGCGCTCTTTACCGACGAGCTGCCGCCCGCTGCGCTAAAATCACGCTTACAAGACAGCCAAATTGAAATCATTCTTCCCCAGGAAGACGCGTAAATTTTtgt                                                                                                                         <  1:153086/139‑1 (MQ=255)       cggcggcggTTGAAATTGGTAACGTGGCACAGGTCACTGCGCTCTTTACCGACGAGCTGCCGCCCGCTGCGCTAAAATCACGCTTACAAGACAGCCAAATTGAAATCATTCTTCCCCAGGAAGACGCGTAAATTTTtgt                                                                                                                         <  2:704584/139‑1 (MQ=255)       cggcggcggTTGAAATTGGTAACGTGGCACAGGTCACTGCGCTCTTTACCGACGAGCTGCCGCCCGCTGCGCTAAAATCACGCTTACAAGACAGCCAAATTGAAATCATTCTTCCCCAGGAAGACGCGTAAATTTTtgt                                                                                                                         >  1:369290/1‑139 (MQ=255)             cggTTGAAATTGGTAACGTGGCACAGGTCACTGCGCTCTTTACCGACGAGCTGCCGCCCGCTGCGCTAAAATCACGCTTACAAGACAGCCAAATTGAAATCATTCTTCCCCAGGAAGACGCGTAAATTTTTGTGACCGt                                                                                                                   >  2:496894/1‑139 (MQ=255)                  gAAATTGGTCACGTGGCACAGGTCACTGCGCTCTTTACCGACGAGCTGCCGCCCGCTGCGCTAAAATCACGCTTACAAGACAGCCAAATTGAAATCATTCTTCCCCAGGAAGACGCGTAAATTTTTGTGACCGTTAAcc                                                                                                              <  2:515255/139‑1 (MQ=255)                  gAAATTGGTAACGTGGCACAGGTCACTGCGCTCTTTACCGACGAGCTGCCGCCCGCTGCGCTAAAATCACGCTTACAAGACAGCCAAATTGAAATCATTCTTCCCCAGGAAGACGCGTAAATTTTTGTGACCGTTAAcc                                                                                                              <  1:96348/139‑1 (MQ=255)                   aaaTTGGTAACGTGGCACAGGTCACTGCGCTCTTTACCGACGAGCTGCCGCCCGCTGCGCTAAAATCACGCTTACAAGACAGCCAAATTGAAATCATTCTTCCCCAGGAAGACGCGTAAATTTTTGTGACCGTTAACCt                                                                                                             <  2:628498/139‑1 (MQ=255)                   aaaTTGGTAACGTGGCACAGGTCACTGCGCTCTTTACCGACGAGCTGCCGCCCGCTGCGCTAAAATCACGCTTACAAGACAGCCAAATTGAAATCATTCTTCCCCAGGAAGACGCGTAAATTTTTGTGACCGTTAACCt                                                                                                             <  2:33031/139‑1 (MQ=255)                                      ggTCACTGCGCTCTTTACCGACGAGCTGCCGCCCGCTGCGCTAAAATCACGCTTACAAGACAGCCAAATTGAAATCATTCTTCCCCAGGAAGACGCGTAAATTTTTGTGACCGTTAACCTGGCTTCATACCTTGCcaca                                                                                          <  1:496894/139‑1 (MQ=255)                                       gTCACTGCGCTCTTTACCGACGAGCTGCCGCCCGCTGCGCTAAAATCACGCTTACAAGACAGCCAAATTGAAATCATTCTTCCCCAGGAAGACGCGTAAATTTTTGTGACCGTTAACCTGGCTTCATACCTTGCCACAt                                                                                         >  2:305969/1‑139 (MQ=255)                                          aCTGCGCTCTTTACCGACGAGCTGCCGCCCGCTGCGCTAAAATCACGCTTACAAGACAGCCAAATTGAAATCATTCTTCCCCAGGAAGACGCGTAAATTTTTGTGACCGTTAACCTGGCTTCATACCTTGCCACATAGc                                                                                      <  2:395074/139‑1 (MQ=255)                                            tGCGCTCTTTACCGACGAGCTGCCGCCCGCTGCGCTAAAATCACGCTTACAAGACAGCCAAATTGAAATCATTCTTCCCCAGGAAGACGCGTAAATTTTTGTGACCGTTAACCTGGCTTCATACCTTGCCACATAGCCa                                                                                    >  2:485113/1‑139 (MQ=255)                                                   tttACCGACGAGCTGCCGCCCGCTGCGCTAAAATCACGCTTACAAGACAGCCAAATTGAAATCATTCTTCCCCAGGAAGACGCGTAAATTTTTGTGACCGTTAACCTGGCTTCATACCTTg                                                                                               <  2:571961/121‑1 (MQ=255)                                                   tttACCGACGAGCTGCCGCCCGCTGCGCTAAAATCACGCTTACAAGACAGCCAAATTGAAATCATTCTTCCCCAGGAAGACGCGTAAATTTTTGTGACCGTTAACCTGGCTTCATACCTTg                                                                                               >  1:571961/1‑121 (MQ=255)                                                          acgaGCTGCCGCCCGCTGCGCTAAAATCACGCTTACAAGACAGCCAAATTGAAATCATTCTTCCCCAGGAAGACGCGTAAATTTTTGTGACCGTTAACCTGGCTTCATACCTTGCCACATAGCCAAACCCATCCTTTcc                                                                      >  2:773315/1‑139 (MQ=255)                                                           cgaGCTGCCGCCCGCTGCGCTAAAATCACGCTTACAAGACAGCCAAATTGAAATCATTCTTCCCCAGGAAGACGCGTAAATTTTTGTGACCGTTAACCTGGCTTCATACCTTGCCACATAGCCAAACCCATCCTTTccc                                                                     >  1:392486/1‑139 (MQ=255)                                                               cTGCCGCCCGCTGCGCTAAAATCACGCTTACAAGACAGCCAAATTGAAATCATTCTTCCCCAGGAAGACGCGTAAATTTTTGTGACCGTTAACCTGGCTTCATACCTTGCCACATAGCCAAACCCATCCTTTCCCGCTa                                                                 <  1:275596/139‑1 (MQ=255)                                                                tGCCGCCCGCTGCGCTAAAATCACGCTTACAAGACAGCCAAATTGAAATCATTCTTCCCCAGGAAGACGCGTAAATTTTTGTGACCGTTAACCTGGCTTCATACCTTGCCACATAGCCAAACCCATCCTTTCCCGCTAc                                                                <  1:773315/139‑1 (MQ=255)                                                                          tGCGCTAAAATCACGCTTACAAGACAGCCAAATTGAAATCATTCTTCCCCAGGAAGACGCGTAAATTTTTGTGACCGTTAACCTGGCTTCATACCTTGCCACATAGCCAAACCCATCCTTTCCCg                                                                    <  2:284607/125‑1 (MQ=255)                                                                          tGCGCTAAAATCACGCTTACAAGACAGCCAAATTGAAATCATTCTTCCCCAGGAAGACGCGTAAATTTTTGTGACCGTTAACCTGGCTTCATACCTTGCCACATAGCCAAACCCATCCTTTCCCg                                                                    >  1:284607/1‑125 (MQ=255)                                                                              cTAAAATCACGCTTACAAGACAGCCAAATTGAAATCATTCTTCCCCAGGAAGACGCGTAAATTTTTGTGACCGTTAACCTGGCTTCATACCTTGCCACATAGCCAAACCCATCCTTTCCCGCTACAGTTAATTTCTtgt                                                  <  2:392486/139‑1 (MQ=255)                                                                                      aCGCTTACAAGACAGCCAAATTGAAATCATTCTTCCCCAGGAAGACGCGTAAATTTTTGTGACCGTTAACCTGGCTTCATACCTTGCCACATAGCCAAACCCATCCTTTCCCGCTACAGTTAATTTCTTGTGGCGCGaa                                          <  1:106992/139‑1 (MQ=255)                                                                                       cGCTTACAAGACAGCCAAATTGAAATCATTCTTCCCCAGGAAGACGCGTAAATTTTTGTGACCGTTAACCTGGCTTCATACCTTGCCACATAGCCAAACCCATCCTTTCCCGCTACAGTTAATTTCTTGTGGCGCGaaa                                         <  2:727174/139‑1 (MQ=255)                                                                                            aCAAGACAGCCAAATTGAAATCATTCTTCCCCAGGAAGACGCGTAAATTTTTGTGACCGTTAACCTGGCTTCATACCTTGCCACATAGCCAAACCCATCCTTTCCCGCTACAGTTAATTTCTTGTGGCGCGAAaggagg                                    <  1:485113/139‑1 (MQ=255)                                                                                                  cAGCCAAATTGAAATCATTCTTCCCCAGGAAGACGCGTAAATTTTTGTGACCGTTAACCTGGCTTCATACCTTGCCACATAGCCAAACCCATCCTTTCCCGCTACAGTTAATTTCTTGTGGCGCGAAAGGAGGCaaaaa                              >  2:146616/1‑139 (MQ=255)                                                                                                  cAGCCAAATTGAAATCATTCTTCCCCAGGAAGACGCGTAAATTTTTGTGACCGTTAACCTGGCTTCATACCTTGCCACATAGCCAAACCCATCCTTTCCCGCTACAGTTAATTTCTTGTGGCGCGAAAGGAGGCaaaaa                              >  1:288469/1‑139 (MQ=255)                                                                                                  cAGCCAAATTGAAATCATTCTTCCCCAGGAAGACGCGTAAATTTTTGTGACCGTTAACCTGGCTTCATACCTTGCCACATAGCCAAACCCATCCTTTCCCGCTACAGTTAATTTCTTGTGGCGCGAAAGGAGGCaaaaa                              <  2:288469/139‑1 (MQ=255)                                                                                                    gCCAAATTGAAATCATTCTTCCCCAGGAAGACGCGTAAATTTTTGTGACCGTTAACCTGGCTTCATACCTTGCcaca                                                                                          <  1:363335/77‑1 (MQ=255)                                                                                                    gCCAAATTGAAATCATTCTTCCCCAGGAAGACGCGTAAATTTTTGTGACCGTTAACCTGGCTTCATACCTTGCcaca                                                                                          >  2:363335/1‑77 (MQ=255)                                                                                                    accAAATTGAAATCATTCTTCCCCAGGAAGACGCGTAAATTTTTGTGACCGTTAACCTGGCTTCATACCTTGCCACATAGCCAAACCCATCCTTTCCCGCTACAGTTAATTTCTTGTGGCGCGAAAGGAGGCAAAAATg                            >  2:418990/2‑139 (MQ=255)                                                                                                      cAAATTGAAATCATTCTTCCCCAGGAAGACGCGTAAATTTTTGTGACCGTTAACCTGGCTTCATACCTTGCCACATAGCCAAACCCATCCTTTCCCGCTACAGTTAATTTCTTGTGGCGCGAAAGGAGGCAAAAATGct                          <  1:418990/139‑1 (MQ=255)                                                                                                          ttGAAATCATTCTTCCCCAGGAAGACGCGTAAATTTTTGTGACCGTTAACCTGGCTTCATACCTTGCCACATAGCCAAACCCATCCTTTCCCGCTACAGTTAATTTCTTGTGGCGCGAAAGGAGGCAAAAATGCTCtat                      <  1:305969/139‑1 (MQ=255)                                                                                                              aaTCATTCTTCCCCAGGAAGACGCGTAAATTTTTGTGACCGTTAACCTGGCTTCATACCTTGCCACATAGCCAAACCCATCCTTTCCCGCTACAGTTAATTTCTTGTGGCGCGAAAGGAGGCAAAAATGCTCTATATCt                  <  1:493937/139‑1 (MQ=255)                                                                                                              aaTCATTCTTCCCCAGGAAGACGCGTAAATTTTTGTGACCGTTAACCTGGCTTCATACCTTGCCACATAGCCAAACCCATCCTTTCCCGCTACAGTTAATTTCTTGTGGCGCGAAAGGAGGCAAAAATGCTCTATATCt                  <  1:502644/139‑1 (MQ=255)                                                                                                              aaTCATTCTTCCCCAGGAAGACGCGTAAATTTTTGTGACCGTTAACCTGGCTTCATACCTTGCCACATAGCCAAACCCATCCTTTCCCGCTACAGTTAATTTCTTGTGGCGCGAAAGGAGGCAAAAATGCTCTATATCt                  <  1:191854/139‑1 (MQ=255)                                                                                                              aaTCATTCTTCCCCAGGAAGACGCGTAAATTTTTGTGACCGTTAACCTGGCTTCATACCTTGCCACATAGCCAAACCCATCCTTTCCCGCTACAGTTAATTTCTTGTGGCGCGAAAGGAGGCAAAAATGCCCTATATCt                  <  2:324772/139‑1 (MQ=255)                                                                                                                       tcCCCAGGAAGACGCGTAAATTTTTGTGACCGTTAACCTGGCTTCATACCTTGCCACATAGCCAAACCCATCCTTTCCCGCTACAGTTAATTTCTTGTGGCGCGAAAGGAGGCAAAAATg                            >  2:31226/1‑120 (MQ=255)                                                                                                                       tcCCCAGGAAGACGCGTAAATTTTTGTGACCGTTAACCTGGCTTCATACCTTGCCACATAGCCAAACCCATCCTTTCCCGCTACAGTTAATTTCTTGTGGCGCGAAAGGAGGCAAAAATg                            <  1:31226/120‑1 (MQ=255)                                                                                                                           cAGGAAGACGCGTAAATTTTTGTGACCGTTAACCTGGCTTCATACCTTGCCACATAGCCAAACCCATCCTTTCCCGCTACAGTTAATTTCTTGTGGCGCGAAAGGAGGCAAAAATGCTCTATATCTTTGATTTAGGTaa     <  1:34677/139‑1 (MQ=255)                                                                                                                              gAAGACGCGTAAATTTTTGTGACCGTTAACCTGGCTTCATACCTTGCCACATAGCCAAACCCATCCTTTCCCGCTACAGTTAATTTCTTGTGGCGCGAAAGGAGGCaaaaa                              >  2:781717/1‑111 (MQ=255)                                                                                                                              gAAGACGCGTAAATTTTTGTGACCGTTAACCTGGCTTCATACCTTGCCACATAGCCAAACCCATCCTTTCCCGCTACAGTTAATTTCTTGTGGCGCGAAAGGAGGCaaaaa                              <  1:781717/111‑1 (MQ=255)                                                                                                                              gAAGACGCGTAAATTTTTGTGACCGTTAACCTGGCTTCATACCTTGCCACATAGCCAAACCCATCCTTTCCCGCTACAGTTAATTTCTTGTGGCGCGAAAGGAGGCAAAAATGCTCTATATCTTTGATTTAGGTAAtgt  >  1:181980/1‑139 (MQ=255)                                                                                                                                   cgcgTAAATTTTTGTGACCGTTAACCTGGCTTCATACCTTGCCACATAGCCAAACCCATCCTTTCCCGCTACAGTTAATTTCTTGTGGCGCGAAAGGAGGCAAAAATGCTCTATATCt                  >  2:550277/1‑118 (MQ=255)                                                                                                                                   cgcgTAAATTTTTGTGACCGTTAACCTGGCTTCATACCTTGCCACATAGCCAAACCCATCCTTTCCCGCTACAGTTAATTTCTTGTGGCGCGAAAGGAGGCAAAAATGCTCTATATCt                  <  1:550277/118‑1 (MQ=255)                                                                                                                                         |                                                                                                                                ATGCTTCGGCGGCGGTTGAAATTGGTAACGTGGCACAGGTCACTGCGCTCTTTACCGACGAGCTGCCGCCCGCTGCGCTAAAATCACGCTTACAAGACAGCCAAATTGAAATCATTCTTCCCCAGGAAGACGCGTAGATTTTTGTGACCGTTAACCTGGCTTCATACCTTGCCACATAGCCAAACCCATCCTTTCCCGCTACAGTTAATTTCTTGTGGCGCGAAAGGAGGCAAAAATGCTCTATATCTTTGATTTAGGTAATGT  >  NC_000913/4075318‑4075581

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