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

ATTGGTATTTCGATGGCCGTCCTTGACGGCGCAATCGCCAACGTCG‑CCCTGCCAACAATCGCCACGGACCTTCATGCCACGCCAGCCAGTTCCATCTGGGTAGTGAACGCCTATCAAATCGCCATTGTCATCTCCCTGCTCTCGTTTTCGTTTCTGGGCGATATGTTTGGCTATCGACGTATTTATAAATGCGGTCTGGTCGTTTTTCTGTTGTCTTCACTGTTCTGCGCCCTTTCTGATTCGCTGCAAATGCTCACCCTTGC  >  NC_000913/1910342‑1910604                                                                                                                                      |                                                                                                                                   aTTGGTATTTCGATGGCCGTCCTTGACGGCGCAATCGCCAACGTCG‑CCCTGCCAACAATCGCCACGGACCTTCATGCCACGCCAGCCAGTTCCATCTGGGTAGTGAACGCCTATCAAATCGCCATTGTCATCCCCCTGc                                                                                                                              <  1:307973/139‑1 (MQ=255) aTTGGTATTTCGATGGCCGTCCTTGACGGCGCAATCGCCAACGTCG‑CCCTGCCAACAATCGCCACGGACCTTCATGCCACGCCAGCCAGTTCCATCTGGGTAGTGAACGCCTATCAAATCGCCATTGTCATCCCCCTGc                                                                                                                              >  2:404816/1‑139 (MQ=255)  ttGGTATTTCGATGGCCGTCCTTGACGGCGCAATCGCCAACGTCG‑CCCTGCCAACAATCGCCACGGACCTTCATGCCACGCCAGCCAGTTCCATCTGGGTAGTGAACGCCTATCAAATCGCCATTGTCATCCCCCTGct                                                                                                                             >  2:407281/1‑139 (MQ=255)         ttCGATGGCCGTCCTTGATGGCGCAATCGCCAACGTCG‑CCCTGCCAACAATCGCCACGGACCTTCATGCCACGCCAGCCAGTTCCATCTGGGTAGTGAACGCCTATCAAATCGCCATTGTCATCCCCCTGCTCtcgttt                                                                                                                      >  2:426793/1‑139 (MQ=255)         ttCGATGGCCGTCCTTGACGGCGCAATCGCCAACGTCG‑CCCTGCCAACAATCGCCACGGACCTTCATGCCACGCCAGCCAGTTCCATCTGGGTAGTGAACGCCTATCAAATCGCCATTGTCATCCCCCTGCTCtcgttt                                                                                                                      >  1:57048/1‑139 (MQ=255)           cGATGGCCGTCCTTGACGGCGCAATCGCCAACGTCG‑CCCTGCCAACAATCGCCACGGACCTTCATGCCACGCCAGCCAGTTCCATCTGGGTAGTGAACGCCTATCAAATCGCCATTGTCATCCCCCTGCTCtcgttttc                                                                                                                    >  1:372517/1‑139 (MQ=255)                gCCGTCCTTGACGGCGCAATCGCCAACGTCG‑CCCTGCCAACAATCGCCACGGACCTTCATGCCACGCCAGCCAGTTCCATCTGGGTAGTGAACGCCTATCAAATCGCCATTGTCATCCCCCTGCTCTCGTTTTCGTTTc                                                                                                               <  2:651136/139‑1 (MQ=255)                    tCCTTGACGGCGCAATCGCCAACGTCG‑CCCTGCCAACAATCGCCACGGACCTTCATGCCACGCCAGCCAGTTCCATCTGGGTAGTGAACGCCTATCAAATCGCCATTGTCATCCCCCTGCTCTCGTTTTCGTTTCTggg                                                                                                           >  2:61300/1‑139 (MQ=255)                              cgcAATCGCCAACGTCGCCCCTGCCAACAATCGCCACGGACCTTCATGCCACGCCAGCCAGTTCCATCTGGGTAGTGAACGCCTATCAAATCGCCATTGTCATCCCCCTGCTCTCGTTTTCGTTTCTGGGCGATATGtt                                                                                                  >  2:593503/1‑139 (MQ=255)                                  aTCGCCAACGCCG‑CCCTGCCAACAATCGCCACGGACCTTCATGCCACGCCAGCCAGTTCCATCTGGGTAGTGAACTCCTATCAAATCTCCATTGTCATCCCCCTGCTCTCGTTTTCGTTTCTGGGCGATATGTTTGGCt                                                                                             <  1:737362/139‑1 (MQ=255)                                    cGCCAACGTCG‑CCCTGCCAACAATCGCCACGGACCTTCATGCCACGCCAGCCAGTTCCATCTGGGTAGTGAACGCCTATCAAATCGCCATTGTCATCCCCCTGCTCTCGTTTTCGTTTCTGGGCGATATGTTTGGCTAt                                                                                           <  1:61300/139‑1 (MQ=255)                                         aCGTCG‑CCCTGCCAACAATCGCCACGGACCTTCATGCCACGCCAGCCAGTTCCATCTGGGTAGTGAACGCCTATCAAATCGCCATTGTCATCCCCCTGCTCTCGTTTTCGTTTCTGGGCGATATGTTTGGCTATCGACg                                                                                      <  1:33068/139‑1 (MQ=255)                                                                    gACCTTCATGCCACGCCAGCCAGTTCCATCTGGGTAGTGAACGCCTATCAAATCGCCATTGTCATCCCCCTGCTCTCGTTTTCGTTTCTGGGCGATATGTTTGGCTATCGACGTATTTATAAATGCGGTCTGGTCGttt                                                            <  2:57048/139‑1 (MQ=255)                                                                      ccTTCATGCCACGCCAGCCAGTTCCATCTGGGTAGTGAACGCCTATCAAATCGCCATTGTCATCCCCCTGCTCTCGTTTTCGTTTCTGGGCGATATGTTTGGCTATCGACGTATTTATAAATGCGGTCTGGTCGttttt                                                          <  2:749538/139‑1 (MQ=255)                                                                             gCCACGCCAGCCAGTTCCATCTGGGTAGTGAACGCCTATCAAATCGCCATTGTCATCCCCCTGCTCTTGTTTTCGTTTCTGGGCGATATGTTTGGCTATCGACGTATTTATAAATGCGGTCTGGTCGTTTTTCtgttgt                                                   <  1:458665/139‑1 (MQ=255)                                                                             gCCACGCCAGCCAGTTCCATCTGGGTAGTGAACGCCTATCAAATCGCCATTGTCATCCCCCTGCTCTCGTTTTCGTTTCTGGGCGATATGTTTGGCTATCGACGTATTTATAAATGCGGTCTGGTCGTTTTTCtgttgt                                                   <  1:339103/139‑1 (MQ=255)                                                                              ccACGCCAGCCAGTTCCATCTGGGTAGTGAACGCCTATCAAATCGCCATTGTCATCCCCCTGCTCTCGTTTTCGTTTCTGGGCGATATGTTTGGCTATCGACGTATTTATAAATGCGGTCTGGTCGTTTTTCTGTTGTc                                                  >  1:385536/1‑139 (MQ=255)                                                                                 cGCCAGCCAGTTCCATCTGGGTAGTGAACGCCTATCAAATCGCCATTGTCATCCCCCTGCTCTCGTTTTCGTTTCTGGGCGATATGTTTGGCTATCGACGTATTTATAAATGCGGTCTGGTCGTTTTTCTGTTGTCTTc                                               >  1:65733/1‑139 (MQ=255)                                                                                   ccagccagTTCCATCTGGGTAGTGAACGCCTATCAAATCGCCATTGTCATCCCCCTGCTCTCGTTTTCGTTTCTGGGCGATATGTTTGGCTATCGACGTATTTATAAATGCGGTCTGGTCGTTTTTCTGTTGTCTTCAc                                             >  2:755153/1‑139 (MQ=255)                                                                                       ccagTTCCATCTGGGTAGTGAACGCCTATCAAATCGCCATTGTCATCCCCCTGCTCTCGTTTTCGTTTCTGGGCGATATGTTTgg                                                                                               <  2:767620/85‑1 (MQ=255)                                                                                       ccagTTCCATCTGGGTAGTGAACGCCTATCAAATCGCCATTGTCATCCCCCTGCTCTCGTTTTCGTTTCTGGGCGATATGTTTgg                                                                                               >  1:767620/1‑85 (MQ=255)                                                                                       ccagTTCCATCTGGGTAGTGAACGCCTATCAAATCGCCATTGTCATCCCCCTGCTCTCGTTTTCGTTTCTGGGCGATATGTTTGGCTATCGACGTATTTATAAATGCGGTCTGGTCGTTTTTCTGTTGTCTTCACTGtt                                         <  1:426793/139‑1 (MQ=255)                                                                                             ccATCTGGGTAGTGAACGCCTATCAAATCGCCATTGTCATCCCCCTGCTCTCGTTTTCGTTTCTGGGCGATATGTTTGGCTATCGACGTATTTATAAATGCGGTCTGGTCGTTTTTCTGTTGTCTTCACTGTTCTgcgc                                   <  1:599858/139‑1 (MQ=255)                                                                                                        gTGAACGCCTATCAAATCGCCATTGTCATCCCCCTGCTCTCGTTTTCGTTTCTGGGCGATATGTTTGGCTATCGa                                                                                        <  2:473245/75‑1 (MQ=255)                                                                                                        gTGAACGCCTATCAAATCGCCATTGTCATCCCCCTGCTCTCGTTTTCGTTTCTGGGCGATATGTTTGGCTATCGa                                                                                        >  1:473245/1‑75 (MQ=255)                                                                                                         tGAACGCCTATCAAATCGCCATTGTCATCCCCCTGCTCTCGTTTTCGTTTCTGGGCGATATGTTTGGCTATCGACGTATTTATAAATGCGGTCTGGTCGTTTTTCTGTTGTCTTCACTGTTCTGCGCCCTTTCTGATTc                       <  2:372517/139‑1 (MQ=255)                                                                                                           aaCGCCTATCAAATCGCCATTGTCATCCCCCTGCTCTCGTTTTCGTTTCTGGGCGATATGTTTGGCTATCGACGTATTTATAAATGCGGTCTGGTCGTTTTTCTGTTGTCTTCACTGTTCTGCGCCCTTTCTGATTCGc                     >  1:609679/1‑139 (MQ=255)                                                                                                              gtcTATCAAATCGCCATTGTCATCCCCCTGCTCTCGTTTTCGTTTCTGGGCGATATGTTTGGCTATCGACGTATTTATAAATGCGGTCTGGTCGTTTTTCtgttgt                                                   <  2:102972/104‑1 (MQ=255)                                                                                                              gCCTATCAAATCGCCATTGTCATCCCCCTGCTCTCGTTTTCGTTTCTGGGCGATATGTTTGGCTATCGACGTATTTATAAATGCGGTCTGGTCGTTTTTCtgttgt                                                   >  1:102972/1‑106 (MQ=255)                                                                                                                   tCAAATCGCCATTGTCATCCCCCTGCTCTCGTTTTCGTTTCTGGGCGATATGTTTGGCTATCGACGTATTTATAAATGCGGTCTGGTCGTTTTTCTGTTGTCTTCACTGTTCTGCGCCCTTTCTGATTCGCTGCAAATg             >  2:660013/1‑139 (MQ=255)                                                                                                                     gacaggCCATTGTTATACCCCTCTTCTCGTTTTCGTTTCTGGGCGCTATGTTTGGCTATCGACTCATTTCTAAATGCGGTCGGGTCGTGTTTCTGTTGCCTTCACCGTTCTCCGCCCTTTCTGATTCGCTGCAAATGCt           <  1:478080/134‑1 (MQ=255)                                                                                                                          gCCATTGTCATCCCCCTGCTCTCGTTTTCGTTTCCGGGCGATATGTTTGGCTATCGACGTATTTATAAATGCGGTCTGGTCGTTTTTCTGTTGTCTTCACTGTTCTGCGCCCTTTCTGATTCGCTGCAAATGCTCAccc      <  2:14849/139‑1 (MQ=255)                                                                                                                              ttGTCATCCCCCTGCTCTCGTTTTCGTTTCTGGGCGATATGTTTGGCTATCGACGTATTTATAAATGCGGTCTGGTCGTTTTTCTGTTGTCTTCACTGTTCTGCGCCCTTTCTGATTCGCTGCAAATGCTCACCCTTgc  <  2:65733/139‑1 (MQ=255)                                                                                                                                      |                                                                                                                                   ATTGGTATTTCGATGGCCGTCCTTGACGGCGCAATCGCCAACGTCG‑CCCTGCCAACAATCGCCACGGACCTTCATGCCACGCCAGCCAGTTCCATCTGGGTAGTGAACGCCTATCAAATCGCCATTGTCATCTCCCTGCTCTCGTTTTCGTTTCTGGGCGATATGTTTGGCTATCGACGTATTTATAAATGCGGTCTGGTCGTTTTTCTGTTGTCTTCACTGTTCTGCGCCCTTTCTGATTCGCTGCAAATGCTCACCCTTGC  >  NC_000913/1910342‑1910604

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