Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	REL606	1,916,109	A→G	intergenic (+12/‑116)	yebK → / → pykA	predicted DNA‑binding transcriptional regulator/pyruvate kinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	REL606	1,916,109	0	A	G	100.0%	236.2 / NA	62	intergenic (+12/‑116)	yebK/pykA	predicted DNA‑binding transcriptional regulator/pyruvate kinase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (27/35);  total (27/35)

GATAAGCAGTTACTTAATTTAAGTGACGATCGCTAAAAACGACTGTCACTGTCCTAATCTTATACGACATCCGAATGAGATTAATTTATCGCCAT  >  REL606/1916062‑1916156                                                |                                                gATAAGCAGTTACTTAATTTAAGTGACGATCGCTAAAAACGactgtcgctg                                              <  1:4138029/51‑1 (MQ=255) gATAAGCAGTTACTTAATTTAAGTGACGATCGCTAAAAACGactgtcgctg                                              >  1:1668581/1‑51 (MQ=255)  aTAAGCAGTTACTTAATTTAAGTGACGATCGCTAAAAACGactgtcgctgt                                             >  1:3403427/1‑51 (MQ=255)  aTAAGCAGTTACTTAATTTAAGTGACGATCGCTAAAAACGactgtcgctgt                                             >  1:6237519/1‑51 (MQ=255)  aTAAGCAGTTACTTAATTTAAGTGACGATCGCTAAAAACGactgtcgctgt                                             <  1:6952239/51‑1 (MQ=255)   tAAGCAGTTACTTAATTTAAGTGACGATCGCTAAAAACGactgtcgctgtc                                            <  1:6773472/51‑1 (MQ=255)    aaGCAGTTACTTAATTTAAGTGACGATCGCTAAAAACGactgtcgctgt                                             <  1:3949991/49‑1 (MQ=255)    aaGCAGTTACTTAATTTAAGTGACGATCGCTAAAAACGACTGTCGCTGTcc                                           >  1:6569683/1‑51 (MQ=255)     aGCAGTTACTTAATTTAAGTGACGATCGCTAAAAACGACTGTCGCTGTCCt                                          <  1:3716047/51‑1 (MQ=255)     aGCAGTTACTTAATTTAAGTGACGATCGCTAAAAACGACTGTCGCTGTCCt                                          <  1:3633398/51‑1 (MQ=255)       cAGTTACTTAATTTAAGTGACGATCGCTAAAAACGACTGTCGCTGTCCTaa                                        >  1:4610596/1‑51 (MQ=255)       cAGTTACTTAATTTAAGTGACGATCGCTAAAAACGACTGTCGCTGTCCTaa                                        >  1:6706043/1‑51 (MQ=255)        aGTTACTTAATTTAAGTGACGATCGCTAAAAACGACTGTCGCTGTCCTAAt                                       >  1:1184595/1‑51 (MQ=255)          ttACTTAATTTAAGTGACGATCGCTAAAAACGACTGTCGCTGTCCTAATCt                                     <  1:6449959/51‑1 (MQ=255)          ttACTTAATTTAAGTGACGATCGCTAAAAACGACTGTCGCTGTCCTAATCt                                     <  1:3222095/51‑1 (MQ=255)          ttACTTAATTTAAGTGACGATCGCTAAAAACGACTGTCGCTGTCCTAATCt                                     >  1:1093415/1‑51 (MQ=255)            aCTTAATTTAAGTGACGATCGCTAAAAACGACTGTCGCTGTCCTAATCTta                                   <  1:3329374/51‑1 (MQ=255)            aCTTAATTTAAGTGACGATCGCTAAAAACGACTGTCGCTGTCCTAATCTta                                   >  1:817644/1‑51 (MQ=255)            aCTTAATTTAAGTGACGATCGCTAAAAACGACTGTCGCTGTCCTAATCTta                                   >  1:4568166/1‑51 (MQ=255)             cTTAATTTAAGTGACGATCGCTAAAAACGACTGTCGCTGTCCTAATCTtat                                  >  1:1572731/1‑51 (MQ=255)              ttAATTTAAGTGACGATCGCTAAAAACGACTGTCGCTGTCCTAATCTtata                                 <  1:6886091/51‑1 (MQ=255)              ttAATTTAAGTGACGATCGCTAAAAACGACTGTCGCTGTCCTAATCTtata                                 >  1:1197804/1‑51 (MQ=255)               tAATTTAAGTGACGATCGCTAAAAACGACTGTCGCTGTCCTAATCTTATAc                                >  1:1680300/1‑51 (MQ=255)               tAATTTAAGTGACGATCGCTAAAAACGACTGTCGCTGTCCTAATCTTATAc                                <  1:6764802/51‑1 (MQ=255)               tAATTTAAGTGACGATCGCTAAAAACGACTGTCGCTGTCCTAATCTTATAc                                <  1:2873792/51‑1 (MQ=255)                aaTTTAAGTGACGATCGCTAAAAACGACTGTCGCTGTCCTAATCTTATACg                               >  1:5229148/1‑51 (MQ=255)                aaTTTAAGTGACGATCGCTAAAAACGACTGTCGCTGTCCTAATCTTATACg                               <  1:1905398/51‑1 (MQ=255)                 aTTTAAGTGACGATCGCTAAAAACGACTGTCGCTGTCCTAATCTTATACGa                              <  1:5640485/51‑1 (MQ=255)                 aTTTAAGTGACGATCGCTAAAAACGACTGTCGCTGTCCTAATCTTATACGa                              >  1:3720037/1‑51 (MQ=255)                  tttAAGTGACGATCGCTAAAAACGACTGTCGCTGTCCTAATCTTATACGAc                             <  1:4038471/51‑1 (MQ=255)                  tttAAGTGACGATCGCTAAAAACGACTGTCGCTGTCCTAATCTTATACGAc                             >  1:69981/1‑51 (MQ=255)                   ttAAGTGACGATCGCTAAAAACGACTGTCGCTGTCCTAATCTTATACGACa                            <  1:5759295/51‑1 (MQ=255)                   ttAAGTGACGATCGCTAAAAACGACTGTCGCTGTCCTAATCTTATACGACa                            <  1:1185371/51‑1 (MQ=255)                   ttAAGTGACGATCGCTAAAAACGACTGTCGCTGTCCTAATCTTATACGACa                            >  1:6575930/1‑51 (MQ=255)                   ttAAGTGACGATCGCTAAAAACGACTGTCGCTGTCCTAATCTTATACGACa                            >  1:6804774/1‑51 (MQ=255)                   ttAAGTGACGATCGCTAAAAACGACTGTCGCTGTCCTAATCTTATACGACa                            <  1:91188/51‑1 (MQ=255)                    tAAGTGACGATCGCTAAAAACGACTGTCGCTGTCCTAATCTTATACGACAt                           <  1:1520426/51‑1 (MQ=255)                    tAAGTGACGATCGCTAAAAACGACTGTCGCTGTCCTAATCTTATACGACAt                           >  1:3379292/1‑51 (MQ=255)                       gTGACGATCGCTAAAAACGACTGTCGCTGTCCTAATCTTATACGACATCCg                        <  1:204315/51‑1 (MQ=255)                         gACGATCGCTAAAAACGACTGTCGCTGTCCTAATCTTATACGACATCCGaa                      <  1:1611855/51‑1 (MQ=255)                         gACGATCGCTAAAAACGACTGTCGCTGTCCTAATCTTATACGACATCCGaa                      <  1:2781167/51‑1 (MQ=255)                          aCGATCGCTAAAAACGACTGTCGCTGTCCTAATCTTATACGACATCCGAAt                     <  1:6731845/51‑1 (MQ=255)                           cGATCGCTAAAAACGACTGTCGCTGTCCTAATCTTATACGACATCCGAATg                    <  1:6856472/51‑1 (MQ=255)                            gATCGCTAAAAACGACTGTCGCTGTCCTAATCTTATACGACATCCGAATga                   <  1:411733/51‑1 (MQ=255)                            gATCGCTAAAAACGACTGTCGCTGTCCTAATCTTATACGACATCCGAATga                   <  1:6484999/51‑1 (MQ=255)                            gATCGCTAAAAACGACTGTCGCTGTCCTAATCTTATACGACATCCGAATga                   <  1:192138/51‑1 (MQ=255)                              tCGCTAAAAACGACTGTCGCTGTCCTAATCTTATACGACATCCGAATgaga                 >  1:5793092/1‑51 (MQ=255)                              tCGCTAAAAACGACTGTCGCTGTCCTAATCTTATACGACATCCGAATgaga                 >  1:2568732/1‑51 (MQ=255)                                gCTAAAAACGACTGTCGCTGTCCTAATCTTATACGACATCCGAATGAGAtt               <  1:2890752/51‑1 (MQ=255)                                gCTAAAAACGACTGTCGCTGTCCTAATCTTATACGACATCCGAATGAGAtt               <  1:2487856/51‑1 (MQ=255)                                 cTAAAAACGACTGTCGCTGTCCTAATCTTATACGACATCCGAATGAGATTa              <  1:4269308/51‑1 (MQ=255)                                  tAAAAACGACTGTCGCTGTCCTAATCTTATACGACATCCGAATGAGATTaa             >  1:2948718/1‑51 (MQ=255)                                  tAAAAACGACTGTCGCTGTCCTAATCTTATACGACATCCGAATGAGATTaa             >  1:1830470/1‑51 (MQ=255)                                  tAAAAACGACTGTCGCTGTCCTAATCTTATACGACATCCGAATGAGATTaa             <  1:2078019/51‑1 (MQ=255)                                   aaaaaCGACTGTCGCTGTCCTAATCTTATACGACATCCGAATGAGATTAAt            >  1:2354615/1‑51 (MQ=255)                                    aaaaCGACTGTCGCTGTCCTAATCTTATACGACATCCGAATGAGATTAAtt           <  1:2371441/51‑1 (MQ=255)                                    aaaaCGACTGTCGCTGTCCTAATCTTATACGACATCCGAATGAGATTAAtt           >  1:6162286/1‑51 (MQ=255)                                     aaaCGACTGTCGCTGTCCTAATCTTATACGACATCCGAATGAGATTAAttt          <  1:7018357/51‑1 (MQ=255)                                     aaaCGACTGTCGCTGTCCTAATCTTATACGACATCCGAATGAGATTAAttt          >  1:3824374/1‑51 (MQ=255)                                      aaCGACTGTCGCTGTCCTAATCTTATACGACATCCGAATGAGATTAATTTa         <  1:6205891/51‑1 (MQ=255)                                      aaCGACTGTCGCTGTCCTAATCTTATACGACATCCGAATGAGATTAATTTa         >  1:2085181/1‑51 (MQ=255)                                      aaCGACTGTCGCTGTCCTAATCTTATACGACATCCGAATGAGATTAATTTa         <  1:2209145/51‑1 (MQ=255)                                       aCGACTGTCGCTGTCCTAATCTTATACGACATCCGAATGAGATTAatatat        <  1:5267719/51‑1 (MQ=255)                                       aCGACTGTCGCTGTCCTAATCTTATACGACATCCGAATGAGATTAATTTAt        >  1:6336874/1‑51 (MQ=255)                                       aCGACTGTCGCTGTCCTAATCTTATACGACATCCGAATGAGATTAATTTAt        <  1:2405661/51‑1 (MQ=255)                                        cGACTGTCGCTGTCCTAATCTTATACGACATCCGAATGAGATTAATTTATc       <  1:6583681/51‑1 (MQ=255)                                        cGACTGTCGCTGTCCTAATCTTATACGACATCCGAATGAGATTAATTTATc       >  1:7110818/1‑51 (MQ=255)                                         gACTGTCGCTGTCCTAATCTTATACGACATCCGAATGAGATTAATTTATCg      >  1:2388572/1‑51 (MQ=255)                                         gACTGTCGCTGTCCTAATCTTATACGACATCCGAATGAGATTAATTTATCg      <  1:3396759/51‑1 (MQ=255)                                           ctgtcgctgtcCTAATCTTATACGACATCCGAATGAGATTAATTTATCGcc    <  1:7241975/51‑1 (MQ=255)                                            tgtcgctgtcCTAATCTTATACGACATCCGAATGAGATTAATTTAt        <  1:6330077/46‑1 (MQ=255)                                            tgtcgctgtcCTAATCTTATACGACATCCGAATGAGATTAATTTATCGCCa   >  1:6944490/1‑51 (MQ=255)                                            tgtcgctgtcCTAATCTTATACGACATCCGAATGAGATTAATTTATCGCCa   <  1:6750754/51‑1 (MQ=255)                                            tgtcgctgtcCTAATCTTATACGACATCCGAATGAGATTAATTTATCGCCa   <  1:7117019/51‑1 (MQ=255)                                             gtcgctgtcCTAATCTTATACGACATCCGAATGAGATTAATTTATCGCCAt  >  1:1471898/1‑51 (MQ=255)                                                |                                                GATAAGCAGTTACTTAATTTAAGTGACGATCGCTAAAAACGACTGTCACTGTCCTAATCTTATACGACATCCGAATGAGATTAATTTATCGCCAT  >  REL606/1916062‑1916156

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