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

TGTGATTGGTATTTCGATGGCCGTCCTTGACGGCGCAATCGCCAACGTCGCCCTGCCAACAATCGCCACGGACCTTCATGCCACGCCAGCCAGTTCCATCTGGGTAGTGAACGCCTATCAAATCGCCATTGTCATCTCCCTGCTCTCGTTTTCGTTTCTGGGCGATATGTTTGGCTATCGACGTATTTATAAATGCGGTCTGGTCGTTTTTCTGTTGTCTTCACTGTTCTGCGCCCTTTCTGATTCGCTGCAAATGCTC  >  NC_000913/1910338‑1910596                                                                                                                                         |                                                                                                                           tgtgATTGGTATTTCGATGGCCGTCCTTGACGGCGCAATCGCCAACGTCGCCCTGCCAACAATCGCCACGGACCTTCATGCCACGCCAGCCAGTTCCATCTGGGTAGTGAACGCCTATCAAATCGCCATTGTCATcccc                                                                                                                          >  2:1020049/1‑139 (MQ=255)     aTTGGTATTTCGATGGCCGTCCTTGACGGCGCAATCGCCGACGTCGCCCTGCCAACAATCGCCACGGACCTTCATGCCACGCCAGCCAGTTCCATCTGGGTAGTGAACGCCTATCAAATCGCCATTGTCATCCCCCTGc                                                                                                                      <  1:1557521/139‑1 (MQ=255)                   ggCCGTCCTTGACGGCGCAATCGCCAACGTCGCCCTGCCAACAATCGCCACGGACCTTCATGCCACGCCAGCCAGTTCCATCTGGGTAGTGAACGCCTATCAAATCGCCATTGTCATCCCCCTGCTCtcgttttcgttt                                                                                                        >  2:154908/1‑139 (MQ=255)                    gCCGTCCTTGACGGCGCAATCGCCAACGTCGCCCTGCCAACAATCGCCACGGACCTTCATGCCACGCCAGCCAGTTCCATCTGGGTAGTGAACGCCTATCAAATCGCCATTGTCATCCCCCTGCTCTCGTTTTCGTTTc                                                                                                       >  2:1790796/1‑139 (MQ=255)                          cTTGACGGCGCAATCGCCAACGTCGCCCTGCCAACAATCGCCACGGACCTTCATGCCACGCCAGCCAGTTCCATCTGGGTAGTGAACGCCTATCAAATCGCCATTGTCATCCCCCTGCTCTCGTTTTCGTTTCTGGGCg                                                                                                 <  1:1519843/139‑1 (MQ=255)                                   gcAATCGCCAACGTCGCCCTGCCAACAATCGCCACGGACCTTCATGCCACGCCAGCCAGTTCCATCTGGGTAGTGAACGCCTATCAAATCGCCATTGTCATCCCCCTGCTCtcgttttcgttt                                                                                                        >  2:829884/1‑123 (MQ=255)                                   gcAATCGCCAACGTCGCCCTGCCAACAATCGCCACGGACCTTCATGCCACGCCAGCCAGTTCCATCTGGGTAGTGAACGCCTATCAAATCGCCATTGTCATCCCCCTGCTCtcgttttcgttt                                                                                                        <  1:829884/123‑1 (MQ=255)                                    cAATCGCCAACGTCGCCCTGCCAACAATCGCCACGGACCTTCATGCCACGCCAGCCAGTTCCATCTGGGTAGTGAACGCCTATCAAATCGCCATTGTCATCCCCCTGCTCTCGTTTTCGTTTCTGGGCGATATGTTTgg                                                                                       <  2:515828/139‑1 (MQ=255)                                        cGCCAACGTCGCCCTGCCAACAATCGCCACGGACCTTCATGCCACGCCAGCCAGTTCCATCTGGGTAGTGAACGCCTATCAAATCGCCATTGTCATCCCCCTGCTCTCGTTTTCGTTTCTGGGCGATATGTTTGGCTAt                                                                                   <  1:1020049/139‑1 (MQ=255)                                           cAACGTCGCCCTGCCAACAATCGCCACGGACCTTCATGCCACGCCAGCCAGTTCCATCTGGGTAGTGAACGCCTATCAAATCGCCATTGTCATCCCCCTGCTCTCGTTTTCGTTTCTGGGCGATATGTTTGGCTATCGa                                                                                <  1:826220/139‑1 (MQ=255)                                            aaCGTCGCCCTGCCAACAATCGCCACGGACCTTCATGCCACGCCAGCCAGTTCCATCTGGGTAGTGAACGCCTATCAAATCGCCATTGTCATCCCCCTGctct                                                                                                                   <  2:423501/103‑1 (MQ=255)                                            aaCGTCGCCCTGCCAACAATCGCCACGGACCTTCATGCCACGCCAGCCAGTTCCATCTGGGTAGTGAACGCCTATCAAATCGCCATTGTCATCCCCCTGctct                                                                                                                   >  1:423501/1‑103 (MQ=255)                                            aaCGTCGCCCTGCCAACAATCGCCACGGACCTTCATGCCACGCCAGCCAGTTCCATCTGGGTAGTGAACGCCTATCAAATCGCCATTGTCATCCCCCTGCTCTCGTTTTCGTTTCTGGGCGATATGTTTGGCTATCGAc                                                                               <  2:866193/139‑1 (MQ=255)                                            aaCGTCGCCCTGCCAACAATCGCCACGGACCTTCATGCCACGCCAGCCAGTTCCATCTGGGTAGTGAACGCCTATCAAATCGCCATTGTCATCCCCCTGCTCTCGTTTTCGTTTCTGGGCGATATGTTTGGCTATCGAc                                                                               <  2:1647565/139‑1 (MQ=255)                                            aaCGTCGCCCTGCCAACAATCGCCACGGACCTTCATGCCACGCCAGCCAGTTCCATCTGGGTAGTGAACGCCTATCAAATCGCCATTGTCATCCCCCTGCTCTCGTTTTCGTTTCTGGGCGATATGTTTGGCTATCGAc                                                                               <  1:597638/139‑1 (MQ=255)                                            aaCGTCGCCCTGCCAACAATCGCCACGGACCTTCATGCCACGCCAGCCAGTTCCATCTGGGTAGTGAACGCCTATCAAATCGCCATTGTCATCCCCCTGCTCTCGTTTTCGTTTCTGGGCGATATGTTTGGCTATCGAc                                                                               <  1:154908/139‑1 (MQ=255)                                                  gCCCTGCCAACAATCGCCACGGACCTTCATGCCACGCCAGCCAGTTCCATCTGGGTAGTGAACGCCTATCAAATCGCCATTGTCATCCCCCTGCTCTCGTTTTCGTTTCTGGGCGATATGTTTGGCTATCGACGTAttt                                                                         >  1:1543483/1‑139 (MQ=255)                                                   cccTGCCAACAATCGCCACGGACCTTCATGCCACGCCAGCCAGTTCCATCTGGGTAGTGAACGCCTATCAAATCGCCATTGTCATCCCCCTGCTCTCGTTTTCGTTTCTGGGCGa                                                                                                >  2:65834/1‑115 (MQ=255)                                                   cccTGCCAACAATCGCCACGGACCTTCATGCCACGCCAGCCAGTTCCATCTGGGTAGTGAACGCCTATCAAATCGCCATTGTCATCCCCCTGCTCTCGTTTTCGTTTCTGGGCGa                                                                                                >  1:1668073/1‑115 (MQ=255)                                                   cccTGCCAACAATCGCCACGGACCTTCATGCCACGCCAGCCAGTTCCATCTGGGTAGTGAACGCCTATCAAATCGCCATTGTCATCCCCCTGCTCTCGTTTTCGTTTCTGGGCGa                                                                                                <  2:1668073/115‑1 (MQ=255)                                                   cccTGCCAACAATCGCCACGGACCTTCATGCCACGCCAGCCAGTTCCATCTGGGTAGTGAACGCCTATCAAATCGCCATTGTCATCCCCCTGCTCTCGTTTTCGTTTCTGGGCGa                                                                                                <  1:65834/115‑1 (MQ=255)                                                               tCGCCACGGACCTTCATGCCACGCCAGCCAGTTCCATCTGGGTAGTGAACGCCTATCAAATCGCCATTGTCATCCCCCTGCTCTCGTTTTCGTTTCTGGGCGATATGTTTGGCTATCGACGTATTta                                                                        <  1:1065780/127‑1 (MQ=255)                                                               tCGCCACGGACCTTCATGCCACGCCAGCCAGTTCCATCTGGGTAGTGAACGCCTATCAAATCGCCATTGTCATCCCCCTGCTCTCGTTTTCGTTTCTGGGCGATATGTTTGGCTATCGACGTATTta                                                                        >  2:1065780/1‑127 (MQ=255)                                                                 gCCACGGACCTTCATGCCACGCCAGCCAGTTCCATCTGGGTAGTGAACGCCTATCAAATCGCCATTGTCATCCCCCTGCTCTCGTTTTCGTTTCTGGGCGATATGTTTGGCTATCGACGTATTTATAAATGCGGTCTgg                                                          <  1:1790796/139‑1 (MQ=255)                                                                  ccACGGACCTTCATGCCACGCCAGCCAGTTCCATCTGGGTAGTGAACGCCTATCAAATCGCCATTGTCGTCCCCCTGCTCTCGTTTTCGTTTCTGGGCGat                                                                                               >  2:1278693/1‑101 (MQ=255)                                                                  ccACGGACCTTCATGCCACGCCAGCCAGTTCCATCTGGGTAGTGAACGCCTATCAAATCGCCATTGTCATCCCCCTGCTCTCGTTTTCGTTTCTGGGCGat                                                                                               <  1:1278693/101‑1 (MQ=255)                                                                  ccACGGACCTTCATGCCACGCCAGCCAGTTCCATCTGGGTAGTGAACGCCTATCAAATCGCCATTGTCATCCCCCTGCTCTCGTTTTCGTTTCTGGGCGATATGTTTGGCTATCGACGTATTTATAAATGCGGTCTGGt                                                         >  1:1269272/1‑139 (MQ=255)                                                                       gACCTTCATGCCACGCCAGCCAGTTCCATCTGGGTAGTGAACGCCTATCAAATCGCCATTGTCATCCCCCTGCTCTCGTTTTCGTTTCTGGGCGATATGTTTGGCTATCGAc                                                                               <  1:1200725/112‑1 (MQ=255)                                                                       gACCTTCATGCCACGCCAGCCAGTTCCATCTGGGTAGTGAACGCCTATCAAATCGCCATTGTCATCCCCCTGCTCTCGTTTTCGTTTCTGGGCGATATGTTTGGCTATCGAc                                                                               >  2:1200725/1‑112 (MQ=255)                                                                                    cGCCAGCCAGTTCCATCTGGGTAGTGAACGCCTATCAAATCGCCATTGTCATCCCCCTGCTCTCGTTTTCGTTTCTGGGCg                                                                                                 >  1:1022168/1‑81 (MQ=255)                                                                                    cGCCAGCCAGTTCCATCTGGGTAGTGAACGCCTATCAAATCGCCATTGTCATCCCCCTGCTCTCGTTTTCGTTTCTGGGCg                                                                                                 <  2:1022168/81‑1 (MQ=255)                                                                                         gccagTTCCATCTGGGTAGTGAACGCCTATCAAATCGCCATTGTCATCCCCCTGCTCTCGTTTTCGTTTCTGGGCGATATGTTTGGCTATCGACGTATTTATAAATGCGGTCTgg                                                          <  1:1055426/115‑1 (MQ=255)                                                                                         gccagTTCCATCTGGGTAGTGAACGCCTATCAAATCGCCATTGTCATCCCCCTGCTCTCGTTTTCGTTTCTGGGCGATATGTTTGGCTATCGACGTATTTATAAATGCGGTCTgg                                                          >  2:1055426/1‑115 (MQ=255)                                                                                                   tCTGGGTAGTGAACGCCTATCAAATCGCCATTGTCATCCCCCTGCTCtcg                                                                                                                 >  2:1548750/1‑50 (MQ=255)                                                                                                   tCTGGGTAGTGAACGCCTATCAAATCGCCATTGTCATCCCCCTGCTCtcg                                                                                                                 <  1:1548750/50‑1 (MQ=255)                                                                                                   tCTGGGTAGTGAACGCCTATCAAATCGCCATTGTCATCCCCCTGCTCTCGTTTTCGTTTCTGGGCGATATGTTTGGCTATCg                                                                                 >  2:555852/1‑82 (MQ=255)                                                                                                   tCTGGGTAGTGAACGCCTATCAAATCGCCATTGTCATCCCCCTGCTCTCGTTTTCGTTTCTGGGCGATATGTTTGGCTATCg                                                                                 <  1:555852/82‑1 (MQ=255)                                                                                                      gggTAGTGAACGCCTATCAAATCGCCATTGTCATCCCCCTGCTCTCGTTTTCGTTTCTGGGCGATATGTTTGGCTATCGAc                                                                               <  1:1613503/81‑1 (MQ=255)                                                                                                      gggTAGTGAACGCCTATCAAATCGCCATTGTCATCCCCCTGCTCTCGTTTTCGTTTCTGGGCGATATGTTTGGCTATCGAc                                                                               >  2:1613503/1‑81 (MQ=255)                                                                                                       ggTAGTGAACGCCTATCAAATCGCCATTGTCATCCCCCTGCTCTCGTTTTCGTTTCTGGGCGATATGTTTGGCTATCGACGTATTTATAAATGCGGTCTGGTCGTTTTTCTGTTGTCTTCACTGTTCTGCGCCCTTTCt                    >  2:499744/1‑139 (MQ=255)                                                                                                        gTAGTGAACGCCTATCAAATCGCCATTGTCATCCCCCTGCTCTCGTTTTCGTTTCTGGGCGATATGTTTGGCTATCGACGTATTTATAAATGCGGTCTGGTCGTTTTTCTGTTGTCTTCACTGTTCTGCGCCCTTTCTg                   >  2:1041223/1‑139 (MQ=255)                                                                                                              aaCGCCTATCAAATCGCCATTGTCATCCCCCTGCTCTCGTTTTCGTTTCTGGGCg                                                                                                 <  2:1742175/55‑1 (MQ=255)                                                                                                              aaCGCCTATCAAATCGCCATTGTCATCCCCCTGCTCTCGTTTTCGTTTCTGGGCg                                                                                                 >  1:1742175/1‑55 (MQ=255)                                                                                                                 gCCTATCAAATCGCCATTGTCATCCCCCTGCTCTCGTTTTCGTTTCTGGGCGATATGTTTGGCTATCGACGt                                                                             >  2:1969524/1‑72 (MQ=255)                                                                                                                 gCCTATCAAATCGCCATTGTCATCCCCCTGCTCTCGTTTTCGTTTCTGGGCGATATGTTTGGCTATCGACGt                                                                             <  1:1969524/72‑1 (MQ=255)                                                                                                                 gCCTATCAAATCGCCATTGTCATCCCCCTGCTCTCGTTTTCGTTTCTGGGCGATATGTTTGGCTATCGACGTATTTATAAATGCGGTCTGGTCGTTTTTCTGTTGTCTTCACTGTTCTGCGCCCTTTCTGATTCGCTGc          <  1:499744/139‑1 (MQ=255)                                                                                                                      tCAAATCGCCATTGTCATCCCCCTGCTCTCGTTTTCGTTTCTGGGCGATATGTTTgg                                                                                       <  2:1810554/57‑1 (MQ=255)                                                                                                                      tCAAATCGCCATTGTCATCCCCCTGCTCTCGTTTTCGTTTCTGGGCGATATGTTTgg                                                                                       >  1:1810554/1‑57 (MQ=255)                                                                                                                      tCAAATCGCCATTGTCATCCCCCTGCTCTCGTTTTCGTTTCTGGGCGATATGTTTGGCTATCg                                                                                 >  1:441604/1‑63 (MQ=255)                                                                                                                      tCAAATCGCCATTGTCATCCCCCTGCTCTCGTTTTCGTTTCTGGGCGATATGTTTGGCTATCg                                                                                 <  2:441604/63‑1 (MQ=255)                                                                                                                         aagcGCCATTGTCATCCCCCTGCTCTCGTTTTCGTTTCTGGGCGATATGTTTGGCTATCGACGTATTTATACATGCGGTCTGGTCGTTTTTCTGTTGTCTTCACTGTTCTGCGCCCTTTCTGATTCGCTGCAAATGCTc  <  2:1543483/136‑1 (MQ=255)                                                                                                                             gCCATTGTCATCCCCCTGCTCTCGTTTTCGTTTCTGGGCGATATGTTTGGCTATCGACg                                                                              <  1:1149628/59‑1 (MQ=255)                                                                                                                             gCCATTGTCATCCCCCTGCTCTCGTTTTCGTTTCTGGGCGATATGTTTGGCTATCGACg                                                                              >  2:1149628/1‑59 (MQ=255)                                                                                                                                         |                                                                                                                           TGTGATTGGTATTTCGATGGCCGTCCTTGACGGCGCAATCGCCAACGTCGCCCTGCCAACAATCGCCACGGACCTTCATGCCACGCCAGCCAGTTCCATCTGGGTAGTGAACGCCTATCAAATCGCCATTGTCATCTCCCTGCTCTCGTTTTCGTTTCTGGGCGATATGTTTGGCTATCGACGTATTTATAAATGCGGTCTGGTCGTTTTTCTGTTGTCTTCACTGTTCTGCGCCCTTTCTGATTCGCTGCAAATGCTC  >  NC_000913/1910338‑1910596

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