Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	minE	1,314,637	A→G	I298V (ATC→GTC)	hisB →	fused histidinol‑phosphatase and imidazoleglycerol‑phosphate dehydratase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	1,314,637	0	A	G	100.0%	313.0 / NA	84	I298V (ATC→GTC)	hisB	fused histidinol‑phosphatase and imidazoleglycerol‑phosphate dehydratase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (43/41);  total (43/41)

GTCGCCCGCACCTGGAATATAAAGCCGAGTTTACCTACCAGCGCGTGGGCGATCTCAGCACCGAAATGATCGAGCACTTCTTCCGTTCGCTCTCATACACCATGGGCGTGACGCTACACCTGAAAACCAAAG  >  minE/1314569‑1314700                                                                     |                                                                gTCGCCCGCACCTGGAATATAAAGCCGAGTTTACCTACCAGCGCGTGGGCGATCTCAGCACCGAAATGGTc                                                               <  1:3955607/71‑1 (MQ=255)      ccGAACCTGGAATATAAAGCCGAGTTTACCTACCAGCGCGTGGGCGATCTCAGCACCGAAATGGTCGAGCa                                                          <  1:2950627/71‑1 (MQ=255)        gCACCTGGAATATAAAGCCGAGTTTACCTACCAGCGCGTGGGCGATCTCAGCACCGAAATGGTCGAGCAct                                                        <  1:3571369/71‑1 (MQ=255)         cACCTGGAATATAAAGCCGAGTTTACCTACCAGCGCGTGGGCGATCTCAGCACCGAAATGGTCGAGCActt                                                       >  1:380408/1‑71 (MQ=255)           ccTGGAATATAAAGCCGAGTTTACCTACCAGCGCGTGGGCGATCTCAGCACCGAAATGGTCGAGCActtct                                                     <  1:417306/71‑1 (MQ=255)              ggAATATAAAGCCGAGTTTACCTACCAGCGCGTGGGCGATCTCAGCACCGAAATGGTCGAGCACTTCTTcc                                                  >  1:223860/1‑71 (MQ=255)               gAATATAAAGCCGAGTTTACCTACCAGCGCGTGGGCGATCTCAGCACCGAAATGGTCg                                                              <  1:847717/58‑1 (MQ=255)               gAATATAAAGCCGAGTTTACCTACCAGCGCGTGGGCGATCTCAGCACCGAAATGGTCGAGCACttcttc                                                   <  1:2599144/69‑1 (MQ=255)               gAATATAAAGCCGAGTTTACCTACCAGCGCGTGGGCGATCTCAGCACCGAAATGGTCGAGCACttcttc                                                   <  1:2304708/69‑1 (MQ=255)               gAATATAAAGCCGAGTTTACCTACCAGCGCGTGGGCGATCTCAGCACCGAAATGGTCGAGCACTTCTTCCg                                                 >  1:955556/1‑71 (MQ=255)                aaTATAAAGCCGAGTTTACCTACCAGCGCGTGGGCGATCTCAGCACCGAAATGGTCGAGCACTTCTTCCGt                                                <  1:1254998/71‑1 (MQ=255)                aaTATAAAGCCGAGTTTACCTACCAGCGCGTGGGCGATCTCAGCACCGAAATGGTCGAGCACTTCTTCCGt                                                <  1:2040106/71‑1 (MQ=255)                 atatAAAGCCGAGTTTACCTACCAGCGCGTGGGCGATCTCAGCACCGAAATGGTCGAGCACTTCTTCCGtt                                               >  1:2740167/1‑71 (MQ=255)                 atatAAAGCCGAGTTTACCTACCAGCGCGTGGGCGATCTCAGCACCGAAATGGTCGAGCACTTCTTCCGtt                                               <  1:2529155/71‑1 (MQ=255)                 atatAAAGCCGAGTTTACCTACCAGCGCGTGGGCGATCTCAGCACCGAAATGGTCGAGCACTTCTTCCGtt                                               <  1:840998/71‑1 (MQ=255)                 atatAAAGCCGAGTTTACCTACCAGCGCGTGGGCGATCTCAGCACCGAAATGGTCGAGCACTTCTTCCGtt                                               <  1:1154314/71‑1 (MQ=255)                   ataAAGCCGAGTTTACCTACCAGCGCGTGGGCGATCTCAGCACCGAAATGGTCGAGCACTTCTTCCGTTc                                              >  1:1163471/1‑70 (MQ=255)                   ataAAGCCGAGTTTACCTACCAGCGCGTGGGCGATCTCAGCACCGAAATGGTCGAGCACTTCTTCCGTTCg                                             >  1:872673/1‑71 (MQ=255)                    taAATCCGAGTTTACCTACCAGCGCGTGGGCGATCTCAGCACCGAAATGGTCGAGc                                                           <  1:1523828/56‑1 (MQ=255)                    taAAGCCGAGTTTACCTACCAGCGCGTGGGCGATCTCAGCACCGAAATGGTCGAGc                                                           <  1:563481/56‑1 (MQ=255)                     aaaGCCGAGTTTACCTACCAGCGCGTGGGCGATCTCAGCACCGAAATGGTCGAGCACTTCTTCCGtt                                               >  1:1519724/1‑67 (MQ=255)                      aaGCCGAGTTTACCTACCAGCGCGTGGGCGATCTCAGCACCGAAATGGTCGAGCACTTCTTcc                                                  <  1:3065811/63‑1 (MQ=255)                      aaGCCGAGTTTACCTACCAGCGCGTGGGCGATCTCAGCACCGAAATGGTCGAGCACTTCTTCCGTTCGctc                                          <  1:3793177/71‑1 (MQ=255)                      aaGCCGAGTTTACCTACCAGCGCGTGGGCGATCTCAGCACCGAAATGGTCGAGCACTTCTTCCGTTCGctc                                          <  1:2535519/71‑1 (MQ=255)                       aGCCGAGTTTACCTACCAGCGCGTGGGCGATCTCAGCACCGAAATGGTCGAGCACTTCTTCCGTTCGctct                                         >  1:1887178/1‑71 (MQ=255)                        gCCGAGTTTACCTACCAGCGCGTGGGCGATCTCAGCACCGAAATGGTCGAGCACTTCTTCCGTTCGctc                                          >  1:866166/1‑69 (MQ=255)                         ccGAGTTTACCTACCAGCGCGTGGGCGATCTCAGCACCGAAATGGTCGAGCACTTCTTCCGTTCGCTCTCa                                       >  1:1081001/1‑71 (MQ=255)                         ccGAGTTTACCTACCAGAGCGTGGGCGATCTCAGCACCGAAATGGTCGAGCACTTCTTCCGTTCGCTCTCa                                       >  1:3152290/1‑71 (MQ=255)                           gAGTTTACCTACCAGCGCGTGGGCGATCTCAGCACCGAAATGGTCGAGCACTTCTTCCGTTCGCTCTCATa                                     >  1:368484/1‑71 (MQ=255)                           gAGTTTACCTACCAGCGCGTGGGCGATCTCAGCACCGAAATGGTCGAGCACTTCTTCCGTTCGCTCTCATa                                     <  1:761202/71‑1 (MQ=255)                             gTTTACCTACCAGCGCGTGGGCGATCTCAGCACCGAAATGGTCGAGCACTTCTTCCGTTCGCTCTCATaca                                   >  1:2054571/1‑71 (MQ=255)                             gTTTACCTACCAGCGCGTGGGCGATCTCAGCACCGAAATGGTCGAGCACTTCTTCCGTTCGCTCTCATaca                                   >  1:3122416/1‑71 (MQ=255)                             gTTTACCTACCAGCGCGTGGGCGATCTCAGCACCGAAATGGTCGAGCACTTCTTCCGTTCGCTCTCATaca                                   >  1:989875/1‑71 (MQ=255)                               ttACCTACCAGCGCGTGGGCGATCTCAGCACCGAAATGGTCGAGc                                                           >  1:1809283/1‑45 (MQ=255)                                tacctaccAGCGCGTGGGCGATCTCAGCACCGAAATGGTCGAGCACTTCTTCCGTTCGCTCTCa                                       >  1:414156/1‑64 (MQ=255)                                tacctaccAGCGCGTGGGCGATCTCAGCACCGAAATGGTCGAGCACTTCTTCCGTTCGCTCTCATACACCa                                >  1:1988263/1‑71 (MQ=255)                                tacctaccAGCGCGTGGGCGATCTCAGCACCGAAATGGTCGAGCACTTCTTCCGTTCGCTCTCATACACCa                                >  1:2546840/1‑71 (MQ=255)                                  cctaccAGCGCGTGGGCGATCTCAGCACCGAAATGGTCGAGCACTTCTTCCGTTCGCTCTCATACACCATg                              >  1:2212450/1‑71 (MQ=255)                                   ctaccAGCGCGTGGGCGATCTCAGCACCGAAATGGTCGAGCACTTCTTCCGTTCGCTCTCATaca                                   <  1:3789317/65‑1 (MQ=255)                                    taccAGCGCGTGGGCGATCTCAGCACCGAAATGGTCGAGCACTTCATCCGTTCg                                             >  1:569949/1‑54 (MQ=39)                                     accAGCGCGTGGGCGATCTCAGCACCGAAATGGTCGAGCACTTCTTCCGTTCGCTCTCATACACCATGGGc                           >  1:319267/1‑71 (MQ=255)                                     accAGCGCGTGGGCGATCTCAGCACCGAAATGGTCGAGCACTTCTTCCGTTCGCTCTCATACACCATGGGc                           >  1:1237998/1‑71 (MQ=255)                                     accAGCGCGTGGGCGATCTCAGCACCGAAATGGTCGAGCACTTCTTCCGTTCGCTCTCATACACCATGGGc                           >  1:4109506/1‑71 (MQ=255)                                     accAGCGCGTGGGCGATCTCAGCACCGAAATGGTCGAGCACTTCTTCCGTTCGCTCTCATACACCATGGGc                           >  1:3532436/1‑71 (MQ=255)                                      ccAGCGCGTGGGCGATCTCAGCACCGAAATGGTCGAGCACTTCTTCCGTTCGctct                                         >  1:3239844/1‑56 (MQ=255)                                      ccAGCGCGTGGGCGATCTCAGCACCGAAATGGTCGAGCACTTCTTCCGTTCGCTCTCATACACCATGGGc                           <  1:1568361/70‑1 (MQ=255)                                       cAGCGCGTGGGCGATCTCAGCACCGAAATGGTCGAGCACTTCTTCCGTTCGCTCTCATACACCATGGGCGt                         <  1:1006655/71‑1 (MQ=255)                                       cAGCGCGTGGGCGATCTCAGCACCGAAATGGTCGAGCACTTCTTCCGTTCGCTCTCATACACCATGGGCGt                         <  1:519130/71‑1 (MQ=255)                                        aGCGCGTGGGCGATCTCAGCACCGAAATGGTCGAGCACTTCTTCCGTTCGCTCTCATACACCATGGGCGt                         >  1:3370652/1‑70 (MQ=255)                                        aGCGCGTGGGCGATCTCAGCACCGAAATGGTCGAGCACTTCTTCCGTTCGCTCTCATACACCATGGGCGTg                        <  1:2726949/71‑1 (MQ=255)                                         gcgcGTGGGCGATCTCAGCACCGAAATGGTCGAGCACTTCTTCCGt                                                <  1:2698460/46‑1 (MQ=255)                                           gcgTGGGCGATCTCAGCACCGAAATGGTCGAGCACTTCTTCCGTTCGCTCTCATACACCATGGGCGTGACg                     >  1:581528/1‑71 (MQ=255)                                           gcgTGGGCGATCTCAGCACCGAAATGGTCGAGCACTTCTTCCGTTCGCTCTCATACACCATGGGCGTGACg                     >  1:714002/1‑71 (MQ=255)                                           gcgTGGGCGATCTCAGCACCGAAATGGTCGAGCACTTCTTCCGTTCGCTCTCATACACCATGGGCGTGACg                     >  1:1028214/1‑71 (MQ=255)                                             gTGGGCGATCTCAGCACCGAAATGGTCGAGCACTTCTTCCGTTCGCTCTCATACACCATGGGCGTGACGc                    <  1:1823142/70‑1 (MQ=255)                                              tGGGCGATCTCAGCACCGAAATGGTCGAGCACTTCTTCCGTTCGCTCTCATACACCATGGGCGTGACGCt                   >  1:3479822/1‑70 (MQ=255)                                                  cGATCTCAGCACCGAAATGGTCGAGCACTTCTTCCGTTCGCTCTCATACACCATGGGCGTGACGCTACAcc              <  1:2960270/71‑1 (MQ=255)                                                   gATCTCAGCACCGAAATGGTCGAGCACTTCTTCCGTTTGCtctc                                        >  1:3225240/1‑44 (MQ=38)                                                   gATCTCAGCACCGAAATGGTCGAGCACTTCTTCCGTTCGCtctc                                        >  1:1607953/1‑44 (MQ=255)                                                   gATCTCAGCACCGAAATGGTCGAGCACTTCTTCCGTTCGCTCTCATACACCATGGGCGTGACGCTACACCt             >  1:3420985/1‑71 (MQ=255)                                                    aTCTCAGCACCGAAATGGTCGAGCACTTCTTCCGTTCGct                                           <  1:2073165/40‑1 (MQ=255)                                                    aTCTCAGCACCGAAATGGTCGAGCACTTCTTCCGTTCGCTCTCATACACCATGGGCGTGACGCTac                 <  1:140816/66‑1 (MQ=255)                                                        cAGCACCGAAATGGTCGAGCACTTCTTCCGTTCGCTCTCATACACCATGGGCGTGACGCTACACCTGaaaa        >  1:3325364/1‑71 (MQ=255)                                                        cAGCACCGAAATGGTCGAGCACTTCTTCCGTTCGCTCTCATACACCATGGGCGTGACGCTACACCTGaaaa        >  1:1464577/1‑71 (MQ=255)                                                        cAGCACCGAAATGGTCGAGCACTTCTTCCGTTCGCTCTCATACACCATGGGCGTGACGCTACACCTGaaaa        >  1:2210318/1‑71 (MQ=255)                                                         aGCACCGAAATGGTCGAGCACTTCTTCCGTTCGCTCTCATACACCATGGGCGTGACGCTACACCTGaa          <  1:1295823/68‑1 (MQ=255)                                                         aGCACCGAAATGGTCGAGCACTTCTTCCGTTCGCTCTCATACACCATGGGCGTGACGCTACACCTGa           <  1:1357701/67‑1 (MQ=255)                                                         aGCACCGAAATGGTCGAGCACTTCTTCCGTTCGCTCTCATACACCATGGGCGTGACGCTACACCTGa           <  1:1792023/67‑1 (MQ=255)                                                         aGCACCGAAATGGTCGAGCACTTCTTCCGTTCGCTCTCATACACCATGGGCGTGACGCTACACCTGAAAAc       >  1:161462/1‑71 (MQ=255)                                                         aGCACCGAAATGGTCGAGCACTTCTTCCGTTCGCTCTCATACACCATGGGCGTGACGCTACACCTGAAAAc       >  1:2874393/1‑71 (MQ=255)                                                          gCACCGAAATGGTCGAGCACTTCTTCCGTTCGCTCTCATACACCATGGGCGt                         <  1:989518/52‑1 (MQ=255)                                                          gCACCGAAATGGTCGAGCACTTCTTCCGTTCGCTCTCATACACCATGGGCGTGACGCTACACCt             <  1:372597/64‑1 (MQ=255)                                                           cACCGAAATGGTCGAGCACTTCTTCCGTTCGCTCTCATACACCATGGGCGTGACGCTACACCTGAAAACCa     <  1:3173532/71‑1 (MQ=255)                                                            aCCGAAATGGTCGAGCACTTCTTCCGTTCGCTCTCATACACCATg                              <  1:1875064/45‑1 (MQ=255)                                                             ccGAAATGGTCGAGCACTTCTTCCGTTCGCTCTCATACACCATGGGCGTGACGCTACACCTGaa          <  1:3014062/64‑1 (MQ=255)                                                             ccGAAATGGTCGAGCACTTCTTCCGTTCGCTCTCATACACCATGGGCGTGACGCTACACCTGaa          <  1:394915/64‑1 (MQ=255)                                                             ccGAAATGGTCGAGCACTTCTTCCGTTCGCTCTCATACACCATGGGCGTGACGCTACACCTGAAAACCaaa   <  1:1993106/71‑1 (MQ=255)                                                             ccGAAATGGTCGAGCACTTCTTCCGTTCGCTCTCATACACCATGGGCGTGACGCTACACCTGAAAACCaaa   <  1:2870837/71‑1 (MQ=255)                                                             ccGAAATGGTCGAGCACTTCTTCCGTTCGCTCTCATACACCATGGGCGTGACGCTACACCTGAAAACCaa    <  1:2502780/70‑1 (MQ=255)                                                             ccGAAATGGTCGAGCACTTCTTCCGTTCGCTCTCATACACCATGGGCGTGACGCTACACCTGAAAACCa     >  1:868059/1‑69 (MQ=255)                                                             ccGAAATGGTCGAGCACTTCTTCCGTTCGCTCTCATACACCATGGGCGTGACGCTACACCTGAAAACCa     >  1:1396197/1‑69 (MQ=255)                                                             ccGAAATGGTCGAGCACTTCTTCCGTTCGATCTCATACACCATGGGCGTGACGCTACACCTGAAAACCaaa   <  1:3825410/71‑1 (MQ=255)                                                              cGAAATGGTCGAGCACTTCTTCCGTTCGCTCTCATACACCATGGGCGTGACGCTACACCTGAAAACCAAAg  >  1:2942698/1‑71 (MQ=255)                                                              cGAAATGGTCGAGCACTTCTTCCGTTCGCTCTCATACACCATGGGCGTGACGCTACACCTGAAAACCAAAg  >  1:219451/1‑71 (MQ=255)                                                                     |                                                                GTCGCCCGCACCTGGAATATAAAGCCGAGTTTACCTACCAGCGCGTGGGCGATCTCAGCACCGAAATGATCGAGCACTTCTTCCGTTCGCTCTCATACACCATGGGCGTGACGCTACACCTGAAAACCAAAG  >  minE/1314569‑1314700

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