BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	REL606	3,292,285	A→G	K807E (AAG→GAG)	gltB →	glutamate synthase, large subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	REL606	3,292,285	0	A	G	100.0%	219.7 / NA	56	K807E (AAG→GAG)	gltB	glutamate synthase, large subunit
Reads supporting (aligned to +/- strand): ref base A (0/0); new base G (27/29); total (27/29)

CGCATTGGTGGGGCAAGCTTTGAAGACTTCCAGCAGGATCTGCTGAATCTGTCGAAACGTGCCTGGCTGGCGCGTAAGCCCATCAGCCAGGGCGGCCTGCTGAAATACGTCCACGGCGGCGAATACCACGCCTACAACCCGGACGT  >  REL606/3292210‑3292355

cGCATTGGTGGGGCAAGCTTTGAAGACTTCCAGCAGGATCTGCTGAATCTGTCGAAACGTGCCTGGCTGGCGCGTa                                                                        <  1:20674/76‑1 (MQ=255)

        tggGGCAAGCTTTGAAGACTTCCAGCAGGATCTGCTGAATCTGTCGAAACGTGCCTGGCTGGCGCGTGAGCCCATc                                                                >  2:611267/1‑76 (MQ=255)

            gcgaGCTTTGAAGTCTGCCAGAAGGATCTGCTGAATCTGTCGAAACGTGCCTGGCTGGCGCGGGAGCCCATCAGcc                                                            <  2:468506/73‑1 (MQ=255)

            gCAAGCTTTGAAGACTTCCAGCAGGATCTGCTGAATCTGTCGAAACGTGCCTGGCTGGCGCGTGAGCCCATCAGcc                                                            <  2:458776/76‑1 (MQ=255)

            gCAAGCTTTGAAGACTTCCAGCAGGATCTGCTGAATCTGTCGAAACGTGCCTGGCTGGCGCGTGAGCCCATCAGcc                                                            <  2:978499/76‑1 (MQ=255)

            gCAAGCTTTGAAGACTTCCAGCAGGATCTGCTGAATCTGTCGAAACGTGCCTGGCTGGCGCGTGAGCCCATCAGcc                                                            <  2:935892/76‑1 (MQ=255)

            gCAAGCTTTGAAGACTTCCAGCAGGATCTGCTGAATCTGTCGAAACGTGCCTGGCTGGCGCGTGAGCCCATCAGcc                                                            <  2:46897/76‑1 (MQ=255)

             cAAGCTTTGAAGACTTCCAGCAGGATCTGCTGAATCTGTCGAAACGTGCCTGGCTGGCGCGTGAGCCCATCAGCCa                                                           <  2:28874/76‑1 (MQ=255)

               aGCTTTGAAGACTTCCAGCAGGATCTGCTGAATCTGTCGAAACGTGCCTGGCTGGCGCGTGAGCCCATCAGCCAgg                                                         <  1:598284/76‑1 (MQ=255)

                   ttGAAGACTTCCAGCAGGATCTGCTGAATCTGTCGAAACGTGCCTGGCTGGCGCGTGAGCCCATCAGCCAcggcgg                                                     >  2:1490600/1‑76 (MQ=255)

                   ttGAAGACTTCCAGCAGGATCTGCTGAATCTGTCGAAACGTGCCTGGCTGGCGCGTGAGCCCATCAGCCAGggcgg                                                     <  1:1324525/76‑1 (MQ=255)

                     gAAGACTTCCAGCAGGATCTGCTGAATCTGTCGAAACGTGCCTGGCTGGCGCGTGAGCCCATCAGCCAGGGCGGcc                                                   >  1:551117/1‑76 (MQ=255)

                       aGACTTCCAGCAGGATCTGCTGAATCTGTCGAAACGTGCCTGGCTGGCGCGTGAGCCCATCAGCCAGGGCGGCctg                                                 <  2:461385/76‑1 (MQ=255)

                       aGACTTCCACCAGGATTTGCTAAATCTGTCGAAGCGTGCCTGGCTGGCGCGTGAGCCCATCAGCCAGGGCGGCctg                                                 <  2:531315/76‑1 (MQ=255)

                            tCCAGCAGGATCTGCTGAATCTGTCGAAACGTGCCTGGCTGGCGCGTGAGCCCATCAGCCAGGGCGGCCTGCTGaa                                            <  2:1363411/76‑1 (MQ=255)

                              cagcagGATCTGCTGAATCTGTCGAAACGTGCCTGGCTGGCGCGTGAGCCCATCAGCCAGGGCGGCCTGCTGAAAt                                          <  1:998816/76‑1 (MQ=255)

                               agcagGATCTGCTGAATCTGTCGAAACGTGCCTGGCTGGCGCGTGAGCCCATCAGCCAGGGCGGCCTGCTGAAATa                                         >  1:10567/1‑76 (MQ=255)

                               agcagGATCTGCTGAATCTGTCGAAACGTGCCTGGCTGGCGCGTGAGCCCATCAGCCAGGGCGGCCTGCTGAAATa                                         <  1:1120202/76‑1 (MQ=255)

                                 cagGATCTGCTGAATCTGTCGAAACGTGCCTGGCTGGCGCGTGAGCCCATCAGCCAGGGCGGCCTGCTGAAATACg                                       >  1:917937/1‑76 (MQ=255)

                                    gATCTGCTGAATCTGTCGAAACGTGCCTGGCTGGCGCGTGAGCCCATCAGCCAGGGCGGCCTGCTGAAATACGTcc                                    >  2:1384523/1‑76 (MQ=255)

                                     aTCTGCTGAATCTGTCGAAACGTGCCTGGCTGGCGCGTGAGCCCATCAGCCAGGGCGGCCTGCTGAAATACGTCCa                                   <  1:388541/76‑1 (MQ=255)

                                       ctgctgAATCTGTCGAAACGTGCCTGGCTGGCGCGTGAGCCCATCAGCCAGGGCGGCCTGCTGAAATACGTCCAcg                                 <  1:1210907/76‑1 (MQ=255)

                                       ctgctgAATCTGTCGAAACGTGCCTGGCTGGCGCGTGAGCCCATCAGCCAGGGCGGCCTGCTGAAATACGTCCAcg                                 >  2:1511064/1‑76 (MQ=255)

                                        tgctgAATCTGTCGAAACGTGCCTGGCTGGCGCGTGAGCCCATCAGCCAGGGCGGCCTGCTGAAATACGTCCAcgg                                >  2:489371/1‑76 (MQ=255)

                                        tgctgAATCTGTCGAAACGTGCCTGGCTGGCGCGTGAGCCCATCAGCCAGGGCGGCCTGCTGAAATACGTCCAcgg                                <  1:1213893/76‑1 (MQ=255)

                                         gctgAATCTGTCGAAACGTGCCTGGCTGGCGCGTGAGCCCATCAGCCAGGGCGGCCTGCTGAAATACGTCCAcggc                               >  2:160885/1‑76 (MQ=255)

                                           tgAATCTGTCGAAACGTGCCTGGCTGGCGCGTGAGCCCATCAGCCAGGGCGGCCTGCTGAAATACGTCCAcggcgg                             >  2:649728/1‑76 (MQ=255)

                                            gAATCTGTCGAAGCGTGCCTGGCTGGCGCGTGAGCCCATCAGCCAGGGCGGCCTGCTGAAATACGTCCACggcggc                            <  2:848375/76‑1 (MQ=255)

                                            gAATCTGTCGAAACGTGCCTGGCTGGCGCGTGAGCTCAACAGCCAGGGCGGCCTGCTGAAATACGTCCACggcggc                            >  1:596188/1‑76 (MQ=255)

                                            gAATCTGTCGAAACGTGCCTGGCTGGCGCGTGAGCCCATCAGCCAGGGCGGCCTGCTGAAATACGTCCACggcggc                            <  1:1071882/76‑1 (MQ=255)

                                            gAATCTGTCGAAACGTGCCTGGCTGGCGCGTGAGCCCATCAGCCAGGGCGGCCTGCTGAAATACGTCCACggcggc                            >  1:686885/1‑76 (MQ=255)

                                             aaTCTGTCGAAACGTGCCTGGCTGGCGCGTGAGCCCATCAGCCAGGGCGGCCTGCTGAAATACGTCCACGgcggcg                           >  2:455025/1‑76 (MQ=255)

                                                cTGTCGAAACGTGCCTGGCTGGCGCGTGAGCCCATCAGCCAGGGCGGCCTGCTGAAATACGTCCACGGCGGCGAAt                        >  1:1057061/1‑76 (MQ=255)

                                                 tGTCGAAACGGGCCTGGCTGGCGCGTGAGCCCATCAGCCAGGGCGGCCTGCTGAAATACGTCCACGGCGGCGAATa                       <  1:1468367/76‑1 (MQ=255)

                                                  gTCGAAACGTGCCTGGCTGGCGCGTGAGCCCATCAGCCAGGGCGGCCTGCTGAAATACGTCCACGGCGGGGAATAc                      >  1:1436211/1‑76 (MQ=255)

                                                  gTCGAAACGTGCCTGGCTGGCGCGTGAGCCCATCAGCCAGGGCGGCCTGCTGAAATACGTCCACGGCGGCGAATaa                      >  2:183165/1‑75 (MQ=255)

                                                    cGAAACGTGCCTGGCTGGCGCGTGAGCCCATCAGCCAGGGCGGCCTGCTGAAATACGTCCACGGCGGCGAATACCa                    <  2:1058283/76‑1 (MQ=255)

                                                    cGAAACGTGCCTGGCTGGCGCGTGAGCCCATCAGCCAGGGCGGCCTGCTGAAATACGTCCACGGCGGCGAATACCa                    <  1:1506352/76‑1 (MQ=255)

                                                    cGAAACGTGCCTGGCTGGCGCGTGAGCCCATCAGCCAGGGCGGCCTGCTGAAATACGTCCACGGCGGCGAATACCa                    >  1:1279715/1‑76 (MQ=255)

                                                     gAAACGTGCCTGGCTGGCGCGTGAGCCCATCAGCCAGGGCGGCCTGCTGAAATACGTCCACGGCGGCGAATACCAc                   >  2:1338746/1‑76 (MQ=255)

                                                     gAAACGTGCCTGGCTGGCGCGTGAGCCCATCAGCCAGGGCGGCCTGCTGAAATACGTCCACGGCGGCGAATACCAc                   <  1:1412429/76‑1 (MQ=255)

                                                      aaaCGTGCCTGGCTGGCGCGTGAGCCCATCAGCCAGGGCGGCCTGCTGAAATACGTCCACGGCGGCGAATACCaag                  >  2:869342/1‑74 (MQ=255)

                                                      aaaCGTGCCTGGCTGGCGCGTGAGCCCATCAGCCAGGGCGGCCTGCTGAAATACGTCCACGGCGGCGAATACCACg                  >  2:888879/1‑76 (MQ=255)

                                                      aaaCGTGCCTGGCTGGCGCGTGAGCCCATCAGCCAGGGCGGCCTGCTGAAATACGTCCACGGCGGCGAATACCACg                  >  2:547171/1‑76 (MQ=255)

                                                      aaaCGTGCCTGGCTGGCGCGTGAGCCCATCAGCCAGGGCGGCCTGCTGAAATACGTCCACGGCGGCGAATACCACg                  <  1:235655/76‑1 (MQ=255)

                                                       aaCGTGCCTGGCTGGCGCGTGAGCCCATCAGCCAGGGCGGCCTGCTGAAATACGTCCACGGCGGCGAATACCACGc                 >  2:1040023/1‑76 (MQ=255)

                                                         cGTGCCTGGCTGGCGCGTGAGCCCATCCGCCAGGGCGGCCTGCTGAAATACGTCCACGGCGGCGAATACCACGCCt               <  2:128471/76‑1 (MQ=255)

                                                         cGTGCCTGGCTGGCGCGTGAGCCCATCAGCCAGGGCGGCCTGCTGAAATACGTCCACGGCGGCGAATACCACGCCt               <  2:1573286/76‑1 (MQ=255)

                                                             ccTGGCTGGCGCGTGAGCCCATCAGCCAGGGCGGCCTGCTGAAATACGTCCACGGCGGCGAATACCACGCCTACaa           <  2:306783/76‑1 (MQ=255)

                                                              ctggctggCGCGTGAGCCCATCAGCCAGGGCGGCCTGCTGAAATACGTCCACGGCGGCGAATACCACGCCTACAAc          <  2:686884/76‑1 (MQ=255)

                                                               tggctggcGCGTGAGCCCATCAGCCAGGGCGGCCTGCTGAAATACGTCCACGGCGGCGAATACCACGCCTACAAcc         >  1:1214972/1‑76 (MQ=255)

                                                                ggctggcGCGTGAGCCCATCAGCCAGGGCGGCCTGCTGAAATACGTCCACGGCGGCGAATACCACGGCCACAAacc        >  2:799207/1‑76 (MQ=255)

                                                                ggctggcGCGTGAGCCCATCAGCCAGGGCGGCCTGCTGAAATACGTCCACGGCGGCGAATACCACGCCTACAAccc        >  1:562152/1‑76 (MQ=255)

                                                                ggctggcGCGTGAGCCCATCAGCCAGGGCGGCCTGCTGAAATACGTCCACGGCGGCGAATACCACGCCTACAAccc        >  2:1091795/1‑76 (MQ=255)

                                                                  ctggcGCGTGAGCCCATCAGCCAGGGCGGCCTGCTGAAATACGTCCACGGCGGCGAATACCACGCCTACAACCCgg      <  1:1421802/76‑1 (MQ=255)

                                                                  ctggcGCGTGAGCCCATCAGCCAGGGCGGCCTGCTGAAATACGTCCACGGCGGCGAATACCACGCCTACAACCCgg      >  2:10329/1‑76 (MQ=255)

                                                                  ctggcGCGTGAGCCCATCAGCCAGGGCGGCCTGCTGAAATACGTCCACGGCGGCGAATACCACGCCTACAACCCgg      >  2:319806/1‑76 (MQ=255)

                                                                    ggcGCGTGAGCCCATCAGCCAGGGCGGCCTGCTGAAATACGTCCACGGCGGCGAATACCACGCCTACAACCCGGAc    <  1:455026/76‑1 (MQ=255)

                                                                      cgcgTGAGCCCATCAGCCAGGGCGGCCTGCTGAAATACGTCCACGGCGGCGAATACCACGCCTACAACCCGGACgt  <  1:1551642/76‑1 (MQ=255)

CGCATTGGTGGGGCAAGCTTTGAAGACTTCCAGCAGGATCTGCTGAATCTGTCGAAACGTGCCTGGCTGGCGCGTAAGCCCATCAGCCAGGGCGGCCTGCTGAAATACGTCCACGGCGGCGAATACCACGCCTACAACCCGGACGT  >  REL606/3292210‑3292355

Alignment Legend
Aligned base mismatch/match (shaded by quality score): `ATCG/ATCG` `< 3 ≤` `ATCG/ATCG` `< 20 ≤` `ATCG/ATCG` `< 27 ≤` `ATCG/ATCG` `< 33 ≤` `ATCG/ATCG` `< 38 ≤` `ATCG/ATCG`
Unaligned base: `atcg` Masked matching base: `atcg` Alignment gap: `‑` Deleted base: `‑`