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%	148.3 / NA	43	E403* (GAG→TAG)	gntT	high‑affinity gluconate transporter
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base T (28/15);  total (28/15)

AGCTGATGGTTATTGCGGTTGGTTCCGGTAGTGTGATTTTCTCTCATGTGAACGATCCGGGCTTCTGGCTGTTCAAAGAGTACTTTAACCTGACTATCGGCGAGACCATCAAATCCTGGTCGATGCTGGAAACCATTATCTCGGTGTGCGGTCTG  >  NC_000913/3547688‑3547842                                                                              |                                                                              aGCTGATGGTTATTGCGGTTGGTTCCGGTAGTGTGATTTTCTCTCATGTGAACGATCCGGGCTTCTGGCTGTTCAAatag                                                                             >  1:630155/1‑80 (MQ=255)   cTGATGGTTATTGCGGTTGGTTCCGGTAGTGTGATTTTCTCTCATGTGAACGATCCGGGCTTCTGGCTGTTCAAATAGTAc                                                                          >  2:630158/1‑81 (MQ=255)    tGATGGTTATTGCGGTTGGTTCCGGTAGTGTGATTTTCTCTCATGTGAACGATCCGGGCTTCTGGCTGTTCAAATAGTACt                                                                         >  2:630153/1‑81 (MQ=255)     gATGGTTATTGCGGTTGGTTCCGGTAGTGTGATTTTCTCTCATGTGAACGATCCGGGCTTCTGGCTGTTCAAATAGTACtt                                                                        >  2:630142/1‑81 (MQ=255)            aTTGCGGTTGGTTCCGGTAGTGTGATTTTCTCTCATGTGAACGATCCGGGCTTCTGGCTGTTCAAATAGTACTTTAACCt                                                                  >  2:630157/1‑80 (MQ=255)             ttGCGGTTGGTTCCGGTAGTGTGATTTTCTCTCATGTGAACGATCCGGGCTTCTGGCTGTTCAAATAGTACTTTAACCTGa                                                                <  2:630110/81‑1 (MQ=255)               gCGGTTGGTTCCGGTAGTGTGATTTTCTCTCATGTGAACGATCCGGGCTTCTGGCTGTTCAAATAGTACTTTAACCTGAc                                                               >  1:630164/1‑80 (MQ=255)                cGGTTGGTTCCGGTAGTGTGATTTTCTCTCATGTGAACGATCCGGGCTTCTGGCTGTTCAAATAGTACTTTAACCTGACt                                                              <  1:630112/80‑1 (MQ=255)                cGGTTGGTTCCGGTAGTGTGATTTTCTCTCATGTGAACGATCCGGGCTTCTGGCTGTTCAAATAGTACTTTAACCTGACt                                                              <  1:630111/80‑1 (MQ=255)                 ggttggttCCGGTAGTGTGATTTTCTCTCATGTGAACGATCCGGGCTTCTGGCTGTTCAAATAGTACTTTAACCTGACTAt                                                            <  2:630113/81‑1 (MQ=255)                  gttggttCCGGTAGTGTGATTTTCTCTCATGTGAACGATCCGGGCTTCTGGCTGTTCAAATAGTACTTTAACCTGACTATc                                                           >  2:630162/1‑81 (MQ=255)                    tggttCCGGTAGTGTGATTTTCTCTCATGTGAACGATCCGGGCTTCTGGCTGTTCAAATAGTACTTTAACCTGACTATCgg                                                         >  1:630161/1‑81 (MQ=255)                    tggttCCGGTAGTGTGATTTTCTCTCATGTGAACGATCCGGGCTTCTGGCTGTTCAAATAGTACTTTAACCTGACTATCgg                                                         >  1:630181/1‑81 (MQ=255)                     ggttCCGGTAGTGTGATTTTCTCTCATGTGAACGATCCGGGCTTCTGGCTGTTCAAATAGTACTTTAACCTGACTATCGGc                                                        >  2:630154/1‑81 (MQ=255)                      gttCCGGTAGTGTGATTTTCTCTCATGTGAACGATCCGGGCTTCTGGCTGTTCAAATAGTACTTTAACCTGACTATCGGc                                                        >  1:630170/1‑80 (MQ=255)                               gtgtGATTTTCTCTCATGTGAACGATCCGGGCTTCTGGCTGTTCAAATAGTACTTTAACCTGACTATCGGCGAGACCATCa                                              <  2:630114/81‑1 (MQ=255)                                    aTTTTCTCTCATGTGAACGATCCGGGCTTCTGGCTGTTCAAATAGTACTTTAACCTGACTATCGGCGAGACCATCAAATc                                          <  1:630115/80‑1 (MQ=255)                                        tctctcATGTGAACGATCCGGGCTTCTGGCTGTTCAAATAGTACTTTAACCTGACTATCGGCGAGACCATCAAATCCTgg                                      <  1:630116/80‑1 (MQ=255)                                         ctctcATGTGAACGATCCGGGCTTCTGGCTGTTCAAATAGTACTTTAACCTGACTATCGGCGAGACCATCAAATCCTGGTc                                    >  1:630166/1‑81 (MQ=255)                                            tcATGTGAACGATCCGGGCTTCTGGCTGTTCAAATAGTACTTTAACCTGACTATCGGCGAGACCATCAAATCCTGGTc                                    <  1:630117/78‑1 (MQ=255)                                             cATGTGAACGATCCGGGCTTCTGGCTGTTCAAATAGTACTTTAACCTGATTATCGGCGAGACCATCAAATCCTGGTCGATg                                >  2:630185/1‑81 (MQ=255)                                             cATGTGAACGATCCGGGCTTCTGGCTGTTCAAATAGTACTTTAACCTGACTATCGGCGAGACCATCAAATCCTGGTCGATg                                >  2:630163/1‑81 (MQ=255)                                             cATGTGAACGATCCGGGCTTCTGGCTGTTCAAATAGTACTTTAACCTGACTATCGGCGAGACCATCAAATCCTGGTCGATg                                >  2:630171/1‑81 (MQ=255)                                              aTGTGAACGATCCGGGCTTCTGGCTGTTCAAATAGTACTTTAACCTGACTATCGGCGAGACCATCAAATCCTGGTCGATGc                               <  2:630118/81‑1 (MQ=255)                                               tgtgAACGATCCGGGCTTCTGGCTGTTCAAATAGTACTTTAACCTGACTATCGGCGAGACCATCAAATCCTGGTCGATGCt                              <  1:630119/81‑1 (MQ=255)                                                 tgAACGATCCGGGCTTCTGGCTGTTCAAATAGTACTTTAACCTGACTATCGGCGAGACCATCAAATCCTGGTCGATGCTg                             >  2:630208/1‑80 (MQ=255)                                                  gAACGATCCGGGCTTCTGGCTGTTCAAATAGTACTTTAACCTGACTATCGGCGAGACCATCAAATCCTGGTCGATGCTgg                            >  2:630168/1‑80 (MQ=255)                                                     cGATCCGGGCTTCTGGCTGTTCAAATAGTACTTTAACCTGACTATCGGCGAGACCATCAAATCCTGGTCGATGCTGGAAAc                        >  2:630159/1‑81 (MQ=255)                                                        tCCGGGCTTCTGGCTGTTCAAATAGTACTTTAACCTGACTATCGGCGAGACCATCAAATCCTGGTCGATGCTGGAAACCAt                     <  2:630120/81‑1 (MQ=255)                                                          cGGGCTTCTGGCTGTTCAAATAGTACTTTAACCTGACTATCGGCGAGACCATCAAATCCTGGTCGATGCTGGAAACCATTa                   >  2:630186/1‑81 (MQ=255)                                                           gggCTTCTGGCTGTTCAAATAGTACTTTAACCTGACTATCGGCGAGACCATCAAATCCTGGTCGATGCTGGAAACCATTAt                  >  2:630167/1‑81 (MQ=255)                                                              cTTCTGGCTGTTCAAATAGTACTTTAACCTGACTATCGGCGAGACCATCAAATCCAGGTCGATGCTGGAAACCATTAtctt               <  1:630121/81‑2 (MQ=255)                                                                 cTGGCTGTTCAAATAGTACTTTAACCTGACTATCGGCGAGACCATCAAATCCTGGTCGATGCTGGAAACCATTATCTCgg             <  2:630122/80‑1 (MQ=255)                                                                 cTGGCTGTTCAAATAGTACTTTAACCTGACTATCGGCGAGACCATCAAATCCTGGTCGATGCTGGAAACCATTATCTCGgt            >  1:630172/1‑81 (MQ=255)                                                                   ggCTGTTCAAATAGTACTTTAACCTGACTATCGGCGAGACCATCAAATCCTGGTCGATGCTGGAAACCATTATCTCGgtgt          >  1:630175/1‑81 (MQ=255)                                                                    gCTGTTCAAATAGTACTTTAACCTGACTATCGGCGAGACCATCAAATCCTGGTCGATGCTGGAAACCATTATCTCGgtg           <  1:630124/79‑1 (MQ=255)                                                                    gCTGTTCAAATAGTACTTTAACCTGACTATCGGCGAGACCATCAAATCCTGGTCGATGCTGGAAACCATTATCTCGgtg           <  2:630123/79‑1 (MQ=255)                                                                    gCTGTTCAAATAGTACTTTAACCTGACTATCGGCGAGACCATCAAATCCTGGTCGATGCTGGAAACCATTATCTCGGtgtg         >  2:630183/1‑81 (MQ=255)                                                                       gTTCAAATAGTACTTTAACCTGACTATCGGCGAGACCATCAAATCCTGGTCGATGCTGGAAACCATTATCTCGGTGTGCgg      >  2:630160/1‑81 (MQ=255)                                                                        ttCAAATAGTACTTTAACCTGACTATCGGCGAGACCATCAAATCCTGGTCGATGCTGGAAACCATTATCTCGGTGTGCGGt     >  2:630176/1‑81 (MQ=255)                                                                        ttCAAATAGTACTTTAACCTGACTATCGGCGAGACCATCAAATCCTGGTCGATGCTGGAAACCATTATCTCGGTGTGCGGt     >  2:630178/1‑81 (MQ=255)                                                                          cAAATAGTACTTTAACCTGACTATCGGCGAGACCATCAAATCCTGGTCGATGCTGGAAACCATTATCTCGGTGTGCGGTCt   >  2:630192/1‑81 (MQ=255)                                                                           aaaTAGTACTTTAACCTGACTATCGGCGAGACCATCAAATCCTGGTCGATGCTGGAAACCATTATCTCGGTGTGCGGTCTg  >  1:630165/1‑81 (MQ=255)                                                                           aaaTAGTACTTTAACCTGACTATCGGCGAGACCATCAAATCCTGGTCGATGCTGGAAACCATTATCTCGGTGTGCGGTCTg  >  2:630187/1‑81 (MQ=255)                                                                              |                                                                              AGCTGATGGTTATTGCGGTTGGTTCCGGTAGTGTGATTTTCTCTCATGTGAACGATCCGGGCTTCTGGCTGTTCAAAGAGTACTTTAACCTGACTATCGGCGAGACCATCAAATCCTGGTCGATGCTGGAAACCATTATCTCGGTGTGCGGTCTG  >  NC_000913/3547688‑3547842

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: ‑