Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	REL606	3,271,597	G→A	P45S (CCG→TCG)	murA ←	UDP‑N‑acetylglucosamine 1‑carboxyvinyltransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	REL606	3,271,597	0	G	A	100.0%	115.1 / NA	31	P45S (CCG→TCG)	murA	UDP‑N‑acetylglucosamine 1‑carboxyvinyltransferase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (13/18);  total (13/18)

GAACCATTACGTTCTACTTTCGCACCCAGCTGGCTTAGCAGCTTCATTGATGTATCGACGTCTTTCAGTTTCGGGACGTTCTGGATCTCTACCGGTTCTTCCGCCAGTAGAGCGGCAAAAAGGATAGGCAGAGCAGCATTTTTAG  >  REL606/3271524‑3271668                                                                          |                                                                        gAACCATTAAGTGCTACTTTCGCGCCCAGCTGGCTTAGCAGCTTCATTGATGTATCGANGTCTTTCAGTTTCgaga                                                                       <  1:153240/76‑1 (MQ=255)  aaCCATTTCGATCTACTTGCGCACACATCTGGCCTAGCTNNTTCATTGNTGTATCGACGTCTTTGAGTTTCGAGAc                                                                      <  1:1248236/76‑1 (MQ=255)    ccATTACGTTCTACTTTCGCACCCAGCTGGCTTAGCAGCTTCATTGATGTATCGACGTCTTTCAGTTTCGAGACGt                                                                    >  2:64417/1‑76 (MQ=255)      aTTACGTTCTACTTTCGCACCCAGCTGGCTTAGCAGCTTCANTGATGTATCGACGTCTTTNAGNTTCGAGACGTTc                                                                  >  1:1039100/1‑76 (MQ=255)          cGTTCTTCTATCGCACCGTGCTGGCTTAGCAGCTTCATTGATGTATCGAGGTCTTTCAGTANCGNGACGTTCTGGa                                                              >  2:195350/1‑76 (MQ=255)                  tttCGCACCCAGCTGGCTTAGCAGCTTCATTGATGTATCGACGTCTTTCAGTTTCGAGACGTTCTGGATCTCTAcc                                                      >  2:799721/1‑76 (MQ=255)                    tCGCACCCAGCTGGCTTAGCAGCTTCATTGATGTATCGACGTCTTTCAGTTTCGAGACGTTCTGGATCTCTACCgg                                                    >  2:407342/1‑76 (MQ=255)                    tCGCACCCAGCTGGCTTAGCAGCTTCATTGATGTATCGACGTCTTTCAGTTTCGAGACGTTCTGGATCTCTACCgg                                                    >  2:796678/1‑76 (MQ=255)                      gCACCCAGCTGGCTTAGCAGCTTCACTGATGTATCGACGTCTTTCAGTTTCGAGACGTTCTGGATCTCTACCGGtt                                                  >  2:905847/1‑76 (MQ=255)                        ccccAGCTGGCTTAGCAGCTTCATTGATGTATCGACGTCTTTCAGTTTCGAGACGTTCTGGATCTCTACCGGttct                                                <  2:112504/75‑1 (MQ=255)                        aCCCAGCTGGCTTAGCAGCTTCATTGATGTATCGACGTCTTTCAGTTTCGAGACGTTCTGGATCTCTACCGGttct                                                >  2:888447/1‑76 (MQ=255)                           cAGCTGGCTTAGCAGCTTCATTGATGTATCGACGTCTTTCAGTTTCGAGACGTTCTGGATCTCTACCGGTTCTTcc                                             >  1:556895/1‑76 (MQ=255)                            aGCTGGCTTAGCAGCTTCATTGATGTATCGACGTCTTTCAGTTTCGAGACGTTCTGGATCTCTACCGGTTCTTCCg                                            >  1:689340/1‑76 (MQ=255)                                ggCTTAGCAGCTTCATTGATGTATCGACGTCTTTCAGTTTCGAGACGTTCTGTATCTCTACCGGTTCTTCCGCCAg                                        <  1:510910/76‑1 (MQ=255)                                 gCTTAGCAGCTTCATTGATGTATCGACGTCTTTCAGTTTCGAGACGTTCTGGATCTCTACCGGTTCTTCCGCCAGt                                       <  2:735323/76‑1 (MQ=255)                                  cTTAGCAGCTTCATTGATGTATCGACGTCTTTCAGTTTCGAGACGTTCTGGATCTCTACCGGTTCTTCCGCCAGTa                                      <  1:986332/76‑1 (MQ=255)                                   ttAGCAGCTTCATTGATGTATCGACGTCTTTCAGTTTCGAGACGTTCTGGATCTCTACCGGTTCTTCCGCCAGTag                                     <  2:960961/76‑1 (MQ=255)                                    tAGCAGCTTCATTGATGTATCGACGTCTTTCAGTTTCGAGACGTTCTGGATCTCTACCGGTTCTTCCGCCAGTaga                                    >  2:30039/1‑76 (MQ=255)                                       cagcTTCATTGATGTATCGACGTCTTTCAGTTTCGAGACGTTCTGGATCTCTACCGGTTCTTCCGCCAGTAGAGCg                                 <  1:734507/76‑1 (MQ=255)                                         gcTTCATTGATGTATCGACGTCTTTCAGTTTCGAGACGTTCTGGATCTCTACCGGTTCTTCCGCCAGTAGAGCGGc                               >  1:956398/1‑76 (MQ=255)                                          cTTCATTGATGTATCGACGTCTTTCAGTTTCGAGACGTTCTGGATCTCTACCGGTTCTTCCGCCAGTAGAGCGGCa                              <  2:1082089/76‑1 (MQ=255)                                           ttCATTGATGTATCGACGTCTTTCAGTTTCGAGACGTTCTGGATCTCTACCGGTTCTTCCGCCAGTAGAGCGGCaa                             <  2:829088/76‑1 (MQ=255)                                                  aTGTATCGACGTCTTTCAGTTTCGAGACGTTCTGGATCTCTACCGGTTCTTCCGCCAGTAGAGCGGCAAAAAGGAt                      <  1:744030/76‑1 (MQ=255)                                                   tGTATCGACGTCTTTCAGTTTCGAGACGTTCTGGATCTCTACCGGTTCTTCCGCCAGTAGAGCGGCAAAAAGGATa                     <  1:389881/76‑1 (MQ=255)                                                        cGACGTCTTTCAGTTTCGAGACGTTCTGGATCTCTACCGGTTCTTCCGCCAGTAGAGCGGCAAAAAGGATAGGCag                <  2:761473/76‑1 (MQ=255)                                                         gACGTCTTTCAGTTTCGAGACGTTCTGGATCTCTACCGGTTCTTCCGCCAGTAGAGCGGCAAAAAGGATAGGCaga               <  1:464996/76‑1 (MQ=255)                                                          aCGTCTTTCAGTTTCGAGACGTACTGGNGCNCTNCCAGTTCNTCCGCCCGTAGAGCGGCAAAAAGGATAGGCagag              <  1:1103296/76‑1 (MQ=255)                                                              cTTTCAGTTTCGAGACGTTCTGGATCTCTACCGGTTCTTCCGCCAGTAGAGCGGCAAAAAGGATAGGCAGagcagc          <  2:1008864/76‑1 (MQ=255)                                                              cTTTCAGTTTCGAGACGTTCTGGATCTCTACCGGTTCTTCCGCCAGTAGAGCGGCAAAAAGGATAGGCAGagcagc          >  2:487117/1‑76 (MQ=255)                                                                ttCAGTTTCGAGACGTTCTGGATCTCTACCGGTTCTTCCGCCAGTAGAGCGGCAAAAAGGATAGGCAGAGCAGCAt        <  2:239659/76‑1 (MQ=255)                                                                ttCAGTTTCGAGACGTTCTGGATCTCTACCGGTTCTTCCGCCAGTAGAGCGGCAAAAAGGATAGGCAGAGCAGCAt        >  1:924913/1‑76 (MQ=255)                                                                 tCAGTTTCGAGACGTGCTGGATCTCTGCCGGTTCGTCCGCCAGTAGAGCGGCAAAAAGGATAGGCAGAGCAGCAtt       <  2:1128022/76‑1 (MQ=255)                                                                     tttCGAGACGTTCTGGATCTCTACCGGTTCTTCCGCCAGTAGAGCGGCAAAAAGGATAGGCAGAGCAGCATTTTTa   <  1:888448/76‑1 (MQ=255)                                                                      ttCGAGACGTTCTGGATCTCTACCGGTTCTTCCGCCAGTAGAGCGGCAAAAAGGATAGGCAGAGCAGCATTTTTAg  <  2:847587/76‑1 (MQ=255)                                                                          |                                                                        GAACCATTACGTTCTACTTTCGCACCCAGCTGGCTTAGCAGCTTCATTGATGTATCGACGTCTTTCAGTTTCGGGACGTTCTGGATCTCTACCGGTTCTTCCGCCAGTAGAGCGGCAAAAAGGATAGGCAGAGCAGCATTTTTAG  >  REL606/3271524‑3271668

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