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

TGTGATTGGTATTTCGATGGCCGTCCTTGACGGCGCAATCGCCAACGTCGCCCTGCCAACAATCGCCACGGACCTTCATGCCACGCCAGCCAGTTCCATCTGGGTAGTGAACGCCTATCAAATCGCCATTGTCATCTCCCTGCTCTCGTTTTCGTTTCTGGGCGATATGTTTGGCTATCGACGTATTTATAAATGCGGTCTGGTCGTTTTTCTGTTGTCTTCACTGTTCTGCGCCCTTTCTGATTCGCTGCAAATGCTCACCCTTGCGCGT  >  NC_000913/1910338‑1910608                                                                                                                                         |                                                                                                                                       tgtgATTGGTATTTCGATGGCCGTCCTTGACGGCGCAATCGCCAACGTCGCCCTGCCAACAATCGCCACGGACCTTCATGCCACGCCAGCCAGTTCCATCTGGGTAGTGAACGCCTATCAAATCGCCATTGTCATcccc                                                                                                                                      >  1:650960/1‑139 (MQ=255) tgtgATTGGTATTTCGATGGCCGTCCTTGACGGCGCAATCGCCAACGTCGCCCTGCCAACAATCGCCACGGACCTTCATGCCACGCCAGCCAGTTCCATCTGGGTAGTGAACGCCTATCAAATCGCCATTGTCATcccc                                                                                                                                      <  1:406676/139‑1 (MQ=255)     aTTGGTATTTCGATGGCCGTCCTTGACGGCGCAATCGCCAACGTCGCCCTGCCAACAATCGCCACGGACCTTCATGCCACGCCAGCCAGTTCCATCTGGGTAGTGAACGCCTATCAAATCGCCATTGTCACCCCCCTGc                                                                                                                                  <  1:495596/139‑1 (MQ=255)       tGGTATTTCGATGGCCGTCCTTGACGGCGCAATCGCCAACGTCGCCCTGCCAACAATCGCCACGGACCTTCATGCCACGCCAGCCAGTTCCATCTGGGTAGTGAACGCCTATCAAATCGCCATTGTCATCCCCCTGctc                                                                                                                                >  2:1111877/1‑139 (MQ=255)        ggTATTTCGATGGCCGTCCTTGACGGCGCAATCGCCAACGTCGCCCTGCCAACAATCGCCACGGACCTTCATGCCACGCCAGCCAGTTCCATCTGGGTAGTGAACGCCTATCAAATCGCCATTGTCATCCCCCTGctct                                                                                                                               >  1:358054/1‑139 (MQ=255)             ttCGATGGCCGTCCTTGACGGCGCAATCGCCAACGTCGCCCTGCCAACAATCGCCACGGACCTTCATGCCACGCCAGCCAGTTCCATCTGGGTAGTGAACGCCTATCAAATCGCCATTGTCATCCCCCTGCTCtcgttt                                                                                                                          >  1:1323071/1‑139 (MQ=255)               cGATGGCCGTCCTTGACGGCGCAATCGCCAACGTCGCCCTGCCAACAATCGCCACGGACCTTCATGCCACGCCAGCCAGTTCCATCTGGGTAGTGAACGCCTATCAAATCGCCATTGTCATCCCCCTGCTCtcgttttc                                                                                                                        >  1:698839/1‑139 (MQ=255)                 aTGGCCGTCCTTGACGGCGCAATCGCCAACGTCGCCCTGCCAACAATCGCCACGGACCTTCATGCCACGCCAGCCAGTTCCATCTGGGTAGTGAACGCCTATCAAATCGCCATTGTCATCCCCCTGCTCtcgttttcgt                                                                                                                      >  1:864553/1‑139 (MQ=255)                 aTGGCCGTCCTTGACGGCGCAATCGCCAACGTCGCCCTGCCAACAATCGCCACGGACCTTCATGCCACGCCAGCCAGTTCCATCTGGGTAGTGAACGCCTATCAAATCGACATTGTCATCCCCCTGCTCtcgttttcgt                                                                                                                      >  1:1392779/1‑139 (MQ=255)                        tCCTTGACGGCGCAATCGCCAACGTCGCCCTGCCAACAATCGCCACGGACCTTCATGCCACGCCAGCCAGTTCCATCTGGGTAGTGAACGCCTATCAAATCGCCATTGTCATCCCCCTGCTCTCGTTTTCGTTTCTggg                                                                                                               >  1:402992/1‑139 (MQ=255)                          cTTGACGGCGCAATCGCCAACGTCGCCCTGCCAACAATCGCCACGGACCTTCATGCCACGCCAGCCAGTTCCATCTGGGTAGTGAACGCCTATCAAATCGCCATTGTCATCCCCCTGCTCTCGTTTTCGTTTCTGGGCg                                                                                                             <  2:1323071/139‑1 (MQ=255)                          cTTGACGGCGCAATCGCCAACGTCGCCCTGCCAACAATCGCCACGGACCTTCATGCCACGCCAGCCAGTTCCATCTGGGTAGTGAACGCCTATCAAATCGCCATTGTCATCCCCCTGCTCTCGTTTTCGTTTCTGGGCg                                                                                                             <  1:1111877/139‑1 (MQ=255)                          cTTGACGGCGCAATCGCCAACGTCGCCCTGCCAACAATCGCCACGGACCTTCATGCCACGCCAGCCAGTTCCATCTGGGTAGTGAACGCCTATCAAATCGCCATTGTCATCCCCCTGCTCTCGTTTTCGTTTCTGGGCg                                                                                                             <  2:358054/139‑1 (MQ=255)                          cTTGACGGCGCAATCGCCAACGTCGCCCTGCCAACAATCGCCACGGACCTTCATGCCACGCCAGCCAGTTCCATCTGGGTAGTGAACGCCTATCAAATCGCCATTGTCATCCCCCTGCTCTCGTTTTCGTTTCTGGGCg                                                                                                             <  2:698839/139‑1 (MQ=255)                                  cgcAATCGCCAACGTCGCCCTGCCAACAATCGCCACGGACCTTCATGCCACGCCAGCCAGTTCCATCTGGGTAGTGAACGCCTATCAAATCGCCATTGTCATCCCCCTGCTCTCGTTTTCGTTTCTGGGCg                                                                                                             <  2:129472/131‑1 (MQ=255)                                  cgcAATCGCCAACGTCGCCCTGCCAACAATCGCCACGGACCTTCATGCCACGCCAGCCAGTTCCATCTGGGTAGTGAACGCCTATCAAATCGCCATTGTCATCCCCCTGCTCTCGTTTTCGTTTCTGGGCg                                                                                                             >  1:129472/1‑131 (MQ=255)                                            aaCGTCGCCCTGCCAACAATCGCCACGGACCTTCATGCCACGCCAGCCAGTTCCATCTGGGTAGTGAACGCCTATCAAATCGCCATTGTCATCCCCCTGCTCTCGTTTTCGTTTCTGGGCGATATGTTTGGCTATCGAc                                                                                           <  2:864553/139‑1 (MQ=255)                                                         caacaaTCGCCACGGACCTTCATGCCACGCCAGCCAGTTCCATCTGGGTAGTGAACGCCTATCAAATCGCCATTGTCATCCCCCTGCTCTCGTTTTCGTTTCTGGGCGATATGTTTGGCTATCGACGTATTTATAAATg                                                                              <  2:230072/139‑1 (MQ=255)                                                               tCGCCACGGACCTTCATGCCACGCCAGCCAGTTCCATCTGGGTAGTGAACGCCTATCAAATCGCCATTGTCATCCCCCTGCTCTCGTTTTCGTTTCTGGGCGATATGTTTGGCTATCGACGTATTTATAAATGCGGTCt                                                                        <  2:402992/139‑1 (MQ=255)                                                                  ccACGGACCTTCATGCCACGCCAGCCAGTTCCATCTGGGTAGTGAACGCCTATCAAATCGCCATTGTCATCCCCCTGCTCTCGTTTTCGTTTCTGGGCGATATGTTTGGCTATCGa                                                                                            <  1:735628/116‑1 (MQ=255)                                                                  ccACGGACCTTCATGCCACGCCAGCCAGTTCCATCTGGGTAGTGAACGCCTATCAAATCGCCATTGTCATCCCCCTGCTCTCGTTTTCGTTTCTGGGCGATATGTTTGGCTATCGa                                                                                            >  2:735628/1‑116 (MQ=255)                                                                  ccACGGACCTTCATGCCACGCCAGCCAGTTCCATCTGGGTAGTGAACGCCTATCAAATCGCCATTGTCATCCCCCTGCTCTCGTTTTCGTTTCTGGGCGATATGTTTGGCTATCGACGTATTTATAAATGCGGTCTGGt                                                                     <  2:353738/139‑1 (MQ=255)                                                                  ccACGGACCTTCATGCCACGCCAGCCAGTTCCATCTGGGTAGTGAACGCCTATCAAATCGCCATTGTCATCCCCCTGCTCTCGTTTTCGTTTCTGGGCGATATGTTTGGCTATCGACGTATTTATAAATGAGGTCTGGt                                                                     >  2:864635/1‑139 (MQ=255)                                                                           ttCATGCCACGCCAGCCAGTTCCATCTGGGTAGTGAACGCCTATCAAATCGCCATTGTCATCCCCCTGCTCTCGTTTTCGTTTCTGGGCGATATGTTTGGCTATCGACGTATTTATAAATGCGGTCTGGTCGTTTTTCt                                                            <  1:368591/139‑1 (MQ=255)                                                                              aTGCCACGCCAGCCAGTTCCATCTGGGTGGTGAACGCCTATCAAATCGCCATTGTCATCCCCCTGCTCTCGTTTTCGTTTCTGGGCGATATGTTTGGCTATCGACGTATTTATAAATGCGGTCTGGTCGTTTTTCtgtt                                                         <  1:276890/139‑1 (MQ=255)                                                                                                       ggTAGTGAACGCCTATCAAATCGCCATTGTCATCCCCCTGCTCTCGTTTTCGTTTCTGGGCGATATGTTTGGCTATCGACGTATTTATAAATGCGGTCTGGTCGTTTTTCTGTTGTCTTCACTGTTCTgcg                                        >  1:247033/1‑131 (MQ=255)                                                                                                       ggTAGTGAACGCCTATCAAATCGCCATTGTCATCCCCCTGCTCTCGTTTTCGTTTCTGGGCGATATGTTTGGCTATCGACGTATTTATAAATGCGGTCTGGTCGTTTTTCTGTTGTCTTCACTGTTCTgcg                                        <  2:247033/131‑1 (MQ=255)                                                                                                        gTAGTGAACGCCTATCAAATCGCCATTGTCATCCCCCTGCTCTCGTTTTCGTTTCTGGGCGATATGTTTGGCTATCGACGTATTTATAAATGCGGTCTGGTCGTTTTTCTGTTGTCTTCACTGTTCTGCGCCCTTTCTg                               >  1:231650/1‑139 (MQ=255)                                                                                                             gAACGCCTATCAAATCGCCATTGTCATCCCCCTGCTCTCGTTTTCGTTTCTGGGCGATATGTTTGGCTATCGACGTATTTATAAATGCGGTCTGGTCGTTTTTCTGTTGTCTTCACTGTTCTgcg                                        <  2:1200731/125‑1 (MQ=255)                                                                                                             gAACGCCTATCAAATCGCCATTGTCATCCCCCTGCTCTCGTTTTCGTTTCTGGGCGATATGTTTGGCTATCGACGTATTTATAAATGCGGTCTGGTCGTTTTTCTGTTGTCTTCACTGTTCTgcg                                        >  1:1200731/1‑125 (MQ=255)                                                                                                              aaCGCCTATCAAATCGCCATTGTCATCCCCCTGCTCTCGTTTTCGTTTCTGGGCGATATGTTTGGCTATCGACGTATTTATAAATGCGGTCTGGTCGTTTTTCTGTTGTCTTCACTGTTCTGCGCCCTTTCTGATTCGc                         >  2:220188/1‑139 (MQ=255)                                                                                                                 gCCTATCAAATCTCCATTGTCATCCCCCTGCTCTCGTTTTCGTTTCTGGGCGATATGTTTGGCTATCGACGTATTTATAAATGCGGTCTGGTCGTTTTTCtgttgt                                                       <  2:708795/106‑1 (MQ=255)                                                                                                                 gCCTATCAAATCGCCATTGTCATCCCCCTGCTCTCGTTTTCGTTTCTGGGCGATATGTTTGGCTATCGACGTATTTATAAATGCGGTCTGGTCGTTTTTCtgttgt                                                       >  1:708795/1‑106 (MQ=255)                                                                                                                  ccTATCAAATCGCCATTGTCATCCCCCTGCTCTCGTTTTCGTTTCTGGGCGATATGTTTGGCTATCGACGTATTTATAAATGCGGTCTGGTCGTTTTTCTGTTGTCTTCACTGTTCTGCGCCCTTTCTGATTCGCTGCa                     <  2:1392779/139‑1 (MQ=255)                                                                                                                      tCAAATCGCCATTGTCATCCCCCTGCTCTCGTTTTCGTTTCTGGGCGATATGTTTGGCTATCGACGTATTTATAAATGCGGTCt                                                                        >  1:149169/1‑84 (MQ=255)                                                                                                                      tCAAATCGCCATTGTCATCCCCCTGCTCTCGTTTTCGTTTCTGGGCGATATGTTTGGCTATCGACGTATTTATAAATGCGGTCt                                                                        <  2:149169/84‑1 (MQ=255)                                                                                                                                aTTGTCATCCCCCTGCTCTCGTTTTCGTTTCTGGGCGATATGTTTGGCTATCGACGTATTTATAAATGCGGTCTGGTCGTTTTTCt                                                            >  2:992408/1‑86 (MQ=255)                                                                                                                                aTTGTCATCCCCCTGCTCTCGTTTTCGTTTCTGGGCGATATGTTTGGCTATCGACGTATTTATAAATGCGGTCTGGTCGTTTTTCt                                                            <  1:992408/86‑1 (MQ=255)                                                                                                                                     cATCCCCCTGCTCTCGTTTTCGTTTCTGGGCGATATGTTTGGCTATCGACGTATTTATAAATGCGGTCTGGTCGTTTTTCTGTTGTCTTCACTGTTCTGCGCCCTTTCTGATTCGCTGCAAATGCTCACCCTTGCGCgt  <  2:650960/139‑1 (MQ=255)                                                                                                                                         |                                                                                                                                       TGTGATTGGTATTTCGATGGCCGTCCTTGACGGCGCAATCGCCAACGTCGCCCTGCCAACAATCGCCACGGACCTTCATGCCACGCCAGCCAGTTCCATCTGGGTAGTGAACGCCTATCAAATCGCCATTGTCATCTCCCTGCTCTCGTTTTCGTTTCTGGGCGATATGTTTGGCTATCGACGTATTTATAAATGCGGTCTGGTCGTTTTTCTGTTGTCTTCACTGTTCTGCGCCCTTTCTGATTCGCTGCAAATGCTCACCCTTGCGCGT  >  NC_000913/1910338‑1910608

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