Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	REL606	433,711	T→G	I69S (ATC→AGC)	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,711	0	T	G	100.0%	146.0 / NA	37	I69S (ATC→AGC)	hupB	HU, DNA‑binding transcriptional regulator, beta subunit
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base G (25/12);  total (25/12)

GTAGGTTTTGGTACTTTTGCCGTTAAAGAGCGTGCTGCCCGTACTGGCCGCAACCCGCAGACCGGTAAAGAGATCACCATCGCTGCTGCTAAAGTACCGAGCTTCCGTGCAGGTAAAGCACTGAAAGACGCGGTAAACTAAG  >  REL606/433638‑433779                                                                          |                                                                     gTAGGTTTTGGTACTTTTGCCGTTAAAGAGCGTGCTGCCCGTACTGGCCGCAACCCGCAGACCGGTAAAGAGAGCa                                                                    <  2:481602/76‑1 (MQ=255) gTAGGTTTTGGTACTTTTGCCGTTAAAGAGCGTGCTGCCCGTACTGGCCGCAACCCGCAGACCGGTAAAGAGAGCa                                                                    >  2:1041711/1‑76 (MQ=255)   aGGTTTTGGTACTTTTGCCGTTAAAGAGCGTGCTGCCCGTACTGGCCGCAACCCGCAGACCGGTAAAGAGAGCAcc                                                                  >  2:983705/1‑76 (MQ=255)      ttttGGTACTTTTGCCGTTAAAGAGCGTGCTGCCCGTACTGGCCGCAACCCGCAGACCGGTAAAGAGAGCACCATc                                                               <  1:878465/76‑1 (MQ=255)          ggTACTTTTGCCGTTAAAGAGCGTGCTGCCCGTACTGGCCGCAACCCGCAGACCGGTAAAGAGAGCACCATCgctg                                                           <  2:1273149/76‑1 (MQ=255)          ggTACTTTTGCCGTTAAAGAGCGTGCTGCCCGTACTGGCCGCAACCCGCAGACCGGTAAAGAGAGCACCATCgctg                                                           >  2:926233/1‑76 (MQ=255)          ggTACTTTTGCCGTTAAAGAGCGTGCTGCCCGTACTGGCCGCAACCCGCAGACCGGTAAAGAGAGCACCATCgctg                                                           <  2:510511/76‑1 (MQ=255)            tACNTTTGCCGTTAAAGAGCGTGCTGCCCGTACTGGCCGCAACCCGCAGACCGGTAAAGAGAGCACCATCgctgct                                                         >  1:1200335/1‑76 (MQ=255)              cTTTTGCCGTTAAAGAGCGTGCTGCCCGTACTGGCCGCAACCCGCAGACCGGTAAAGAGAGCACCATCgctgctgc                                                       >  2:163153/1‑76 (MQ=255)               ttttGCCGTTAAAGAGCGTGCTGCCCGTACTGGCCGCAACCCGCAGACCGGTAAAGAGAGCACCATCgctgctgct                                                      >  2:894575/1‑76 (MQ=255)               ttttGCCGTTAAAGAGCGTGCTGCCCGTACTGGCCGCAACCCGCAGACCGGTAAAGAGAGCACCATCgctgctgct                                                      >  1:610561/1‑76 (MQ=255)                 ttGCCGTTAAAGAGCGTGCTGCCCGGACTGGCCGCTGCCCGCAGACGGGTAAAGAAAGGACACGCGCTGCTGCTaa                                                    >  2:20026/1‑76 (MQ=255)                   gCCGTTAAAGAGCGTGCTGCCCGTACTGGCCGCAACCCGCAGACCGGTAAAGAGAGCACCATCGCTGCTGCTAAAg                                                  >  1:623373/1‑76 (MQ=255)                      gTTAAAGAGCGTGCTGCCCGTACTGGCCGCAACCCGCAGACCGGTAAAGAGAGCACCATCGCTGCTGCTAAAGTAc                                               >  1:25583/1‑76 (MQ=255)                        tAAAGAGCGTGCTGCCCGTACTGGCCGCATCCCGCAGACCGGTAAAGAGAGCACCATCGCTGCTGCTAAAGTACCg                                             <  2:36898/76‑1 (MQ=37)                        tAAAGAGCGTGCTGCCCGTACTGGCCGCAACCCGCAGACCGGTAAAGAGAGCACCATCGCTGCTGCTAAAGTACCg                                             >  2:335171/1‑76 (MQ=255)                        tAAAGAGCGTGCTGCCCGTACTGGCCGCAACCCGCAGACCGGTAAAGAGAGCACCATCGCTGCTGCTAAAGTACCg                                             >  1:306925/1‑76 (MQ=255)                         aaaGAGCGTGCTGCCCGTACTGGCAGAAACCCGCTGACCTGAACAGAGGGCTCCCTCGCTCGCGATAAAGTACCGa                                            >  1:1344781/1‑76 (MQ=255)                          aaGAGCGTGCTGCCCGTACTGGCCGCAACCCGCAGACCGGTAAAGAGAGCACCATCGCTGCTGCTAAAGTACCGAg                                           <  2:596356/76‑1 (MQ=255)                              gCGTGCTGCCCGTACTGGCCGCAACCCGCAGACCGGTAAAGAGAGCACCATCGCTGCTGCTAAAGTACCGAGCTTc                                       >  2:1259770/1‑76 (MQ=255)                               cGTGCTGCCCNAGGTGGCCGCNACCCGCAGACCGATAAAGAGAGCACCATCGCTGCTGCTAAAGTACCGAGCTTcc                                      >  2:1324109/1‑76 (MQ=255)                                gTGCTGCCCGTACTGGCCGCAACCCGCAGACCGGTAAAGAGAGCACCATCGCTGCTGCTAAAGTACCGAGCTTCCg                                     >  1:897377/1‑76 (MQ=255)                                 tgctgcCCGTACTGGCCGCAACCCGCAGACCGGTAAAGAGAGCACCATNACTGCTGCTAAAGTACCGAGCTTCCGt                                    <  2:4917/76‑1 (MQ=255)                                  gctgcCCGTACTGGCCGCAACCCGCAGACTGGTAAAGAGAGCACCATCGCTGCTGCTAAAGTACCGAGCTTCCGTg                                   >  2:765776/1‑76 (MQ=255)                                    tgcCCGTACTGGCCGCAACCCGCAGACCGGTAAAGAGAGCACCATCGCTGCTGCTAAAGTACCGAGCTTCCGTGCa                                 >  1:564759/1‑76 (MQ=255)                                    tgcCCGTACTGGCCGCAACCCGCAGACCGGTAAAGAGAGCACCATCGCTGCTGCTAAAGTACCGAGCTTCCGTGCa                                 <  2:138435/76‑1 (MQ=255)                                         gTACTGGCCGCAACCCGCAGACCGGTAAAGAGAGCACCATCGCTGCTGCTAAAGTACCGAGCTTCCGTGCAGGTaa                            >  2:710335/1‑76 (MQ=255)                                         gTACTGGCCGCAACCCGCAGACCGGTAAAGAGAGCACCATCGCTGCTGCTAAAGTACCGAGCTTCCGTGCAGGTaa                            >  1:323458/1‑76 (MQ=255)                                           aCTGGCCGCAACCCGCAGACCGGTAAAGAGAGCACCATCGCTGCTGCTAAAGTACCGAGCTTCCGTGCCGGTAAAg                          >  2:240959/1‑76 (MQ=255)                                             tGGCCGGAAGCCGCAGATCGGTAAAGAGAGCACCATCGCTGCTGCTAAAGTACCGAGCTTCNGTNCAGGTAAAGCa                        >  2:206355/1‑76 (MQ=255)                                             tGGCCGCAACCCGCAGACCGGTAAAGAGAGCACCATCGCTGCTGCTAAAGTACCGAGCTTCCGTGCAGGTAAAGCa                        >  1:897337/1‑76 (MQ=255)                                              ggCCGCAACCCGCAGACCGGTAAAGAGAGCACCATCGCTGCTGCTAAAGTACCGAGCTTCCGTGCAGGTAAAGCAc                       <  2:924290/76‑1 (MQ=255)                                                       ccGCAGACCGGTAAAGAGAGCACCATCGCTGCTGCTAAAGTACCGAGCTTCCGTGCAGGTAAAGCACTGAAAGAcg              >  1:992232/1‑76 (MQ=255)                                                        cGCAGACCGGTAAAGAGAGCACCATCGCTGCTGCTAAATTACCGAGCTTCCGTGCAGGTAAAGCACTGAAAGAcgc             <  1:1122309/76‑1 (MQ=255)                                                          cAGACCGGTAAAGAGAGCACCATCGCTGCTGCTAAAGTACCGAGCTTCCGTGCAGGTAAAGNACTGAAAGACGCgg           <  1:3586/76‑1 (MQ=255)                                                          cAGACCGGTAAAGAGAGCACCATCGCTGCTGCTAAAGTACCGAGCTTCCGTGCAGGTAAAGCACTGAAAGACGCgg           >  2:1038726/1‑76 (MQ=255)                                                            gACCGGTAAAGAGAGCACCATCGCTGCTGCTAAAGTACCGAGCTTCCGTGCAGGTAAAGCACTGAAAGACGCGGTa         >  1:1061243/1‑76 (MQ=255)                                                                 gTAAAGAGAGCACCATCGCTGCTGCTAAAGTACCGAGCTTCCGTGCAGGTAAAGCACTGAAAGACGCGGTAAACTa    >  1:599094/1‑76 (MQ=255)                                                                   aaaGAGAGCACCATCGCTGCTGCTAAAGTACCGAGCTTCCGTGCAGGTAAAGCACTGAAAGACGCGGTAAACTAAg  <  2:49279/76‑1 (MQ=255)                                                                          |                                                                     GTAGGTTTTGGTACTTTTGCCGTTAAAGAGCGTGCTGCCCGTACTGGCCGCAACCCGCAGACCGGTAAAGAGATCACCATCGCTGCTGCTAAAGTACCGAGCTTCCGTGCAGGTAAAGCACTGAAAGACGCGGTAAACTAAG  >  REL606/433638‑433779

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