Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	REL606	4,421,946	T→A	D66E (GAT→GAA)	cysQ →	PAPS (adenosine 3'‑phosphate 5'‑phosphosulfate) 3'(2'),5'‑bisphosphate nucleotidase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	REL606	4,421,946	0	T	A	100.0%	117.4 / NA	33	D66E (GAT→GAA)	cysQ	PAPS (adenosine 3'‑phosphate 5'‑phosphosulfate) 3'(2'),5'‑bisphosphate nucleotidase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base A (14/19);  total (14/19)

TGCCGCTCACACCGTTATCATGGACGGTTTACGTACGCTGACACCGGATATTCCGGTCCTTTCTGAAGAAGATCCTCCCGGTTGGGAAGTCCGTCAGCACTGGCAGCGTTACTGGCTGGTAGACCCGCTGGATGGTACTAAAGA  >  REL606/4421874‑4422017                                                                         |                                                                        tGCCGCTCACACCGTTATCATGGACGGTTTACGTACGCTGACACCGGATATTCCGGTCCTTTCTGAAGAAGAacct                                                                      <  2:500164/76‑1 (MQ=255)   ccGCTCACACCGTTATCATGGACGGTTTACGTACGCTGACACCGGATATTCCGGTCCTTTCTGAAGAAGAacctcc                                                                    <  1:536218/76‑1 (MQ=255)     gCTCACTCCGTTATCATGGACGGTTTACGTACGCTGACACCGGATATTCCGGTCCTTTCTGAAGAAGAACCTCCCg                                                                  <  2:1176303/76‑1 (MQ=255)          cacCGTTATCATGGACGGTTTACGTACGCTGACACCGGATATTCCGGTCCTTTCTGAAGAAGAACCTCCCGGTTgg                                                             >  2:402305/1‑76 (MQ=255)           acCGTTATNATGGACGGTTTACGTACGCTGACACCGGATATTCCGGTCCTTTCTGAAGAAGAACCTCCCGGTTggg                                                            <  1:1064985/76‑1 (MQ=255)           acCGTTATCATGGACGGTTTACGTACGCTGACACCGGATATTCCGGTCCTTTCTGAAGAAGAACCTCCCGGTTggg                                                            <  2:1069253/76‑1 (MQ=255)           acCGTTATCATGGACGGTTTACGTACGCTGACACCGGATATTCCGGTCCTTTCTGAAGAAGAACCTCCCGGTTggg                                                            >  1:946169/1‑76 (MQ=255)           acCGTTATCATGGACGGTTTACGTACGCTGACACCGGATATTCCGGTCCTTTCTGAAGAAGAACCTCCCGGTTggg                                                            >  2:834624/1‑76 (MQ=255)              gTTATCATGGACGGTTTACGTACGCTCACACCGGATATTCCGGTCCTTTCTGAAGAAGAACCTCCCGGTTGGGAAg                                                         >  2:580469/1‑76 (MQ=255)                 aTCATGGACGGTTTACGTACGCTGACACCGGANATTCCGGTCCTTTCTGAAGAAGAACCTCCCGGTTGGGAANGcc                                                      >  2:1133585/1‑76 (MQ=255)                  tcctGGACGGATTACGTACGCTGACACCGGATATTCCGGTCCTTTCTGACGAAGAACCTCCCGGTTGGACAGTCCg                                                     >  2:524716/4‑76 (MQ=255)                    aTGGACGGTTTACGTACGCTGACACCGGATATTCCGGTCCTTTCTGAAGAAGAACCTCCCGGTTGGGAAGTCCGTc                                                   <  2:907190/76‑1 (MQ=255)                     tGGACGGTTTACGTACGCTGACACCGGATATTCCGGTCCTTTCTGAAGAAGAACCTCCCGGTTGGGAAGTCCGTCa                                                  <  2:810668/76‑1 (MQ=255)                       gACGGTTTACGTACGCTGACACCGGATATTCCGGTCCTTTCTGAAGAAGAACCTCCCGGTTGGGAAGTCCGTCAGc                                                >  2:1172523/1‑76 (MQ=255)                         cGGTTTACGTACGCTGACACCGGATATTCCGGTCCTTTCTGAAGAAGAACCTCCCGGTTGGGAAGTCCGTCAGCAc                                              >  2:854591/1‑76 (MQ=255)                             ttACGTACGCTGACACCGGATATTCCGGTCCTTTCTGAAGAAGAACCTCCCGGTTGGGAAGTCCGTCAGCACTGGc                                          <  2:490039/76‑1 (MQ=255)                             ttACGTACGCTGACACCGGATATTCCGGTCCTTTCTGAAGAAGAACCTCCCGGTTGGGAAGTCCGTCAGCACTGGc                                          <  2:490056/76‑1 (MQ=255)                                      cTGACACCNGATATTCCGGNCCTTTCTGAAGAAGAACCTCCCNGTTGGGAAGTCCGTCAGCACTGGCAGCGTTAct                                 <  2:1008463/76‑1 (MQ=255)                                      cTGACACCGGATATTCCGGTCCTTTCTGAAGAAGAACCTCCCGGTTGGGAAGTCCGTCAGCACTGGCAGCGTTAct                                 <  2:541885/76‑1 (MQ=255)                                       tGACACCGGATATTCCGGTCCTTTCTGAAGAAGAACCTCCCGGTTGGGAAGTCCGTCAGCACTGGCAGCGTTActg                                >  1:319772/1‑76 (MQ=255)                                       tGACACCGGATATTCCGGTCCTTTCTGAAGAAGAACCTCCCGGTTGGGAAGTCCGTCAGCACTGGCAGCGGTActg                                >  2:930125/1‑76 (MQ=255)                                           acCGGATATTCCGGTCCTTTCTGAAGAAGAACCTCCCGGTTGGGAAGTCCGTCAGCACTGGCAGCGTTActggctg                            >  2:976606/1‑76 (MQ=255)                                           acCGGATATTCCGGTCCTTTCTGAAGAAGAACCTCCCGGTTGGGAAGTCCGTCAGCACTGGCAGCGTTActggctg                            >  2:53904/1‑76 (MQ=255)                                             cGGATATTCCGGTCCTTTCTGAAGAAGAACCTCCCGGTTGGGAAGTCCGTCAGCACTGGCAGCNTTACTGGCTGGt                          <  1:13465/76‑1 (MQ=255)                                                 tatTCCGGTCCTTTCTGAAGAAGAACCTCCCGGTTGGGAAGTCCGTCAGCACTGGCAGCGTTACTGGCTGGGAGAc                      >  2:576910/1‑76 (MQ=255)                                                    tCCGGTCCTTTCTGAAGAAGAACTTCCCGGTTGGGAAGTCCGTCAGCACTGGCAGCGTTACTGGCTGGTAGACCCg                   <  1:1043880/76‑1 (MQ=255)                                                       ggTCCTTTCTGAAGAAGAACCTCCCGGTTGGGAAGTCCGTCAGCACTGGCAGCGTTACTGGCTGGTAGACCCGCTg                <  2:1066087/76‑1 (MQ=255)                                                          ccTTTCTGAAGAAGAACCTCCCGGTTGGGAAGTCCGTCAGCACTGGCAGCGTTACTGGCTGGTAGACCCGCTGGAt             <  2:642927/76‑1 (MQ=255)                                                           cTTTCTGAAGAAGAACCTCCCGGTTGGGAAGTCCGTCAGCACTGGCAGCGTTACTGGCTGGTAGACCCGCTGGATg            <  2:27114/76‑1 (MQ=255)                                                           cTTTCTGAAGAAGAACCTCCCGGTTGGGAAGTCCGTCAGCACTGGCAGCGTTACTGGCTGGTAGACCCGCTGGATg            <  2:726200/76‑1 (MQ=255)                                                            tttCTGAAGAAGAACCTCCCGGTTGGGAAGTCCGTCAGCACTGGCAGCGTTACTGGCTGGGAGACCCGCTGGGTgg           >  2:447120/1‑76 (MQ=255)                                                              tCTGAAGAAGAACCTCCCGGTTGGGAAGTCCGTCAGCACTGGCAGCGTTACTGGCTGGTAGACCCGCTGGATGGTa         <  2:428138/76‑1 (MQ=255)                                                              tCTGAAGAAGAACCTCCCGGTTGGGAAGTCCGTCAGCACTGGCAGCGTTACTGGCTGGTAGACCCGCTGGATGGTa         <  1:854592/76‑1 (MQ=255)                                                               cTGAAGAAGAACCTCCCGGTTGGGAAGTCCGTCAGCACTGGCAGCGTTACTGGCTGGTAGACCCGCTGGATGGTAc        <  1:376586/76‑1 (MQ=255)                                                                   agaagaACCTCCCGGTTGGGAAGTCCGTCAGCACTGGCAGCGTTACTGGCTGGTAGACCCGCTGGATGGTACTaaa    <  1:516879/76‑1 (MQ=255)                                                                     aagaACCTCCCGGTTGGGAAGTCCGTCAGCACTGGCAGCGTTACTGGCTGGGAGACCCGCTGGGTGGTACTAAaga  >  2:9170/1‑76 (MQ=255)                                                                         |                                                                        TGCCGCTCACACCGTTATCATGGACGGTTTACGTACGCTGACACCGGATATTCCGGTCCTTTCTGAAGAAGATCCTCCCGGTTGGGAAGTCCGTCAGCACTGGCAGCGTTACTGGCTGGTAGACCCGCTGGATGGTACTAAAGA  >  REL606/4421874‑4422017

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