Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	REL606	3,462,730	G→A	G410D (GGC→GAC)	pckA →	phosphoenolpyruvate carboxykinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	REL606	3,462,730	0	G	A	100.0%	136.3 / NA	37	G410D (GGC→GAC)	pckA	phosphoenolpyruvate carboxykinase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (15/22);  total (15/22)

GGCTTCACCGCCAAACTGGCCGGTACTGAGCGTGGCATCACCGAACCGACGCCAACCTTCTCCGCTTGCTTCGGCGCGGCATTCCTGTCGCTGCACCCGACTCAGTACGCAGAAGTGCTGGTGAAACGTATGCAGGCGGCGGGCGC  >  REL606/3462657‑3462802                                                                          |                                                                         ggCTTCACCGCCAAACTGGCCGGTACTGAGCGTGGCATCACCGAACCGACGCCAACCTTCTCCGCTTGCTTCGacg                                                                        <  2:13540/76‑1 (MQ=255) ggCTTCACCGCCAAACTGGCCGGTACTGAGCGTGGCATCACCGAACCGACACCAACCTTCTCCCCTTGTTTCGacg                                                                        >  2:1368816/1‑76 (MQ=255)     tCACCGCCAAACTGGCCGGTACTGAGCGTGGCATCACCGAACCGACGCCAACCTTCTCCGCTTGCTTCGACGCGGc                                                                    >  1:981717/1‑76 (MQ=255)         cGCCAAACTGGCCGGTACTGAGCGTGGCATCACCGAACCGACGCCAACCTTCTCCGCTTGCTTCGACGCGGCATTc                                                                >  2:360442/1‑76 (MQ=255)          gCCAAACTGGCCGGTACTGAGCGTGGCATCACCGAACCGACGCCAACCTTCTCCGCTTGCTTCGACGCGGCATTcc                                                               >  1:1489595/1‑76 (MQ=255)          gCCAAACTGGCCGGTACTGAGCGTGGCATCACCGAACCGACGCCAACCTTCTCCGCTTGCTTCGACGCGGCATTcc                                                               >  1:849348/1‑76 (MQ=255)           ccAAACTGGCCGGTACTGAGCGTGGCATCACCGAACCGACGCCAACCTTCTCCGCTTGCTTCGACGCGGCATTCCt                                                              <  1:549889/76‑1 (MQ=255)              aaCTGGCCGGTACTGAGCGTGGCATCACCGAACCGACGCCAACCTTCTCCGCTTGCTTCGACGCGGCATTCCTGTc                                                           >  2:994130/1‑76 (MQ=255)                cTGGCCGGTACTGAGCGTGGCATCACCGAACCGACGCCAACCTTCTCCGCTTGCTTCGACGCGGCATTCCTGTCGc                                                         >  1:813509/1‑76 (MQ=255)                 tGGCCGGTACTGAGCGTGGCATCACCGAACCGACGCCAACCTTCTCCGCTTGCTTCGACGCGGCATTCCTGTCGCt                                                        <  2:1354895/76‑1 (MQ=255)                        tACTGAGCGTGGCATCACCGAACCGACGCCAACCTTCTCCGCTTGCTTCGACGCGGCATTCCTGTCGCTGCACCCg                                                 <  1:1238404/76‑1 (MQ=255)                        tACTGAGCGTGGCATCACCGAACCGACGCAAACCTTCTCCGCTTGCTTCGACGCGGCATTCCTGTCGCTGCACCCg                                                 <  2:313264/76‑1 (MQ=255)                          cTGAGCGTGGCATCACCGAACNGACGCCAACCTTCTCCTCTTGCTTCGACGCGGCATTCCTGTCGCTGCACCCGAc                                               <  2:295684/76‑1 (MQ=255)                            gcgCGTGGCATCACCGAACCGACGCCAACCTTCTCCGCTTGCTTCGACGCGGCATTCCTGTCGCTGCACCCGACTc                                             <  2:905001/74‑1 (MQ=255)                               cGTGGCATCACCGAACCGACGCCAACCTTCTCCGCTTGCTTCGACGCGGCATTCCTGTCGCTGCACCCGACTCAGt                                          <  1:629656/76‑1 (MQ=255)                                gTGGCATCACCGAACCGACGCCAACCTTCTCCGCTTGCTTCGACGCGGCATTCCTGTCGCTGCACCCGACTCAGTa                                         <  2:1209558/76‑1 (MQ=255)                                 tgccATCACCGAACCGACGCCAACCTTCTCCGCTTGCTTCGACGCGGCATTCCTGTCGCTGCACCCGACTCAGTAc                                        <  1:226325/73‑1 (MQ=255)                                 tGGCATCACCGAACCGACGCCAACCTTCTCCGCTTGCTTCGACGCGGCATTCCTGTCGCTGCACCCGACTCAGTAc                                        <  2:843876/76‑1 (MQ=255)                                 tGGCATCACCGAACCGACGCCAACCTTCTCCGCTTGCTTCGACGCGGCATTCCTGTCGCTGCACCCGACTCAGTAc                                        <  1:244671/76‑1 (MQ=255)                                   gCATCACCGAACCGACGCCAACCTTCTCCGCTTGCTTCGACGCGGCATTCCTGTCGCTGCACCCGACTCAGTACGc                                      <  2:148819/76‑1 (MQ=255)                                    cATCACCGAACCGACGCCAACCTTCTCCGCTTGCTTCGACGCGGCATTCCTGTCGCTGCACCCGACTCAGTACGCa                                     >  1:778966/1‑76 (MQ=255)                                        accgaaccgaCGCCAACCTTCTCCGCTTGCTTCGACGCGGCATTCCTGTCGCTGCACCCGACTCAGTACGCAGAAg                                 <  1:65925/76‑1 (MQ=255)                                         ccgaaccgaCGCCAACCTTCTCCGCTTGCTTCGACGCGGCATTCCTGTCGCTGCACCCGACTCAGTACGCAGAAGt                                <  2:75633/76‑1 (MQ=255)                                            aaccgaCGCCAACCTTCTCCGCTTGCTTCGACGCGGCATTCCTGTCGCTGCACCCGACTCAGTACGCAGAAGTGcg                             >  1:306925/1‑75 (MQ=255)                                            aaccgaCGCCAACCTTCTCCGCTTGCTTCGACGCGGCATTCCTGTCGCTGCACCCGACTCAGTACGCAGAAGTGCt                             <  2:1011607/76‑1 (MQ=255)                                              ccgaCGCCAACCTTCTCCGCTTGCTTCGACGCGGCATTCCTGTCGCTGCACCCGACTCAGTACGCAGAAGTGCTgg                           <  2:1435055/76‑1 (MQ=255)                                                gaCGCCAACCTTCTCCGCTTGCTTCGACGCGGCATTCCTGTCGCTGCACCCGACTCAGTACGCAGAAGTGCTGGTg                         <  1:1069520/76‑1 (MQ=255)                                                   gCCAACCTTCTCCGCTTGCTTCGACGCGGCATTCCTGTCGCTGCACCCGACTCAGTACGCAGAAGTGCTGGTGaaa                      <  1:681199/76‑1 (MQ=255)                                                   gCCAACCTTCTCCGCTTGCTTCGACGCGGCATTCCTGTCGCTGCACCCGACTCAGTACGCAGAAGTGCTGGTGaaa                      <  1:994130/76‑1 (MQ=255)                                                    ccAACCTTCTCCGCTTGCTTCGACGCGGCATTCCTGTCGCTGCACCCGACTCAGTACGCAGAAGTGCTGGTGAAAc                     >  2:868142/1‑76 (MQ=255)                                                     cAACCTCCTCCGCTTGCTTCGACGCGGCATTCCTGTCGCTGCACCCGACTCAGTACGCAGAAGTGCTGGTGAAACg                    <  2:441164/76‑1 (MQ=255)                                                         cTTCTCCGCTTGCTTCGACGCGGCATTCCTGTCGCTGCACCCGACTCAGTACGCAGAAGTGCTGGTGAAACGTATg                >  2:1375172/1‑76 (MQ=255)                                                          ttCTCCGCTTGCTTCGACGCGGCATTCCTGTCGCTGCACCCGACTCAGTACGCAGAAGTGCTGGTGAAACGTATGc               <  2:691512/76‑1 (MQ=255)                                                          ttCTCCGCTTGCTTCGACGCGGCATTCCTGTCGCTGCACCCGACTCAGTACGCAGAAGTGCTGGTGAAACGTATGc               <  2:1160126/76‑1 (MQ=255)                                                            ctcCGCTTGCTTCGACGCGGCATTCCTGTCGCTGCACCCGACTCAGTACGCAGAAGTGCTGGTGAAACGTATGCAg             >  2:1484930/1‑76 (MQ=255)                                                              ccGCTTGCTTCGACGCGGCATTCCTGTCGCTGCACCCGACTCAGTACGCAGAAGTGCTGGTGAAACGTATGCAggc           >  2:771326/1‑76 (MQ=255)                                                                   tgcttCGACGCGGCATTCCTGTCGCTGCACCCGACTCAGTACGCAGAAGCGCTGGTGAAACGTATGCAGGaggcgg      >  1:1598866/1‑76 (MQ=255)                                                                       tCGACGCGGCATTCCTGTCGCTGCACCCGACTCAGTACGCAGAAGTGCTGGTGAAACGTATGCAGGCGGCGGgcgc  >  2:1230763/1‑76 (MQ=255)                                                                       tCGACGCGGCATTCCTGTCGCTGCACCCGACTCAGTACGCAGAAGTGCTGGTGAAACGTATGCAGGCGGCGGgcgc  >  1:603814/1‑76 (MQ=255)                                                                          |                                                                         GGCTTCACCGCCAAACTGGCCGGTACTGAGCGTGGCATCACCGAACCGACGCCAACCTTCTCCGCTTGCTTCGGCGCGGCATTCCTGTCGCTGCACCCGACTCAGTACGCAGAAGTGCTGGTGAAACGTATGCAGGCGGCGGGCGC  >  REL606/3462657‑3462802

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