Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,612,018	C→T	G245D (GGT→GAT)	b1525 ←	putative aldehyde dehydrogenase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,612,018	0	C	T	100.0%	269.0 / NA	74	G245D (GGT→GAT)	b1525	putative aldehyde dehydrogenase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (47/27);  total (47/27)

CTGTCCGGTATTCTGATAACGTCCGGCTACCGCCGCTTTCACCGCCAGTTCCAGATCGGCATCGTTAAGCACAATAAACGGATCCGAACCGCCCAGTTCCAGTACGCATTTTTTCAGTGCCGCTCCAGCCTGTGCGCCAATAGCCGCTCCCGCACGAACACTTCCGGTCACCGTGACA  >  NC_000913/1611930‑1612107                                                                                         |                                                                                          cTGTCCGGTATTCTGATAACGTCCGGCTACCGCCGCCTTCACCGCCCGTTCCAGATCGGCATCGTTAAGAACAATACACGGATCCGAATCGcc                                                                                       >  2:1944124/1‑93 (MQ=255)   gTCCGGTATTCTGATAACGTCCGGCTACCGCCGCTTTCACCGCCAGTTCCAGATCGGCATCGTTAAGCACAATAAACGGATCCGAATCGCCCa                                                                                     >  2:301232/1‑93 (MQ=255)   gTCCGGTATTCTGATAACGTCCGGCTACCGCCGCTTTCACCGCCAGTTCCAGATCGGCATCGTTAAGCACAATAAACGGATCCGAATCGCCCa                                                                                     >  2:301233/1‑93 (MQ=255)   gTCCGGTATTCTGATAACGTCCGGCTACCGCCGCTTTCACCGCCAGTTCCAGATCGGCATCGTTAAGCACAATAAACGGATCCGAATCGCCCa                                                                                     >  1:211267/1‑93 (MQ=255)   gTCCGGTATTCTGATAACGTCCGGCTACCGCCGCTTTCACCGCCAGTTCCAGATCGGCATCGTTAAGCACAATAAACGGATCCGAATCGCCCa                                                                                     <  1:1593255/93‑1 (MQ=255)        gTATTCTGATAACGTCCGGCTACCGCCGCTTTCACCGCCAGTTCCAGATCGGCATCGTTAAGCACAATAAACGGATCCGAATCGCCCAGTTcc                                                                                >  1:1175296/1‑93 (MQ=255)         tATTCTGATAACGTCCGGCTACCGCCGCTTTCACCGCCAGTTCCAGATCGGCATCGTTAAGCACAATAAACGGATCCGAATCGCCCAGTTCCa                                                                               <  2:3022624/93‑1 (MQ=255)           ttCTGATAACGTCCGGCTACCGCCGCTTTCACCGCCAGTTCCAGATCGGCATCGTTAAGCACAATAAACGGATCCGAATCGCCCAGTTCCAGt                                                                             >  1:3957688/1‑93 (MQ=255)             cTGATAACGTCCGGCTACCGCCGCTTTCACCGCCAGTTCCAGATCGGCATCGTTAAGCACAATAAACGGATCCGAATCGCCCAGTTCCAGTAc                                                                           <  1:1177271/93‑1 (MQ=255)             cTGAAAACGTCCGGACCCCGCCACTTTCACCGCCAGGTCCAGATCGGCATCGTTAAGCACAATAAACGGATCCGAATCGCCCAGTTCCAGTAc                                                                           >  2:3222286/1‑93 (MQ=255)              tGATAACGTCCGGCTACCGCCGCTTTCACCGCCAGTTCCAGATCGGCATCGTTAAGCACAATAAACGGATCCGAATCGCCCAGTTCCAGTACg                                                                          >  2:529824/1‑93 (MQ=255)              tGATAACGTCCGGCTACCGCCGCTTTCACCGCCAGTTCCAGATCGGCATCGTTAAGCACAATAAACGGATCCGAATCGCCCAGTTCCAGTACg                                                                          >  1:3036625/1‑93 (MQ=255)               gATAACGTCCGGCTACCGCCGCTTTCACCGCCAGTTCCAGATCGGCATCGTTAAGCACAATAAACGGATCCGAATCGCCCCGTTCCAGTACGc                                                                         >  1:4324380/1‑93 (MQ=255)               gATAACGTCCGGCTACCGCCGCTTTCACCGCCAGTTCCAGATCGGCATCGTTAAGCACAATAAACGGATCCGAATCGCCCAGTTCCAGTACGc                                                                         <  2:3473520/93‑1 (MQ=255)               gATAACGTCCGGCTACCGCCGCTTTCACCGCCAGTTCCAGATCGGCATCGTTAAGCACAATAAACGGATCCGAATCGCCCAGTTCCAGTACGc                                                                         >  1:2734258/1‑93 (MQ=255)                aTAACGTCCGGCTACCGCCGCTTTCACCTCCAGTTCCAGATCGGCATCGTTAAGCACAATAAACGGATCCGAATCGCCCAGTTCCAGTACGCa                                                                        >  1:263541/1‑93 (MQ=255)                 tAACGTCCGGCTACCGCCGCTTTTACCGCCAGTTCCAGATCGGCATCGATAAGCACAATAAACGGATTCGAATTGCCCAGTTACAGTACGCAt                                                                       >  2:2923891/1‑93 (MQ=255)                 tAACGTCCGGCTACCGCCGCTTTCACCGCCAGTTCCAGATCGGCATCGTTAAGCACAATAAACGGATCCGAATCGCCCAGTTCCAGTACGCAt                                                                       >  1:2168609/1‑93 (MQ=255)                   aCGTCCGGCTACCGCCGCTTTCACCGCCAGTTCCAGATCGGCATCGTTAAGCACAATAAACGGATCCGAATCGCCCAGTTCCAGTACGCAttt                                                                     >  1:461528/1‑93 (MQ=255)                    cgcccGGCTACCGCCGCTTTCACCGCCAGTTCCAGATCGGCATCGTTAAGCACAATAAACGGATCCGAATCGCCCAGTTCCAGTACGCAtttt                                                                    <  1:3567563/90‑1 (MQ=255)                      tCCGGCTACCGCCGCTTTCACCGCCAGTTCCAGATCGGCATCGTTAAGCACAATAAACGGATCCGAATCGCCCAGTTCCAGTACGCAtttttt                                                                  <  1:814668/93‑1 (MQ=255)                        cGGCTACCGCCGCTTTCACCGCCAGTTCCAGATCGGCATCGTTAAGCACAATAAACGGATCCGAATCGCCCAGTTCCAGTACGCATTTTTTCa                                                                >  2:3986439/1‑93 (MQ=255)                         ggCTACCGCCGCTTTCACCGCCAGTTCCAGATCGGCATCGTTAAGCACAATAAACGGATCCGAATCGCCCAGTTCCAGTACGCATTTTTTCAg                                                               <  1:693201/93‑1 (MQ=255)                          gCTACCGCCGCTTTCACCGCCAGTTCCAGATCGGCATCGTTAAGCACAATAAACGGATCCGAATCGCCCAGTTCCAGTACGCATTTTTTCAGt                                                              <  1:993693/93‑1 (MQ=255)                           cTACCGCCGCTTTCACCGCCAGTTCCAGATCGGCATCGTTAAGCACAATAAACGGATCCGAATCGCCCAGTTCCAGTGCGCATTTTTTCAGTg                                                             <  1:882292/93‑1 (MQ=255)                             aCCGCCGCTTTCACCGCCAGTTCCAGATCGGCATCGTTAAGCACAATAAACGGATCCGAATCGCCCAGTTCCAGTACGCATTTTTTCAGTGcc                                                           <  1:765066/93‑1 (MQ=255)                             aCCGCCGCTTTCACCGCCAGTTCCAGATCGGCATCGTTAAGCACAATAAACGGATCCGAATCGCCCAGTTCCAGTACGCATTTTTTCAGTGcc                                                           <  2:4430565/93‑1 (MQ=255)                               cgccgcTTTCACCGCCAGTTCCAGATCGGCATCGTTAAGCACAATAAACGGATCCGAATCGCCCAGTTCCAGTACGCATTTTTTCAGTGCCGc                                                         <  1:4137490/93‑1 (MQ=255)                                gccgcTTTCACCGCCAGTTCCAGATCGGCATCGTTAAGCACAATAAACGGATCCGAATCGCCCAGTTCCAGTACGCATTTTTTCAGTGCCGCt                                                        >  1:3575248/1‑93 (MQ=255)                                  cgcTTTCACCGCCAGTTCCAGATCGGCATCGTTAAGCACAATAAACGGATCCGAATCGCCCAGTTCCAGTACGCATTTTTTCAGTGCCGCTcc                                                      >  2:4312698/1‑93 (MQ=255)                                    cTTTCACCGCCAGTTCCAGATCGGCATCGTTAAGCACAATAAACGGATCCGAATCGCCCAGTTCCAGTACGCATTTTTTCAGTGCCGCTCCAg                                                    >  1:1518852/1‑93 (MQ=255)                                    cTTTCACCGCCAGTTCCAGATCGGCATCGTTAAGCACAATAAACGGATCCGAATCGCCCAGTTCCAGTACGCATTTTTTCAGTGCCGCTCCAg                                                    >  1:767257/1‑93 (MQ=255)                                      ttCACCGCCAGTTCCAGATCGGCATCGTTAAGCACAATAAACGGATCCGAATCGCCCAGTTCCAGTACGCATTTTTTCAGTGCCGCTCCAGcc                                                  <  1:1254671/93‑1 (MQ=255)                                      ttCACCGCCAGTTCCAGATCGGCATCGTTAAGCACAATAAACGGATCCGAATCGCCCAGTTCCAGTACGCATTTTTTCAGTGCCGCTCCAGcc                                                  >  1:1992600/1‑93 (MQ=255)                                       tCACCGCCAGTTCCAGATCGGGATCGTTAAGATCAATAAACGGAGCCGAATCGCCCAGTTCCAGGACGCATTTTTTCAGTGCCGCTCCAGCCt                                                 >  2:4228694/1‑93 (MQ=255)                                            gCCAGTTCCAGATGGGCATCGTTAAGCACAATAAACGGATCCGAATCGCCCAGTTCCAGTACGCATTTTTTCAGTGCCGCTCCAGCCTGTgcg                                            >  2:4605938/1‑93 (MQ=255)                                                gTTCCAGATCGGCATCGTTAAGCACAATAAACGGATCCGAATCGCCCAGTTCCAGTACGCATTTTTTCAGTGCCGCTCCAGCCTGTGCGCCaa                                        <  1:468199/93‑1 (MQ=255)                                                gTTCCAGATCGGCATCGTTAAGCACAATAAACGGATCCGAATCGCCCAGTTCCAGTACGCATTTTTTCAGTGCCGCTCCAGCCTGTGCGCCaa                                        >  2:4267614/1‑93 (MQ=255)                                                gTTCCAGATCGGCATCGTTAAGCACAATAAACGGATCCGAATCGCCCAGTTCCAGTACGCATTTTTTCAGTGCCGCTCCAGCCTGTGCGCCaa                                        >  2:202304/1‑93 (MQ=255)                                                gTTCCAGATCGGCATCGTTAAGCACAATAAACGGATCCGAATCGCACAGTTCCAGTACGCATTTTTTCAGTGCCGCTCCAGCCTGTGCGCCaa                                        >  2:457500/1‑93 (MQ=255)                                                 ttCCAGATCGGCATCGTTAAGCACAATAAACGGATCCGAATCGCCCAGTTCCAGTACGCATTTTTTCAGTGCCGCTCCAGCCTGTGCGCCAAt                                       >  1:2527531/1‑93 (MQ=255)                                                 ttCCAGATCGGCATCGTTAAGCACAATAAACGGATCCGAATCGCCCAGTTCCAGTACGCATTTTTTCAGTGCCGCTCCAGCCTGTGCGCCAAt                                       >  1:1293877/1‑93 (MQ=255)                                                   ccAGATCGGCATCGTTAAGCACAATAAACGGATCCGAATCGCCCAGTTCCAGTACGCATTTTTTCAGTGCCGCTCCAGCCTGTGCGCCAATAg                                     <  2:335043/93‑1 (MQ=255)                                                    cAGATCGGCATCGTTAAGCACAATAAACGGATCCGAATCGCCCAGTTCCAGTACGCATTTTTTCAGTGCCGCTCCAGCCTGTGCGCCAATAGc                                    >  2:751324/1‑93 (MQ=255)                                                      gATCGGCATCGTTAAGCACAATAAACGGATCCGAATCGCCCAGTTCCAGTACGCATTTTTTCAGTGCCGCTCCAGCCTGTGCGCCAATAGCCg                                  >  2:1402805/1‑93 (MQ=255)                                                       aTCGGCATCGTTAAGCACAATAAACGGATCCGAATCGCCCAGTTCCAGTACGCATTTTTTCAGTGCCGCTCCAGCCTGTGCGCCAATAGCCGc                                 <  2:3013400/93‑1 (MQ=255)                                                       aTCGGCATCGTTAAGCACAATAAACGGATCCGAATCGCCCAGTTCCAGTACGCATTTTTTCAGTGCCGCTCCAGCCTGTGCGCCAATAGCCGc                                 >  2:3047729/1‑93 (MQ=255)                                                         cGGCATCGTTAATCACAATAAACGGATCCGAATCGCCCAGTTCCAGTACGCATTTTTTCAGTGCCGCTCCAGCCTGTGCGCCAATAGCCGCTc                               <  2:3276697/93‑1 (MQ=255)                                                         cGGCATCGTTAAGCACAATAAACGGATCCGAATCGCCCAGTTCCAGTACGCATTTTTTCAGTGCCGCTCCAGCCTGTGCGCCAATAGCCGCTc                               >  2:4167583/1‑93 (MQ=255)                                                           gCATCGTTAAGCACAATAAACGGATCCGAATCGCCCAGTTCCAGTACGCATTTTTTCAGTGCCGCTCCAGCCTGTGCGCCAATAGCCGCTccc                             >  2:1948567/1‑93 (MQ=255)                                                             aTCGTTAAGCACAATAAACGGATCCGAATCGCCCAGTTCCAGTACGCATTTTTTCAGTGCCGCTCCAGCCTGTGCGCCAATAGCCGCTCCCGc                           >  1:1244234/1‑93 (MQ=255)                                                             aTCGTTAAGCACAATAAACGGATCCGAATCGCCCAGTTCCAGTACGCATTTTTTCAGTGCCGCTCCAGCCTGTGCGCCAATAGCCGCTCCCGc                           >  2:2607535/1‑93 (MQ=255)                                                              tCGTTAAGCACAATAAACGGATCCGAATCGCCCAGTTCCAGTACGCATTTTTTCAGTGCCGCTCCAGCCTGTGCGCCAATAGCCGCTCCCGCa                          >  2:4433330/1‑93 (MQ=255)                                                                gTTAAGCACAATAAACGGATCCGAATCGCCCAGTTCCAGTACGCATTTTTTCAGTGCCGCTCCAGCCTGTGCGCCAATAGCCGCTCCCGCACg                        >  2:2108934/1‑93 (MQ=255)                                                                gTTAAGCACAATAAACGGATCCGAATCGCCCAGTTCCAGTACGCATTTTTTCAGTGCCGCTCCAGCCTGTGCGCCAATAGCCGCTCCCGCACg                        >  2:2489038/1‑93 (MQ=255)                                                                 ttAAGCACAATAAACGGATCCGAATCGCCCAGTTCCAGTACGCATTTTTTCAGTGCCGCTCCAGCCTGTGCGCCAATAGCCGCTCCCGCACGa                       <  1:1755837/93‑1 (MQ=255)                                                                  tAAGCACAATAAACGGATCCGAATCGCCCAGTTCCAGTACGCATTTTTTCAGTGCCGCTCCAGCCTGTGCGCCAATAGCCGCTCCCGCACGaa                      >  2:1317940/1‑93 (MQ=255)                                                                  tAAGCACAATAAACGGATCCGAATCGCCCAGTTCCAGTACGCATTTTTTCAGTGCCGCTCCAGCCTGTGCGCCAATAGCCGCTCCCGCACGaa                      <  1:1768945/93‑1 (MQ=255)                                                                  tAAGCACAATAAACGGACCCGAATCTTCCAGTTCCAGTACGCATTTTTTCAGTGCCGCCCCAGCCTGTGCGCCAATAGCCGCTCCCGCACGaa                      <  1:3876216/93‑1 (MQ=255)                                                                   aaGCACAATAAACGGATCCGAATCGCCCAGTTCCAGTACGCATTTTTTCAGTGCCGCTCCAGCCTGTGCGCCAATAGCCGCTCCCGCACGAac                     <  2:1451080/93‑1 (MQ=255)                                                                     gCACAATAAACGGATCCGAATCGCCCAGTTCCAGTACGCATTTTTTCAGTGCCGCTCCAGCCTGTGCGCCAATAGCCGCTCCCGCACGAacac                   >  1:1134338/1‑93 (MQ=255)                                                                       acaATAAACGGATCCGAATCGCCCAGTTCCAGTACGCATTTTTTCAGTGCCGCTCCAGCCTGTGCGCCAATAGCCGCTCCCGCACGAACACtt                 <  2:3158344/93‑1 (MQ=255)                                                                       acaATAAACGGATCCGAATCGCCCAGTTCCAGTACGCATTTTTTCAGTGCCGCTCCAGCCTGTGCGCCAATAGCCGCTCCCGCACGAACACtt                 >  1:3453121/1‑93 (MQ=255)                                                                        caATAAACGGATCCGAATCGCCCAGTTCCAGTACGCATTTTTTCAGTGCCGCTCCAGCCTGTGCGCCAATAGCCGCTCCCGCACGAACACTTc                <  2:572849/93‑1 (MQ=255)                                                                         aaTAAACGGATCCGAATCGCCCAGTTCCAGTACGCATTTTTTCAGTGCCGCTCCAGCCTGTGCGCCAATAGCCGCTCCCGCACGAACACTTcc               >  2:1883991/1‑93 (MQ=255)                                                                         aaTAAACGGATCCGAATCGCCCAGTTCCAGTACGCATTTTTTCAGTGCCGCTCCAGCCTGTGCGCCAATAGCCGCTCCCGCACGAACACTTcc               <  1:3491784/93‑1 (MQ=255)                                                                          aTAAACGGATCCGAATCGCCCAGTTCCAGTACGCATTTTTTCAGTGCCGCTCCAGCCTGTGCGCCAATAGCCGCTCCCGCACGAACACTTCCg              >  1:3240121/1‑93 (MQ=255)                                                                          aTAAACGGATCCGAATCGCCCAGTTCCAGTACGCATTTTTTCAGTGCCGCTCCAGCCTGTGCGCCAATAGCCGCTCCCGCACGAACACTTCCg              >  1:1517567/1‑93 (MQ=255)                                                                            aaaCGGATCCGAATCGCCCAGTTCCAGTACGCATTTTTTCAGTGCCGCTCCAGCCTGTGCGCCAATAGCCGCTCCCGCACCAACACTTCCGGt            >  1:4076188/1‑93 (MQ=255)                                                                             aaCGGATCCGAATCGCCCAGTTCCAGTACGCATTTTTTCAGTGCCGCTCCAGCCTGTGCGCCAATAGCCGCTCCCGCACGAACACTTCCGGTc           >  1:2944772/1‑93 (MQ=255)                                                                              aCGGATCCGAATCGCCCAGTTCCAGTACGCATTTTTTCAGTGCCGCTCCAGCCTGTGCGCCAATAGCCGCTCCCGCACGAACACTTCCGGTCa          >  2:3915006/1‑93 (MQ=255)                                                                                gtaTCCGAATCGCCCAGTTCCAGTACGCATTTTTTCAGTGCCGCTCCAGCCTGTGCGCCAATAGCCGCTCCCGCACGAACACTTCCGGTCAcc        <  1:4294211/91‑1 (MQ=255)                                                                                ggATCCGAATCGCCCAGTTCCAGTACGCATTTTTTCAGTGCCGCTCCAGCCTGTGCGCCAATAGCCGCTCCCGCACGAACACTTCCGGTCAcc        >  1:863318/1‑93 (MQ=255)                                                                                 gATCCGAATCGCCCAGTTCCAGTACGCATTTTTTCAGTGCCGCTCCAGCCTGTGCGCCAATAGCCGCTCCCGCACGAACACTTCCGGTCACCg       <  2:1055951/93‑1 (MQ=255)                                                                                  aTCCGAATCGCCCAGTTCCAGTACGCATTTTTTCAGTGCCGCTCCAGCCTGTGCGCCAATAGCCGCTCCCGCACGAACACTTCCGGTCACCGt      <  1:162311/93‑1 (MQ=255)                                                                                  aTCCGAATCGCCCAGTTCCAGTACGCATTTTTTCAGTGCCGCTCCAGCCTGTGCGCCAATAGCCGCTCCCGCACGAACACTTCCGGTCACCGt      >  2:4611729/1‑93 (MQ=255)                                                                                   tCCGAATCGCCCAGTTCCAGTACGCATTTTTTCAGTGCCGCTCCAGCCTGTGCGCCAATAGCCGCTCCCGCACGAACACTTCCGGTCACCGTg     >  1:3084881/1‑93 (MQ=255)                                                                                   tCCGAATCGCCCAGTTCCAGTACGCATTTTTTCAGTGCCGCTCCAGCCTGTGCGCCAATAGCCGCTCCCGCACGAACACTTCCGGTCACCGTg     <  1:302403/93‑1 (MQ=255)                                                                                      gAATCGCCCAGTTCCAGTACGCATTTTTTCAGTGCCGCTCCAGCCTGTGCGCCAATAGCCGCTCCCGCACGAACACTTCCGGTCACCGTGAca  >  2:3286361/1‑93 (MQ=255)                                                                                         |                                                                                          CTGTCCGGTATTCTGATAACGTCCGGCTACCGCCGCTTTCACCGCCAGTTCCAGATCGGCATCGTTAAGCACAATAAACGGATCCGAACCGCCCAGTTCCAGTACGCATTTTTTCAGTGCCGCTCCAGCCTGTGCGCCAATAGCCGCTCCCGCACGAACACTTCCGGTCACCGTGACA  >  NC_000913/1611930‑1612107

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