Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	1	1,623,191	A→G	intergenic (+127/‑145)	unknown → / → unknown	Diacylglycerol kinase (EC 2.7.1.107)/hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	1	1,623,191	0	A	G	100.0%	100.8 / NA	37	intergenic (+127/‑145)	unknown/unknown	Diacylglycerol kinase (EC 2.7.1.107)/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (14/23);  total (14/23)

CGCATCGGCACCGAACGCCACGAACTCTCCGGGCGTGCCAAAGATCAAGGCTCCGCTGCCGTGTTTGTCGCCTTAGTTATCGTAGTGGTCACTTGGGGATTAATGTTGTTTAACTAAAAAACGTCGGCAAAAATGACCGCACTTTGCTCTGCACCCCAAAACCTGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTATCAATTTCGTTAAAAAACTATCAAACAGGTCACCGAGCAGAATTTTGGTATCCGTGAGGTGGCTAAATTTCATCAGATTTCTCGTTCTCAAGTCATTTATTGGAAAAGCCTTTCGTGAAAGAGGGCTTAATGGCGTAAAATCCCCTTATATAAACCCTCAACGCCCTAAAATAGTG  >  1/1622948‑1623326                                                                                                                                                                                                                                                    |                                                                                                                                        cGCATCGGCACCGAACGCCACGAACTCTCCGGGCGTGCCAAAGATCAAGGCTCCGCTGCCGTGTTTGTCGCCTTAGTTATCGTAGTGGTCACTTGGGGATTAATGTTGTTTAACTAAAAAACGTCGGCAAAAATGACCGCACTTTGCTCTGCACCCCAAAACCTGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTATCAATTTCGTTAAAAAACTATCAAACAGGTCACCGAGCAGGATTTTgg                                                                                                                                  <  1:12202/251‑1 (MQ=255)   cATCGGCACCGAACGCCACGAACTCTCCGGGCGTGCCAAAGATCAAGGCTCCGCTGCCGTGTTTGTCGCCTTAGTTATCGTAGTGGTCACTTGGGGATTAATGTTGTTTAACTAAAAAACGTCGGCAAAAATGACCGCACTTTGCTCTGCACCCCAAAACCTGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTATCAATTTCGTTAAAAAACTATCAAACAGGTCACCGAGCAGGATTTTGGTa                                                                                                                                >  1:471524/1‑251 (MQ=255)                       aaCTCTCCGGGCGTGCCAAAGATCAAGGCTCCGCTGCCGTGTTTGTCGCCTTAGTTATCGTAGTGGTCACTTGGGGATTAATGTTGTTTAACTAAAAAACGTCGGCAAAAATGACCGCACTTTGCTCTGCACCCCAAAACCTGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTATCAATTTCGTTAAAAAACTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCTAAAtt                                                                                                             <  1:6153/250‑1 (MQ=255)                         ctctCCGGGCGTGCCAAAGATCAAGGCTCCGCTGCCGTGTTTGTCGCCTTAGTTATCGTAGTGGTCACTTGGGGATTAATGTTGTTTAACTAAAAAACGTCGGCAAAAATGACCGCACTTTGCTCTGCACCCCAAAACCTGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTATCAATTTCGTTAAAAAACTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCTAAATTtca                                                                                                          >  1:656662/1‑251 (MQ=255)                                                        cTGCCGTGTTTGTCGCCTTAGTTATCGTAGTGGTCACTTGGGGATTAATGTTGTTTAACTAAAAAACGTCGGCAAAAATGACCGCACTTTGCTCTGCACCCCAAAACCTGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTATCAATTTCGTTAAAAAACTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCTAAATTTCATCAGATTTCTCGTTCTCAAGTCATTTATTgg                                                                           <  1:415799/251‑1 (MQ=255)                                                          gCCGTGTTTGTCGCCTTAGTTATCGTAGTGGTCACTTGGGGATTAATGTTGTTTAACTAAAAAACGTCGGCAAAAATGACCGCACTTTGCTCTGCACCCCAAAACCTGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTATCAATTTCGTTAAAAAACTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCGaaac                                                                                                              >  1:770060/1‑213 (MQ=255)                                                          gCCGTGTTTGTCGCCTTAGTTATCGTAGTGGTCACTTGGGGATTAATGTTGTTTAACTAAAAAACGTCGGCAAAAATGACCGCACTTTGCTCTGCACCCCAAAACCTGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTATCAATTTCGTTAAAAAACTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCGaaac                                                                                                              <  1:45632/214‑2 (MQ=255)                                                          gCCGTGTTTGTCGCCTTAGTTATCGTAGTGGTCACTTGGGGATTAATGTTGTTTAACTAAAAAACGTCGGCAAAAATGACCGCACTTTGCTCTGCACCCCAAAACCTGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTATCAATTTCGTTAAAAAACTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCGaaac                                                                                                              <  1:668595/214‑2 (MQ=255)                                                                         ttAGTTATCGTAGTGGTCACTTGGGGATTAATGTTGTTTAACTAAAAAACGTCGGCAAAAATGACCGCACTTTGCTCTGCACCCCAAAACCTGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTATCAATTTCGTTAAAAAACTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCTAAATTTCATCAGATTTCTCGTTCTCAAGTCATTTATTGGAAAAGCCTTTCGTGaaa                                                          >  1:524790/1‑251 (MQ=255)                                                                          tAGTTATCGTAGTGGTCACTTGGGGATTAATGTTGTTTAACTAAAAAACGTCGGCAAAAATGACCGCACTTTGCTCTGCACCCCAAAACCTGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTATCAATTTCGTTAAAAAACTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCTAAATTTCATCAGATTTCTCGTTCTCAAGTCATTTATTGGAAAAGCCTTTCGTGAAag                                                         <  1:509496/251‑1 (MQ=255)                                                                           aGTTATCGTAGTGGTCACTTGGGGATTAATGTTGTTTAACTAAAAAACGTCGGCAAAAATGACCGCACTTTGCTCTGCACCCCAAAACCTGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTATCAATTTCGTTAAAAAACTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCTAAATTTCATCAGATTTCTCGTTCTCAAGTCATTTATTGGAAAAGCCTTTCGTGAAaga                                                        >  1:766715/1‑251 (MQ=255)                                                                            gTTATCGTAGTGGTCACTTGGGGATTAATGTTGTTTAACTAAAAAACGTCGGCAAAAATGACCGCACTTTGCTCTGCACCCCAAAACCTGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTATCAATCTCGTTAAAAAACTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCTAAATTTCATCAGATTTCTCGTTCTCAAGTCATTTATTGGAAAAGCCTTTCGTGAAag                                                         <  1:877314/249‑1 (MQ=255)                                                                             ttATCGTAGTGGTCACTTGGGGATTAATGTTGTTTAACTAAAAAACGTCGGCAAAAATGACCGCACTTTGCTCTGCACCCCAAAACCTGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTATCAATTTCGTTAAAAAACTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCTAAATTTCATCAGATTTCTCGTTCTCAAGTCATTTATTGGAAAAGCCTTTCGTGAAAGAgg                                                      >  1:305195/1‑251 (MQ=255)                                                                               aTCGTAGTGGTCACTTGGGGATTAATGTTGTTTAACTAAAAAACGTCGGCAAAAATGACCGCACTTTGCTCTGCACCCCAAAACCTGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTATCAATTTCGTTAAAAAACTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCTAAATTTCATCAGATTTCTCGTTCTCAAGTCATTTATTGGAAAAGCCTTTCGTGAAAGAGGGc                                                    <  1:485290/251‑1 (MQ=255)                                                                                tCGTAGTGGTCACTTGGGGATTAATGTTGTTTAACTAAAAAACGTCGGCAAAAATGACCGCACTTTGCTCTGCACCCCAAAACCTGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTATCAATTTCGTTAAAAAACTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCTAAAtt                                                                                                             <  1:455492/193‑1 (MQ=255)                                                                                      tGGTCACTTGGGGATTAATGTTGTTTAACTAAAAAACGTCGGCAAAAATGACCGCACTTTGCTCTGCACCCCAAAACCTGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTATCAATTTCGTTAAAAAACTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCTAAATTTCATCAGATTTCTCGTTCTCAAGTCATTTATTGGAAAAGCCTTTCGTGAAAGAGGGCTTAATgg                                             <  1:193772/251‑1 (MQ=255)                                                                                        gTCACTTGGGGATTAATGTTGTTTAACTAAAAAACGTCGGCAAAAATGACCGCACTTTGCTCTGCACCCCAAAACCTGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTATCAATTTCGTTAAAAAACTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCTAAATTTCATCAGATTTCTCGTTCTCAAGTCATTTATTGGAAAAGCCTTTCGTGAAAGAGGGCTTAATGGCg                                           >  1:480213/1‑251 (MQ=255)                                                                                                 ggATTAATGTTGTTTAACTAAAAAACGTCGGCAAAAATGACCGCACTTTGCTCTGCACCCCAAAACCTGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTATCAATTTCGTTAAAAAACTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCTAAATTTCATCAGATTTCTCGTTCTCAAGTCATTTATTGGAAAAGCCTTTCGTGAAAGAGGGCTTAATGGCGTAAAATcc                                   >  1:785495/1‑250 (MQ=255)                                                                                                  gATTAATGTTGTTTAACTAAAAAACGTCGGCAAAAATGACCGCACTTTGCTCTGCACCCCAAAACCTGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTATCAATTTCGTTAAAAAACTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAgg                                                                                                                       <  1:716176/165‑1 (MQ=255)                                                                                                  gATTAATGTTGTTTAACTAAAAAACGTCGGCAAAAATGACCGCACTTTGCTCTGCACCCCAAAACCTGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTATCAATTTCGTTAAAAAACTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCTAAATTTCATCAGATTTCTCGTTCTCAAGTCATTTATTgg                                                                           >  1:737416/1‑209 (MQ=255)                                                                                                  gATTAATGTTGTTTAACTAAAAAACGTCGGCAAAAATGACCGCACTTTGCTCTGCACCCCAAAACCTGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTATCAATTTCGTTAAAAAACTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCTAAATTTCATCAGATTTCTCGTTCTCAAGTCATTTATTGGAAAAGCCTTTCGTGAAAGAGGGCTTAATGGCGTAAAATccc                                  >  1:50358/1‑250 (MQ=255)                                                                                                  gATTAATGTTGTTTAACTAAAAAACGTCGGCAAAAATGACCGCACTTTGCTCTGCACCCCAAAACCTGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTATCAATTTCGTTAAAAAACTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCTAAATTTCATCAGATTTCTCGTTCTCAAGTCATTTATTGGAAAAGCCTTTCGTGAAAGAGGGCTTAATGGCGTAAAATccc                                  >  1:896022/1‑250 (MQ=255)                                                                                                   aTTAATGTTGTTTAACTAAAAAACGTCGGCAAAAATGACCGCACTTTGCTCTGCACCCCAAAACCTGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTATCAATTTCGTTAAAAAACTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCTAAATTTCATCAGATTTCTCGTTCTCAAGTCATTTATTGGAAAAGCCTTTCGTGAAAGAGGGCTTAATGGCGTAAAATcccc                                 <  1:472858/250‑1 (MQ=255)                                                                                                      aaTGTTGTTTAACTAAAAAACGTCGGCAAAAATGACCGCACTTTGCTCTGCACCCCAAAACCTGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTATCAATTTCGTTAAAAAACTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCTAAATTTCATCAGATTTCTCGTTCTCAAGTCATTTATTGGAAAAGCCTTTCGTGAAAGAGGGCTTAATGGCGTAAAATCCCCTtat                             >  1:534749/1‑251 (MQ=255)                                                                                                        tgttgtTTAACTAAAAAACGTCGGCAAAAATGACCGCACTTTGCTCTGCACCCCAAAACCTGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTATCAATTTCGTTAAAAAACTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCTAAATTTCATCAGATTTCTCGTTCTCAAGTCATTTATTGGAAAAGCCTTTCGTGAAAGAGGGCTTAATGGCGTAAAATCCCCTtatat                           <  1:487212/251‑1 (MQ=255)                                                                                                        tgttgtTTAACTAAAAAACGTCGGCAAAAATGACCGCACTTTGCTCTGCACCCCAAAACCTGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTATCAATTTCGTTAAAAAACTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCTAAATTTCATCAGATTTCTCGTTCTCAAGTCATTTATTGGAAAAGCCTTTCGTGAAAGAGGGCTTAATGGCGTAAAATCCCCTtatat                           <  1:379174/251‑1 (MQ=255)                                                                                                        tgttgtTTAACTAAAAAACGTCGGCAAAAATGACCGCACTTTGCTCTGCACCCCAAAACCTGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTATCAATTTCGTTAAAAAACTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCTAAATTTCATCAGATTTCTCGTTCTCAAGTCATTTATTGGAAAAGCCTTTCGTGAAAGAGGGCTTAATGGCGTAAAATCCCCTtatat                           <  1:301352/251‑1 (MQ=255)                                                                                                         gttgttTAACTAAAAAACGTCGGCAAAAATGACCGCACTTTGCTCTGCACCCCAAAACCTGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTATCAATTTCGTTAAAAAACTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGt                                                                                                                      <  1:739342/159‑1 (MQ=255)                                                                                                          ttgttTAACTAAAAAACGTCGGCAAAAATGACCGCACTTTGCTCTGCACCCCAAAACCTGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTATCAATTTCGTTAAAAAACTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCTAAATTTCATCAGATTTCTCGTTCTCAAGTCATTTATTGGAAAAGCCTTTCGTGAAAGAGGGCTTAATGGCGTAAAATCCCCTTATATaa                         >  1:168261/1‑251 (MQ=255)                                                                                                           tgttTAACTAAAAAACGTCGGCAAAAATGACCGCACTTTGCTCTGCACCCCAAAACCTGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTATCAATTTCGTTAAAAAACTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCTAAATTTCATCAGATTTCTCg                                                                                              >  1:622019/1‑181 (MQ=255)                                                                                                             tttAACTAAAAAACGTCGGCAAAAATGACCGCACTTTGCTCTGCACCCCAAAACCTGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTATCAATTTCGTTAAAAAACTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCTAAATTTCATCAGATTTCTCGTTCTCAAGTCATTTATTGGAAAAGCCTTTCGTGAAAGAGGGCTTAATGGCGTAAAATCCCCTTATATAAAcc                      <  1:357001/251‑1 (MQ=255)                                                                                                              ttAACTAAAAAACGTCGGCAAAAATGACCGCACTTTGCTCTGCACCCCAAAACCTGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTATCAATTTCGTTAAAAAACTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCTAAATTTCATCAGATTTCTCGTTCTCAAGTCATTTATTGGAAAAGCCTTTCGTGAAAGAGGGCTTAATGGCGTAAAATCCCCTTATATAAAccc                     <  1:880833/251‑1 (MQ=255)                                                                                                                 aCTAAAAAACGTCGGCAAAAATGACCGCACTTTGCTCTGCACCCCAAAACCTGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTATCAATTTCGTTAAAAAACTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCTAAATTTCATCAGATTTCTCGTTCTCAAGTCATTTATTGGAAAAGCCTTTCGTGAAAGAGGGCTTAATGGCGTAAAATCCCCTtatat                           <  1:187441/242‑1 (MQ=255)                                                                                                                   tAAAAAACGTCGGCAAAAATGACCGCACTTTGCTCTGCACCCCAAAACCTGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTATCAATTTCGTTAAAAAACTACCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCTAAATTTCATCAGAATTCTCGTTCTCAAGTCATTTATTGGAAAAGCCTTTCGTGAAAGAGGGCTTAATGGCGTAAAATCCCCTTATATAAACCCTCAAc                <  1:608973/251‑1 (MQ=255)                                                                                                                     aaaaaCGTCGGCAAAAATGACCGCACTTTGCTCTGCACCCCAAAACCTGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTATCAATTTCGTTAAAAAACTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCGaaac                                                                                                              <  1:275424/155‑2 (MQ=255)                                                                                                                             cGGCAAAAATGACCGCACTTTGCTCTGCACCCCAAAACCTGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTATCAATTTCGTTAAAAAACTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCTAAATTTCATCAGATTTCTCGTTCTCAAGTCATTTATTGGAAAAGCCTTTCGTGAAAGAGGGCTTAATGGCGTAAAATCCCCTTATATAAACCCTCAACGCCCTAAAAt      <  1:807172/251‑1 (MQ=21)                                                                                                                                 aaaaaTGACCGCACTTTGCTCTGCACCCCAAAACCTGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTATCAATTTCGTTAAAAAACTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCTAAATTTCATCAGATTTCTCGTTCTCAAGTCATTTATTGGAAAAGCCTTTCGTGAAAGAGGGCTTAATGGCGTAAAATCCCCTTATATAAACCCTCAACGCCCTAAAATAGTg  <  1:242576/251‑1 (MQ=17)                                                                                                                                                                                                                                                    |                                                                                                                                        CGCATCGGCACCGAACGCCACGAACTCTCCGGGCGTGCCAAAGATCAAGGCTCCGCTGCCGTGTTTGTCGCCTTAGTTATCGTAGTGGTCACTTGGGGATTAATGTTGTTTAACTAAAAAACGTCGGCAAAAATGACCGCACTTTGCTCTGCACCCCAAAACCTGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTATCAATTTCGTTAAAAAACTATCAAACAGGTCACCGAGCAGAATTTTGGTATCCGTGAGGTGGCTAAATTTCATCAGATTTCTCGTTCTCAAGTCATTTATTGGAAAAGCCTTTCGTGAAAGAGGGCTTAATGGCGTAAAATCCCCTTATATAAACCCTCAACGCCCTAAAATAGTG  >  1/1622948‑1623326

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