Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	REL606	3,499,711	A→C	S365A (TCA→GCA)	glgC ←	glucose‑1‑phosphate adenylyltransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	REL606	3,499,711	0	A	C	100.0%	130.2 / NA	32	S365A (TCA→GCA)	glgC	glucose‑1‑phosphate adenylyltransferase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base C (15/17);  total (15/17)

CAGCGGCGCAGACGGCACGAGCGACCTACCCATACTTCCGGTAACAATACGGCGGAATCAATGTTGCAGAATGAATTCACGCGAACGCGCGAGAACAGAACGGACTGCACCACCACCGAACCGGAGATCACACAACCGCCGGAAAC  >  REL606/3499638‑3499783                                                                          |                                                                         cAGCGGCGCAGACGGCACGAGCGACCTACCCATACTTCCGGTAACAATACGGCGGAATCAATGTTGCAgaatgcat                                                                        <  1:139611/76‑1 (MQ=255)      gcgcAGACGGCACGAGCGACCTACCCATACTTCCGGTAACAATACGGCGGAATCAATGTTGCAGAATGCATTCAcg                                                                   <  2:497349/76‑1 (MQ=255)      gcgcAGACGGCACGAGCGACCTACCCATACTTCCGGTAACAATACGGCGGAATCAATGTTGCAGAATGCATTCAcg                                                                   <  1:794464/76‑1 (MQ=255)       cgcAGACGGCACGAGCGACCTACCCATACTTCCGGTAACAATACGGCGGAATCAATGTTGCAGAATGCATTCAcgc                                                                  >  1:748502/1‑76 (MQ=255)        gcAGACGGCACGAGCGACCTACCCATACTTCCGGTAACAATACGGCGGAATCAATGTTGCAGAATGCATTCAcgcg                                                                 <  1:1055792/76‑1 (MQ=255)         cggACGGCACGAGCGACCTACCCATACTTCCGGTAACAATACGGCGGAATCAATGTTGCAGAATGCATTCACGCGa                                                                >  1:1307994/3‑76 (MQ=255)          aGACGGCACGAGCGACCTACCCATACTTCCGGTAACAATACGGCGGAATCAATGTTGCAGAATGAATTCACGCGaa                                                               >  1:1336208/1‑76 (MQ=255)           gACGGCACGAGCGACCTACCCATACTTCCGGTAACAATACGGCGGAATCAATGTTGCAGAATGCATTCACGCGAAc                                                              >  1:1461808/1‑76 (MQ=255)                cACGAGCGACCTACCCATACTTCCGGTAACAATACGGCGGAATCAATGTTGCAGAATGCATTCACGCGAAcgcgcg                                                         >  2:1617341/1‑76 (MQ=255)                cACGAGCGACCTACCCATACTTCCGGTAACAATACGGCGGAATCAATGTTGCAGAATGCATTCACGCGAAcgcgcg                                                         >  2:1272530/1‑76 (MQ=255)                  cGAGCGACCTACCCATACTTCCGGTAACAATACGGCGGAATCAATGTTGCAGAATGCATTCACGCGAACGCGCgag                                                       <  1:1345340/76‑1 (MQ=255)                  cGAGCGACCTACCCATACTTCCGGTAACAATACGGCGGAATCAATGTTGCAGAATGCATTCACGCGAACGCGCgag                                                       <  2:1456114/76‑1 (MQ=255)                     gCGACCTACCCATACTTCCGGTAACAATACGGCGGAATCAATGTTGCAGAATGCATTCACGCGAACGCGCGagaac                                                    <  1:1570331/76‑1 (MQ=255)                      cGACCTACCCATACTTCCGGTAACAATACGGCGGAATCAATGTTGCAGAATGCATTCACGCGAACGCGCGagaaca                                                   <  1:344218/76‑1 (MQ=255)                       gACCTACCCATACTTCCGGTAACAATACGGCGGAATCAATGTTGCAGAATGCATTCACGCGAACGCGCGagaacag                                                  >  2:970817/1‑76 (MQ=255)                       gACCTACCCATACTTCCGGTAACAATACGGCGGAATCAATGTTGCAGAATGCATTCACGCGAACGCGCGagaacag                                                  <  2:941577/76‑1 (MQ=255)                        aCCTACCCATACTTCCGGTAACAATACGGCGGAATCAATGTTGCAGAATGCATTCACGCGAACGCGCGagaacaga                                                 <  1:1548344/76‑1 (MQ=255)                          cTACCCATACTTCCGGTAACAATACGGCGGAATCAATGTTGCAGAATGCATTCACGCGAACGCGCGagaacagaac                                               <  1:679506/76‑1 (MQ=255)                            aCCCATACTTCCGGTAACAATACGGCGGAATCAATGTTGCAGAATGCATTCACGCGAACGCGCGAGAACAGAACgg                                             >  2:581319/1‑76 (MQ=255)                             cccATACTTCCGGTAACAATACGGCGGAATCAATGTTGCAGAATGCATTCACGCGAACGCGCGAGAACAGAACGGa                                            <  2:1212483/76‑1 (MQ=255)                                 caCTTCCGGTAACAATACGGCGGAATCAATGTTGAAGAATGCATTCACGAGATCGCGCGAGAACAGCACGGACTGc                                        >  2:416656/2‑76 (MQ=255)                                  aCTTCCGGTAACAATACGGCGGAATCAATGTTGCAGAATGCATTCACGCGAACGCGCGAGAACAGAACGGACTGca                                       <  2:884874/76‑1 (MQ=255)                                     tCCGGTAACAATACGGCGGAATCAATGTTGCAGAATGCATTCACGCGAACGCGCGAGAACAGAACGGACTGcacca                                    >  1:1593598/1‑76 (MQ=255)                                     tCCGGTAACAATACGGCGGAATCAATGTTGCAGAATGCATTCACGCGAACGCGCGAGAACAGAACGGACTGcacca                                    <  1:1588904/76‑1 (MQ=255)                                              aaTACGGCGGAATCAATGTTGCAGAATGCATTCACGCGAACGCGCGAGAACAGAACGGACTGCACCACCACCGAAc                           <  1:413617/76‑1 (MQ=255)                                                  cggcggAATCAATGTTGCAGAATGCATTCACGCGAACGCGCGAGAACAGAACGGACTGCACCACCACCGAACCGga                       <  1:272771/76‑1 (MQ=255)                                                        aaTCAATGTTGCAGAATGCATTCACGCGAACGCGCGAGAACAGAACGGACTGCACCACCACCGAACCGGAGATcac                 <  2:132723/76‑1 (MQ=255)                                                         aTCAATGTTGCAGAATGCATTCACGCGAACGCGCGAGAACAGAACGGACTGCACCACCACCGAACCGGAGATcaca                >  2:946962/1‑76 (MQ=255)                                                          tCAATGTTGCAGAATGCATTCACGCGAACGCGCGAGAACAGAACGGACTGCACCACCACCGAACCGGAGATcacac               <  2:770918/76‑1 (MQ=255)                                                            aaTGTTGCAGAATGCATTCACGCGAACGCGCGAGAACAGAACGGACTGCACCACCACCGAACCGGAGATCACACaa             >  1:665029/1‑76 (MQ=255)                                                              tGTTGCAGAATGCATTCACGCGAACGCGCGAGAACAGAACGGACTGCACCACCACCGAACCGGAGATCACACAAcc           >  1:814188/1‑76 (MQ=255)                                                                ttGCAGAATGCATTCACGCGAACGCGCGAGAACAGAACGGACTGCACCACCACCGAACCGGAGATCACACAAccgc         >  2:370580/1‑76 (MQ=255)                                                                   cAGAATGCATTCACGCGAACGCGCGAGAACAGAACGGACTGCACCACCACCGAACCGGAGATCACACAACCGCCgg      >  2:1587713/1‑76 (MQ=255)                                                                    aGAATGCATTCACGCGAACGCGCGAGAACAGAACGGACTGCACCACCACCGAACCGGAGATCACACAACCGCCGGa     >  1:1560724/1‑76 (MQ=255)                                                                      aatgcatTCACGCGAACGCGCGAGAACAGAACGGACTGCACCACCACCGAACTGGAGATCACACAACCGCCGGaaa   >  2:1361308/1‑76 (MQ=255)                                                                      aatgcatTCACGCGAACGCGCGAGAACAGAACGGACTGCACCACCACCGAACCGGAGATCACACAACCGCCGGaaa   >  2:1102054/1‑76 (MQ=255)                                                                       atgcatTCACGCGAACGCGCGAGAACAGAACGGACTGCACCACCACCGAACCGGAGATCACACAACCGCCGGAAAc  <  1:1184695/76‑1 (MQ=255)                                                                          |                                                                         CAGCGGCGCAGACGGCACGAGCGACCTACCCATACTTCCGGTAACAATACGGCGGAATCAATGTTGCAGAATGAATTCACGCGAACGCGCGAGAACAGAACGGACTGCACCACCACCGAACCGGAGATCACACAACCGCCGGAAAC  >  REL606/3499638‑3499783

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