Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	REL606	3,875,632	(T)7→8	intergenic (‑66/+287)	glmU ← / ← atpC	bifunctional N‑acetylglucosamine‑1‑phosphate uridyltransferase/glucosamine‑1‑phosphate acetyltransferase/F0F1 ATP synthase subunit epsilon

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	REL606	3,875,625	1	.	T	100.0%	105.9 / NA	29	intergenic (‑59/+294)	glmU/atpC	bifunctional N‑acetylglucosamine‑1‑phosphate uridyltransferase/glucosamine‑1‑phosphate acetyltransferase/F0F1 ATP synthase subunit epsilon
Reads supporting (aligned to +/- strand):  ref base . (0/0);  new base T (11/18);  total (11/18)

ATAGCATTATTCAACATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACA‑TTTTT‑TTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGACGAAAGC  >  REL606/3875550‑3875692                                                                             |                                                                     aTAGCATTATTCAACATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACa                                                                       >  2:244913/1‑76 (MQ=255)  tAGCATTATTCAACATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACA‑t                                                                     >  1:362391/1‑76 (MQ=255)         attCAACATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACA‑ttttt‑ttt                                                             >  1:25182/1‑75 (MQ=255)          ttCAACATACGCGTCCTGNCTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACA‑ttttt‑tttc                                                            <  1:845381/76‑3 (MQ=255)              aCATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTT‑TTCATCg                                                         >  2:505630/1‑76 (MQ=255)                aTACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTT‑TTCATCGTg                                                       >  2:105207/1‑76 (MQ=255)                 naCGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTT‑TTCATCGTGa                                                      <  1:939838/75‑1 (MQ=255)                  aCGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTT‑TTCATCGTGaa                                                     <  1:719918/76‑1 (MQ=255)                  aCGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTT‑TTCATCGTGaa                                                     >  1:546064/1‑76 (MQ=255)                    gcgTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTT‑TTCATCGTGAAAt                                                   <  1:60958/76‑1 (MQ=255)                       tCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTT‑TTCATCGTGAAATGGa                                                >  1:332211/1‑76 (MQ=255)                         cTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTT‑TTCATCGTGAAATGGACa                                              <  1:663655/76‑1 (MQ=255)                           gACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTT‑TTCATCGTGAAATGGACaga                                            <  1:482918/76‑1 (MQ=255)                            aCTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTT‑TTCATCGTGAAATGGACagag                                           >  1:1033643/1‑76 (MQ=255)                              tGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTT‑TTCATCGTGAAATGGACAGAGGa                                         <  2:142954/76‑1 (MQ=255)                               gTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTT‑TTCATCGTGAAATGGACAGAGGAt                                        <  1:750895/76‑1 (MQ=255)                               gTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTT‑TTCATCGTGAAATGGACAGAGGAt                                        <  2:522090/76‑1 (MQ=255)                                tAATTTGAGAAGGAATTAAAACCGCTTCACGATGACAAAACTNCN‑TTTTTCTNCATCGTGAAATGGACAGAGGATa                                       <  1:999073/76‑1 (MQ=255)                                tAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTT‑TTCATCGTGAAATGGACAGAGGATa                                       <  2:362391/76‑1 (MQ=255)                                 aacttGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTT‑TTCATCGTGAAATGGACAGAGGATaa                                      <  1:518408/73‑1 (MQ=255)                                    ntGNGNACTAAGTGNAATCGCTTGACATGGAANAAACTACAGTTTTT‑TTCATCGTGAAATGGACAGAGGATaaagt                                   <  1:916139/75‑3 (MQ=255)                                            gAATTTAAACCGCTTCACCTTGAAAAAACTACANTTTTT‑TTCATCGTGAAATGGACAGAGGATAAATTGTTCAATa                           <  1:875866/76‑1 (MQ=255)                                             aaTTTAAACCGCTTCACCTTGAAAAANCTACATTTTTT‑TTCATCGTGAAATGGACAGAGGATAAATTGTTCACTaa                          >  2:933697/1‑76 (MQ=255)                                             aaTTTAAACCGCTTCACCTTGAAAAAACTACATTTTTT‑TTCATCGTGAAATGGACAGAGGATAAATTGTTCAATan                          >  1:947380/1‑75 (MQ=255)                                                ttAAACCGCTTCACCTTGAAAAAACTACATTTTTT‑TTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACgg                       <  1:39016/76‑1 (MQ=255)                                                      cGCTTCACCTTGAAAAAACTACATTTTTT‑TTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATc                 >  1:249034/1‑76 (MQ=255)                                                       gCTTCACCTTGAAAAAACTACATTTTTT‑TTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATcc                <  2:662308/76‑1 (MQ=255)                                                            aCCTTGAAAAAACTACATTTTTT‑TTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCt           <  2:22609/76‑1 (MQ=255)                                                            aCCTTGAAAAAACTACATTTTTT‑TTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCt           <  1:807648/76‑1 (MQ=255)                                                             ccTTGAAAAAACTACATTTTTT‑TTCATCGTGAAATGGATAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTg          <  2:198308/76‑1 (MQ=255)                                                             ccTTGAAAAAACTACATTTTTT‑TTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTg          >  2:247421/1‑76 (MQ=255)                                                             ccTTGAAAAAACTACATTTTTT‑TTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTg          <  2:609291/76‑1 (MQ=255)                                                              cTTGAAAAAACTACATTTTTT‑TTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGa         <  2:435344/76‑1 (MQ=255)                                                               ttGAAAAAACTACATTTTTT‑TTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGAc        >  1:896360/1‑76 (MQ=255)                                                                    aaaaCTACATTTTTT‑TTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGACGAAAg   <  1:393499/76‑1 (MQ=255)                                                                     aaaCTACATTTTTT‑TTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGACGAAAgc  >  1:405938/1‑76 (MQ=255)                                                                             |                                                                     ATAGCATTATTCAACATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACA‑TTTTT‑TTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGACGAAAGC  >  REL606/3875550‑3875692

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