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%	108.2 / NA	30	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 (14/16);  total (14/16)

TAGCATTATTCAACATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACA‑TT‑TTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGACGAA  >  REL606/3875551‑3875689                                                                            |                                                                  tAGCATTATTCAACATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACA‑t                                                                  >  1:944028/1‑76 (MQ=255) tAGCATTATTCAACATACGCGTCCTGACTCTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACA‑t                                                                  >  1:1143407/1‑76 (MQ=255)    cATTATTCAACATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACA‑tt‑tt                                                              <  1:703828/76‑1 (MQ=255)    cATTATTCAACATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACA‑tt‑tt                                                              <  2:808388/76‑1 (MQ=255)       tattCAACATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACA‑tt‑ttttt                                                           <  2:329849/76‑1 (MQ=255)       tattCAACATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACA‑tt‑ttttt                                                           >  1:302532/1‑76 (MQ=255)       tattCAACATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCACTTCACCTTGAAAAAACTACA‑tt‑ttttt                                                           <  2:1123599/76‑1 (MQ=255)          tCAACATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACA‑tt‑ttttttca                                                        >  2:759240/1‑73 (MQ=255)                 aCGCGTCCTGACTGTATTTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTT‑TTTTTCATCGTGaa                                                  <  2:765358/76‑1 (MQ=255)                 aCGCGTCCTGACTGTAATTTGAGAACGCATTTAAACCGCTTCACCTTGAAAAAACTACATTT‑TTTTTCATCGTGaa                                                  <  1:655728/76‑1 (MQ=255)                 aCGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTT‑TTTTTCATCGTGan                                                  <  1:413748/76‑2 (MQ=255)                 aCGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTT‑TTTTTCATCGTGaa                                                  >  1:344618/1‑76 (MQ=255)                 aCGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTT‑TTTTTCATCGTGaa                                                  <  1:860056/76‑1 (MQ=255)                  cgcgTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTT‑TTTTTCATCGTGaaa                                                 >  2:451223/1‑76 (MQ=255)                  cgcgTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTT‑TTTTTCATCGTGaaa                                                 <  2:832880/76‑1 (MQ=255)                      tCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTT‑TTTTTCATCGTGAAATGGa                                             <  1:317854/76‑1 (MQ=255)                        cTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTT‑TTTTTCATCGTGAAATGGACa                                           <  2:944028/76‑1 (MQ=255)                        cTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACGACATTT‑TTTTTCATCGTGAAATGGACa                                           <  2:1143405/76‑1 (MQ=255)                         tGACTGTAATTTAAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTT‑TTTTTCATCGTGAAATGGACag                                          >  2:633597/1‑76 (MQ=255)                             tGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTT‑TTTTTCATCGTGAAATGGACAGAGGa                                      >  1:372689/1‑76 (MQ=255)                              gTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTT‑TTTTTCATCGTGAAATGGACAGAGGAt                                     >  2:83841/1‑76 (MQ=255)                                aaTTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTT‑TTTTTCAGCGTGAAATGGACAGAGGATaa                                   >  2:969747/1‑76 (MQ=255)                                        aaCGAATTTAAACCGCTTCACCTTGAAAAAACTACATTT‑TTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCa                           <  2:794192/76‑1 (MQ=255)                                        aaCGAATTTAAACCGCTTCACCTTGAAAAAACTACATTT‑TTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCa                           <  1:541708/76‑1 (MQ=255)                                        aaCGAATTTAAACCGCTTCACCTTGAAAAAACTACATTT‑TTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCa                           >  2:53500/1‑76 (MQ=255)                                         aCGAATTTAAACCGCTTCACCTTGAAAAAACTACATTT‑TTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCaa                          <  2:108776/76‑1 (MQ=255)                                           gtaTTTAAACCACTTGAGCTTGAAAAAACTGCA‑TTCTTTTTCNTNGTGAACTNGACAGAGGATAAATTGTTCAATa                        <  1:1077110/74‑1 (MQ=255)                                           gAATTTAAACCGCTTCACCTTGAAAAAACTACATTT‑TTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATa                        <  2:303890/76‑1 (MQ=255)                                            aaTTTAAACCGCTTCACCTTGAAAAAACTACATTT‑TTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATaa                       <  2:40562/76‑1 (MQ=255)                                             aTTTAAACCGCTTCACCTTGAAAAAACTACATTT‑TTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAAc                      <  1:859986/76‑1 (MQ=255)                                              tttAAACCGCTTCACCTTGAAAAAACTACATTT‑TTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACg                     >  2:359076/1‑76 (MQ=255)                                                 aaaCCGCTTCACCTTGAAAAAACTACATTT‑TTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGAt                  >  1:645822/1‑76 (MQ=255)                                                   aCCGCTTCACCTTGAAAAAACTACATTT‑TTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTa                >  1:747041/1‑76 (MQ=255)                                                     cGCTTCACCTTGAAAAAACTACATTT‑TTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATc              <  1:83841/76‑1 (MQ=255)                                                           aCCCTGAAAAAACTACATTT‑TTTTTCATCGTGAAATGGACAGAGGATAANTTGTTCAATAACGGATTATCCCCGCt        <  2:994937/76‑1 (MQ=255)                                                             ctagAAAAAACTACATTT‑TTTTTCATNGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGATGa      >  2:999319/4‑76 (MQ=255)                                                             cTTGAAAAAACTACATTT‑TTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGa      >  1:1155585/1‑76 (MQ=255)                                                               tGAAAAAACTACATTT‑TTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGACg    >  1:862068/1‑76 (MQ=255)                                                                 aaaaaaCTACATTT‑TTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGACGaa  >  1:524042/1‑76 (MQ=255)                                                                            |                                                                  TAGCATTATTCAACATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACA‑TT‑TTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGACGAA  >  REL606/3875551‑3875689

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