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%	124.9 / 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 (20/19);  total (20/19)

CGGCGGCGGTTGAAATTGGTAACGTGGCACAGGTCACTGCGCTCTTTACCGACGAGCTGCCGCCCGCTGCGCTAAAATCACGCTTACAAGACAGCCAAATTGAAATCATTCTTCCCCAGGAAGACGCGTAGATTTTTGTGACCGTTAACCTGGCTTCATACCTTGCCACATAGCCAAACCCATCCTTTCCCGCTACAGTTAATTTCTTGTGGCGCGAAAGGAGGCAAAAATGCTCTATATCTTTGATTTAGGTAATGTGATTG  >  NC_000913/4075324‑4075586                                                                                                                                   |                                                                                                                                     cggcggcggTTGAAATTGGTAACGTGGCACAGGTCACTGCGCTCTTTACCGACGAGCTGCCGCCCGCTGCGCTAAAATCACGCTTACAAGACAGCCAAATTGAAATCATTCTTCCCCAGGAAGACGCGTAAATTTTtgt                                                                                                                              >  2:714833/1‑139 (MQ=255) cggcggcggTTGAAATTGGTAACGTGGCACAGGTCACTGCGCTCTTTACCGACGAGCTGCCGCCCGCTGCGCTAAAATCACGCTTACAAGACAGCCAAATTGAAATCATTCTTCCCCAGGAAGACGCGTAAATTTTtgt                                                                                                                              >  1:91161/1‑139 (MQ=255)         gTTGAAATTGGTAACGTGGCACAGGTCACTGCGCTCTTTACCGACGAGCTGCCGCCCGCTGCGCTAAAATCACGCTTACAAGACAGCCAAATTGAAATCATTCTTCCCCAGGAAGACGCGTAAATTTTtgtg                                                                                                                             <  2:530798/132‑1 (MQ=255)         gTTGAAATTGGTAACGTGGCACAGGTCACTGCGCTCTTTACCGACGAGCTGCCGCCCGCTGCGCTAAAATCACGCTTACAAGACAGCCAAATTGAAATCATTCTTCCCCAGGAAGACGCGTAAATTTTtgtg                                                                                                                             >  1:530798/1‑132 (MQ=255)             aaaTTGGTAACGTGGCACAGGTCACTGCGCTCTTTACCGACGAGCTGCCGCCCGCTGCGCTAAAATCACGCTTACAAGACAGCCAAATTGAAATCATTCTTCCCCAGGAAGACGCGTAAATTTTTGTGACCGTTAACCt                                                                                                                  <  2:91161/139‑1 (MQ=255)                   gTAACGTGGCACAGGTCACTGCGCTCTTTACCGACGAGCTGCCGCCCGCTGCGCTAAAATCACGATTACAAGACAGCCAAATTGAAATCATTCTTCCCCAGGAAGACGCGTAAATTTTTGTGACCGTTAACCTGGCTTc                                                                                                            >  2:586124/1‑139 (MQ=255)                       cGTGGCACAGGTCACTGCGCTCTTTACCGACGAGCTGCCGCCCGCTGCGCTAAAATCACGCTTACAAGACAGCCAAATTGAAATCATTCTTCCCCAGGAAGACGCGTAAATTTTTGTGACCGTTAACCTGGCTTCATAc                                                                                                        <  1:113424/139‑1 (MQ=255)                           gCACAGGTCACTGCGCTCTTTACCGACGAGCTGCCGCCCGCTGCGCTAAAATCACGCTTACAAGACAGCCAAATTGAAATCATTCTTCCCCAGGAAGACGCGTAAATTTTTGTGACCGTTAACCTGGCTTCATACCTTg                                                                                                    >  2:226740/1‑139 (MQ=255)                               aGGTCACTGCGCTCTTTACCGACGAGCTGCCGCCCGCTGCGCTAAAATCACGCTTACAAGACAGCCAAATTGAAATCATTCTTCCCCAGGAAGACGCGTAAATTTTTGTGACCGTTAACCTGGCTTCATATCTTGCcac                                                                                                <  1:714833/139‑1 (MQ=255)                                 gTCACTGCGCTCTTTACCGACGAGCTGCCGCCCGCTGCGCTAAAATCACGCTTACAAGACAGCCAAATTGAAATCATTCTTCCCCAGGAAGACGCGTAAATTTTTGTGACCGTTAACCTGGCTTCATACCTTGCCACAt                                                                                              >  1:92852/1‑139 (MQ=255)                                   cACTGCGCTCTTTACCGACGAGCTGCCGCCCGCTGCGCTAAAATCACGCTTACAAGACAGCCAAATTGAAATCATTCTTCCCCAGGAAGACGCGTAAATTTTTGTGACCGTTAACCTGGCTTCATa                                                                                                         >  1:425039/1‑126 (MQ=255)                                   cACTGCGCTCTTTACCGACGAGCTGCCGCCCGCTGCGCTAAAATCACGCTTACAAGACAGCCAAATTGAAATCATTCTTCCCCAGGAAGACGCGTAAATTTTTGTGACCGTTAACCTGGCTTCATa                                                                                                         <  2:425039/126‑1 (MQ=255)                                        cgcTCTTTACCGACGAGCTGCCGCCCGCTGCGCTAAAATCACGCTTACAAGACAGCCAAATTGAAATCATTCTTCCCCAGGAAGACGCGTAAATTTTTGTGACCGTTAACCTGGCTTCATACCTTGCCACATAGCCaaa                                                                                       >  2:509255/1‑139 (MQ=255)                                           tctTTACCGACGAGCTGCCGCCCGCTGCGCTAAAATCACGCTTACAAGACAGCCAAATTGAAATCATTCTTCCCCAGGAAGACGCGTAAATTTTTGTGACCGTTAACCTGGCTTCATACCTTGCCACATAGCCAAAccc                                                                                    >  2:14486/1‑139 (MQ=255)                                               tACCGACGAGCTGCCGCCCGCTGCGCTAAAATCACGCTTACAAGACAGCCAAATTGAAATCATTCTTCCCCAGGAAGACGCGTAAATTTTTGTGACCGTTAACCTGGCTTCATACCTTGCCACATAGCCAAACCCATcc                                                                                <  2:248009/139‑1 (MQ=255)                                                       aGCTGCCGCCCGCTGCGCTAAAATCACGCTTACAAGACAGCCAAATTGAAATCATTCTTCCCCAGGAAGACGCGTAAATTTTTGTGACCGTTAACCTGGCTTCATACCTTGCCACATAGCCAAACCCATCCTTTCCCGc                                                                        >  2:586208/1‑139 (MQ=255)                                                       aGCTGCCGCCCGCTGCGCTAAAATCACGCTTACAAGACAGCCAAATTGAAATCATTCTTCCCCAGGAAGACGCGTAAATTTTTGTGACCGTTAACCTGGCTTCATACCTTGCCACATAGCCAAACCCATCCTTTCCCGc                                                                        >  2:118636/1‑139 (MQ=255)                                                       aGCTGCCGCCCGCTGCGCTAAAATCACGCTTACAAGACAGCCAAATTGAAATCATTCTTCCCCAGGAAGACGCGTAAATTTTTGTGACCGTTAACCTAGCTTCATACCTTGCCACATAGCCAAACCCATCCTTTCCCGc                                                                        >  2:462109/1‑139 (MQ=255)                                                         cTGCCGCCCGCTGCGCTAAAATCACGCTTACAAGACAGCCAAATTGAAATCATTCTTCCCCAGGAAGACGCGTAAATTTTTGTGACCGTTAACCTGGCTTCATACCTTGCCACATAGCCAAACCCATCCTTTCCCGCTa                                                                      >  1:598309/1‑139 (MQ=255)                                                                        cTAAAATCACGCTTACAAGACAGCCAAATTGAAATCATTCTTCCCCAGGAAGACGCGTAAATTTTTGTGACCGTTAACCTGGCTTCATACCTTGCCACATAGCCAAACCCATCCTTTCCCGCTACAGTTAATTTCTtgt                                                       <  2:598309/139‑1 (MQ=255)                                                                                aCGCTTACAAGACAGCCAAATTGAAATCATTCTTCCCCAGGAAGACGCGTAAATTTTTGTGACCGTTAACCTGGCTTCATACCTTGCCACATAGCCAAACCCATCCTTTCCCGCTa                                                                      <  2:588890/116‑1 (MQ=255)                                                                                aCGCTTACAAGACAGCCAAATTGAAATCATTCTTCCCCAGGAAGACGCGTAAATTTTTGTGACCGTTAACCTGGCTTCATACCTTGCCACATAGCCAAACCCATCCTTTCCCGCTa                                                                      >  1:588890/1‑116 (MQ=255)                                                                                       cAAGACAGCCAAATTGAAATCATTCTTCCCCAGGAAGACGCGTAAATTTTTGTGACCGTTAACCTGGCTTCATACCTTGCCACATAGCCAAACCCATCCTTTCCCGCTACAGTTAATTTCTTGTGGCGCGAAAGGAGGc                                        <  2:700183/139‑1 (MQ=255)                                                                                           aCAGCCAAATTGAAATCATTCTTCCCCAGGAAGACGCGTAAATTTTTGTGACCGTTAACCTGGCTTCATACCTTGTCACATAGCCAAACCCATCCTTTCCCGCTACAGTTAATTTCTTGTGGCGCGAAAGGAGGCaaaa                                    <  2:377045/139‑1 (MQ=255)                                                                                           aCAGCCAAATTGAAATCATTCTTCCCCAGGAAGACGCGTAAATTTTTGTGACCGTTAACCTGGCTTCATACCTTGCCACATAGCCAAACCCATCCTTTCCCGCTACAGTTAATTTCTTGTGGCGCGAAAGGAGGCaaaa                                    <  1:671727/139‑1 (MQ=255)                                                                                           aCAGCCAAATTGAAATCATTCTTCCCCAGGAAGACGCGTAAATTTTTGTGACCGTTAACCTGGCTTCATACCTTGCCACATAGCCAAACCCATCCTTTCCCGCTACAGTTAATTTCTTGTGGCGCGAAAGGAGGCaaaa                                    <  1:586208/139‑1 (MQ=255)                                                                                            cAGCCAAATTGAAATCATTCTTCCCCAGGAAGACGCGTAAATTTTTGTGACCGTTAACCTGGCTTCATACCTTGCCACATAGCCAAACCCATCCTTTCCCGCTACAGTTAATTTCTTGTGGCGCGAAAGGAGGCaaaaa                                   <  1:14486/139‑1 (MQ=255)                                                                                                  aaTTGAAATCATTCTTCCCCAGGAAGACGCGTAAATTTTTGTGACCGTTAACCTGGCTTCATACCTTGCCACATAgtca                                                                                         <  2:396067/79‑4 (MQ=255)                                                                                                  aaTTGAAATCATTCTTCCCCAGGAAGACGCGTAAATTTTTGTGACCGTTAACCTGGCTTCATACCTTGCCACATAgtca                                                                                         >  1:396067/1‑76 (MQ=255)                                                                                                        aaTCATTCTTCCCCAGGAAGACGCGTAAATTTTTGTGACCGTTAACCTGGCTTCATACCTTGCCACATAGCCAAACCCATCCTTTCCCGCTACAGTTAATTTCTTGTGGCGCGAAAGGAGGCAAAAATGCTCTATATCt                       <  1:226740/139‑1 (MQ=255)                                                                                                        aaTCATTCTTCCCCAGGAAGACGCGTAAATTTTTGTGACCGTTAACCTGGCTTCATACCTTGCCACATAGCCAAACCCATCCTTTCCCGCTACAGTTAATTTCTTGTGGCGCGAAAGGAGGCAAAAATGCTCTATATCt                       <  2:749608/139‑1 (MQ=255)                                                                                                              tcttcCCCAGGAAGACGCGTAAATTTTTCTGACCGTTAACCTGGCTTCATACCTTGCCACATAGCCAAACCCATCCTTTCCCGCTACAGTTAATTTCTTGTGGCGCGAAAGGAGGCAAAAATGCTCTATATCTTTGAtt                 <  1:462109/139‑1 (MQ=255)                                                                                                                 tcCCCAGGAAGACGCGTAAATTTTTGTGACCGTTAACCTGGCTTCATACCTTGCCACATAGCCAAACCCATc                                                                                 >  2:601554/1‑72 (MQ=255)                                                                                                                 tcCCCAGGAAGACGCGTAAATTTTTGTGACCGTTAACCTGGCTTCATACCTTGCCACATAGCCAAACCCATc                                                                                 <  1:601554/72‑1 (MQ=255)                                                                                                                        gAAGACGCGTAAATTTTTGTGACCGTTAACCTGGCTTCATACCTTGCCACATAGTCAAACCCATCCTTTCCCGCTACAGTTAATTTCTTGTGGCGCGAAAGGAGGCAAAAATGCTCTATATCTTTGATTTAGGTAAtgt       >  1:619280/1‑139 (MQ=255)                                                                                                                        gAAGACGCGTAAATTTTTGTGACCGTTAACCTGGCTTCATACCTTGCCACATAGCCAAACCCATCCTTTCCCGCTACAGTTAATTTCTTGTGGCGCGAAAGGAGGCAAAAATGCTCTATATCTTTGATTTAGGTAAtgt       >  2:330693/1‑139 (MQ=255)                                                                                                                            aCGCGTAAATTTTTGTGACCGTTAACCTGGCTTCATACCTTGCCACATAGCCAAACCCAt                                                                                  >  1:325860/1‑60 (MQ=255)                                                                                                                            aCGCGTAAATTTTTGTGACCGTTAACCTGGCTTCATACCTTGCCACATAGCCAAACCCAt                                                                                  <  2:325860/60‑1 (MQ=255)                                                                                                                             cgcgTAAATTTTTGTGACCGTTAACCTGGCTTCATACCTTGCCACATAGCCAAACCCATCCTTTCCCGCTACAGTTAATTTCTTGTGGCGCGAAAGGAGGCAAAAATGCTCTATATCTTTGATTTAGGTAATGTGATTg  >  1:267638/1‑139 (MQ=255)                                                                                                                                   |                                                                                                                                     CGGCGGCGGTTGAAATTGGTAACGTGGCACAGGTCACTGCGCTCTTTACCGACGAGCTGCCGCCCGCTGCGCTAAAATCACGCTTACAAGACAGCCAAATTGAAATCATTCTTCCCCAGGAAGACGCGTAGATTTTTGTGACCGTTAACCTGGCTTCATACCTTGCCACATAGCCAAACCCATCCTTTCCCGCTACAGTTAATTTCTTGTGGCGCGAAAGGAGGCAAAAATGCTCTATATCTTTGATTTAGGTAATGTGATTG  >  NC_000913/4075324‑4075586

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