Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	REL606	433,668	C→T	R55C (CGT→TGT)	hupB →	HU, DNA‑binding transcriptional regulator, beta subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	REL606	433,668	0	C	T	100.0%	175.3 / NA	52	R55C (CGT→TGT)	hupB	HU, DNA‑binding transcriptional regulator, beta subunit
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (29/23);  total (29/23)

GGGGATGATGTAGCACTGGTAGGTTTTGGTACTTTTGCCGTTAAAGAGCGTGCTGCCCGTACTGGCCGCAACCCGCAGACCGGTAAAGAGATCAC  >  REL606/433620‑433714                                                 |                                               ggggATGATGTAGCACTGGTAGGTTTTGGTACTTTTGCCGTTAAAGAgtgt                                              <  1:2240439/51‑1 (MQ=255)  gggATGATGTAGCACTGGTAGGTTTTGGTACTTTTGCCGTTAAAGAGtgtg                                             >  1:7085997/1‑51 (MQ=255)   ggATGATGTAGCACTGGTAGGTTTTGGTACTTTTGCCGTTAAAGAGTGtgc                                            <  1:1357150/51‑1 (MQ=255)     atgatgTAGCACTGGTAGGTTTTGGTACTTTTGCCGTTAAAGAGTGtgct                                           <  1:1174076/50‑1 (MQ=255)      tgatgTAGCACTGGTAGGTTTTGGTACTTTTGCCGTTAAAGAGTGtgctgc                                         >  1:6363718/1‑51 (MQ=255)      tgatgTAGCACTGGTAGGTTTTGGTACTTTTGCCGTTAAAGAGTGtgctgc                                         >  1:1947379/1‑51 (MQ=255)      tgatgTAGCACTGGTAGGTTTTGGTACTTTTGCCGTTAAAGAGTGtgctgc                                         >  1:6695412/1‑51 (MQ=255)       gatgTAGCACTGGTAGGTTTTGGTACTTTTGCCGTTAAAGAGTGTGCTGcc                                        >  1:3354298/1‑51 (MQ=255)       gatgTAGCACTGGTAGGTTTTGGTACTTTTGCCGTTAAAGAGTGTGCTGcc                                        >  1:1809454/1‑51 (MQ=255)        atgTAGCACTGGTAGGTTTTGGTACTTTTGCCGTTAAAGAGTGTGCTGccc                                       <  1:2652207/51‑1 (MQ=255)        atgTAGCACTGGTAGGTTTTGGTACTTTTGCCGTTAAAGAGTGTGCTGccc                                       >  1:3076715/1‑51 (MQ=255)        atgTAGCACTGGTAGGTTTTGGTACTTTTGCCGTTAAAGAGTGTGCTGccc                                       <  1:1333190/51‑1 (MQ=255)         tgTAGCACTGGTAGGTTTTGGTACTTTTGCCGTTAAAGAGTGTGCTGCCCg                                      >  1:1459699/1‑51 (MQ=255)          gTAGCACTGGTAGGTTTTGGTACTTTTGCCGTTAAAGAGTGTGCTGCCCGt                                     <  1:3048945/51‑1 (MQ=255)           tAGCACTGGTAGGTTTTGGTACTTTTGCCGTTAAAGAGTGTGCTGCCCGTa                                    <  1:6005412/51‑1 (MQ=255)            aGCACTGGTAGGTTTTGGTACTTTTGCCGTTAAAGAGTGTGCTGCCCGTAc                                   >  1:3718498/1‑51 (MQ=255)             gCACTGGTAGGTTTTGGTACTTTTGCCGTTAAAGAGTGTGCTGCCCGTACt                                  <  1:4245605/51‑1 (MQ=255)                 tGGTAGGTTTTGGTACTTTTGCCGTTAAAGAGTGTGCTGCCCGTACTGGcc                              >  1:2317134/1‑51 (MQ=255)                 tGGTAGGTTTTGGTACTTTTGCCGTTAAAGAGTGTGCTGCCCGTACTGGcc                              <  1:5029057/51‑1 (MQ=255)                 tGGTAGGTTTTGGTACTTTTGCCGTTAAAGAGTGTGCTGCCCGTACTGGcc                              <  1:2793128/51‑1 (MQ=255)                  ggTAGGTTTTGGTACTTTTGCCGTTAAAGAGTGTGCTGCCCGTACTGGCCg                             >  1:5574504/1‑51 (MQ=255)                   gTAGGTTTTGGTACTTTTGCCGTTAAAGAGTGTGCTGCCCGTACTGGCCGc                            <  1:5635032/51‑1 (MQ=255)                     aGGTTTTGGTACTTTTGCCGTTAAAGAGTGTGCTGCCCGTACTGGCCGCaa                          <  1:5631264/51‑1 (MQ=255)                       gTTTTGGTACTTTTGCCGTTAAAGAGTGTGCTGCCCGTACTGGCCGCAAcc                        >  1:7239597/1‑51 (MQ=255)                       gTTTTGGTACTTTTGCCGTTAAAGAGTGTGCTGCCCGTACTGGCCGCAAcc                        <  1:1211948/51‑1 (MQ=255)                       gTTTTGGTACTTTTGCCGTTAAAGAGTGTGCTGCCCGTACTGGCCGCAAcc                        <  1:7196513/51‑1 (MQ=255)                       gTTTTGGTACTTTTGCCGTTAAAGAGTGTGCTGCCCGTACTGGCCGCAAcc                        >  1:6474118/1‑51 (MQ=255)                        ttttGGTACTTTTGCCGTTAAAGAGTGTGCTGCCCGTACTGGCCGCAAccc                       >  1:2377045/1‑51 (MQ=255)                        ttttGGTACTTTTGCCGTTAAAGAGTGTGCTGCCCGTACTGGCCGCAAccc                       >  1:1505425/1‑51 (MQ=255)                         tttGGTACTTTTGCCGTTAAAGAGTGTGCTGCCCGTACTGGCCGCAACCCg                      >  1:3386107/1‑51 (MQ=255)                          ttGGTACTTTTGCCGTTAAAGAGTGTGCTGCCCGTACTGGCCGCAACCCGc                     >  1:2098535/1‑51 (MQ=255)                           tGGTACTTTTGCCGTTAAAGAGTGTGCTGCCCGTACTGGCCGCAACCCGCa                    <  1:6829674/51‑1 (MQ=255)                             gTACTTTTGCCGTTAAAGAGTGTGCTGCCCGTACTGGCCGCAACCCGCAGa                  >  1:6354949/1‑51 (MQ=255)                              tACTTTTGCCGTTAAAGAGTGTGCTGCCCGTACTGGCCGCAACCCGCAGAc                 <  1:1885912/51‑1 (MQ=255)                               aCTTTTGCCGTTAAAGAGTGTGCTGCCCGTACTGGCCGCAACCCGCAGAcc                <  1:140024/51‑1 (MQ=255)                                cTTTTGCCGTTAAAGAGTGTGCTGCCCGTACTGGCCGCAACCCGCAGAcc                <  1:7193907/50‑1 (MQ=255)                                 ttttGCCGTTAAAGAGTGTGCTGCCCGTACTGGCCGCAACCCGCAGACCgg              >  1:4723916/1‑51 (MQ=255)                                 ttttGCCGTTAAAGAGTGTGCTGCCCGTACTGGCCGCAACCCGCAGACCgg              >  1:2255600/1‑51 (MQ=255)                                  tttGCCGTTAAAGAGTGTGCTGCCCGTACTGGCCGCAACCCGCAGACCGGt             >  1:3710637/1‑51 (MQ=255)                                   ttGCCGTTAAAGAGTGTGCTGCCCGTACTGGCCGCAACCCGCAGACCGGTa            >  1:6463545/1‑51 (MQ=255)                                   ttGCCGTTAAAGAGTGTGCTGCCCGTACTGGCCGCAACCCGCAGACCGGTa            <  1:6920751/51‑1 (MQ=255)                                   ttGCCGTTAAAGAGTGTGCTGCCCGTACTGGCCGCAACCCGCAGACCGGTa            <  1:4059886/51‑1 (MQ=255)                                    tGCCGTTAAAGAGTGTGCTGCCCGTACTGGCCGCAACCCGCAGACCGGTaa           >  1:3767191/1‑51 (MQ=255)                                    tGCCGTTAAAGAGTGTGCTGCCCGTACTGGCCGCAACCCGCAGACCGGTaa           <  1:1138980/51‑1 (MQ=255)                                     gCCGTTAAAGAGTGTGCTGCCCGTACTGGCCGCAACCCGCAGACCGGTaaa          >  1:2957922/1‑51 (MQ=255)                                      ccGTTAAAGAGTGTGCTGCCCGTACTGGCCGCAACCCGCAGACCGGTAAag         >  1:6119421/1‑51 (MQ=255)                                      ccGTTAAAGAGTGTGCTGCCCGTACTGGCCGCAACCCGCAGACCGGTAAag         <  1:3069226/51‑1 (MQ=255)                                       cGTTAAAGAGTGTGCTGCCCGTACTGGCCGCAACCCGCAGACCGGTAAaga        <  1:4180672/51‑1 (MQ=255)                                       cGTTAAAGAGTGTGCTGCCCGTACTGGCCGCAACCCGCAGACCGGTAAaga        >  1:1319526/1‑51 (MQ=255)                                         ttAAAGAGTGTGCTGCCCGTACTGGCCGCAACCCGCAGACCGGTAAAgaga      >  1:1170244/1‑51 (MQ=255)                                          tAAAGAGTGTGCTGCCCGTACTGGCCGCAACCCGCAGACCGGTAAAGAGAt     <  1:4555544/51‑1 (MQ=255)                                            aaGAGTGTGCTGCCCGTACTGGCCGCAACCCGCAGACCGGTAAAGAGATCa   >  1:1317413/1‑51 (MQ=255)                                             agagTGTGCTGCCCGTACTGGCCGCAACCCGCAGACCGGTAAAGAGATCAc  >  1:6771611/1‑51 (MQ=255)                                             agagTGTGCTGCCCGTACTGGCCGCAACCCGCAGACCGGTAAAGAGATCAc  >  1:6700984/1‑51 (MQ=255)                                                 |                                               GGGGATGATGTAGCACTGGTAGGTTTTGGTACTTTTGCCGTTAAAGAGCGTGCTGCCCGTACTGGCCGCAACCCGCAGACCGGTAAAGAGATCAC  >  REL606/433620‑433714

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