Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	REL606	3,875,625	+T	intergenic (‑59/+294)	glmU/atpC	bifunctional N‑acetylglucosamine‑1‑phosphate uridyltransferase/glucosamine‑1‑phosphate acetyltransferase/F0F1 ATP synthase subunit epsilon

Read alignment evidence...
	seq id	position		change	freq	score	reads	annotation	genes	product
*	REL606	3,875,625	1	.→T	100.0%	78.8	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):  new base (16/11):  ref base (1/1):  total (17/12)

AACGAATTTAAACCGCTTCACCTTGAAAAAACTACA.TTT.T.TTTCATCGTGAAATGGACAGAGGATA  >  REL606/3875590‑3875655                                     |                                 aaCGAATTTAAACCGCTTCACCTTGAAAAAACTACa                                   >  1:5988844/1‑36 (MQ=255) aaCGAATTTAAACCGCTTCACCTTGAAAAAACTACa                                   >  1:6720278/1‑36 (MQ=255)  aCGAATTTAAACCGCTTCACCTTGAAAAAACTACA.t                                 >  1:595072/1‑36 (MQ=255)   gaaaTTTAACCCGCTTCACCTTGAAAAAACTACA.tt                                <  1:885242/34‑1 (MQ=25)   cGAATTTAAACCGCTTCACCTTGAAAAAACTACA.tt                                <  1:590454/36‑1 (MQ=255)   cGAATTTAAACCGCTTCACCTTGAAAAAACTACA.tt                                >  1:2769620/1‑36 (MQ=255)    gAATTTAGACCGCTTCACCTTGAAAAAACTACA.ttt                               <  1:6131610/36‑1 (MQ=255)    gAATTTAAACCGCTTCACCTTGAAAAAACTACA.ttt                               <  1:3336338/36‑1 (MQ=255)    gAATATAAACCGCTTCACCTTGAAAAAACTACA.ttt                               >  1:6941160/1‑36 (MQ=255)     aaTTTAAACCGCTTCACCTTGAAAAAACTACA.ttt.t                             >  1:1546860/1‑36 (MQ=255)      aTTTAAACCGCTTCACCTTGAAAAAACTACA.ttt.t.t                           >  1:1141827/1‑36 (MQ=255)      aTTTAAACCGCTTCACCTTGAAAAAACTACA.ttt.t.t                           >  1:3566197/1‑36 (MQ=255)      aTTTAAACAGCATCACCTTGAAAAAACTACA.ttg.t.t                           <  1:6097407/36‑1 (MQ=21)       tttAAACCGCTTCACCTTGAAAAAACTACA.ttt.t.tt                          <  1:6098064/36‑1 (MQ=255)        ttAAACCGCTTCACCTTGAAAAAACTACA.tat.a.ttt                         >  1:3537153/1‑36 (MQ=25)          aaaCCGCTTCACCTTGAAAAAACTACA.ttt.t.ttttc                       >  1:2036068/1‑34 (MQ=255)           aaCCGCTTCACCTTGAAAAAACTACA.ttt.t.ttttca                      <  1:5922249/36‑4 (MQ=255)           aaCCGCTTCACCTTGAAAAAACTACA.ttt.t.ttttca                      >  1:3760410/1‑33 (MQ=255)            aCCGCTTCACCTTGAAAAAACTACA.TTTAT.TTTCAt                      >  1:5755026/1‑36 (MQ=25)              cGCTTCACCTTGAAAAAACTACATTTT.T.TTTCATCg                    >  1:411427/1‑36 (MQ=255)               gCTTCACCTTGAAAAAACTACATTTT.T.TTTCATCGt                   <  1:4443459/36‑1 (MQ=25)               gCTTCACCTTGAAAAAACTACATTTT.T.TTTCATCGt                   <  1:6770695/36‑1 (MQ=25)               gCTTCACCTTGAAAAAACTACATTTT.T.TTTCATCGt                   <  1:173527/36‑1 (MQ=25)               ccTTCACCTTGAAAAAACTACATTTT.T.TTTCATCGt                   >  1:758935/2‑36 (MQ=21)                 ttCACCTTGAAAAAACTACATTTT.T.TTTCATCGTGa                 >  1:6108888/1‑36 (MQ=25)                 ttCACCTTGAAAAAACTACATTTT.T.TTTCATCGTGa                 >  1:2823563/1‑36 (MQ=25)                  tCACCTTGAAAAAACTACATTTT.T.TTTCATCGTGaa                <  1:2744376/36‑1 (MQ=25)                   cACCTTGAAAAAACTACATTTT.T.TTTCATCGTGaaa               <  1:6197661/36‑1 (MQ=25)                   cACCTTGAAAAAACTACATTTT.T.TTTCATCGTGaaa               <  1:1653252/36‑1 (MQ=25)                    aCCTTGAAAAAACTACATTTT.T.TTTCATCGTGAAAt              >  1:6484497/1‑36 (MQ=25)                    aCCTTGAAAAAACTACATTTT.T.TTTCATCGTGAAAt              >  1:3453267/1‑36 (MQ=25)                     ccTTGAAAAAACTACATTTT.T.TTTCATCGTGAAATg             <  1:1138005/36‑1 (MQ=25)                      cTTGAAAAAACTACATTTT.T.TTTCATCGTGAAATgg            >  1:5853925/1‑36 (MQ=25)                      cTTGAAAAAACTACATTTT.T.TTTCATCGTGAAATgg            >  1:5610010/1‑36 (MQ=25)                        tGAAAAAACTACATTTT.T.TTTCATCGTGAAATGGAc          <  1:4391506/36‑1 (MQ=25)                          aaaaaaCTACATTTT.T.TTTCATCGTGAACTGGCCag        >  1:5560754/1‑36 (MQ=18)                          aaaaaaCTACATTTT.T.TTTCATCGTGAAATGGACag        >  1:2992163/1‑36 (MQ=25)                          aaaaaaCTACA.TTT.TATTTCATCGTGAAATGGACag        <  1:6809251/36‑1 (MQ=25)                           aaaaaCTACATTTT.T.TTTCATCGTGAAATGGACaga       >  1:5106536/1‑36 (MQ=25)                           aaaaaCTACATTTT.T.TTTCATCGTGAAATGGACaga       >  1:6673486/1‑36 (MQ=25)                           aaaaaCTACATTTT.T.TTTCATCGTGAAATGGACaga       <  1:3550583/36‑1 (MQ=25)                           aaaaaCTACATTTT.T.TTTCATCGTGAAATGCACaga       >  1:3537482/1‑36 (MQ=21)                            aaaaCTACATTTT.T.TTTCATCGTGAAATGGACagag      <  1:4837616/36‑1 (MQ=25)                            aaaaCTACATTTT.T.TTTCATCGTGAAATGGACagag      >  1:3295603/1‑36 (MQ=25)                             aaaCTACATTTT.T.TTTCATCGTGAATTGGACAGAgg     >  1:562490/1‑36 (MQ=21)                               aCTACATTTT.T.TTTCATCGTGAAATGGACAGAGGAt   <  1:7083489/36‑1 (MQ=25)                                cTACATTTT.T.TTTCATCGTGAAATGGACAGAGGATa  >  1:4190902/1‑36 (MQ=25)                                     |                                 AACGAATTTAAACCGCTTCACCTTGAAAAAACTACA.TTT.T.TTTCATCGTGAAATGGACAGAGGATA  >  REL606/3875590‑3875655 Base quality scores:  ATCG  < 3 ≤  ATCG  < 11 ≤  ATCG  < 19 ≤  ATCG  < 30 ≤  ATCG  < 32 ≤  ATCG