Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,736,520	A→G	Q218R (CAG→CGG)	purR →	PurR

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,736,520	0	A	G	100.0%	245.9 / NA	79	Q218R (CAG→CGG)	purR	PurR
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (34/45);  total (34/45)

TGGAACGTAACACCGGCGCAGGCCGCCTTGCCGGTTTTATGAAGGCGATGGAAGAAGCGATGATCAAGGTGCCGGAAAGCTGGATTGTGCAGGGTGACTTTGAACCTGAATCCGGTTATCGCGCCATGCAGCAAATCCTGTCGCAGCCGCATCGCCCTACTGCCGTCTTCTGTGGTGG  >  NC_000913/1736430‑1736607                                                                                           |                                                                                        tGGAACGTAACACCGGCGCAGGCCGCCTTGCCGGTTTTATGAAGGCGATGGAAGAAGCGATGATCAAGGTGCCGGAAAGCTGGATTGTGCggg                                                                                       <  2:1500893/93‑1 (MQ=255)  ggAACGTAACACCGGCGCAGGCCGCCTTGCCGGTTTTATGAAGGCGATGGAAGAAGCGATGATCAAGGTGCCGGAAAGCTGGATTGTGCgggg                                                                                      >  2:4024850/1‑93 (MQ=255)   gAACGTAACACCGGCGCAGGCCGCCTTGCCGGTTTTATGAAGGCGATGGAAGAAGCGATGATCAAGGTGCCGGAAAGCTGGATTGTGCGGGGt                                                                                     <  2:3767205/93‑1 (MQ=255)   gAACGTAACACCGGCGCAGGCCGCCTTGCCGGTTTTATGAAGGCGATGGAAGAAGCGATGATCAAGGTGCCGGAAAGCTGGATTGTGCGGGGt                                                                                     <  1:2007475/93‑1 (MQ=255)    aaCGTAACACCGGCGCAGGCCGCCTTGCCGGTTTTATGAAGGCGATGGAAGAAGCGATGATCAAGGTGCCGGAAAGCTGGATTGTGCGGGGTg                                                                                    >  2:1003837/1‑93 (MQ=255)    aaCGTAACACCGGCGCAGGCCGCCTTGCCGGTTTTATGAAGGCGATGGAAGAAGCGATGATCAAGGTGCCGGAAAGCTGGATTGTGCGGGGTg                                                                                    >  2:1842461/1‑93 (MQ=255)     aCGTAACACCGGCGCAGGCCGCCTTGCCGGTTTTATGAAGGCGATGGAAGAAGCGATGATCAAGGTGCCGGAAAGCTGGATTGTGCGGGGTGa                                                                                   >  1:1590127/1‑93 (MQ=255)     aCGTAACACCGGCGCAGGCCGCCTTGCCGGTTTTATGAAGGCGATGGAAGAAGCGATGATCAAGGTGCCGGAAAGCTGGATTGTGCGGGGTGa                                                                                   >  1:193825/1‑93 (MQ=255)      cGTAACACCGGCGCAGGCCGCCTTGCCGGTTTTATGAAGGCGATGGAAGAAGCGATGATCAAGGTGCCGGAAAGCTGGATTGTGCGGGGTGAc                                                                                  <  1:2783350/93‑1 (MQ=255)        tAACACCGGCGCAGGCCGCCTTGCCGGTTTTATGAAGGCGATGGAAGAAGCGATGATCAAGGTGCCGGAAAGCTGGATTGTGCGGGGTGACtt                                                                                <  1:3699638/93‑1 (MQ=255)         aaCACCGGCGCAGGCCGCCTTGCCGGTTTTATGAAGGCGATGGAAGAAGCGATGATCAAGGTGCCGGAAAGCTGGATTGTGCGGGGTGACttt                                                                               <  2:4306998/93‑1 (MQ=255)          acacCGGCGCAGGCCGCCTTGCCGGTTTTATGAAGGCGATGGAAGAAGCGATGATCAAGGTGCCGGAAAGCTGGATTGTGCGGGGTGACTTTg                                                                              >  1:528314/1‑93 (MQ=255)           cacCGGCGCAGGCCGCCTTGCCGGTTTTATGAAGGCGATGGAAGAAGCGATGATCAAGGTGCCGGAAAGCTGGATTGTGCGGGGTGACTTTGa                                                                             <  1:4410849/93‑1 (MQ=255)              cGGCGCAGGCCGCCTTGCCGGTTTTATGAAGGCGATGGAAGAAGCGATGATCAAGGTGCCGGAAAGCTGGATTGTGCGGGGTGACTTTGAAcc                                                                          <  1:1163016/93‑1 (MQ=255)              cGGCGCAGGCCGCCTTGCCGGTTTTATGAAGGCGATGGAAGAAGCGATGATCAAGGTGCCGGAAAGCTGGATTGTGCGGGGTGACTTTGAAcc                                                                          <  2:2418641/93‑1 (MQ=255)               ggCGCAGGCCGCCTTGCCGGTTTTATGAAGGCGATGGAAGAAGCGATGATCAAGGTGCCGGAAAGCTGGATTGTGCGGGGTGACTTTGAACCt                                                                         >  1:1367481/1‑93 (MQ=255)                 ggcAGGTCGCCTTGCTGGTTTTATGAAGGCGATGGAAGAAGCGATGATCAAGGTGCCGGAAAGCTGGATTGTGCGGGGTGACTTTGAACCTGa                                                                       <  1:4447425/92‑1 (MQ=255)                 cgcAGGCCGCCTTGCCGGTTTTATGAAGGCGATGGAAGAAGCGATGATCAAGGTGCCGGAAAGCTGGATTGTGCGGGGTGACTTTGAACCTGa                                                                       <  2:4162395/93‑1 (MQ=255)                  gcAGTCCTCCTTGCCGGTTTTATGAAGGCGATGGAAGAAGCGATGATCAAGGTGCCGGAAAGCTGGATTGTGCGGGGTGCCTTTGAACCTGaa                                                                      <  2:631020/93‑1 (MQ=255)                  gcAGGCCGCCTTGCCGGTTTTATGAAGGCGATGGAAGAAGCGATGATCAAGGTGCCGGAAAGCTGGATTGTGCGGGGTGACTTTGAACCTGaa                                                                      >  1:522088/1‑93 (MQ=255)                  gcAGGCCGCCTTGCCGGTTTTATGAAGGCGATGGAAGAAGCGATGATCAAGGTGCCGGAAAGCTGGATTGTGCGGGGTGACTTTGAACCTGaa                                                                      <  1:302806/93‑1 (MQ=255)                    cggCCGCCTTGCCGGTTTTATGAAGGCGATGGAAGAAGCGATGATCAAGGTGCCGGAAAGCTGGATTGTGCGGGGTGACTTTGAACCTGAATc                                                                    <  1:3103015/92‑1 (MQ=255)                    aGGCCGCCTTGCCGGTTTTATGAAGGCGATGGAAGAAGCGATGATCAAGGTGCCGGAAAGCTGGATTGTGCGGGGTGACTTTGAACCTGAATc                                                                    <  1:2726736/93‑1 (MQ=255)                      gccgccTTGCCGGTTTTATGAAGGCGATGGAAGAAGCGATGATCAAGGTGCCGGAAAGCTGGATTGTGCGGGGTGACTTTGAACCTGAATCCg                                                                  <  1:297123/93‑1 (MQ=255)                           cTTGCCGGTTTTATGAAGGCGATGGAAGAAGCGATGATCAAGGTGCCGGAAAGCTGGATTGTGCGGGGTGACTTTGAACCTGAATCCGGTTAt                                                             >  2:3736921/1‑93 (MQ=255)                           cTTGCCGGTTTTATGAAGGCGACGGAAGAAGCGATGATCAAGGTGCCGGAAAGCTGTATTGTGCGGGGTGACTTTGAACCTGAATCCGGTTAt                                                             <  1:4547384/93‑1 (MQ=255)                            ttGCCGGTTTTATGAAGGCGATGGAAGAAGCGATGATCAAGGTGCCGGAAAGCTGGATTGTGCGGGGTGACTTTGAACCTGAATCCGGTTATc                                                            >  2:2790916/1‑93 (MQ=255)                            ttGCCGGTTTTATGAAGGCGATGGAAGAAGCGATGATCAAGGTGCCGGAAAGCTGGATTGTGCGGGGTGACTTTGAACCTGAATCCGGTTATc                                                            <  1:2999433/93‑1 (MQ=255)                            ttGCCGGTTTTATGAAGGCGATGGAAGAAGCGATGATCAAGGTGCCGGAAAGCTGGATTGTGCGGGGTGACTTTGAACCTGAATCCGGTTATc                                                            >  2:1119731/1‑93 (MQ=255)                              gCCGGTTTTATGAAGGCGATGGAAGAAGCGATGATCAAGGTGCCGGAAAGCTGGATTGTGCGGGGTGACTTTGAACCTGAATCCGGTTATcgc                                                          <  1:4204892/93‑1 (MQ=255)                                cGGTTTTATGAAGGCTATGGAAGAAGCGATGATCAAGGTGCCGGAAAGCTGGATTGTGCGGGGTGACTTTGAACCTGAATCCGGTTATCgcgc                                                        >  1:3704219/1‑93 (MQ=255)                                cGGTTTTATGAAGGCGATGGAAGAAGCGATGATCAAGGTGCCGGAAAGCTGGATTGTGCGGGGTGACTTTGAACCTGAATCCGGTTATCgcgc                                                        >  1:442299/1‑93 (MQ=255)                                  gTTTTATGAAGGCGATGGAAGAAGCGATGATCAAGGTGCCGGAAAGCTGGATTGTGCGGGGTGACTTTGAACCTGAATCCGGTTATCGCGCCa                                                      >  1:4650617/1‑93 (MQ=255)                                   ttttATGAAGGCGATGGAAGAAGCGATGATCAAGGTGCCGGAAAGCTGGATTGTGCGGGGTGACTTTGAACCTGAATCCGGTTATCGCGCCAt                                                     >  1:204810/1‑93 (MQ=255)                                   ttttATGAAGGCGATGGAAGAAGCGATGATCAAGGTGCCGGAAAGCTGGATTGTGCGGGGTGACTTTGAACCTGAATCCGGTTATCGCGCCAt                                                     >  2:2571558/1‑93 (MQ=255)                                     ttATGAAGGCGATGGAAGAAGCGATGATCAAGGTGCCGGAAAGCTGGATTGTGCGGGGCGACTTTGAACCTGAATCCGGTTATCGCGCCATgc                                                   >  1:3878355/1‑93 (MQ=255)                                       aTGAAGGCGATGGAAGAAGCGATGATCAAGGTGCCGGAAAGCTGGATTGTGCGGGGTGACTTTGAACCTGAATCCGGTTATCGCGCCATgcag                                                 >  1:1492963/1‑93 (MQ=255)                                        tGAAGGCGATGGAAGAAGCGATGATCAAGGTGCCGGAAAGCTGGATTGTGCGGGGTGACTTTGAACCTGAATCCGGTTATCGCGCCATgcagc                                                <  1:193097/93‑1 (MQ=255)                                           aGGCGATGGAAGAAGCGATGATCAAGGTGCCGGAAAGCTGGATTGTGCGGGGTGACTTTGAACCTGAATCCGGTTATCGCGCCATGCAGCaaa                                             <  1:1844456/93‑1 (MQ=255)                                           aGGCGATGGAAGAAGCGATGATCAAGGTGCCGGAAAGCTGGATTGTGCGGGGTGACTTTGAACCTGAATCCGGTTATCGCGCCATGCAGCaaa                                             <  2:1016400/93‑1 (MQ=255)                                            ggCGATGGAAGAAGCGATGATCAAGGTGCCGGAAAGCTGGATTGTGCGGGGTGACTTTGAACCTGAATCCGGTTATCGCGCCATGCAGCAAAt                                            <  2:3235635/93‑1 (MQ=255)                                              cGATGGAAGAAGCGATGATCAAGGTGCCGGAAAGCTGGATTGTGCGGGGTGACTTTGAACCTGAATCCGGTTATCGCGCCATGCAGCAAATcc                                          <  2:2914979/93‑1 (MQ=255)                                               gATGGAAGAAGCGATGATCAAGGTGCCGGAAAGCTGGATTGTGCGGGGTGACTTTGAACCTGAATCCGGTTATCGCGCCATGCAGCAAATCCt                                         >  1:3600753/1‑93 (MQ=255)                                                aTGGAAGAAGCGATGATCAAGGTACCGGAAAGCTGGATTGTGCGGGGTGACTTTGAACCTGAATCCGGTTATCGCGCCATGCAGCAAATCCTg                                        <  2:1429044/93‑1 (MQ=255)                                                 tGGAAGAAGCGATGATCAAGGTGCCGGAAAGCTGGATTGTGCGGGGTGACTTTGAACCTGAATCCGGTTATCGCGCCATGCAGCAAATCCTGt                                       >  2:2984485/1‑93 (MQ=255)                                                 tGGAAGAAGCGATGATCAAGGTGCCGGAAAGCTGGATTGTGCGGGGTGACTTTGAACCTGAATCCGGTTATCGCGCCATGCAGCAAATCCTGt                                       >  1:330853/1‑93 (MQ=255)                                                 tGGAAGAAGCGATGATCAAGGTGCCGGAAAGCTGGATTGTGCGGGGTGACTTTGAACCTGAATCCGGTTATCGCGCCATGCAGCAAATCCTGt                                       >  2:3668167/1‑93 (MQ=255)                                                 tGGAAGAAGCGATGATCAAGGTGCCGGAAAGCTGGATTGTGCGGGGTGACTTTGAACCTGAATCCGGTTATCGCGCCATGCAGCAAATCCTGt                                       <  2:446701/93‑1 (MQ=255)                                                  ggAAGAAGCGATGATCAAGGTGCCGGAAAGCTGGATTGTGCGGGGTGACTTTGAACCTGAATCCGGTTATCGCGCCATGCAGCAAATCCTGTc                                      >  2:1514542/1‑93 (MQ=255)                                                    aagaagCGATGATCAAGGTGCCGGAAAGCTGGATTGTGCGGGGTGACTTTGAACCTGAATCCGGTTATCGCGCCATGCAGCAAATCCTGTCGc                                    <  1:1253960/93‑1 (MQ=255)                                                      gaagCGATGATCAAGGTGCCGGAAAGCTGGATTGTGCGGGGTGACTTTGAACCTGAATCCGGTTATCGCGCCATGCAGCAAATCCTGTCGCAg                                  >  1:893465/1‑93 (MQ=255)                                                      gaagCGATGATCAAGGTGCCGGAAAGCTGGATTGTGCGGGGTGACTTTGAACCTGAATCCGGTTATCGCGCCATGCAGCAAATCCTGTCGCAg                                  >  1:2458209/1‑93 (MQ=255)                                                      gaagCGATGATCAAGGTGCCGGAAAGCTGGATTGTGCGGGGTGACTTTGAACCTGAATCCGGTTATCGCGCCATGCAGCAAATCCTGTCGCAg                                  >  1:2694564/1‑93 (MQ=255)                                                           gatgatCAAGGTGCCGGAAAGCTGGATTGTGCGGGGTGACTTTGAACCTGAATCCGGTTATCGCGCCATGCAGCAAATCCTGTCGCAGCCGCa                             <  2:4457230/93‑1 (MQ=255)                                                           gatgatCAAGGTGCCGGAAAGCTGGATTGTGCGGGGTGACTTTGAACCTGAATCCGGTTATCGCGCCATGCAGCAAATCCTGTCGCAGCCGCa                             <  2:1984322/93‑1 (MQ=255)                                                            atgatCAAGGTGCCGGAAAGCTGGATTGTGCGGGGTGACTTTGAACCTGAATCCGGTTATCGCGCCATGCAGCAAATCCTGTCGCAGCCGCAt                            <  1:3118381/93‑1 (MQ=255)                                                             tgatCAAGGTGCCGGAAAGCTGGATTGTGCGGGGTGACTTTGAACCTGAATCCGGTTATCGCGCCATGCAGCAAATCCTGTCGCAGCCGCATc                           <  1:172310/93‑1 (MQ=255)                                                             tgatCAAGGTGCCGGAAAGCTGGATTGTGCGAGGTGACTTTGAACCTGAATCCGGTTATCGCGCCATGCAGCAAATCCTGTCGCAGCCGCATc                           <  1:2192629/93‑1 (MQ=255)                                                             tgaaCAAGGTGCCGGAAAGCTGGATTGTGCGGGGTGACTTTGAACCTGAATCCGGTTATCGCGCCATGCAGCAAATCCTGTCGCAGCCGCATc                           <  1:2570466/93‑1 (MQ=255)                                                              gatCAAGGTGCCGGAAAGCTGGATTGTGCGGGGTGACTTTGAACCTGAATCCGGTTATCGCGCCATGCAGCAAATCCTGTCGCAGCCGCATCg                          >  1:1005835/1‑93 (MQ=255)                                                              gatCAAGGTGCCGGAAAGCTGGATTGTGCGGGGTGACTTTGAACCTGAATCCGGTTATCGCGCCATGCAGCAAATCCTGTCGCAGCCGCATCg                          <  2:3084689/93‑1 (MQ=255)                                                               atCAAGGTGCCGGAAAGCTGGATTGTGCGGGGTGACTTTGAACCTGAATCCGGTTATCGCGCCATGCAGCAAATCCTGTCGCAGCCGCATCGc                         >  2:2920497/1‑93 (MQ=255)                                                               atCAAGGTGCCGGAAAGCTGGATTGTGCGGGGTGACTTTGAACCTGAATCCGGTTATCGCGCCATGCAGCAAATCCTGTCGCAGCCGCATCGc                         >  1:3478092/1‑93 (MQ=255)                                                                  aaGGTGCCGGAAAGCTGGATTGTGCGGGGTGACTTTGAACCTGAATCCGGTTATCGCGCCATGCAGCAAATCCTGTCGCAGCCGCATCGCCCt                      <  1:3872896/93‑1 (MQ=255)                                                                   aGGTGCCGGAAAGCTGGATTGTGCGGGGTGACTTTGAACCTGAATCCGGTTATCGCGCCATGCAGCAAATCCTGTGGCAGCCGCATCGCCCTa                     <  2:3607886/93‑1 (MQ=255)                                                                   aGGTGCCGGAAAGCTGGATTGTGCGGGGTGACTTTGAACCTGAATCCGGTTATCGCGCCATGCAGCAAATCCTGTCGCAGCCGCATCGCCCTa                     >  1:3896799/1‑93 (MQ=255)                                                                    ggTGCCGGAAAGCTGGATTGTGCGGGGTGACTTTGAACCTGAATCCGGTTATCGCGCCATGCAGCAAATCCTGTCGCAGCCGCATCGCCCTAc                    <  1:2913318/93‑1 (MQ=255)                                                                    ggTGCCGGAAAGCTGGATTGTGCGGGGTGACTTTGAACCTGAATCCGGTTATCGCGCCATGCAGCAAATCCTGTCGCAGCCGCATCGCCCTAc                    <  1:2375239/93‑1 (MQ=255)                                                                     gTGCCGGAAAGCTGGATTGTGCGGGGTGACTTTGAACCTGAATCCGGTTATCGCGCCATGCAGCAAATCCTGTCGCAGCCGCATCGCCCTACt                   <  1:3443462/93‑1 (MQ=255)                                                                        ccGGAAAGCTGGATTGTGCGGGGTGACTTTGAACCTGAATCCGGTTATCGCGCCATGCAGCAAATCCTGTCGCAGGCGCATCGCCCTACTGcc                >  2:3503817/1‑93 (MQ=255)                                                                        ccGGAAAGCTGGATTGTGCGGGGTGACTTTGAACCTGAATCCGGTTATCGCGCCATGCAGCAAATCCTGTCGCAGCCGCATCGCCCTACTGcc                <  2:2446957/93‑1 (MQ=255)                                                                          ggAAAGCTGGATTGTGCGGGGTGACTTTGAACCTGAATCCGGTTATCGCGCCATGCAGCAAATCCTGTCGCAGCCGCATCGCCCTACTGCCGt              >  1:2357710/1‑93 (MQ=255)                                                                           gAAAGCTGGATTGTGCGGGGTGACTTTGAACCTGAATCCGGTTATCGCGCCATGCAGCAAATCCTGTCGCAGCCGCATCGCCCTACTGCCGtc             >  2:1581849/1‑93 (MQ=255)                                                                            aaaGCTGGATTGTGCGGGGTGACTTTGAACCTGAATCCGGTTATCGCGCCATGCAGCAAATCCTGTCGCAGCCGCATCGCCCTACTGCCGtct            <  2:908668/93‑1 (MQ=255)                                                                              aGCTGGATTGTGCGGGTTGACTTTGAACCTGAATCCGGTTATCGCGCCATGCAGCAAATCCTGTCGCAGCCGCATCGCCCTACTGCCGtcttc          <  1:218990/93‑1 (MQ=255)                                                                              aGCTGGATTGTGCGGGGTGACTTTGAACCTGAATCCGGTTATCGCGCCATGCAGCAAATCCTGTCGCAGCCGCATCGCCCTACTGCCGtcttc          <  1:4305630/93‑1 (MQ=255)                                                                               gCTGGATTGTGCGGGGTGACTTTGAACCTGAATCCGGTTATCGCGCCATGCAGCAAATCCTGTCGCAGCCGCATCGCCCTACTGCCGtcttct         >  1:2131261/1‑93 (MQ=255)                                                                               gCTGGATTGTGCGGGGTGACTTTGAACCTGAATCCGGTTATCGCGCCATGCAGCAAATCCTGTCGCAGCCGCATCGCCCTACTGCCGtcttct         >  2:2380282/1‑93 (MQ=255)                                                                                  ggATTGTGCGGGGTGACTTTGAACCTGAATCCGGTTATCGCGCCATGCAGCAAATCCTGTCGCAGCCGCATCGCCCTACTGCCGTCTTCtgtg      <  2:3421646/93‑1 (MQ=255)                                                                                   gATTGTGCGGGGTGACTTTGAACCTGAATCCGGTTATCGCGCCATGCAGCAAATCCTGTCGCAGCCGCATCGCCCTACTGCCGTCTTCTgtgg     <  1:3261796/93‑1 (MQ=255)                                                                                     ttGTGCGGGGTGACTTTGAACCTGAATCCGGTTATCGCGCCATGCAGCAAATCCTGTCGCAGCCGCATCGCCCTACTGCCGTCTTCTgtggtg   >  2:3964987/1‑93 (MQ=255)                                                                                      tgtgCGGGGTGACTTTGAACCTGAATCCGGTTATCGCGCCATGCAGCAAATCCTGTCGCAGCCGCATCGCCCTACTGCCGTCTTCTGaggtgg  <  2:4530956/93‑1 (MQ=255)                                                                                           |                                                                                        TGGAACGTAACACCGGCGCAGGCCGCCTTGCCGGTTTTATGAAGGCGATGGAAGAAGCGATGATCAAGGTGCCGGAAAGCTGGATTGTGCAGGGTGACTTTGAACCTGAATCCGGTTATCGCGCCATGCAGCAAATCCTGTCGCAGCCGCATCGCCCTACTGCCGTCTTCTGTGGTGG  >  NC_000913/1736430‑1736607

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