Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	REL606	3,959,951	C→T	A44A (GCC→GCT)	uvrD →	DNA‑dependent ATPase I and helicase II

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	REL606	3,959,951	0	C	T	100.0%	224.1 / NA	57	A44A (GCC→GCT)	uvrD	DNA‑dependent ATPase I and helicase II
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (18/39);  total (18/39)

CGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGTGGTAAGACGCGCGTACTGGTGCATCGTATCGCCTGGCTGATGAGCGTGGAAAACTGCTCTCCATACTCGATTATGGCGGTGACGTTTACCAACAAAGCGGCG  >  REL606/3959883‑3960020                                                                     |                                                                      cGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGTGGTAAGACGCGCGTACTGGTGCATCGTATCGCTTGGCtga                                                                <  2:1662213/76‑1 (MQ=255) cGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGTGGTAAGACGCGCGTACTGGTGCATCGTATCGCTTGGCtga                                                                <  2:136413/76‑1 (MQ=255) cGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGTGGTAAGACGCGCGTACTGGTGCATCGTATCGCTTGGCtga                                                                >  2:1147306/1‑76 (MQ=255) cGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGTGGTAAGACGCGCGTACTGGTGCATCGTATCGCTTGGCtga                                                                >  1:1865583/1‑76 (MQ=255)   cagcaACCTTCTGGTGCTGGCGGGCGCGGGCAGTGGTAAGACGCGCGTACTGGTGCATCGTATCGCTTGGCtgatg                                                              <  1:1933046/76‑1 (MQ=255)   caacaACCTTCTGGTGCTGGCGGGCGCGGGCAGTGGTAAGACGCGCGTACTGGTGCATCGTATCGCTTGGCtgatg                                                              >  2:220285/4‑76 (MQ=255)     gcaACCTTCTGGTGCTGGCGGGCGCGGGCTGTGGTAAGACTCGCGTACTGGTGGATCGTATCGCTTGGCtgatggg                                                            >  2:951705/1‑74 (MQ=38)       aaCCTTCTGGTGCTGGCGGGCGCGGGCAGTGGTAAGACGCGCGTACTGGTGCATCGTATCGCTTGGCTGATGGGCg                                                          >  2:914574/1‑76 (MQ=255)       aaCCTTCTGGTGCTGGCGGGCGCGGGCAGTGGTAAGACGCGCGTACTGGTGCATCGTATCGCTTGGCTGATGAGCg                                                          >  2:482565/1‑76 (MQ=255)       aaCCTTCTGGGGCTGGCGGGCGCGGGCAGAGGTAAGACGCGCGTCCTGGTGCATCGTATCGTTTGGGTGCTGAGCg                                                          >  2:218070/1‑76 (MQ=255)        aCCTTCTGGTGCTGGCGGGCGCGGGCAGTGGTAAGACGCGCGTACTGGTGCATCGTATCGCTTGGCTGATGAGCGt                                                         >  2:1618268/1‑76 (MQ=255)        aCCTTCTGGTGCTGGCGGGCGCGGGCAGTGGTAAGACGCGCGTACTGGTGCATCGTATCGCTTGGCTGATGAGCGt                                                         >  1:1643756/1‑76 (MQ=255)         cgttCTGGTGCGGGGGGGCGGGGGCAGTGGTAAGACGTGCGTACTGGTGCATCGTATCGCTTGGCTGATGAGCGTg                                                        <  2:520854/74‑1 (MQ=255)             ggggTGCGGGGGGGCGGGGGCAGTGGTAAGACGCGCGTACTGGTGCATCGTATCGCTTGGCTGATGAGCGTGGaaa                                                    <  2:990741/74‑1 (MQ=255)             cTGGTGCTGGCGGGCGCGGGCAGTGGTAAGACGCGCGTACTGGTGCATCGTATCGCTTGGCTGATGAGCGTGGaaa                                                    <  2:65075/76‑1 (MQ=255)                     ggggggcgggggcAGTGGTAAGACGCGCGTACTGGTGCATCGTATCGCTTGGCTGATGAGCGTGGAAAACTGctct                                            <  2:799131/73‑1 (MQ=255)                         ggcgcgggcAGTGGTAAGACGCGCGTACTGGTGCATCGTATCGCTTGGCTGATGAGCGTGGAAAACTGCTCTCCAt                                        <  1:1829872/76‑1 (MQ=255)                          gcgggggcAGTGGTAAGACGCGCGTACTGGTGCATCGTATCGCTTGGCTGATGAGCGTGGAAAACTGCTCTCCATa                                       <  1:1268200/76‑1 (MQ=255)                          gcgcgggcAGTGGTAAGACGCGCGTACTGGTGCATCGTATCGCTTGGCTGATGAGCGTGGAAAACTGCTCTCCATa                                       <  1:1306344/76‑1 (MQ=255)                          gcgcgggcAGTGGTAAGACGCGCGTACTGGTGCATCGTATCGCTTGGCTGATGAGCGTGGAAAACTGCTCTCCATa                                       <  2:306773/76‑1 (MQ=255)                           ggggggcAGTGGTAAGACGCGCGTACTGGTGCATCGTATCGCTTGGCTGATGAGCGTGGAAAACTGCTCTCCATAc                                      <  1:910921/73‑1 (MQ=255)                            gggggcAGTGGTAAGACGCGCGTACTGGTGCATCGTATCGCTTGGCTGATGAGCGTGGAAAACTGCTCTCCATACt                                     <  2:320116/74‑1 (MQ=255)                            gcgggcAGTGGTAAGACGCGCGTACTGGTGCATCGTATCGCTTGGCTGATGAGCGTGGAAAACTGCTCTCCATACt                                     <  1:2064645/76‑1 (MQ=255)                             cgggcAGTGGTAAGACGCGCGTACTGGTGCATCGTATCGCTTGGCTGATGAGCGTGGAAAACTGCTCTCCATACTc                                    >  2:286983/1‑76 (MQ=255)                              gggcAGTGGTAAGACGCGCGTACTGGTGCATCGTATCGCTTGGCTGATGAGCGTGGAAAACTGCTCTCCATACTCg                                   <  1:100229/76‑1 (MQ=255)                                gcAGTGGTAAGACGCGCGTACTGGTGCATCGTATCGCTTGGCTGATGAGCGTGGAAAACTGCTCTCCATACTCGAt                                 >  2:1617129/1‑76 (MQ=255)                                 cAGTGGTAAGACGCGCGTACTGGTGCATCGTATCGCTTGGCTGATGAGCGTGGAAAACTGCTCTCCATACTCGAtt                                <  1:1367035/76‑1 (MQ=255)                                  aGTGGTAAGACGCGCGTACTGGTGCATCGTATCGCTTGGCTGATGAGCGTGGAAAACTGCTCTCCATACTCGATTa                               <  2:365329/76‑1 (MQ=255)                                  aGTGGTAAGACGCGCGTACTGGTGCATCGTATCGCTTGGCTGATGAGCGTGGAAAACTGCTCTCCATACTCGATTa                               <  2:858338/76‑1 (MQ=255)                                   gTGGTAAGACGCGCGTACTGGTGCATCGTATCGCTTGGCTGATGAGCGTGGAAAACTGCTCTCCATACTCGATTAt                              <  1:1251065/76‑1 (MQ=255)                                   gTGGTAAGACGCGCGTACTGGTGCATCGTATCGCTTGGCTGATGAGCGTGGAAAACTGCTCTCCATACTCGATTAt                              <  2:1269850/76‑1 (MQ=255)                                    tGGTAAGACGCGCGTACTGGTGCATCGTATCGCTTGGCTGATGAGCGTGGAAAACTGCTCTCCATACTCGATTATg                             >  2:1311588/1‑76 (MQ=255)                                         aGACGCGCGTACTGGTGCATCGTATCGCTTGGCTGATGAGCGTGGAAAACTGCTCTCCATACTCGATTATGGCggg                        >  1:1545832/1‑75 (MQ=255)                                          gACGCGCGTACTGGTGCATCGTATCGCTTGGCTGATGAGCGTGGAAAACTGCTCTCCATACTCGATTATGGCGGTg                       <  2:851994/76‑1 (MQ=255)                                          gACGCGCGTACTGGTGCATCGTATCGCTTGGCTGATGAGCGTGGAAAACTGCTCTCCATACTCGATTATGGCGGTg                       <  2:1848942/76‑1 (MQ=255)                                           aCGCGCGTACTGGTGCATCGTATCGCTTGGCTGATGAGCGTGGAAAACTGCTCTCCATACTCGATTATGGCGGTGa                      <  1:568530/76‑1 (MQ=255)                                            cgcgcgTACTGGTGCATCGTATCGCTTGGCTGATGAGCGTGGAAAACTGCTCTCCATACTCGATTATGGCGGTGAc                     <  1:502279/76‑1 (MQ=255)                                            cgcgcgTACTGGTGCATCGTATCGCTTGGCTGATGAGCGTGGAAAACTGCTCTCCATACTCGATTATGGCGGTGAc                     <  2:1952827/76‑1 (MQ=255)                                            cgcgcgTACTGGTGCATCGTATCGCTTGGCTGATGAGCGTGGAAAACTGCTCTCCATACTCGATTATGGCGGTGAc                     <  1:703766/76‑1 (MQ=255)                                              cgcgTACTGGTGCATCGTATCGCTTGGCTGATGAGCGTGGAAAACTGCTCTCCATACTCGATTATGGCGGGGACGt                   >  2:2025557/1‑76 (MQ=255)                                               gcgTACTGGTGCATCGTATCGCTTGGCTGATGAGCGTGGAAAACTGCTCTCCATACTCGATTATGGCGGTGACGtt                  <  1:1455004/76‑1 (MQ=255)                                               gcgTACTGGTGCATCGTATCGCTTGGCTGATGAGCGTGGAAAACTGCTCTCCATACTCGATTATGGCGGTGACGtt                  <  2:503517/76‑1 (MQ=255)                                                cgTACTGGTGCATCGTATCGCTTGGCTGATGAGCGTGGAAAACTGCTCTCCATACTCGATTATGGCGGTGACGttt                 <  2:1142392/76‑1 (MQ=255)                                                cgTACTGGTGCATCGTATCGCTTGGCTGATGAGCGTGGAAAACTGCTCTCCATACTCGATTATGGCGGGGACGctt                 >  1:905925/1‑76 (MQ=255)                                                 gTACTGGTGCATCGTATCGCTTGGCTGATGAGCGTGGAAAACTGCTCTCCATACTCGATTATGGCGGTGAGGTTTa                >  1:1132691/1‑76 (MQ=255)                                                 gTACTGGTGCATCGTATCGCTTGGCTGATGAGCGTGGAAAACTGCTCTCCATACTCGATTATGGCGGGGGCGtttt                >  1:1197379/1‑75 (MQ=255)                                                  ttcTGGTGCATCGTATCGCTTGGCTGATGAGCGTGGAAAACTGCTCTCCATACTCGATTATGGCGGTGACGTTTAc               <  1:1582537/74‑1 (MQ=255)                                                    cTGGTGCATCGTATCGCTTGGCTGATGAGCGTGGAAAACTGCTCTCCATACTCGATTATGGCGGTGAGGTTTACca             >  1:29946/1‑76 (MQ=255)                                                    cTGGTGCATCGTATCGCTTGGCTGATGAGCGTGGAAAACTGCTCTCCATACTCGATTATGGCGGTGACGTTTACca             <  1:1914809/76‑1 (MQ=255)                                                     tGGTGCATCGTATCGCTTGGCTGATGAGCGTGGAAAACTGCTCTCCATACTCGATTATGGCGGTGACGTTTACcaa            <  2:1046871/76‑1 (MQ=255)                                                     tGGTGCATCGTATCGCTTGGCTGATGAGCGTGGAAAACTGCTCTCCATACGCGATTATGGCGGTGACGTTTACcaa            <  1:1324055/76‑1 (MQ=255)                                                      ggTGCATCGTATCGCTTGGCTGATGAGCGTGGAAAACTGCTCTCCATACTCGATTATGGCGGTGACGTTTACcaaa           >  1:1998661/1‑75 (MQ=255)                                                       gTGCATCGTATCGCTTGGCTGATGAGCGTGGAAAACTGCTCTCCATACTCGATTATGGCGGTGACGTTTACcaaca          <  2:24387/76‑1 (MQ=255)                                                         gCATCGTATCGCTTGGCTGATGAGCGTGGAAAACTGCTCTCCATACTCGATTATGGCGGTGACGTTTACCAACaaa        <  1:1259766/76‑1 (MQ=255)                                                         gCATCGTATCGCTTGGCTGATGAGCGTGGAAAACTGCTCTCCATACTCGATTATGGCGGTGACGTTTACCAACaaa        <  1:298885/76‑1 (MQ=255)                                                         gCATCGTATCGCTTGGCTGATGAGCGTGGAAAACTGCTCTCCATACTCGATTATGGCGGTGACGTTTACCAACaaa        <  1:130369/76‑1 (MQ=255)                                                          cATCGTATCGCTTGGCTGATGAGCGTGGAAAACTGCTCTCCATACTCGATTATGGCGGTGACGTTTTCCAACaaaa       >  2:779833/1‑75 (MQ=255)                                                           aTCGTATCGCTTGGCTGATGAGCGTGGAAAACTGCTCTCCATACTCGATTATGGCGGTGACGTTTACCAACAAAgc      <  2:1899215/76‑1 (MQ=255)                                                              gTATCGCTTGGCTGATGAGCGTGGAAAACTGCTCTCCATACTCGATTATGGCGGTGACGTTTACCAACAAAgcggc   <  1:866268/76‑1 (MQ=255)                                                               tATCGCTTGGCTGATGAGCGTGGAAAACTGCTCTCCATACTCGATTATGGCGGTGACGTTTACCAACAAAgcggcg  >  2:756062/1‑76 (MQ=255)                                                                     |                                                                      CGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGTGGTAAGACGCGCGTACTGGTGCATCGTATCGCCTGGCTGATGAGCGTGGAAAACTGCTCTCCATACTCGATTATGGCGGTGACGTTTACCAACAAAGCGGCG  >  REL606/3959883‑3960020

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