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%	123.1 / NA	39	*262* (TAG→TAA)	yihW	putative transcriptional regulator for sulphoquinovose utilization
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (19/20);  total (19/20)

ATGCTTCGGCGGCGGTTGAAATTGGTAACGTGGCACAGGTCACTGCGCTCTTTACCGACGAGCTGCCGCCCGCTGCGCTAAAATCACGCTTACAAGACAGCCAAATTGAAATCATTCTTCCCCAGGAAGACGCGTAGATTTTTGTGACCGTTAACCTGGCTTCATACCTTGCCACATAGCCAAACCCATCCTTTCCCGCTACAGTTAATTTCTTGTGGCGCGAAAGGAGGCAAAAATGCTCTATATCTTTGATTTAGGTAATGTGA  >  NC_000913/4075318‑4075583                                                                                                                                         |                                                                                                                                  aTGCTTCGGCGGCGGTTGAAATTGGTAACGTGGCACAGGTCACTGCGCTCTTTACCGACGAGCTGCCGCCCGCTGCGCTAAAATCACGCTTACAAGACAGCCAAATTGAAATCATTCTTCCCCAGGAAGACGCGTAAAt                                                                                                                                 >  1:698289/1‑139 (MQ=255)     ttCGGCGGCGGTTGAAATTGGTAACGTGGCACAGGTCACTGCGCTCTTTACCGACGAGCTGCCGCCCGCTGCGCTAAAATCACGCTTACAAGACAGCCAAATTGAAATCATTCTTCCCCAGGAAGACGCGTAAAttttt                                                                                                                             >  2:129224/1‑139 (MQ=255)       cggcggcggTTGAAATTGGTAACGTGGCACAGGTCACTGCGCTCTTTACCGACGAGCTGCCGCCCGCTGCGCTAAAATCACGCTTACAAGACAGCCAAATTGAAATCATTCTTCCCCAGGAAGACGCGTAAATTTTtgt                                                                                                                           <  2:999052/139‑1 (MQ=255)             cggTTGAAATTGGTAACGTGGCACAGGTCACTGCGCTCTTTACCGACGAGCTGCCGCCCGCTGCGCTAAAATCACGCTTACAAGACAGCCAAATTGAAATCATTCTTCCCCAGGAAGACGCGTAAATTTTTGTGACCGt                                                                                                                     >  2:1012677/1‑139 (MQ=255)                  gAAATTGGTAACGTGGCACAGGTCACTGCGCTCTTTACCGACGAGCTGCCGCCCGCTGCGCTAAAATCACGCTTACAAGACAGCCAAATTGAAATCATTCTTCCCCAGGAAGACGCGTAAATTTTTGTGACCGTTAAcc                                                                                                                <  1:811515/139‑1 (MQ=255)                  gAAATTGGTAACGTGGCACAGGTCACTGCGCTCTTTACCGACGAGCTGCCGCCCGCTGCGCTAAAATCACGCTTACAAGACAGCCAAATTGAAATCATTCTTCCCCAGGAAGACGCGTAAATTTTTGTGACCGTTAAcc                                                                                                                <  1:1023606/139‑1 (MQ=255)                   aaaTTGGTAACGTGGCACAGGTCACTGCGCTCTTTACCGACGAGCTGCCGCCCGCTGCGCTAAAATCACGCTTACAAGACAGCCAAATTGAAATCATTCTTCCCCAGGAAGACGCGTAAATTTTTGTGACCGTTAACCt                                                                                                               <  2:159105/139‑1 (MQ=255)                          tAACGTGGCACAGGTCACTGCGCTCTTTACCGACGAGCTGCCGCCCGCTGCGCTAAAATCACGCTTACAAGACAGCCAAATTGAAATCATTCTTCCCCAGGAAGACGCGTAAATTTTTGTGACCGTTAACCTGGCTTCa                                                                                                        <  1:795198/139‑1 (MQ=255)                          tAACGTGGCACAGGTCACTGCGCTCTTTACCGACGAGCTGCCGCCCGCTGCGCTAAAATCACGCTTACAAGACAGCCAAATTGAAATCATTCTTCCCCAGGAAGACGCGTAAATTTTTGTGACCGTTAACCTGGCTTCa                                                                                                        >  2:305060/1‑139 (MQ=255)                             cGTGGCACAGGTCACTGCGCTCTTTACCGACGAGCTGCCGCCCGCTGCGCTAAAATCACGCTTACAAGACAGCCAAATTGAAATCATTCTTCCCCAGGAAGACGCGTAAATTTTTGTGACCGTTAACCTGGCTTCATAc                                                                                                     >  2:866147/1‑139 (MQ=255)                               tGGCACAGGTCACTGCGCTCTTTACCGACGAGCTGCCGCCCGCTGCGCTAAAATCACGCTTACAAGACAGCCAAATTGAAATCATTCTTCCCCAGGAAGACGCGTAAATTTTTGTGACCGTTAACCTGGCTTCATACCt                                                                                                   >  1:124881/1‑139 (MQ=255)                                 gCACAGGTCACTGCGCTCTTTACCGACGAGCTGCCGCCCGCTGCGCTAAAATCACGCTTACAAGACAGCCAAATTGAAATCATTCTTCCCCAGGAAGACGCGTAAATTTTTGTGACCGTTAACCTGGCTTCATACCTTg                                                                                                 >  1:39635/1‑139 (MQ=255)                                                     tACCGACGAGCTGCCGCCCGCTGCGCTAAAATCACGCTTACAAGACAGCCAAATTGAAATCATTCTTCCCCAGGAAGACGCGTAAATTTTTGTGACCGTTAACCTGGCTTCATACCTTGCCACATAGCCAAACCCATcc                                                                             <  1:736425/139‑1 (MQ=255)                                                      aCCGACGAGCTGCCGCCCGCTGCGCTAAAATCACGCTTACAAGACAGCCAAATTGAAATCATTCTTCCCCAGGAAGACGCGTAAATTTTTGTGACCGTTAACCTGGCTTCATACCTTGCCACATAGCCAAACCCATCCt                                                                            >  2:881391/1‑139 (MQ=255)                                                             aGCTGCCGCCCGCTGCGCTAAAATCACGCTTACAAGACAGCCAAATTGAAATCATTCTTCCCCAGGAAGACGCGTAAATTTTTGTGACCGTTAACCTGGCTTCATACCTTGCca                                                                                              >  2:453336/1‑114 (MQ=255)                                                             aGCTGCCGCCCGCTGCGCTAAAATCACGCTTACAAGACAGCCAAATTGAAATCATTCTTCCCCAGGAAGACGCGTAAATTTTTGTGACCGTTAACCTGGCTTCATACCTTGCca                                                                                              <  1:453336/114‑1 (MQ=255)                                                             aGCTGCCGCCCGCTGCGCTAAAATCACGCTTACAAGACAGCCAAATTGAAATCATTCTTCCCCAGGAAGACGCGTAAATTTTTGTGACCGTTAACCTGGCTTCATACCTTGCCACATAGCCAAACCCATCCTTTCCCGc                                                                     >  2:190668/1‑139 (MQ=255)                                                              gCTGCCGCCCGCTGCGCTAAAATCACGCTTACAAGACAGCCAAATTGAAATCATTCTTCCCCAGGAAGACGCGTAAATTTTTGTGACCGTTAACCTGGCTTCATACCTTGCCACATAGCCAAACCCATCCTTTCCCGCt                                                                    <  1:190668/139‑1 (MQ=255)                                                               cTGCCGCCCGCTGCGCTAAAATCACGCTTACAAGACAGCCAAATTGAAATCATTCTTCCCCAGGAAGACGCGTAAATTTTTGTGACCGTTAACCTGGCTTCATACCTTGCCACATAGCCAAACCCATCCTTTCCCGCTa                                                                   <  2:336042/139‑1 (MQ=255)                                                                tGCCGCCCGCTGCGCTAAAATCACGCTTACAAGACAGCCAAATTGAAATCATTCTTCCCCAGGAAGACGCGTAAATTTTTGTGACCGTTAACCTGGCTTCATACCTTGCCACATAGCCAAACCCATCCTTTCCCGCTAc                                                                  <  1:382284/139‑1 (MQ=255)                                                                         cTGCGCTAAAATCACGCTTACAAGACAGCCAAATTGAAATCATTCTTCCCCAGGAAGACGCGTAAATTTTTGTGACCGTTAACCTGGCTTCATTCCTTGCCACATAGCCAAACCCACCCTTTCCCGCTACAGTTAAttt                                                         >  2:474855/1‑139 (MQ=255)                                                                                      aCGCTTACAAGACAGCCAAATTGAAATCATTCTTCCCCAGGAAGACGCGTAAATTTTTGTGACCGTTAACCTGGCTTCATACCTTGCCACATAGCCAAACCCATCCTTTCCCGCTACAGTTAATTTCTTGTGGCGCGaa                                            >  2:593494/1‑139 (MQ=255)                                                                                         cTTACAAGACAGCCAAATTGAAATCATTCTTCCCCAGGAAGACGCGTAAATTTTTGTGACCGTTAACCTGGCTTCATACCTTGCCACATAGCCAAACCCATCCTTTCCCGCTACAGTTAATTTCTTGTGGCGCGAAagg                                         <  2:124881/139‑1 (MQ=255)                                                                                                   aGCCAAATTGAAATCATTCTTCCCCAGGAAGACGCGTAAATTTTTGTGACCGTTAACCTGGCTTCATACCTTGCCACATAGCCAAACCCATCCTTTCCCGCTACAGTTAATTTCTTGTGGCGCGAAAGGAGGCAAAAAt                               >  1:142195/1‑139 (MQ=255)                                                                                                    gCCAAATTGAAATCATTCTTCCCCAGGAAGACGCGTAAATTTTTGTGACCGTTAACCTGGCTTCATACCTTGCCACATAGCCAAACCCATCCTTTCCCGCTACAGTTAATTTCTTGTGGCGCGAAAGGAGGCAAAAATg                              >  1:209460/1‑139 (MQ=255)                                                                                                    gCCAAATTGAAATCATTCTTCCCCAGGAAGACGCGTAAATTTTTGTGACCGTTAACCTGGCTTCATACCTTGCCACATAGCCAAACCCATCCTTTCCCGCTACAGTTAATTTCTTGTGGCGCGAAAGGAGGCAAAAATg                              >  1:554298/1‑139 (MQ=255)                                                                                                     ccAAATTGAAATCATTCTTCCCCAGGAAGACGCGTAAATTTTTGTGACCGTTAACCTGGCTTCATACCTTGCCACATAGc                                                                                        >  1:508355/1‑80 (MQ=255)                                                                                                     ccAAATTGAAATCATTCTTCCCCAGGAAGACGCGTAAATTTTTGTGACCGTTAACCTGGCTTCATACCTTGCCACATAGc                                                                                        <  2:508355/80‑1 (MQ=255)                                                                                                       aaaTTGAAATCATTCTTCCCCAGGAAGACGCGTAAATTTTTGTGACCGTTAACCTGGCTTCATACCTTGCCACATAGCCAAACCCATCCTTTCCCGCTACAGTTAATTTCTTGTGGCGCGAAAGGAGGCAAAAATGctc                           <  1:345042/139‑1 (MQ=255)                                                                                                          ttGAAATCATTCTTCCCCAGGAAGACGCGTAAATTTTTGTGACCGTTAACCTGGCTTCATACCTTGCCACATAGCCAAACCCATCCTTTCCCGCTACAGTTAATTTCTTGTGGCGCGAAAGGAGGCAAAAATGCTCtat                        <  1:305060/139‑1 (MQ=255)                                                                                                              aaTCATTCTTCCCCAGGAAGACGCGTAAATTTTTGTGACCGTTAACCTGGCTTCATACCTTGCCACATAGCCAAACCCATCCTTTCCCGCTACAGTTAATTTCTTGTGGCGCGAAAGGAGGCAAAAATGCTCTATATCt                    <  1:691642/139‑1 (MQ=255)                                                                                                              aaTCATTCTTCCCCAGGAAGACGCGTAAATTTTTGTGACCGTTAACCTGGCTTCATACCTTGCCACATAGCCAAACCCATCCTTTCCCGCTACAGTTAATTTCTTGTGGCGCGAAAGGAGGCAAAAATGCTCTATATCt                    <  1:1040956/139‑1 (MQ=255)                                                                                                                        ccccAGGAAGACGCGTAAATTTTTGTGACCGTTAACCTGGCTTCATACCTTGCCACATAGCCAAACCCATCCTTTCCCGCTACAGTTAATTTCTTGTGGCGCGAAAGGAGGCAAAAATGCTCTATATCTTTGATTTAgg          <  1:1012677/139‑1 (MQ=255)                                                                                                                            aGGAAGACGCGTAAATTTTTGTGACCGTTAACCTGGCTTCATACCTTGCCACATAGCCAAACCCATCCTTTCCCGCTACAGTTAATTTCTTGTGGCGCGAAAGGAGGCAAAAATGCTCTATATCTTTGATTTAGGTAAt      <  2:209460/139‑1 (MQ=255)                                                                                                                             ggAAGACGCGTAAATTTTTGTGACCGTTAACCTGGCTTCATACCTTGCCACATAGCCAAACCCATCCTTTCCCGCTACAGTTAATTTCTTGTGGCGCGAAAGGAGGCAAAAATGCTCTATATCTTTGATTTAGGTAAtg     >  2:837429/1‑139 (MQ=255)                                                                                                                                aggcgcgGAAATTGTTGTGACCGTTAACCTGGCTTCATACCTTGCCACATCGCCAATCCCTTCCTTTCCCGCTACATTTAATTTCTTGTGGCGCGAAAGGATTCAAGATGGCTCTATATCTTTGATTTAGGTAATGTGa  <  2:237960/136‑1 (MQ=255)                                                                                                                                   cgcgTAAATTTTTGTGACCGTTAACCTGGCTTCATACCTTGCCACATAGCCAAACCCAt                                                                               <  2:475970/59‑1 (MQ=255)                                                                                                                                   cgcgTAAATTTTTGTGACCGTTAACCTGGCTTCATACCTTGCCACATAGCCAAACCCAt                                                                               >  1:475970/1‑59 (MQ=255)                                                                                                                                   cgcgTAAATTTTTGTGACCGTTAACCTGGCTTCATACCTTGCCACATAGCCAAACCCATCCTTTCCCGCTACAGTTAATTTCTTGTGGCGCGAAAGGAGGCaaaa                                 >  1:64686/1‑105 (MQ=255)                                                                                                                                   cgcgTAAATTTTTGTGACCGTTAACCTGGCTTCATACCTTGCCACATAGCCAAACCCATCCTTTCCCGCTACAGTTAATTTCTTGTGGCGCGAAAGGAGGCaaaa                                 <  2:64686/105‑1 (MQ=255)                                                                                                                                         |                                                                                                                                  ATGCTTCGGCGGCGGTTGAAATTGGTAACGTGGCACAGGTCACTGCGCTCTTTACCGACGAGCTGCCGCCCGCTGCGCTAAAATCACGCTTACAAGACAGCCAAATTGAAATCATTCTTCCCCAGGAAGACGCGTAGATTTTTGTGACCGTTAACCTGGCTTCATACCTTGCCACATAGCCAAACCCATCCTTTCCCGCTACAGTTAATTTCTTGTGGCGCGAAAGGAGGCAAAAATGCTCTATATCTTTGATTTAGGTAATGTGA  >  NC_000913/4075318‑4075583

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