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

GGTAATCAGGTGAGTTTGCACTGGTCAATGTCGATGATTGGTTATGTGGGCGCGGTGCCTGCTTTTGCCGCGTTAATTTGGTCAATCATTATATTTCGCCGCTGGCCAGTGACACTCGAAGAACAGACGCAATAGTTGAAAGGCCCATTCGGGCCTTTTTTAATGGTACGTTTTAATGATTTCCAGGATGCCGTTAATAATAAACTGCACACCCATACA  >  NC_000913/1618084‑1618302                                                                                                                                       |                                                                                     ggtaATCAGGTGAGTTTGCACTGGTCAATGTCGATGATTGGTTATGTGGGCGCGGTGCCTGCTTTTGCCGCGTTAATTTGGTCAATCATTATATTTCGCCGCTGGCCAGTGACACTCGAAGAACAGACGCAATAATTGa                                                                                  <  1:1395342/139‑1 (MQ=255)         ggTGAGTTTGCACTGGTCAATGTCGATGATTGGTTATGTGGGCGCGGTGCCTGCTTTTGCCGCGTTAATTTGGTCAATCATTATATTTCGCCGCTGGCCAGTGACACTCGAAGAACAGACGCAATAATTGAAAGGCCCa                                                                          >  1:755787/1‑139 (MQ=255)                    aCTGGTCAATGTCGATGATTGGTTATGTGGGCGCGGTGCCTGCTTTTGCCGCGTTAATTTGGTCAATCATTATATTTCGCCGCTGGCCAGTGACACTCGAAGAACAGACGCAATAATTGAAAGGCCCATTCGGGCCttt                                                               <  1:961603/139‑1 (MQ=255)                     cTGGTCAATGTCGATGATTGGTTATGTGGGCGCGGTGCCTGCTTTTGCCGCGTTAATTTGGTCAATCATTATATTTCGCCGCTGGCCAGTGACACTCGAAGAACAGACGCAATAATTGAAAGGCCCATTCGGGCCtttt                                                              >  2:232513/1‑139 (MQ=255)                         tCAATGTCGATGATTGGTTATGTGGGCGCGGTGCCTGCTTTTGCCGCGTTAATTTGGTCAATCATTATATTTCGCCGCTGGCCAGTGACACTCGAAGAACAGACGCAATAATTGAAAGGCCCATTCGGGCCTTTTTTaa                                                          >  2:788577/1‑139 (MQ=255)                             tGTCGATGATTGGTTATGTGGGCGCGGTGCCTGCTTTTGCCGCGTTAATTTGGTCAATCATTATATTTCGCCGCTGGCCAGTGACACTCGAAGAACAGACGCAATAATTGAAAGGCCCATTCGGGcc                                                                  <  2:1513760/127‑1 (MQ=255)                             tGTCGATGATTGGTTATGTGGGCGCGGTGCCTGCTTTTGCCGCGTTAATTTGGTCAATCATTATATTTCGCCGCTGGCCAGTGACACTCGAAGAACAGACGCAATAATTGAAAGGCCCATTCGGGcc                                                                  >  1:1513760/1‑127 (MQ=255)                                  atgatTGGTTATGTGGGCGCGGTGCCTGCTTTTGCCGCGTTAATTTGGTCAATCATTATATTTCGCCGCTGGCCAGTGACACTCGAAGAACAGACGCAATAATTGAAAGGCCCATTCGGGCCTTTTTTAATGGTACGtt                                                 >  2:939969/1‑139 (MQ=255)                                       tGGTTATGTGGGCGCGGTGCCTGCTTTTGCCGCGTTAATTTGGTCAATCATTATATTTCGCCGCTGGCCAGTGACACTCGAAGAACAGACGCAATAATTGAAAGGCCCATTCGGGCCTTTTTTAATGGTACGTTTTAAt                                            <  2:1226226/139‑1 (MQ=255)                                        ggTTATGTGGGCGCGGTGCCTGCTTTTGCCGCGTTAATTTGGTCAATCATTATATTTCGCCGCTGGCCAGTGACACTCGAAGAACAGACGCAATAATTGAAAGGCCCATTCGGGCCTTTTTTAATGGTACGTTTTAATg                                           <  2:435639/139‑1 (MQ=255)                                             tgtgGGCGCGGTGCCTGCTTTTGCCGCGTTAATTTGGTCAATCATTATATTTCGCCGCTGGCCAGTGACACTCGAAGAACAGACGCAATAATTGAAAGGCCCATTCGGGCCTTTTTTAATGGTACGTTTTAATGATTTc                                      >  1:631226/1‑139 (MQ=255)                                                   cgcgGTGCCTGCTTTTGCCGCGTTAATTTGGTCAATCATTATATTTCGCCGCTGGCCAGTGACACTCGAAGAACAGACGCAATAATTGAAAGGCCCATTCGGGcc                                                                  >  2:880141/1‑105 (MQ=255)                                                   cgcgGTGCCTGCTTTTGCCGCGTTAATTTGGTCAATCATTATATTTCGCCGCTGGCCAGTGACACTCGAAGAACAGACGCAATAATTGAAAGGCCCATTCGGGcc                                                                  <  1:880141/105‑1 (MQ=255)                                                   cgcgGTGCCTGCTTTTGCCGCGTTAATTTGGTCAATCATTATATTTCGCCGCTGGCCAGTGACACTCGAAGAACAGACGCAATAATTGAAAGGCCCATTCGGGCCTTTTTTAATGGTACGTTTTAATGATTTCCAGGAt                                >  1:1222524/1‑139 (MQ=255)                                                   cgcgGTGCCTGCTTTTGCCGCGTTAATTTGGTCAATCATTATATTTCGCCGCTGGCCAGTGACACTCGAAGAACAGACGCAATAATTGAAAGGCCCATTCGGGCCTTTTTTAATGGTACGTTTTAATGATTTCCAGGAt                                >  1:222360/1‑139 (MQ=255)                                                   cgcgGTGCCTGCTTTTGCCGCGTTAATTTGGTCAATCATTATATCTCGCCGCTGGCCAGTGACACTCGAAGAACAGACGCAATAATTGAAAGGCCCATTCGGGCCTTTTTTAATGGTACGTTTTaa                                             >  1:1442375/1‑126 (MQ=255)                                                   cgcgGTGCCTGCTTTTGCCGCGTTAATTTGGTCAATCATTATATCTCGCCGCTGGCCAGTGACACTCGAAGAACAGACGCAATAATTGAAAGGCCCATTCGGGCCTTTTTTAATGGTACGTTTTaa                                             <  2:1442375/126‑1 (MQ=255)                                                    gcgGTGCCTGCTTTTGCCGCGTTAATTTGGTCAATCATTATATTTCGCCGCTGGCCAGTGACACTCGAAGAACAGACGCAATAATTGAAAGGCCCATTCGGGCCTTTTTTAATGGTACGTTTTAATGATTTCCAGGATg                               >  2:1000735/1‑139 (MQ=255)                                                                  tGCCGCGTTAATTTGGTCAATCATTATATTTCGCCGCTGGCCAGTGACACTCGAAGAACAGACGCAATAATTGAAAGGCCCATTCGGGCCTTTTTTAATGGTACGTTTTAATGATTTCCAGGATGCCGt                           <  1:425018/129‑1 (MQ=255)                                                                  tGCCGCGTTAATTTGGTCAATCATTATATTTCGCCGCTGGCCAGTGACACTCGAAGAACAGACGCAATAATTGAAAGGCCCATTCGGGCCTTTTTTAATGGTACGTTTTAATGATTTCCAGGATGCCGt                           >  2:425018/1‑129 (MQ=255)                                                                  tGCCGCGTTAATTTGGTCAATCATTATATTTCGCCGCTGGCCAGTGACACTCGAAGAACAGACGCAATAATTGAAAGGCCCATTCGGGCCTTTTTTAATGGTACGTTTTAATGATTTCCAGGATGCCGTTAATAATaaa                 <  1:177112/139‑1 (MQ=255)                                                                    ccGCGTTAATTTGGTCAATCATTATATTTCGCCGCTGGCCAGTGACACTCGAAGAACAGACGCAATAATTGAAAGGc                                                                             <  2:232109/77‑1 (MQ=255)                                                                    ccGCGTTAATTTGGTCAATCATTATATTTCGCCGCTGGCCAGTGACACTCGAAGAACAGACGCAATAATTGAAAGGc                                                                             >  1:232109/1‑77 (MQ=255)                                                                    ccGCGTTAATTTGGTCAATCATTATATTTCGCCGCTGGCCAGTGACACTCGAAGAACAGACGCAATAATTGAAAGGCCCATTCGGGCCTTTTTTAATGGTACGTTTTAATGATTTCCAGGATGCCGTTAATAATAAACt               >  2:839243/1‑139 (MQ=255)                                                                     cgcgTTAATTTGGTCAATCATTATATTTCGCCGCTGGCCAGTGACACTCGAAGAACAGACGCAATAATTGAAAGGCCCATTCGGGCCTTTTTTAATGGTACGTTTTAATGATTTCCAGGATGCCGTTAATAATAAACTg              <  2:755787/139‑1 (MQ=255)                                                                     cgcgTTAATTTGGTCAATCATTATATTTCGCCGCTGGCCAGTGACACTCGAAGAACAGACGCAATAATTGAAAGGCCCATTCGGGCCTTTTTTAATGGTACGTTTTAATGATTTCCAGGATGCCGTTAATAATAAACTg              <  1:1351311/139‑1 (MQ=255)                                                                       cgTTAATTTGGTCAATCATTATATTTCGCCGCTGGCCAGTGACACTCGAAGAACAGACGCAATAATTGAAAGGCCCATTCGGGCCTTTTTTAATGGTACGTTTTAATGATTTCCAgg                                  >  2:189175/1‑117 (MQ=255)                                                                       cgTTAATTTGGTCAATCATTATATTTCGCCGCTGGCCAGTGACACTCGAAGAACAGACGCAATAATTGAAAGGCCCATTCGGGCCTTTTTTAATGGTACGTTTTAATGATTTCCAgg                                  <  1:189175/117‑1 (MQ=255)                                                                       cgTTAATTTGGTCAATCATTATATTTCGCCGCTGGCCAGTGACACTCGAAGAACAGACGCAATAATTGAAAGGCCCATTCGGGCCTTTTTTAATGGTACGTTTTAATGATTTCCAGGATGCCGTTAATAATAAACTGca            <  2:943121/139‑1 (MQ=255)                                                                       cgTTAATTTGGTCAATCATTATATTTCGCCGCTGGCCAGTGACACTCGAAGAACAGACGCAATAATTGAAAGGCCCATTCGGGCCTTTTTTAATGGTACGTTTTAATGATTTCCAGGATGCCGTTAATAATAAACTGca            <  1:232513/139‑1 (MQ=255)                                                                       cgTTAATTTGGTCAATCATTATATTTCGCCGCTGGCCAGTGACACTCGAAGAACAGACGCAATAATTGAAAGGCCCATTCGGGCCTTTTTTAATGGTACGTTTTAATGATTTCCAGGATGCCGTTAATAATAAACTGca            <  1:939969/139‑1 (MQ=255)                                                                       cgTTAATTTGGTCAATCATTATATTTCGCCGCTGGCCAGTGACACTCGAAGAACAGACGCAATAATTGAAAGGCCCATTCGGGCCTTTTTTAATGGTACGTTTTAATGATTTCCAGGATGCCGTTAATAATAAACTGca            <  1:473352/139‑1 (MQ=255)                                                                       cgTTAATTTGGTCAATCATTATATTTCGCCGCTGGCCAGTGACACTCGAAGAACAGACGCAATAATTGAAAGGCCCATTCGGGCCTTTTTTAATGGTACGTTTTAATGATTTCCAGGATGCCGTTAATAATAAACTGca            <  1:839243/139‑1 (MQ=255)                                                                       cgTTAATTTGGTCAATCATTATATTTCGCCGCTGGCCAGTGACACTCGAAGAACAGACGCAATAATTGAAAGGCCCATTCGGGCCTTTTTTAATGGTACGTTTTAATGATTTCCAGGATGCCGTTAATAATAAACTGca            <  1:788577/139‑1 (MQ=255)                                                                          tAATTTGGTCAATCATTATATTTCGCCGCTGGCCAGTGACACTCGAAGAACAGACGCAATAATTGAAAGGCCCATTCGGGCCTTTTTTAATGGTACGTTTTAATGATTTCCAGGATGCCGTTAATAATAAACTGCacac         <  1:963408/139‑1 (MQ=255)                                                                          tAATTTGGTCAATCATTATATTTCGCCGCTGGCCAGTGACACTCGAAGAACAGACGCAATAATTGAAAGGCCCATTCGGGCCTTTTTTAATGGTACGTTTTAATGATTTCCAGGATGCCGTTAATAATAAACTGCacac         <  2:222360/139‑1 (MQ=255)                                                                                           atatTTCGCCGCTGGCCAGTGACACTCGAAGAACAGACGCAATAATTGAAAGGCCCATTCGGGCCTTTTTTAATGGTACGTTTTaa                                             >  2:1477653/1‑86 (MQ=255)                                                                                           atatTTCGCCGCTGGCCAGTGACACTCGAAGAACAGACGCAATAATTGAAAGGCCCATTCGGGCCTTTTTTAATGGTACGTTTTaa                                             <  1:1477653/86‑1 (MQ=255)                                                                                                tCGCCGCTGGCCAGTGACGCTCGAAGAACAGACGCAATAATTGAAAGGCCCATTCGGGcg                                                                  <  1:942560/60‑2 (MQ=255)                                                                                                tCGCCGCTGGCCAGTGACACTCGAAGAACAGACGCAATAATTGAAAGGCCCATTCGGGcc                                                                  >  2:1369127/1‑60 (MQ=255)                                                                                                tCGCCGCTGGCCAGTGACACTCGAAGAACAGACGCAATAATTGAAAGGCCCATTCGGGcc                                                                  >  2:942560/1‑60 (MQ=255)                                                                                                tCGCCGCTGGCCAGTGACACTCGAAGAACAGACGCAATAATTGAAAGGCCCATTCGGGcc                                                                  <  1:1369127/60‑1 (MQ=255)                                                                                                 cgccgcTGGCCAGTGACACTCGAAGAACAGACGCAATAATTGAAAGGCCCATTCGGGcc                                                                  <  2:1264297/59‑1 (MQ=255)                                                                                                 cgccgcTGGCCAGTGACACTCGAAGAACAGACGCAATAATTGAAAGGCCCATTCGGGcc                                                                  >  1:1264297/1‑59 (MQ=255)                                                                                                 cgccgcTGGCCAGTGACACTCGAAGAACAGACGCAATAATTGAAAGGCCCATTCGGGCCTTTTTTAATGGTACGTTTTAATGATTTCCAGGATGCCGt                           >  1:114915/1‑98 (MQ=255)                                                                                                 cgccgcTGGCCAGTGACACTCGAAGAACAGACGCAATAATTGAAAGGCCCATTCGGGCCTTTTTTAATGGTACGTTTTAATGATTTCCAGGATGCCGt                           <  2:114915/98‑1 (MQ=255)                                                                                                    cgcTGGCTAGTGACACTCGAAAAACAGACGCAATAATTGAAAGGCCCATTCGGGCCTTTTTTAATGGTACGTTTTAATGATTTCCAGGATGCCGt                           <  2:171781/95‑1 (MQ=255)                                                                                                    cgcTGGCCAGTGACACTCGAAGAACAGACGCAATAATTGAAAGGCCCATTCGGGCCTTTTTTAATGGTACGTTTTAATGATTTCCAGGATGCCGt                           >  1:171781/1‑95 (MQ=255)                                                                                                            aGTGACACTCGAAGAACAGACGCAATAATTGAAAGGCCCATTCggg                                                                    <  2:290266/46‑1 (MQ=255)                                                                                                            aGTGACACTCGAAGAACAGACGCAATAATTGAAAGGCCCATTCggg                                                                    >  1:290266/1‑46 (MQ=255)                                                                                                                   cTCGAAGAACAGACGCAATAATTGAAAGGCCCATTCGGGCCTTTTTTAATGGt                                                      <  1:1525964/53‑1 (MQ=255)                                                                                                                   cTCGAAGAACAGACGCAATAATTGAAAGGCCCATTCGGGCCTTTTTTAATGGt                                                      >  2:1525964/1‑53 (MQ=255)                                                                                                                   cTCGAAGAACAGACGCAATAATTGAAAGGCCCATTCGGGCCTTTTTTAATGGTACGTTTTAATGATTTCCAGGATGCCGTtaataata                   <  2:1347894/88‑1 (MQ=255)                                                                                                                   cTCGAAGAACAGACGCAATAATTGAAAGGCCCATTCGGGCCTTTTTTAATGGTACGTTTTAATGATTTCCAGGATGCCGTtaataata                   >  1:1347894/1‑88 (MQ=255)                                                                                                                              gACGCAATAATTGAAAGGCCCATTCGGGCCTTTTTTAATGGTACGTTTTAATGATTTCCAGGATGCCGTTAATAATAAACTGCACACCcataca  <  2:1123725/94‑1 (MQ=255)                                                                                                                              gACGCAATAATTGAAAGGCCCATTCGGGCCTTTTTTAATGGTACGTTTTAATGATTTCCAGGATGCCGTTAATAATAAACTGCACACCcataca  >  1:1123725/1‑94 (MQ=255)                                                                                                                               aCGCAATAATTGAAAGGCCCATTCGGGCCTTTTTTAATGGTACGTTTTAATg                                           <  2:530533/52‑1 (MQ=255)                                                                                                                               aCGCAATAATTGAAAGGCCCATTCGGGCCTTTTTTAATGGTACGTTTTAATg                                           >  1:530533/1‑52 (MQ=255)                                                                                                                                       |                                                                                     GGTAATCAGGTGAGTTTGCACTGGTCAATGTCGATGATTGGTTATGTGGGCGCGGTGCCTGCTTTTGCCGCGTTAATTTGGTCAATCATTATATTTCGCCGCTGGCCAGTGACACTCGAAGAACAGACGCAATAGTTGAAAGGCCCATTCGGGCCTTTTTTAATGGTACGTTTTAATGATTTCCAGGATGCCGTTAATAATAAACTGCACACCCATACA  >  NC_000913/1618084‑1618302

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