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%	110.7 / NA	35	*397* (TAG→TAA)	ydeA	arabinose efflux transporter, arabinose‑inducible
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (16/19);  total (16/19)

AGTTTGCACTGGTCAATGTCGATGATTGGTTATGTGGGCGCGGTGCCTGCTTTTGCCGCGTTAATTTGGTCAATCATTATATTTCGCCGCTGGCCAGTGACACTCGAAGAACAGACGCAATAGTTGAAAGGCCCATTCGGGCCTTTTTTAATGGTACGTTTTAATGATTTCCAGGATGCCGTTAATAATAAACTGCACACCCATACATACCAGCAGGAATCCCATCAGACGGGAGATCGCTTCAATGCCA  >  NC_000913/1618096‑1618345                                                                                                                           |                                                                                                                                aGTTTGCACTGGTCAATGTCGATGATTGGTTATGTGGGCGCGGTGCCTGCTTTTGCCGCGTTAATTTGGTCAATCATTATATTTCGCCGCTGGCCAGTGACACTCGAAGAACAGACGCAATAATTGAAAGGCCCATTCg                                                                                                                 >  1:505800/1‑139 (MQ=255)   tttGCACTGGTCAATGTCGATGATTGGTTATGTGGGCGCGGTGCCTGCTTTTGCCGCGTTAATTTGGTCAATCATTATATTTCGCCGCTGGCCAGTGACACTCGAAGAACAGACGCAATAATTGAAAGGCCCATTCggg                                                                                                               <  1:1020898/139‑1 (MQ=255)     tGCACTGGTAAATGTCGATGATTGGTTATGTGGGCGCGGTGCCTGCTTTTGCCGCGTTAATTTGGTCAATCATTATATTTCGCCGCTGGCCAGTGACACTCGAAGAACAGACGCAATAATTGAAAGGCCCATTCGGGcc                                                                                                             <  2:505800/139‑1 (MQ=255)             tCAATGTCGATGATTGGTTATGTGGGCGCGGTGCCTGCTTTTGCCGCGTTAATTTGGTCAATCATTATATTTCGCCGCTGGCCAGTGACACTCGAAGAACAGACGCAATAATTGAAAGGCCCATTCGGGCCTTTTTTaa                                                                                                     <  2:565188/139‑1 (MQ=255)                      atgatTGGTTATGTGGGCGCGGTGCCTGCTTTTGCCGCGTTAATTTGGTCAATCATTATATTTCGCCGCTGGCCAGTGACACTCGAAGAACAGACGCAATAATTGAAAGGCCCATTCgg                                                                                                                >  1:859677/1‑119 (MQ=255)                      atgatTGGTTATGTGGGCGCGGTGCCTGCTTTTGCCGCGTTAATTTGGTCAATCATTATATTTCGCCGCTGGCCAGTGACACTCGAAGAACAGACGCAATAATTGAAAGGCCCATTCgg                                                                                                                <  2:859677/119‑1 (MQ=255)                                     ggCGCGGTGCCTGCTTTTGCCGCGTTAATTTGGTCAATCATTATATTTCGCCGCTGGCCAGTGACACTCGAAGAACAGACGCAATAATTGAAAGGCCCATTCGGGCCTTTTTTAATGGTACGTTTTAATGATTTCCAgg                                                                             >  2:1057949/1‑139 (MQ=255)                                       cgcgGTGCCTGCTTTTGCCGCGTTAATTTGGTCAATCATTATATTTCGCCGCTGGCCAGTGACACTCGTAGAACAGACGCAATAATTGAAAGGCCCAt                                                                                                                    <  2:620342/98‑1 (MQ=255)                                       cgcgGTGCCTGCTTTTGCCGCGTTAATTTGGTCAATCATTATATTTCGCCGCTGGCCAGTGACACTCGAAGAACAGACGCAATAATTGAAAGGCCCAt                                                                                                                    >  1:620342/1‑98 (MQ=255)                                                         cgcgTTAATTTGGTCAATCATTATATTTCGCCGCTGGCCAGTGACACTCGAAGAACAGACGCAATAATTGAAAGGCCCATTCGGGCCTTTTTTAATGGTACGTTTTAATGATTTCCAGGATGCCGTTAATAATAAACTg                                                         <  1:1057949/139‑1 (MQ=255)                                                         cgcgTTAATTTGGTCAATCATTATATTTCGCCGCTGGCCAGTGACACTCGAAGAACAGACGCAATAATTGAAAGGCCCATTCGGGCCTTTTTTAATGGTACGTTTTAATGATTTCCAGGATGCCGTTAATAATAAACTg                                                         <  1:454206/139‑1 (MQ=255)                                                                     gTCAATCATTATATTTCGCCGCTGGCCAGTGACACTCGAAGAACAGACGCAATAATTGAAAGGCCCATTCGGGCCTTTTTTAATGGTACGTTTTAATGATTTCCAGGATGCCGTTAATAATAAACTGCACACCcataca                                             >  1:245393/1‑139 (MQ=255)                                                                                    tCGCCGCTGGCCAGTGACACTCGAAGAACAGACGCAATAATTGAAAGGCCCATTCGGGCCTTTTTTAATGGTACGTTTTAATGATTTCCAGGATGCCGTTAATAATAAACTGCACACCCATACATACCAGCAGGAATcc                              >  1:674969/1‑139 (MQ=255)                                                                                    tCGCCGCTGGCCAGTGACACTCGAAGAACAGACGCAATAATTGAAAGGCCCATTCGGGCCTTTTTTAATGGTACGTTTTAATGATTTCCAGGATGCCGTTAATAATAAACTGCACACCCATACATACCAGCAGGAATcc                              <  2:674969/139‑1 (MQ=255)                                                                                     cgccgcTGGCCAGTGACACTCGAAGAACAGACGCAATAATTGAAAGGCCCATTCGGGCCTTTTTTAATGGTa                                                                                                <  2:346836/72‑1 (MQ=255)                                                                                     cgccgcTGGCCAGTGACACTCGAAGAACAGACGCAATAATTGAAAGGCCCATTCGGGCCTTTTTTAATGGTa                                                                                                >  1:346836/1‑72 (MQ=255)                                                                                     cgccgcTGGCCAGTGACACTCGAAGAACAGACGCAATAATTGAAAGGCCCATTCGGGCCTTTTTTAATGGTACGTTTTAATGATTTCCAgg                                                                             >  2:525718/1‑91 (MQ=255)                                                                                     cgccgcTGGCCAGTGACACTCGAAGAACAGACGCAATAATTGAAAGGCCCATTCGGGCCTTTTTTAATGGTACGTTTTAATGATTTCCAgg                                                                             <  1:525718/91‑1 (MQ=255)                                                                                     cgccgcTGGCCAGTGACACTCGAAGAACAGACGCAATAATTGAAAGGCCCATTCGGGCCTTTTTTAATGGTACGTTTTAATGATTTCCAGGATGCCGTTAATAATAAACTGCACACCCATACATACCAGCAGGAATccc                             >  2:777324/1‑139 (MQ=255)                                                                                       ccgcTGGCCAGTGACACTCGAAGAACAGACGCAATAATTGAAAGGCCCATTCGGGCCTTTTTTAATGGTACGTTTTAATGATTTCCAGGATGCCGTTAATAATAAACTGCACACCCATACATACCAGCAGGAATCCCAt                           >  2:569747/1‑139 (MQ=255)                                                                                            ggCCAGTGACACTCGAAGAACAGACGCAATAATTGAAAGGCCCATTCGGGCCTTTTTTAATGGTACGTTTTAATGATTTCCAGGATGCCGTTAATAATAAACTGCACACCCATACATACCAGCAGGAATCCCATCAGAc                      <  1:777324/139‑1 (MQ=255)                                                                                                    acacTCGAAGAACAGACGCAATAATTGAAAGGCCCATTCGGGCCTTTTTTAATGGTACGTTTTAATGATTTCCAGGATGCCGTTAATAATAAACTGca                                                       >  1:884997/1‑98 (MQ=255)                                                                                                    acacTCGAAGAACAGACGCAATAATTGAAAGGCCCATTCGGGCCTTTTTTAATGGTACGTTTTAATGATTTCCAGGATGCCGTTAATAATAAACTGca                                                       <  2:884997/98‑1 (MQ=255)                                                                                                       cTCGAAGAACAGACGCAATAATTGAAAGGCCCATTCGGGCCTTTTTTAATGGTACGTTTTAATGATTTCCAGGATGCCGTTAATAATAAACTGca                                                       >  2:169675/1‑95 (MQ=255)                                                                                                       cTCGAAGAACAGACGCAATAATTGAAAGGCCCATTCGGGCCTTTTTTAATGGTACGTTTTAATGATTTCCAGGATGCCGTTAATAATAAACTGca                                                       <  1:169675/95‑1 (MQ=255)                                                                                                       cTCGAAGAACAGACGCAATAATTGAAAGGCCCATTCGGGCCTTTTTTAATGGTACGTTTTAATGATTTCCAGGATGCCGTTAATAATAAACTGCacac                                                    <  2:915959/98‑1 (MQ=255)                                                                                                       cTCGAAGAACAGACGCAATAATTGAAAGGCCCATTCGGGCCTTTTTTAATGGTACGTTTTAATGATTTCCAGGATGCCGTTAATAATAAACTGCacac                                                    >  1:915959/1‑98 (MQ=255)                                                                                                        tCGAAGAACAGACGCAATAATTGAAAGGCCCATTCGGGCCTTTTTTAATGGTACGTTTTAATGATTTCCAGGATGCCGTTAATAATAAACTGCACACCCATACATACCAGCAGGAATCCCATCAGACGGGAGATCGCtt          <  1:569747/139‑1 (MQ=255)                                                                                                          gaagaaCAGACGCAATAATTGAAAGGCCCATTCGGGCCTTTTTTAATGGTACGTTTTAATGATTTCCAGGATGCCGTTAATAATAAACTGCACACCca                                                 <  2:468267/98‑1 (MQ=255)                                                                                                          gaagaaCAGACGCAATAATTGAAAGGCCCATTCGGGCCTTTTTTAATGGTACGTTTTAATGATTTCCAGGATGCCGTTAATAATAAACTGCACACCca                                                 >  1:468267/1‑98 (MQ=255)                                                                                                            agaaCAGACGCAATAATTGAAAGGCCCATTCGGGCCTTTTTTAATGGTACGTTTTAATGATTTCCAGGATGCCGTTAATAATAAACTGca                                                       >  1:760272/1‑90 (MQ=255)                                                                                                            agaaCAGACGCAATAATTGAAAGGCCCATTCGGGCCTTTTTTAATAGTACGTTTTAATGATTTCCAGGATGCCGTTAATAATAAACTGca                                                       <  2:760272/90‑1 (MQ=255)                                                                                                                cAGACGCAATAATTGAAAGGCCCATTCGGGCCTTTTTTAATGGTACGTTTTAATGATTTCCAGGATGCCGTTAATAATAAACTGCACACCCATACATACCAGCAGGAATCCCATCAGACGGGAGATCGCTTCAATGCCa  <  2:245393/139‑1 (MQ=255)                                                                                                                    cGCAATAATTGAAAGGCCCATTCGGGCCTTTTTTAATGGTACGTTTTAATGATTTCCAGGATGCCGt                                                                      <  2:15793/67‑1 (MQ=255)                                                                                                                    cGCAATAATTGAAAGGCCCATTCGGGCCTTTTTTAATGGTACGTTTTAATGATTTCCAGGATGCCGt                                                                      >  1:15793/1‑67 (MQ=255)                                                                                                                           |                                                                                                                                AGTTTGCACTGGTCAATGTCGATGATTGGTTATGTGGGCGCGGTGCCTGCTTTTGCCGCGTTAATTTGGTCAATCATTATATTTCGCCGCTGGCCAGTGACACTCGAAGAACAGACGCAATAGTTGAAAGGCCCATTCGGGCCTTTTTTAATGGTACGTTTTAATGATTTCCAGGATGCCGTTAATAATAAACTGCACACCCATACATACCAGCAGGAATCCCATCAGACGGGAGATCGCTTCAATGCCA  >  NC_000913/1618096‑1618345

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