Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	REL606	3,067,926	G→A	A40V (GCA→GTA)	yghS ←	predicted protein with nucleoside triphosphate hydrolase domain

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	REL606	3,067,926	0	G	A	100.0%	233.9 / NA	68	A40V (GCA→GTA)	yghS	predicted protein with nucleoside triphosphate hydrolase domain
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (50/18);  total (50/18)

ATTCGCCCGGACGATTGCCCGAGATAAATGTGTTCTGTTGGCATTCTTGCTGCCAGTTCATTTACCAGGCTTGCCGTGAGGGTCGATTTACCTGAA  >  REL606/3067878‑3067973                                                 |                                                aTTCGCCCGGACGATTGCCCGAGATAAATGTGTTCTGTTGGCATTCTtact                                               <  1:6112183/51‑1 (MQ=255)  ttCGCCCGGACGATTGCCCGAGATAAATGTGTTCTGTTGGCATTCTtactg                                              <  1:6090288/51‑1 (MQ=255)  ttCGCCCGGACGATTGCCCGAGATAAATGTGTTCTGTTGGCATTCTtactg                                              <  1:3409742/51‑1 (MQ=255)  ttCGCCCGGACGATTGCCCGAGATAAATGTGTTCTGTTGGCATTCTtactg                                              >  1:4097589/1‑51 (MQ=255)   tCGCCCGGACGATTGCCCGAGATAAATGTGTTCTGTTGGCATTCTtactgc                                             >  1:762062/1‑51 (MQ=255)     gCCCGGACGATTGCCCGAGATAAATGTGTTCTGTTGGCATTCTTACTGCCa                                           >  1:1099011/1‑51 (MQ=255)     gCCCGGACGATTGCCCGAGATAAATGTGTTCTGTTGGCATTCTTACTGCCa                                           >  1:3181213/1‑51 (MQ=255)      cccGGACGATTGCCCGAGATAAATGTGTTCTGTTGGCATTCTTACTGCCAg                                          >  1:5331681/1‑51 (MQ=255)      cccGGACGATTGCCCGAGATAAATGTGTTCTGTTGGCATTCTTACTGCCAg                                          >  1:1256494/1‑51 (MQ=255)      cccGGACGATTGCCCGAGATAAATGTGTTCTGTTGGCATTCTTACTGCCAg                                          <  1:2713634/51‑1 (MQ=255)       ccGGACGATTGCCCGAGATAAATGTGTTCTGTTGGCATTCTTACTGCCAGt                                         >  1:1260505/1‑51 (MQ=255)        cNGACGATTGCCCGAGATAAATGTGTTCTGTTGGCATTCTTACTGCCAGtt                                        >  1:5197768/1‑51 (MQ=255)        cGGACGATTGCCCGAGATAAATGTGTTCTGTTGGCATTCTTACTGCCAGtt                                        >  1:5170635/1‑51 (MQ=255)         ggACGATTGCCCGAGATAAATGTGTTCTGTTGGCATTCTTACTGCCAGTTc                                       <  1:6509822/51‑1 (MQ=255)          gACGATTGCCCGAGATAAATGTGTTCTGTTGGCATTCTTACTGCCAGTTCa                                      <  1:7291065/51‑1 (MQ=255)            cGATTGCCCGAGATAAATGTGTTCTGTTGGCATTCTTACTGCCAGTTCAtt                                    <  1:6112565/51‑1 (MQ=255)              aTTGCCCGAGATAAATGTGTTCTGTTGGCATTCTTACTGCCAGTTCATTTa                                  >  1:6480181/1‑51 (MQ=255)               ttGCCCGAGATAAATGTGTTCTGTTGGCATTCTTACTGCCAGTTCATTTAc                                 >  1:2299604/1‑51 (MQ=255)                tGCCCGAGATAAATGTGTTCTGTTGGCATTCTTACTGCCAGTTCATTTAcc                                >  1:32116/1‑51 (MQ=255)                 gCCCGAGATAAATGTGTTCTGTTGGCATTCTTACTGCCAGTTCATTTACCa                               <  1:5144940/51‑1 (MQ=255)                 gCCCGAGATAAATGTGTTCTGTTGGCATTCTTACTGCCAGTTCATTTACCa                               <  1:2694689/51‑1 (MQ=255)                 gCCCGAGATAAATGTGTTCTGTTGGCATTCTTACCGCCAGTTCATTTACCa                               >  1:6536822/1‑51 (MQ=255)                  cccGAGATAAATGTGTTCTGTTGGCATTCTTACTGCCAGTTCATTTACCAg                              >  1:6886929/1‑51 (MQ=255)                   ccGAGATAAATGTGTTCTGTTGGCATTCTTACTGCCAGTTCATTTACCAgg                             <  1:7007762/51‑1 (MQ=255)                   ccGAGATAAATGTGTTCTGTTGGCATTCTTACTGCCAGTTCATTTACCAgg                             >  1:3076228/1‑51 (MQ=255)                   ccGAGATAAATGTGTTCTGTTGGCATTCTTACTGCCAGTTCATTTACCAgg                             <  1:3779506/51‑1 (MQ=255)                   ccGAGATAAATGTGTTCTGTTGGCATTCTTACTGCCAGTTCATTTACCAgg                             <  1:3183008/51‑1 (MQ=255)                   ccGAGATAAATGTGTTCTGTTGGCATTCTTACTGCCAGTTCATTTACCAgg                             >  1:3997107/1‑51 (MQ=255)                   ccGAGATAAATGTGTTCTGTTGGCATTCTTACTGCCAGTTCATTTACCAgg                             >  1:503045/1‑51 (MQ=255)                    cGAGATAAATGTGTTCTGTTGGCATTCTTACTGCCAGTTCATTTACCa                               >  1:685817/1‑48 (MQ=255)                    cGAGATAAATGTGTTCTGTTGGCATTCTTACTGCCAGTTCATTTACCAGGc                            >  1:1581594/1‑51 (MQ=255)                     gagaTAAATGTGTTCTGTTGGCATTCTTACTGCCAGTTCATTTACCAGGCt                           >  1:19787/1‑51 (MQ=255)                     gagaTAAATGTGTTCTGTTGGCATTCTTACTGCCAGTTCATTTACCAGGCt                           >  1:2068505/1‑51 (MQ=255)                      agaTAAATGTGTTCTGTTGGCATTCTTACTGCCAGTTCATTTACCAGGCtt                          <  1:5465507/51‑1 (MQ=255)                      agaTAAATGTGTTCTGTTGGCATTCTTACTGCCAGTTCATTTACCAGGCtt                          >  1:2582798/1‑51 (MQ=255)                        aTAAATGTGTTCTGTTGGCATTCTTACTGCCAGTTCATTTACCAGGCTTGc                        >  1:3221740/1‑51 (MQ=255)                        aTAAATGTGTTCTGTTGGCATTCTTACTGCCAGTTCATTTACCAGGCTTGc                        >  1:1595859/1‑51 (MQ=255)                        aTAAATGTGTTCTGTTGGCATTCTTACTGCCAGTTCATTTACCAGGCTTGc                        >  1:3887086/1‑51 (MQ=255)                        aTAAATGTGTTCTGTTGGCATTCTTACTGCCAGTTCATTTACCAGGCTTGc                        >  1:4931556/1‑51 (MQ=255)                         tAAATGTGTTCTGTTGGCATTCTTACTGCCAGTTCATTTACCAGGCTTGcc                       >  1:4945320/1‑51 (MQ=255)                          aaaTGTGTTCTGTTGGCATTCTTACTGCCAGTTCATTTACCAGGc                            <  1:4740010/45‑1 (MQ=255)                          aaaTGTGTTCTGTTGGCATTCTTACTGCCAGTTCATTTACCAGGc                            <  1:5151043/45‑1 (MQ=255)                          aaaTGTGTTCTGTTGGCATTCTTACTGCCAGTTCATTTACCAGGCTTGCCg                      >  1:4740665/1‑51 (MQ=255)                          aaaTGTGTTCTGTTGGCATTCTTACTGCCAGTTCATTTACCAGGCTTGCCg                      >  1:3973470/1‑51 (MQ=255)                          aaaTGTGTTCTGTTGGCATTCTTACTGCCAGTTCAATTACCAGGCTTGCCg                      >  1:2366290/1‑51 (MQ=255)                           aaTGTGTTCTGTTGGCATTCTTACTGCCAGTTCATTTACCAGGCTTGCCGt                     >  1:2053874/1‑51 (MQ=255)                            aTGTGTTCTGTTGGCATTCTTACTGCCAGTTCATTTACCAGGCTTGCCGTg                    >  1:5495256/1‑51 (MQ=255)                            aTGTGTTCTGTTGGCATTCTTACTGCCAGTTCATTTACCAGGCTTGCCGTg                    >  1:5213986/1‑51 (MQ=255)                             tgtgTTCTGTTGGCATTCTTACTGCCAGTTCATTTACCAGGCTTGCCGTGa                   >  1:5172924/1‑51 (MQ=255)                              gtgtTCTGTTGGCATTCTTACTGCCAGTTCATTTACCAGGCTTGCCGTGAg                  >  1:6137849/1‑51 (MQ=255)                              gtgtTCTGTTGGCATTCTTACTGCCAGTTCATTTACCAGGCTTGCCGTGAg                  >  1:2437602/1‑51 (MQ=255)                               tgtTCTGTTGGCATTCTTACTGCCAGTTCATTTACCAGGCTTGCCGTGAgg                 >  1:6133247/1‑51 (MQ=255)                                gtTCTGTTGGCATTCTTACTGCCAGTTCATTTACCAGGCTTGCCGTGAggg                >  1:1112455/1‑51 (MQ=255)                                  tCTGTTGGCATTCTTACTGCCAGTTCATTTACCAGGCTTGCCGTGAGGGTc              <  1:4784104/51‑1 (MQ=255)                                   cTGTTGGCATTCTTACTGCCAGTTCATTTACCAGGCTTGCCGTGAGGGTCg             >  1:4503790/1‑51 (MQ=255)                                   cTGTTGGCATTCTTACTGCCAGTTCATTTACCAGGCTTGCCGTGAGGGTCg             >  1:5392741/1‑51 (MQ=255)                                    tGTTGGCATTCTTACTGCCAGTTCATTTACCAGGCTTGCCGTGAGGGTCGa            >  1:6962258/1‑51 (MQ=255)                                     gTTGGCATTCTTACTGCCAGTTCATTTACCAGGCTTGCCGTGAGGGTCGAt           >  1:6046415/1‑51 (MQ=255)                                     gTTGGCATTCTTACTGCCAGTTCATTTACCAGGCTTGCCGTGAGGGTCGAt           <  1:3829046/51‑1 (MQ=255)                                     gTTGGCATTCTTACTGCCAGTTCATTTACCAGGCTTGCCGTGAGGGTCGAt           >  1:4266834/1‑51 (MQ=255)                                      ttGGCATTCTTACTGCCAGTTCATTTACCAGGCTTGCCGTGAGGGTCGAtt          >  1:906839/1‑51 (MQ=255)                                      ttGGCATTCTTACTGCCAGTTCATTTACCAGGCTTGCCGTGAGGGTCGAtt          >  1:1819669/1‑51 (MQ=255)                                      ttGGCATTCTTACTGCCAGTTCATTTACCAGGCTTGCCGTGAGGGTCGAtt          >  1:3859866/1‑51 (MQ=255)                                      ttGGCATTCTTACTGCCAATTCATTTACCAGGCTTGCCGTGAGGGTCGAtt          >  1:3284694/1‑51 (MQ=255)                                       tGGCATTCTTACTGCCAGTTCATTTACCAGGCTTGCCGTGAGGGTCGAttt         >  1:5585135/1‑51 (MQ=255)                                       tGGCATTCTTACTGCCAGTTCATTTACCAGGCTTGCCGTGAGGGTCGAttt         >  1:489772/1‑51 (MQ=255)                                          cATTCTTACTGCCAGTTCATTTACCAGGCTTGCCGTGAGGGTCGATTTAcc      <  1:5770520/51‑1 (MQ=255)                                            ttCTTACTGCCAGTTCATTTACCAGGCTTGCCGTGAGGGTCGATTTACCTg    >  1:3392004/1‑51 (MQ=255)                                            ttCTTACTGCCAGTTCATTTACCAGGCTTGCCGTGAGGGTCGATTTACCTg    >  1:6439857/1‑51 (MQ=255)                                             tCTTACTGCCAGTTCATTTACCAGGCTTGCCGTGAGGGTCGATTTACCTGa   <  1:1400330/51‑1 (MQ=255)                                             tCTTACTGCCAGTTCATTTACCAGGCTTGCCGTGAGGGTCGATTTACCTGa   <  1:1160159/51‑1 (MQ=255)                                             tCTTACTGCCAGTTCATTTACCAGGCTTGCCGTGAGGGTCGATTTACCTGa   >  1:4084602/1‑51 (MQ=255)                                              cTTACTGCCAGTTCATTTACCAGGCTTGCCGTGAGGGTCGATTTACCTGaa  <  1:4054430/51‑1 (MQ=255)                                                 |                                                ATTCGCCCGGACGATTGCCCGAGATAAATGTGTTCTGTTGGCATTCTTGCTGCCAGTTCATTTACCAGGCTTGCCGTGAGGGTCGATTTACCTGAA  >  REL606/3067878‑3067973

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