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

ACGCGTCGGCTGGATTATCGCCGGCTCTGTTGACGGAATTATGGTTGAAATCTGGAGCTACCATGCGTTGTTCTGGCTGGCGATAGGGATGTTGGGTATTGCGATGATTTGCCTGCTGTTTATTAAAGATATTTAGCTCACCAGTCTGGGCTAAGAAAAATCGGCGTGGCGCAAACAGAGCCATCCTCGCCATTTTTGCTATCATGCCTGCATACATAAACGACAAAACAGTATGCAGAGGGAAAAATGGGTTCCACCAG  >  NC_000913/3838002‑3838261                                                                                                                                        |                                                                                                                             aCGCGTCGGCTGGATTATCGCCGGCTCTGTTGACGGAATTATGGTTGAAATCTGGAGCTACCATGCGTTGTTCTGGCTGGCGATAGGGATGTT‑GGTATTGCGATGATTTGCCTGCTGTTTATTAAAGATATTTAACTCa                                                                                                                          <  2:1387646/139‑1 (MQ=255)  cgcgTCGGCTGGATTATCGCCGGCTCTGTTGACGGAATTATGGTTGAAATCTGGAGCTACCATGCGTTGTTCTGGCTGGCGATAGGGATGTTGGGTATTGCGATGATTTGCCTGCTGTTTATTAAAGATATTTAACTCa                                                                                                                          <  1:199997/139‑1 (MQ=255)  cgcgTCGGCTGGATTATCGCCGGCTCTGTTGACGGAATTATGGTTGAAATCTGGAGCTACCATGCGTTGTTCTGGCTGGCGATAGGGATGTTGGGTATTGCGATGATTTGCCTGCTGTTTATTAAAGATATTTAACTCa                                                                                                                          >  2:199997/1‑139 (MQ=255)                         ctctGTTGACGGAATTATGGTTGAAATCTGGAGCTACCATGCGTTGTTCTGGCTGGCGATAGGGATGTTGGGTATTGCGATGATTTGCCTGCTGTTTATTAAAGATATTTAACTCACCAGTCTGGGCTAAGAAAAATCg                                                                                                   >  1:787957/1‑139 (MQ=255)                             gTTGACGGAATTATGGTTGAAATCTGGAGCTACCATGCGTTGTTCTGGCTGGCGATAGGGATGTTGGGTATTGCGATGATTTGCCTGCTGTTTACTAAAGATATTTAACTCACCAGTCTGGGCTAAGAAAAATCGGCGt                                                                                               <  2:1012218/139‑1 (MQ=255)                                               gAAATCTGGAGCTACCATGCGTTGTTCTGGCTGGCGATAGGGATGTTGGGTATTGCGATGATTTGCCTGCTGTTTATTAAAGATATTTAACTCACCAGTCTGGGCTAAGAAAAATCGGc                                                                                                 <  1:385012/119‑1 (MQ=255)                                               gAAATCTGGAGCTACCATGCGTTGTTCTGGCTGGCGATAGGGATGTTGGGTATTGCGATGATTTGCCTGCTGTTTATTAAAGATATTTAACTCACCAGTCTGGGCTAAGAAAAATCGGc                                                                                                 >  2:385012/1‑119 (MQ=255)                                                 aaTCTGGAGCTACCATGCGTTGTTCTGGCTGGCGATAGGGATGTTGGGTATTGCGATGATTTGCCTGCTGTTTATTAAAGATATTTAACTCACCAGt                                                                                                                     <  1:767742/97‑1 (MQ=255)                                                 aaTCTGGAGCTACCATGCGTTGTTCTGGCTGGCGATAGGGATGTTGGGTATTGCGATGATTTGCCTGCTGTTTATTAAAGATATTTAACTCACCAGt                                                                                                                     >  2:767742/1‑97 (MQ=255)                                                            aCCATGCGTTGTTCTGGCTGGCGATAGGGATGTTGGGTATTGCGATGATTTGCCTGCTGTTTATTAAAGATATTTAACTCACCAGTCTGGGCTAAGAAAAATCGGCGTGGCGCAAACa                                                                                     <  1:1365377/118‑1 (MQ=255)                                                            aCCATGCGTTGTTCTGGCTGGCGATAGGGATGTTGGGTATTGCGATGATTTGCCTGCTGTTTATTAAAGATATTTAACTCACCAGTCTGGGCTAAGAAAAATCGGCGTGGCGCAAACa                                                                                     >  2:1365377/1‑118 (MQ=255)                                                            aCCATGCGTTGTTCTGGCTGGCGATAGGGATGTTGGGTATTGCGATGATTTGCCTGCTGTTTATTAAAGATATTTAACTCACCAGTCTGGGCTAAGAAAAATCGGCGTGGCGCAAACAGAGCCATCCTCGCCATTTTTg                                                                >  2:1419820/1‑139 (MQ=255)                                                              cATGCGTTGTTCTGGCTGGCGATAGGGATGTTGGGTATTGCGATGATTTGCCTGCTGTTTATTAAAGATATTTAACTCACCAGTCTGGGCTAAGAAAAATCGGCGTGGCGCAAACAGAGCCATCCTCGCCATTTTTGCt                                                              >  1:563590/1‑139 (MQ=255)                                                               aTGCGTTGTTCTGGCTGGCGATAGGGATGTTGGGTATTGCGATGATTTGCCTGCTGTTTATTAAAGATATTTAACTCACCAGTCTGGGCTAAGAAAAATCGGCGTGGCGCAAACAGAGCCATCCTCGCCATTTTTGCTa                                                             >  2:690086/1‑139 (MQ=255)                                                                    ttgttCTGGCTGGCGATAGGGATGTTGGGTATTGCGATGATTTGCCTGCTGTTTATTAAAGATATTTAACTCACCAGTCTGGGCTAAGAAAAATCGGCGTGGCGCAAACAGAGCCATCCTCGCCATTTTTGCTATCATg                                                        >  2:597148/1‑139 (MQ=255)                                                                       ttCTGGCTGGCGATAGGGATGTTGGGTATTGCGATGATTTGCCTGCTGTTTATTAAAGATATTTAACTCAc                                                                                                                         >  1:980725/1‑71 (MQ=255)                                                                       ttCTGGCTGGCGATAGGGATGTTGGGTATTGCGATGATTTGCCTGCTGTTTATTAAAGATATTTAACTCAc                                                                                                                         <  2:980725/71‑1 (MQ=255)                                                                       ttCTGGCTGGCGATAGGGATGTTGGGTATTGCGATGATTTGCCTGCTGTTTATTAAAGATATTTAACTCACCAGTCTGGGCTAAGAAAAATCGGCGTGGCGCAAACAGAGCCATCCTc                                                                          <  2:1305239/118‑1 (MQ=255)                                                                       ttCTGGCTGGCGATAGGGATGTTGGGTATTGCGATGATTTGCCTGCTGTTTATTAAAGATATTTAACTCACCAGTCTGGGCTAAGAAAAATCGGCGTGGCGCAAACAGAGCCATCCTc                                                                          >  1:1305239/1‑118 (MQ=255)                                                                       ttCTGGCTGGCGATAGGGATGTTGGGTATTGCGATGATTTGCCTGCTGTTTATTAAAGATATTTAACTCACCAGTCTGGGCTAAGAAAAATCGGCGTGGCGCAAACAGAGCCATCCTCGCCATTTTTGCTATCATGc                                                       <  2:149649/137‑1 (MQ=255)                                                                       ttCTGGCTGGCGATAGGGATGTTGGGTATTGCGATGATTTGCCTGCTGTTTATTAAAGATATTTAACTCACCAGTCTGGGCTAAGAAAAATCGGCGTGGCGCAAACAGAGCCATCCTCGCCATTTTTGCTATCATGc                                                       >  1:149649/1‑137 (MQ=255)                                                                       ttCTGGCTGGCGATAGGGATGTTGGGTATTGCGATGATTTGCCTGCTGTTTATTAAAGATATTTAACTCACCAGTCTGGGCTAAGAAAAATCGGCGTGGCGCAAACAGAGCCATCCTCGCCATTTTTGCTATCATGCCt                                                     <  1:690086/139‑1 (MQ=255)                                                                       ttCTGGCTGGCGATAGGGATGTTGGGTATTGCGATGATTTGCCTGCTGTTTATTAAAGATATTTAACTCACCAGTCTGGGCTAAGAAAAATCGGCGTGGCGCAAACAGAGCCATCCTCGCCATTTTTACTATCATTCCt                                                     >  1:926537/1‑139 (MQ=255)                                                                        tCTGGCTGGCGATAGGGATGTTGGGTATTGCGATGATTTGCCTGCTGTTTATTAAAGATATTTAACTCACCAGTCTGGGCTAAGAAAAATCGGCGTGGCGCAAACAGAGCTATCCTCGCCATTTTTGCTATCATGCCTg                                                    <  2:926537/139‑1 (MQ=255)                                                                        tCTGGCTGGCGATAGGGATGTTGGGTATTGCGATGATTTGCCTGCTGTTTATTAAAGATATTTAACTCACCAGTCTGGGCTAAGAAAAATCGGCGTGGCGCAAACAGAGCCATCCTCGCCATTTTTGCTATCATGCCTg                                                    <  1:1419820/139‑1 (MQ=255)                                                                            gctggcGATAGGGATGTTGGGTATTGCGATGATTTGCCTGCTGTTTATTAAAGATATTTAACTCAc                                                                                                                         <  1:1158089/66‑1 (MQ=255)                                                                            gctggcGATAGGGATGTTGGGTATTGCGATGATTTGCCTGCTGTTTATTAAAGATATTTAACTCAc                                                                                                                         >  2:1158089/1‑66 (MQ=255)                                                                                 cGATAGGGATGTTGGGTATTGCGATGATTTGCCTGCTGTTTATTAAAGATATTTAACTCACCAGTCTGGGCTAAGAAAAATCGGCGTGGCGCAAACAGAGCCATCCTCGCCATTTTTGCTAt                                                            >  2:1141413/1‑122 (MQ=255)                                                                                 cGATAGGGATGTTGGGTATTGCGATGATTTGCCTGCTGTTTATTAAAGATATTTAACTCACCAGTCTGGGCTAAGAAAAATCGGCGTGGCGCAAACAGAGCCATCCTCGCCATTTTTGCTAt                                                            <  1:1141413/122‑1 (MQ=255)                                                                                      gggATGTTGGGTATTGCGATGATTTGCCTGCTGTTTATTAAAGATATTTAACTCACCAGTCTGGGCTAAGAAAAATCGGCGTGGCGCAAACAGAGCCATCCTCGCCATTTTTg                                                                >  2:290837/1‑113 (MQ=255)                                                                                      gggATGTTGGGTATTGCGATGATTTGCCTGCTGTTTATTAAAGATATTTAACTCACCAGTCTGGGCTAAGAAAAATCGGCGTGGCGCAAACAGAGCCATCCTCGCCATTTTTg                                                                <  1:290837/113‑1 (MQ=255)                                                                                       ggATGTTGGGTATTGCGATGATTTGCCTGCTGTTTATTAAAGATATTTAACTCACCAGTCTGGGCTAAGAAAAATCGGCGTGGCGCAAACAGAGCCATCCTCGCCATTTTTGCTATCATGCCTGCATACATAAACGACa                                     <  1:703883/139‑1 (MQ=255)                                                                                          tGTTGGGTATTGCGATGATTTGCCTGCTGTTTATTAAAGATATTTAACTCACCAGTCTGGGCTAAGAAAAATCGGCGTGGCGCAAACAGAGCCATCCTCGCCATTTTTGCTATCATGCCTGCATACATAAACGACaaaa                                  <  2:563590/139‑1 (MQ=255)                                                                                                    tGCGATGATTTGCCTGCTGTTTATTAAAGATATTTAACTCACCAGTCTGGGCTAAGAAAAATCGGCGTGGCGCAAACAGAGCCATCCTCGCCATTTTTGCTATCATGCCTGCATACATAAACGACAAAACAGTATGCag                        >  2:1120408/1‑139 (MQ=255)                                                                                                    tGCGATGATTTGCCTGCTGTTTATTAAAGATATTTAACTCACCAGTCTGGGCTAAGAAAAATCGGCGTGGCGCAAACAGAGCCATCCTCGCCATTTTTGCTATCATGCCTGCATACATAAACGACAAAACAGTATGCag                        >  1:253803/1‑139 (MQ=255)                                                                                                         tgatTTGCCTGCTGTTTCTTAAAGATATTTAACTCACCAGTCTGGGCTAAGAAAAATCGGCGTGGCGCAAAGAGAGCCATCCTCGCCATTTTTGCTAt                                                            >  2:821212/1‑98 (MQ=255)                                                                                                         tgatTTGCCTGCTGTTTATTAAAGATATTTAACTCACCAGTCTGGGCTAAGAAAAATCGGCGTGGCGCAAACAGAGCCATCCTCGCCATTTTTGCTAt                                                            <  1:821212/98‑1 (MQ=255)                                                                                                           atTTGCCTGCTGTTTATTAAAGATATTTAACTCACCAGTCTGGGCTAAGAAAAATCGGCGTGGCGCAAACAGAGCCATCCTCGCCATTTTTGCTATCATGCCTGCATACATAAACGACAAAACAGTATGCAGAGGGaaa                 >  2:320373/1‑139 (MQ=255)                                                                                                                ccTGCTGTTTATTAAAGATATTTAACTCACCAGTCTGGGCTAAGAAAAATCGGCGTGGCGCAAACAGAGCCATCCTCGCCAttttt                                                                 >  1:974947/1‑86 (MQ=255)                                                                                                                ccTGCTGTTTATTAAAGATATTTAACTCACCAGTCTGGGCTAAGAAAAATCGGCGTGGCGCAAACAGAGCCATCCTCGCCAttttt                                                                 <  2:974947/86‑1 (MQ=255)                                                                                                                   gctgTTTATTAAAGATATTTAACTCACCAGTCTGGGCTAAGAAAAATCGGCGTGGCGCAAACAGAGCCATCCTCGCCATTTTTGCTATCATGCCTGCATACATAAACGACAAAACAGTATGCAGAGGGAAAAATGGGtt         <  2:253803/139‑1 (MQ=255)                                                                                                                    ctgTTTATTAAAGATATTTAACTCACCAGTCTGGGCTAAGAAAAATCGGCGTGGCGCAAACAGAGCCATCCTCGCCATTTTTGCTATCATg                                                        >  2:179766/1‑91 (MQ=255)                                                                                                                    ctgTTTATTAAAGATATTTAACTCACCAGTCTGGGCTAAGAAAAATCGGCGTGGCGCAAACAGAGCCATCCTCGCCATTTTTGCTATCATg                                                        <  1:179766/91‑1 (MQ=255)                                                                                                                          attaAAGATATTTAACTCACCAGTCTGGGCTAAGAAAAATCGGCGTGGCGCAAACAGAGCCATCCTCGCCATTTTTGCTATCATGCCTGCATACATAAACGACAAAACAGTATGCAGAGGGAAAAATGGGTTCCACCAg  >  1:290871/1‑139 (MQ=255)                                                                                                                                 atatTTAACTCACCAGTCTGGGCTAAGAAAAATCGGc                                                                                                 <  1:286526/37‑1 (MQ=255)                                                                                                                                 atatTTAACTCACCAGTCTGGGCTAAGAAAAATCGGc                                                                                                 >  2:286526/1‑37 (MQ=255)                                                                                                                                 atatTTAACTCACCAGTCTGGGCTAAGAAAAATCGGCGTGGCGCAAACAGAGCCATCCTCGCCATTTTTGCTATCATg                                                        >  1:1423350/1‑78 (MQ=255)                                                                                                                                 atatTTAACTCACCAGTCTGGGCTAAGAAAAATCGGCGTGGCGCAAACAGAGCCATCCTCGCCATTTTTGCTATCATg                                                        <  2:1423350/78‑1 (MQ=255)                                                                                                                                 atatTTAACTCACCAGTCTGGGCTAAGAAAAATCGGCGTGGCGCAAACAGAGCCATCCTCGCCATTTTTGCTATCATGCCTGCATACATAAACGACAAAACAGTATGCAGAGGGAAAAATGGGTTc        <  2:333194/126‑1 (MQ=255)                                                                                                                                 atatTTAACTCACCAGTCTGGGCTAAGAAAAATCGGCGTGGCGCAAACAGAGCCATCCTCGCCATTTTTGCTATCATGCCTGCATACATAAACGACAAAACAGTATGCAGAGGGAAAAATGGGTTc        >  1:333194/1‑126 (MQ=255)                                                                                                                                   atTTAACTCACCAGTCTGGGCTAAGAAAAATCGGCGTGGCGCAAACAGAGCCATCCTCGCCATTTTTGc                                                               <  2:1315464/69‑1 (MQ=255)                                                                                                                                   atTTAACTCACCAGTCTGGGCTAAGAAAAATCGGCGTGGCGCAAACAGAGCCATCCTCGCCATTTTTGc                                                               >  1:1315464/1‑69 (MQ=255)                                                                                                                                        |                                                                                                                             ACGCGTCGGCTGGATTATCGCCGGCTCTGTTGACGGAATTATGGTTGAAATCTGGAGCTACCATGCGTTGTTCTGGCTGGCGATAGGGATGTTGGGTATTGCGATGATTTGCCTGCTGTTTATTAAAGATATTTAGCTCACCAGTCTGGGCTAAGAAAAATCGGCGTGGCGCAAACAGAGCCATCCTCGCCATTTTTGCTATCATGCCTGCATACATAAACGACAAAACAGTATGCAGAGGGAAAAATGGGTTCCACCAG  >  NC_000913/3838002‑3838261

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