Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	minE	558,552	G→A	intergenic (+102/‑202)	iaaA → / → cmr	L‑asparaginase/multidrug efflux system protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	558,552	0	G	A	100.0%	314.1 / NA	96	intergenic (+102/‑202)	iaaA/cmr	L‑asparaginase/multidrug efflux system protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (40/56);  total (40/56)

AAGTGATTAAAAGGGAGCCAATACAGGCAAGTCGTTGAGAATAAAGTGCAGGTTAAACTGGGTAAAGCGGCATCGTCTTATTTCCCTCAAGCGGCCTGTTTACGGTGGGTGATTGTAACGGGCATAGGTTAAA  >  minE/558484‑558616                                                                     |                                                                 aaGTGATTAAAAGGGAGCCAATACAGGCAAGTCGTTGAGAATAAAGTGCAGGTTAAACTGGGTAAagcagc                                                                <  1:388560/71‑1 (MQ=255) aaGTGATTAAAAGGGAGCCAATACAGGCAAGTCGTTGAGAATAAAGTGCAGGTTAAACTGGGTAAagcagc                                                                <  1:3788149/71‑1 (MQ=255) aaGTGATTAAAAGGGAGCCAATACAGGCAAGTCGTTGAGAATAAAGTGCAGGTTAAACTGGGTAAagcagc                                                                <  1:3410910/71‑1 (MQ=255)      aTTAAAAGGGAGCCAATACAGGCAAGTCGTTGAGAATAAAGTGCAGGTTAAACTGGGTAAAGCAGCATCGt                                                           >  1:4783657/1‑71 (MQ=255)         aaaaGGGAGCCAATACAGGCAAGTCGTTGAGAATAAAGTGCAGGTTAAACTGGGTAAAGCAGCATCGTCt                                                         >  1:4561682/1‑70 (MQ=255)         aaaaGGGAGCCAATACAGGCAAGTCGTTGAGAATAAAGTGCAGGTTAAACTGGGTAAAGCAGCATCGTCt                                                         >  1:2413833/1‑70 (MQ=255)         aaaaGGGAGCCAATACAGGCAAGTCGTTGAGAATAAAGTGCAGGTTAAACTGGGTAAAGCAGCATCGTCt                                                         >  1:2664801/1‑70 (MQ=255)         aaaaGGGAGCCAATACAGGCAAGTCGTTGAGAATAAAGTGCAGGTTAAACTGGGTAAAGCAGCATCGTCt                                                         >  1:3712689/1‑70 (MQ=255)         aaaaGGGAGCCAATACAGGCAAGTCGTTGAGAATAAAGTGCAGGTTAAACTGGGTAAAGCAGCATCGTCt                                                         >  1:1778380/1‑70 (MQ=255)         aaaaGGGAGCCAATACAGGCAAGTCGTTGAGAATAAAGTGCAGGTTAAACTGGGTAAAGCAGCATCGTCt                                                         >  1:3610227/1‑70 (MQ=255)         aaaaGGGAGCCAATACAGGCAAGTCGTTGAGAATAAAGTGCAGGTTAAACTGGGTAAAGCAGCATCGTCt                                                         >  1:3002389/1‑70 (MQ=255)          aaaGGGAGCCAATACAGGCAAGTCGTTGAGAATAAAGTGCAGGTTAAACTGGGTAAAGCAGCATCGTCtt                                                        <  1:3989949/70‑1 (MQ=255)          aaaGGGAGCCAATACAGGCAAGTCGTTGAGAATAAAGTGCAGGTTAAACTGGGTAAAGCAGCATCGTCtt                                                        <  1:4404425/70‑1 (MQ=255)          aaaGGGAGCCAATACAGGCAAGTCGTTGAGAATAAAGTGCAGGTTAAACTGGGTAAAGCAGCATCGTCtt                                                        <  1:2870759/70‑1 (MQ=255)          aaaGGGAGCCAATACAGGCAAGTCGTTGAGAATAAAGTGCAGGTTAAACTGGGTAAAGCAGCATCGTCtt                                                        <  1:2906442/70‑1 (MQ=255)          aaaGGGAGCCAATACAGGCAAGTCGTTGAGAATAAAGTGCAGGTTAAACTGGGTAAAGCAGCATCGTCtt                                                        <  1:1776294/70‑1 (MQ=255)                  ccAATACAGGCAAGTCGTTGAGAATAAAGTGCAGGTTAAACTGGGTAAAGCAGCATCGTCTTATTTCCCTc                                               <  1:3779433/71‑1 (MQ=255)                  ccAATACAGGCAAGTCGTTGAGAATAAAGTGCAGGTTAAACTGGGTAAAGCAGCATCGTCTTATTTCCCTc                                               <  1:4444373/71‑1 (MQ=255)                   cAATACAGGCAAGTCGTTGAGAATAAAGTGCAGGTTAAACTGGGTAAAGCAGCATCGTCTTATTTCCCTCa                                              <  1:3024652/71‑1 (MQ=255)                   cAATACAGGCAAGTCGTTGAGAATAAAGTGCAGGTTAAACTGGGTAAAGCAGCATCGTCTTATTTCCCTCa                                              <  1:3349627/71‑1 (MQ=255)                   cAATACAGGCAAGTCGTTGAGAATAAAGTGCAGGTTAAACTGGGTAAAGCAGCATCGTCTTATTTCCCTCa                                              <  1:2810605/71‑1 (MQ=255)                   cAATACAGGCAAGTCGTTGAGAATAAAGTGCAGGTTAAACTGGGTAAAGCAGCATCGTCTTATTTCCCTCa                                              <  1:4345456/71‑1 (MQ=255)                   cAATACAGGCAAGTCGTTGAGAATAAAGTGCAGGTTAAACTGGGTAAAGCAGCATCGTCTTATTTCCCTCa                                              <  1:2672092/71‑1 (MQ=255)                   cAATACAGGCAAGTCGTTGAGAATAAAGTGCAGGTTAAACTGGGTAAAGCAGCATCGTCTTATTTCCCTCa                                              <  1:3714138/71‑1 (MQ=255)                   cAATACAGGCAAGTCGTTGAGAATAAAGTGCAGGTTAAACTGGGTAAAGCAGCATCGTCTTATTTCCCTCa                                              <  1:3738944/71‑1 (MQ=255)                   cAATACAGGCAAGTCGTTGAGAATAAAGTGCAGGTTAAACTGGGTAAAGCAGCATCGTCTTATTTCCCTCa                                              <  1:4102058/71‑1 (MQ=255)                   cAATACAGGCAAGTCGTTGAGAATAAAGTGCAGGTTAAACTGGGTAAAGCAGCATCGTCTTATTTCCCTCa                                              <  1:4351145/71‑1 (MQ=255)                   cAATACAGGCAAGTCGTTGAGAATAAAGTGCAGGTTAAACTGGGTAAAGCAGCATCGTCTTATTTCCCTCa                                              <  1:2215550/71‑1 (MQ=255)                   cAATACAGGCAAGTCGTTGAGAATAAAGTGCAGGTTAAACTGGGTAAAGCAGCATCGTCTTATTTCCCTCa                                              <  1:4767182/71‑1 (MQ=255)                   cAATACAGGCAAGTCGTTGAGAATAAAGTGCAGGTTAAACTGGGTAAAGCAGCATCGTCTTATTTCCCTCa                                              <  1:1360694/71‑1 (MQ=255)                   cAATACAGGCAAGTCGTTGAGAATAAAGTGCAGGTTAAACTGGGTAAAGCAGCATCGTCTTATTTCCCTCa                                              <  1:1822164/71‑1 (MQ=255)                   cAATACAGGCAAGTCGTTGAGAATAAAGTGCAGGTTAAACTGGGTAAAGCAGCATCGTCTTATTTCCCTCa                                              <  1:1967860/71‑1 (MQ=255)                    aaTACAGGCAAGTCGTTGAGAATAAAGTGCAGGTTAAACTGGGTAAAGCAGCATCGTCTTATTTCCCTCa                                              <  1:838164/70‑1 (MQ=255)                    aaTACAGGCAAGTCGTTGAGAATAAAGTGCAGGGTAAACTGGGTAAAGCAGCATCGTCtt                                                        <  1:2884192/60‑1 (MQ=255)                          ggCAAGTCGTTGAGAATAAAGTGCAGGTTAAACT‑GGTAAAGCAGCATCGTCTTATTTCCCTCAAGCGGCCt                                      <  1:653067/71‑1 (MQ=255)                           gCAAGTCGTTGAGAATAAAGTGCAGGTTAAACTGGGTAAAGCAGCATCGTCTTATTTCCCTCAAGCGGCCt                                      <  1:1798082/71‑1 (MQ=255)                           gCAAGTCGTTGAGAATAAAGTGCAGGTTAAACTGGGTAAAGCAGCATCGTCTTATTTCCCTCAAGCGGCCt                                      <  1:1655310/71‑1 (MQ=255)                           gCAAGTCGTTGAGAATAAAGTGCAGGTTAAACTGGGTAAAGCAGCATCGTCTTATTTCCCTCAAGCGGCCt                                      <  1:3826872/71‑1 (MQ=255)                           gCAAGTCGTTGAGAATAAAGTGCAGGTTAAACTGGGTAAAGCAGCATCGTCTTATTTCCCTCAAGCGGCCt                                      <  1:4812228/71‑1 (MQ=255)                           gCAAGTCGTTGAGAATAAAGTGCAGGTTAAACTGGGTAAAGCAGCATCGTCTTATTTCCCTCAAGCGGCCt                                      <  1:1941865/71‑1 (MQ=255)                               gTCGTTGAGAATAAAGTGCAGGTTAAACTGGGTAAAGCAGCATCGTCTTATTTCCCTCAAGCgg                                         <  1:2295214/64‑1 (MQ=255)                               gTCGTTGAGAATAAAGTGCAGGTTAAACTGGGTAAAGCAGCATCGTCTTATTTCCCTCAAGCgg                                         <  1:1949797/64‑1 (MQ=255)                               gTCGTTGAGAATAAAGTGCAGGTTAAACTGGGTAAAGCAGCATCGTCTTATTTCCCTCAAGCgg                                         <  1:948081/64‑1 (MQ=255)                               gTCGTTGAGAATAAAGTGCAGGTTAAACTGGGTAAAGCAGCATCGTCTTATTTCCCTCAAGCgg                                         <  1:2607971/64‑1 (MQ=255)                                 cGTTGAGAATAAAGTGCAGGTTAAACTGGGTAAAGCAGCATCGTCTTATTTCCCTCAAGCGGcc                                       >  1:4331649/1‑64 (MQ=255)                                 cGTTGAGAATAAAGTGCAGGTTAAACTGGGTAAAGCAGCATCGTCTTATTTCCCTCAAGCGGcc                                       >  1:2915913/1‑64 (MQ=255)                                 cGTTGAGAATAAAGTGCAGGTTAAACTGGGTAAAGCAGCATCGTCTTATTTCCCTCAAGCGGcc                                       >  1:3753914/1‑64 (MQ=255)                                 cGTTGAGAATAAAGTGCAGGTTAAACTGGGTAAAGCAGCATCGTCTTATTTCCCTCAAGCGGcc                                       >  1:3106658/1‑64 (MQ=255)                                 cGTTGAGAATAAAGTGCAGGTTAAACTGGGTAAAGCAGCATCGTCTTATTTCCCTCAAGCGGcc                                       >  1:3482932/1‑64 (MQ=255)                                 cGTTGAGAATAAAGTGCAGGTTAAACTGGGTAAAGCAGCATCGTCTTATTTCCCTCAAGCGGcc                                       >  1:964496/1‑64 (MQ=255)                                 cGTTGAGAATAAAGTGCAGGTTAAACTGGGTAAAGCAGCATCGTCTTATTTCCCTCAAGCGGCCTGTTTAc                                >  1:3766411/1‑71 (MQ=255)                                 cGTTGAGAATAAAGTGCAGGTTAAACTGGGTAAAGCAGCATCGTCTTATTTCCCTCAAGCGGCCTGTTTAc                                >  1:58188/1‑71 (MQ=255)                                 cGTTGAGAATAAAGTGCAGGTTAAACTGGGTAAAGCAGCATCGTCTTATTTCCCTCAAGCGGCCTGTTTAc                                >  1:2745671/1‑71 (MQ=255)                                 cGTTGAGAATAAAGTGCAGGTTAAACTGGGTAAAGCAGCATCGTCTTATTTCCCTCAAGCGGCCTGTTTAc                                >  1:4793303/1‑71 (MQ=255)                                  gTTGAGAATAAAGTGCAGGTTAAACTGGGTAAAGCAGCATCGTCTTATTTCCCTCAAGCGGCCTGTTTAc                                <  1:376241/70‑1 (MQ=255)                                  gTTGAGAATAAAGTGCAGGTTAAACTGGGTAAAGCAGCATCGTCTTATTTCCCTCAAGCGGCCTGTTTAc                                <  1:678256/70‑1 (MQ=255)                                  gTTGAGAATAAAGTGCAGGTTAAACTGGGTAAAGCAGCATCGTCTTATTTCCCTCAAGCGGCCTGTTTAc                                <  1:2302561/70‑1 (MQ=255)                                  gTTGAGAATAAAGTGCAGGTTAAACTGGGTAAAGCAGCATCGTCTTATTTCCCTCAAGCGGCCTGTTTAc                                <  1:1999093/70‑1 (MQ=255)                                  gTTGAGAATAAAGTGCAGGTTAAACTGGGTAAAGCAGCATCGTCTTATTTCCCTCAAGCGGCCTGTTTAc                                <  1:2011322/70‑1 (MQ=255)                                       gaATAAAGTGCAGGTTAAACTGGGTAAAGCAGCATCGTCTTATTTc                                                   >  1:3365706/1‑46 (MQ=255)                                       gaATAAAGTGCAGGTTAAACTGGGTAAAGCAGCATCGTCTTATTTCCCTCAAGCGGCCTGTTTACggtggg                          >  1:166959/1‑71 (MQ=255)                                       gaATAAAGTGCAGGTTAAACTGGGTAAAGCAGCATCGTCTTATTTCCCTCAAGCGGCCTGTTTACggtggg                          >  1:4322536/1‑71 (MQ=255)                                       gaATAAAGTGCAGGTTAAACTGGGTAAAGCAGCATCGTCTTATTTCCCTCAAGCGGCCTGTTTACggtggg                          >  1:3035794/1‑71 (MQ=255)                                       gaATAAAGTGCAGGTTAAACTGGGTAAAGCAGCATCGTCTTATTTCCCTCAAGCGGCCTGTTTACggtggg                          >  1:294881/1‑71 (MQ=255)                                       gaATAAAGTGCAGGTTAAACTGGGTAAAGCAGCATCGTCTTATTTCCCTCAAGCGGCCTGTTTACggtggg                          >  1:4835390/1‑71 (MQ=255)                                       gaATAAAGTGCAGGTTAAACTGGGTAAAGCAGCATCGTCTTATTTCCCTCAAGCGGCCTGTTTACggtggg                          >  1:247157/1‑71 (MQ=255)                                       gaATAAAGTGCAGGTTAAACTGGGTAAAGCAGCATCGTCTTATTTCCCTCAAGCGGCCTGTTTACggtggg                          >  1:2394126/1‑71 (MQ=255)                                       gaATAAAGTGCAGGTTAAACTGGGTAAAGCAGCATCGTCTTATTTCCCTCAAGCGGCCTGTTTACggtggg                          >  1:1112631/1‑71 (MQ=255)                                       gaATAAAGTGCAGGTTAAACTGGGTAAAGCAGCATCGTCTTATTTCCCTCAAGCGGCCTGTTTACggtgg                           >  1:4318033/1‑70 (MQ=255)                                       gaATAAAGTGCAGGTTAAACTGGGTAAAGCAGCATCGTCTTATTTCCCTCAAGCGGCCTGTTTACggtgg                           >  1:979225/1‑70 (MQ=255)                                         aTAAAGTGCAGGTTAAACTGGGTAAAGCAGCATCGTCTTATTTCCCTCAAGCGGCCTGTTTACggtgggt                         <  1:1129104/70‑1 (MQ=255)                                         aTAAAGTGCAGGTTAAACTGGGTAAAGCAGCATCGTCTTATTTCCCTCAAGCGGCCTGTTTACggtgggt                         <  1:2621773/70‑1 (MQ=255)                                          tAAAGTGCAGGTTAAACTGGGTAAAGCAGCATCGTCTTATTTCCCTCAAGCGGCCTGTTTATggtgggt                         <  1:4192371/69‑1 (MQ=255)                                           aaaGTGCAGGTTAAACTGGGTAAAGCAGCATCGTCTTATTTCCCTCAAGCgg                                         <  1:2009413/52‑1 (MQ=255)                                            aaGTGCAGGTTAAACTGGGTAAAGCAGCATCGTCTTATTTCCCTCAAgag                                          >  1:3593718/1‑48 (MQ=255)                                            aaGTGCAGGTTAAACTGGGTAAAGCAGCATCGTCTTATTTCCCTCAAGCGGCCTGTTTACGGTGGGTGAtt                     >  1:1152658/1‑71 (MQ=255)                                            aaGTGCAGGTTAAACTGGGTAAAGCAGCATCGTCTTATTTCCCTCAAGCGGCCTGTTTACGGTGGGTGAtt                     >  1:4714436/1‑71 (MQ=255)                                            aaGTGCAGGTTAAACTGGGTAAAGAAGCATCGTCTTATTTCCCTCAAGCGGCCTGTTTACGGTGGGTGAtt                     >  1:2100504/1‑71 (MQ=255)                                               tGCAGGTTAAACTGGGTAAAGCAGCATCGTCTTATTTCCCTCAAGCGGCCTGTTTACGGTGGGTGATTGTa                  <  1:4845074/71‑1 (MQ=255)                                               tGCAGGTTAAACTGGGTAAAGCAGCATCGTCTTATTTCCCTCAAGCGGCCTGTTTACGGTGGGTGATTGTa                  <  1:2417079/71‑1 (MQ=255)                                                    gTTAAACTGGGTAAAGCAGCATCGTCTTATTTCCCTCATGCGGCCTGtt                                   <  1:3555559/49‑1 (MQ=255)                                                     ttAAACTGGGTAAAGCAGCATCGTCTTATTTCCCTCAAGCGGCCTGTTTACGGTGGGTGATTGTAACGGGc            <  1:4756636/71‑1 (MQ=255)                                                      tAAACTGGGTAAAGCAGCATCGTCTTATTTCCCTCAAGCGGCCTGTTTACGGTGGGTGATTGTAACgg              >  1:4325348/1‑68 (MQ=255)                                                      tAAACTGGGTAAAGCAGCATCGTCTTATTTCCCTCAAGCGGCCTGTTTACGGTGGGTGATTGTAACgg              >  1:200356/1‑68 (MQ=255)                                                      tAAACTGGGTAAAGCAGCATCGTCTTATTTCCCTCAAGCGGCCTGTTTACGGTGGGTGATTGTAACgg              >  1:215808/1‑68 (MQ=255)                                                      tAAACTGGGTAAAGCAGCATCGTCTTATTTCCCTCAAGCGGCCTGTTTACGGTGGGTGATTGTAACgg              >  1:2687865/1‑68 (MQ=255)                                                      tAAACTGGGTAAAGCAGCATCGTCTTATTTCCCTCAAGCGGCCTGTTTACGGTGGGTGATTGTAACgg              >  1:2809473/1‑68 (MQ=255)                                                       aaaCTGGGTAAAGCAGCATCGTCTTATTTCCCTCAAGCGGCCt                                      <  1:1802913/43‑1 (MQ=255)                                                       aaaCTGGGTAAAGCAGCATCGTCTTATTTCCCTCAAGCGGCCTGTTTACGGTGGGTGATTGTAACGGGCa           <  1:2263652/70‑1 (MQ=255)                                                       aaaCTGGGTAAAGCAGCATCGTCTTATTTCCCTCAAGCGGCCTGTTTACGGTGGGTGATTGTAACGGGCa           <  1:209924/70‑1 (MQ=255)                                                       aaaCTGGGTAAAGCAGCATCGTCTTATTTCCCTCAAGCGGCCTGTTTACGGTGGGTGATTGTAACGGGCa           <  1:4379577/70‑1 (MQ=255)                                                       aaaCTGGGTAAAGCAGCATCGTCTTATTTCCCTCAAGCGGCCTGTTTACGGTGGGTGATTGTAACGGGCa           <  1:619438/70‑1 (MQ=255)                                                       aaaCTGGGTAAAGCAGCATCGTCTTATTTCCCTCAAGCGGCCTGTTTACGGTGGGTGATTGTAACGGGCa           <  1:1640229/70‑1 (MQ=255)                                                       aaaCTGGGTAAAGCAGCATCGTCTTATTTCCCTCAAGCGGCCTGTTTACGGTGGGTGATTGTAACGGGCa           <  1:1621412/70‑1 (MQ=255)                                                       aaaCTGGGTAAAGCAGCATCGTCTTATTTCCCTCAAGCGGCCTGTTTACGGTGGGTGATTGTAACGGGCa           <  1:15479/70‑1 (MQ=255)                                                         aCTGGGTAAAGCAGCATCGTCTTATTTCCCTCAAGCgg                                         <  1:4906404/38‑1 (MQ=255)                                                             ggTAAAGCAGCATCGTCTTATTTCCCTCAAGCGGCCt                                      <  1:2102303/37‑1 (MQ=255)                                                                aaaGCAGCATCGTCTTATTTCCCTCAAGCGGCCTGTTTACGGTGGGTGATTGTAACGGGCATAGGTtaaa  >  1:2409951/1‑70 (MQ=255)                                                                aaaGCAGCATCGTCTTATTTCCCTCAAGCGGCCTGTTTACGGTGGGTGATTGTAACGGGCATAGGTtaaa  >  1:3920748/1‑70 (MQ=255)                                                                     |                                                                 AAGTGATTAAAAGGGAGCCAATACAGGCAAGTCGTTGAGAATAAAGTGCAGGTTAAACTGGGTAAAGCGGCATCGTCTTATTTCCCTCAAGCGGCCTGTTTACGGTGGGTGATTGTAACGGGCATAGGTTAAA  >  minE/558484‑558616

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