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

ATGCTTCGGCGGCGGTTGAAATTGGTAACGTGGCACAGGTCACTGCGCTCTTTACCGACGAGCTGCCGCCCGCTGCGCTAAAATCACGCTTACAAGACAGCCAAATTGAAATCATTCTTCCCCAGGAAGACGCGTAGATTTTTGTGACCGTTAACCTGGCTTCATACCTTGCCACATAGCCAAACCCATCCTTTCCCGCTACAGTTAATTTCTTGTGGCGCGAAAGGAGGCAAAAATGCTCTATATCTTTGATTTAGGTAATGTGATTG  >  NC_000913/4075318‑4075586                                                                                                                                         |                                                                                                                                     aTGCTTCGGCGGCGGTTGAAATTGGTAACGTGGCACAGGTCACTGCGCTCTTTACCGACGAGCTGCCGCCCGCTGCGCTAAAATCACGCTTACAAGACAGCCAAATTGAAATCATTCTTCCCCAGGAAGACGCGTAAAt                                                                                                                                    >  2:1454610/1‑139 (MQ=255) aTGCTTCGGCGGCGGTTGAAATTGGTAACGTGGCACAGGTCACTGCGCTCTTTACCGACGAGCTGCCGCCCGCTGCGCTAAAATCACGCTTACAAGACAGCCAAATTGAAATCATTCTTCCCCAGGAAGACGCGTAAAt                                                                                                                                    >  1:619224/1‑139 (MQ=255)   gCTTCGGCGGCGGTTGAAATTGGTAACGTGGCACAGGTCACTGCGCTCTTTACCGACGAGCTGCCGCCCGCTGCGCTAAAATCACGCTTACAAGACAGCCAAATTGAAATCATTCTTCCCCAGGAAGACGCGTAAAttt                                                                                                                                  >  2:386946/1‑139 (MQ=255)     ttCGGCGGCGGTTGAAATTGGTAACGTGGCACAGGTCACTGCGCTCTTTACCGACGAGCTGCCGCCCGCTGCGCTAAAATCACGCTTACAAGACAGCCAAATTGAAATCATTCTTCCCCAGGAAGACGCGTAAAttttt                                                                                                                                >  1:759538/1‑139 (MQ=255)       cggcggcggTTGAAATTGGTAACGTGGCACAGGTCACTGCGCTCTTTACCGACGAGCTGCCGCCCGCTGCGCTAAAATCACGCTTACAAGACAGCCAAATTGAAATCATTCTTCCCCAGGAAGACGCGTAAATTTTtgt                                                                                                                              <  2:1226713/139‑1 (MQ=255)       cggcggcggTTGAAATTGGTAACGTGGCACAGGTCACTGCGCTCTTTACCGACGAGCTGCCGCCCGCTGCGCTAAAATCACGCTTACAAGACAGCCAAATTGAAATCATTCTTCCCCAGGAAGACGCGTAAATTTTtgt                                                                                                                              <  2:759538/139‑1 (MQ=255)       cggcggcggTTGAAATTGGTAACGTGGCACAGGTCACTGCGCTCTTTACCGACGAGCTGCCGCCCGCTGCGCTAAAATCACGCTTACAAGACAGCCAAATTGAAATCATTCTTCCCCAGGAAGACGCGTAAATTTTtgt                                                                                                                              >  1:371619/1‑139 (MQ=255)                  gAAATTGGTAACGTGGCACAGGTCACTGCGCTCTTTACCGACGAGCTGCCGCCCGCTGCGCTAAAATCACGCTTACAAGACAGCCAGATTGAAATCATTCTTCCCCAGGAAGACGCGTAAATTTTTGCGACCGTTAccc                                                                                                                   <  1:850653/139‑1 (MQ=255)                  gAAATTGGTAACGTGGCACAGGTCACTGCGCTCTTTACCGACGAGCTGCCGCCCGCTGCGCTAAAATCACGCTTACAAGACAGCCAAATTGAAATCATTCTTCCCCAGGAAGACGCGTAAATTTTTGTGACCGTTAAcc                                                                                                                   <  1:1091131/139‑1 (MQ=255)                  gAAATTGGTAACGTGGCACAGGTCACTGCGCTCTTTACCGACGAGCTGCCGCCCGCTGCGCTAAAATCACGCTTACAAGACAGCCAAATTGAAATCATTCTTCCCCAGGAAGACGCGTAAATTTTTGTGACCGTTAAcc                                                                                                                   <  1:1454610/139‑1 (MQ=255)                                     aGGTCACTGCGCTCTTTACCGACGAGCTGCCGCCCGCTGCGCTAAAATCACGCTTACAAGACAGCCAAATTGAAATCATTCTTCCCCAGGAAGACGCGTAAATTTTTGTGACCGTTAACCTGGCTTCATACCTTGCcac                                                                                                <  1:1392947/139‑1 (MQ=255)                                          aCTGCGCTCTTTACCGACGAGCTGCCGCCCGCTGCGCTAAAATCACGCTTACAAGACAGCCAAATTGAAATCATTCTTCCCCAGGAAGACGCGTAAATTTTTGTGACCGTTAACCTGGCTTCATACCTTGCCACATAGc                                                                                           >  1:965413/1‑139 (MQ=255)                                          aCTGCGCTCTTTACCGACGAGCTGCCGCCCGCTGCGCTAAAATCACGCTTACAAGACAGCCAAATTGAAATCATTCTTCCCCAGGAAGACGCGTAAATTTTTGTGACCGTTAACCTGGCTTCATACCTTGCCACATAGc                                                                                           <  2:371619/139‑1 (MQ=255)                                          aCTGCGCTCTTTACCGACGAGCTGCCGCCCGCTGCGCTAAAATCACGCTTACAAGACAGCCAAATTGAAATCATTCTTCCCCAGGAAGACGCGTAAATTTTTGTGACCGTTAACCTGGCTTCATACCTTGCCACATAGc                                                                                           <  2:965413/139‑1 (MQ=255)                                           cTGCGCTCTTTACCGACGAGCTGCCGCCCGCTGCGCTAAAATCACGCTTACAAGACAGCCAAATTGAAATCATTCTTCCCCAGGAAGACGCGTAAATTTTtgt                                                                                                                              <  1:195189/103‑1 (MQ=255)                                           cTGCGCTCTTTACCGACGAGCTGCCGCCCGCTGCGCTAAAATCACGCTTACAAGACAGCCAAATTGAAATCATTCTTCCCCAGGAAGACGCGTAAATTTTtgt                                                                                                                              >  2:195189/1‑103 (MQ=255)                                              cgcTCTTTACCGACGAGCTGCCGCCCGCTGCGCTAAAATCACGCTTACAAGACAGCCAAATTGAAATCATTCTTCCCCAGGAAGACGCGTAAATTTTTGTGACCGTTAACCTGGCTTCATACCTTGCCACATAGc                                                                                           <  2:768597/135‑1 (MQ=255)                                              cgcTCTTTACCGACGAGCTGCCGCCCGCTGCGCTAAAATCACGCTTACAAGACAGCCAAATTGAAATCATTCTTCCCCAGGAAGACGCGTAAATTTTTGTGACCGTTAACCTGGCTTCATACCTTGCCACATAGc                                                                                           >  1:768597/1‑135 (MQ=255)                                                     tACCGACGAGCTGCCGCCCGCTGCGCTAAAATCACGCTTACAAGACAGCCAAATTGAAATCATTCTTCCCCAGGAAGACGCGTAAATTTTTGTGACCGTTAACCTGGCTTCATACCTTGCCACATAGCCAAACCCATcc                                                                                <  1:218799/139‑1 (MQ=255)                                                     tACCGACGAGCTGCCGCCCGCTGCGCTAAAATCACGCTTACAAGACAGCCAAATTGAAATCATTCTTCCCCAGGAAGACGCGTAAATTTTTGTGACCGTTAACCTGGCTTCATACCTTGCCACATAGCCAAACCCATcc                                                                                >  2:218799/1‑139 (MQ=255)                                                        cgacgaGCTGCCGCCCGCTGCGCTAAAATCACGCTTACAAGACAGCCAAATTGAAATCATTCTTCCCCAGGAAGACGCGTAAATTTTTGTGACCGTTAACCTGGCTTCATACCTTGCCACATAGCCAAACCCATCCttt                                                                             >  2:670364/1‑139 (MQ=255)                                                                  ccgcccgcTGCGCTAAAATCACGCTTACAAGACAGCCAAATTGAAATCATTCTTCCCCAGGAAGACGCGTAAATTTTTGTGACCGTTAACCTGGCTTCATACCTTGCCACATAGCCAAACCCATCCTTTCCCGCTACAg                                                                   <  2:158701/139‑1 (MQ=255)                                                                                            aCAAGACAGCCAAATTGAAATCATTCTTCCCCAGGAAGACGCGTAAATTTTTGTGACCGTTAACCTGGCTTCATACCTTGCCACATAGCCAAACCCATCCTTTCCCGCTACAGTTAATTTCTTGTGGCGCGAAaggagg                                         <  1:1194056/139‑1 (MQ=255)                                                                                                  cAGCCAAATTGAAATCATTCTTCCCCAGGAAGACGCGTAAATTTTTGTGACCGTTAACCTGGCTTCATACCTTGCCACATAGCCAAACCCATCCTTTCCCGCTACAGTTAATTTCTTGTGGCGCGAAAGGAGGCaaaaa                                   <  1:386946/139‑1 (MQ=255)                                                                                                  cAGCCAAATTGAAATCATTCTTCCCCAGGAAGACGCGTAAATTTTTGTGACCGTTAACCTGGCTTCATACCTTGCCACATAGCCAAACCCATCCTTTCCCGCTACAGTTAATTTCTTGTGGCGCGAAAGGAGGCaaaaa                                   >  2:760182/1‑139 (MQ=255)                                                                                                  cAGCCAAATTGAAATCATTCTTCCCCAGGAAGACGCGTAAATTTTTGTGACCGTTAACCTGGCTTCATACCTTGCCACATAGCCAAACCCATCCTTTCCCGCTACAGTTAATTTCTTGTGGCGCGAAAGGAGGCaaaaa                                   >  2:1069919/1‑139 (MQ=255)                                                                                                    gCCAAATTGAAATCATTCTTCCCCAGGAAGACGCGTAAATTTTtgt                                                                                                                              >  2:1140406/1‑46 (MQ=255)                                                                                                    gCCAAATTGAAATCATTCTTCCCCAGGAAGACGCGTAAATTTTtgt                                                                                                                              <  1:1140406/46‑1 (MQ=255)                                                                                                     ccAAATTGAAATCATTCTTCCCCAGGAAGACGCGTAAATTTTtgt                                                                                                                              >  1:944487/1‑45 (MQ=255)                                                                                                     ccAAATTGAAATCATTCTTCCCCAGGAAGACGCGTAAATTTTtgt                                                                                                                              <  2:944487/45‑1 (MQ=255)                                                                                                     ccAAATTGAAATCATTCTTCCCCAGGAAGACGCGTAAATTTTTGTGACCGTTAACCTGGCTTCATACCTTGCcac                                                                                                >  1:968649/1‑75 (MQ=255)                                                                                                     ccAAATTGAAATCATTCTTCCCCAGGAAGACGCGTAAATTTTTGTGACCGTTAACCTGGCTTCATACCTTGCcac                                                                                                <  2:968649/75‑1 (MQ=255)                                                                                                     ccAAATTGAAATCATTCTTCCCCAGGAAGACGCGTAAATTTTTGTGACCGTTAACCTGGCTTCATACCTTGCCACATAGCCAAACCCATCCTTTCCCGCTACAGTTAATTTCTTGTGGCGCGAAAGGAGGCAAAAATGc                                >  2:339610/1‑139 (MQ=255)                                                                                                      cAAATTGAAATCATTCTTCCCCAGGAAGACGCGTAAATTTTTGTGACCGTTAACCTGGCTTCATACCTTGCCACATAGCCAAACCCATCCTTTCCCGCTACAGTTAATTTCTTGTGGCGCGAAAGGAGGCAAAAATGct                               <  1:760182/139‑1 (MQ=255)                                                                                                         aTTGAAATCATTCTTCCCCAGGAAGACGCGTAAATTTTTGTGACCGTTAACCTGGCTTCATACCTTGCCACATAGCCAAACCCATCCTTTCCCGCTACAGTTAATTTCTTGTGGCGCGAAAGGAGGCAAAAATGCTCta                            <  2:505217/139‑1 (MQ=255)                                                                                                              aaTCATTCTTCCCCAGGAAGACGCGTAAATTTTTGTGACCGTTAACCTGGCTTCATACCTTGCCACATAGCCAAACCCATCCTTTCCCGCTACAGTTAATTTCTTGTGGCGCGAAAGGAGGCAAAAATGCTCTATATCt                       <  2:639466/139‑1 (MQ=255)                                                                                                              aaTCATTCTTCCCCAGGAAGACGCGTAAATTTTTGTGACCGTTAACCTGGCTTCATACCTTGCCACATAGCCAAACCCATCCTTTCCCGCTACAGTTAATTTCTTGTGGCGCGAAAGGAGGCAAAAATGCTCTATATCt                       <  1:1069919/139‑1 (MQ=255)                                                                                                               aTCATTCTTCCCCAGGAAGACGCGTAAATTTTTGTGACCGTTAACCTGGCTTCATACCTTGCCACATAGCCAAACCCATCCTTTCCCGCTACAGTTAATTTCTTGTGGCGCGAAaggagg                                         >  2:1407557/1‑120 (MQ=255)                                                                                                               aTCATTCTTCCCCAGGAAGACGCGTAAATTTTTGTGACCGTTAACCTGGCTTCATACCTTGCCACATAGCCAAACCCATCCTTTCCCGCTACAGTTAATTTCTTGTGGCGCGAAaggagg                                         <  1:1407557/120‑1 (MQ=255)                                                                                                                    tcttcCCCAGGAAGACGCGTAAATTTTTGTGACCGTTAACCTGGCTTCATACCTTGCCACATAGCCAAACCCATCCTTTCCCGCTACAGTTAATTTCTTGTGGCGCGAAAGGAGGCAAAAATGCTCTATATCTTTGAtt                 >  2:176618/1‑139 (MQ=255)                                                                                                                        ccccAGGAAGACGCGTAAATTTTTGTGACCGTTAACCTGGCTTCATACCTTGCCACATAGCCAAACCCATCCTTTCCCGCTACAGTTAATTTCTTGTGGCGCGAAAGGAGGCAAAAATGCTCTATATCTTTGATTTAgg             <  1:176618/139‑1 (MQ=255)                                                                                                                         cccAGGAAGACGCGTAAATTTTTGTGACCGTTAACCTGGCTTCATACCTTGCCACATAGCCAAACCCa                                                                                   <  1:859123/68‑1 (MQ=255)                                                                                                                         cccAGGAAGACGCGTAAATTTTTGTGACCGTTAACCTGGCTTCATACCTTGCCACATAGCCAAACCCa                                                                                   >  2:859123/1‑68 (MQ=255)                                                                                                                           cAGGAAGACGCGTAAATTTTTGTGACCGTTAACCTGGCTTCATACCTTGCCACATAGCCAAACCCATCCTTTCCCGCTACAGTTAATTTCTTGTGGCGCGAAAGGAGGCAAAAATGCTCTATATCt                       <  2:1073081/126‑1 (MQ=255)                                                                                                                           cAGGAAGACGCGTAAATTTTTGTGACCGTTAACCTGGCTTCATACCTTGCCACATAGCCAAACCCATCCTTTCCCGCTACAGTTAATTTCTTGTGGCGCGAAAGGAGGCAAAAATGCTCTATATCt                       >  1:1073081/1‑126 (MQ=255)                                                                                                                              gAAGACGCGTAAATTTTTGTGACCGTTAACCTGGCTTCATACCTTGCCACATAGCCAAACCCATCCTTTCCCGCTACAGTTAATTTCTTGTGGCGCGAAaggagg                                         <  2:316672/105‑1 (MQ=255)                                                                                                                              gAAGACGCGTAAATTTTTGTGACCGTTAACCTGGCTTCATACCTTGCCACATAGCCAAACCCATCCTTTCCCGCTACAGTTAATTTCTTGTGGCGCGAAaggagg                                         >  1:316672/1‑105 (MQ=255)                                                                                                                              gAAGACGCGTAAATTTTTGTGACCGTTAACCTGGCTTCATACCTTGCCACATAGCCAAACCCATCCTTTCCCGCTACAGTTAATTTCTTGTGGCGCGAAAGGAGGCAAAAATg                                 <  2:431467/113‑1 (MQ=255)                                                                                                                              gAAGACGCGTAAATTTTTGTGACCGTTAACCTGGCTTCATACCTTGCCACATAGCCAAACCCATCCTTTCCCGCTACAGTTAATTTCTTGTGGCGCGAAAGGAGGCAAAAATg                                 >  1:431467/1‑113 (MQ=255)                                                                                                                                  aCGCGTAAATTTTTGTGACCGTTAACCTGGCTTCATACCTTGCCACATAGCCAAACCCATCCTTTCCCGCTACAGTTAATTTCTTGTGGCGCGAAaggagg                                         <  2:180295/101‑1 (MQ=255)                                                                                                                                  aCGCGTAAATTTTTGTGACCGTTAACCTGGCTTCATACCTTGCCACATAGCCAAACCCATCCTTTCCCGCTACAGTTAATTTCTTGTGGCGCGAAaggagg                                         >  1:180295/1‑101 (MQ=255)                                                                                                                                   cgcgTAAATTTTTGTGACCGTTAACCTGGCTTCATACCtt                                                                                                     <  1:438054/40‑1 (MQ=255)                                                                                                                                   cgcgTAAATTTTTGTGACCGTTAACCTGGCTTCATACCtt                                                                                                     >  2:438054/1‑40 (MQ=255)                                                                                                                                   cgcgTAAATTTTTGTGACCGTTAACCTGGCTTCATACCTTGCCACATAGCCAAACCCATCCTTTCCCGCt                                                                       <  1:482892/70‑1 (MQ=255)                                                                                                                                   cgcgTAAATTTTTGTGACCGTTAACCTGGCTTCATACCTTGCCACATAGCCAAACCCATCCTTTCCCGCt                                                                       >  2:482892/1‑70 (MQ=255)                                                                                                                                   cgcgTAAATTTTTGTGACCGTTAACCTGGCTTCATACCTTGCCACATAGCCAAACCCATCCTTTCCCGCTACAGt                                                                  <  1:1366672/75‑1 (MQ=255)                                                                                                                                   cgcgTAAATTTTTGTGACCGTTAACCTGGCTTCATACCTTGCCACATAGCCAAACCCATCCTTTCCCGCTACAGt                                                                  >  2:1366672/1‑75 (MQ=255)                                                                                                                                   cgcgTAAATTTTTGTGACCGTTAACCTGGCTTCATACCTTGCCACATAGCCAAACCCATCCTTTCCCGCTACAGTTAATTTCTTGTGGCGCGAAAGGAGGCAAAAATGCTCTATATCTTTGATTTAGGTAATGTGATTg  >  2:1164895/1‑139 (MQ=255)                                                                                                                                   cgcgTAAATCTTTGTGACCGTTAACCTGGCTTCATACCTTGCCACATAGCCAAACCCATCCTTTCCCGCTACAGTTAATTTCTTGTGGCGCGAAAGGAGGCAAAAATGCTCTATATCTTTGATTTAGGTAATGTGATTg  >  2:113341/1‑139 (MQ=255)                                                                                                                                     cgTAAATTTTTGTGACCGTTAACCTGGCTTCATACCTTGCCACATAGCCAAAcc                                                                                     <  2:626386/54‑1 (MQ=255)                                                                                                                                     cgTAAATTTTTGTGACCGTTAACCTGGCTTCATACCTTGCCACATAGCCAAAcc                                                                                     >  1:626386/1‑54 (MQ=255)                                                                                                                                         |                                                                                                                                     ATGCTTCGGCGGCGGTTGAAATTGGTAACGTGGCACAGGTCACTGCGCTCTTTACCGACGAGCTGCCGCCCGCTGCGCTAAAATCACGCTTACAAGACAGCCAAATTGAAATCATTCTTCCCCAGGAAGACGCGTAGATTTTTGTGACCGTTAACCTGGCTTCATACCTTGCCACATAGCCAAACCCATCCTTTCCCGCTACAGTTAATTTCTTGTGGCGCGAAAGGAGGCAAAAATGCTCTATATCTTTGATTTAGGTAATGTGATTG  >  NC_000913/4075318‑4075586

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