Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,175,815	A→G	G64G (GGT→GGC)	yqiA ←	esterase YqiA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,175,815	0	A	G	100.0%	135.8 / NA	42	G64G (GGT→GGC)	yqiA	esterase YqiA
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (23/19);  total (23/19)

CTCGTTCTGACCGAGATAGTCCGTCAGCAGTTCAAACGGGCGCACCGCCGGGTTTACCACCACTGCGGGCAGCATAAAACATTGTGACAACCAGGTGGCGTAATATCCCCCCAGTGACGAACCGACAATACCCAGCGAATCACCGCCATGTTCCAGGACAATGGATTCCAGCAGCTCTGCCGCGTCGGAAGGATACGGCGGCAACTGCGGAATGATCATCTCAACGTCAGGGTGATGTTCCGCCAGCCAGTTTTTTAACAAGCT  >  NC_000913/3175686‑3175949                                                                                                                                  |                                                                                                                                       ctcGTTCTGACCGAGATAGTCCGTCAGCAGTTCAAACGGGCGCACCGCCGGGTTTACCACCACTGCGGGCAGCATAAAACATTGTGACAACCAGGTGGCGTAATATCCCCCCAGTGACGAACCGACAATGCCCAGCGAAt                                                                                                                              >  2:635397/1‑140 (MQ=255) ctcGTTCTGACCGAGATAGTCCGTCAGCAGTTCAAACGGGCGCACCGCCGGGTTTACCACCACTGCGGGCAGCATAAAACATTGTGACAACCAGGTGGCGTAATATCCCCCCAGTGACGAACCGACAATGCCCAGCGAAt                                                                                                                              >  2:613354/1‑140 (MQ=255) ctcGTTCTGACCGAGATAGTCCGTCAGCAGTTCAAACGGGCGCACCGCCGGGTTTACCACCACTGCGGGCAGCATAAAACATTGTGACAACCAGGTGGCGTAATATCCCCCCAGTGACGAACCGACAATGCCCAGCGAAt                                                                                                                              >  1:509861/1‑140 (MQ=255)   cGTTCTGACCGAGATAGTCCGTCAGCAGTTCAAACGGGCGCACCGCCGGGTTTACCACCACTGCGGGCAGCATAAAACATTGTGACAACCAGGTGGCGTAATATCCCCCCAGTGACGAACCGACAATGCCCAGCGAATCa                                                                                                                            >  2:650076/1‑140 (MQ=255)         gACCGAGATAGTCCGTCAGCAGTTCAAACGGGCGCACCGCCGGGTTTACCACCACTGCGGGCAGCATAAAACATTGTGACAACCAGGTGGCGTAATATCCCCCCAGTGACGAACCGACAATGCCCAGCGAATCACCGCCa                                                                                                                      <  1:119637/140‑1 (MQ=255)         gACCGAGATAGTCCGTCAGCAGTTCAAACGGGCGCACCGCCGGGTTTACCACCACTGCGGGCAGCATAAAACATTGTGACAACCAGGTGGCGTAATATCCCCCCAGTGACGAACCGACAATGCCCAGCGAATCACCGCCa                                                                                                                      <  2:509861/140‑1 (MQ=255)         gACCGAGATAGTCCGTCAGCAGTTCAAACGGGCGCACCGCCGGGTTTACCACCACTGCGGGCAGCATAAAACATTGTGACAACCAGGTGGCGTAATATCCCCCCAGTGACGAACCGACAATGCCCAGCGAATCACCGCCa                                                                                                                      <  1:773057/140‑1 (MQ=255)         gACCGAGATAGTCCGTCAGCAGTTCAAACGGGCGCACCGCCGGGTTTACCACCACTGCGGGCAGCATAAAACATTGTGACAACCAGGTGGCGTAATATCCCCCCAGTGACGAACCGACAATGCCCAGCGAATCACCGCCa                                                                                                                      >  1:697204/1‑140 (MQ=255)         gACCGAGATAGTCCGTCAGCAGTTCAAACGGGCGCACCGCCGGGTTTACCACCACTGCGGGCAGCATAAAACATTGTGACAACCAGGTGGCGTAATATCCCCCCAGTGACGAACCGACAATGCCCAGCGAATCACCGCCa                                                                                                                      <  1:650076/140‑1 (MQ=255)          aCCGAGATAGTCCGTCAGCAGTTCAAACGGGCGCACCGCCGGGTTTACCACCACTGCGGGCAGCATAAAACATTGTGACAACCAGGTGGCGTAATATCCCCCCAGTGACGAACCGACAATGCCCAGCGaa                                                                                                                               <  2:231845/130‑1 (MQ=255)          aCCGAGATAGTCCGTCAGCAGTTCAAACGGGCGCACCGCCGGGTTTACCACCACTGCGGGCAGCATAAAACATTGTGACAACCAGGTGGCGTAATATCCCCCCAGTGACGAACCGACAATGCCCAGCGaa                                                                                                                               >  1:231845/1‑130 (MQ=255)                    tCCGTCAGCAGTTCAAACGGGCGCACCGCCGGGTTTACCACCACTGCGGGCAGCATAAAACATTGTGACAACCAGGTGGCGTAATATCCCCCCAGTGACGAACCGACAATGCCCAGCGAATCACCGCCATGTTCCAGGAc                                                                                                           >  2:170905/1‑140 (MQ=255)                        tCAGCAGTTCAAACGGGCGCACCGCCGGGTTTACCACCACTGCGGGCAGCATAAAACATTGTGACAACCAGGTGGCGTAATATCCCCCCAGTGACGAACCGACAATGCCCAGCGAATCACCGCCATGTTCCAGGACAATg                                                                                                       >  1:540333/1‑140 (MQ=255)                                      gggCGCACCGCCGGGTTTACCACCACTGCGGGCAGCATAAAACATTGTGACAACCAGGTGGCGTAATATCCCCCCAGTGACGAACCGACAATGCCCAGCGAATCACCGCCATGTTCCAGGACAATGGATTCCAGCAGctc                                                                                         >  1:657236/1‑140 (MQ=255)                                         cgcACCGCCGGGTTTACCACCACTGCGGGCAGCATAAAACATTGTGACAACCAGGTGGCGTAATATCCCCCCAGTGACGAACCGACAATGCCCAGCGAATCACCGCCATGTTCCAGGACAATGGATTCCAGCAGCTCTGc                                                                                      >  2:118560/1‑140 (MQ=255)                                                 ggggTTTACCACCACTGCGGGCAGCAAAAAACATTGTGAAAACCAGGTGGCGTAATACCCCCCCAGTGACGAACCGACAATGCCCGGCGAATCACCGCCATGTTCCAGGACAATGGATTCCAGCAGCTCTGCCGCGTCgg                                                                              <  2:628299/139‑1 (MQ=255)                                                   ggTTTACCACCACTGCGGGCAGCATAAAACATTGTGACAACCAGGTGGCGTAATATCCCCCCAGTGACGAACCGACAATGCCCAGCGAATCACCGCCATGTTCCAGGACAATGGATTCCAGCAGCTCTGCCGCGTCGGaa                                                                            >  1:615028/1‑140 (MQ=255)                                                      ttACCACCACTGCGGGCAGCATAAAACATTGTGACAACCAGGTGGCGTAATATCCCCCCAGTGACGAACCGACAATGCCCAGCGAATCACCGCCATGTTCCAGGACAATGGATTCCAGCAGCTCTGCCGCGTCGGAAGGa                                                                         <  1:118560/140‑1 (MQ=255)                                                       tACCACCACTGCGGGCAGCATAAAACATTGTGACAACCAGGTGGCGTAATATCCCCCCAGTGACGAACCGACAATGCCCAGCGAATCACCGCCATGTTCCAGGACAATGGATTCCAGCAGCTCTGCCGCGTCGGAAGGAt                                                                        <  2:697204/140‑1 (MQ=255)                                                       tACCACCACTGCGGGCAGCATAAAACATTGTGACAACCAGGTGGCGTAATATCCCCCCAGTGACGAACCGACAATGCCCAGCGAATCACCGCCATGTTCCAGGACAATGGATTCCAGCAGCTCTGCCGCGTCGGAAGGAt                                                                        <  2:480440/140‑1 (MQ=255)                                                         ccaccaCTGCGGGCAGCATAAAACATTGTGACAACCAGGTGGCGTAATATCCCCCCAGTGACGAACCGACAATGCCCAGCGAATCACCGCCATGTTCCAGGACAATGGATTCCAGCAGCTCTGCCGCGTCGGAAGGATAc                                                                      <  1:170905/140‑1 (MQ=255)                                                               cTGCGGGCAGCATAAAACATTGTGACAACCAGGTGGCGTAATATCCCCCCAGTGACGAACCGACAATGCCCAGCGAATCACCGCCATGTTCCAGGACAATGGATTCCAGCAGCTCTGCCGCGTCGGAAGGATACggcggc                                                                >  1:122912/1‑140 (MQ=255)                                                               cTGCGGGCAGCATAAAACATTGTGACAACCAGGTGGCGTAATATCCCCCCAGTGACGAACCGACAATGCCCAGCGAATCACCGCCATGTTCCAGGACAATGGATTCCAGCAGCTCTGCCGCGTCGGAAGGATACggcggc                                                                <  2:615028/140‑1 (MQ=255)                                                                  cGGGCAGCATAAAACATTGTGACAACCAGGTGGCGTAATATCCCCCCAGTGACGAACCGACAATGCCCAGCGAATCACCGCCATGTTCCAGGACAATGGATTCCAGCAGCTCTGCCGCGTCGGAAGGATACGGCGGCAAc                                                             >  1:469353/1‑140 (MQ=255)                                                                      cagcaTAAAACATTGTGACAACCAGGTGGCGTAATATCCCCCCAGTGACGAACCGACAATGCCCAGCGAATCACCGCCATGTTCCAGGACAATGGATTCCAGCAGCTCTGCCGCGTCGGAAGGATACGGCGGCAACTGCg                                                         <  1:191461/140‑1 (MQ=255)                                                                            aaaaCATTGTGACAACCAGGTGGCGTAATACCCCCCCAGTGACGAACCGACAATGCCCAGCGAATCACCGCCATGTTCCAGGACAATGGATTCCAGCAGCTCTGCCGCGTCGGAAGGATACGGCGGCAACTGCGGAATGa                                                   <  1:606852/140‑1 (MQ=255)                                                                              aaCATTGTGACAACCAGGTGGCGTAATATCCCCCCAGTGACGAACCGACAATGCCCAGCGAATCACCGCCATGTTCCAGGACAATGGATTCCAGCAGCTCTGCCGCGTCGGAAGGATACGGCGGCAACTGCGGAATGatc                                                 <  2:122912/140‑1 (MQ=255)                                                                                  ttGTGACAACCAGGTGGCGTAATATCCCCCCAGTGACGAACCGACAATGCCCAGCGAATCACCGCCATGTTCCAGGACAATGGATTCCAGCAGCTCTGCCGCGTCGGAAGGATACGGCGGCAACTGCGGAATGATCAtct                                             <  2:596627/140‑1 (MQ=255)                                                                                   tgtgACAACCAGGTGGCGTAATATCCCCCCAGTGACGAACCGACAATGCCCAGCGAATCACCGCCATGTTCCAGGACAATGGATTCCAGCAGCTCTGCCGCGTCGGAAGGATACGGCGGCAACTGCg                                                         >  2:375581/1‑127 (MQ=255)                                                                                   tgtgACAACCAGGTGGCGTAATATCCCCCCAGTGACGAACCGACAATGCCCAGCGAATCACCGCCATGTTCCAGGACAATGGATTCCAGCAGCTCTGCCGCGTCGGAAGGATACGGCGGCAACTGCg                                                         <  1:375581/127‑1 (MQ=255)                                                                                   tgtgACAACCAGGTGGCGTAATATCCCCCCAGTGACGAACCGACAATGCCCAGCGAATCACCGCCATGTTCCAGGACAATGGATTCCAGCAGCTCTGCCGCGTCGGAAGGATACGGCGGCAACTGCGGAATGATCAtctc                                            <  1:613354/140‑1 (MQ=255)                                                                                   tgtgACAACCAGGTGGCGTAATATCCCCCCAGTGACGAACCGACAATGCCCAGCGAATCACCGCCATGTTCCAGGACAATGGATTCCAGCAGCTCTGCCGCGTCGGAAGGATACGGCGGCAACTGCGGAATGATCAtctc                                            >  1:440921/1‑140 (MQ=255)                                                                                        cAACCAGGTGGCGTAATACCCCCCCAGTGACGAACCGACAATGCCCAGCGAATCACCGCCATGTTCCAGGACAATGGATTCCAGCAGCTCTGCCGCGTCGGAAGGATACGGCGGCAACTGCGGAATGATCATCTCAACGt                                       <  2:469353/140‑1 (MQ=255)                                                                                          aCCAGGTGGCGTAATATCCCCCCAGTGACGAACCGACAATGCCCAGCGAATCACCGCCATGTTCCAGGACAATGGATTCCAGCAGCTCTGCCGCGTCGGAAGGATACGGCGGCAACTGCGGAATGATCATCTCAACGTCa                                     >  1:676336/1‑140 (MQ=255)                                                                                               gTGGCGTAATATCCCCCCAGTGACGAACCGACAATGCCCAGCGAATCACCGCCATGTTCCAGGACAATGGATTCCAGCAGCTCTGCCGCGTCGGAAGGATACGGCGGCAACTGCGGAATGATCATCTCAACGTCAggggg                                >  1:258302/1‑138 (MQ=255)                                                                                                          tCCCCCCAGTGACGAACCGACAATGCCCAGCGAATCACCGCCATGTTCCAGGACAATGGATTCCAGCAGCTCTGCCGCGTCGGAAGGATACGGCGGCAACTGCGGAATGATCATCTCAACGTCAGGGTGATGTTCCgcca                     >  1:769399/1‑140 (MQ=255)                                                                                                             ccccAGTGACGAACCGACAATGCCCAGCGAATCACCGCCATGTTCCAGGACAATGGATTCCAGCAGCTCTGCCGCGTCGGAAGGATACGGCGGCAACTGCGGAATGATCATCTCAACGTCAGGGTGATGTTCCgccagcc                  >  2:280120/1‑140 (MQ=255)                                                                                                               ccAGTGACGAACCGACAATGCCCAGCGAATCACCGCCATGTTCCAGGACAATGGATTCCAGCAGCTCTGCCGCGTCGGAAGGATACGGCGGCAACTGCGGAATGATCATCTCAACGTCAGGGTGATGTTCCGccagccag                >  2:381911/1‑140 (MQ=255)                                                                                                                   tGACGAACCGACAATGCCCAGCGAATCACCGCCATGTTCCAGGACAATGGATTCCAGCAGCTCTGCCGCGTCGGAAGGATACGGCGGCAACTGCGGAATGATCATCTCAACGTCAGGGTGATg                             >  2:499147/1‑123 (MQ=255)                                                                                                                   tGACGAACCGACAATGCCCAGCGAATCACCGCCATGTTCCAGGACAATGGATTCCAGCAGCTCTGCCGCGTCGGAAGGATACGGCGGCAACTGCGGAATGATCATCTCAACGTCAGGGTGATg                             <  1:499147/123‑1 (MQ=255)                                                                                                                         aCCGACAATGCCCAGCGAATCACCGCCATGTTCCAGGACAATGGATTCCAGCAGCTCTGCCGCGTCGGAAGGATACGGCGGCAACTGCGGAATGATCATCTCAACGTCAGGGTGATGTTCCGCCAGCCAGTTTTTTAACa      >  2:37034/1‑140 (MQ=255)                                                                                                                             aCAATGCCCAGCGAATCACCGCCATGTTCCAGGACAATGGATTCCAGCAGCTCTGCCGCGTCGGAAGGATACGGCGGCAACTGCGGAATGATCATCTCAACGTCAGGGTGATGTTCCGCCAGCCAGTTTTTTAACAAGCt  >  2:554775/1‑140 (MQ=255)                                                                                                                                  |                                                                                                                                       CTCGTTCTGACCGAGATAGTCCGTCAGCAGTTCAAACGGGCGCACCGCCGGGTTTACCACCACTGCGGGCAGCATAAAACATTGTGACAACCAGGTGGCGTAATATCCCCCCAGTGACGAACCGACAATACCCAGCGAATCACCGCCATGTTCCAGGACAATGGATTCCAGCAGCTCTGCCGCGTCGGAAGGATACGGCGGCAACTGCGGAATGATCATCTCAACGTCAGGGTGATGTTCCGCCAGCCAGTTTTTTAACAAGCT  >  NC_000913/3175686‑3175949

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