Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,175,651	G→A	S119L (TCA→TTA)	yqiA ←	acyl CoA esterase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,175,651	0	G	A	100.0%	181.8 / NA	45	S119L (TCA→TTA)	yqiA	acyl CoA esterase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (26/19);  total (26/19)

TTCTATGACAGTCTGGCGGCAGGAAGCGTAGTACGCCACCGCCTGGCGGTAATCCAGCACTTCATCTCCCGTCTGTTGCAGCAGCCAGATCAAATCCGGCGCTTCCAGCGGGTCAATCTGCATGACTTTAAGATCGTAAATATGGCGTGACTCTAGCACATATTGCTGCCCGGTGTAGGGGTTCTCGTTCTGACCGAGATAGTCCGTCAGCAGTTCAAACGGGCGCACCGCCGGGTTTACCACCACTGCGGGCAGCATAAAACATTGTGACAACCAGGTGGCGTAATAT  >  NC_000913/3175503‑3175791                                                                                                                                                     |                                                                                                                                             ttCTATGACAGTCTGGCGGCAGGAAGCGTAGTACGCCACCGCCTGGCGGTAATCCAGCACTTCATCTCCCGTCTGTTGCAGCAGCCAGATCAAATCCGGCGCTTCCAGCGGGTCAATCTGCATGACTTTAAGATCGTAAATATGGCGTAAc                                                                                                                                            >  2:423507/1‑151 (MQ=255)        aCAGTCTGGCGGCAGGAAGCGTAGTACGCCACCGCCTGGCGGTAATCCAGCACTTCATCTCCCGTCTGTTGCAGCAGCCAGATCAAATCCGGCGCTTCCAGCGGGTCAATCTGCATGACTTTAAGATCGTAAATATGGCGTAACTCTAGca                                                                                                                                     <  2:554165/151‑1 (MQ=255)            tCTGGCGGCAGGAAGCGTAGTACGCCACCGCCTGGCGGTAATCCAGCACTTCATCTCCCGTCTGTTGCAGCAGCCAGATCAAATCCGGCGCTTCCAGCGGGTCAATCTGCATGACTTTAAGATCGTAAATATGGCGTAACTCTAGCACat                                                                                                                                  >  1:449308/1‑150 (MQ=255)              tGGCGGCAGGAAGCGTAGTACGCCACCGCCTGGCGGTAATCCAGCACTTCATCTCCCGTCTGTTGCAGCAGCCAGATCAAATCCGGCGCTTCCAGCGGGTCAATCTGCATGACTTTAAGATCGTAAATATGGCGTAACTCTAGCACATAtt                                                                                                                               >  1:109049/1‑151 (MQ=255)               ggcggcAGGAAGCGTAGTACGCCACCGCCTGGCGGTAATCCAGCACTTCATCTCCCGTCTGTTGCAGCAGCCAGATCAAATCCGGCGCTTGCAGCGGGTCAATCTGCATGACGTTAAGATCGTAAATATGGCGTAACTCTAGCACatat                                                                                                                                <  1:266575/149‑1 (MQ=255)               ggcggcAGGAAGCGTAGTACGCCACCGCCTGGCGGTAATCCAGCACTTCATCTCCCGTCTGTTGCAGCAGCCAGATCAAATCCGGCGCTTCCAGCGGGTCAATCTGCATGACTTTAAGATCGTAAATATGGCGTAACTCTAGCACatat                                                                                                                                <  1:42190/149‑1 (MQ=255)                            gtagtaCGCCACCGCCTGGCGGTAATCCAGCACTTCATCTCCCGTCTGTTGCAGCAGCCAGATCAAATCCGGCGCTTCCAGCGGGTCAATCTGCATGACTTTAAGATCGTAAATATGGCGTAACTCTAGCACATATTGCTGCCCGGTGTAg                                                                                                                 <  2:333068/151‑1 (MQ=255)                                        cGCCTGGCGGTAATCCAGCACTTCATCTCCCGTCTGTTGCAGCAGCCAGATCAAATCCGGCGCTTCCAGCGGGTCAATCTGCATGACTTTAAGATCGTAAATATGGCGTAACTCTAGCACATATTGCTGCCCGGTGTAGGGGTTCTCGTTc                                                                                                     <  2:482136/151‑1 (MQ=255)                                                 gTAATCCAGCACTTCATCTCCCGTCTGTTGCAGCAGCCAGATCAAATCCGGCGCTTCCAGCGGGTCAATCTGCATGACTTTAAGATCGTAAATATGGCGTAACTCTAGCACATATTGCTGCCCGGTGTAGGGGTTCTCGTTCTGACCg                                                                                               >  1:641362/1‑148 (MQ=255)                                                       cAGCACTTCATCTCCCGTCTGTTGCAGCAGCCAGATCAAATCCGGCGCTTCCAGCGGGTCAATCTGCATGACTTTAAGATCGTAAATATGGCGTAACTCTAGCACATATTGCTGCCCGGAGTAGGGGTTCTCGTTCTGACCGAGATAGt                                                                                        >  1:665749/1‑149 (MQ=255)                                                                aTCTCCCGTCTGTTGCAGCAGCCAGATCAAATCCGGCGCTTCCAGCGGGTCAATCTGCATGACTTTAAGATCGTAAATATGGCGTAACTCTAGCACATATTGCTGCCCGGTGTAGGGGTTCTCGTTCTGACCGAGATAGTCCGTCAGCAGt                                                                             >  1:252690/1‑151 (MQ=255)                                                                    cccGTCTGTTGCAGCAGCCAGATCAAATCCGGCGCTTCCAGCGGGTCAATCTGCATGACTTTAAGATCGTAAATATGGCGTAACTCTAGCACATATTGCTGCCCGGTGTAGGGGTTCTCGTTCTGACCGAGATAGTCCGTCAGCAGTTCa                                                                          <  2:250781/150‑1 (MQ=255)                                                                    cccGTCTGTTGCAGCAGCCAGATCAAATCCGGCGCTTCCAGCGGGTCAATCTGCATGACTTTAAGATCGTAAATATGGCGTAACTCTAGCACATATTGCTGCCCGGTGTAGGGGTTCTCGTTCTGACCGAGATAGTCCGTCAGCAGTTCa                                                                          >  1:650461/1‑150 (MQ=255)                                                                        tCTGTTGCAGCAGCCAGATCAAATCCGGCGCTTCCAGCGGGTCAATCTGCATGACTTTAAGATCGTAAATATGGCGTAACTCTAGCACATATTGCTGCCCGGTGTAGGGGTTCTCGTTCTGACCGAGATAGTCCGTCAGCAGTTCAAACgg                                                                     >  2:71912/1‑151 (MQ=255)                                                                           gTTGCAGCAGCCAGATCAAATCCGGCGCTTCCAGCGGGTCAATCTGCATGACTTTAAGATCGTAAATATGGCGTAACTCTAGCACATATTGCTGCCCGGTGTAGGGGTTCTCGTTCTGACCGAGATAGTCCGTCAGCAGTTCAAACGGgcg                                                                  <  2:107514/151‑1 (MQ=255)                                                                              gcagcaGCCAGATCAAATCCGGCGCTTCCAGCGGGTCAATCTGCATGACTTTAAGATCGTAAATATGGCGTAACTCTAGCACATATTGCTGCCCGGTGTAGGGGTTCTCGTTCTGACCGAGATAGTCCGTCAGCAGTTCAAACGGGCGCa                                                                <  1:153576/150‑1 (MQ=255)                                                                                agcagcCAGATCAAATCCGGCGCTTCCAGCGGGTCAATCTGCATGACTTTAAGATCGTAAATATGGCGTAACTCTAGCACATATTGCTGCCCGGTGTAGGGGTTCTCGTTCTGACCGAGATAGTCCGTCAGCAGTTCAAACGGGCGCAccg                                                             <  1:310533/151‑1 (MQ=255)                                                                                agcagcCAGATCAAATCCGGCGCTTCCAGCGGGTCAATCTGCATGACTTTAAGATCGTAAATATGGCGTAACTCTAGCACATATTGCTGCCCGGTGTAGGGGTTCTCGTTCTGACCGAGATAGTCCGTCAGCAGTTCAAACGGGCGCAccg                                                             <  1:261520/151‑1 (MQ=255)                                                                                    gcCAGATCAAATCCGGCGCTTCCAGCGGGTCAATCTGCATGACTTTAAGATCGTAAATATGGCGTAACTCTAGCACATATTGCTGCCCGGTGTAGGGGTTCTCGTTCTGACCGAGATAGTCCGTCAGCAGTTCAAACGGGCGCACCGCCgg                                                         <  2:637703/151‑1 (MQ=255)                                                                                         aTCAAATCCGGCGCTTCCAGCGGGTCAATCTGCATGACTTTAAGATCGTAAATATGGCGTAACTCTAGCACATATTGCTGCCCGGTGTAGGGGTTCTCGTTCTGACCGAGATAGTCCGTCAGCAGTTCAAACGGGCGCACCGCCGGGTTTa                                                    <  2:242295/151‑1 (MQ=255)                                                                                            aaaTCCGGCGCTTCCAGCGGGTCAATCTGCATGACTTTAAGATCGTAAATATGGCGTAACTCTAGCACATATTGCTGCCCGGTGTAGGGGTTCTCGTTCTGACCGAGATAGTCCGTCAGCAGTTCAAACGGGCGCACCGCCGGGTTTacc                                                  >  2:294738/1‑150 (MQ=255)                                                                                                 cGGCGCTTCCAGCGGGTCAATCTGCATGACTTTAAGATCGTAAATATGGCGTAACTCTAGCACATATTGCTGCCCGGTGTAGGGGTTCTCGTTCTGACCGAGATAGTCCGTCAGCAGTTCAAACGGGCGCACCGCCGGGTTTACCACCACt                                            >  2:371871/1‑151 (MQ=255)                                                                                                 cGGCGCTTCCAGCGGGTCAATCTGCATGACTTTAAGATCGTAAATATGGCGTAACTCTAGCACATATTGCTGCCCGGTGTAGGGGTTCTCGTTCTGACCGAGATAGTCCGTCAGCAGTTCAAACGGGCGCACCGCCGGGTTTACCACCACt                                            >  1:138606/1‑151 (MQ=255)                                                                                                 cGGCGCTTCCAGCGGGTCAATCTGCATGACTTTAAGATCGTAAATATGGCGTAACTCTAGCACATATTGCTGCCCGGTGTAGGGGTTCTCGTTCTGACCGAGATAGTCCGTCAGCAGTTCAAACGGGCGCACCGCCGGGTTTACCACCACt                                            >  1:139555/1‑151 (MQ=255)                                                                                                    cgcTTCCAGCGGGTCAATCTGCATGACTTTAAGATCGTAAATATGGCGTAACTCTAGCACATATTGCTGCCCGGTGTAGGGGTTCTCGTTCTGACCGAGATAGTCCGTCAGCAGTTCAAACGGGCGCACCGCCGGGTTTACCACCACTGCg                                         <  1:24117/151‑1 (MQ=255)                                                                                                        tCCAGCGGGTCAATCTGCATGACTTTAAGATCGTAAATATGGCGTAACTCTAGCACATATTGCTGCCCGGTGTAGGGGTTCTCGTTCTGACCGAGATAGTCCGTCAGCAGTTCAAACGGGCGCACCGCCGGGTTTACCACCACTGCggg                                       >  1:76508/1‑149 (MQ=255)                                                                                                        tCCAGCGGGTCAATCTGCATGACTTTAAGATCGTAAATATGGCGTAACTCTAGCACATATTGCTGCCCGGTGTAGGGGTTCTCGTTCTGACCGAGATAGTCCGTCAGCAGTTCAAACGGGCGCACCGCCGGGTTTACCACCACTGCGGgc                                      >  1:432304/1‑150 (MQ=255)                                                                                                        tCCAGCGGGTCAATCTGCATGACTTTAAGATCGTAAATATGGCGTAACTCTAGCACATATTGCTGCCCGGTGTAGGGGTTCTCGTTCTGACCGAGATAGTCCGTCAGCAGTTCAAACGGGCGCACCGCCGGGTTTACCACCACTGCGGgc                                      <  1:404822/150‑1 (MQ=255)                                                                                                           aGCGGGTCAATCTGCATGACTTTAAGATCGTAAATATGGCGTAACTCTAGCACATATTGCTGCCCGGTGTAGGGGTTCTCGTTCTGACCGAGATAGTCCGTCAGCAGTTCAAACGGGCGCACCGCCGGGTTTACCACCACTGCGGgcagc                                   >  2:415325/1‑150 (MQ=255)                                                                                                            gCGGGTCAATCTGCATGACTTTAAGATCGTAAATATGGCGTAACTCTAGCACATATTGCTGCCCGGTGTAGGGGTTCTCGTTCTGACCGAGATAGTCCGTCAGCAGTTCAAACGGGCGCACCGCCGGGTTTACCACCACTGCGGGCAGCAt                                 <  1:54899/151‑1 (MQ=255)                                                                                                             cGGGTCAATCTGCATGACTTTAAGATCGTAAATATGGCGTAACTCTAGCACATATTGCTGCCCGGTGTAGGGGTTCTCGTTCTGACCGAGATAGTCCGTCAGCAGTTCAAACGGGCGCACCGCCGGGTTTACCACCACTGCGGGCAGCAt                                 >  1:547506/1‑150 (MQ=255)                                                                                                             cGGGTCAATCTGCATGACTTTAAGATCGTAAATATGGCGTAACTCTAGCACATATTGCTGCCCGGTGTAGGGGTTCTCGTTCTGACCGAGATAGTCCGTCAGCAGTTCAAACGGGCGCACCGCCGGGTTTACCACCACTGCGGGCAGCATa                                <  2:384723/151‑1 (MQ=255)                                                                                                               ggTCAATCTGCATGACTTTAAGATCGTAAATATGGCGTAACTCTAGCACATATTGCTGCCCGGTGTAGGGGTTCTCGTTCTGACCGAGATAGTCCGTCAGCAGTTCAAACGGGCGCACCGCCGGGTTTACCACCACTGCGGGCAGCATaa                               >  1:497515/1‑150 (MQ=255)                                                                                                                  cAATCTGCATGACTTTAAGATCGTAAATATGGCGTAACTCTAGCACATATTGCTGCCCGGTGTAGGGGTTCTCGTTCTGACCGAGATAGTCCGTCAGCAGTTCAAACGGGCGCACCGCCGGGTTTACCACCACTGCGGGCAGCATAAAACa                           <  1:433765/151‑1 (MQ=255)                                                                                                                        gCATGACTTTAAGATCGTAAATATGGCGTAACTCTAGCACATATTGCTGCCCGGTGTAGGGGTTCTCGTTCTGACCGAGATAGTCCGTCAGCAGTTCAAACGGGCGCACCGCCGGGTTTACCACCACTGCGGGCAGCATAAAACATtgtg                      >  1:392315/1‑150 (MQ=255)                                                                                                                         cATGACTTTAAGATCGTAAATATGGCGTAACTCTAGCACATATTGCTGCCCGGTGTAGGGGTTCTCGTTCTGACCGAGATAGTCCGTCAGCAGTTCAAACGGGCGCACCGCCGGGTTTACCACCACTGCGGGCAGCATAAAACATTGTGAc                    >  1:621032/1‑151 (MQ=255)                                                                                                                          aTGACTTTAAGATCGTAAATATGGCGTAACTCTAGCACATATTGCTGCCCGGTGTAGGGGTTCTCGTTCTGACCGAGATAGTCCGTCAGCAGTTCAAACGGGCGCACCGCCGGGTTTACCACCACTGCGGGCAGCATAAAACATTGTGAc                    >  2:386845/1‑150 (MQ=255)                                                                                                                          aTGACTTTAAGATCGTAAATATGGCGTAACTCTAGCACATATTGCTGCCCGGTGTAGGGGTTCTCGTTCTGACCGAGATAGTCCGTCAGCAGTTCAAACGGGCGCACCGCCGGGTTTACCACCACTGCGGGCAGCATAAAACATTGTGACa                   <  2:1953/151‑1 (MQ=255)                                                                                                                               tttAAGATCGTAAATATGGCGTAACTCTAGCACATATTGCTGCCCGGTGTAGGGGTTCTCGTTCTGACCGAGATAGTCCGTCAGCAGTTCAAACGGGCGCACCGCCGGGTTTACCACCACTGCGGGCAGCATAAAACATTGTGACAAcc                >  2:453991/1‑149 (MQ=255)                                                                                                                                  aaGATCGTAAATATGGCGTAACTCTAGCACATATTGCTGCCCGGTGTAGGGGTTCTCGTTCTGACCGAGATAGTCCGTCAGCAGTTCAAACGGGCGCACCGCCGGGTTTACCACCACTGCGGGCAGCATa                                <  1:261791/130‑1 (MQ=255)                                                                                                                                  aaGATCGTAAATATGGCGTAACTCTAGCACATATTGCTGCCCGGTGTAGGGGTTCTCGTTCTGACCGAGATAGTCCGTCAGCAGTTCAAACGGGCGCACCGCCGGGTTTACCACCACTGCGGGCAGCATa                                >  2:261786/1‑130 (MQ=255)                                                                                                                                   aGATCGTAAATATGGCGTAACTCTAGCACATATTGCTGCCCGGTGTAGGGGTTCTCGTTCTGACCGAGATAGTCCGTCAGCAGTTCAAACGGGCGCACCGCCGGGTTTACCACCACTGCGGGCAGCATAAAACATTGTGACAACCAGGTgg          >  1:668299/1‑151 (MQ=255)                                                                                                                                       cGTAAATATGGCGTAACTCTAGCACATATTGCTGCCCGGTGTAGGGGTTCTCGTTCTGACCGAGATAGTCCGTCAGCAGTTCAAACGGGCGCACCGCCGGGTTTACCACCACTGCGGGCAGCATAAAACATTGTGACAACCAGGTGGCGt       >  2:235398/1‑150 (MQ=255)                                                                                                                                       cGTAAATATGGCGTAACTCTAGCACATATTGCTGCCCGGTGTAGGGGTTCTCGTTCTGACCGAGATAGTCCGTCAGCAGTTCAAACGGGCGCACCGCCGGGTTTACCACCACTGCGGGCAGCATAAAACATTGTGACAACCAGGTGGCGt       >  2:87136/1‑150 (MQ=255)                                                                                                                                           aaTATGGCGTAACTCTAGCACATATTGCTGCCCGGTGTAGGGGTTCTCGTTCTGACCGAGATAGTCCGTCAGCAGTTCAAACGGGCGCACCGCCGGGTTTACCACCACTGCGGGCAGCATAAAACATTGTGACAACCAGGTGGCGTAatat  >  2:200995/1‑151 (MQ=255)                                                                                                                                                     |                                                                                                                                             TTCTATGACAGTCTGGCGGCAGGAAGCGTAGTACGCCACCGCCTGGCGGTAATCCAGCACTTCATCTCCCGTCTGTTGCAGCAGCCAGATCAAATCCGGCGCTTCCAGCGGGTCAATCTGCATGACTTTAAGATCGTAAATATGGCGTGACTCTAGCACATATTGCTGCCCGGTGTAGGGGTTCTCGTTCTGACCGAGATAGTCCGTCAGCAGTTCAAACGGGCGCACCGCCGGGTTTACCACCACTGCGGGCAGCATAAAACATTGTGACAACCAGGTGGCGTAATAT  >  NC_000913/3175503‑3175791

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