Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	REL606	3,893,551	+G	intergenic (+6/‑50)	kup → / → insJ‑5	potassium transporter/IS150 hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	REL606	3,893,550	1	.	G	100.0%	249.1 / NA	65	intergenic (+5/‑51)	kup/insJ‑5	potassium transporter/IS150 hypothetical protein
Reads supporting (aligned to +/- strand):  ref base . (0/0);  new base G (23/42);  total (23/42)

GCGTGCACCAGATCAATTTGAAATCCCGCCAAACAGGGTTATCGAACTGGGTACTCAGGTCGAAATCTAACGCCA‑GACTGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCTAATATGTCAAAGCCAAAATACCC  >  REL606/3893476‑3893620                                                                            |                                                                       gCGTGCACCAGATCAATTTGAAATCCCGCCAAACAGGGTTATCGAACTGGGTACTCAGGTCGAAATCTAACGCCA‑g                                                                       >  2:1571632/1‑76 (MQ=255)  cGTGCACCAGATCAATTTGAAATCCCGCCAAACAGGGTTATCGAACTGGGTACTCAGGTCGAAATCTAACGCCA‑gg                                                                      <  1:1232689/76‑2 (MQ=255)    tGCACCAGATCAATTTGAAATCCCGCCNAACAGGGTAATCGAACTGGGTACTCAGGTCGAAATCTAACGCCA‑ggac                                                                    >  1:1554433/1‑73 (MQ=255)       aCCAGATCAATTTGAAATCCCGCCAAACAGGGTTATCGAACTGGGTACTCAGGTCGAAATCTAACGCCAGGACTGt                                                                  <  1:116581/76‑1 (MQ=255)       aCCAGATCAATTTGAAATCCCGCCAAACAGGGTTATCGAACTGGGTACTCAGGTCGAAATCTAACGCCAGGACTGt                                                                  >  1:1305453/1‑76 (MQ=255)         cAGATCAATTTGAAATCCCGCCAAACAGGGTTATCGAACTGGGTACTCAGGTCGAAATCTAACGCCAGGACTGTAc                                                                >  1:1382420/1‑76 (MQ=255)           gATCAATTTGAAATCCCGCCAAACAGGGTTATCGAACTGGGTACTCAGGTCGAAATCTAACGCCAGGACTGTACTg                                                              <  1:1060008/76‑1 (MQ=255)            aTCAATTTGAAATCCCGCCAAACAGGGTTATCGAACTGGGTACTCAGGTCGAAATCTAACGCCAGGACTGTACTGc                                                             <  2:125306/76‑1 (MQ=255)                   tGAAATCCCGCCAAACAGGGTTATCGAATTGGGTACTCAGGTCGAAATCTAACGCCAGGACTGTACTGCACCCAtt                                                      >  1:107002/1‑76 (MQ=255)                    gAAATCCCGCCAAACAGGGTTATCGAACTGGGTACTCAGGTCGAAATCTAACGCCAGGACTGTACTGCACCCAttt                                                     <  1:1463130/76‑1 (MQ=255)                      aaTCCCGCCAAACAGGGTTATCGAACTGGGTACTCAGGTCGAAATCTAACGCCAGGACTGTACTGCACCCATTttg                                                   <  1:813455/76‑1 (MQ=255)                      aaTCCCGCCAAACAGGGTTATCGAACTGGGTACTCAGGCCGAAATCTAACGTCAGGACTGTNCTGCACCCATTttg                                                   <  1:846/76‑1 (MQ=255)                       aTCCCGCCAAACAGGGTTATCGAACTGGGTACTCAGGTCGAAATCTAACGCCAGGACTGTACTGCACCCATTttgt                                                  <  1:592025/76‑1 (MQ=255)                         cccGCCAAACAGGGTTATCGAACTGGGTACTCAGGTCGAAATCTAACGCCAGGACTGTACTGCACCCATTttgttg                                                <  2:1382416/76‑1 (MQ=255)                          ccGCCAAACAGGGTTATCGAACTGGGTACTCAGGTCGAAATCTAACGCCAGGACTGTACTGCACCCATTTTGTTgg                                               <  1:335534/76‑1 (MQ=255)                           cGCCAAACAGGGTTATCGAACTGGGTACTCAGTTCGAAATCTAACGCCAGGACTGTACTGCACCCATTTTGTTGGa                                              <  2:1591186/76‑1 (MQ=38)                             ccAAACAGGGTTATCGAACTGGGTACTCAGGTCGAAATCTAACGCCAGGACTGTACTGCACCCATTTTGTTGGACg                                            >  1:159647/1‑76 (MQ=255)                             ccAAACAGGGTTATCGAACTGGGTACTCAGGTCGAAATCTAACGCCAGGACTGTACTGCACCCATTTTGTTGGACg                                            <  1:310900/76‑1 (MQ=255)                             ccAAACAGGGTTATCGAACTGGGTACTCAGGTCGAAATCTAACGCCAGGACTGTACTGCACCCATTTTGTTGGACg                                            <  2:981817/76‑1 (MQ=255)                              cAAACAGGGTTATCGAACTGGGTACTCAGGTCGAAATCTAACGCCAGGACTGTACTGCACCCATTTTGTTGGACGa                                           <  2:1554429/76‑1 (MQ=255)                               aaaCAGGGTTATCGAACTGGGTACTCAGGTCGAAATCTAACGCCAGGACTGTACTGCACCCATGTTGTTGGTCGAt                                          <  1:1650342/76‑1 (MQ=38)                               aaaCAGGCGTATCGAACTGGGTACTCAGGTCGAAATCTAACGCCAGGACTGTACTGCACCCATTTTGTTGGACGAt                                          <  1:1424498/76‑1 (MQ=37)                                aaCAGGGTTATCGAACTGGGTACTCAGGTCGAAATCTAACGCCAGGACTGTACTGCACCCATTTTGTTGGTCGATg                                         <  2:1145699/76‑1 (MQ=38)                                aaCAGGGTTATCGAACTGGGTACTCAGGTCGAAATCTAACGCCAGGACTGTACTGCACCCATTTTGTTGGACGATg                                         >  1:1291676/1‑76 (MQ=255)                                 aCAGGGTTATCGAACTGGGTACTCAGGTCGAAATCTAACGCCAGGACTGTACTGCACCCATTTTGTTGGACGATGa                                        >  1:388095/1‑76 (MQ=255)                                 aCAGGGTTATCGAACTGGGTACTCAGGTCGAAATCTAACGCCAGGACTGTACTGCACCCATTTTGTTGGACGATGa                                        >  1:744250/1‑76 (MQ=255)                                  cAGGGTTATCGAACTGGGTACTCAGGTCGAAATCTAACGCCAGGACTGTACTGCACCCATTTTGTTGGACGATGaa                                       >  2:495799/1‑76 (MQ=255)                                        tATCGAACTGGGTACTCAGGTCGAAATCTAACGCCAGGACTGTACTGCACCCATTTTGTTGGACGATGAAATGGaa                                 >  1:1333583/1‑76 (MQ=255)                                        tATCGAACTGGGTACTCAGGTCGAAATATAAGGCNAGGAGTGNACNGCNCCCACTTTGTTGGACGATGAAATGGaa                                 >  1:1536087/1‑76 (MQ=21)                                         aTCGAACTGGGTACTCAGGTCGAAATCTAACGCCAGGACTGTACTGCACCCATTTTGTTGGACGATGAAATGGAAt                                <  2:128602/76‑1 (MQ=255)                                         aTCGAACTGGGTACTCAGGTCGAAATCTAACGCCAGGACTGTACTGCACCCATTTTGTTGGACGATGAAATGGAAt                                <  2:900660/76‑1 (MQ=255)                                         aTCGAACTGGGTACTCAGGTCGAAATCTAACGCCAGGACTGTACTGCACCCATTTTGTTGGACGATGAAATGGAAt                                >  2:1300133/1‑76 (MQ=255)                                          tCGAACTGGGTACTCAGGTCGAAATCTAACGCCAGGACTGTACTGCACCCATTTTGTTGGACGATGAAATGGAATa                               <  2:579826/76‑1 (MQ=255)                                          tCGAACTGGGTACTCAGGTCGAAATCTAACGCCAGGACTGTACTGCACCCATTTTGCTGGACGATGAAATGGAATa                               <  1:1499929/76‑1 (MQ=25)                                           cGAACTGGGTACTCAGGTCGAAATCTAACGCCAGGACTGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAg                              <  2:1591243/76‑1 (MQ=255)                                           cGAACTGGGTACTCAGGTCGAAATCTAACGCCAGGACTGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAg                              <  1:1025583/76‑1 (MQ=255)                                            gAACTGGGTACTCAGGTCGAAATCTAACGCCAGGACTGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGc                             >  2:902863/1‑76 (MQ=255)                                             aaCTGGGTACTCAGGTCGAATTCTAACGCCAGGACTGTACTGCACCCATTTTGTTGGACGATGAAATGGAACAGcc                            >  2:137580/1‑76 (MQ=21)                                             aaCTGGGTACTCAGGTCGAAATCTAACGCCAGGACTGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGcc                            >  2:1240663/1‑76 (MQ=255)                                              aCTGGGTACTCAGGTCGAAATCTAACGCCAGGACTGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGccc                           <  2:1004450/76‑1 (MQ=255)                                               cTGGGTACTCAGGTCGAAATCTAACGCCAGGACTGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGcccc                          >  1:1284934/1‑76 (MQ=255)                                               cTGGGTACTCAGGTCGAAATCTAACGCCAGGACTGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGcccc                          <  1:765437/76‑1 (MQ=255)                                                tGGGTACTCAGGTCGAAATCTAACGCCAGGACTGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCt                         <  2:1221442/76‑1 (MQ=255)                                                 gggTACTCAGGTCGAAATCTAACGCCAGGACTGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCTa                        >  2:1222293/1‑76 (MQ=255)                                                 gggTACTCAGGCCGAAATCTAACGCCAGGACTGTACTTCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCTa                        >  1:52252/1‑76 (MQ=18)                                                 gggTAATCAGTTCGAAATCTAACGCCAGGACTGTACTGCACCCATTTTGTTGGACGATGAAATCGAATAGCCCCTa                        <  2:63997/76‑1 (MQ=18)                                                   gTACTCAGGTCGAAATCTAACGCCAGGACTGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCTAat                      <  2:1661169/76‑1 (MQ=255)                                                    tACTCAGGTCGGAATCTAACGCCAGGACTGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCTAata                     <  2:934444/76‑1 (MQ=18)                                                    tACTCAGGTCGAAATCTAACGCCAGGACTGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCTAata                     <  2:1197440/76‑1 (MQ=255)                                                    tACTCAGGTCGAAATCTAACGCCAGGACTGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCTAata                     <  2:1656020/76‑1 (MQ=255)                                                     aCTCAGGTCGAAATCTAACGCCAGGACTGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCTAatat                    <  1:856307/76‑1 (MQ=255)                                                     aCTCAGGTCGAAATCTAACGCCAGGACTGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCTAatat                    <  2:1382501/76‑1 (MQ=255)                                                     aCTCAGGTCGAAATCTAACGCCAGGACTGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCTAatat                    <  1:1055039/76‑1 (MQ=255)                                                     aCTCAGGTCGAAATCCAACGCCAGGACTGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCTAatat                    >  1:160774/1‑76 (MQ=18)                                                        cAGGTCGAAATCTAACGCCAGGACTGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCTNATAGGTc                 <  1:137253/76‑1 (MQ=18)                                                        cAGGTCGAAATCTAACGCCAGGACTGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCTAATATGTc                 >  2:1312852/1‑76 (MQ=255)                                                           gTCGAAATCTAATGACAGGACAGTACTGCCCCCATTTTGTCGGGCGATGAAATGGAATAGCCCCTAAAATGTCaaa              <  2:31230/76‑1 (MQ=14)                                                            tCGAAATCTAACGCCAGGACTGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCTAATATGTCAAAg             <  1:1407987/76‑1 (MQ=255)                                                            tCGAAATCTAACGCCAGGACTGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCTAATATGTCAAAg             <  2:1504924/76‑1 (MQ=255)                                                               aaaTCTGACGNCAGGACTGTANTGCACCCATTTTGTTGGACGATGAACTGGAATAGCCCCTAATATGTCAAAGCCa          >  2:1743880/1‑76 (MQ=14)                                                                aaTCTAACGCCAGGACTGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCTAATATGTCAAAGCCaa         <  2:1251725/76‑1 (MQ=255)                                                                aaTCTAACGCCAGGACTGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCTAATATGTCAAAGCCaa         <  2:324295/76‑1 (MQ=255)                                                                aaTCTAACGCCAGGACTGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCTAATATGTCAAAGCCaa         <  2:1489889/76‑1 (MQ=255)                                                                 aTCTAACGCCAGGACTGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCTAATATGTCAAAGCCaaa        <  1:1433652/76‑1 (MQ=255)                                                                 aTCTAACGCCAGGACCGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCTAATATGTCAAAGCCaaa        >  2:791186/1‑76 (MQ=14)                                                                   cTAACGCCAGGACTGTACTGCACCCATTTTGTTGGGCGATGAAATGGAATAGCCCCTAATATGTCAAAGCCAAAAt      >  2:1213329/1‑76 (MQ=14)                                                                     aaCGCCAGGACTGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCTAATATGTCAAAGCCAAAATAc    >  2:1534034/1‑76 (MQ=255)                                                                       cGCCAGGACTGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCTAATATGTCAAAGCCAAAATAccc  <  1:314458/76‑1 (MQ=255)                                                                       cGCCAGGACTGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCTAATATGTCAAAGCCAAAATAccc  >  1:464144/1‑76 (MQ=255)                                                                            |                                                                       GCGTGCACCAGATCAATTTGAAATCCCGCCAAACAGGGTTATCGAACTGGGTACTCAGGTCGAAATCTAACGCCA‑GACTGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCTAATATGTCAAAGCCAAAATACCC  >  REL606/3893476‑3893620

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