Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	REL606	4,438,531	T→G	D110A (GAC→GCC)	fbp ←	fructose‑1,6‑bisphosphatase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	REL606	4,438,531	0	T	G	100.0%	141.6 / NA	35	D110A (GAC→GCC)	fbp	fructose‑1,6‑bisphosphatase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base G (13/22);  total (13/22)

CGCGGCGGTAGATGGAGAAAATGGTACCGACAGAGACGTTAACATCGATGTTGGACGAGCCATCCAGGGGGTCCATCAGCACCACGTATTTTGCGTGTTCACAGCCTTCAAAGACGACAATCTCATCTTCTTCTTCAGAGGC  >  REL606/4438460‑4438601                                                                        |                                                                       cgcgGCGGTAGATGGAGAAAATGGTACCGACAGAGACGTTAACATCGATGTTGGACGAGCCATCCAGGGGGGCCAt                                                                    <  2:367096/76‑1 (MQ=255)  gcggcgGTAGATGGAGAAAATGGTACCGACAGAGACGTTAACATCGATGTTGGACGAGCCATCCAGGGGGGCCATc                                                                   <  2:1371033/76‑1 (MQ=255)  gcggcgGTAGATGGAGAAAATGGTACCGACAGAGACGTTAACATCGATGTTGGACGAGCCATCCAGGGGGGCCATc                                                                   <  2:219529/76‑1 (MQ=255)   cggcggTAGATGGAGAAAATGGTACCGACAGAGACGTTAACATCGATGTTGGACGAGCCATCCAGGGGGGCCATCa                                                                  <  1:1130560/76‑1 (MQ=255)       ggTAGATGGAGAAAATGGTACCGACAGAGACGTTAACATCGATGTTGGACGAGCCATCCAGGGGGGCCATCAGcac                                                              <  1:1147888/76‑1 (MQ=255)         tAGATGGAGAAAATGGTACCGACAGAGACGTTAACATCGATGTTGGACGAGCCATCCAGGGGGGCCATCAGcacca                                                            >  1:1196417/1‑76 (MQ=255)          ngATGGAGAAAATGGTACCGACAGAGACGTTAACATCGATGTTGGACGAGCCATCCAGGGGGGCCATCAGcaccac                                                           >  1:596352/2‑76 (MQ=255)                  aaaaTGGTACCGACAGAGACGTTAACATCGATGTTGGACGAGCCATCCAGGGGGGCCATCAGCACCACGTATTTTg                                                   >  2:119870/1‑76 (MQ=255)                    aaTGGTACCGACAGAGACGTTAACATCGATGTTGGACGAGCCATCCAGGGGGGCCATCAGCACCACGTATTTTGCg                                                 >  2:498976/1‑76 (MQ=255)                    aaTGGTACCGACAGAGACGTTAACATCGATGTTGGACGAGCCATCCAGGGGGGCCATCAGCACCACGTATTTTGCg                                                 >  2:1190087/1‑76 (MQ=255)                     aTGGTACCGACAGAGACGTTAACATCGATGTTGGACGAGCCATCCAGGGGGGCCATCAGCACCACGTATTTTGCgt                                                >  1:1457843/1‑76 (MQ=255)                     aTGGTACCGACAGAGACGTTAACATCGATGTTGGACGAGCCATCCAGGGGGGCCATCAGCACCACGTATTTTGCgt                                                >  1:955896/1‑76 (MQ=255)                       ggTACCGACAGAGACGTTAACATCGATGTTGGACGAGCCATCCAGGGGGGCCATCAGCACCACGTATTTTGCgtgt                                              >  1:1421066/1‑76 (MQ=255)                       ggTACCGACAGAGACGTTAACATCGATGTTGGACGAGCCATCCAGGGGGGCCATCAGCACCACGTATTTTGCgtgt                                              >  1:1059577/1‑76 (MQ=255)                        gTACCGACAGAGACGTTAACATCGATGTTGGACGAGCCATCCAGGGGGGCCATCAGCACCACGTATTTTGCGTGtt                                             <  2:779322/76‑1 (MQ=255)                           ccGACAGAGACGTTAACATCGATGTTGGACGAGCCATCCAGGGGGGCCATCAGCACCACGTATTTTGCGTGTTcac                                          <  1:202662/76‑1 (MQ=255)                                agagACGTTAACATCGATGTTGGACGAGCCATCCAGGGGGGCCATCAGCACCACGTATTTTGCGTGTTCACAGCCt                                     <  1:673804/76‑1 (MQ=255)                                agagACGTTAACATCGATGTTGGACGAGCCATCCAGGGGGGCCATCAGCACCACGTATTTTGCGTGTTCACAGCCt                                     >  1:236929/1‑76 (MQ=255)                                    aCGTTAACATCGATGTTGGACGAGCCATCCAGGGGGGCCATCAGCACCACGTATTTTGCGTGTTCACAGCCTTCaa                                 <  2:303140/76‑1 (MQ=255)                                      gTTAACATCGATGTTGGACGAGCCATCCAGGGGGGCCATCAGCACCACGTATTTTGCGTGTTCACAGCCTTCAAAg                               <  1:359715/76‑1 (MQ=255)                                         aaCATCGATGTTGGACGAGCCATCCAGGGGGGCCATCAGCACCACGTATTTTGCGTGTTCACAGCCTTCAAAgacg                            <  1:1404848/76‑1 (MQ=255)                                             tCGATGTTGGACGAGCCATCCAGGGGGGCCATCAGCACCACGTATTTTGCGTGTTCACAGCCTTCAAAGACGACaa                        >  2:394432/1‑76 (MQ=255)                                              cGATGTTGGACGAGCCATCCAGGGGGGCCATCAGCACCACGTATTTTGCGTGTTCACAGCCTTCAAAGACGACAAt                       <  1:1016336/76‑1 (MQ=255)                                              cGATGTTGGACGAGCCATCCAGGGGGGCCATCAGCACCACGTATTTTGCGTGTTCACAGCCTTCAAAGACGACAAt                       <  1:1068580/76‑1 (MQ=255)                                              cGATGTTGGACGAGCCATCCAGGGGGGCCATCAGCACCACGTATTTTGCGTGTTCACAGCCTTCAAAGACGACAAt                       <  1:1108854/76‑1 (MQ=255)                                              cGATGTTGGACGAGCCATCCAGGGGGGCCATCAGCACCACGTATTTTGCGTGTTCACAGCCTTCAAAGACGACAAt                       <  1:920000/76‑1 (MQ=255)                                              cGATGCTGGACGAGCCATCCAGGGGGGCCATCAGCACCACGTATTTTGCGTGTTCACAGCCTTCAAAGACGACAAt                       <  1:1193106/76‑1 (MQ=255)                                                 tGTTGGACGAGCCATCCAGGGGGGCCATCAGCACCACGTATTTTGCGTGTTCACAGCCTTCAAAGACGACAAtctc                    <  2:443683/76‑1 (MQ=255)                                                   ttGGACGAGCCATCCAGGGGGGCCATCAGCACCACGTATTTTGCGTGTTCACAGCCTTCAAAGACGACAATCTCAt                  <  1:553638/76‑1 (MQ=255)                                                   ttGGACGAGCCATCCAGGGGGGCCATCAGCACCACGTATTTTGCGTGTTCACAGCCTTCAAAGACGACAATCTCAt                  <  1:581547/76‑1 (MQ=255)                                                    tGGACGAGCCATCCAGGGGGGCCATCAGCACCACGTATTTTGCGTGTTCACAGCCTTCAAAGACGACAATCTCAtc                 >  2:814642/1‑76 (MQ=255)                                                     ggACGAGCCATCCAGGGGGGCCATCAGCACCACGTATTTTGCGTGTTCACAGCCTTCAAAGACGACAATCTCAtct                >  2:194755/1‑76 (MQ=255)                                                     ggACGAGCCATCCAGGGGGGCCATCAGCACCACGTATTTTGCGTGTTCACAGCCTTCAAAGACGACAATCTCAtct                >  2:955889/1‑76 (MQ=255)                                                      gACGAGCCATCCAGGGGGGCCATCAGCACCACGTATTTTGCGTGTTCACAGCCTTCAAAGACGACAATCTCAtctt               <  2:827145/76‑1 (MQ=255)                                                        cGAGCCATCCAGGGGGGCCATCAGCACCACGTATTTTGCGTGTTCACAGCCTTCAAAGACGACAATCTCAtcttct             <  2:1430031/76‑1 (MQ=255)                                                             cATCCAGGGGGGCCATCAGCACCACGTATTTTGCGTGTTCACAGCCTTCAAAGACGACAATCTCATCttcttcttc        <  2:509257/76‑1 (MQ=255)                                                                   ggggggCCATCAGCACCACGTATTTTGCGTGTTCACAGCCTTCAAAGACGACAATCTCATCTTCTTCTTCAGAGGc  >  1:339008/1‑76 (MQ=255)                                                                        |                                                                       CGCGGCGGTAGATGGAGAAAATGGTACCGACAGAGACGTTAACATCGATGTTGGACGAGCCATCCAGGGGGTCCATCAGCACCACGTATTTTGCGTGTTCACAGCCTTCAAAGACGACAATCTCATCTTCTTCTTCAGAGGC  >  REL606/4438460‑4438601

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