Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,871,692	C→T	R13R (CGG→CGA)	dgoT ←	putative D‑galactonate transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,871,692	0	C	T	100.0%	149.2 / NA	37	R13R (CGG→CGA)	dgoT	putative D‑galactonate transporter
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (22/15);  total (22/15)

ACAAATGACTACCGTAATAAAGATCATCACCAGCGTCAGATAACGCCGACGCCCCGGCTTTGCTGCATTAACGGGAATATCCATAGCG  >  NC_000913/3871647‑3871734                                              |                                           aCAAATGACTACCGTAATAAAGATCATCACCAGCGTCAGATAACGTCGa                                         >  1:687502/1‑49 (MQ=255)  cAAATGACTACCGTAATAAAGATCATCGCCAGCGTCAGATAACGTCGAc                                        >  1:1766232/1‑49 (MQ=38)    aaTGACTACCGTAATAAAGATCATCACCAGCGTCAGATAACGTCGACGc                                      <  1:1209317/49‑1 (MQ=255)      tGACTACCGTAATAAAGATCATCACCAGCGTCAGATAACGTCGACGccc                                    >  1:1035435/1‑49 (MQ=255)      tGACTACCGTAATAAAGATCATCACCAGCGTCAGATAACGTCGACGccc                                    >  1:1897835/1‑49 (MQ=255)      tGACTACCGTAATAAAGATCATCACCAGCGTCAGATAACGTCGACGccc                                    >  1:444904/1‑49 (MQ=255)       gACTACCGTAATAAAGATCATCACCAGCGTCAGATAACGTCGACGcccc                                   >  1:1555242/1‑49 (MQ=255)        aCTACCGTAATAAAGATCATCACCAGCGTCAGATAACGTCGACGcccc                                   <  1:1548746/48‑1 (MQ=255)        aCTACCGTAATAAAGATCATCACCAGCGTCAGATAACGTCGACGCCCCg                                  <  1:614370/49‑1 (MQ=255)         cTACCGTAATAAAGATCATCACCAGCGTCAGATAACGTCGACGCCCCgg                                 >  1:1137060/1‑49 (MQ=255)          tACCGTAATAAAGATCATCACCAGCGTCAGATAACGTCGACGCCCCGGc                                >  1:1829877/1‑49 (MQ=255)          tACCGTAATAAAGATCATCACCAGCGTCAGATAACGTCGACGCCCCGGc                                >  1:165652/1‑49 (MQ=255)             cGTAATAAAGATCATCACCAGCGTCAGATAACGTCGACGCCCCGGCttt                             <  1:310297/49‑1 (MQ=255)             cGTAATAAAGATCATCACCAGCGTCAGATAACGTCGACGCCCCGGCttt                             >  1:1554940/1‑49 (MQ=255)              gTAATAAAGATCATCACCAGCGTCAGATAACGTCGACGCCCCGGCTTtg                            <  1:2335396/49‑1 (MQ=255)                aataaAGATCATCACCAGCGTCAGATAACGTCGACGCCCCGGCTTtgct                          >  1:1485547/1‑49 (MQ=255)                  taaAGATCATCACCAGCGTCAGATAACGTCGACGCCCCGGCTTtgctg                         >  1:903475/1‑48 (MQ=255)                    aaGATCATCACCAGCGTCAGATAACGTCGACGCCCCGGCTTTGCTGCAt                      >  1:1663088/1‑49 (MQ=255)                     aGATCATCACCAGCGTCAGATAACGTCGACGCCCCGGCTTTGCTGCAt                      >  1:336157/1‑48 (MQ=255)                      gATCATCACCAGCGTCAGATAACGTCGACGCCCCGGCTTTGCTGCAtt                     >  1:1421473/1‑48 (MQ=255)                      gATCATCACCAGCGTCAGATAACGTCGACGCCCCGGCTTTGCTGCATTa                    <  1:628698/49‑1 (MQ=255)                       atcatcACCAGCGTCAGATAACGTCGACGCCCCGGCTTTGCTGCATTaa                   <  1:1506892/49‑1 (MQ=255)                        tcatcaCCAGCGTCAGATAACGTCGACGCCCCGGCTTTGCTGCATTAAc                  <  1:1103833/49‑1 (MQ=255)                        tcatcaCCAGCGTCAGATAACGTCGACGCCCCGGCTTTGCTGCATTAAc                  <  1:2337387/49‑1 (MQ=255)                        tcatcaCCAGCGTCAGATAACGTCGACGCCCCGGCTTTGCTGCATTAAc                  >  1:329214/1‑49 (MQ=255)                          atcaCCAGCGTCAGATAACGTCGACGCCCCGGCTTTGCTGCATTAACgg                >  1:1097472/1‑49 (MQ=255)                              ccAGCGTCAGATAACGTCGACGCCCCGGCTTTGCTGCATTAACGGGAat            <  1:1692111/49‑1 (MQ=255)                               cAGCGTCAGATAACGTCGACGCCCCGGCTTTGCTGCATTAACGGGAata           >  1:682400/1‑49 (MQ=255)                                aGCGTCAGATAACGTCGACGCCCCGGCTTTGCTGCATTAACGGGAatac          >  1:1709053/1‑48 (MQ=255)                                  cGTCAGATAACGTCGACGCCCCGGCTTTGCTGCATTAACGGGAATATcc        <  1:1350507/49‑1 (MQ=255)                                    tCAGATAACGTCGACGCCCCGGCTTTGCTGCATTAACGGGAATATCCAt      <  1:1656035/49‑1 (MQ=255)                                    tCAGATAACGTCGACGCCCCGGCTTTGCTGCATTAACGGGAATATCCAt      <  1:170942/49‑1 (MQ=255)                                     cAGATAACGTCGACGCCCCGGCTTTGCTGCATTAACGGGAATATCCa       >  1:815467/1‑47 (MQ=255)                                      aGATAACGTCGACGCCCCGGCTTTGCTGCATTAACGGGAATATCCATAg    >  1:202839/1‑49 (MQ=255)                                       gATAACGTCGACGCACCGGCTTTGCTGCATTAACGGGAATATCCATAGc   <  1:102592/49‑1 (MQ=38)                                        aTAACGTCGACGCCCCGGCTTTGCTGCATTAACGGGAATATCCATAGCg  <  1:1499797/49‑1 (MQ=255)                                         tAACGTCGACGCCCCGGCTTTGCTGCATTAACGGGAATATCCATAGc   >  1:590921/1‑47 (MQ=255)                                              |                                           ACAAATGACTACCGTAATAAAGATCATCACCAGCGTCAGATAACGCCGACGCCCCGGCTTTGCTGCATTAACGGGAATATCCATAGCG  >  NC_000913/3871647‑3871734

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