Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	73,641	A→G	L294S (TTA→TCA)	thiP ←	thiamine ABC transporter membrane subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	73,641	0	A	G	100.0%	156.8 / NA	45	L294S (TTA→TCA)	thiP	thiamine ABC transporter membrane subunit
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (18/27);  total (18/27)

GGCAACTGGCGATTTACCCCATCGACGATCACCGCCAGTAACGGTGGCAGCAACAGCAGCAGCGCCAGCACAATTAACACCGTGTCGCAAATGCGGCTATGCAGACGATCGTCCGGGTCGCGCCAGCCTTGCAGCAGCGTGGTGCCGGGCGCA  >  NC_000913/73566‑73718                                                                            |                                                                              ggCAACTGGCGATTTACCCCATCGACGATCACCGCCAGTAACGGTGGCAGCAACAGCAGCAGCGCCAGCACAATTGACAcc                                                                          >  2:16360/1‑81 (MQ=255)  gCAACTGGCGATTTACCCCATCGACGATCACCGCCAGTAACGGTGGCAGCAACAGCAGCAGCGCCAGCACAATTGACACCg                                                                         >  1:16356/1‑81 (MQ=255)          cGATTTACCCCATCGACGATCACCGCCAGTAACGGTGGCAGCAACAGCAGCAGCGCCAGCACAATTGACACCGTGTCGCaa                                                                 <  1:16308/81‑1 (MQ=255)           gATTTACCCCATCGACGATCACCGCCAGTAACGGTGGCAGCAACAGCAGCAGCGCCAGCACAATTGACACCGTGTCGCaaa                                                                <  2:16309/81‑1 (MQ=255)            aTTTACCCCATCGACGATCACCGCCAGTAACGGTGGCAGCAACAGCAGCAGCGCCAGCACAATTGACACCGTGTCGCAAAt                                                               <  1:16310/81‑1 (MQ=255)            aTTTACCCCATCGACGATCACCGCCAGTAACGGTGGCAGCAACAGCAGCAGCGCCAGCACAATTGACACCGTGTCGCAAAt                                                               <  1:16311/81‑1 (MQ=255)            aTTTACCCCATCGACGATCACCGCCAGTAACGGTGGCAGCAACAGCAGCAGCGCCAGCACAATTGACACCGTGTCGCAAAt                                                               <  1:16312/81‑1 (MQ=255)             tttACCCCATCGACGATCACCGCCAGTAACGGTGGCAGCAACAGCAGCAGCGCCAGCACAATTGACACCGTGTCGCAAATg                                                              <  1:16313/81‑1 (MQ=255)                 ccccATCGAGGATCACCGCCAGTAAAGGTGGCAGAAACAGCAGCAGCGCCAGCACAATTGACACCGTGTCGCAAATGCGGc                                                          <  2:16314/81‑1 (MQ=255)                    cATCGACGATCACCGCCAGTAACGGTGGCAGCAACAGCAGCAGCGCCAGCACAATTGACACCGTGTCGCAAATGCGGCTAt                                                       >  1:16371/1‑81 (MQ=255)                          cgaTCACCGCCAGTAACGGTGGCAGCAACAGCAGCAGCGCCAGCACAATTGACACCGTGTCGCAAATGCGGCTATGCAGAc                                                 <  1:16315/81‑1 (MQ=255)                            aTCACCGCCAGTAACGGTGGCAGCAACAGCAGCAGCGCCAGCACAATTGACACCGTGTCGCAAATGCGGCTATGCAGACGa                                               >  2:16363/1‑81 (MQ=255)                             tCACCGCCAGTAACGGTGGCAGCAACAGCAGCAGCGCCAGCACAATTGACACCGTGTCGCAAATGCGGCTATGCAGACGAt                                              <  1:16317/81‑1 (MQ=255)                             tCACCGCCAGTAACGGTGGCAGCAACAGCAGCAGCGCCAGCACAATTGACACCGTGTCGCAAATGCGGCTATGCAGACGAt                                              <  1:16316/81‑1 (MQ=255)                                  gCCAGTAACGGTGGCAGCAACAGCAGCAGCGCCAGCACAATTGACACCGTGTCGCAAATGCGGCTATGCAGACGATCGTcc                                         <  1:16319/81‑1 (MQ=255)                                  gCCAGTAACGGTGGCAGCAACAGCAGCAGCGCCAGCACAATTGACACCGTGTCGCAAATGCGGCTATGCAGACGATCGTcc                                         <  1:16318/81‑1 (MQ=255)                                     aGTAACGGTGGCAGCAACAGCAGCAGCGCCAGCACAATTGACACCGTGTCGCAAATGCGGCTATGCAGACGATCGTCCggg                                      >  2:16359/1‑81 (MQ=255)                                       tAACGGTGGCAGCAACAGCAGCAGCGCCAGCACAATTGACACCGTGTCGCAAATGCGGCTATGCAGACGATCGTCCGGGTc                                    <  2:16320/81‑1 (MQ=255)                                       tAACGGTGGCAGCAACAGCAGCAGCGCCAGCACAATTGACACCGTGTCGCAAATGCGGCTATGCAGACGATCGTCCGGGTc                                    <  2:736580/81‑1 (MQ=255)                                          cGGTGGCAGCAACAGCAGCAGCGCCAGCACAATTGACACCGTGTCGCAAATGCGGCTATGCAGACGATCGTCCGGGTcgcg                                 <  2:16321/81‑1 (MQ=255)                                           ggTGGCAGCAACAGCAGCAGCGCCAGCACAATTGACACCGTGTCGCAAATGCGGCTATGCAGACGATCGTCCGGGTCgcgc                                >  2:16362/1‑81 (MQ=255)                                             tGGCAGCAACAGCAGCAGCGCCAGCACAATTGACACCGTGTCGCAAATGCGGCTATGCAGACGATCGTCCGGGTCgcgc                                >  2:16358/1‑79 (MQ=255)                                             tGGCAGCAACAGCAGCAGCGCCAGCACAATTGACACCGTGTCGCAAATGCGGCTATGCAGACGATCGTCCGGGTCgcgc                                >  2:16357/1‑79 (MQ=255)                                              ggCAGCAACAGCAGCAGCGCCAGCACAATTGACACCGTGTCGCAAATGCGGCTATGCAGACGATCGTCCGGGTCGCGCCAg                             >  2:16369/1‑81 (MQ=255)                                                cagcaACAGCAGCAGCGCCAGCACAATTGACACCGTGTCGCAAATGCGGCTATGCAGACGATCGTCCGGGTCGCGCCAGcc                           <  2:16322/81‑1 (MQ=255)                                                 agcaACAGCAGCAGCGCCAGCACAATTGACACCGTGTCGCAAATGCGGCTATGCAGACGATCGTCCGGGTCGCGCCAGCCt                          >  2:16374/1‑81 (MQ=255)                                                 agcaACAGCAGCAGCGCCAGCACAATTGACACCGTGTCGCAAATGCGGCTATGCAGACGATCGTCCGGGTCGCGCCAGCCt                          >  2:16373/1‑81 (MQ=255)                                                  gcaACAGCAGCAGCGCCAGCACAATTGACACCGTGTCGCAAATGCGGCTATGCAGACGATCGTCCGGGTCGCGCCAGCCtt                         <  1:16324/81‑1 (MQ=255)                                                  gcaACAGCAGCAGCGCCAGCACAATTGACACCGTGTCGCAAATGCGGCTATGCAGACGATCGTCCGGGTCGCGCCAGCCtt                         <  2:16323/81‑1 (MQ=255)                                                    aaCAGCAGCAGCGCCAGCACAATTGACACCGTGTCGCAAATGCGGCTATGCAGACGATCGTCCGGGTCGCGCCAGCCTTgc                       <  1:16325/81‑1 (MQ=255)                                                     aCAGCAGCAGCGCCAGCACAATTGACACCGTGTCGCAAATGCGGCTATGCAGACGATCGTCCGGGTCGCGCCAGCCTTgca                      <  1:16326/81‑1 (MQ=255)                                                        gcagcagcGCCAGCACAATTGACACCGTGTCGCAAATGCGGCTATGCAGACGATCGTCCGGGTCGCGCCAGCCTTgcag                     >  2:16368/1‑79 (MQ=255)                                                            cagcGCCAGCACAATTGACACCGTGTCGCAAATGCGGCTATGCAGACGATCGTCCGGGTCGCGCCAGCCTTGCAGCAGCgt               >  2:16390/1‑81 (MQ=255)                                                               cgcCAGCACAATTGACACCGTGTCGCAAATGCGGCTATGCAGACGATCGTCCGGGTCGCGCCAGCCTTGCAGCAGCgtggt            >  1:16379/1‑81 (MQ=255)                                                               cgcCAGCACAATTGACACCGTGTCGCAAATGCGGCTATGCAGACGATCGTCCGGGTCGCGCCAGCCTTGCAGCAGCgtggt            <  2:16327/81‑1 (MQ=255)                                                                  cAGCACAATTGACACCGTGTCGCAAATGCGGCTATGCAGACGATCGTCCGGGTCGCGCCAGCCTTGCAGCAGCGTGGTGcc         <  2:16329/81‑1 (MQ=255)                                                                  cAGCACAATTGACACCGTGTCGCAAATGCGGCTATGCAGACGATCGTCCGGGTCGCGCCAGCCTTGCAGCAGCGTGGTGcc         <  1:16328/81‑1 (MQ=255)                                                                  cAGCACAATTGACACCGTGTCGCAAATGCGGCTATGCAGACGATCGTCCGGGTCGCGCCAGCCTTGCAGCAGCGTGGTGcc         >  2:16376/1‑81 (MQ=255)                                                                  cAGCACAATTGACACCGTGTCGCAAATGCGGCTATGCAGACGATCGTCCGGGTCGCGCCAGCCTTGCAGCAGCGTGGTGcc         >  2:16387/1‑81 (MQ=255)                                                                    gCACAATTGACACCGTGTCGCAAATGCGGCTATGCAGACGATCGTCCGGGTCGCGCCAGCCTTGCAGCAGCGTGGTGCCgg       >  2:16365/1‑81 (MQ=255)                                                                     cacaATTGACACCGTGTCGCAAATGCGGCTATGCAGACGATCGTCCGGGTCGCGCCAGCCTTGCAGCAGCGTGGTGCCggg      <  2:16330/81‑1 (MQ=255)                                                                      acaATTGACACCGTGTCGCAAATGCGGCTATGCAGACGATCGTCCGGGTCGCGCCAGCCTTGCAGCAGCGTGGTGCCGGgc     <  2:16331/81‑1 (MQ=255)                                                                         aTTGACACCGTGTCGCAAATGCGGCTATGCAGACGATCGTCCGGGTCGCGCCAGCCTTGCAGCAGCGTGGTGCCGGgcgc   >  2:16382/1‑80 (MQ=255)                                                                         aTTGACACCGTGTCGCAAATGCGGCTATGCAGACGATCGTCCGGGTCGCGCCAGCCTTGCAGCAGCGTGGTGCCGGGCGCa  <  2:16332/81‑1 (MQ=255)                                                                         aTTGACACCGTGTCGCAAATGCGGCTATGCAGACGATCGTCCGGGTCGCGCCAGCCTTGCAGCAGCGTGGTGCCGGGCGCa  <  2:16333/81‑1 (MQ=255)                                                                            |                                                                              GGCAACTGGCGATTTACCCCATCGACGATCACCGCCAGTAACGGTGGCAGCAACAGCAGCAGCGCCAGCACAATTAACACCGTGTCGCAAATGCGGCTATGCAGACGATCGTCCGGGTCGCGCCAGCCTTGCAGCAGCGTGGTGCCGGGCGCA  >  NC_000913/73566‑73718

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