Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	REL606	183,037	G→T	R320L (CGC→CTC)	dgt →	deoxyguanosinetriphosphate triphosphohydrolase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	REL606	183,037	0	G	T	100.0%	234.2 / NA	62	R320L (CGC→CTC)	dgt	deoxyguanosinetriphosphate triphosphohydrolase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base T (25/37);  total (25/37)

CACGAAGCGTGGGGCCAGCATGAGAAAGGTTCGCTCTTTTCGCTAGTGGTTGAAAATGCCTGGGAAAAATCACGCTCAAATAGTTTAAGCCGCAGTACGGAAGATCAGTTTTTTATGTATTTACGAGTAAACACCCTAAATAA  >  REL606/182964‑183106                                                                          |                                                                      cACGAAGCGTGGGGCCAGCATGAGAAAGGTTCGCTCTTTTCGCTAGTGGTTAAAAATGCCTGNGAAAAATCACGcg                                                                     >  1:201750/1‑75 (MQ=255)   cGAAGCGTGGGGCCAGCATGAGAAAGGTTCGCTCTTTTCGCTAGTGGTTGAAAATGCCTGGGAAAAATCACTCTCa                                                                   <  1:848143/76‑1 (MQ=255)    gAAGCGTGGGGCCAGCATGAGAAAGGTTCGCTCTTTTCGCTAGTGGTTGAAAATGCCTGGGAANAATCACTCTCaa                                                                  <  1:14123/76‑1 (MQ=255)    gAAGCGTGGGGCCAGCATGAGAAAGGTTCGCTCTTTTCGCTAGTGGTTGAAAATGCCTGGGAAAAATCACTCTCaa                                                                  >  1:614949/1‑76 (MQ=255)         gTGGGGCCAGCATGAGAAAGGTTCGCTCTTTTCGCTAGTGGTTGAAAATGCCTGGGAAAAATCACTCTCAAATAGt                                                             >  1:225285/1‑76 (MQ=255)          tGGGGCCAGCATGAGAAAGGTTCGCTCTTTTCGCTAGTGGTTGAAAATGCCTGGGAAAAATCACTCTCAAATAGtt                                                            >  2:69235/1‑76 (MQ=255)          tGGGGCCAGCATGAGAAAGGTTCGCTCTTTTCGCTAGTGGTTGAAAATGCCTGGGAAAAATCACTCTCAAATAGtt                                                            <  1:220290/76‑1 (MQ=255)          tGGGGCCAGCATGAGAAAGGTTCGCTCTTTTCGCTAGGGGTTGAAAATGCCTGGGAAAAATCACTCTCAAATAGtt                                                            <  1:220285/76‑1 (MQ=255)           ggggCCAGCATGAGAAAGGTTCGCTCTTTTCGCTAGTGGTTGAAAATGCCTGGGAAAAATCACTCTCAAATAGttt                                                           <  2:1002360/76‑1 (MQ=255)            gggCCAGCATGAGAAAGGTTCGCTCTTTTCGCTAGTGGTTGAAAATGCCTGGGAAAAATCACTCTCAAATAGTTTa                                                          <  2:883746/76‑1 (MQ=255)               ccAGCATGAGAAAGGTTCGCTCTTTTCGCTAGTAGTTGAAAATGCCTGGGAAAAATCACTCTCAAATAGTTTAAGc                                                       >  2:1312494/1‑76 (MQ=255)                 aGCATGAGAAAGGTTCGCTCTTTTCGCTAGTGGTTGAAAATGCCTGGGATAAATCACTCTCAAATAGTTTAAGCCg                                                     >  1:429417/1‑76 (MQ=255)                   cATGAGAANGGTTCGCTCTTTTCGCTAGTGGTTGAAAATGCCTGGGAAAAATCACTCTCAAATAGTTTAAGCCGCa                                                   <  1:1340379/76‑1 (MQ=255)                   cATGAGAAAGGTTCGCTCTTTTCGCTAGTGGTTGAAAATGCCTGGGAAAAATCACTCTCAAATAGTTTAAGCCGCa                                                   <  2:536441/76‑1 (MQ=255)                    aTGAGAAAGGTTCGCTCTTTTCGCTAGTGGTTGAAAATGCCTGGGAAAAATCACTCTCAAATAGTTTAAGCCGCAg                                                  >  2:580498/1‑76 (MQ=255)                     tGAGNAAGGTTCGCTCTTTTCGCTAGTGATTAAAAATGCNTGGGAAAAATCACTCTCAAATAGTTTAAGCCGCAGt                                                 <  1:187269/76‑1 (MQ=255)                     tGAGAAAGGTTCGCTCTTTTCGCTAGTGGTTGAAAATGCCTGGGAAAAATCACTCTCAAATAGTTTAAGCCGCAGt                                                 >  1:297829/1‑76 (MQ=255)                      gagaAAGGTTCGCTCTTTTCGCTAGTGGTTGAAAATGCCTGGGACAAATCACTCTCAAATAGTTTAAGCCGCAGTa                                                >  2:244762/1‑76 (MQ=255)                      gagaAAGGTTCGCTCTTTTCGCTAGTGGTTGAAAATGCCTGGGAAAAATCACTCTCAAATAGTTTAAGCCGCAGTa                                                <  2:1060959/76‑1 (MQ=255)                      gagaAAGGTTCGCTCTTTTCGCTAGTGGTTGAAAATGCCTGGGAAAAATCACTCTCAAATAGTTTAAGCCGCAGTa                                                <  1:1131892/76‑1 (MQ=255)                      gagaAAGGTTCGCTCTTTTCGCTAGTGGTTGAAAATGCCTGGGAAAAATCACTCTCAAATAGTTTAAGCCGCAGTa                                                <  1:1094467/76‑1 (MQ=255)                      gagaAAGGTTCGCTCTTTTCGCTAGTGGTTGAAAATGCCTGGGAAAAATCACTCTCAAATAGTTTAAGCCGCAGTa                                                >  1:581965/1‑76 (MQ=255)                       agaAAGGTTCGCTCTTTTCGCTAGTGGTTGAAAATGCCTGGGAAAAATCACTCTCAAATAGTTTAAGCCGCAGTAc                                               >  1:212857/1‑76 (MQ=255)                       agaAAGGTTCGCTCTTTTCGCTAGTGGTTGAAAATGCCTGGGAAAAATCACTCTCAAATAGTTTAAGCCGCAGTAc                                               >  2:1009841/1‑76 (MQ=255)                         aaaGGTTCGCTCTTTTCGCTAGTGGTTGAAAATGCCTGGGAAAAATCACTCTCAAATAGTTTAAGCCGCAGTACgg                                             >  2:1421059/1‑76 (MQ=255)                          aaGGTTCGCTCTTTTCGCTAGTGGTTGAAAATGCCTGGGAAAAATCACTCTCAAATAGTTTAAGCCGCAGTACGGa                                            >  2:354525/1‑76 (MQ=255)                            ggTTCGCTCTTTTCGCTAGTGGTTGAAAATGCCTGGGAAAAATCACTCTCAAATAGTTTAAGCCGCAGTACGGAAg                                          <  1:680739/76‑1 (MQ=255)                             gTTCGCTCTTTTCGCTAGTGGTTGAAAATGCCTGGGAAAAATCCCTCTCAAATAGTTTAAGCCGCAGTACGGAAGa                                         <  2:1155487/76‑1 (MQ=255)                             gTTCGCTCTTTTCGCTAGTGGTTGAAAATGCCTGGGAAAAATCACTCTCAAATAGTTTAAGCCGCAGTACGGAAGa                                         <  2:1155474/76‑1 (MQ=255)                              ttCGCTCTTTTCGCTAGTGGTTGAAAATGCCTGGGAAAAATCACTCTCAAATAGTTTAAGCCGCAGTACGGAAGAt                                        <  1:946714/76‑1 (MQ=255)                                cGCTCTTTTCGCTAGTGGTTGAAAATGCCTGGGAAAAATCACTCTCAANTAGTTTAAGCCGCAGTACGGAAGATCa                                      <  1:1524238/76‑1 (MQ=255)                                cGCTCTTTTCGCTAGTGGTTGAAAATGCCTGGGAAAAATCACTCTCAAATAGTTTAAGCCGCAGTACGGAAGATCa                                      <  2:384815/76‑1 (MQ=255)                                 gCTCTTTTCGCTAGTGGTTGAAAATGCCTGGCAAAAATCACTCTCAAATAGTTTAAGCCGCAGTACGGAAGATCAg                                     <  2:249369/76‑1 (MQ=255)                                   tctTTTCGCTAGTGGTTGAAAATGCCTGGGAAAAATCACTCTCAAATAGTTTAAGCCGCAGTACGGAAGATCAGtt                                   <  2:575318/76‑1 (MQ=255)                                   tctTTTCGCTAGTGGTTGAAAATGCCTGGGAAAAATCACTCTCAAATAGTTTAAGCCGCAGTACGGAAGATCAGtt                                   >  1:1332300/1‑76 (MQ=255)                                    ctTTTCGCTAGTGGTTGAAAATGCCTGGGAAAAATCACTCTCAAATAGTTTAAGCCGCAGTACGGAAGATCAGttt                                  >  1:872835/1‑76 (MQ=255)                                     ttttCGCTAGTGGTTGAAAATGCCTGGGAAAAATAACTCTCAAATAGTTTAAGCCGCAGTACGGAAGATCAGtttt                                 >  2:1533922/1‑76 (MQ=255)                                      tttCGCTAGTGGTTGAAAATCCTTGGGAAAAATCACTCTCAAATAGTTTAAGCCGCAGTACGGAAGATCAGttttt                                >  1:8349/1‑76 (MQ=255)                                      tttCGCTAGTGGTTGAAAATACCTGGGAAAAATCACTCTCAAATAGTTTAAGCCGCAGTACGGAAGATCAGttttt                                >  1:8358/1‑76 (MQ=255)                                         cGCTAGTGGTTGAAAATGCCTGGGAAAAATCACTCTCAAATAGTTTAAGCCGCAGTACGGAAGATCAGTTTTTTAt                             <  2:653557/76‑1 (MQ=255)                                           cTAGTGGTTGAAAATGCCTGGGAAAAATCACTCTCAAATAGTTTAAGCCGCAGTACGGAAGATCAGTTTTTTATGt                           <  2:497838/76‑1 (MQ=255)                                             aGTGGTTGAAAATGCCTGGGAAAAATCACTCTCAAATAGTTTAAGCCGCAGTACGGAAGATCAGTTTTTTATGTAt                         >  1:1098720/1‑76 (MQ=255)                                               tGGTTGAAAATGCCTGGGAAAAATCACTCTCAAATAGTTTAAGCCGCAGTACGGAAGATCAGTTTTTTATGTAttt                       <  1:1006581/76‑1 (MQ=255)                                                ggTTGAAAATGCCTGGGAAAAATCACTCTCAAATAGTTTAAGCCGCAGTACGGAAGATCAGTTTTTTATGTATTTa                      <  2:439019/76‑1 (MQ=255)                                                   tGAAAATGCCTGGGAAAAATCACTCTCAAATAGTTTAAGCCGCAGTACGGAAGATCAGTTTTTTATGTATTTACGa                   <  1:1193083/76‑1 (MQ=255)                                                     aaaaTGCCTGGGAAAAATCACTCTCAAATAGTTTAAGCCGCAGTACGGAAGATCAGTTTTTTATGTATTTACGAGt                 <  2:603246/76‑1 (MQ=255)                                                     aaaaTGCCTGGGAAAAATCACTCTCAAATAGTTTAAGCCGCAGTACGGAAGATCAGTTTTTTATGTATTTACGAGt                 <  1:1297360/76‑1 (MQ=255)                                                      aaaTGCCTGGGAAAAATCACTCTCAAATAGTTTAAGCCGCAGTACGGAAGATCAGTTTTTTATGTATTTACGAGTa                <  1:810242/76‑1 (MQ=255)                                                       aaTGCCTGTGAAAAATCACTCTCAAATAGCTAAAGCCGCAGTACGGAAGATCAGCTTTTTCTGTATCTCCGCGTat               >  2:365307/1‑75 (MQ=255)                                                       aaTGCCTGGGAAAAATCACTCTCAAATAGTTTAAGCCGCAGTACGGAAGATCAGTTTTTTATGTATTTACGAGTaa               >  2:365276/1‑76 (MQ=255)                                                         tGCCTGGGAAAAATCACTCTCAAATAGATTAAGCCGCAGTACGGAAGATCAGTTTTTTATGTATTTACGCGTAAac             >  1:1506769/1‑76 (MQ=255)                                                            cTGGGAAAAATCACTCTCAAATAGTTTAAGCCGCAGTACGGAAGATCAGTTTTTTATGTATTTACGAGTAAACAcc          <  2:1163687/76‑1 (MQ=255)                                                              gggAAAAATCACTCTCAAATAGTTTAAGCCGCAGTACGGAAGATCAGTTTTTTATGTATTTACGAGTAAACACCCt        <  1:1214519/76‑1 (MQ=255)                                                               ggAAAAATCACTCTCAAATAGTTTAAGCCGCAGTACGGAAGATCAGTTTTTTATGTATTTACGAGTAAACACCCta       <  2:479776/76‑1 (MQ=255)                                                               ggAAAAATCACTCTCAAATAGTTTAAGCCGCAGTACGGAAGATCAGTTTTTTATGTATTTACGAGTAAACACCCta       <  2:479791/76‑1 (MQ=255)                                                               ggAAAAATCACTCTCAAATAGTTTAAGCCGCAGTACGGAAGATCAGTTTTTTATGTATTTACGAGTAAACACCCta       <  2:16922/76‑1 (MQ=255)                                                               ggAAAAATCACTCTCAAATAGTTTAAGCCGCAGTACGGAAGATCAGTTTTTTATGTATTTACGAGTAAACACCCta       <  2:1523892/76‑1 (MQ=255)                                                               ggAAAAATCACTCTCAAATAGTTTAAGCCGAAGTACGGAAGATCAGTTTTTTATGTATTTACGAGTAAAGACCCta       >  1:1397278/1‑76 (MQ=255)                                                                 aaaaaTCACTCTCAAATAGTTTAAGTCGCAGTACGGAAGATCAGTTTTTTATGTATTTACGAGNAAACACCCtaaa     <  1:17090/76‑1 (MQ=255)                                                                   aaaTCACTCTCAAATAGTTTAAGCCGCAGTACGGAAGATCAGTTTTTTATGTATTTACGAGTAAACACCCtaaata   >  2:413106/1‑76 (MQ=255)                                                                   aaaTCACTCTCAAATAGTTTAAGCCGCAGTACGGAAGATCAGTTTTTTATGTATTTACGAGTAAACACCCtaaata   >  2:1185440/1‑76 (MQ=255)                                                                    aaTCACTCTCAAATAGTTTAAGCCGCAGTACGGAAGATCAGTTTTTTATGTATTTACGAGTAGACGCCCtaaataa  <  2:560526/76‑1 (MQ=255)                                                                    aaTCACTCTCAAATAGTTTAAGCCGCAGTACGGAAGATCAGTTTTTTATGTATTTACGAGTAAACACCCtaaataa  <  2:560525/76‑1 (MQ=255)                                                                          |                                                                      CACGAAGCGTGGGGCCAGCATGAGAAAGGTTCGCTCTTTTCGCTAGTGGTTGAAAATGCCTGGGAAAAATCACGCTCAAATAGTTTAAGCCGCAGTACGGAAGATCAGTTTTTTATGTATTTACGAGTAAACACCCTAAATAA  >  REL606/182964‑183106

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