Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,547,765	G→T	E403* (GAG→TAG)	gntT →	high‑affinity gluconate transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,547,765	0	G	T	100.0%	159.5 / NA	44	E403* (GAG→TAG)	gntT	high‑affinity gluconate transporter
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base T (30/14);  total (30/14)

GCTGATGGTTATTGCGGTTGGTTCCGGTAGTGTGATTTTCTCTCATGTGAACGATCCGGGCTTCTGGCTGTTCAAAGAGTACTTTAACCTGACTATCGGCGAGACCATCAAATCCTGGTCGATGCTGGAAACCATTATCTCGGTGTGCGGTCT  >  NC_000913/3547689‑3547841                                                                             |                                                                             gCTGATGGTTATTGCGGTTGGTTCCGGTAGTGTGATTTTCTCTCATGTGAACGATCCGGGCTTCTGGCTGTTCAAATAGTa                                                                          >  2:680589/1‑81 (MQ=255) gCTGATGGTTATTGCGGTTGGTTCCGGTAGTGTGATTTTCTCTCATGTGAACGATCCGGGCTTCTGGCTGTTCAAATAGTa                                                                          >  2:680601/1‑81 (MQ=255)  cTGATGGTTATTGCGGTTGGTTCCGGTAGTGTGATTTTCTCTCATGTGAACGATCCGGGCTTCTGGCTGTTCAAATAGTAc                                                                         >  2:680587/1‑81 (MQ=255)   tGATGGTTATTGCGGTTGGTTCCGGTAGTGTGATTTTCTCTCATGTGAACGATCCGGGCTTCTGGCTGTTCAAATAGTACt                                                                        >  1:680583/1‑81 (MQ=255)    gATGGTTATTGCGGTTGGTTCCGGTAGTGTGATTTTCTCTCATGTGAACGATCCGGGCTTCTGGCTGTTCAAATAGTACtt                                                                       >  1:680614/1‑81 (MQ=255)     aTGGTTATTGCGGTTGGTTCCGGTAGTGTGATTTTCTCTCATGTGAACGATCCGGGCTTCTGGCTGTTCAAATAGTACtt                                                                       >  1:680577/1‑80 (MQ=255)        gTTATTGCGGTTGGTTCCGGTAGTGTGATTTTCTCTCATGTGAACGATCCGGGCTTCTGGCTGTTCAAATAGTACTTTaa                                                                    <  1:680536/80‑1 (MQ=255)         ttATTGCGGTTGGTTCCGGTAGTGTGATTTTCTCTCATGTGAACGATCCGGGCTTCTGGCTGTTCAAATAGTACTTTAAcc                                                                  <  1:680537/81‑1 (MQ=255)         ttATTGCGGTTGGTTCCGGTAGTGTGATTTTCTCTCATGTGAACGATCCGGGCTTCTGGCTGTTCAAATAGTACTTTAAcc                                                                  >  1:680610/1‑81 (MQ=255)          tATTGCGGTTGGTTCCGGTAGTGTGATTTTCTCTCATGTGAACGATCCGGGCTTCTGGCTGTTCAAATAGTACTTTAACCt                                                                 >  2:680609/1‑81 (MQ=255)            ttGCGGTTGGTTCCGGTAGTGTGATTTTCTCTCATGTGAACGATCCGGGCTTCTGGCTGTTCAAATAGTACTTTAACCTGa                                                               >  2:680564/1‑81 (MQ=255)              gCGGTTGGTTCCGGTAGTGTGATTTTCTCTCATGTGAACGATCCGGGCTTCTGGCTGTTCAAATAGTACTTTAACCTGAc                                                              >  1:680538/1‑80 (MQ=255)               cGGTTGGTTCCGGTAGTGTGATTTTCTCTCATGTGAACGATCCGGGCTTCTGGCTGTTCAAATAGTACTTTAACCTGACt                                                             <  1:680539/80‑1 (MQ=255)               cGGTTGGTTCCGGTAGTGTGATTTTCTCTCATGTGAACGATCCGGGCTTCTGGCTGTTCAAATAGTACTTTAACCTGACTa                                                            >  2:680620/1‑81 (MQ=255)                ggttggttCCGGTAGTGTGATTTTCTCTCATGTGAACGATCCGGGCTTCTGGCTGTTCAAATAGTACTTTAACCTGACTAt                                                           >  1:680613/1‑81 (MQ=255)                    ggttCCGGTAGTGTGATTTTCTCTCATGTGAACGATCCGGGCTTCTGGCTGTTCAAATAGTACTTTAACCTGACTATCGGc                                                       >  2:680584/1‑81 (MQ=255)                    ggttCCGGTAGTGTGATTTTCTCTCATGTGAACGATCCGGGCTTCTGGCTGTTCAAATAGTACTTTAACCTGACTATCGGc                                                       >  2:680582/1‑81 (MQ=255)                          ggTAGTGTGATTTTCTCTCATGTGAACGATCCGGGCTTCTGGCTGTTCAAATAGTACTTTAACCTGACTATCGGCGAGAc                                                  >  1:680594/1‑80 (MQ=255)                           gTAGTGTGATTTTCTCTCATGTGAACGATCCGGGCTTCTGGCTGTTCAAATAGTACTTTAACCTGACTATCGGCGAGACCa                                                <  2:680540/81‑1 (MQ=255)                              gtgtGATTTTCTCTCATGTGAACGATCCGGGCTTCTGGCTGTTCAAATAGTACTTTAACCTGACTATCGGCGAGACCATc                                              >  2:680604/1‑80 (MQ=255)                                gtgATTTTATCTAATGTGAACGATCCGGGCTTCTGGCTGTTCAAATAGTACTTTAACCTGACTATCGGCGAGACCATCaaa                                           >  2:680588/1‑81 (MQ=255)                                     tttCTCTCATGTGAACGATCCGGGCTTCTGGCTGTTCAAATAGTACTTTAACCTGACTATCGGCGAGACCATCAAATCCTg                                      >  1:680585/1‑81 (MQ=255)                                       tctctcATGTGAACGATCCGGGCTTCTGGCTGTTCAAATAGTACTTTAACCTGACTATCGGCGAGACCATCAAATCCTGGt                                    >  1:680595/1‑81 (MQ=255)                                          ctcATGTGAACGATCCGGGCTTCTGGCTGTTCAAATAGTACTTTAACCTGACTATCGGCGAGACCATCAAATCCTGGTCGa                                 >  2:680599/1‑81 (MQ=255)                                          ctcATGTGAACGATCCGGGCTTCTGGCTGTTCAAATAGTACTTTAACCTGACTATCGGCGAGACCATCAAATCCTGGTCGa                                 >  2:680631/1‑81 (MQ=255)                                           tcATGTGAACGATCCGGGCTTCTGGCTGTTCAAATAGTACTTTAACCTGACTATCGGCGAGACCATCAAATCCTGGTCGAt                                >  2:680622/1‑81 (MQ=255)                                             aTGTGAACGATCCGGGCTTCTGGCTGTTCAAATAGTACTTTAACCTGACTATCGGCGAGACCATCAAATCCTGGTCGATGc                              >  2:680611/1‑81 (MQ=255)                                              tgtgAACGATCCGGGCTTCTGGCTGTTCAAATAGTACTTTAACCTGACTATCGGCGAGACCATCAAATCCTGGTCGATGCt                             <  1:680541/81‑1 (MQ=255)                                               gtgAACGATCCGGGCTTCTGGCTGTTCAAATAGTACTTTAACCTGACTATCGGCGAGACCATCAAATCCTGGTCGATGCTg                            >  2:680615/1‑81 (MQ=255)                                                tgAACGATCCGGGCTTCTGGCTGTTCAAATAGTACTTTAACCTGACTATCGGCGAGACCATCAAATCCTGGTCGATGCTgg                           <  1:680543/81‑1 (MQ=255)                                                tgAACGATCCGGGCTTCTGGCTGTTCAAATAGTACTTTAACCTGACTATCGGCGAGACCATCAAATCCTGGTCGATGCTg                            <  2:680542/80‑1 (MQ=255)                                                 gAACGATCCGGGCTTCTGGCTGTTCAAATAGTACTTTAACCTGACTATCGGCGAGACCATCAAATCCTGGTCGATGCTgg                           >  1:680621/1‑80 (MQ=255)                                                 gAACGATCCGGGCTTCTGGCTGTTCAAATAGTACTTTAACCTGACTATCGGCGAGACCATCAAATCCTGGTCGATGCTgg                           >  1:680597/1‑80 (MQ=255)                                                   aCGATCCGGGCTTCTGGCTGTTCAAATAGTACTTTAACCTGACTATCGGCGAGACCATCAAATCCTGGTCGATGCTGGaaa                        >  2:680602/1‑81 (MQ=255)                                                      aTCCGGGCTTCTGGCTGTTCAAATAGTACTTTAACCTGACTATCGGCGAGACCATCAAATCCTGGTCGATGCTGGAAACCa                     <  1:680544/81‑1 (MQ=255)                                                       tCCGGGCTTCTGGCTGTTCAAATAGTACTTTAACCTGACTATCGGCGAGACCATCAAATCCTGGTCGATGCTGGAAACCAt                    >  2:680623/1‑81 (MQ=255)                                                         cGGGCTTCTGGCTGTTCAAATAGTACTTTAACCTGACTATCGGCGAGACCATCAAATCCTGGTCGATGCTGGAAACCAtt                   <  2:680545/80‑1 (MQ=255)                                                          gggCTTCTGGCTGTTCAAATAGTACTTTAACCTGACTATCGGCGAGACCATCAAATCCTGGTCGATGCTGGAAACCATTAt                 <  1:680546/81‑1 (MQ=255)                                                           ggCTTCTGGCTGTTCAAATAGTACTTTAACCTGACTATCGGCGAGACCATCAAATCCTGGTCGATGCTGGAAACCATTAtc                >  2:680607/1‑81 (MQ=255)                                                                cTGGCTGTTCAAATAGTACTTTAACCTGACTATCGGCGAGACCATCAAATCCTGGTCGATGCTGGAAACCATTATCTCGgt           >  2:680612/1‑81 (MQ=255)                                                                   gCTGTTCAAATAGTACTTTAACCTGACTGTCGGCGAGACCATCAAATCCTGGTCGATGCTGGAAACCATTATCTCGgt           <  2:680547/78‑1 (MQ=255)                                                                   gCTGTTCAAATAGTACTTTAACCTGACTATCGGCGAGACCATCAAATCCTGGTCGATGCTGGAAACCATTATCTCGgt           <  1:680548/78‑1 (MQ=255)                                                                          aaaTAGTACTTTAACCTGACTATCGGCGAGACCGTCAAATCCTGGTCGATGCTGGAAACCATTATCTCGGTGTGCGGTCt  <  1:680550/80‑1 (MQ=255)                                                                          aaaTAGTACTTTAACCTGACTATCGGCGAGACCATCAAATCCTGGTCGATGCTGGAAACCATTATCTCGGTGTGCGGTCt  <  2:680549/80‑1 (MQ=255)                                                                             |                                                                             GCTGATGGTTATTGCGGTTGGTTCCGGTAGTGTGATTTTCTCTCATGTGAACGATCCGGGCTTCTGGCTGTTCAAAGAGTACTTTAACCTGACTATCGGCGAGACCATCAAATCCTGGTCGATGCTGGAAACCATTATCTCGGTGTGCGGTCT  >  NC_000913/3547689‑3547841

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