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%	104.9 / NA	30	S67P (TCC→CCC)	yebQ	putative transporter
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (14/16);  total (14/16)

TGTGATTGGTATTTCGATGGCCGTCCTTGACGGCGCAATCGCCAACGTCGCCCTGCCAACAATCGCCACGGACCTTCATGCCACGCCAGCCAGTTCCATCTGGGTAGTGAACGCCTATCAAATCGCCATTGTCATCTCCCTGCTCTCG‑TTTTCGTTTCTGGGCGATATGTTTGGCTATCGACGTATTTATAAATGCGGTCTGGTCGTTTTTCTGTTGTCTTCACTGTTCTGCGCCCTTTCTGATTCGCTGCAAATGCTCACCCTTGCGCGTG  >  NC_000913/1910338‑1910609                                                                                                                                         |                                                                                                                                         tgtgATTGGTATTTCGATGGCCGTCCTTGACGGCGCAATCGCCAACGTCGCCCTGCCAACAATCGCCACGGACCTTCATGCCACGCCAGCCAGTTCCATCTGGGTAGTGAACGCCTATCAAATCGCCATTGTCATcccc                                                                                                                                        >  1:901951/1‑139 (MQ=255)     aTTGGTATTTCGATGGCCGTCCTTGACGGCGCAATCGCCAACGTCGCCCTGCCAACAATCGCCACGGACCTTCATGCCACGCCAGCCAGTTCCATCTGGGTAGTGAACGCCTATCAAATCGCCATTGTCATcccc                                                                                                                                        >  1:732859/1‑135 (MQ=255)        ggTATTTCGATGGCCGTCCTTGACGGCGCAATCGCCAACGTCGCCCTGCCAACAATCGCCACGGACCTTCATGCCACGCCAGCCAGTTCCATCTGGGTAGTGAACGCCTATCAAATCGCCATTGTCATCCCCCTGctct                                                                                                                                 >  2:802813/1‑139 (MQ=255)        ggTATTTCGATGGCCGTCCTTGACGGCGCAATCGCCAACGTCGCCCTGCCAACAATCGCCACGGACCTTCATGCCACGCCAGCCAGTTCCATCTGGGTAGTGAACGCCTATCAAATCGCCATTGTCATCCCCCTGctct                                                                                                                                 >  1:601586/1‑139 (MQ=255)           aTTTCGATGGCCGTCCTTGACGGCGCAATCGCCAACGTCGCCCTGCCAACAATCGCCACGGACCTTCATGCCACGCCAGCCAGTTCCATCTGGGTAGTGAACGCCTATCAAATCGCCATTGTCATCCCCCTGCTCtcg‑t                                                                                                                             >  1:1281583/1‑139 (MQ=255)             ttCGATGGCCGTCCTTGACGGCGCAATCGCCAACGTCGCCCTGCCAACAATCGCCACGGACCTTCATGCCACGCCAGCCAGTTCCATCTGGGTAGTGAACGCCTATCAAATCGCCATTGTCATCCCCCTGCTCtcg‑ttt                                                                                                                           <  2:1281583/139‑1 (MQ=255)                   ggCCGTCCTTGACGGCGCAATCGCCAACGTCGCCCTGCCAACAATCGCCACGGACCTTCATGCCACGCCAGCCAGTTCCATCTGGGTAGTGAACGCCTATCAAATCGCCATTGTCATCCCCCTGCTCtcg‑ttttcgttt                                                                                                                     >  1:174907/1‑139 (MQ=255)                    gCCGTCCTTGACGGCGCAATCGCCAACGTCGCCCTGCCAACAATCGCCACGGACCTTCATGCCACGCCAGCCAGTTCCATCTGGGTAGTGAACGCCTATCAAATCGCCATTGTCATCCCCCTGCTCTCG‑TTTTCGTTTc                                                                                                                    <  2:901951/139‑1 (MQ=255)                    gCCGTCCTTGACGGCGCAATCGCCAACGTCGCCCTGCCAACAATCGCCACGGACCTTCATGCCACGCCAGCCAGTTCCATCTGGGTAGTGAACGCCTATCAAATCGCCATTGTCATCCCCCTGCTCTCG‑TTTTCGTTTc                                                                                                                    <  1:772517/139‑1 (MQ=255)                    gCCGTCCTTGACGGCGCAATCGCCAACGTCGCCCTGCCAACAATCGCCACGGACCTTCATGCCACGCCAGCCAGTTCCATCTGGGTAGTGAACGCCTATCAAATCGCCATCGTCATCCCCCTGCTCTCG‑TTTTCGTTTc                                                                                                                    >  2:820181/1‑139 (MQ=255)                     ccGTCCTTGACGGCGCAATCGCCAACGTCGCCCTGCCAACAATCGCCACGGACCTTCATGCCACGCCAGCCAGTTCCATCTGGGTAGTGAACGCCTATCAAATCGCCATTGTCATCCCCCTGCTCTCG‑TTTTCGTTTCt                                                                                                                   <  2:601586/139‑1 (MQ=255)                            tGACGGCGCAATCGCCAACGTCGCCCTGCCAACAATCGCCACGGACCTTCATGCCACGCCAGCCAGTTCCATCTGGGTAGTGAACGCCTATCAAATCGCCATTGTCATCCCCCTGCTCTCG‑TTTTCGTTTCTGGGCGat                                                                                                            >  2:275505/1‑139 (MQ=255)                                  cgcAATCGCCAACGTCGTCCTGCCAACAATCGCCACGGACCTTCATGCCACGCCAGCCAGTTCCATCTGGGTAGTGAACGCCTATCAAATCGCCATTGTCATCCCCCTGCTCTCGTTTTTCGTTTCTGGCCGATATGtt                                                                                                       >  1:935506/1‑139 (MQ=255)                                      aTCGCCAACGTCGCCCTGCCAACAATCGCCACGGACCTTCATGCCACGCCAGCCAGTTCCATCTGGGTAGTGAACGCCTATCAAATCGCCATTGTCATCCCCCTGCTCTCG‑TTTTCGTTTCTGGGCGATATGTTTGGCt                                                                                                  <  1:491849/139‑1 (MQ=255)                                          ccAACGTCGCCCTGCCAACAATCGCCACGGACCTTCATGCCACGCCAGCCAGTTCCATCTGGGTAGTGAACGCCTATCAAATCGCCATTGTCATCCCCCTGCTCTCG‑TTTTCGTTTCTGGGCGATATGTTTGGCTATCg                                                                                              <  1:275505/139‑1 (MQ=255)                                            aaCGTCGCCCTGCCAACAATCGCCACGGACCTTCATGCCACGCCAGCCAGTTCCATCTGGGTAGTGAACGCCTATCAAATCGCCATTGTCATCCCCCTGCTCTCG‑TTTTCGTTTCTGGGCGATATGTTTGGCTATCGAc                                                                                            <  2:703442/139‑1 (MQ=255)                                             aCGTCGTCCTGCCAACAATCGCCACGGACCTTCATGCCACGCCAGCCAGTTCCATCTGGGTAGTGAACGCCTATCAAATCGCCATTGTCATCCCCCTGCTCTCGTTTTTCGTTTCTGGGCGATATGTTTGGCTATCGAc                                                                                            <  2:935506/139‑1 (MQ=255)                                                 cGCCCTGCCAACAATCGCCACGGACCTTCATGCCACGCCAGCCAGTTCCATCTGGGTAGTGAACGCCTATCAAATCGCCATTGTCATCCCCCTGCTCTCG‑TTTTCGTTTCTGGGCGATATGTTTgg                                                                                                    >  1:203985/1‑126 (MQ=255)                                                 cGCCCTGCCAACAATCGCCACGGACCTTCATGCCACGCCAGCCAGTTCCATCTGGGTAGTGAACGCCTATCAAATCGCCATTGTCATCCCCCTGCTCTCG‑TTTTCGTTTCTGGGCGATATGTTTgg                                                                                                    <  2:203985/126‑1 (MQ=255)                                                        ccAACAATCGCCACGGACCTTCATGCCACGCCAGCCAGTTCCATCTGGGTAGTGAACGCCTATCAAATCGCCATTGTCATCCCCCTGCTCTCG‑TTTTCGTTTCTGGGCGATATGTTTGGCTATCGACGTATTTATAAAt                                                                                >  2:221421/1‑139 (MQ=255)                                                                    aCGGACCTTCATGCCACGCCAGCCAGTTCCATCTGGGTAGTGAACGCCTATCAAATCGCCATTGTCATCCCCCTGCTCTCG‑TTTTCGTTTCTGGGCGATATGTTTGGCTATCGACGTATTTATAAATGCGGTCTGGTCg                                                                    <  1:221421/139‑1 (MQ=255)                                                                       gACCTTCATGCCACGCCAGCCAGTTCCATCTGGGTAGTGAACGCCTATCAAATCGCCATTGTCATCCCCCTGCTCTCG‑TTTTCGTTTCTGGGCGATATGTTTGGCTATCGACGTATTTATAAATGCGGTCTGGTCGttt                                                                 >  1:1207348/1‑139 (MQ=255)                                                                                  cACGCCAGCCAGTTCCATCTGGGTAGTGAACGCCTATCAAATCGCCATTGTCATCCCCCTGCTCTCG‑TTTTCGTTTCTGGGCGATATGTTTGGCTATCGACGTATTTATAAATGCGGTCTGGTCGTTTTTCTGTTGTCt                                                      >  2:565232/1‑139 (MQ=255)                                                                                   aagccagccaGTTCCATCTGGGTAGTGAACGCCTATCAAATCGCCATTGTCATCCCCCTGCTCTCG‑TTTTCGTTTCTGGGCGATATGTTTGGCTATCGACGTATTTATAAATGCGGTCCGGTCGTTTTTCTGTTGTCtt                                                     <  1:802813/137‑1 (MQ=255)                                                                                    cGCCAGCCAGTTCCATCTGGGTAGTGAACGCCTATCAAATCGCCATTGTCATCCCCCTGCTCTCG‑TTTTCGTTTCTGGGCGATATGTTTGGCTATCGACGTATTTATAAATGCGGTCTGGTCGTTTTTCTGTTGTCTTc                                                    >  2:410688/1‑139 (MQ=255)                                                                                                 cATCTGGGTAGTGAACGCCTATCAAATCGCCATTGTCATCCCCCTGCTCTCG‑TTTTCGTTTCTGGGCGATATGTTTGGCTATCGACGTATTTATAAATGCGGTCTGGTCGTTTTTCTGTTGTCTTCACTGTTCTGCGcc                                       <  2:320366/139‑1 (MQ=255)                                                                                                   tCTTGGTAGTGAACGCCTATCAAATCGCCATTGTCATCCCCCTGCTCTCG‑TTTTCGTTTCTGGGCGATATGTTTGGCTATCGACGTATTTACAAATGCGGTCTGGTCGTTTTTCTGCTGTCTTCACTGTTCTGCGCCCt                                     <  1:254694/139‑1 (MQ=255)                                                                                                                           tCGCCATTGTCATCCCCCTGCTCTCG‑TTTTCGTTTCTGGGCGATATGTTTGGCTATCGACGTATTTATAAATGCGGTCTGGTCGTTTTTCTGTTGTCTTCACTGTTCTgcgc                                        >  2:454400/1‑112 (MQ=255)                                                                                                                           tCGCCATTGTCATCCCCCTGCTCTCG‑TTTTCGTTTCTGGGCGATATGTTTGGCTATCGACGTATTTATAAATGCGGTCTGGTCGTTTTTCTGTTGTCTTCACTGTTCTgcgc                                        <  1:454400/112‑1 (MQ=255)                                                                                                                                    tCATCCCCCTGCTCTCG‑TTTTCGTTTCTGGGCGATATGTTTGGCTATCGACGTATTTATAAATGCGGTCTGGTCGTTTTTCTGTTGTCTTCACTGTTCTGCGCCCTTTCTGATTCGCTGCAAATGCTCACCCTTGcgcg    >  2:496994/1‑139 (MQ=255)                                                                                                                                     cATCCCCCTGCTCTCG‑TTTTCGTTTCTGGGCGATATGTTTGCCTATCGACGTATTTATAAATGCGGTCTGGTCGTTTTTCTGTTGTCTTCACTGTTCTGCGCCCTTTCTGATTCGCTGCAAATGCTCACCCTTgcgctt   <  1:820181/139‑3 (MQ=255)                                                                                                                                      aTCCCCCTGCTCTCG‑TTTTCGTTTCTGGGCGATATGTTTGGCTATCGACGTATTTATAAATGCGGTCTGGTCGTTTTTCTGTTGTCTTCACTGTTCTGCGCCCTTTCTGATTCGCTGCAAATGCTCACCCTTGCGCgtg  <  2:174907/139‑1 (MQ=255)                                                                                                                                         |                                                                                                                                         TGTGATTGGTATTTCGATGGCCGTCCTTGACGGCGCAATCGCCAACGTCGCCCTGCCAACAATCGCCACGGACCTTCATGCCACGCCAGCCAGTTCCATCTGGGTAGTGAACGCCTATCAAATCGCCATTGTCATCTCCCTGCTCTCG‑TTTTCGTTTCTGGGCGATATGTTTGGCTATCGACGTATTTATAAATGCGGTCTGGTCGTTTTTCTGTTGTCTTCACTGTTCTGCGCCCTTTCTGATTCGCTGCAAATGCTCACCCTTGCGCGTG  >  NC_000913/1910338‑1910609

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