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%	282.3 / NA	69	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/39);  total (30/39)

ATAGCATTATTCAACATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTAC‑A‑TTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGACGAAAGCGCC  >  REL606/3875550‑3875695                                                                              |                                                                       aTAGCATTATTCAACATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTAC‑a                                                                         >  1:144434/1‑76 (MQ=255)   aGCATTATTCAACATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTAC‑A‑tt                                                                      >  2:460689/1‑76 (MQ=255)   aGCATTATTCAACATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTAC‑A‑tt                                                                      <  1:1873076/76‑1 (MQ=255)   aGCATTATTCAACATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTAC‑A‑tt                                                                      >  2:1190951/1‑76 (MQ=255)   aGCATTATTCAACATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTAC‑A‑tt                                                                      >  2:1190936/1‑76 (MQ=255)    gCATTATTCAACATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTAC‑A‑ttt                                                                     >  2:1756621/1‑76 (MQ=255)         attCAACATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTAC‑A‑tttttttt                                                                >  1:1792884/1‑75 (MQ=255)         attCAACATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTAC‑A‑tttttttt                                                                <  2:1411754/76‑2 (MQ=255)               cATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTAC‑ATTTTTTTTCATCGt                                                           <  1:366059/76‑1 (MQ=255)                aTACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTAC‑ATTTTTTTTCATCGTg                                                          <  1:1802588/76‑1 (MQ=255)                aTACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTAC‑ATTTTTTTTCATCGTg                                                          <  1:720784/76‑1 (MQ=255)                  aCGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTAC‑ATTTTTTTTCATCGTGaa                                                        <  1:1822713/76‑1 (MQ=255)                  aCGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTAC‑ATTTTTTTTCATCGTGaa                                                        <  1:1270481/76‑1 (MQ=255)                  aCGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTAC‑ATTTTTTTTCATCGTGaa                                                        >  1:1600020/1‑76 (MQ=255)                    gcgTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTAC‑ATTTTTTTTCATCGTGAAAt                                                      <  2:839940/76‑1 (MQ=255)                     cgTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTAC‑ATTTTTTTTCATCGTGAAATg                                                     >  2:1136274/1‑76 (MQ=255)                       tCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTAC‑ATTTTTTTTCATCGTGAAATGGa                                                   <  1:207838/76‑1 (MQ=255)                        ccTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTAC‑ATTTTTTTTCATTGTGAAATGGAc                                                  <  1:90861/76‑1 (MQ=255)                        ccTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTAC‑ATTTTTTTTCATCGTGAAATGGAc                                                  >  1:697333/1‑76 (MQ=255)                        ccTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTAC‑ATTTTTTTTCATCGTGAAATGGAc                                                  >  1:697315/1‑76 (MQ=255)                        ccTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTAC‑ATTTTTTTTCATCGTGAAATGGAc                                                  <  2:373547/76‑1 (MQ=255)                        ccTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTAC‑ATTTTTTTTCATCGTGAAATGGAc                                                  >  2:412038/1‑76 (MQ=255)                         cTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTAC‑ATTTTTTTTCATCGTGAAATGGACa                                                 <  2:1536105/76‑1 (MQ=255)                         cTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTAC‑ATTTTTTTTCATCGTGAAATGGACa                                                 <  2:46679/76‑1 (MQ=255)                         cTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTAC‑ATTTTTTTTCATCGTGAAATGGACa                                                 <  2:554016/76‑1 (MQ=255)                         cTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTAC‑ATTTTTTTTCATCGTGAAATGGACa                                                 <  2:987371/76‑1 (MQ=255)                         cTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTAC‑ATTTTTTTTCATCGTGAAATGGACa                                                 <  1:1611118/76‑1 (MQ=255)                          tGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTAC‑ATTTTTTTTCATCGTGAAATGGACag                                                <  2:688875/76‑1 (MQ=255)                          tGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTAC‑ATTTTTTTTCATCGTGAAATGGACag                                                <  1:1109642/76‑1 (MQ=255)                          tGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTAC‑ATTTTTTTTCATCGTGAAATGGACag                                                <  2:1267723/76‑1 (MQ=255)                          tGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTAC‑ATTTTTTTTCATCGTGAAATGGACag                                                >  1:896306/1‑76 (MQ=255)                          tGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTAC‑ATTTTTTTTCATCGTGAAATGGACag                                                <  1:1623008/76‑1 (MQ=255)                            aCTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTAC‑ATTTTTTTTCATCGTGAAATGGACagag                                              >  2:1399016/1‑76 (MQ=255)                             cTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTAC‑ATTTTTTTTTATCGTGAAATGGACAGAgg                                             <  1:751114/76‑1 (MQ=255)                              tGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTAC‑ATTTTTTTTCATCGTGAAATGGACAGAGGa                                            <  1:1278215/76‑1 (MQ=255)                               gTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTAC‑ATTTTTTTTCATCGTGAAATGGACAGAGGAt                                           >  2:1599773/1‑76 (MQ=255)                                 aaTTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTAC‑ATTTTTTTTCATCGTGAAATGGACAGAGGATaa                                         <  1:212655/76‑1 (MQ=255)                                   tttGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTAC‑ATTTTTTTTCATCGTGAAATGGACAGAGGATAAAt                                       <  1:701937/76‑1 (MQ=255)                                   tttGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTAC‑ATTTTTTTTCATCGTGAAATGGACAGAGGATAAAt                                       <  1:701935/76‑1 (MQ=255)                                    ttGAGGACGAATTTAAACCGCTTCACCTTGAAAAAACTAC‑ATTTTTTTTCATCGTGAAATGGACAGAGGATAAAtt                                      >  1:1371219/1‑76 (MQ=255)                                    ttGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTAC‑ATTTTTTTTCATCGTGAAATGGACAGAGGATAAAtt                                      <  1:1722566/76‑1 (MQ=255)                                     tGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTAC‑ATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTg                                     <  2:504215/76‑1 (MQ=255)                                      gagaACGAATTTAAACCGCTTCACCTTGAAAAAACTAC‑ATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGt                                    >  1:1055490/1‑76 (MQ=255)                                        gaACGAATTTAAACCGCTTCACCTTGAAAAAACTAC‑ATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTc                                  <  1:1059410/76‑1 (MQ=255)                                         atcGAATTTAACCCGCTTCACCTTGAAAAAACTACGT‑TTTTTTTCAACGTGATCCGGACTGTGGACAAATAGTCCa                                 >  2:1503332/3‑76 (MQ=255)                                         aaCGAATTTAAACCGCTTCACCTTGAAAAAACTAC‑ATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCa                                 <  2:969222/76‑1 (MQ=255)                                         aaCGAATTTAAACCGCTTCACCTTGAAAAAACTAC‑ATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCa                                 <  1:470614/76‑1 (MQ=255)                                         aaCGAATTTAAACCGCTTCACCTTGAAAAAACTAC‑ATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCa                                 <  2:1762856/76‑1 (MQ=255)                                          aCGAATTTAAACCGCTTCACCTTGAAAAAACTAC‑ATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCaa                                <  1:1044870/76‑1 (MQ=255)                                          aCGAATTTAAACCGCTTCACCTTGAAAAAACTAC‑ATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCaa                                <  1:447048/76‑1 (MQ=255)                                           cGAATTTAAACCGCTTCACCTTGAAAAAACTAC‑ATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAAt                               <  1:1165738/76‑1 (MQ=255)                                            gAATTTAAACCGCTTCACCTTGAAAAAACTAC‑ATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATa                              <  1:1676171/76‑1 (MQ=255)                                             aaTTTAAACCGCTTCACCTTGAAAAAACTAC‑ATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATaa                             >  1:1881184/1‑76 (MQ=255)                                              aTTTAAACCGCTTCACCTTGAAAAAACTAC‑ATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAAc                            >  1:891807/1‑76 (MQ=255)                                               tttAAACCGCTTCACCTTGAAAAAACTAC‑ATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACg                           >  1:460974/1‑76 (MQ=255)                                               tttAAACCGCTTCACCTTGAAAAAACTAC‑ATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACg                           <  2:1413107/76‑1 (MQ=255)                                                ttAAACCGCTTCACCTTGAAAAAACTAC‑ATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACgg                          >  2:1684472/1‑76 (MQ=255)                                                  aaaCCGCTTCACCTTGAAAAAACTAC‑ATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGAt                        <  1:1367659/76‑1 (MQ=255)                                                  aaaCCGCTTCACCTTGAAAAAACTAC‑ATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGAt                        >  1:1444348/1‑76 (MQ=255)                                                    aCCGCTTCACCTTGAAAAAACTAC‑ATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTa                      <  1:1691188/76‑1 (MQ=255)                                                       gCTTCACCTTGAAAAAACTAC‑ATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATcc                   >  2:1879224/1‑76 (MQ=255)                                                        cTTCACCTTGAAAAAACTAC‑ATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAAGGGATTATccc                  >  1:1358223/1‑76 (MQ=255)                                                        cTTCACCTTGAAAAAACTAC‑ATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATccc                  >  2:1291000/1‑76 (MQ=255)                                                          tCACCTTGAAAAAACTAC‑ATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCg                >  2:146586/1‑76 (MQ=255)                                                           cACCTTGAAAAAACTAC‑ATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGc               >  1:1296069/1‑76 (MQ=255)                                                           cACCTTGAAAAAACTAC‑ATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGc               >  1:1307928/1‑76 (MQ=255)                                                           cACCTTGAAAAAACTAC‑ATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGc               >  1:599504/1‑76 (MQ=255)                                                               ttGAAAAAACTAC‑ATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGAc           <  1:1281046/76‑1 (MQ=255)                                                               ttGAAAAAACTAC‑ATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGAc           >  1:1417856/1‑76 (MQ=255)                                                               ttGAAAAAACTAC‑ATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGAc           >  1:82925/1‑76 (MQ=255)                                                                 gAAAAAACTAC‑ATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGACGa         >  1:982685/1‑76 (MQ=255)                                                                    aaaaCTAC‑ATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGACGAAAg      <  1:1516515/76‑1 (MQ=255)                                                                     aaaCTAC‑ATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGACGAAAgc     >  2:1854876/1‑76 (MQ=255)                                                                      aaCTAC‑ATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGACGAAAgcg    <  1:75228/76‑1 (MQ=255)                                                                       aCTAC‑ATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGACGAAAgcgc   >  1:237744/1‑76 (MQ=255)                                                                        cTAC‑ATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGACGAAAGCGcc  >  1:609575/1‑76 (MQ=255)                                                                        cTAC‑ATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGACGAAAGCGcc  >  1:180667/1‑76 (MQ=255)                                                                        cTAC‑ATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGACGAAAGCGcc  >  1:477695/1‑76 (MQ=255)                                                                              |                                                                       ATAGCATTATTCAACATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTAC‑A‑TTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGACGAAAGCGCC  >  REL606/3875550‑3875695

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