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%	195.4 / NA	50	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 (30/20);  total (30/20)

ATAGCATTATTCAACATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACA‑TTTTTTT‑CATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGACGAAAGCG  >  REL606/3875550‑3875693                                                                             |                                                                      aTAGCATTATTCAACATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACa                                                                        >  1:268206/1‑76 (MQ=255)   aGCATTATTCAACATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACA‑tt                                                                     <  1:1055818/76‑1 (MQ=255)      attattCAACATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACA‑ttttt                                                                  <  2:189986/76‑1 (MQ=255)        tattCAACATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACA‑ttttttt                                                                >  1:1387910/1‑76 (MQ=255)         attCAACATTCGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACA‑tttttttt                                                               >  2:1639927/1‑75 (MQ=255)         attCAACATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACA‑tttttttt                                                               >  1:483232/1‑75 (MQ=255)         attCAACATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACA‑tttttttt                                                               >  1:388298/1‑75 (MQ=255)               cATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTT‑CATCGt                                                         >  2:1597819/1‑76 (MQ=255)                 tACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTT‑CATCGTGa                                                       <  2:1624562/76‑1 (MQ=255)                   cgcgTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTT‑CATCGTGaaa                                                     <  1:506483/76‑1 (MQ=255)                   cgcgTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTT‑CATCGTGaaa                                                     >  1:731487/1‑76 (MQ=255)                     cgTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTT‑CATCGTGAAATg                                                   <  2:623306/76‑1 (MQ=255)                      gTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTT‑CATCGTGAAATgg                                                  >  1:939658/1‑76 (MQ=255)                      gTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTT‑CATCGTGAAATgg                                                  >  1:939651/1‑76 (MQ=255)                        ccTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTT‑CATCGTGAAATGGAc                                                >  2:1770951/1‑76 (MQ=255)                        ccTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTT‑CATCGTGAAATGGAc                                                <  2:1052582/76‑1 (MQ=255)                           gACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTT‑CATCGTGAAATGGACaga                                             >  2:564592/1‑76 (MQ=255)                           gACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTT‑CATCGTGAAATGGACaga                                             <  1:589774/76‑1 (MQ=255)                               gTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTT‑CATCGTGAAATGGACAGAGGAt                                         <  2:239817/76‑1 (MQ=255)                               gTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTT‑CATCGTGAAATGGACAGAGGAt                                         <  1:768726/76‑1 (MQ=255)                                tAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTT‑CATCGTGAAATGGACAGAGGATa                                        >  2:962443/1‑76 (MQ=255)                                 aaTTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTT‑CATCGTGAAATGGACAGAGGATaa                                       >  1:533389/1‑76 (MQ=255)                                 aaTTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTT‑CATCGTGAAATGGACAGAGGATaa                                       >  2:1123137/1‑76 (MQ=255)                                  aTTTGAGAACGAATTTAAACCGCTTTACCTTGAAAAAACTACATTTTTTTT‑CATCGTGAAATGGACAGAGGATaaa                                      >  2:1489474/1‑76 (MQ=255)                                  aTTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTT‑CATCGTGAAATGGACAGAGGATaaa                                      <  1:686995/76‑1 (MQ=255)                                      gagaACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTT‑CATCGTGAAATGGACAGAGGATAAATTGt                                  >  2:1490861/1‑76 (MQ=255)                                       agaACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTT‑CATCGTGAAATGGACAGAGGATAAATTGtt                                 >  2:833905/1‑76 (MQ=255)                                        gaACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTT‑CATCGTGAAATGGACAGAGGATAAATTGTTc                                <  1:92974/76‑1 (MQ=255)                                         aaCGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTT‑CATCGTGAAATGGACAGAGGATAAATTGTTCa                               <  1:174409/76‑1 (MQ=255)                                         aaCGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTT‑CATCGTGAAATGGACAGAGGATAAATTGTTCa                               <  2:1301791/76‑1 (MQ=255)                                           cGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTT‑CATCGTGAAATGGACAGAGGATAAATTGTTCAAt                             >  1:393318/1‑76 (MQ=255)                                           cGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTT‑CATCGTGAAATGGACAGAGGATAAATTGTTCAAt                             >  1:1349782/1‑76 (MQ=255)                                             aaTTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTT‑CATCGTGAAATGGACAGAGGATAAATTGTTCAATaa                           >  1:381175/1‑76 (MQ=255)                                              aTTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTT‑CATCGTGAAATGGACAGAGGATAAATTGTTCAATAAc                          <  1:901052/76‑1 (MQ=255)                                              aTTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTT‑CATCGTGAAATGGACAGAGGATAAATTGTTCAATAAc                          <  2:1712070/76‑1 (MQ=255)                                               tttAAACCGCTTCACCTTGAAAAAACTACATTTTTTTT‑CATCGTGAAATGGACAGAGGATAAATTGTTCAATAACg                         <  2:42346/76‑1 (MQ=255)                                                 tAAACCGCTTCACCTTGAAAAAACTACATTTTTTTT‑CATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGa                       <  2:299197/76‑1 (MQ=255)                                                 tAAACCGCTTCACCTTGAAAAAACTACATTTTTTTT‑CATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGa                       >  1:1274173/1‑76 (MQ=255)                                                  aaaCCGCTTCACCTTGAAAAAACTACATTTTTTTT‑CATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGAt                      >  1:1555791/1‑76 (MQ=255)                                                    aCCGCTTCACCTTGAAAAAACTACATTTTTTTT‑CATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTa                    <  1:1135237/76‑1 (MQ=255)                                                        cTTCACCTTGAAAAAACTACATTTTTTTT‑CATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATccc                <  2:1815383/76‑1 (MQ=255)                                                        cTTCACCTTGAAAAAACTACATTTTTTTT‑CATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATccc                >  1:268683/1‑76 (MQ=255)                                                          tCACCTTGAAAAAACTTCATTTTTTTT‑CATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCg              >  2:407839/1‑76 (MQ=255)                                                          tCACCTTGAAAAAACTACATTTTTTTT‑CATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCg              >  2:196027/1‑76 (MQ=255)                                                           cACCTTGAAAAAACTACATTTTTTTT‑CATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGc             >  1:1746715/1‑76 (MQ=255)                                                             ccTTGAAAAAACTACATTTTTTTT‑CATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTg           >  2:1514396/1‑76 (MQ=255)                                                             ccTTGAAAAAACTACATTTTTTTT‑CATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTg           >  1:1272507/1‑76 (MQ=255)                                                               ttGAAAAAACTACATTTTTTTT‑CATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGAc         >  2:190011/1‑76 (MQ=255)                                                                tGAAAAAACTACATTTTTTTT‑CATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGACg        <  1:1701725/76‑1 (MQ=255)                                                                tGAAAAAACTACATTTTTTTT‑CATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGACg        <  2:179809/76‑1 (MQ=255)                                                                  aaaaaaCTACATTTTTTTT‑CATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGACGaa      >  1:125585/1‑76 (MQ=255)                                                                   aaaaaCTACATTTTTTTT‑CATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGACGaaa     <  1:775268/76‑1 (MQ=255)                                                                    aaaaCTACATTTTTTTT‑CATCGTGAAATGGACAGATGATAAATTGTTCAATAACGGATTATCCCCGCTGACGAAAg    >  2:1623441/1‑76 (MQ=255)                                                                    aaaaCTACATTTTTTTT‑CATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGACGAAAg    >  2:1095223/1‑76 (MQ=255)                                                                    aaaaCTACATTTTTTTT‑CATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGACGAAAg    >  2:152415/1‑76 (MQ=255)                                                                    aaaaCTACATTTTTTTT‑CATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGACGAAAg    >  1:84082/1‑76 (MQ=255)                                                                      aaCTACATTTTTTTT‑CATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGACGAAAgcg  >  2:1455521/1‑76 (MQ=255)                                                                      aaCTACA‑TTTTTTTACATCGTGAAATGGACAGAGGATAAATTGTTCAATAGCGGATTATCCCCGCTGACGAAAgcg  <  2:347433/76‑1 (MQ=255)                                                                             |                                                                      ATAGCATTATTCAACATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACA‑TTTTTTT‑CATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGACGAAAGCG  >  REL606/3875550‑3875693

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