Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	REL606	433,524	A→T	I7F (ATC→TTC)	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,524	0	A	T	100.0%	223.3 / NA	56	I7F (ATC→TTC)	hupB	HU, DNA‑binding transcriptional regulator, beta subunit
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base T (27/29);  total (27/29)

GTTCGTACCTTGAAGGATTCAAGTGCGATATAAATTATAAAGAGGAAGAGAAGAGTGAATAAATCTCAATTGATCGACAAGATTGCTGCAGGGGCTGATATCTCTAAAGCTGCGGCTGGCCGTGCGTTAGATGCTATTATTGCTT  >  REL606/433452‑433596                                                                         |                                                                         gTTCGTACCTTGAAGGATTCAAGTGCGATATAAATTATAAAGAGGAAGAGAAGAGTGAATAAATCTCAATTGTTCg                                                                       >  2:1472017/1‑76 (MQ=255) gTTCGTACCTTGAAGGATTCAAGTGCGATATAAATTATAAAGAGGAAGAGAAGAGTGAATAAATCTCAATTGTTCg                                                                       <  1:1588804/76‑1 (MQ=255)  ttCGTACCTTGAAGGATTCAAGTGCGATATAAATTATAAAGAGGAAGATAAGAGTGAATAAATCTCAATTGTTCGa                                                                      <  2:1095106/76‑1 (MQ=255)  ttCGTACCTTGAAGGATTCAAGTGCGATATAAATTATAAAGAGGAAGAGAAGAGTGAATAAATCTCAATTGTTCGa                                                                      >  1:595121/1‑76 (MQ=255)  ttCGTACCTTGAAGGATTCAAGTGCGATATAAATTATAAAGAGGAAGAGAAGAGTGAATAAATCTCAATTGTTCGa                                                                      <  2:1114127/76‑1 (MQ=255)   tCGTACCTTGAAGGATTCAAGTGCGATATAAATTATAAAGAGGAAGAGAAGAGTGAATAAATCTCAATTGTTCGAc                                                                     >  1:240719/1‑76 (MQ=255)   tCGTACCTTGAAGGATTCAAGTGCGATATAAATTATAAAGAGGAAGAGAAGAGTGAATAAATCTCAATTGTTCGAc                                                                     >  1:1810370/1‑76 (MQ=255)      tACCTTGAAGGATTCAAGTGCGATATAAATTATAAAGAGGAAGAGAAGAGTGAATAAATCTCAATTGTTCGACAAg                                                                  >  1:952749/1‑76 (MQ=255)             aaGGATTCAAGTGCGATATAAATTATAAAGAGGAAGAGAAGAGTGAATAAATCTCAATTGTTCGACAAGATtgctg                                                           >  1:1480837/1‑76 (MQ=255)                gATTCAAGTGCGATATAAATTATAAAGAGGAAGAGAAGAGTGAATAAATCTCAATTGTTCGACAAGATTGCTGCAg                                                        <  1:1007371/76‑1 (MQ=255)                gATTCAAGTGCGATATAAATTATAAAGAGGAAGAGAAGAGTGAATAAATCTCAATTGTTCGACAAGATTGCTGCAg                                                        <  2:851093/76‑1 (MQ=255)                 aTTCAAGTGCGATATAAATTATAAAGAGGAAGAGAAGANTGAATAAATCTCAATTGTTCGACAAGATTGCTGCAgg                                                       >  1:1599702/1‑76 (MQ=255)                 aTTCAAGTGCGATATAAATTATAAAGAGGAAGAGAAGAGTGAATAAATCTCAATTGTTCGACAAGATTGCTGCAgg                                                       <  2:35459/76‑1 (MQ=255)                     aaGTGCGATATAAATTATAAAGAGGAAGAGAAGAGTGAATAAATCTCAATTGTTCGACAAGATTGCTGCAGGGGCt                                                   >  2:1528724/1‑76 (MQ=255)                        tGCGATATAAATTATAAAGAGGAAGAGAAGAGTGAATAAATCTCAATTGTTCGACAAGATTGCTGCAGGGGCTGat                                                >  2:532888/1‑76 (MQ=255)                          cGATATAAATTATAAAGAGGAAGAGAAGAGTGAATAAATCTCAATTGTTCGACAAGATTGCTGCAGGGACTGatat                                              >  1:1871252/1‑76 (MQ=255)                           ggnataAATTATAAAGAGGAAGAGAAGAGTGAATAAATCTCAATTGTTCGACAAGATTGCTGCAGGGGCTGATAtc                                             >  2:1412831/4‑76 (MQ=255)                             tataAATTATAAAGAGGAAGAGAAGAGTGAATAAATCTCAATTGTTCGACAAGATTGCTGCAGGGGCTGATAtctc                                           <  2:13861/76‑1 (MQ=255)                             tataAATTATAAAGAGGAAGAGAAGAGTGAATAAATCTCAATTGTTCGACAAGATTGCTGCAGGGGCTGATAtctc                                           >  2:13125/1‑76 (MQ=255)                              ataAATTATAAAGAGGAAGAGAAGAGTGAATAAATCTCAATTGTTCGACAAGATTGCTGCAGGGGCTGATATctct                                          >  1:1248524/1‑76 (MQ=255)                              ataAATTATAAAGAGGAAGAGAAGAGTGAATAAATCTCAATTGTTCGACAAGATTGCTGCAGGGGCTGATATctct                                          <  1:217548/76‑1 (MQ=255)                                aaaTTATAAAGAGGAAGAGAAGAGTGAATAAATCTCAATTGTTCGACAAGATTGCTGCAGGGGCTGATATCTCTaa                                        >  2:995960/1‑76 (MQ=255)                                aaaTTATAAAGAGGAAGAGAAGAGTGAATAAATCTCAATTGTTCGACAAGATTGCTGCAGGGGCTGATATCTCTaa                                        <  1:361842/76‑1 (MQ=255)                                 aaTTATAAAGAGGAAGAGAAGAGTGAATAAATCTCAATTGTTCGACAAGATTGCTGCAGGGGCTGATATCTCTaaa                                       <  2:1984282/76‑1 (MQ=255)                                  aTTATAAAGAGGAAGAGAAGAGTGAATAAATCTCAATTGTTCGACAAGATTGCTGCAGGGGCTGATATCTCTAAAg                                      >  1:892682/1‑76 (MQ=255)                                  aTTATAAAGAGGAAGAGAAGAGAGAATAAATCTCAATTGTTCGACAAGATTGCTGCAGGGGCTGATATCTCTAAAg                                      >  1:1609487/1‑76 (MQ=255)                                    tataAAGAGGAAGAGAAGAGTGAATAAATCTCAATTGTTCGACAAGATTGCTGCAGGGGCTGATATCTCTAAAGCt                                    <  2:739088/76‑1 (MQ=255)                                      taAAGAGGAAGAGAAGAGTGAATAAATCTCAATTGTTCGACAAGATTGCTGCAGGGGCTGATATCTCTAAAGCTGc                                  <  2:1690155/76‑1 (MQ=255)                                       aaaGAGGAAGAGAAGAGTGAATAAATCTCAATTGTTCGACAAGATTGCTGCAGGGGCTGATATCTCTAAAGCTGCg                                 >  2:1704338/1‑76 (MQ=255)                                       aaaGAGGAAGAGAAGAGTGAATAAATCTCAATTGTTCGACAAGATTGCTGCAGGGGCTGATATCTCTAAAGCTGCg                                 >  2:270740/1‑76 (MQ=255)                                        aaGAGGAAGAGAAGAGTGAATAAATCTCAATTGTTCGACAAGATTGCTGCAGGGGCTGATATCTCTAAAGCTGCgg                                <  1:1042800/76‑1 (MQ=255)                                         agagGAAGAGAAGAGTGAATAAATCTCAATTGTTCGACAAGATTGCTGCAGGGGCTGATATCTCTAAAGCTGCGGc                               >  2:113110/1‑76 (MQ=255)                                          gagGAAGAGAAGAGTGAATAAATCTCAATTGTTCGACAAGATTGCTGCAGGGGCTGATATCTCTAAAGCTGCGGCt                              <  1:1742833/76‑1 (MQ=255)                                           agGAAGAGAAGAGTGAATAAATCTCAATTGTTCGACAAGATTGCTGCAGGGGCTGATATCTCTAAAGCTGCGGCTg                             >  2:22542/1‑76 (MQ=255)                                           agGAAGAGAAGAGTGAATAAATCTCAATTGTTCGACAAGATTGCTGCAGGGGCTGATATCTCTAAAGCTGCGGCTg                             <  1:5377/76‑1 (MQ=255)                                                gagaagagTGAATAAATCTCAATTGTTCGACAAGATTGCTGCAGGGGCTGATATCTCTAAAGCTGCGGCTGGCCGt                        >  1:1532944/1‑76 (MQ=255)                                                gagaagagTGAATAAATCTCAATTGTTCGACAAGATTGCTGCAGGGGCTGATATCTCTAAAGCTGCGGCTGGCCGt                        <  2:163634/76‑1 (MQ=255)                                                 agaagagTGAATAAATCTCAATTGTTCGACAAGATTGCTGCAGGGGCTGATATCTCTAAAGCTGCGGCTGGCCGTg                       >  2:107995/1‑76 (MQ=255)                                                    agagTGAATAAATCTCAATTGTTCGACAAGATTGCTGCAGGGGCTGATATCTCTAAAGCTGCGGCTGGCCGTGCGt                    <  2:771430/76‑1 (MQ=255)                                                    agagTGAATAAATCTCAATTGTTCGACAAGATTGCTGCAGGGGCTGATATCTCTAAAGCTGCGGCTGGCCGTGCGt                    <  1:1062500/76‑1 (MQ=255)                                                      agTGAATAAATCTCAATTGTTCGACAAGATTGCTGCAGGGGCTGATATCTCTAAAGCTGCGGCTGGCCGTGCGTTa                  <  2:1146405/76‑1 (MQ=255)                                                       gTGAATAAATCTCAATTGTTCGACAAGATTGCTGCAGGGGCTGATATCTCTAAAGCTGCGGCTGGCCGTGCGTTAg                 <  1:554433/76‑1 (MQ=255)                                                       gTGAATAAATCTCAATTGTTCGACAAGATTGCTGCAGGGGCTGATATCTCTAAAGCTGCGGCTGGCCGTGCGTTAg                 <  1:180638/76‑1 (MQ=255)                                                        tGAATAAATCTCAATTGTTCGACAAGATTGCTGCAGGGGCTGATATCTCTAAAGCTGCGGCTGGCCGTGCGTTAGa                >  2:1242552/1‑76 (MQ=255)                                                         gAATAAATCTCAATTGTTCGACAAGATTGCTGCAGGGGCTGATATCTCTAAAGCTGCGGCTGGCCGTGCGTTAGAt               >  2:309930/1‑76 (MQ=255)                                                         gAATAAATCTCAATTGTTCGACAAGATTGCTGCAGGGGCTGATATCTCTAAAGCTGCGGCTGGCCGTGCGGTAGAt               >  1:1110598/1‑76 (MQ=255)                                                            tAAATCTCAATTGTTCGACAAGATTGCTGCAGGGGCTGATATCTCTAAAGCTGCGGCTGGCCGTGCGTTAGATGCt            <  2:240719/76‑1 (MQ=255)                                                              aaTCTCAATTGTTCGACAAGATTGCTGCAGGGGCTGATATCTCTAAAGCTGCGGCTGGCCGTGCGTTAGATGCtat          <  1:1354559/76‑1 (MQ=255)                                                              aaTCTCAATTGTTCGACAAGATTGCTGCAGGGGCTGATATCTCTAAAGCTGCGGCTGGCCGTGCGTTAGATGCtat          <  1:461257/76‑1 (MQ=255)                                                               aTCTCAATTGTTCGACAAGATTGCTGCAGGGGCTGATATCTCTAAAGCTGCGGCTGGCCGTGCGTTAGATGCtatt         <  1:1892869/76‑1 (MQ=255)                                                               aTCTCAATTGTTCGACAAGATTGCTGCAGGGGCTGATATCTCTAAAGCTGCGGCTGGCCGTGCGTTAGATGCtatt         <  1:870261/76‑1 (MQ=255)                                                                tctcAATTGTTCGACAAGATTGCTGCAGGGGCTGATATCTCTAAAGCTGCGGCTGGCCGTGCGTTAGATGCtatta        <  2:298458/76‑1 (MQ=255)                                                                   cAATTGTTCGACAAGATTGCTGCAGGGGCTGATATCTCTAAAGCTGCGGCTGGCCGTGCGTTAGATGCTATTATTg     <  1:1370445/76‑1 (MQ=255)                                                                   cAATTGTTCGACAAGATTGCTGCAGGGGCTGATATCTCTAAAGCTGCGGCTGGCCGTGCGTTAGATGCTATTATTg     <  1:127963/76‑1 (MQ=255)                                                                    aaTTGTTCGACAAGATTGCTGCAGGGGCTGATATCTCTAAAGCTGCGGCTGGCCGTGCGTTAGATGCTATTATTGc    >  2:539714/1‑76 (MQ=255)                                                                     aTTGTTCGACAAGATTGCTGCAGGGGCTGATATCTCTAAAGCTGCGGCTGGCCGTGCGTTAGATGCTATTATTGCt   >  1:646805/1‑76 (MQ=255)                                                                      ttGTTCGACAAGATTGCTGCAGGGGCTGATATCTCTAAAGCTGCGGCTGGCCGTGCGTTAGATGCTATTATTGCtt  >  2:253923/1‑76 (MQ=255)                                                                         |                                                                         GTTCGTACCTTGAAGGATTCAAGTGCGATATAAATTATAAAGAGGAAGAGAAGAGTGAATAAATCTCAATTGATCGACAAGATTGCTGCAGGGGCTGATATCTCTAAAGCTGCGGCTGGCCGTGCGTTAGATGCTATTATTGCTT  >  REL606/433452‑433596

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