Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,618,218	G→A	*397* (TAG→TAA)	ydeA →	arabinose efflux transporter, arabinose‑inducible

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,618,218	0	G	A	100.0%	132.5 / NA	42	*397* (TAG→TAA)	ydeA	arabinose efflux transporter, arabinose‑inducible
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (20/22);  total (20/22)

GGTGAGTTTGCACTGGTCAATGTCGATGATTGGTTATGTGGGCGCGGTGCCTGCTTTTGCCGCGTTAATTTGGTCAATCATTATATTTCGCCGCTGGCCAGTGACACTCGAAGAACAGACGCAATAGTTGAAAGGCCCATTCGGGCCTTTTTTAATGGTACGTTTTAATGATTTCCAGGATGCCGTTAATAATAAACTGCACACCCATACATACCAGCAGGAATCCCATCAGACGGGAGATCGCTTCAATGC  >  NC_000913/1618092‑1618343                                                                                                                               |                                                                                                                              ggTGAGTTTGCACTGGTCAATGTCGATGATTGGTTATGTGGGCGCGGTGCCTGCTTTTGCCGCGTTAATTTGGTCAATCATTATATTTCGCCGCTGGCCAGTGACACTCGAAGAACAGACGCAATAATTGAAAGGCCCa                                                                                                                   >  1:46239/1‑139 (MQ=255)     aGTTTGCACTGGTCAATGTCGATGATTGGTTATGTGGGCGCGGTGCCTGCTTTTGCCGCGTTAATTTGGTCAATCATTATATTTCGCCGCTGGCCAGTGACACTCGAAGAACAGACGCAATAATTGAAAGGCCCATTCg                                                                                                               >  1:309731/1‑139 (MQ=255)         tGCACTGGTCAATGTCGATGATTGGTTATGTGGGCGCGGTGCCTGCTTTTGCCGCGTTAATTTGGTCAATCATTATATTTCGCCGCTGGCCAGTGACACTCGAAGAACAGACGCAATAATTGAAAGGCCCATTCGGGcc                                                                                                           >  2:419857/1‑139 (MQ=255)                gTCAATGTCGATGATTGGTTATGTGGGCGCGGTGCCTGCTTTTGCCGCGTTAATTTGGTCAATCATTATATTTCGCCGCTGGCCAGTGACACTCGAAGAACAGACGCAATAATTGAAAGGCCCATTCGGGcc                                                                                                           >  1:105867/1‑132 (MQ=255)                gTCAATGTCGATGATTGGTTATGTGGGCGCGGTGCCTGCTTTTGCCGCGTTAATTTGGTCAATCATTATATTTCGCCGCTGGCCAGTGACACTCGAAGAACAGACGCAATAATTGAAAGGCCCATTCGGGcc                                                                                                           <  2:105867/132‑1 (MQ=255)                     tGTCGATGATTGGTTATGTGGGCGCGGTGCCTGCTTTTGCCGCGTTAATTTGGTCAATCATTATATTTCGCCGCTGGCCAGTGACACTCGAAGAACAGACGCAATAATTGAAAGGCCCATTCGGGCCTTTTTTAATGGt                                                                                               >  1:511146/1‑139 (MQ=255)                              ttGGTTATGTGGGCGCGGTGCCTGCTTTTGCCGCGTTAATTTGGTCAATCATTATATTTCGCCGCTGGCCAGTGACACTCGAAGAACAGACGCAATAATTGAAAGGCCCATTCGGGCCTTTTTTAATGGTACGTTTTaa                                                                                      <  2:46239/139‑1 (MQ=255)                                 gTTATGTGGGCGCGGTGCCTGCTTTTGCCGCGTTAATTTGGTCAATCATTATATTTCGCCGCTGGCCAGTGACACTCGAAGAACAGACGCAATAATTGAAAGGCCCATTCGGGCCttt                                                                                                        >  2:465725/1‑118 (MQ=255)                                 gTTATGTGGGCGCGGTGCCTGCTTTTGCCGCGTTAATTTGGTCAATCATTATATTTCGCCGCTGGCCAGTGACACTCGAAGAACAGACGCAATAATTGAAAGGCCCATTCGGGCCttt                                                                                                        <  1:465725/118‑1 (MQ=255)                                         ggCGCGGTGCCTGCTTTTGCCGCGTTAATTTGGTCAATCATTATATTTCGCCGCTGGCCAGTGACACTCGAAGAACAGACGCAATAATTGAAAGGCCCATTCGGGCCTTTTTTAATGGTACGTTTTAATGATTTCCAgg                                                                           <  2:52269/139‑1 (MQ=255)                                           cgcgGTGCCTGCTTTTGCCGCGTTAATTTGGTCAATCATTATATTTCGCCGCTGGCCAGTGACACTCGAAGAACAGACGCAATAATTGAAAGGCCCATTCgg                                                                                                              >  2:38581/1‑102 (MQ=255)                                           cgcgGTGCCTGCTTTTGCCGCGTTAATTTGGTCAATCATTATATTTCGCCGCTGGCCAGTGACACTCGAAGAACAGACGCAATAATTGAAAGGCCCATTCgg                                                                                                              <  1:38581/102‑1 (MQ=255)                                           cgcgGTGCCTGCTTTTGCCGCGTTAATTTGGTCAATCATTATATTTCGCCGCTGGCCAGTGACACTCGAAGAACAGACGCAATAATTGAAAGGCCCATTCGGGCCTTTTTTAATGGTACGTTTTAATGATTTCCAGGAt                                                                         >  2:283405/1‑139 (MQ=255)                                                             cgcgTTAATTTGCTCAATCATTATATTTCGCCGCTGGCCAGTGACACTCGAAGCACCGACTCAATAATTTAAAGGCCCATTCGCGCCTTTTTTAATGGTACGTTTTAATGATTTCCAGGATGCCGTTAATAATAAACTg                                                       <  1:392876/139‑1 (MQ=255)                                                              gcgTTAATTTGGTCAATCATTATATTTCGCCGCTGGCCAGTGACACTCGAAGAACAGACGCAATAATTGAAAGGCCCATTCGGGCCTTTTTTAATGGTACGTTTTAATGATTTCCAGGATGCCGTTAATAATAAACTGc                                                      <  2:309731/139‑1 (MQ=255)                                                               cgTTAATTTGGTCAATCATTATATTTCGCCGCTGGCCAGTGACACTCGAAGAACAGACGCAATAATTGAAAGGCCCATTCGGGCCTTTTTTAATGGTACGTTTTAATGATTTCCAGGATGCCGTTAATAATAAACTGca                                                     <  2:102713/139‑1 (MQ=255)                                                                                        tCGCCGCTGGCCAGTGACACTCGAAGAACAGACGCAATAATTGAAAGGCCCATTCGGGCCTTTTTTAATGGt                                                                                               >  1:304790/1‑72 (MQ=255)                                                                                        tCGCCGCTGGCCAGTGACACTCGAAGAACAGACGCAATAATTGAAAGGCCCATTCGGGCCTTTTTTAATGGt                                                                                               <  2:304790/72‑1 (MQ=255)                                                                                        tCGCCGCTGGCCAGTGACACTCGAAGAACAGACGCAATAATTGAAAGGCCCATTCGGGCCTTTTTTAATGGTACGTTTTAATGATTTCCAGGATGCCGt                                                                    <  1:354457/99‑1 (MQ=255)                                                                                        tCGCCGCTGGCCAGTGACACTCGAAGAACAGACGCAATAATTGAAAGGCCCATTCGGGCCTTTTTTAATGGTACGTTTTAATGATTTCCAGGATGCCGt                                                                    >  2:354457/1‑99 (MQ=255)                                                                                        tCGCCGCTGGCCAGTGACACTCGAAGAACAGACGCAATAATTGAAAGGCCCATTCGGGCCTTTTTTAATGGTACGTTTTAATGATTTCCAGGATGCCGTTAATAATAAACTGCACACCCATACATACCAGCAGGAATcc                            <  2:5462/139‑1 (MQ=255)                                                                                         cgccgcTGGCCAGTGACACTCGAAGAACAGACGCAATAATTGAAAGGCCCATTCGGGcc                                                                                                           >  2:645663/1‑59 (MQ=255)                                                                                         cgccgcTGGCCAGTGACACTCGAAGAACAGACGCAATAATTGAAAGGCCCATTCGGGcc                                                                                                           <  1:645663/59‑1 (MQ=255)                                                                                         cgccgcTGGCCAGTGACACTCGAAGAACAGACGCAATAATTGAAAGGCCCATTCGGGCCTTTTTTAATGGTACGTTTTAAt                                                                                     <  2:281868/81‑1 (MQ=255)                                                                                         cgccgcTGGCCAGTGACACTCGAAGAACAGACGCAATAATTGAAAGGCCCATTCGGGCCTTTTTTAATGGTACGTTTTAAt                                                                                     >  1:281868/1‑81 (MQ=255)                                                                                         cgccgcTGGCCAGTGACACTCGAAGAACAGACGCAATAATTGAAAGGCCCATTCGGGCCTTTTTTAATGGTACGTTTTAATGATTTCCAGGATg                                                                        >  2:387073/1‑94 (MQ=255)                                                                                         cgccgcTGGCCAGTGACACTCGAAGAACAGACGCAATAATTGAAAGGCCCATTCGGGCCTTTTTTAATGGTACGTTTTAATGATTTCCAGGATg                                                                        <  1:387073/94‑1 (MQ=255)                                                                                         cgccgcTGGCCAGTGACACTCGAAGAACAGACGCAATAATTGAAAGGCCCATTCGGGCCTTTTTTAATGGTACGTTTTAATGATTTCCAGGATGCCGTTAATAATAAACTg                                                       >  2:174012/1‑111 (MQ=255)                                                                                         cgccgcTGGCCAGTGACACTCGAAGAACAGACGCAATAATTGAAAGGCCCATTCGGGCCTTTTTTAATGGTACGTTTTAATGATTTCCAGGATGCCGTTAATAATAAACTg                                                       <  1:174012/111‑1 (MQ=255)                                                                                          gccgcTGGCCAGTGACACTCGAAGAACAGACGCAATAATTGAAAGGCCCATTCGGGCCTTTTTTAATGGTACGTTTTAATGATTTCCAGGATGcc                                                                      >  2:205719/1‑95 (MQ=255)                                                                                          gccgcTGGCCAGTGACACTCGAAGAACAGACGCAATAATTGAAAGGCCCATTCGGGCCTTTTTTAATGGTACGTTTTAATGATTTCCAGGATGcc                                                                      <  1:205719/95‑1 (MQ=255)                                                                                           ccgcTGGCCAGTGACACTCGAAGAACAGAGGCAATAATTGAAAGGCCCATTCGGGCCTTTTTTAATGGTa                                                                                              <  2:205952/70‑1 (MQ=255)                                                                                           ccgcTGGCCAGTGACACTCGAAGAACAGACGCAATAATTGAAAGGCCCATTCGGGCCTTTTTTAATGGTa                                                                                              >  1:205952/1‑70 (MQ=255)                                                                                           ccgcTGGCCAGTGACACTCGAAGAACAGACGCAATAATTGAAAGGCCCATTCGGGCCTTTTTTAATGGTACGTTTTAATGATTTCCAGGATGCCGTTAATAATAAACTGCACACCCATACATACCAGCAGGAATCCCAt                         <  2:511146/139‑1 (MQ=255)                                                                                           ccgcTGGCCAGTGACACTCGAAGAACAGACGCAATAATTGAAAGGCCCATTCGGGCCTTTTTTAATGGTACGTTTTAATGATTTCCAGGATGCCGTTAATAATAAACTGCACACCCATACATACCAGCAGGAATCCCAt                         >  2:514168/1‑139 (MQ=255)                                                                                              cTGGCCAGTGACACTCGAAGAACAGACGCAATAATTGAAAGGCCCATTCGGGCCTTTTTTAATGGTACGTTTTAATGATTTCCAGGATGCCGTTAATAATaaa                                                          >  2:273091/1‑103 (MQ=255)                                                                                              cTGGCCAGTGACACTCGAAGAACAGACGCAATAATTGAAAGGCCCATTCGGGCCTTTTTTAATGGTACGTTTTAATGATTTCCAGGATGCCGTTAATAATaaa                                                          <  1:273091/103‑1 (MQ=255)                                                                                              cTGGCCAGTGACACTCGAAGAACAGACGCAATAATTGAAAGGCCCATTCGGGCCTTTTTTAATGGTACGTTTTAATGATTTCCAGGATGCCGTTAATAATAAACTGCACACCCATACATACCAGCAGGAATCCCATCAg                      >  1:493385/1‑139 (MQ=255)                                                                                                           cTCGAAGAACAGACGCAATAATTGAAAGGCCCATTCGGGCCTTTTTTAATGGTACGTTTTAATGATTTCCAGGATGCCGTTAATAATAAACTGca                                                     >  1:515635/1‑95 (MQ=255)                                                                                                           cTCGAAGAACAGACGCAATAATTGAAAGGCCCATTCGGGCCTTTTTTAATGGTACGTTTTAATGATTTCCAGGATGCCGTTAATAATAAACTGca                                                     <  2:515635/95‑1 (MQ=255)                                                                                                                  aaCAGACGCAATAATTGAAAGGCCCATTCGGGCCTTTTTTAATGGTACGTTTTAATGATTTCCAGGATGCCGTTAATAATAAACTGCACACCCATACATACCAGCAGGAATCCCATCAGACGGGAGATCGCTTCAATGc  <  2:493385/139‑1 (MQ=255)                                                                                                                            aTAATTGAAAGGCCCATTCGGGCCTTTTTTAATGGTACGTTTTAAt                                                                                     <  2:658411/46‑1 (MQ=255)                                                                                                                               |                                                                                                                              GGTGAGTTTGCACTGGTCAATGTCGATGATTGGTTATGTGGGCGCGGTGCCTGCTTTTGCCGCGTTAATTTGGTCAATCATTATATTTCGCCGCTGGCCAGTGACACTCGAAGAACAGACGCAATAGTTGAAAGGCCCATTCGGGCCTTTTTTAATGGTACGTTTTAATGATTTCCAGGATGCCGTTAATAATAAACTGCACACCCATACATACCAGCAGGAATCCCATCAGACGGGAGATCGCTTCAATGC  >  NC_000913/1618092‑1618343

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