Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,838,137	G→A	*395* (TAG→TAA)	setC →	putative arabinose efflux transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,838,137	0	G	A	100.0%	107.1 / NA	35	*395* (TAG→TAA)	setC	putative arabinose efflux transporter
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (18/17);  total (18/17)

GCGTCGGCTGGATTATCGCCGGCTCTGTTGACGGAATTATGGTTGAAATCTGGAGCTACCATGCGTTGTTCTGGCTGGCGATAGGGATGTTGGGTATTGCGATGATTTGCCTGCTGTTTATTAAAGATATTTAGCTCACCAGTCTGGGCTAAGAAAAATCGGCGTGGCGCAAACAGAGCCATCCTCGCCATTTTTGCTATCATGCCTGCATACATAAACGACAAAACAGTATGCAGAGGGAAAAATGGGTTCCACCAGAAAGG  >  NC_000913/3838004‑3838266                                                                                                                                      |                                                                                                                                  gcgTCGGCTGGATTATCGCCGGCTCTGTTTACGGAATTATGGTTGAAATCTGGAGCTACCATGCGTTGTTCTGGCTGGCGATAGGGCTGTTGGGTATTGCGATGATTTGCCTGCTGTTTATTAAAGATATTTAACTCAc                                                                                                                              <  1:629461/139‑1 (MQ=255)  cgTCGGCTGGATTATCGCCGGCTCTGTTGACGGAATTATGGTTGAAATCTGGAGCTACCATGCGTTGCTCTGGCTGGCGATAGGGATGTTGGGTATTGCGATGATTTGCCTGCTGTTTATTAAAGATATTTAACTCAcc                                                                                                                             >  2:197908/1‑139 (MQ=255)      ggCTGGATTATCGCCGGCTCTGTTGACGGAATTATGGTTGAAATCTGGAGCTACCATGCGTTGTTCTGGCTGGCGATAGGGATGTTGGGTATTGCGATGATTTGCCTGCTGTTTATTAAAGATATTTAACTCACCAGTc                                                                                                                         <  1:698383/139‑1 (MQ=255)          ggATTATCGCCGGCTCTGTTGACGGAATTATGGTTGAAATCTGGAGCTACCATGCGTTGTTCTGGCTGGCGATAGGGATGTTGGGTATTGCGATGATTTGCCTGCTGTTTATTAAAGATATTTAACTCACCAGTCTggg                                                                                                                     <  1:61532/139‑1 (MQ=255)                  gCCGGCTCTGTTGACGGAATTATGGTTGAAATCTGGAGCTACCATGCGTTGTTCTGGCTGGCGATAGGGATGTTGGGTATTGCGATGATTTGCCTGCTGTTTATTAAAGATATTTAACTCACCAGTCTGGGCTAAGaaa                                                                                                             >  2:900862/1‑139 (MQ=255)                         ctGTTGACGGAATTATGGTTGAAATCTGGAGCTACCATGCGTTGTTCTGGCTGGCGATAGGGATGTTGGGTATTGCGATGATTTGCCTGCTGTTTATTAAAGATATTTAACTCACCAGTCTGGGCTAAGAAAAATTGGc                                                                                                      <  2:638134/139‑1 (MQ=255)                                   aaTTATGGTTGAAATCTGGAGCTACCATGCGTTGTTCTGGCTGGCGATAGGGATGTTGGGTATTGCGATGATTTGCCTGCTGTTTATTAAAGATATTTAACTCACCAGTCTGGGCTAAGAAAAATCGGCGTGGCGCaaa                                                                                            <  2:861594/139‑1 (MQ=255)                                             gAAATCTGGAGCTACCATGCGTTGTTCTGGCTGGCGATAGGGATGTTGGGTATTGCGATGATTTGCCTGCTGTTTATTAAAGATATTTAACTCACCAGTCTGGGCTAAGAAAAATCGGCGTGGCGCAAACAGAGCCATc                                                                                  >  2:915238/1‑139 (MQ=255)                                             gAAATCTGGAGCTACCATGCGTTGTTCTGGCTGGCGATAGGGATGTTGGGTATTGCGATGATTTGCCTGCTGTTTATTAAAGATATTTAACTCACCAGTCTGGGCTAAGAAAAATCGGCGTGGCGCAAACAGAGCCATc                                                                                  >  1:736389/1‑139 (MQ=255)                                                                cGTTGTTCTGGCTGGCGATAGGGATGTTGGGTATTGCGATGATTTGCCTGCTGTTTATTAAAGATATTTAACTCACCAGTCTGGGCTAAGAAAAATCGGCGTGGCGCAAACAGAGCCATCCTCGCCATTTTTGCTATCa                                                               <  1:184994/139‑1 (MQ=255)                                                                     ttCTGGCTGGCGATAGGGATGTTGGGTATTGCGATGATTTGCCTGCTGTTTATTAAAGATATTTAACTCACCAGTCTGGGCTAAGAAAAATCGGCGTGGCGCAAACAGAGCCATCCTCGCCATTTTTGCTATCATGCCt                                                          >  1:630010/1‑139 (MQ=255)                                                                     ttCTGGCTGGCGATAGGGATGTTGGGTATTGCGATGATTTGCCTGCTGTTTATTAAAGATATTTAACTCACCAGTCTGGGCTAAGAAAAATCGGCGTGGCGCAAACAGAGCCATCCTCGCCATTTTTGCTATCATGCCt                                                          >  2:207803/1‑139 (MQ=255)                                                                      tCTGGCTGGCGATAGGGATGTTGGGTATTGCGATGATTTGCCTGCTGTTTATTAAAGATATTTAACTCACCAGTCTGGGCTAAGAAAAATCGGCGTGGCGCAAACAGAGCCATCCTCGCCATTTTTGCTATCATGCCTg                                                         <  2:630010/139‑1 (MQ=255)                                                                               cGATAGGGATGTTGGGTATTGCGATGATTTGCCTGCTGTTTATTAAAGATATTTAACTCACCAGTCTGGGCTAAGAAAAATCGGCGTGGTGCAAACAGAGCCATCCTCGCCATTTTTGCTATCATGCCTGCATACATaa                                                >  2:832390/1‑139 (MQ=255)                                                                               cGATAGGGATGTTGGGTATTGCGATGATTTGCCTGCTGTTTATTAAAGATATTTAACTCACCAGTCTGGGCTAAGAAAAATCGGCGTGGCGCAAACAGAGCCATCCTCGCCATTTTTGCTATCATGCCTGCATACATaa                                                >  2:131463/1‑139 (MQ=255)                                                                                gATAGGGATGTTGGGTATTGCGATGATTTGCCTGCTGTTTATTAAAGATATTTAACTCACCAGTCTGGGCTAAGAAAAATCGGCGTGGCGCAAACAGAGCCATCCTCGCCATTTTTGCTATCATGCCTGCATACATaaa                                               <  1:131463/139‑1 (MQ=255)                                                                                       aTGTTGGGTATTGCGATGATTTGCCTGCTGTTTATTAAAGATATTTAACTCACCAGTCTGGGCTAAGAAAAATCGGCGTGGTGCAAACAGAGCCATCCTCGCCATTTTTGCTATCATGCCTGCATACATAAACGACaaa                                        <  1:832390/139‑1 (MQ=255)                                                                                       aTGTTGGGTATTGCGATGATTTGCCTGCTGTTTATTAAAGATATTTAACTCACCAGTCTGGGCTAAGAAAAATCGGCGTGGCGCAAACAGAGCCATCCTCGCCATTTTTGCTATCATGCCTGCATACATAAACGACaaa                                        <  1:207803/139‑1 (MQ=255)                                                                                        tGTTGGGTATTGCGATGATTTGCCTGGTGTTTATTAAAGATATTTAACTCACCAGTCTGGGCTAAGAAAACTCTGCGGGGCGCAAACAGAGCCATCCTCGCCATTGTTGCTATCATGCCTGCATACATAAACGACaaaa                                       >  2:232591/1‑139 (MQ=255)                                                                                        tGTTGGGTATTGCGATGATTTGCCTGCTGTTTATTAAAGATATTTAACTCACCAGTCTGGGCTAAGAAAAATCGGCGTGGCGCAAACAGAGCCATCCTCGCCATTTTTGCTATCATGCCTGCATACATAAACGACaaaa                                       >  1:397577/1‑139 (MQ=255)                                                                                         gTTGGGTATTGCGATGATTTGCCTGCTGTTTATTAAAGATATTTAACTCAc                                                                                                                              <  1:572924/51‑1 (MQ=255)                                                                                         gTTGGGTATTGCGATGATTTGCCTGCTGTTTATTAAAGATATTTAACTCAc                                                                                                                              >  2:572924/1‑51 (MQ=255)                                                                                             ggTATTGCGATGATTTGCCTGCTGTTTATTAAAGATATTTAACTCACCAGTCTGGGCTAAGAAAAATCGGCGTGGCGCAAACAGAGCCATCCTCGCCATTTTTGCTATCATGCCTGCATACATAAACGa                                            >  2:330844/1‑129 (MQ=255)                                                                                             ggTATTGCGATGATTTGCCTGCTGTTTATTAAAGATATTTAACTCACCAGTCTGGGCTAAGAAAAATCGGCGTGGCGCAAACAGAGCCATCCTCGCCATTTTTGCTATCATGCCTGCATACATAAACGa                                            <  1:330844/129‑1 (MQ=255)                                                                                                  tGCGATGATTTGCCTGCTGTTTATTAAAGATATTTAACTCACCAGTCTGGGCTAAGAAAAATCGGCGTGGCGCAAACAGAGCCATCCTCGCCATTTTTGCTATCATGCCTGca                                                       <  1:797388/113‑1 (MQ=255)                                                                                                  tGCGATGATTTGCCTGCTGTTTATTAAAGATATTTAACTCACCAGTCTGGGCTAAGAAAAATCGGCGTGGCGCAAACAGAGCCATCCTCGCCATTTTTGCTATCATGCCTGca                                                       >  2:797388/1‑113 (MQ=255)                                                                                                    cGATGATTTGCCTGCTGTTTATTAAAGATATTTAACTCACCAGTCTGGGCTAAGAAAAATCGGCGTGGCGCAAACAGAGCCATCCTCGCCATTTTTGCTATCACGCCTGCATACATAAACGACACAACAGAATGCagag                           >  2:246456/1‑139 (MQ=255)                                                                                                            tGCCTGCTGTTTATTAAAGATATTTAACTCACCAGTCTGGGCTAAGAAAAATCGGCGTGGCGCAAACAGAGCCATCCTCGCCATTTTTGCTATCATGCCTg                                                         <  1:854679/101‑1 (MQ=255)                                                                                                            tGCCTGCTGTTTATTAAAGATATTTAACTCACCAGTCTGGGCTAAGAAAAATCGGCGTGGCGCAAACAGAGCCATCCTCGCCATTTTTGCTATCATGCCTg                                                         >  2:854679/1‑101 (MQ=255)                                                                                                                  ctgTTTATTAAAGATATTTAACTCACCAGTCTGGGCTAAGAAAAATCGGCGTGGCGCAAACAGAGCCATCCTCGCCATTTTTGCTATCATGCCTGCATACATAAACGACAAAACAGTATGCAGAGGGAAAAATGGGTTc             <  1:72265/139‑1 (MQ=255)                                                                                                                   tgTTTATTAAAGATATTTAACTCACCAGTCTGGGCTAAGAAAAATCGGCGTGGCGCAAACAGAGCCATCCTCGCCATTTTTGCTATCATGCCTGCATACATAAACGACAAAACAGTATGCAGAGGGAAAAATGGGTTcc            >  1:699554/1‑139 (MQ=255)                                                                                                                     tttATTAAAGATATTTAACTCACCAGTCTGGGCTAAGAAAAATCGGCGTGGCGCAAACAGAGCCATCCTCGCCATTTTTGCTATCATGCCTGCATACATAAACGACAAAACAGTATGCAGAGGGAAAAATGGGTTccac          <  2:699554/139‑1 (MQ=255)                                                                                                                            aaGATATTTAACTCACCAGTCTGGGCTAAGAAAAATCGGCGt                                                                                                    <  2:447362/42‑1 (MQ=255)                                                                                                                            aaGATATTTAACTCACCAGTCTGGGCTAAGAAAAATCGGCGt                                                                                                    >  1:447362/1‑42 (MQ=255)                                                                                                                             aGATATTTAACTCACCAGTCTGGGCTAAGAAAAATCGGCGTGGCGCAAACAGAGCCATCCTCGCCATTTTTGCTATCATGCCTGCATACATAAACGACAAAACAGTATGCAGAGGGAAAAATGGGTTCCACCAGAAAgg  >  1:187155/1‑139 (MQ=255)                                                                                                                                      |                                                                                                                                  GCGTCGGCTGGATTATCGCCGGCTCTGTTGACGGAATTATGGTTGAAATCTGGAGCTACCATGCGTTGTTCTGGCTGGCGATAGGGATGTTGGGTATTGCGATGATTTGCCTGCTGTTTATTAAAGATATTTAGCTCACCAGTCTGGGCTAAGAAAAATCGGCGTGGCGCAAACAGAGCCATCCTCGCCATTTTTGCTATCATGCCTGCATACATAAACGACAAAACAGTATGCAGAGGGAAAAATGGGTTCCACCAGAAAGG  >  NC_000913/3838004‑3838266

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: ‑