Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,910,474	T→C	S67P (TCC→CCC)	yebQ →	putative transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,910,474	0	T	C	100.0%	129.9 / NA	38	S67P (TCC→CCC)	yebQ	putative transporter
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (14/24);  total (14/24)

TGTGATTGGTATTTCGATGGCCGTCCTTGACGGCGCAATCGCCAACGTCGCCCTGCCAACAATCGCCACGGACCTTCATGCCACGCCAGCCAGTTCCATCTGGGTAGTGAACGCCTATCAAATCGCCATTGTCAT‑CTCCCTGCTCTCGTTTTCGTTTCTGGGCGATATGTTTGGCTATCGACGTATTTATAAATGCGGTCTGGTCGTTTTTCTGTTGTCTTCACTGTTCTGCGCCCTTTCTGATTCGCTGCAAATGCTCACC  >  NC_000913/1910338‑1910599                                                                                                                                          |                                                                                                                              tgtgATTGGTATTTCGATGGCCGTCCTTGACGGCGCAATCGCCAACGTCGCCCTGCCAACAATCGCCACGGACCTTCATGCCACGCCAGCCAGTTCCATCTGGGTAGTGAACGCCTATCAAATCGCCATTGTCAT‑cccc                                                                                                                             >  1:11030/1‑139 (MQ=255)   tgATTGGTATTTCGATGGCCGCCCTTGACGGCGCAATCGCCAACGTCGCCCTGCCAACAATCGCCACGGACCTTCATGCCACGCCAGCCAGTTCCATCTGGGTAGTGAACGCCTATCAAATCGCCATTGTCAT‑CCCCCt                                                                                                                           >  2:1127635/1‑139 (MQ=255)       tGGTATTTCGATGGCCGTCCTTGACGGCGCAATCGCCAACGTCGCCCTGCCAACAATCGCCACGGACCTTCATGCCACGCCAGCCAGTTCCATCTGGGTAGTGAACGCCTATCAAATCGCCATTGTCAT‑CCCCCTGctc                                                                                                                       <  1:820089/139‑1 (MQ=255)            tttCGATGGCCGTCCTTGACGGCGCAATCGCCAACGTCGCCCTGCCAACAATCGCCACGGACCTTCATGCCACGCCAGCCAGTTCCATCTGGGTAGTGAACGCCTATCAAATCGCCATTGTCAT‑CCCCCTGCTCtcgtt                                                                                                                  >  1:699985/1‑139 (MQ=255)              tCGATGGCCGTCCTTGACGGCGCAATCGCCAACGTCGCCCTGCCAACAATCGCCACGGACCTTCATGCCACGCCAGCCAGTTCCATCTGGGTAGTGAACGCCTATCAAATCGCCATTGTCAT‑CCCCCTGCTCtcgtttt                                                                                                                <  2:699985/139‑1 (MQ=255)                          cTTGACGGCGCAATCGCCAACGTCGCCCTGCCAACAATCGCCACGGACCTTCATGCCACGCCAGCCAGTTCCATCTGGGTAGTGAACGCCTATCAAATCGCCATTGTCAT‑CCCCCTGCTCTCGTTTTCGTTTCTGGGCg                                                                                                    <  1:1294457/139‑1 (MQ=255)                          cTTGACGGCGCAATCGCCAACGTCGCCCTGCCAACAATCGCCACGGACCTTCATGCCACGCCAGCCAGTTCCATCTGGGTAGTGAACGCCTATCAAATCGCCATTGTCAT‑CCCCCTGCTCTCGTTTTCGTTTCTGGGCg                                                                                                    <  2:978157/139‑1 (MQ=255)                                          ccAACGTCGCCCTGCCAACAAGCGCCACGGACCTTCATGCCACGCCAGCCAGTTCCATCTGGGTAGTGAACGCCTATCAAATCGCCATTGTCAT‑CCCCCTGCTCTCGTTTTCGTTTCTGGGCGATATGTTTGGCTATCg                                                                                    <  1:630985/139‑1 (MQ=255)                                               gTCGCCCTGCCAACAATCGCCACGGACCTTCATGCCACGCCAGCCAGTTCCATCTGGGTAGTGAACGCCTATCAAATCGCCATTGTCAT‑CCCCCTGCTCTCGTTTTCGTTTCTGGGCGATATGTTTGGCTATCGACGTa                                                                               <  2:203269/139‑1 (MQ=255)                                                         caacaaTCGCCACGGACCTTCATGCCACGCCAGCCAGTTCCATCTGGGTAGTGAACGCCTATCAAATCGCCATTGTCAT‑CCCCCTGCTCTCGTTTTCGTTTCTGGGCGATATGTTTGGCTATCGACGTATTTATaa                                                                        >  2:1139825/1‑136 (MQ=255)                                                         caacaaTCGCCACGGACCTTCATGCCACGCCAGCCAGTTCCATCTGGGTAGTGAACGCCTATCAAATCGCCATTGTCAT‑CCCCCTGCTCTCGTTTTCGTTTCTGGGCGATATGTTTGGCTATCGACGTATTTATaa                                                                        <  1:1139825/136‑1 (MQ=255)                                                              aTCGCCACGGACCTTCATGCCACGCCAGCCAGTTCCATCTGGGTAGTGAACGCCTATCAAATCGCCATTGTCAT‑CCCCCTGCTCTCGTTTTCGTTTCt                                                                                                         <  1:295014/98‑1 (MQ=255)                                                              aTCGCCACGGACCTTCATGCCACGCCAGCCAGTTCCATCTGGGTAGTGAACGCCTATCAAATCGCCATTGTCAT‑CCCCCTGCTCTCGTTTTCGTTTCt                                                                                                         >  2:295014/1‑98 (MQ=255)                                                               tCGCCACGGACCTTCATGCCACGCCAGCCAGTTCCATCTGGGTAGTGAACGCCTATCAAATCGCCATTGTCAT‑CCCCCTGCTCTCGTTTTCGTTTCTGGGCGATATGTTTGGCTATCGACGTATTTATAAATGCGGTCt                                                               <  1:1127635/139‑1 (MQ=255)                                                                  ccACGGACCTTCATGCCACGCCAGCCAGTTCCATCTGGGTAGTGAACGCCTATCAAATCGCCATTGTCAT‑CCCCCTGCTCTCGTTTTCGTTTCTGGGCg                                                                                                    >  2:894768/1‑99 (MQ=255)                                                                  ccACGGACCTTCATGCCACGCCAGCCAGTTCCATCTGGGTAGTGAACGCCTATCAAATCGCCATTGTCAT‑CCCCCTGCTCTCGTTTTCGTTTCTGGGCg                                                                                                    <  1:894768/99‑1 (MQ=255)                                                                  ccACGGACCTTCATGCCACGCCAGCCAGTTCCATCTGGGTAGTGAACGCCTATCAAATCGCCATTGTCAT‑CCCCCTGCTCTCGTTTTCGTTTCTGGGCGATATGTTTGGCTATCGACGTATTTATAAATGCGGTCTGGt                                                            <  2:898038/139‑1 (MQ=255)                                                                      ggACCTTCATGCCACGCCAGCCAGTTCCATCTGGGTAGTGAACGCCTATCAAATCGCCATTGTCAT‑CCCCCTGCTCTCGTTTTCGTTTCTGGGCGa                                                                                                   <  2:759559/96‑1 (MQ=255)                                                                      ggACCTTCATGCCACGCCAGCCAGTTCCATCTGGGTAGTGAACGCCTATCAAATCGCCATTGTCAT‑CCCCCTGCTCTCGTTTTCGTTTCTGGGCGa                                                                                                   >  1:759559/1‑96 (MQ=255)                                                                             cATGCCACGCCAGCCAGTTCCATCTGGGTAGTGAACGCCTATCAAATCGCCATTGTCAT‑CCCCCTGCTCTCGTTTTCGTTTCTGGGCGATATGTTTGGCTATCGACGTATTTATAAATGCGGTCTGGTCGTTTTTCtgt                                                 >  1:1147981/1‑139 (MQ=255)                                                                              aTGCCACGCCAGCCAGTTCCATCTGGGTAGTGAACGCCTATCAAATCGCCATTGTCAT‑CCCCCTGCTCTCGTTTTCGTTTCTGGGCGATATGTTTGGCTATCGACGTATTTATAAATGCGGTCTGGTCGTTTTTCtgtt                                                <  2:979119/139‑1 (MQ=255)                                                                                    cGCCAGCCAGTTCCATCTGGGTAGTGAACGCCTATCAAATCGCCATTGTCAT‑CCCCCTGCTCTCGTTTTCGTTTCTGGGCGa                                                                                                   <  1:1260151/82‑1 (MQ=255)                                                                                    cGCCAGCCAGTTCCATCTGGGTAGTGAACGCCTATCAAATCGCCATTGTCAT‑CCCCCTGCTCTCGTTTTCGTTTCTGGGCGa                                                                                                   >  2:1260151/1‑82 (MQ=255)                                                                                         gccagTTCCATCTGGGTAGTGAACGCCTATCAAATCGCCATTGTCAT‑CCCCCTGCTCTCGTTTTCGTTTCTGGGCGATATGTTTGGCTATCGACGTATTTATAAATGCGGTCTGGTCGTTTTTCTGTTGTCTTCACTGt                                     >  1:987109/1‑139 (MQ=255)                                                                                          ccagTTCCATCTGGGTAGTGAACGCCTATCAAATCGCCATTGTCAT‑CCCCCTGCTCTCGTTTTCGTTTCTGGGCGATATGTTTGGCTATCGACGTATTTATAAATGCGGTCTGGTCGTTTTTCTGTTGTCTTCACTGtt                                    <  2:11030/139‑1 (MQ=255)                                                                                               tCCATCTGGGTAGTGAACGCCTATCAAATCGCCATTGTCAT‑CCCCCTGCTCTCGTTTTCGTTTCTGGGCGATATGTTTGGCTATCGACGTATTTATAAATGCGGTCTGGTCGTTTTTCTGTTGTCTTCACTGTTCTgcg                               <  1:684749/139‑1 (MQ=255)                                                                                               tCCATCTGGGTAGTGAACGCCTATCAAATCGCCATTGTCAT‑CCCCCTGCTCTCGTTTTCGTTTCTGGGCGATATGTTTGGCTATCGACGTATTTATAAATGCGGTCTGGTCGTTTTTCTGTTGTCTTCACTGTTCTgcg                               <  1:180180/139‑1 (MQ=255)                                                                                                   tCTGGGTAGTGAACGCCTATCAAATCGCCATTGTCATCCCCCCTGCTCTCGTTTTCGTTTCTGGGCGATATGt                                                                                              >  2:1036688/1‑73 (MQ=255)                                                                                                   tCTGGGTAGTGAACGCCTATCAAATCGCCATTGTCATCCCCCCTGCTCTCGTTTTCGTTTCTGGGCGATATGt                                                                                              <  1:1036688/73‑1 (MQ=255)                                                                                                               aCGCCTATCAAATCGCCATTGTCAT‑CCCCCTGCTCTCGTTTTCGTTTCTGGGCGATATGTTTGGCTATCg                                                                                    >  2:1257308/1‑70 (MQ=255)                                                                                                               aCGCCTATCAAATCGCCATTGTCAT‑CCCCCTGCTCTCGTTTTCGTTTCTGGGCGATATGTTTGGCTATCg                                                                                    <  1:1257308/70‑1 (MQ=255)                                                                                                                   cTATCAAATCGCCATTGTCAT‑CCCCCTGCTCTCGTTTTCGTTTCTGGGCGATATGTTTGGCTATCGACGTATTTATAAATGCGGTCTGGTCGTTTTTCTGTTGTCTTCACTGTTCTGCGCCCTTTCTGATTCGCTGCaa           <  2:1147981/139‑1 (MQ=255)                                                                                                                      tCAAATCGCCATTGTCAT‑CCCCCTGCTCTCGTTTTCGTTTCTGGGCGATATGTTTGGCTATCg                                                                                    >  1:820765/1‑63 (MQ=255)                                                                                                                      tCAAATCGCCATTGTCAT‑CCCCCTGCTCTCGTTTTCGTTTCTGGGCGATATGTTTGGCTATCg                                                                                    <  2:820765/63‑1 (MQ=255)                                                                                                                            cGCCATTGTCAT‑CCCCCTGCTCTCGTTTTCGTTTCTGGGCGATATGTTTGGCTATCGACGTATTTATAAATGCGGTCTGGTCGTTTTTCTGTTGTCTTCACTGTTCTGCGCCCTTTCTGATTCGCTGCAAATGCTCAcc  <  2:987109/139‑1 (MQ=255)                                                                                                                             gCCATTGTCAT‑CCCCCTGCTCTCGTTTTCGTTTCTGGGCGATATGTTTGGCTATCGACGTATTTATAAATGCGGTCTGGTCGTTTTTCTGTTGTCTTCACTGTTCTGCGCCCTTTCCGATTCGCTGCAAATGCt      <  2:442949/134‑1 (MQ=255)                                                                                                                             gCCATTGTCAT‑CCCCCTGCTCTCGTTTTCGTTTCTGGGCGATATGTTTGGCTATCGACGTATTTATAAATGCGGTCTGGTCGTTTTTCTGTTGTCTTCACTGTTCTGCGCCCTTTCCGATTCGCTGCAAATGCt      >  1:442949/1‑134 (MQ=255)                                                                                                                               cATTGTCAT‑CCCCCTGCTCTCGTTTTCGTTTCTGGGCGATATGTTTGGCTATCGACGTATTta                                                                           >  2:701764/1‑63 (MQ=255)                                                                                                                               cATTGTCAT‑CCCCCTGCTCTCGTTTTCGTTTCTGGGCGATATGTTTGGCTATCGACGTATTta                                                                           <  1:701764/63‑1 (MQ=255)                                                                                                                                          |                                                                                                                              TGTGATTGGTATTTCGATGGCCGTCCTTGACGGCGCAATCGCCAACGTCGCCCTGCCAACAATCGCCACGGACCTTCATGCCACGCCAGCCAGTTCCATCTGGGTAGTGAACGCCTATCAAATCGCCATTGTCAT‑CTCCCTGCTCTCGTTTTCGTTTCTGGGCGATATGTTTGGCTATCGACGTATTTATAAATGCGGTCTGGTCGTTTTTCTGTTGTCTTCACTGTTCTGCGCCCTTTCTGATTCGCTGCAAATGCTCACC  >  NC_000913/1910338‑1910599

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