Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	REL606	433,564	C→T	A20V (GCG→GTG)	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,564	0	C	T	100.0%	265.1 / NA	67	A20V (GCG→GTG)	hupB	HU, DNA‑binding transcriptional regulator, beta subunit
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (38/29);  total (38/29)

AAGAGGAAGAGAAGAGTGAATAAATCTCAATTGATCGACAAGATTGCTGCAGGGGCTGATATCTCTAAAGCTGCGGCTGGCCGTGCGTTAGATGCTATTATTGCTTCCGTAACTGAATCTCTGAAAGAAGGGGATGATGTAGCA  >  REL606/433491‑433634                                                                          |                                                                       aaGAGGAAGAGAAGAGTGAATAAATCTCAATTGATCGACAAGATTGCTGCAGGGGCTGATATCTCTAAAGCTGTgg                                                                      >  2:1384008/1‑76 (MQ=255) aaGAGGAAGAGAAGAGTGAATAAATCTCAATTGATCGACAAGATTGCTGCAGGGGCTGATATCTCTAAAGCTGTgg                                                                      <  2:1433870/76‑1 (MQ=255)  agagGAAGAGAAGAGTGAATAAATCTCAATTGATCGACAAGATTGCTGCAGGGGCTGATATCTCTAAAGCTGTGGc                                                                     >  2:1153615/1‑76 (MQ=255)  agagGAAGAGAAGAGTGAATAAATCTCAATTGATCGACAAGATTGCTGCAGGGGCTGATATCTCTAAAGCTGTGGc                                                                     >  2:534583/1‑76 (MQ=255)  agagGAAGAGAAGAGTGAATAAATCTCAATTGATCGACAAGATTGCTGCAGGGGCTGATATCTCTAAAGCTGTGGc                                                                     >  2:534592/1‑76 (MQ=255)       aagagaagagTGAATAAATCTCAATTGATCGACAAGATTGCTGCAGGGGCTGATATCTCTAAAGCTGTGGCTGGcc                                                                <  1:378443/76‑1 (MQ=255)       aagagaagagTGAATAAATCTCAATTGATCGACAAGATTGCTGCAGGGGCTGATATCTCTAAAGCTGTGGCTGGcc                                                                >  2:260246/1‑76 (MQ=255)       aagagaagagTGAATAAATCTCAATTGATCGACAAGATTGCTGCAGGGGCTGATATCTCTAAAGCTGTGGCTGGcc                                                                >  1:1815440/1‑76 (MQ=255)        agagaagagTGAATAAATCTCAATTGATCGACAAGATTGCTGCAGGGGCTGATATCTCTAAAGCTGTGGCTGGCCg                                                               >  1:228866/1‑76 (MQ=255)          agaagagTGAATAAATCTCAATTGATCGACAAGATTGCTGCAGGGGCTGATATCTCTAAAGCTGTGGCTGGCCGTg                                                             <  2:846408/76‑1 (MQ=255)            aagagTGAATAAATCTCAATTGATCGACAAGATTGCTGCAGGGGCTGATATCTCTAAAGCTGTGGCTGGCCGTGCg                                                           >  2:1644354/1‑76 (MQ=255)            aagagTGAATAAATCTCAATTGATCGACAAGATTGCTGCAGGGGCTGATATCTCTAAAGCTGTGGCTGGCCGTGCg                                                           >  2:691312/1‑76 (MQ=255)            aagagTGAATAAATCTCAATTGATCGACAAGATTGCTGCAGGGGCTGATATCTCTAAAGCTGTGGCTGGCCGTGCg                                                           >  2:641109/1‑76 (MQ=255)            aagagTGAATAAATCTCAATTGATCGACAAGATTGCTGCAGGGGCTGATATCTCTAAAGCTGTGGCTGGCCGTGCg                                                           >  1:908620/1‑76 (MQ=255)             agagTGAATAAATCTCAATTGATCGACAAGATTGCTGCAGGGGCTGATATCTCTAAAGCTGTGGCTGGCCGTGCGt                                                          >  1:1844201/1‑76 (MQ=255)              gagTGAATAAATCTCAATTGATCGACAAGATTGCTGCAGGGGCTGATATCTCTAAAGCTGTGGCTGGCCGTGCGtt                                                         <  2:713381/76‑1 (MQ=255)                gTGAATAAATCTCAATTGATCGACAAGATTGCTGCAGGGGCTGATATCTCTAAAGCTGTGGCTGGCCGTGCGTTAg                                                       <  1:529652/76‑1 (MQ=255)                gTGAATAAATCCCAATTGAACGACAAGATTGCTGCAGGGGCTGATATCTCTAAAGCTGTGGCTGGCCGTGCGTTAg                                                       <  2:1850236/76‑1 (MQ=255)                 tGAATAAATCTCAATTGATCGACAAGATTGCTGCAGGGGCTGATATCTCTAAAGCTGTGGCTGGCCGTGCGTTAGa                                                      >  2:1560329/1‑76 (MQ=255)                  gAATAAATCTCAATTGATCGACAAGATTGCTGCAGGGGCTGATATCTCTAAAGCTGTGGCTGGCCGTGCGTTAGAt                                                     >  1:19956/1‑76 (MQ=255)                     tAAATCTCAATTGATCGACACGATTGCTGCAGGGGCTGATATCTCTACAGCTGTGGCTGGCCGCGGGTTAGATGCt                                                  >  2:67213/1‑76 (MQ=255)                     tAAATCTCAATTGATCGACAAGATTGCTGCAGGGGCTGATATCTCTAAAGCTGTGGCTGGCCGTGCGTTAGATGCt                                                  >  1:203179/1‑76 (MQ=255)                      aaaTCTCAATTGATCGACAAGNGTGCTGCAGGGGCTGAGATCTCTAAAGCTGTGGCTGGCCGTGCGTTAGATGCta                                                 <  2:307352/76‑1 (MQ=255)                      aaaTCTCAATTGATCGACAAGATTGCTGCAGGGGCTGATATCTCTAAAGCTGTGGCTGGCCGTGCGTTAGATGCta                                                 >  2:975598/1‑76 (MQ=255)                       aaTCTCAATTGATCGACAAGATTGCTGCAGGGGCTGATATCTCTAAAGCTGTGGCTGGCCGTGCGTTAGATGCtat                                                <  1:1745081/76‑1 (MQ=255)                        aTCTCAATTGATCGACAAGATTGCTGCAGGGGCTGATATCTCTAAAGCTGTGGCTGGCCGTGCGTTAGATGCtatt                                               <  2:1184431/76‑1 (MQ=255)                        aTCTCAATTGATCGACAAGATTGCTGCAGGGGCTGATATCTCTAAAGCTGTGGCTGGCCGTGCGTTAGATGCtatt                                               <  1:1336320/76‑1 (MQ=255)                         tctcAATTGATCGACAAGATTGCTGCAGGGGCTGATATCTCTAAAGCTGTGGCTGGCCGTGCGTTAGATGCtatta                                              <  2:427090/76‑1 (MQ=255)                         tctcAATTGATCGACAAGATTGCTGCAGGGGCTGATATCTCTAAAGCTGTGGCTGGCCGTGCGTTAGATGCtatta                                              >  1:496455/1‑76 (MQ=255)                           tcAATTGATCGACAAGATTGCTGCAGGGGCTGATATCTCTAAAGCTGTGGCTGGCCGTGCGTTAGATGCTattatt                                            <  1:1242137/76‑1 (MQ=255)                            cAATTGATCGACAAGATTGCTGCAGGGGCTGATATCTCTAAAGCTGTGGCTGGCCGTGCGTTAGATGCTATTATTg                                           >  1:330046/1‑76 (MQ=255)                            cAATTGATCGACAAGATTGCTGCAGGGGCTGATATCTCTAAAGCTGTGGCTGGCCGTGCGTTAGATGCTATTATTg                                           >  1:669454/1‑76 (MQ=255)                            cAATTGATCGACAAGATTGCTGCAGGGGCTGATATCTCTAAAGCTGTGGCTGGCCGTGCGTTAGATGCTATTATTg                                           >  1:364225/1‑76 (MQ=255)                               ttGATCGACAAGATTGCTGCAGGGGCTGATATCTCTAAAGCTGTGGCTGGCCGTGCGTTAGATGCTATTATTGCtt                                        >  1:89551/1‑76 (MQ=255)                               ttGATCGACAAGATTGCTGCAGGGGCTGATATCTCTAAAGCTGTGGCTGGCCGTGCGTTAGATGCTATTATTGCtt                                        >  1:911684/1‑76 (MQ=255)                                 gATCGACAAGATTGCTGCAGGGGCTGATATCTCTAAAGCTGTGGCTGGCCGTGCGTTAGATGCTATTATTGCTTcc                                      <  2:178215/76‑1 (MQ=255)                                  aTCGACAAGATTGCTGCAGGGGCTGATATCTCTAAAGCTGTGGCTGGCCGTGCGTTAGATGCTATTATTGCTTCCg                                     >  1:1684758/1‑76 (MQ=255)                                   tCGACAAGATTGCTGCAGGGGCTGATATCTCTAAAGCTGTGGCTGGCCGTGCGTTAGATGCTATTATTGCTTCCgn                                    <  1:296576/76‑2 (MQ=255)                                   tCGACAAGATTGCTGCAGGGGCTGATATCTCTAAAGCTGTGGCTGGCCGTGCGTTAGATGCTATTATTGCTTCCGt                                    <  2:820900/76‑1 (MQ=255)                                     gACAAGATTGCTGCAGGGGCTGATATCTCTAAAGCTGTGGCTGGCCGTGCGTTAGATGCTATTATTGCTTCCGTaa                                  <  2:155099/76‑1 (MQ=255)                                     gACAAGATTGCTGCAGGGGCTGATATCTCTAAAGCTGTGGCTGGCCGTGCGTTAGATGCTATTATTGCTTCCGTaa                                  <  2:1326716/76‑1 (MQ=255)                                      aCAAGATTGCTGCAGGGGCTGATATCTCTAAAGCTGTGGCTGGCCGTGCGTTAGATGCTATTATTGCTTCCGTAAc                                 >  1:1130421/1‑76 (MQ=255)                                         aGATTGCTGCAGGGGCTGATATCTCTAAAGCTGTGGCTGGCCGTGCGTTAGATGCTATTATTGCTTCCGTAACTGa                              <  2:70790/76‑1 (MQ=255)                                         aGATTGCTGCAGGGGCTGATATCTCTAAAGCTGTGGCTGGCCGTGCGTTAGATGCTATTATTGCTTCCGTAACTGa                              <  1:677076/76‑1 (MQ=255)                                           aTTGCTGCAGGGGCTGATATCTCTAAAGCTGTGGCTGGCCGTGCGTTAGATGCTATTATTGCTTCCGTAACTGAAt                            >  1:408113/1‑76 (MQ=255)                                               ctgcAGGGGCTGATATCTCTAAAGCTGTGGCTGGCCGTGCGTTAGATGCTATTATTGCTTCCGTAACTGAATctct                        <  1:465857/76‑1 (MQ=255)                                                  cAGGGGCTGATATCTCTAAAGCTGTGGCTGGCCGTGCGTTAGATGCTATTATTGCTTCCGTAACTGAATCTCTGaa                     >  1:1520144/1‑76 (MQ=255)                                                  cAGGGGCTGATATCTCTAAAGCTGTGGCTGGCCGTGCGTTAGATGCTATTATTGCTTCCGTAACTGAATCTCTGaa                     >  2:181219/1‑76 (MQ=255)                                                   aGGGGCTGATATCTCTAAAGCTGTGGCTGGCCGTGCGTTAGATGCTATTATTGCTTCCGTAACTGAATCTCTGaaa                    >  2:157541/1‑76 (MQ=255)                                                      ggCTGATATCTCTAAAGCTGTGGCTGGCCGTGCGTTAGATGCTATTATTGCCTCCGTAACTGAATCTCTGAaagaa                 <  1:1827214/76‑1 (MQ=255)                                                        cTGATATCTCTAAAGCTGTGGCTGGCCGTGCGTTAGATGCTATTATTGCTTCCGTAACTGAATCTCTGAAAGAAgg               <  1:640345/76‑1 (MQ=255)                                                          gATATCTCTAAAGCTGTGGCTGGCCGTGCGTTAGATGCTATTATTGCTTCCGTAACTGAATCTCTGAAAGAAgggg             >  2:835619/1‑76 (MQ=255)                                                          gATATCTCTAAAGCTGTGGCTGGCCGTGCGTTAGATGCTATTATTGCTTCCGTAACTGAATCTCTGAAAGAAgggg             <  1:1088025/76‑1 (MQ=255)                                                          gATATCTCTAAAGCTGTGGCTGGCCGTGCGTTAGATGCTATTATTGCTTCCGTAACTGAATCTCTGAAAGAAgggg             >  1:1689098/1‑76 (MQ=255)                                                          gATATCTCTAAAGCTGTGGCTGGCCGTGCGTTAGATGCTATTATTGCTTCCGTAACTGAATCTCTGAAAGAAgggg             >  1:1777953/1‑76 (MQ=255)                                                            tatCTCTAAAGCTGTGGCTGGCCGTGCGTTAGATGCTATTATTGCTTCCGTAACTGAATCTCTGAAAGAAGGGgat           >  1:401881/1‑76 (MQ=255)                                                            tatCTCTAAAGCTGTGGCTGGCCGTGCGTTAGATGCTATTATTGCTTCCGTAACTGAATCTCTGAAAGAAGGGgat           >  1:215017/1‑76 (MQ=255)                                                             atCTCTAAAGCTGTGGCTGGCCGTGCGTTAGATGCTATTATTGCTTCCGTAACTGAATCTCTGAAAGAAGGGgatg          >  1:1198988/1‑76 (MQ=255)                                                             atCTCTAAAGCTGTGGCTGGCCGTGCGTTAGATGCTATTATTGCTTCCGTAACTGAATCTCTGAAAGAAGGGgatg          <  1:167934/76‑1 (MQ=255)                                                             atCTCTAAAGCTGTGGCTGGCCGTGCGTTAGATGCTATTATTGCTTCCGTAACTGAATCTCTGAAAGAAGGGgatg          >  1:864848/1‑76 (MQ=255)                                                              tctcTAAAGCTGTGGCTGGCCGTGCGTTAGATGCTATTATCGCTTCCGTAACTGAATCTCTGAAAGAAGGGgatga         >  1:457516/1‑76 (MQ=255)                                                               ctctAAAGCTGTGGCTGGCCGTGCGTTAGATGCTATTATTGCTTCCGTAACTGAATCTCTGAAAGAAGGGgatgat        <  1:237725/76‑1 (MQ=255)                                                                tctAAAGCTGTGGCTGGCCGTGCGTTAGATGCTATTATTGCTTCCGTAACTGAATCTCTGAAAGAAGGGGatgatg       >  1:854053/1‑76 (MQ=255)                                                                 ctAAAGCTGTGGCTGGCCGTGCGTTAGATGCTATTATTGCTTCCGTAACTGAATCTCTGAAAGAAGGGGATGATGt      >  2:519592/1‑76 (MQ=255)                                                                  tAAAGCTGTGGCTGGCCGTGCGTTAGATGCTATTATTGCTTCCGTAACTGAATCTCTGAAAGAAGGGGATGATGTa     <  2:1515655/76‑1 (MQ=255)                                                                    aaGCTGTGGCTGGCCGTGCGTTAGATGCTATTATTGCTTCCGTAACTGAATCTCTGAAAGAAGGGGATGATGTAGc   <  1:1089236/76‑1 (MQ=255)                                                                    aaGCTGTGGCTGGCCGTGCGTTAGATGCTATTATTGCTTCCGTAACTGAATCTCTGAAAGAAGGGGATGATGTAGc   <  1:1089218/76‑1 (MQ=255)                                                                     aGCTGTGGCTGGCCGTGCGTTAGATGCTATTATTGCTTCCGTAACTGAATCTCTGAAAGAAGGGGATGATGTAGCa  <  1:1533678/76‑1 (MQ=255)                                                                          |                                                                       AAGAGGAAGAGAAGAGTGAATAAATCTCAATTGATCGACAAGATTGCTGCAGGGGCTGATATCTCTAAAGCTGCGGCTGGCCGTGCGTTAGATGCTATTATTGCTTCCGTAACTGAATCTCTGAAAGAAGGGGATGATGTAGCA  >  REL606/433491‑433634

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: ‑