Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,826,476	T→C	N428S (AAC→AGC)	b1746 ←	putative aldehyde dehydrogenase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,826,476	0	T	C	100.0%	252.6 / NA	79	N428S (AAC→AGC)	b1746	putative aldehyde dehydrogenase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (33/46);  total (33/46)

GCATACCAGGCGCTGGGGCGATGGTTACCGGAAGCACCAATGCCGCCGAATGGCGCGGTACTGGCAGCACCGGTAAGCGGTTTGTTCCAGTTAACAATCCCCGCCCGCGCCTCCAGCAACAGTTGATCGAACTTTTCCCGCTCGGGGGAAACCAGACCGCAAGAGAGTCCGAAGCGAGT  >  NC_000913/1826386‑1826564                                                                                           |                                                                                         gCATACCAGGCGCTGGGGCGATGGTTACCGGAAGCACCAATGCCGCCGAATGGCGCGGTACTGGCAGCACCGGTAAGCGGTTTGTTCCAGCTa                                                                                        >  2:4300518/1‑93 (MQ=255)    tACCAGGCGCTGGGGCGATGGTTACCGGAAGCACCAATGCCGCCGAATGGCGCGGTACTGGCAGCACCGGTAAGCGGTTTGTTCCAGCTAACa                                                                                     <  2:1356114/93‑1 (MQ=255)    tACCAGGCGCTGGGGCGATGGTTACCGGAAGCACCAATGCCGCCGAATGGCGCGGTACTGGCAGCACCGGTAAGCGGTTTGTTCCAGCTAACa                                                                                     >  1:3358768/1‑93 (MQ=255)     aCCAGGCGCTGGGGCGATGGTTACCGGAAGCACCAATGCCGCCGAATGGCGCGGTACTGGCAGCACCGGTAAGCGGTTTGTTCCAGCTAACaa                                                                                    <  1:1548481/93‑1 (MQ=255)     aCCAGGCGCTGGGGCGATGGTTACCGGAAGCACCAATGCCGCCGAATGGCGCGGTACTGGCAGCACCGGTAAGCGGTTTGTTCCAGCTAACaa                                                                                    <  1:4024046/93‑1 (MQ=255)       cAGGCGCTGGGGCGATGGTTACCGGAAGCACCAATGCCGCCGAATGGCGCGGTACTGGCAGCACCGGTAAGCGGTTTGTTCCAGCTAACAATc                                                                                  <  2:3086544/93‑1 (MQ=255)       cAGGCGCTGGGGCGATGGTTACCGGAAGCACCAATGCCGCCGAATGGCGCGGTACTGGCAGCACCGGTAAGCGGTTTGTTCCAGCTAACAATc                                                                                  >  2:637287/1‑93 (MQ=255)          gcgcTGGGGCGATGGTTACCGGAAGCACCAATGCCGCCGAATGGCGCGGTACTGGCAGCACCGGTAAGCGGTTTGTTCCAGCTAACAATcccc                                                                               <  1:464382/93‑1 (MQ=255)            gcTGGGGCGATGGTTACCGGAAGCACCAATGCCGCCGAATGGCGCGGTACTGGCAGCACCGGTAAGCGGTTTGTTCCAGCTAACAATCcccgc                                                                             >  2:2396585/1‑93 (MQ=255)             gtGGGTCGATGGTTACCGGAAGCACCAATGCCGCCGAATGGCGCGGTACTGGCAGCACCGGTAAGCGGTTTGTTCCAGCTAACAACCcccgcc                                                                            <  2:3170295/92‑1 (MQ=255)              tGGGGCGATGGTTACCGGAAGCACCAATGCCGCCGAATGGCGCGGTACTGGCAGCACCGGTAAGCGGTTTGTTCCAGCTAACAATCcccgccc                                                                           >  1:1050240/1‑93 (MQ=255)               ggggCGATGGTTACCGGAAGCACCAATGCCGCCGAATGGCGCGGTACTGGCAGCACCGGTAAGCGGTTTGTTCCAGCTAACAATCcccgcccg                                                                          >  1:1707483/1‑93 (MQ=255)                  gCGATGGTTACCGGAAGCACCAATGCCGCCGAATGGCGCGGTACTGGCAGCACCGGTAAGCGGTTTGTTCCAGCTAACAATCCCCGCCCgcgc                                                                       <  2:2706420/93‑1 (MQ=255)                   cGATGGTTACCGGAAGCACCAATGCCGCCGAATGGCGCGGTACTGGCAGCACCGGTAAGCGGTTTGTTCCAGCTAACAATCCCCGCCCGCGcc                                                                      <  1:4440766/93‑1 (MQ=255)                   cGATGGTTACCGGAAGCACCAATGCCGCCGAATGGCGCGGTACTGGCAGCACCGGTAAGCGGTTTGTTCCAGCTAACAATCCCCGCCCGCGcc                                                                      >  2:196519/1‑93 (MQ=255)                    gATGGTTACCGGAAGCACCAATGCCGCCGAATGGCGCGGTACTGGCAGCACCGGTAAGCGGTTTGTTCCAGCTAACAATCCCCGCCCGCGCCt                                                                     >  2:1882089/1‑93 (MQ=255)                    gATGGTTACCGGAAGCACCAATGCCGCCGAATGGCGCGGTACTGGCAGCACCGGTAAGCGGTTTGTTCCAGCTAACAATCCCCGCCCGCGCCt                                                                     >  1:2032783/1‑93 (MQ=255)                      tGGTTACCGGAAGCACCAATGCCGCCGAATGGCGCGGTACTGGCAGCACCGGTAAGCGGTTTGTTCCAGCTAACAATCCCCGCCCGCGCCtac                                                                   <  2:2077643/93‑3 (MQ=255)                         ttACCGGAAGCACCAATGCCGCCGAATGGCGCGGTACTGGCAGCACCGGTAAGCGGTTTGTTCCAGCTAACAATCCCCGCCCGCGCCTCCAGc                                                                <  2:1602399/93‑1 (MQ=255)                           aCCGGAAGCACCAATGCCGCCGAATGGCGCGGTACTGGCAGCACCGGTAAGCGGTTTGTTCCAGCTAACAATCCCCGCCCGCGCCTCCAGCaa                                                              >  1:912318/1‑93 (MQ=255)                               gAAGCACCAATGCCGCCGAATGGCGCGGTACTGGCAGCACCGGTAAGCGGTTTGTTCCAGCTAACAATCCCCGCCCGCGCCTCCAGCAACAGt                                                          >  1:863040/1‑93 (MQ=255)                                  gCACCAATGCCGCCGAATGGCGCGGTACTGGCAGCACCGGTAAGCGGTTTGTTCCAGCTAACAATCCCCGCCCGCGCCTCCAGCAACAGTTGa                                                       <  2:1071523/93‑1 (MQ=255)                                    aCCAATGCCGCCGAATGGCGCGGTACTGGCAGCACCGGTAAGCGGTTTGTTCCAGCTAACAATCCCCGCCCGCGCCTCCAGCAACAGTTGATc                                                     <  1:432093/93‑1 (MQ=255)                                     ccAATGCCGCCGAATGGCGCGGTACTGGCAGCACCGGTAAGCGGTTTGTTCCAGCTAACAATCCCCGCCCGCGCCTCCAGCAACAGTTGATCg                                                    >  2:3685164/1‑93 (MQ=255)                                       aaTGCCGCCGAATGGCGCGGTACTGGCAGCACCGGTAAGCGGTTTGTTCCAGCTAACAATCCTCGCCGGCGCCTCCAGCAACAGTTGATCGag                                                  >  2:2776948/1‑92 (MQ=255)                                       aaTGCCGCCGAATGGCGCGGTACTGGCAGCACCGGTAAGCGGTTTGTTCCAGCTAACAATCCCCGCCCGCGCCTCCAGCAACAGTTGATCGaa                                                  <  1:4392371/93‑1 (MQ=255)                                          gccgccGAATGGCGCGGTACTGGCAGCACCGGTAAGCGGTTTGTTCCAGCTAACAATCCCCGCCCGCGCCTCCAGCAACAGTTGATCGAACtt                                               <  1:390598/93‑1 (MQ=255)                                           ccgccgAATGGCGCGGTACTGGCAGCACCGGTAAGCGGTTTGTTCCAGCTAACAATCCCCGCCCGCGCCTCCAGCAACAGTTGATCGAACttt                                              <  1:190110/93‑1 (MQ=255)                                           ccgccgAATGGCGCGGTACTGGCAGCACCGGTAAGCGGTTTGTTCCAGCTAACAATCCCCGCCCGCGCCTCCAGCAACAGTTGATCGAACttt                                              <  1:4392052/93‑1 (MQ=255)                                            cgccgAAGGGCGCGGTACTGGCAGCACCGGTAAGCGGTTTGTTCCAGCTAACAATCCCCGCCCGCGCCTCCAGCAACAGTTGATCGAACtttt                                             <  2:3846444/93‑1 (MQ=255)                                             gccgAATGGCGCGGTACTGGCAGCACCGGTAAGCGGTTTGTTCCAGCTAACAATCGCCGCCCGCGCCTCCAGCAACAGTTGATCGAACTTTTc                                            >  2:4077528/1‑93 (MQ=255)                                             gccgAATGGCGCGGTACTGGCAGCACCGGTAAGCGGTTTGTTCCAGCTAACAATCCCCGCCCGCGCCTCCAGCAACAGTTGATCGAACTTTTc                                            >  2:1905249/1‑93 (MQ=255)                                             gccgAATGGCGCGGTACTGGCAGCACCGGTAAGCGGTTTGTTCCAGCTAACAATCCCCGCCCGCGCCTCCAGCAACAGTTGATCGAACTTTTc                                            <  1:916175/93‑1 (MQ=255)                                              ccgAATGGCGCGGTACTGGCAGCACCGGTAAGCGGTTTGTTCCAGCTAACAATCCCCGCCGGGGCCTCCAGCAAAAGTTGATCGAACTTTTcc                                           >  1:3222926/1‑93 (MQ=255)                                              ccgAATGGCGCGGTACTGGCAGCACCGGTAAGCGGTTTGTTCCAGCTAACAATCCCCGCCCGCGCCTCCAGCAACAGTTGATCGAACTTTTcc                                           >  1:2385232/1‑93 (MQ=255)                                              ccgAATGGCGCGGTACTGGCAGCACCGGTAAGCGGTTTGTTCCAGCTAACAATCCCCGCCCGCGCCTCCAGCAACAGTTGATCGAACTTTTcc                                           <  2:3769937/93‑1 (MQ=255)                                               cgAATGGGGCGGTACTGGCAGCACCGGTAAGCGGTTTGTTCCAGCTAACAATCCCCGCCCGCGCCTCCAGCAACAGTTGATCGAACTTTTccc                                          <  1:3702033/93‑1 (MQ=255)                                               cgAATGGCGCGGTACTGGCAGCACCGGTAAGCGGTTTGTTCCAGCTAACAATCCCCGCCCGCGCCTCCAGCAACAGTTGATCGAACTTTTccc                                          >  2:4171527/1‑93 (MQ=255)                                               cgAATGGCGCGGTACTGGCAGCACCGGTAAGCGGTTTGTTCCAGCTAACAATCCCCGCCCGCGCCTCCAGCAACAGTTGATCGAACTTTTccc                                          <  1:3254643/93‑1 (MQ=255)                                                 aagggCGCGGTACTGGCAGCACCGGTAAGCGGTTTGTTCCAGCTAACAATCCCCGCCCGCGCCTCCAGCAACAGTTGATCGAACTTTTCCCGc                                        <  2:4028525/90‑1 (MQ=255)                                                   tGGCGCGGTACTGGCAGCACCGGTAAGCGGTTTGTTCCAGCTAACAATCCCCGCCCGCGCCTCCAGCAACAGTTGATCGAACTTTTCCCGCTc                                      <  2:3293921/93‑1 (MQ=255)                                                    ggCGCGGTACTGGCAGCACCGGTAAGCGGTTTGTTCCAGCTAACAATCCCCGCCCGCGCCTCCAGCAACAGTTGATCGAACTTTTCCCGCTCg                                     <  1:4077423/93‑1 (MQ=255)                                                     gcgcGGTACTGGCAGCACCGGTAAGCGGTTTGTTCCAGCTAACAATCCCCGCCCGCGCCTCCAGCAACAGTTGATCGAACTTTTCCCGCTCgg                                    >  1:3343195/1‑93 (MQ=255)                                                       gcgGTACTGGCAGCACCGGTAAGCGGTTTGTTCCAGCTAACAATCCCCGCCCGCGCCTCCAGCAACAGTTGATCGAACTTTTCCCGCTCgggg                                  >  1:1104203/1‑93 (MQ=255)                                                        cgGTACTGGCAGCACCGGTAAGCGGTTTGTTCCAGCTAACAATCCCCGCCCGCGCCTCCAGCAACAGTTGATCGAACTTTTCCCGCTCggggg                                 <  2:2859430/93‑1 (MQ=255)                                                           tACTGGCATCACCGTTAAGCGGTTTGTTCCAGCAAGCAATCCCCGCCCGCGCCTCCAGCAACAGTTGATCGAGCTTTTCCCGCTCGGGGGaaa                              <  1:2400391/93‑1 (MQ=255)                                                           tACTGGCAGCACCGGTAAGCGGTTTGTTCCAGCTAACAATCCCCGCCCGCGCCTCCAGCAACAGTTGATCGAACTTTTCCCGCTCGGGGGaaa                              <  2:3163789/93‑1 (MQ=255)                                                            aCTGGCAGCACCGGTAAGCGGTTTGTTCCAGCTAACAATCCCCGCCCGCTCCTCCAGCAACAGTTGATCGAACTTTTCCCGCTCGGGGGAAAc                             <  1:324248/93‑1 (MQ=255)                                                            aCTGGCAGCACCGGTAAGCGCTTTGTTCCAGCTAACAATCCCCGCCCGCGCCTCCAGCAACAGTTGATCGAACTTTTCCCGCTCGGGGGAAAc                             <  1:610296/93‑1 (MQ=255)                                                             cTGGCAGCACCGGTAAGCGGTTTGTTCCAGCTAACAATCCCCGCCCGCGCCTCCAGCAACAGTTGATCGAACTTTTCCCGCTCGGGGGAAAcc                            <  2:3400478/93‑1 (MQ=255)                                                               ggCAGCACCGGTAAGCGGTTTTTTCCAGCTAACAATCCCCCCCCGCGCCTCCAGCAACAGTTGATCGAACTTTTCCCGCTCGGGGGAAACCAg                          <  1:2828077/93‑1 (MQ=255)                                                                gcagcaCCGGTAAGCGGTTTGTTCCAGCTAACAATCCCCGCCCGCGCCTCCAGCAACAGTTGATCGAACTTTTCCCGCTCGGGGGAAACCAGa                         <  1:1569000/93‑1 (MQ=255)                                                                 cagcaCCGGTAAGCGGTTTGTTCCAGCTAACAATCCCCGCCCGCGCCTCCAGCAACAGTTGATCGAACTTTTCCCGCTCGGGGGAAACCAGAc                        <  2:4427780/93‑1 (MQ=255)                                                                  cacaCCGGTAAGCGGTTTGTTCCAGCTAACAATCCCCGCCCGCGCCTCCAGCAACAGTTGATCGAACTTTTCCCGCTCGGGGGAAACCAGAcc                       <  2:4008628/91‑1 (MQ=255)                                                                   gcaCCGGTAAGCGGTTTGTTCCAGCTAACAATCCCCGCCCGCGCCTCCAGCAACAGTTGATCGAACTTTTCCCGCTCGGGGGAAACCAGACCg                      >  1:4039482/1‑93 (MQ=255)                                                                    caCCGGTAAGCGGTTTGTTCCAGCTAACAATCCCCGCCCGCGCCTCCAGCAACAGTTGATCGAACTTTTCCCGCTCGGGGGAAACCAGACCGc                     <  2:3398547/93‑1 (MQ=255)                                                                     aCCGGTAAGCGGTTTGTTCCAGCTAACAATCCCCGCCCGCGCCTCCAGCAACAGTTGATCGAACTTTTCCCGCTCGGGGGAAACCAGACCGCa                    >  2:865893/1‑93 (MQ=255)                                                                      ccGGTAAGCGGTTTGTTCCAGCTAACAATCCCCGCCCGCGCCTCCAGCAACAGTTGATCGAACTTTTCCCGCTCGGGGGAAACCAGACCGCaa                   <  1:3435835/93‑1 (MQ=255)                                                                        ggTAAGCGTTTTGTTCCAGCTAACAATCCCCGCCCGCGCCTCCAGCAACAGTTGATCGAACTTTTCCCGCTCGGGGGAAACCAGACCGCAaga                 <  2:2853252/93‑1 (MQ=255)                                                                        ggTAAGCGGTTTGTTCCAGCTAACAATCCCCGCCCGCGCCTCCAGCAACAGTTGATCGAACTTTTCCCGCTCGGGGGAAACCAGACCGCAaga                 >  1:1674523/1‑93 (MQ=255)                                                                        ggTAAGCGGTTTGTTCCAGCTAACAATCCCCGCCCGCGCCTCCAGCAACAGTTGATCGAACTTTTCCCGCTCGGGGGAAACCAGACCGCAaga                 >  2:70077/1‑93 (MQ=255)                                                                        ggTAAGCGGTTTGTTCCAGCTAACAATCCCCGCCCGCGCCTCCAGCAACAGTTGATCGAACTTTTCCCGCTCGGGGGAAACCAGACCGCAaga                 <  2:3488353/93‑1 (MQ=255)                                                                        ggTAAGCGGTTTGTTCCAGCTAACAATCCCCGCCCGCGCCTCCAGCAACAGTTGATCGAACTTTTCCCGCTCGGGGGAAACCAGACCGCAaga                 >  2:3910385/1‑93 (MQ=255)                                                                        ggTAAGCGGTTTGTCCCAGCTAGCAATCCCCGCCCGCGCCGCCAGCAACAGTTGATCGAACTTTTCCCGCTCGGGGGAAACCAGACCGCAaga                 <  1:2629571/93‑1 (MQ=255)                                                                         gTAAGCGGTTTGTTCCAGCTAACAATCCCCGCCCGCGCCTCCAGCAACAGTTGATCGAACTTTTCCCGCTCGGGGGAAACCAGACCGCAagag                <  1:2979833/93‑1 (MQ=255)                                                                          tAAGCTGTTCATTACAGCTACCAAGCCCAACCCGCGCCTCCAGCAGCAGTTGATCGAACTTTTCCCGCTCGGGGGAAACCAGACCGCAagaga               <  2:3092997/93‑1 (MQ=255)                                                                          tAAGCGGTTTGTTCCAGCTAACAATCCCCGCCCGCGCCTCCAGCAACAGTTGATCGAACTTTTCCCGCTCGGGGGAAACCAGACCGCAagaga               <  1:2224235/93‑1 (MQ=255)                                                                          tAAGCGGTTTGTTCCAGCTAACAATCCCCGCCCGCGCCTCCAGCAACAGTTGATCGAAATTTTCCCGCTCGCGGGAAAACAGACCGCAagaga               >  2:1887060/1‑93 (MQ=255)                                                                           aaGCGGTTTGTTCCAGCTAACAATCCCCGCCCGCGCCTCCAGCAACAGTTGATCGAACTTTTCCCGCTCGGGGGAAACCAGACCGCAagagag              >  1:817578/1‑93 (MQ=255)                                                                            aGCGGTTTGTTCCAGCTAACAATCCCCGCCCGCGCCTCCAGCGACAGTTGATCGAACTTTTCCCGCTCGGGGGAAACCAGACCGCAAGAGAGt             >  1:2515519/1‑93 (MQ=255)                                                                            aGCGGTTTGTTCCAGCTAACAATCCCCGCCCGCGCCTCCAGCAACAGTTGATCGAACTTTTCCCGCTCGGGGGAAACCAGACCGCAAGAGAGt             <  1:2054395/93‑1 (MQ=255)                                                                            aGCGGTTTGTTCCAGCTAACAATCCCCGCCCGCGCCTCCAGCAACAGTTGATCGAACTTTTCCCGCTCGGGGGAAACCAGACCGCAAGAGAGt             >  2:2585207/1‑93 (MQ=255)                                                                            aGCGGTTTGTTCCAGCTAACAATCCCCGCCCGCGCCTCCAGCAACAGTTGATCGAACTTTTCCCGCTCGGGGGAAACCAGACCGCAAGAGAGt             <  1:272045/93‑1 (MQ=255)                                                                               ggTTTGTTCCAGCTAACAATCCCCGCCCGCGCCTCCAGCAACAGTTGATCGAACTTTTCCCGCTCGGGGGAAACCAGACCGCAAGAGAGTCCg          >  1:2773474/1‑93 (MQ=255)                                                                                gTTTGTTCCAGCTAACAATCCCCGCCCGCGCCTCCAGCATCAGTTGATCGAACTTTTCCCGCTCGGGGGAAACCAGACCGCAAGAGAGTCCGa         <  1:843040/93‑1 (MQ=255)                                                                                gTTTGTTCCAGCTAACAATCCCCGCCCGCGCCTCCAGCAACAGTTGATCGAACTTTTCCCGCTCGGGGGAAACCAGACCGCAAGAGAGTCCGa         >  2:4571048/1‑93 (MQ=255)                                                                                gTTTGTCCCAGCTAACAATCCCCGCCCGCGCCTCCAGCAACAGTTGATCGAACTTTTCCCGCTCGGGGGAAACCAGACCGCAAGAGAGTCCGa         <  1:555915/93‑1 (MQ=255)                                                                                 tttGTTCCAGCTAACAATCCCCGCCCGCGCCTCCAGCAACAGTTGATCGAACTTTTCCCGCTCGGGGGAAACCAGACCGCAAGAGAGTCCGaa        <  2:3751327/93‑1 (MQ=255)                                                                                 tttGTTCCAGCTAACAATCCCCGCCCGCGCCTCCAGCAACAGTTGATCGAACTTTTCCCGCTCGGGGGAAACCAGACCGCAAGAGAGTCCGaa        <  1:456618/93‑1 (MQ=255)                                                                                 tttGTTCCAGCTAACAATCCCCGCCCGCGCCTCCAGCAACAGTTGATCGAACTTTTCCCGCTCGGGGGAAACCAGACCGCAAGAGAGTCCGaa        <  1:4610268/93‑1 (MQ=255)                                                                                   tGTTCCAGCTAACAATCCCCGCCCGCGCCTCCAGCAACAGTTGATCGAACTTTTCCCGCTCGGGGGAAACCAGACCGCAAGAGAGTCCGAAGc      >  1:3581719/1‑93 (MQ=255)                                                                                     ttacAGCTAACAATCCCCGCCCGCGCCTCCAGCAACAGTTGATCGAACTTTTCCCGCTCGGGGGAAACCAGACCGCAAGAGAGTCCGAAGCGa    >  1:4111393/4‑93 (MQ=255)                                                                                       ccAGCTAACAATCCCCGCCCGCGCCTCCAGCAACAGTTGATCGAACTTTTCCCGCTCGGGGGAAACCAGACCGCAAGAGAGTCCGAAGCGAgt  <  1:3101331/93‑1 (MQ=255)                                                                                           |                                                                                         GCATACCAGGCGCTGGGGCGATGGTTACCGGAAGCACCAATGCCGCCGAATGGCGCGGTACTGGCAGCACCGGTAAGCGGTTTGTTCCAGTTAACAATCCCCGCCCGCGCCTCCAGCAACAGTTGATCGAACTTTTCCCGCTCGGGGGAAACCAGACCGCAAGAGAGTCCGAAGCGAGT  >  NC_000913/1826386‑1826564

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