Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,983,935	C→A	A379S (GCC→TCC)	araG ←	arabinose ABC transporter ATP binding subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,983,935	0	C	A	100.0%	108.5 / NA	30	A379S (GCC→TCC)	araG	arabinose ABC transporter ATP binding subunit
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base A (20/10);  total (20/10)

ATTTCCGCCTGAGAGATTCATGATCAGTTGCTCCGCGCCCGGCGTTTTGATGTTGAGCGAACGAATGTGGTGATCGGCATTGTTTTCTTCCCAACCGTTGTTGATTACACAACCGCCGAGCACATGTTTACGTCTGGCACTGATGTTGATATTGT  >  NC_000913/1983858‑1984012                                                                              |                                                                              aTTTCCGCCTGAGAGATTCATGATCAGTTGCTCCGCGCCCGGCGTTTTGATGTTGAGCGAACGAATGTGGTGATCGGAAtt                                                                            <  1:348641/81‑1 (MQ=255) aTTTCCGCCTGAGAGATTCATGATCAGTTGCTCCGCGCCCGGCGTTTTGATGTTGAGCGAACGAATGTGGTGATCGGAAtt                                                                            >  2:348687/1‑81 (MQ=255)   ttCCGCCTGAGAGATTCATGATCAGTTGCTCCGCGCCCGGCGTTTTGATGTTGAGCGAACGAATGTGGTGATCGGAATTGt                                                                          >  1:348677/1‑81 (MQ=255)     ccGCCTGAGAGATTCATGATCAGTTGCTCCGCGCCCGGCGTTTTGATGTTGAGCGAACGAATGTGGTGATCGGAATTGtt                                                                         <  2:348642/80‑1 (MQ=255)             gagaTTCATGATCAGTTGCTCCGCGCCCGGCGTTTTGATGTTGAGCGAACGAATGTGGTGATCGGAATTGTTTTCTTCCCa                                                                <  2:348643/81‑1 (MQ=255)               gaTTCATGATCAGTTGCTCCGCGCCCGGCGTTTTGATGTTGAGCGAACGAATGTGGTGATCGGAATTGTTTTCTTCCCAAc                                                              >  2:348679/1‑81 (MQ=255)                aTTCATGATCAGTTGCTCCGCGCCCGGCGTTTTGATGTTGAGCGAACGAATGTGGTGATCGGAATTGTTTTCTTCCCAAcc                                                             <  1:348644/81‑1 (MQ=255)                aTTCATGATCAGTTGCTCCGCGCCCGGCGTTTTGATGTTGAGCGAACGAATGTGGTGATCGGAATTGTTTTCTTCCCAAcc                                                             >  1:348676/1‑81 (MQ=255)                 ttCATGATCAGTTGCTCCGCGCCCGGCGTTTTGATGTTGAGCGAACGAATGTGGTGATCGGAATTGTTTTCTTCCCAACCg                                                            >  1:348688/1‑81 (MQ=255)                         cAGTTGCTCCGCGCCCGGCGTTTTGATGTTGAGCGAACGAATGTGGTGATCGGAATTGTTTTCTTCCCAACCGTTGTTGAt                                                    >  1:348685/1‑81 (MQ=255)                           gTTGCTCCGCGCCCGGCGTTTTGATGTTGAGCGAACGAATGTGGTGATCGGAATTGTTTTCTTCCCAACCGTTGTTGATTa                                                  >  1:348690/1‑81 (MQ=255)                               cTCCGCGCCCGGCGTTTTGATGTTGAGCGAACGAATGTGGTGATCGGAATTGTTTTCTTCCCAACCGTTGTTGATTacac                                               <  1:348645/80‑1 (MQ=255)                                  cgcgCCCGGCGTTTTGATGTTGAGCGAACGAATGTGGTGATCGGAATTGTTTTCTTCCCAACCGTTGTTGATTACACAAc                                            <  1:348646/80‑1 (MQ=255)                                   gcgcCCGGCGTTTTGATGTTGAGCGAACGAATGTGGTGATCGGAATTGTTTTCTTCCCAACCGTTGTTGATTACACAAtcg                                          >  2:348678/1‑81 (MQ=255)                                      cccGGCGTTTTGATGTTGAGCGAACGAATGTGGTGATCGGAATTGTTTTCTTCCCAACCGTTGTTGATTACACAAccgcc                                        >  2:348682/1‑80 (MQ=255)                                      cccGGCGTTTTGATGTTGAGCGAACGAATGTGGTGATCGGAATTGTTTTCTTCCCAACCGTTGTTGATTACACAAccgcc                                        >  2:348675/1‑80 (MQ=255)                                       ccGGCGTTTTGATGTTGAGCGAACGAATGTGGTGATCGGAATTGTTTTCTTCCCAACCGTTGTTGATTACACAACCGCCGa                                      <  1:348647/81‑1 (MQ=255)                                               ttGATGTTGAGCGAACGAATGTGGTGATCGGAATTGTTTTCTTCCCAACCGTTGTTGATTACACAACCGCCGAGCACATGt                              >  2:348684/1‑81 (MQ=255)                                                 gATGTTGAGCGAACGAATGTGGTGATCGGAATTGTTTTCTTCCCAACCGTTGTTGATTACACAACCGCCGAGCACATGtt                             >  2:348674/1‑80 (MQ=255)                                                  aTGTTGAGCGAACGAATGTGGTGATCGGAATTGTTTTCTTCCCAACCGTTGTTGATTACACAACCGCCGAGCACATGttt                            <  2:348648/80‑1 (MQ=255)                                                    gTTGAGCGAACGAATGTGGTGATCGGAATTGTTTTCTTCCCAACCGTTGTTGATTACACAACCGCCGAGCACATGTTTACg                         >  1:348689/1‑81 (MQ=255)                                                        aGCGAACGAATGTGGTGATCGGAATTGTTTTCTTCCCAACCGTTGTTGATTACACAACCGCCGAGCACATGTTTACGTCTg                     >  2:348691/1‑81 (MQ=255)                                                           gaacgaaTGTGGTGATCGGAATTGTTTTCTTCCCAACCGTTGTTGATTACACAACCGCCGAGCACATGTTTACGTCTGGCa                  >  2:348681/1‑81 (MQ=255)                                                              cgaaTGTGGTGATCGGAATTGTTTTCTTCCCAACCGTTGTTGATTACACAACCGCCGAGCACATGTTTACGTCTGGCACt                >  2:348680/1‑80 (MQ=255)                                                              cgaaTGTGGTGATCGGAATTGTTTTCTTCCCAACCGTTGTTGATTACACAACCGCCGAGCACATGTTTACGTCTGGCACTg               <  2:348649/81‑1 (MQ=255)                                                               gaaTGTGGTGATCGGAATTGTTTTCTTCCCAACCGTTGTTGATTACACAACCGCCGAGCACATGTTTACGTCTGGCACTGa              >  1:348692/1‑81 (MQ=255)                                                               gaaTGTGGTGATCGGAATTGTTTTCTTCCCAACCGTTGTTGATTACACAACCGCCGAGCACATGTTTACGTCTGGCACTGa              >  2:348701/1‑81 (MQ=255)                                                                      gtgATCGGAATTGTTTTCTTCCCAACCGTTGTTGATTACACAACCGCCGAGCACATGTTTACGTCTGGCACTGATGTTGat       <  2:348650/81‑1 (MQ=255)                                                                       tgATCGGAATTGTTTTCTTCCCAACCGTTGTTGATTACACAACCGCCGAGCACATGTTTACGTCTGGCACTGATGTTGata      >  2:348693/1‑81 (MQ=255)                                                                           cGGAATTGTTTTCTTCCCAACCGTTGTTGATTACACAACCGCCGAGCACATGTTTACGTCTGGCACTGATGTTGATATTGt  >  1:348695/1‑81 (MQ=255)                                                                              |                                                                              ATTTCCGCCTGAGAGATTCATGATCAGTTGCTCCGCGCCCGGCGTTTTGATGTTGAGCGAACGAATGTGGTGATCGGCATTGTTTTCTTCCCAACCGTTGTTGATTACACAACCGCCGAGCACATGTTTACGTCTGGCACTGATGTTGATATTGT  >  NC_000913/1983858‑1984012

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