Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,871,835	C→T	intergenic (‑105/+15)	dgoT ← / ← dgoD	putative D‑galactonate transporter/D‑galactonate dehydratase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,871,835	0	C	T	100.0%	104.6 / NA	30	intergenic (‑105/+15)	dgoT/dgoD	putative D‑galactonate transporter/D‑galactonate dehydratase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (16/14);  total (16/14)

TAGCGAAGCCGCTCACCATGCCCTGTAATTTACAGAGGGTTATATTTTTGTATTGCTGTTTTAGTGCCCGATGAGGGGCTTACGTGGCAGGAATTACCACTCTGCTACGCTGTTATCTTCATGACGCCAGAGCGGATTACGCCAGTCCGGGGCAT  >  NC_000913/3871757‑3871911                                                                               |                                                                             tAGCGAAGCCGCTCACCATGCCCTGTAATTTACAGAGGGTTATATTTTTGTATTGCTGTTTTAGTGCCCGATGAGGGGttt                                                                            <  1:765969/81‑1 (MQ=255)   gCGAAGCCGCTCACCATGCCCTGTAATTTACAGAGGGTTATATTTTTGTATTGCTGTTTTAGTGCCCGATGAGGGGTTTAc                                                                          <  2:765970/81‑1 (MQ=255)    cGAAGCCGCTCACCATGCCCTGTAATTTACAGAGGGTTATATTTTTGTATTGCTGTTTTAGTGCCCGATGAGGGGTTTACg                                                                         >  1:766016/1‑81 (MQ=255)     gAAGCCGCTCACCATGCTCTGTAATTTACAGAGGGTTATATTTTTGTATTGCTGTTTTAGTGCCCGATGAGGGGTTTACGt                                                                        >  2:766002/1‑81 (MQ=255)     gAAGCCGCTCACCATGCCCTGTAATTTACAGAGGGTTATATTTTTGTATTGCTGTTTTAGTGCCCGATGAGGGGTTTACGt                                                                        <  2:765971/81‑1 (MQ=255)     gAAGCCGCTCACCATGCCCTGTAATTTACAGAGGGTTATATTTTTGTATTGCTGTTTTAGTGCCCGATGAGGGGTTTACGt                                                                        >  1:766009/1‑81 (MQ=255)       aGCCGCTCACCATGCCCTGTAATTTACAGAGGGTTATATTTTTGTATTGCTGTTTTAGTGCCCGATGAGGGGTTTACGTgg                                                                      >  1:766001/1‑81 (MQ=255)               aCCATGCCCTGTAATTTACAGAGGGTTATATTTTTGTATTGCTGTTTTAGTGCCCGATGAGGGGTTTACGTGGCAGGaa                                                                <  1:765973/79‑1 (MQ=255)               aCCATGCCCTGTAATTTACAGAGGGTTATATTTTTGTATTGCTGTTTTAGTGCCCGATGAGGGGTTTACGTGGCAGGAAtt                                                              <  2:765972/81‑1 (MQ=255)                    gCCCTGTAATTTACCGAGGGTTATATTCTTGTATTGCTGTTTTAGTGTCCGATGAGGGGTTTACGTGGCAGGAATTACCAc                                                         >  1:766013/1‑81 (MQ=255)                    gCCCTGTAATTTACAGAGGGTTATATTTTTGTATTGCTGTTTTAGTGCCCGATGAGGGGTTTACGTGGCAGGAATTACCAc                                                         >  2:766005/1‑81 (MQ=255)                        tGTAATTTACAGAGGGTTATATTTTTGTATTGCTGTTTTAGTGCCCGATGAGGGGTTTACGTGGCAGGAATTACCACTCTg                                                     >  2:766020/1‑81 (MQ=255)                         gTAATTTACAGAGGGTTATATTTTTGTATTGCTGTTTTAGTGCCCGATGAGGGGTTTACGTGGCAGGAATTACCACTCTGc                                                    >  1:766033/1‑81 (MQ=255)                           aaTTTACAGAGGGTTATATTTTTGTATTGCTGTTTTAGTGCCCGATGAGGGGTTTACGTGGCAGGAATTACCACTCTGCTa                                                  >  1:765999/1‑81 (MQ=255)                                  agagGGTTATATTTTTGTATTGCTGTTTTAGTGCCCGATGAGGGGTTTACGTGGCAGGAATTACCACTCTGCTACGCTGtt                                           <  1:765974/81‑1 (MQ=255)                                        ttATATTTTTGTATTGCTGTTTTAGTGCCCGATGAGGGGTTTACGTGGCAGGAATTACCACTCTGCTACGCTGTTATCtt                                      >  2:766010/1‑80 (MQ=255)                                        ttATATTTTTGTATTGCTGTTTTAGTGCCCGATGAGGGGTTTACGTGGCAGGAATTACCACTCTGCTACGCTGTTATCTTc                                     <  2:765975/81‑1 (MQ=255)                                        ttATATTTTTGTATTGCTGTTTTAGTGCCCGATGAGGGGTTTACGTGGCAGGAATTACCACTCTGCTACGCTGTTATCTTc                                     <  2:765976/81‑1 (MQ=255)                                               tttGTATTGCTGTTTTAGTGCCCGATGAGGGGTTTACGTGGCAGGAATTACCACTCTGCTACGCTGTTATCTTCATGACg                               <  2:765977/80‑1 (MQ=255)                                               tttGTATTGCTGTTTTAGTGCCCGATGAGGGGTTTACGTGGCAGGAATTACCACTCTGCTACGCTGTTATCTTCATGACGc                              >  1:766015/1‑81 (MQ=255)                                                 tGTATTGCTGTTTTAGTGCCCGATGAGGGGTTTACGTGGCAGGAATTACCACTCTGCTACGCTGTTATCTTCATGACGcc                             <  2:765978/80‑1 (MQ=255)                                                  gTATTGCTGTTTTAGTGCCCGATGAGGGGTTTACGTGGCAGGAATTACCACTCTGCTACGCTGTTATCTTCATGACGCCa                            <  2:765979/80‑1 (MQ=255)                                                     ttGCTGTTTTAGTGCCCGATGAGGGGTTTACGTGGCAGGAATTACCACTCTGCTACGCTGTTATCTTCATGACGCCAGAGc                        >  2:766014/1‑81 (MQ=255)                                                       gCTGTTTTAGTGCCCGATGAGGGGTTTACGTGGCAGGAATTACCACTCTGCTACGCTGTTATCTTCATGACGCCAGAGCg                       <  2:765980/80‑1 (MQ=255)                                                        cTGTTTTAGTGCCCGATGAGGGGTTTACGTGGCAGGAATTACCACTCTGCTACGCTGTTATCTTCATGACGCCAGAGCGGa                     <  2:765981/81‑1 (MQ=255)                                                          gTTTTAGTGCCCGATGAGGGGTTTACGTGGCAGGAATTACCACTCTGCTACGCTGTTATCTTCATGACGCCAGAGCGGAt                    >  2:766018/1‑80 (MQ=255)                                                                gTGCCCGATGAGGGGTTTACGTGGCAGGAATTACCACTCTGCTACGCTGTTATCTTCATGACGCCAGAGCGGATTACGCCa             <  1:765982/81‑1 (MQ=255)                                                                 tGCCCGATGAGGGGTTTACGTGGCAGGAATTACCACTCTGCTACGCTGTTATCTTCATGACGCCAGAGCGGATTACGCCAg            >  1:766019/1‑81 (MQ=255)                                                                  gCCCGATGAGGGGTTTACGTGGCAGGAATTACCACTCTGCTACGCTGTTATCTTCATGACGCCAGAGCGGATTACGCCAGt           >  2:766023/1‑81 (MQ=255)                                                                    ccGATGAGGGGTTTACGTGGCAGGAATTACCACTCTGCTACGCTGTTATCTTCATGACGCCAGAGCGGATTACGCCAGTcc         >  2:766028/1‑81 (MQ=255)                                                                           ggggTTTACGTGGCAGGAATTACCACTCTGCTACGCTGTTATCTTCATGACGCCAGAGCGGATTACGCCAGTCCGGGGCAt  <  1:765983/81‑1 (MQ=255)                                                                               |                                                                             TAGCGAAGCCGCTCACCATGCCCTGTAATTTACAGAGGGTTATATTTTTGTATTGCTGTTTTAGTGCCCGATGAGGGGCTTACGTGGCAGGAATTACCACTCTGCTACGCTGTTATCTTCATGACGCCAGAGCGGATTACGCCAGTCCGGGGCAT  >  NC_000913/3871757‑3871911

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