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%	167.3 / NA	42	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 (25/17);  total (25/17)

ATAGCATTATTCAACATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACA‑TTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGACGAAAGC  >  REL606/3875550‑3875692                                                                             |                                                                    aTAGCATTATTCAACATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACa                                                                      >  2:2514/1‑76 (MQ=255)   aGCATTATTCAACATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACA‑tt                                                                   <  2:1334505/76‑1 (MQ=255)     cATTATTCAACATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACA‑tttt                                                                 <  2:760147/76‑1 (MQ=255)        tattCAACATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACA‑ttttttt                                                              >  1:604602/1‑76 (MQ=255)        tattCAACATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACA‑ttttttt                                                              <  2:1494420/76‑1 (MQ=255)         attCAACATACGCGTCCTGACTGTAATTTTAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACA‑tttttttt                                                             >  1:867672/1‑75 (MQ=255)          ttCAACATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACA‑ttttttttc                                                            <  1:276171/76‑3 (MQ=255)              aCATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCg                                                         >  1:777729/1‑76 (MQ=255)              aCATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCg                                                         <  2:312493/76‑1 (MQ=255)               cATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGt                                                        >  1:1525921/1‑76 (MQ=255)                aTACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTg                                                       >  1:247280/1‑76 (MQ=255)                aTACGCGACCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTg                                                       >  1:1047255/1‑76 (MQ=255)                 tACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCgaca                                                      >  2:381215/1‑73 (MQ=255)                   cgcgTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGaaa                                                    >  2:351339/1‑76 (MQ=255)                     cgTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATg                                                  >  1:1066366/1‑76 (MQ=255)                     cgTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATg                                                  <  1:867454/76‑1 (MQ=255)                         cTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACa                                              >  2:294797/1‑76 (MQ=255)                         cTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACa                                              >  1:353459/1‑76 (MQ=255)                           gACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACaga                                            >  1:1037958/1‑76 (MQ=255)                            aCTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACagag                                           <  2:1161821/76‑1 (MQ=255)                              tGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGa                                         <  2:1105046/76‑1 (MQ=255)                               gTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGAt                                        <  1:386839/76‑1 (MQ=255)                                  aTTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATaaa                                     <  2:1319703/76‑1 (MQ=255)                                   tttGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAAt                                    <  1:931223/76‑1 (MQ=255)                                   tttGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAAt                                    <  1:1423619/76‑1 (MQ=255)                                    ttGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAAtt                                   >  2:240130/1‑76 (MQ=255)                                      gagaACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGt                                 <  2:1058512/76‑1 (MQ=255)                                         caCGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCa                              >  1:854993/2‑76 (MQ=255)                                           cGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGNATAAATTGTTCAAt                            <  1:886140/76‑1 (MQ=255)                                           cGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAAt                            <  2:1028053/76‑1 (MQ=255)                                                ttAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACgg                       <  1:1488576/76‑1 (MQ=255)                                                 tAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGa                      <  1:686497/76‑1 (MQ=255)                                                  agaCCGCCTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGAt                     >  1:1061483/3‑76 (MQ=255)                                                  aaaCCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGAt                     >  2:461228/1‑76 (MQ=255)                                                   aaCCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGAtt                    >  2:895684/1‑76 (MQ=255)                                                      cGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATc                 >  2:872815/1‑76 (MQ=255)                                                        cTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATccc               <  1:568389/76‑1 (MQ=255)                                                        cTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATccc               >  2:787156/1‑76 (MQ=255)                                                          tCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCg             >  2:1343930/1‑76 (MQ=255)                                                          tCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCg             >  1:438383/1‑76 (MQ=255)                                                               ttGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGAc        >  2:1612610/1‑76 (MQ=255)                                                                tGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGACg       >  2:706272/1‑76 (MQ=255)                                                                   aaaaaCTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGACGaaa    <  1:287944/76‑1 (MQ=255)                                                                   aaaaaCTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGACGaaa    >  1:1350787/1‑76 (MQ=255)                                                                   aaaaaCTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGACGaaa    >  2:106658/1‑76 (MQ=255)                                                                    aaaaCTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGACGAAAg   >  1:1537284/1‑76 (MQ=255)                                                                    aaaaCTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGACGAAAg   >  1:1358812/1‑76 (MQ=255)                                                                     aaaCTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGACGAAAgc  <  2:166577/76‑1 (MQ=255)                                                                     aaaCTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGACGAAAgc  <  2:1270984/76‑1 (MQ=255)                                                                             |                                                                    ATAGCATTATTCAACATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACA‑TTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGACGAAAGC  >  REL606/3875550‑3875692

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