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%	117.1 / NA	38	*395* (TAG→TAA)	setC	putative arabinose efflux transporter
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (18/20);  total (18/20)

CACGCGTCGGCTGGATTATCGCCGGCTCTGTTGACGGAATTATGGTTGAAATCTGGAGCTACCATGCGTTGTTCTGGCTGGCGATAGGGATGTTGGGTATTGCGATGATTTGCCTGCTGTTTATTAAAGATATTTAGCTCACCAGTCTGGGCTAAGAAAAATCGGCGTGGCGCAAACAGAGCCATCCTCGCCATTTTTGCTATCATGCCTGCATACATAAACGACAAAACAGTATGCAGA  >  NC_000913/3838001‑3838240                                                                                                                                         |                                                                                                        cACGCGTCGGCTGGATTATCGCCGGCTCTGTTGACGGAATTATGGTTGAAATCTGGAGCTACCATGCGTTGTTCTTGCTGGCGATAGGGATGTTGGGTATTGCGATGATTTGCCTGCTGTTTATTAAAGATATTTAACt                                                                                                       >  1:505803/1‑139 (MQ=255)   cgcgTCGGCTGGATTATCGCCGGCTCTGTTGACGGAATTATGGTTGAAATCTGGAGCTACCATGCGTTGTTCTGGCTGGCGATAGGGATGTTGGGTATTGCGATGATTTGCCTGCTGTTTATTAAAGATATTTAACTCa                                                                                                     <  2:533554/139‑1 (MQ=255)   cgcgTCGGCTGGATTATCGCCGGCTCTGTTGACGGAATTATGGTTGAAATCTGGAGCTACCATGCGTTGTTCTGGCTGGCGATAGGGATGTTGGGTATTGCGATGATTTGCCTGCTGTTTATTAAAGATATTTAACTCa                                                                                                     <  1:679696/139‑1 (MQ=255)                            ctGTTGACGGAATTATGGTTGAAATCTGGAGCTACCATGCGTTGTTCTGGCTGGCGATAGGGATGTTGGGTATTGCGATGATTTGCCTGCTGTTTATTAAAGATATTTAACTCACCAGTCTGGGCTAAGAAAAATCGGc                                                                            <  2:518915/139‑1 (MQ=255)                             tGTTGACGGAATTATGGTTGAAATCTGGAGCTACCATGCGTTGTTCTGGCTGGCGATAGGGATGTTGGGTATTGCGATGATTTGCCTGCTGTTTATTAAAGATATTTAACTCACCAGTCTGGGCTAAGAAAAATCGGCg                                                                           >  1:177027/1‑139 (MQ=255)                              gTTGACGGAATTATGGTTGAAATCTGGAGCTACCATGCGTTGTTCTGGCTGGCGATAGGGATGTTGGGTATTGCGATGATTTGCCTGCTGTTTATTAAAGATATTTAACTCACCAGTCTGGGCTAAGAAAAATCGGCGt                                                                          >  2:639116/1‑139 (MQ=255)                                                            tACCATGCGTTGTTCTGGCTGGCGATAGGGATGTTGGGTATTGCGATGATTTGCCTGCTGTTTATTAAAGATATTTAACTCACCAGTCTGGGCTAAGAAAAATCGGCGTGGCGCAAACAGAGCCATCCTCGCCAttttt                                            <  1:472298/139‑1 (MQ=255)                                                              ccATGCGTTGTTCTGGCTGGCGATAGGGATGTTGGGTATTGCGATGATTTGCCTGCTGTTTATTAAAGATATTTAACTCACCAGTCTGGGCTAAGAAAAATCGGCGTGGCGCAAACAGAGCCATCCTCGCCATTTTTGc                                          >  2:244802/1‑139 (MQ=255)                                                                     ttgttCTGGCTGGCGATAGGGATGTTGGGTATTGCGATGATTTGCCTGCTGTTTATTAAAGATATTTAACTCACCAGTCTGGGCTAAGAAAAATCGGCGTGGCGCAAACAGAGCCATCCTCGCCATTTTTGCTATCATg                                   <  1:690546/139‑1 (MQ=255)                                                                      tgttCTGGCTGGCGATAGGGATGTTGGGTATTGCGATGATTTGCCTGCTGTTTATTAAAGATATTTAACTCACCAGTATGGGCTAAGAAAAATCGGCGTGGCGCAAACAGAGCCATCCTCGCCATTTTTGCTATCATGc                                  >  2:188944/1‑139 (MQ=255)                                                                        ttCTGGCTGGCGATAGGGATGTTGGGTATTGCGATGATTTGCCTGCTGTTTATTAAAGATATTTAACTCACCAGTCTGGGCTAAGAAAAATCGGCGTGGCGCAAACAGAGCCATCCTCGCCATTTTTGCTATCATGCCt                                <  1:189493/139‑1 (MQ=255)                                                                                 gcGATAGGGATGTTGGGGATTGCGATGATTTGCCTGCTGTTTATTAAAGATATTTAACTCACCAGTCTGGGCTAAGAAAAATCGGCGTGGCGCAAACAGAGCCATCCTCGCCATTTTTGCTATCATGCCTGcatacata                       <  2:76951/139‑1 (MQ=255)                                                                                  cGATAGGGATGTTGGGTATTGCGATGATTTGCCTGCTGTTTATTAAAGATATTTAACTCACCAGTCTGGGCTAAGAAAAATCGGCGTGGCGCAAACAGAGCCATCCTCGCCATTTTTGCTATCATGCCTGCATACATaa                      >  1:641098/1‑139 (MQ=255)                                                                                  cGATAGGGATGTTGGGTATTGCGATGATTTGCCTGCTGTTTATTAAAGATATTTAACTCACCAGTCTGGGCTAAGAAAAATCGGCGTGGCGCAAACAGAGCCATCCTCGCCATTTTTGCTATCATGCCTGCATACATaa                      >  1:718999/1‑139 (MQ=255)                                                                                      aGGGATGTTGGGTATTGCGATGATTTGCCTGCTGTTTATTAAAGATATTTAACTCACCAGTCTGGGCTAAGAAAAATCGGCGTGGCGCAAACAGAGCCATCCTCGCCATTTTTGCTATCATGCCTGCATACATAAACGa                  <  1:639116/139‑1 (MQ=255)                                                                                        ggATGTTGGGTATTGCGATGATTTGCCTGCTGTTTATTAAAGATATTTAACTCACCAGTCTGGGCTAAGAAAAATCGGCGTGGCGCAAACAGAGCCATCCTCGCCATTTTTGCTATCATGCCTGCATACATAAACGACa                <  2:177027/139‑1 (MQ=255)                                                                                          aTGTTGGGTATTGCGATGATTTGCCTGCTGTTTATTAAAGATATTTAACTCACCAGTCTGGGCTAAGAAAAATCGGCGTGGCGCAAACAGAGCCATCCTCGCCATTTTTGCTATCATGCCTGCATACATAAACGACaaa              <  1:244802/139‑1 (MQ=255)                                                                                             ttGGGTATTGCGATGATTTGCCTGCTGTTTATTAAAGATATTTAACTCACCAGTCTGGGCTAAGAAAAATCGGCGTGGCGCAAACAGAGCCATCCTCGCCATTTTTGCTATCATGCCTGCATACATAAACGACAAAACa           <  2:505803/139‑1 (MQ=255)                                                                                                  tATTGCGATGATTTGCCTGCTGTTTATTAAAGATATTTAACTCACCAGTCTGGGCTAAGAAAAATCGGCGTGGCGCAAACAGAGCCATCCTCGCCATTTTTGCTATCATGCCTGCATACATAAACGACAAAACAGTATg      >  1:783427/1‑139 (MQ=255)                                                                                                     tGCGATGATTTGCCTGCTGTTTATTAAAGATATTTAACTCACCAGTCTGGGCTAAGAAAAATCGGCGTGGCGCAAACAGAGCCATCCTCGCCATTTTTGCTATCATGCCTGCATACATAAACGa                  <  1:66352/124‑1 (MQ=255)                                                                                                     tGCGATGATTTGCCTGCTGTTTATTAAAGATATTTAACTCACCAGTCTGGGCTAAGAAAAATCGGCGTGGCGCAAACAGAGCCATCCTCGCCATTTTTGCTATCATGCCTGCATACATAAACGa                  >  2:66352/1‑124 (MQ=255)                                                                                                      gCGATGATTTGCCTGCTGTTTATTAAAGATATTTAACTCACCAGTCTGGGCTAAGAAAAATCGGCGTGGTGCAAACAGAGCCATCCTCGCCATTTTTGCTATCATGCCTGCATACATAAACGACAAAACAGTATGCaga  >  2:178479/1‑139 (MQ=255)                                                                                                       cGATGATTTGCCTGCTGTTTCTTAAAGATATTTAACTCACCAGTCTGGGCTAAGAAAAATCGGc                                                                            >  1:620561/1‑64 (MQ=255)                                                                                                       cGATGATTTGCCTGCTGTTTATTAAAGATATTTAACTCACCAGTCTGGGCTAAGAAAAATCGGc                                                                            <  2:620561/64‑1 (MQ=255)                                                                                                       cGATGATTTGCCTGCTGTTTATTAAAGATATTTAACTCACCAGTCTGGGCTAAGAAAAATCGGCGTGGCGCAAACAGAGCCATCCTCGCCATTTTTGCTATCa                                     >  1:7849/1‑103 (MQ=255)                                                                                                       cGATGATTTGCCTGCTGTTTATTAAAGATATTTAACTCACCAGTCTGGGCTAAGAAAAATCGGCGTGGCGCAAACAGAGCCATCCTCGCCATTTTTGCTATCa                                     <  2:7849/103‑1 (MQ=255)                                                                                                                      tgTTTATTAAAGATATTTAACTCACCAGTCTGGGCTAAGAAAAATCGGCGTGGCGCAAACAGAGCCATCCTCGCCATTTTTGCTATCATGCCTGCATACATAAACGACaaa              >  2:562990/1‑111 (MQ=255)                                                                                                                      tgTTTATTAAAGATATTTAACTCACCAGTCTGGGCTAAGAAAAATCGGCGTGGCGCAAACAGAGCCATCCTCGCCATTTTTGCTATCATGCCTGCATACATAAACGACaaa              <  1:562990/111‑1 (MQ=255)                                                                                                                                  atatTTAACTCACCAGTCTGGGCTAAGAAAAATCGGCGTGGCGCAAACAGAGCCATCCTCGCCAttttt                                            <  2:273626/69‑1 (MQ=255)                                                                                                                                  atatTTAACTCACCAGTCTGGGCTAAGAAAAATCGGCGTGGCGCAAACAGAGCCATCCTCGCCAttttt                                            >  1:273626/1‑69 (MQ=255)                                                                                                                                  atatTTAACTCACCAGTCTGGGCTAAGAAAAATCGGCGTGGCGCAAACAGAGCCATCCTCGCCATTTTTGCTATCATGCCTGcataca                         >  2:181571/1‑88 (MQ=255)                                                                                                                                  atatTTAACTCACCAGTCTGGGCTAAGAAAAATCGGCGTGGCGCAAACAGAGCCATCCTCGCCATTTTTGCTATCATGCCTGcataca                         <  1:181571/88‑1 (MQ=255)                                                                                                                                  atatTTAACTCACCAGTCTGGGCTAAGAAAAATCGGCGTGGCGCAAACAGAGCCATCCTCGCCATTTTTGCTATCATGCCTGCATACATAAACGACa                <  2:228933/97‑1 (MQ=255)                                                                                                                                  atatTTAACTCACCAGTCTGGGCTAAGAAAAATCGGCGTGGCGCAAACAGAGCCATCCTCGCCATTTTTGCTATCATGCCTGCATACATAAACGACa                >  1:228933/1‑97 (MQ=255)                                                                                                                                    atTTAACTCACCAGTCTGGGCTAAGAAAAATCGGCGTGGCGCAAACTGAGCCATCCTCGCCATTTTTGCTATCATGCCTGCATACATCAACGACAacac            <  2:485410/99‑1 (MQ=255)                                                                                                                                    atTTAACTCACCAGTCTGGGCTAAGAAAAATCGGCGTGGCGCAAACAGAGCCATCCTCGCCATTTTTGCTATCATGCCTGcatacata                       >  1:595541/1‑88 (MQ=255)                                                                                                                                    atTTAACTCACCAGTCTGGGCTAAGAAAAATCGGCGTGGCGCAAACAGAGCCATCCTCGCCATTTTTGCTATCATGCCTGcatacata                       <  2:595541/88‑1 (MQ=255)                                                                                                                                    atTTAACTCACCAGTCTGGGCTAAGAAAAATCGGCGTGGCGCAAACAGAGCCATCCTCGCCATTTTTGCTATCATGCCTGCATACATAAACGACAAAAc            >  1:485410/1‑99 (MQ=255)                                                                                                                                         |                                                                                                        CACGCGTCGGCTGGATTATCGCCGGCTCTGTTGACGGAATTATGGTTGAAATCTGGAGCTACCATGCGTTGTTCTGGCTGGCGATAGGGATGTTGGGTATTGCGATGATTTGCCTGCTGTTTATTAAAGATATTTAGCTCACCAGTCTGGGCTAAGAAAAATCGGCGTGGCGCAAACAGAGCCATCCTCGCCATTTTTGCTATCATGCCTGCATACATAAACGACAAAACAGTATGCAGA  >  NC_000913/3838001‑3838240

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