Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	504,712	(T)5→4	intergenic (‑16/‑202)	fsr ← / → ushA	fosmidomycin efflux pump/5'‑nucleotidase/UDP‑sugar hydrolase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	504,708	0	T	.	100.0%	124.6 / NA	30	intergenic (‑12/‑206)	fsr/ushA	fosmidomycin efflux pump/5'‑nucleotidase/UDP‑sugar hydrolase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base . (12/18);  total (12/18)

GATGTTCGGGCCTTGGTCGTTGACGCAGCCGCGCCCGCCACAGGCTGGGGTTGTTCACTCATTGCCATAGGAAAGTCACTTTTTCAGGGTTGCGATGTAAAGAATGATCTTATTTGTGATTATTACCAGACTAACATACCTGTATGCGTCGTCTG  >  NC_000913/504629‑504783                                                                                |                                                                            gATGTTCGGGCCTTGGTCGTTGACGCAGCCGCGCCCGCCACAGGCTGGGGTTGTTCACTCATTGCCATAGGAAAGTCACtt                                                                            >  1:112492/1‑81 (MQ=255)   tGTTCGGGCCTTGGTCGTTGACGCAGCCGCGCCCGCCACAGGCTGGGGTTGTTCACTCATTGCCATAGGAAAGTCACtttt                                                                          <  2:112457/81‑1 (MQ=255)   tGTTCGGGCCTTGGTCGTTGACGCAGCCGCGCCCGCCACAGGCTGGGGTTGTTCACTCATTGCCATAGGAAAGTCACtttt                                                                          <  1:112458/81‑1 (MQ=255)     ttCGGGCCTTGGTCGTTGACGCAGCCGCGCCCGCCACAGGCTGGGGTTGTTCACTCATTGCCATAGGAAAGTCACtttt                                                                          >  2:112491/1‑79 (MQ=255)       cGGGCCTTGGTCGTTGACGCAGCCGCGCCCGCCACAGGCTGGGGTTGTTCACTCATTGCCATAGGAAAGTCACttttcagg                                                                      >  2:112505/1‑77 (MQ=255)       cGGGCCTTGGTCGTTGACGCAGCCGCGCCCGCCACAGGCTGGGGTTGTTCACTCATTGCCATAGGAAAGTCACttttcagg                                                                      >  2:112511/1‑77 (MQ=255)         ggCCTTGGTCGTTGACGCAGCCGCGCCCGCCACAGGCTGGGGTTGTTCACTCATTGCCATAGGAAAGTCAC‑TTTTCAGGGt                                                                   <  2:112459/81‑1 (MQ=255)           ccTTGGTCGTTGACGCAGCCGCGCCCGCCACAGGCTGGGGTTGTTCACTCATTGCCATAGGAAAGTCAC‑TTTTCAGGGTTg                                                                 <  1:112460/81‑1 (MQ=255)            cTTGGTCGTTGACGCAGCCGCGCCCGCCACAGGCTGGGGTTGTTCACTCATTGCCATAGGAAAGTCAC‑TTTTCAGGGTTGc                                                                <  2:112461/81‑1 (MQ=255)              tGGTCGTTGACGCAGCCGCGCCCGCCACAGGCTGGGGTTGTTCACTCATCGCCATAGGAAAGTCAC‑TTTTCAGGGTTGCGa                                                              >  2:112500/1‑81 (MQ=255)                gTCGTTGACGCAGCCGCGCCCGCCACAGGCTGGGGTTGTTCACTCATTGCCATAGGAAAGTCAC‑TTTTCAGGGTTGCGAt                                                             <  2:112463/80‑1 (MQ=255)                gTCGTTGACGCAGCCGCGCCCGCCACAGGCTGGGGTTGTTCACTCATTGCCATAGGAAAGTCAC‑TTTTCAGGGTTGCGATg                                                            <  1:112462/81‑1 (MQ=255)                   gTTGACGCAGCCGCGCCCGCCACAGGCTGGGGTTGTTCACTCATTGCCATAGGAAAGTCAC‑TTTTCAGGGTTGCGATGTa                                                          <  1:112464/80‑1 (MQ=255)                    ttGACGCAGCCGCGCCCGCCACAGGCTGGGGTTGTTCACTCATTGCCATAGGAAAGTCAC‑TTTTCAGGGTTGCGATGTaaa                                                        <  1:112465/81‑1 (MQ=255)                      gACGCAGCCGCGCCCGCCACAGGCTGGGGTTGTTCACTCATTGCCATAGGAAAGTCAC‑TTTTCAGGGTTGCGATGTAAAg                                                       >  2:112507/1‑80 (MQ=255)                        cGCAGCCGCGCCCGCCACAGGCTGGGGTTGTTCACTCATTGCCATAGGAAAGTCAC‑TTTTCAGGGTTGCGATGTAAAGAAt                                                    <  1:112466/81‑1 (MQ=255)                        cGCAGCCGCGCCCGCCACAGGCTGGGGTTGTTCACTCATTGCCATAGGAAAGTCAC‑TTTTCAGGGTTGCGATGTAAAGAAt                                                    >  2:112503/1‑81 (MQ=255)                            gCCGCGCCCGCCACAGGCTGGGGTTGTTCACTCATTGCCATAGGAAAGTCAC‑TTTTCAGGGTTGCGATGTAAAGAATGATc                                                <  1:112467/81‑1 (MQ=255)                             ccGCGCCCGCCACAGGCTGGGGTTGTTCACTCATTGCCATAGGAAAGTCAC‑TTTTCAGGGTTGCGATGTAAAGAATGATCt                                               <  1:112468/81‑1 (MQ=255)                              cgcgCCCGCCACAGGCTGGGGTTGTTCACTCATTGCCATAGGAAAGTCAC‑TTTTCAGGGTTGCGATGTAAAGAATGATCtt                                              <  1:112469/81‑1 (MQ=255)                              cgcgCCCGCCACAGGCTGGGGTTGTTCACTCATTGCCATAGGAAAGTCAC‑TTTTCAGGGTTGCGATGTAAAGAATGATCt                                               >  1:112501/1‑80 (MQ=255)                               gcgcCCGCCACAGGCTGGGGTTGTTCACTCATTGCCATAGGAAAGTCAC‑TTTTCAGGGTTGCGATGTAAAGAATGATCTTa                                             >  1:112520/1‑81 (MQ=255)                                   ccgccACAGGCTGGGGTTGTTCACTCATTGCCATAGGAAAGTCAC‑TTTTCAGGGTTGCGATGTAAAGAATGATCTTATTtg                                         >  1:112498/1‑81 (MQ=255)                                             cTGGGGTTGTTCACTCATTGCCATAGGAAAGTCAC‑TTTTCAGGGTTGCGATGTAAAGAATGATCTTATTTGTGATTATTAc                               <  1:112470/81‑1 (MQ=255)                                             cTGGGGTTGTTCACTCATTGCCATAGGAAAGTCAC‑TTTTCAGGGTTGCGATGTAAAGAATGATCTTATTTGTGATTATTAc                               >  2:112502/1‑81 (MQ=255)                                                gggTTGTTCACTCATTGCCATAGGAAAGTCAC‑TTTTCAGGGTTGCGATGTAAAGAATGATCTTATTTGTGATTATTACCAg                            <  1:112471/81‑1 (MQ=255)                                                 ggTTGTTCACTCATTGCCATAGGAAAGTCAC‑TTTTCAGGGTTGCGATGTAAAGAATGATCTTATTTGTGATTATTACCAGa                           >  2:112519/1‑81 (MQ=255)                                                 ggTTGTTCACTCATTGCCATAGGAAAGTCAC‑TTTTCAGGGTTGCGATGTAAAGAATGATCTTATTTGTGATTATTACCAGa                           <  1:112472/81‑1 (MQ=255)                                                     gttCACTCATTGCCATAGGAAAGTCAC‑TTTTCAGGGTTGCGATGTAAAGAATGATCTTATTTGTGATTATTACCAGACTa                        <  2:112473/80‑1 (MQ=255)                                                         aCTCATTGCCATAGGGAAGTCAC‑TTTTCAGGGTTGCGATGTAAAGAATGATCTTATTTGTGATTATTACCAGACTAACAt                    <  1:112474/80‑1 (MQ=255)                                                                gCCATAGGAAAGTCAC‑TTTTCAGGGTTGCGATGTAAAGAATGATCTTATTTGTGATTATTACCAGACTAACATACCTGTa             <  2:112475/80‑1 (MQ=255)                                                                  cATAGGAAAGTCAC‑TTTTCAGGGTTGCGATGTAAAGAATGATCTTATTTGTGATTATTACCAGACTAACATACCTGTATGc          >  2:112504/1‑81 (MQ=255)                                                                    tAGGAAAGTCAC‑TTTTCAGGGTTGCGATGTAAAGAATGATCTTATTTGTGATTATTACCAGACTAACATACCTGTATGcgt        >  2:112509/1‑81 (MQ=255)                                                                       gAAAGTCAC‑TTTTCAGGGTTGCGATGTAAAGAATGATCTTATTTGTGATTATTACCAGACTAACATACCTGTATGcgtcgt     >  1:112523/1‑81 (MQ=255)                                                                         aaGTCAC‑TTTTCAGGGTTGCGATGTAAAGAATGATCTTATTTGTGATTATTACCAGACTAACATACCTGTATGCGTCGTCt   <  1:112476/81‑1 (MQ=255)                                                                           gTCAC‑TTTTCAGGGTTGCGATGTAAAGAATGATCTTATTTGTGATTATTACCAGACTAACATACCTGTATGCGTCGTCTg  >  2:112510/1‑80 (MQ=255)                                                                                |                                                                            GATGTTCGGGCCTTGGTCGTTGACGCAGCCGCGCCCGCCACAGGCTGGGGTTGTTCACTCATTGCCATAGGAAAGTCACTTTTTCAGGGTTGCGATGTAAAGAATGATCTTATTTGTGATTATTACCAGACTAACATACCTGTATGCGTCGTCTG  >  NC_000913/504629‑504783

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