Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,859,148	A→G	F361S (TTC→TCC)	ydjK ←	putative MFS sugar transporter, membrane protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,859,148	0	A	G	100.0%	166.2 / NA	44	F361S (TTC→TCC)	ydjK	putative MFS sugar transporter, membrane protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (21/23);  total (21/23)

CAACGGCATAAGGTGCGGCAATACCACTGATTCGCCCTACCGCATTCGCCAGACCGGAGCCACGGAGTTTGGCCTCTGTCGGCCAGATTTCAGGGACATACACTGCCGAGGCATAGCAAACGTACATATAGACGAAAGTAATCAGGAAGAAACCAATCAACGTTATGAGCAACATACTGGTTTGCAGTGAATAGATATATCCGAGCACCGCAATCAGGATTAATAGCCCCACACCCATTGTTTTACGCGGAATTTTATCCATCACCAGCATGGCAATAAAAATACCAAATGGCG  >  NC_000913/1859002‑1859295                                                                                                                                                   |                                                                                                                                                    cAACGGCATAAGGTGCGGCAATACCACTGATTCGCCCTACCGCATTCGCCAGACCGGAGCCACGGAGTTTGGCCTCTGTCGGCCAGATTTCAGGGACATACACTGCCGAGGCATAGCAAACGTACATATAGACGAAAGTAATCAGggagaa                                                                                                                                                 <  2:310247/151‑1 (MQ=255)     ggCATAAGGTGCGGCAATACCACTGATTCGCCCTACCGCATTCGCCAGACCGGAGCCACGGAGTTTGGCCTCTGTCGGCCAGATTTCAGGGACATACACTGCCGAGGCATAGCAAACGTACATATAGACGAAAGTAATCAGGGAGAAAcc                                                                                                                                              >  2:193869/1‑150 (MQ=255)         tAAGGTGCGGCAATACCACTGATTCGCCCTACCGCATTCGCCAGACCGGAGCCACGGAGTTTGGCCTCTGTCGGCCAGATTTCAGGGACATACACTGCCGAGGCATAGCAAACGTACATATAGACGAAAGTAATCAGGGAGAAACCAAt                                                                                                                                           <  1:539977/149‑1 (MQ=255)          aaGGTGCGGCAATACCACTGATTCGCCCTACCGCATTCGCCAGACCGGAGCCACGGAGTTTGGCCTCTGTCGGCCAGATTTCAGGGACATACACTGCCGAGGCATAGCAAACGTACATATAGACGAAAGTAATCAGGGAGAAACCAATCa                                                                                                                                         >  1:645490/1‑150 (MQ=255)            ggTGCGGCAATACCACTGATTCGCCCTACCGCATTCGCCAGACCGGAGCCACGGAGTTTGGCCTCTGTCGGCCAGATTTCAGGGACATACACTGCCGAGGCATAGCAAACGTACATATAGACGAAAGTAATCAGGGAGAAACCAATCAACg                                                                                                                                      <  1:203258/151‑1 (MQ=255)             gTGCGGCAATACCACTGATTCGCCCTACCGCATTCGCCAGACCGGAGCCACGGAGTTTGGCCTCTGTCGGCCAGATTTCAGGGACATACACTGCCGAGGCATAGCAAACGTACATATAGACGAAAGTAATCAGGGAGAAACCAATCAACGt                                                                                                                                     <  2:503062/151‑1 (MQ=255)               gCGGCAATACCACTGATTCGCCCTACCGCATTCGCCAGACCGGAGCCACGGAGTTTGGCCTCTGTCGGCCAGATTTCAGGGACATACACTGCCGAGGCATAGCAAACGTACATATAGACGAAAGTAATCAGGGAGAAACCAATCAACGTTa                                                                                                                                   <  2:470895/151‑1 (MQ=255)                     aTACCACTGATTCGCCCTACCGCATTCGCCAGACCGGAGCCACGGAGTTTGGCCTCTGTCGGCCAGATTTCAGGGACATACACTGCCGAGGCATAGCAAACGTACATATAGACGAAAGTAATCAGGGAGAAACCAATCAACGTTATGAGCa                                                                                                                             <  1:354171/151‑1 (MQ=255)                        ccACTGATTCGCCCTACCGCATTCGCCAGACCGGAGCCACGGAGTTTGGCCTCTGTCGGCCAGATTTCAGGGACATACACTGCCGAGGCATAGCAAACGTACATATAGACGAAAGTAATCAGGGAGAAACCAATCAACGTTATGAGCAACa                                                                                                                          <  2:85514/151‑1 (MQ=255)                                       aCCGCATTCGCCAGACCGGAGCCACGGAGTTTGGCCTCTGTCGGCCAGATTTCAGGGACATACACTGCCGAGGCATAGCAAACGTACATATAGACGAAAGTAATCAGGGAGAAACCAATCAACGTTATGAGCAACATACTGGTTTGCAg                                                                                                             <  1:123163/149‑1 (MQ=255)                                       aCCGCATTCGCCAGACCGGAGCCACGGAGTTTGGCCTCTGTCGGCCAGATTTCAGGGACATACACTGCCGAGGCATAGCAAACGTACATATAGACGAAAGTAATCAGGGAGAAACCAATCAACGTTATGAGCAACATACTGGTTTGCAg                                                                                                             <  1:18132/149‑1 (MQ=255)                                       aCCGCATTCGCCAGACCGGAGCCACGGAGTTTGGCCTCTGTCGGCCAGATTTCAGGGACATACACTGCCGAGGCATAGCAAACGTACATATAGACGAAAGTAATCAGGGAGAAACCAATCAACGTTATGAGCAACATACTGGTTTGCAGTg                                                                                                           <  2:769884/151‑1 (MQ=255)                                       aCCGCATTCGCCAGACCGGAGCCACGGAGTTTGGCCTCTGTCGGCCAGATTTCAGGGACATACACTGCCGAGGCATAGCAAACGTACATATAGACGAAAGTAATCAGGGAGAAACCAATCAACGTTATGAGCAACATACTGGTTTGCAGTg                                                                                                           >  1:218208/1‑151 (MQ=255)                                          gCATTCGCCAGACCGGAGCCACGGAGTTTGGCCTCTGTCGGCCAGATTTCAGGGACATACACTGCCGAGGCATAGCAAACGTACATATAGACGAAAGTAATCAGGGAGAAACCAATCAACGTTATGAGCAACATACTGGTTTGCAGTGAAt                                                                                                        <  2:753712/151‑1 (MQ=255)                                          gCATTCGCCAGACCGGAGCCACGGAGTTTGGCCTCTGTCGGCCAGATTTCAGGGACATACACTGCCGAGGCATAGCAAACGTACATATAGACGAAAGTAATCAGGGAGAAACCAATCAACGTTATGAGCAACATACTGGTTTGCAGTGAAt                                                                                                        <  2:420227/151‑1 (MQ=255)                                           cATTCGCCAGACCGGAGCCACGGAGTTTGGCCTCTGTCGGCCAGATTTCAGGGACATACACTGCCGAGGCATAGCAAACGTACATATAGACGAAAGTAATCAGGGAGAAACCAATCAACGTTATGAGCAACATACTGGTTTGCAGTGAATa                                                                                                       <  1:719361/151‑1 (MQ=255)                                            aTTCGCCAGACCGGAGCCACGGAGTTTGGCCTCTGTCGGCCAGATTTCAGGGACATACACTGCCGAGGCATAGCAAACGTACATATAGACGAAAGTAATCAGGGAGAAACCAATCAACGTTATGAGCAACATACTGGTTTGCAGTGAAt                                                                                                        >  2:448071/1‑149 (MQ=255)                                               cGCCAGACCGGAGCCACGGAGTTTGGCCTCTGTCGGCCAGATTTCAGGGACATACACTGCCGAGGCATAGCAAACGTACATATAGACGAAAGTAATCAGGGAGAAACCATTCAACGTTATGAGCAACATACTGGTTTGCAGTGAATAGat                                                                                                    >  2:222411/1‑150 (MQ=255)                                               cGCCAGACCGGAGCCACGGAGTTTGGCCTCTGTCGGCCAGATTTCAGGGACATACACTGCCGAGGCATAGCAAACGTACATATAGACGAAAGTAATCAGGGAGAAACCAATCAACGTTATGAGCAACATACTGGTTTGCAGTGAATAGat                                                                                                    >  1:701281/1‑150 (MQ=255)                                                 ccAGACCGGAGCCACGGAGTTTGGCCTCTGTCGGCCAGATTTCAGGGACATACACTGCCGAGGCATAGCAAACGTACATATAGACGAAAGTAATCAGGGAGAAACCAATCAACGTTATGAGCAACATACTGGTTTGCAGTGAATAGatata                                                                                                 <  2:243627/151‑1 (MQ=255)                                                  cAGACCGGAGCCACGGAGTTTGGCCTCTGTCGGCCAGATTTCAGGGACATACACTGCCGAGGCATAGCAAACGTACATATAGACGAAAGTAAGCAGGGAGAAACCAATCAACGTTATGAGCAACATACTGGTTTGCAGTGAATAGatatat                                                                                                <  1:118563/151‑1 (MQ=255)                                                         gAGCCACGGAGTTTGGCCTCTGTCGGCCAGATTTCAGGGACATACACTGCCGAGGCATAGCAAACGTACATATAGACGAAAGTAATCAGGGAGAAACCAATCAACGTTATGAGCAACATACTGGTTTGCAGTGAATAGATATATCCGAGc                                                                                          >  2:721640/1‑150 (MQ=255)                                                                ggAGTTTGGCCTCTGTCGGCCAGATTTCAGGGACATACACTGCCGAGGCATAGCAAACGTACATATAGACGAAAGTAATCAGGGAGAAACCAATCAACGTTATGAGCAACATACTGGTTTGCAGTGAATAGATATATCCGAGCACCGCAAt                                                                                  >  1:691211/1‑151 (MQ=255)                                                                ggAGTTTGGCCTCTGTCGGCCAGATTTCAGGGACATACACTGCCGAGGCATAGCAAACGTACATATAGACGAAAGTAATCAGGGAGAAACCAATCAACGTTATGAGCAACATACTGGTTTGCAGTGAATAGATATATCCGAGCACCGCAAt                                                                                  <  2:691285/151‑1 (MQ=255)                                                                      tGGCCTCTGTCGGCCAGATTTCAGGGACATACACTGCCGAGGCATAGCAAACGTACATATAGACGAAAGTAATCAGGGAGAAACCAATCAACGTTATGAGCAACATACTGGTTTGCAGTGAATAGATATATCCGAGCACCGCAATCAGGAt                                                                            >  1:367932/1‑151 (MQ=255)                                                                        gCCTCTGTCGGCCAGATTTCAGGGACATACACTGCCGAGGCATAGCAAACGTACATATAGACGAAAGTAATCAGGGAGAAACCAATCAACGTTATGAGCAACATACTGGTTTGCAGTGAATAGATATATCCGAGCACCGCAATCAGGAtt                                                                           >  2:304390/1‑150 (MQ=255)                                                                                   ccAGATTTCAGGGACATACACTGCCGAGGCATAGCAAACGTACATATAGACGAAAGTAATCAGGGAGAAACCAATCAACGTTATGAGCAACATACTGGTTTGCAGTGAATAGATATATCCGAGCACCGCAATCAGGATTAATAGCCCcac                                                                >  1:446985/1‑150 (MQ=255)                                                                                    cAGATTTCAGGGACATACACTGCCGAGGCATAGCAAACGTACATATAGACGAAAGTAATCAGGGAGAAACCAATCAACGTTATGAGCAACATACTGGTTTGCAGTGAATAGATATATCCGAGCACCGCAATCAGGATTAATAGCCCCacac                                                              >  1:275522/1‑151 (MQ=255)                                                                                              ggACATACACTGCCGAGGCATAGCAAACGTACATATAGACGAAAGTAATCAGGGAGAAACCAATCAACGTTATGAGCAACATACTGGTTTGCAGTGAATAGATATATCCGAGCACCGCAATCAGGATTAATAGCCCCACACCCATTGtttt                                                    <  1:137527/151‑1 (MQ=255)                                                                                                  aTACACTGCCGAGGCATAGCAAACGTACATATAGACGAAAGTAATCAGGGAGAAACCAATCAACGTTATGAGCAACATACTGGTTTGCAGTGAATAGATATATCCGAGCACCGCAATCAGGATTAATAGCCCCACACCCATTGTTTTAcgc                                                >  2:81404/1‑151 (MQ=255)                                                                                                           cGAGGCATAGCAAACGTACATATAGACGAAAGTAATCAGGGAGAAACCAATCAACGTTATGAGCAACATACTGGTTTGCAGTGAATAGATATATCCGAGCACCGCAATCAGGATTAATAGCCCCACACCCATTGTTTTACGCGGAAtttt                                        >  1:569553/1‑150 (MQ=255)                                                                                                           cGAGGCATAGCAAACGTACATATAGACGAAAGTAATCAGGGAGAAACCAATCAACGTTATGAGCAACATACTGGTTTGCAGTGAATAGATATATCCGAGCACCGCAATCAGGATTAATAGCCCCACACCCATTGTTTTACGCGGAATTTTa                                       <  2:486683/151‑1 (MQ=255)                                                                                                            gAGGCATAGCAAACGTACATATAGACGAAAGTAATCAGGGAGAAACCAATCAACGTTATGAGCAACATACTGGTTTGCAGTGAATAGATATATCCGAGCACCGCAATCAGGATTAATAGCCCCACACCCATTGTTTTACGCGGAATTTTAt                                      >  2:593710/1‑151 (MQ=255)                                                                                                            gAGGCATAGCAAACGTACATATAGACGAAAGTAATCAGGGAGAAACCAATCAACGTTATGAGCAACATACTGGTTTGCAGTGAATAGATATATCCGAGCACCGCAATCAGGATTAATAGCCCCACACCCATTGTTTTACGCGGAATTTTAt                                      <  2:599134/151‑1 (MQ=255)                                                                                                              ggCATAGCAAACGTACATATAGACGAAAGTAATCAGGGAGAAACCAATCAACGTTATGAGCAACATACTGGTTTGCAGTGAATAGATATATCCGAGCACCGCAATCAGGATTAATAGCCCCACACCCATTGTTTTACGCGGAATTTTATcc                                    <  2:456488/151‑1 (MQ=255)                                                                                                                 aTAGCAAACGTACATATAGACGAAAGTAATCAGGGAGAAACCAATCAACGTTATGAGCAACATACTGGTTTGCAGTGAATAGATATATCCGAGCACCGCAATCAGGATTAATAGCCCCACACCCATTGTTTTACGTGGAATTTTATCCATc                                 <  2:107167/151‑1 (MQ=255)                                                                                                                                ataGACGAAAGTAATCAGGGAGAAACCAATCAACGTTATGAGCAACATACTGGTTTGCAGTGAATAGATATATCCGAGCACCGCAATCAGGATTAATAGCCCCACACCCATTGTTTTACGCGGAATTTTATCCATCACCAGCATGGCAAt                   <  2:386240/150‑1 (MQ=255)                                                                                                                                ataGACGAAAGTAATCAGGGAGAAACCAATCAACGTTATGAGCAACATACTGGTTTGCAGTGAATAGATATATCCGAGCACCGCAATCAGGATTAATAGCCCCACACCCATTGTTTTACGCGGAATTTTATCCATCACCAGCATGGCAAt                   >  1:190173/1‑150 (MQ=255)                                                                                                                                ataGACGAAAGTAATCAGGGAGAAACCAATCAACGTTATGAGCAACATACTGGTTTGCAGTGAATAGATATATCCGAGCACCGCAATCAGGATTAATAGCCCCACACCCATTGTTTTACGCGGAATTTTATCCATCACCAGCATGGCAATa                  <  2:243580/151‑1 (MQ=255)                                                                                                                                  aGACGAAAGTAATCAGGGAGAAACCAATCAACGTTATGAGCAACATACTGGTTTGCAGTGAATAGATATATCAGAGCACCGCAATCAGGATTAATAGCCCCACACCCATTGTTTTACGCGGAATTTTATCCATCACCAGCATGGCAATaa                 >  1:192162/1‑150 (MQ=255)                                                                                                                                       aaaGTAATCAGGGAGAAACCAATCAACGTTATGAGCAACATACTGGTTTGCAGTGAATAGATATATCCGAGCACCGCAATCAGGATTAATAGCCCCACACCCATTGTTTTACGCGGAATTTTATCCATCACCAGCATGGCAATAAAAATAc           >  2:162590/1‑151 (MQ=255)                                                                                                                                         aGTAATCAGGGAGAAACCAATCAACGTTATGAGCAACATACTGGTTTGCAGTGAATAGATATATCCGAGCACCGCAATCAGGATTAATAGCCCCACACCCATTGTTTTACGCGGAATTTTATCCATCACCAGCATGGCAATAAAAATAcc          >  1:532320/1‑150 (MQ=255)                                                                                                                                         aGTAATCAGGGAGAAACCAATCAACGTTATGAGCAACATACTGGTTTGCAGTGAATAGATATATCCGAGCACCGCAATCAGGATTAATAGCCCCACACCCATTGTTTTACGCGGAATTTTATCCATCACCAGCATGGCAATAAAAATAcc          <  2:532372/150‑1 (MQ=255)                                                                                                                                                aGGGAGAAACCAATCAACGTTATGAGCAACATACTGGTTTGCAGTGAATAGATATATCCGAGCACCGCAATCAGGATTAATAGCCCCACACCCATTGTTTTACGCGGAATTTTATCCATCACCAGCATGGCAATAAAAATACCAAATGgcg  >  2:732613/1‑151 (MQ=255)                                                                                                                                                aGGGAGAAACCAATCAACGTTATGAGCAACATACTGGTTTGCAGTGAATAGATATATCCGAGCACCGCAATCAGGATTAATAGCCCCACACCCATTGTTTTACGCGGAATTTTATCCATCACCAGCATGGCAATAAAAATACCAAATGgc   >  1:601215/1‑150 (MQ=255)                                                                                                                                                   |                                                                                                                                                    CAACGGCATAAGGTGCGGCAATACCACTGATTCGCCCTACCGCATTCGCCAGACCGGAGCCACGGAGTTTGGCCTCTGTCGGCCAGATTTCAGGGACATACACTGCCGAGGCATAGCAAACGTACATATAGACGAAAGTAATCAGGAAGAAACCAATCAACGTTATGAGCAACATACTGGTTTGCAGTGAATAGATATATCCGAGCACCGCAATCAGGATTAATAGCCCCACACCCATTGTTTTACGCGGAATTTTATCCATCACCAGCATGGCAATAAAAATACCAAATGGCG  >  NC_000913/1859002‑1859295

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