Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_017635	3,053,299	(A)8→9	intergenic (‑32/‑286)	ygcW ← / → yqcE	SDR family oxidoreductase/MFS transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_017635	3,053,291	1	.	A	100.0%	134.6 / NA	37	intergenic (‑24/‑294)	ygcW/yqcE	SDR family oxidoreductase/MFS transporter
Reads supporting (aligned to +/- strand):  ref base . (0/0);  new base A (17/20);  total (17/20)

AACGCATTGAGAGATTCGATTGACATAGTTAGCTCCAATTTCCGTAAGGC‑AAAAAAAAGAAAGGTAACTTCCCCCGCGAGAGGGAAGTTACCT  >  NC_017635/3053242‑3053334                                                   |                                            aaCGCATTGAGAGATTCGATTGACATAGTTAGCTCCAATTTCCGTAAGGc                                              <  2:1810405/50‑1 (MQ=255)  aCGCATTGAGAGATTCGATTGACATAGTTAGCTCCAATTTCCGTAAGGc                                              >  1:1075239/1‑49 (MQ=255)   cGCATTGAGAGATTCGATTGACATAGTTAGCTCCAATTTCCGTAAGGC‑aa                                           <  1:1058007/50‑1 (MQ=255)   cGCATTGAGAGATTCGATTGACATAGTTAGCTCCAATTTCCGTAAGGC‑a                                            >  1:1704265/1‑49 (MQ=255)    gCATTGAGAGATTCGATTGACATAGTTAGCTCCAATTTCCGTAAGGC‑aaa                                          <  2:338873/50‑1 (MQ=255)    gCATTGAGAGATTCGATTGACATAGTTAGCTCCAATTTCCGTAAGGC‑aa                                           >  1:2062774/1‑49 (MQ=255)      aTTGAGAGATTCGATTGACATAGTTAGCTCCAATTTCCGTAAGGC‑aaaa                                         >  2:288765/1‑49 (MQ=255)      aTTGAGAGATTCGATTGACATAGTTAGCTCCAATTTCCGGAAGGC‑aaaa                                         >  1:2034919/1‑49 (MQ=255)       ttGAGAGATTCGATTGACATAGTTAGCTCCAATTTCCGTAAGGC‑aaaaa                                        >  1:1734688/1‑49 (MQ=255)       ttGAGAGATTCGATTGACATAGTTAGCTCCAATTTCCGTAAGGC‑aaaaa                                        >  1:509582/1‑49 (MQ=255)        tGAGAGATTCGATTGACATAGTTAGCTCCAATTTCCGTAAGGC‑aaaaaa                                       >  2:1295178/1‑49 (MQ=255)        tGAGAGATTCGATTGACATAGTTAGCTCCAATTTCCGTAAGGC‑aaaaaa                                       >  1:1894099/1‑49 (MQ=255)         gagagaTTCGATTGACATAGTTAGCTCCAATTTCCGTAAGGC‑aaaaaaa                                      >  1:1091445/1‑49 (MQ=255)         gagagaTTCGATTGACATAGTTAGCTCCAATTTCCGTAAGGC‑aaaaaaa                                      <  1:407119/49‑1 (MQ=255)          agagaTTCGATTGACATAGTTAGCTCCAATTTCCGTAAGGC‑aaaaaaa                                      >  2:1399527/1‑48 (MQ=255)            agaTTCGATTGACATAGTTAGCTCCAATTTCCGTAAGGC‑aaaaaaaaag                                   >  1:648616/1‑47 (MQ=255)            agaTTCGATTGACATAGTTAGCTCCAATTTCCGTAAGGC‑aaaaaaaaag                                   >  2:1077418/1‑47 (MQ=255)               ttCGATTGACATAGTTAGCTCCAATTTCCGTAAGGCAAAAAAaaagaaag                                >  1:988114/1‑50 (MQ=255)               ttCGATTGACATAGTTAGCTCCAATTTCCGTAAGGCAAAAAAaaagaaag                                <  1:1750557/50‑1 (MQ=255)               ttCGATTGACATAGTTAGCTCCAATTTCCGTAAGGCAAAAAAaaagaaag                                <  2:441617/50‑1 (MQ=255)                tCGATTGACATAGTTAGCTCCAATTTCCGTAAGGCAAAAAAAAAGAAAgg                               >  2:1398582/1‑50 (MQ=255)                tCGATTGACATAGTTAGCTCCAATTTCCGTAAGGCAAAAAAAAAGAAAgg                               >  1:428093/1‑50 (MQ=255)                 cGATTGACATAGTTAGCTCCAATTTCCGTAAGGCAAAAAAAAAGAAAgg                               <  2:554192/49‑1 (MQ=255)                  gATTGACATAGTTAGCTCCAATTTCCGTAAGGCAAAAAAAAAGAAAGGt                              >  1:218107/1‑49 (MQ=255)                  gATTGACATAGTTAGCTCCAATTTCCGTAAGGCAAAAAAAAAGAAAGGTa                             <  1:1241354/50‑1 (MQ=255)                   aTTGACATAGTTAGCTCCAATTTCCGTAAGGCAAAAAAAAAGAAAGGTaa                            >  2:1767462/1‑50 (MQ=255)                   aTTGACATAGTTAGCTCCAATTTCCGTAAGGCAAAAAAAAAGAAAGGTaa                            <  1:1030315/50‑1 (MQ=255)                     tGACATAGTTAGCTCCAATTTCCGTAAGGCAAAAAAAAAGAAAGGTAACt                          >  1:1342899/1‑50 (MQ=255)                      gACATAGTTAGCTCCAATTTCCGTAAGGCAAAAAAAAAGAAAGGTAACt                          <  1:1613297/49‑1 (MQ=255)                        cATAGTTAGCTCCAATTTCCGTAAGGCAAAAAAAAAGAAAGGTAACTTcc                       <  2:269449/50‑1 (MQ=255)                        cATAGTTAGCTCCAATTTCCGTAAGGCAAAAAAAAAGAAAGGTAACTTcc                       >  2:274066/1‑50 (MQ=255)                        cATAGTTAGCTCCAATTTCCGTAAGGCAAAAAAAAAGAAAGGTAACTTcc                       <  1:314929/50‑1 (MQ=255)                        cATAGTTAGCTCCAATTTCCGTAAGGCAAAAAAAAAGAAAGGTAACTTcc                       <  1:1913649/50‑1 (MQ=255)                        cATAGTTAGCTCCAATTTCCGTAAGGCAAAAAAAAAGAAAGGTAACTTcc                       <  2:716128/50‑1 (MQ=255)                        cATAGTTAGCTCCAATTTCCGTAAGGCAAAAAAAAAGAAAGGTAACTTcc                       >  2:786171/1‑50 (MQ=255)                        cATAGTTAGCTCCAATTTCCGTAAGGCAAAAAAAAAGAAAGGTAACTTcc                       <  2:1920647/50‑1 (MQ=255)                         aTAGTTAGCTCCAATTTCCGTAAGGCAAAAAAAAAGAAAGGTAACTTccc                      <  2:60940/50‑1 (MQ=39)                          tAGTTAGCTCCAATTTCCGTAAGGCAAAAAAAAAGAAAGGTAACTTcccc                     <  2:676825/50‑1 (MQ=39)                          tAGTTAGCTCCAATTTCCGTAAGGCAAAAAAAAAGAAAGGTAACTTcccc                     <  1:1636614/50‑1 (MQ=39)                           aGTTAGCTCCAATTTCCGTAAGGCAAAAAAAAAGAAAGGTAACTTccccc                    >  2:88946/1‑50 (MQ=39)                           aGTTAGCTCCAATTTCCGTAAGGCAAAAAAAAAGAAAGGTAACTTccccc                    >  1:696362/1‑50 (MQ=39)                             ttAGCTCCAATTTCCGTAAGGCAAAAAAAAAGAAAGGTAACTTCCCCcgc                  <  2:1476972/50‑1 (MQ=39)                             ttAGCTCCAATTTCCGTAAGGCAAAAAAAAAGAAAGGTAACTTCCCCcgc                  <  2:446523/50‑1 (MQ=39)                               aGCTCCAATTTCCGTAAGGCAAAAAAAAAGAAAGGTAACTTCCCCcgc                  >  2:1189195/1‑48 (MQ=38)                                 cTCCAATTTCCGTAAGGCAAAAAAAAAGAAAGGTAACTTCCCCCGCga                >  1:948945/1‑48 (MQ=38)                                    cAATTTCCGTAAGGCAAAAAAAAAGAAAGGTAACTTCCCCCGCGAGAggg           <  1:749859/50‑1 (MQ=39)                                    cAATTTCCGTAAGGCAAAAAAAAAGAAAGGTAACTTCCCCCGCGAGAggg           <  1:1798185/50‑1 (MQ=39)                                      aTTTCCGTAAGGCAAAAAAAAAGAAAGGTAACTTCCCCCGCGAGAggg           >  1:2025374/1‑48 (MQ=38)                                      aTTTCCGTAAGGCAAAAAAAAAGAAAGGTAACTTCCCCCGCGAGAggg           >  1:1161830/1‑48 (MQ=38)                                      aTTTCCGTAAGGCAAAAAAAAAGAAAGGTAACTTCCCCCGCGAGAGGGaa         <  1:1646397/50‑1 (MQ=38)                                           cGTAAGGCAAAAAAAAAGAAAGGTAACTTCCCCCGCGAGAGGGAAGTTAc    >  1:346534/1‑50 (MQ=37)                                             tAAGGCAAAAAAAAAGAAAGGTAACTTCCCCCGCGAGAGGGAAGTTACCt  >  2:256657/1‑50 (MQ=25)                                             tAAGGCAAAAAAAAAGAAAGGTAACTTCCCCCGCGAGAGGGAAGTTACCt  <  1:998280/50‑1 (MQ=25)                                             tAAGGCAAAAAAAAAGAAAGGAAACTTCCCCCGCGAGAGGGAAGTTACCt  >  2:1804059/1‑50 (MQ=25)                                                   |                                            AACGCATTGAGAGATTCGATTGACATAGTTAGCTCCAATTTCCGTAAGGC‑AAAAAAAAGAAAGGTAACTTCCCCCGCGAGAGGGAAGTTACCT  >  NC_017635/3053242‑3053334

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