Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,417,219	G→T	E90* (GAA→TAA)	yjfP →	carboxylesterase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,417,219	0	G	T	100.0%	156.1 / NA	48	E90* (GAA→TAA)	yjfP	carboxylesterase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base T (22/26);  total (22/26)

CTATTTTGCCGTTGCGCTGGCGCAGGCTGGTTTGCGGGTGATCATGCCGGATGCGCCCGATCACGGTAGCCGTTTTAGTGGTGACGCAGCGCGGCGGTTAAATCAATTCTGGCAAATCTTGCTACAAAGTATGCAGGAATTCACTACTTTACGTGCGGCAATAGCCGAAGAAAACTGGCTGCTTGATGACCGTCTGGCAGTCGGTGGCGCGTCGATGGGCGCGATGACGGCACTGGGGATTACCGCTCGCCACCCCACGGTGAGATGTACCGC  >  NC_000913/4417083‑4417355                                                                                                                                         |                                                                                                                                         cTATTTTGCCGTTGCGCTGGCGCAGGCTGGTTTGCGGGTGATCATGCCGGATGCGCCCGATCACGGTAGCCGTTTTAGTGGTGACGCAGCGCGGCGGTTAAATCAATTCTGGCAAATCTTGCTACAAAGTATGCAGTAAt                                                                                                                                       >  2:293014/1‑140 (MQ=255) cTATTTTGCCGTTGCGCTGGCACAGGCTGGTTTGCGGGTGATCATGCCGGATGCGCCCGATCACGGTAGCCGTTTTAGTGGTGACGCAGCGCGGCGGTTAAATCAATTCTGGCAAATCTTGCTACAAAGTATGCAGTAAt                                                                                                                                       <  2:981844/140‑1 (MQ=255)  tATTTTGCCGTTGCGCTGGCGCAGGCTGGTTTGCGGGTGATCATGCCGGATGCGCCCGATCACGGTAGCCGTTTTAGTGGTGACGCAGCGCGGCGGTTAAATCAATTCTGGCAAATCTTGCTACAAAGTATGCAGTAAtt                                                                                                                                      <  2:136274/140‑1 (MQ=255)   aTTTTGCCGTTGCGCTGGCGCAGGCTGGTTTGCGGGTGATCATGCCGGATGCGCCCGATCACGGTAGCCGTTTTAGTGGTGACGCAGCGCGGCGGTTAAATCAATTCTGGCAAATCTTGCTACAAAGTATGCAGTAATTc                                                                                                                                     >  2:219180/1‑140 (MQ=255)           gTTGCGCTGGCGCAGGCTGGTTTGCGGGTGATCATGCCGGATGCGCCCGATCACGGTAGCCGTTTTAGTGGTGACGCAGCGCGGCGGTTAAAACAAATCCGGCGAATCTTTCTACCAAGTATGCTATAATTCACTACttt                                                                                                                             >  2:637833/1‑140 (MQ=255)                         ggCTGGTTTGCGGGTGATCATGCCGGATGCGCCCGATCACGGTAGCCGTTTTAGTGGTGACGCAGCGCGGCGGTTAAATCAATTCTGGCAAATCTTGCTACAAAGTATGCAGTAATTCACTACTTTACGTGCGGCAATAg                                                                                                               <  1:539571/140‑1 (MQ=255)                         ggCTGGTTTGCGGGTGATCATGCCGGATGCGCCCGATCACGGTAGCCGTTTTAGTGGTGACGCAGCGCGGCGGTTAAATCAATTCTGGCAAATCTTGCTACAAAGTATGCAGTAATTCACTACTTTACGTGCGGCAATAg                                                                                                               <  1:684773/140‑1 (MQ=255)                          gCTGGTTTGCGGGTGATCATGCCGGATGCGCTCGATCACGGTAGCCGTTTTAGTGGTGACGCAGCGCGGCGGTTAAATCAATTCTGGCAAATCTTGCTACAAAGTATGCAGTAATTCACTACTTTACGTGCGGCAATAGc                                                                                                              >  1:21916/1‑140 (MQ=255)                              gTTTGCGGGTGATCATGCCGGATGCGCCCGATCACGGTAGCCGTTTTAGTGGTGACGCAGCGCGGCGGTTAAATCAATTCTGGCAAATCTTGCTACAAAGTATGCAGTAATTCACTACTTTACGTGCGGCAATAGCCgaa                                                                                                          <  2:338976/140‑1 (MQ=255)                                    gggTGATCATGCCGGATGCGCCCGATCACGGTAGCCGTTTTAGTGGTGACGCAGCGCGGCGGTTTAATCAATTTTGGCAAATCTTGCTACAAAATATGCAGTAATTCACTCCTTTACGTGCGGCCATAGCCGAAGAAAAc                                                                                                    >  2:1227410/1‑140 (MQ=255)                                    gggTGATCATGCCGGATGCGCCCGATCACGGTAGCCGTTTTAGTGGTGACGCAGCGCGGCGGTTAAATCAATTCTGGCAAATCTTGCTACAAAGTATGCAGTAATTCACTACTTTACGTGCGGCAATAGCCGAAGAAAAc                                                                                                    >  2:186783/1‑140 (MQ=255)                                        gATCATGCCGGATGCGCCCGATCACGGTAGCCGTTTTAGTGGTGACGCAGCGCGGCGGTTAAATCAATTCTGGCAAATCTTGCTACAAAGTATGCAGTAATTCACTACTTTACGTGCGGCAATAGCCGAAGAAAACTGgc                                                                                                <  1:1227410/140‑1 (MQ=255)                                        gATCATGCCGGATGCGCCCGATCACGGTAGCCGTTTTAGTGGTGACGCAGCGCGGCGGTTAAATCAATTCTGGCAAATCTTGCTACAAAGTATGCAGTAATTCACTACTTTACGTGCGGCAATAGCCGAAGAAAACTGgc                                                                                                <  1:540406/140‑1 (MQ=255)                                         aTCATGCCGGATGCGCCCGATCACGGTAGCCGTTTTAGTGGTGACGCAGCGCGGCGGTTAAATTCAATCTGGCAAATCTTGCTACAAAATATGCAGTAATTTACTACTTTACCTGCGGCAATAGCCGAAGAAAACAGgct                                                                                               >  2:976320/1‑140 (MQ=255)                                             tGCCGGATGCGCCCGATCACGGTAGCCGTTTTAGTGGTGACGCAGCGCGGCGGTTAAATCAATTCTGGCAAATCTTGCTACAAAGTATGCAGTAATTCACTACTTTACGTGCGGCAATAGCCGAAGAAAACTGGCTGCtt                                                                                           <  2:540234/140‑1 (MQ=255)                                              gCCGGATGCGCCCGATCACGGTAGCCGTTTTAGTGGTGACGCAGCGCGGCGGTTAAATCAATTCTGGCAAATCTTGCTACAAAGTATGCAGTAATTCACTACTTTACGTGCGGCAATAGCCGAAGAAAACTGGCTGCTtg                                                                                          >  1:656083/1‑140 (MQ=255)                                                  gATGCGCCCGATCACGGTAGCCGTTTTAGTGGTGACGCAGCGCGGCGGTTAAATCAATTCTGGCAAATCTTGCTACAAAGTATGCAGTAATTCACTACTTTACGTGCGGCAATAGCCGAAGAAAACTGGCTGCTtgatga                                                                                      <  1:219180/140‑1 (MQ=255)                                                  gATGCGCCCGATCACGGTAGCCGTTTTAGTGGTGACGCAGCGCGGCGGTTAAATCAATTCTGGCAAATCTTGCTACAAAGTATGCAGTAATTCACTACTTTACGTGCGGCAATAGCCGAAGAAAACTGGCTGCTtgatga                                                                                      <  2:754114/140‑1 (MQ=255)                                                   aTGCGCCCGATCACGGTAGCCGTTTTAGTGGTGACGCAGCGCGGCGGTTAAATCAATTCTGGCAAAACTTGCTACAAAGTATGCAGTAATTCACTTCTTTACGTGCGGCAATAGCCGAAGAAAACTGGCTGCTTGTTGAc                                                                                     >  2:627814/1‑140 (MQ=255)                                                    tGCGCCCGATCACGGTAGCCGTTTTAGTGGTGACGCAGCGCGGCGGTTAAATCAATTCTGGCAAATCTTGCTACAAAGTATGCAGTAATTCACTACTTTACGTGCGGCAATAGCCGAAGAAAACTGGCTGCTTGATGAcc                                                                                    >  1:928814/1‑140 (MQ=255)                                                    tGCGCCCGATCACGGTAGCCGTTTTAGTGGTGACGCAGCGCGGCGGTTAAATCAATTCTGGCAAATCTTGCTACAAAGTATGCAGTAATTCACTACTTTACGTGCGGCAATAGCCGAAGAAAACTGGCTGCTTGATGAcc                                                                                    >  2:724557/1‑140 (MQ=255)                                                         ccGATCACGGTAGCCGTTTTAGTGGTGACGCAGCGCGGCGGTTAAATCAATTCTGGCAAATCTTGCTACAAAGTATGCAGTAATTCACTACTTTACGTGCGGCAATAGCCGAAGAAAACTGGCTGCTTGATGACCGTCTg                                                                               >  2:83972/1‑140 (MQ=255)                                                          cGATCACGGTAGCCGTTTTAGTGGTGACGCAGCGCGGCGGTTAAATCAATTCTGGCAAATCTTGCTACAAAGTATGCAGTAATTCACTACTTTACGTGCGGCAATAGCCGAAGAAAACTGGCTGCTTGATGACCGTCTgg                                                                              <  1:293014/140‑1 (MQ=255)                                                                   tAGCCGTTTTAGTGGTGACGCAGCGCGGCGGTTAAATCAATTCTGGCAAATCTTGCTACAAAGTATGCAGTAATTCACTACTTTACGTGCGGCAATAGCCGAAGAAAACTGGCTGCTTGATGACCGTCTGGCAGTCGGTg                                                                     >  1:549712/1‑140 (MQ=255)                                                                   tAGCCGTTTTAGTGGTGACGCAGCGCGGCGGTTAAATCAATTCTGGCAAATCTTGCTACAAAGTATGCAGTAATTCACTACTTTACGTGCGGCAATAGCCGAAGAAAACTGGCTGCTTGATGACCGTCTGGCAGTCGGTg                                                                     <  2:656083/140‑1 (MQ=255)                                                                     gCCGTTTTAGTGGTGACGCAGCGCGGCGGTTAAATCAATTCTGGCAAATCTTGCTACAAAGTATGCAGTAATTCACTACTTTACGTGCGGCAATAGCCGAAGAAAACTGGCTGCTTGATGACCGTCTGGCAGTCGGTGgc                                                                   >  2:1037215/1‑140 (MQ=255)                                                                     gCCGTTTTAGTGGTGACGCAGCGCGGCGGTTAAATCAATTCTGGCAAATCTTGCTACAAAGTATGCAGTAATTCACTACTTTACGTGCGGCAATAGCCGAAGAAAACTGGCTGCTTGATGACCGTCTGGCAGTCGGTGgc                                                                   >  1:88714/1‑140 (MQ=255)                                                                        gTTTTAGTGGTGACGCAGCGCGGCGGTTAAATCAATTCTGGCAAATCTTGCTACAAAGTATGCAGTAATTCACTACTTTACGTGCGGCAATAGCCGAAGAAAACTGGCTGCTTGATGACCGTCTGGCAGTCGGTGGcgcg                                                                >  2:292245/1‑140 (MQ=255)                                                                                ggtgACGCAGCGCGGCGGTTAAATCAATTCTGGCAAATCTTGCTACAAAGTATGCAGTAATTCACTACTTTACGTGCGGCAATAGCCGAAGAAAACTGGCTGCTTGATGACCGTCTGGCAGTCGGTGGCGCGTCGATggg                                                        <  2:21916/140‑1 (MQ=255)                                                                                 gtgACGCAGCGCGGCGGTTAAATCAATTCTGGCAAATCTTGCTACAAAGTATGCAGTAATTCACTACTTTACGTGCGGCAATAGCCGAAGAAAACTGGCTGCTTGATGACCGTCTGGCAGTCGGTGGCGCGTCGATGGgc                                                       <  1:903465/140‑1 (MQ=255)                                                                                 gtgACGCAGCGCGGCGGTTAAATCAATTCTGGCAAATCTTGCTACAAAGTATGCAGTAATTCACTACTTTACGTGCGGCAATAGCCGAAGAAAACTGGCTGCTTGATGACCGTCTGGCAGTCGGTGGCGCGTCGATGGgc                                                       <  1:693190/140‑1 (MQ=255)                                                                                    aCGCAGCGCGGCGGTTAAATCAATTCTGGCAAATCTTGCTACAAAGTATGCAGTAATTCACTACTTTACGTGCGGCAATAGCCGAAGAAAACTGGCTGCTTGATGACCGTCTGGCAGTCGGTGGCGCGTCGATGGGcgcg                                                    >  2:293088/1‑140 (MQ=255)                                                                                    aCGCAGCGCGGCGGTTAAATCAATTCTGGCAAATCTTGCTACAAAGTATGCAGTAATTCACTACTTTACGTGCGGCAATAGCCGAAGAAAACTGGCTGCTTGATGACCGTCTGGCAGTCGGTGGCGCGTCGATGGGcgcg                                                    >  2:999777/1‑140 (MQ=255)                                                                                        aGCGCGGCGGTTAAATCAATTTTGGCAAATCCTGCTACAAAGGAATGAGAAATTCACCACCTTACCGGCGGCCATAGCCGAAGAAAACTGGCTGCCTGTCGACCATCTGGCAGTCGGTGGCGCGGCCATGGGCGCGAAGa                                                >  2:1257113/1‑140 (MQ=255)                                                                                        aGCGCGGCGGTTAAATCAATTCTGGCAAATCTTGCTACAAAGTATGCAGTAATTCACTACTTTACGTGCGGCAATAGCCGAAGAAAACTGGCTGCTTGATGACCGTCTGGCAGTCGGTGGCGCGTCGATGGGCGCGATGa                                                <  1:83972/140‑1 (MQ=255)                                                                                         gcgcGGCGGTTAAATCAATTCTGGCAAATCTTGCTACAAAGTATGCAGTAATTCACTACTTTACGTGCGGCAATAGCCGAAGAAAACTGGCTGCTTGATGACCGTCTGGCAGTCGGTGGCGCGTCGATGGGCGCGATGa                                                >  1:911358/1‑139 (MQ=255)                                                                                         gcgcGGCGGTTAAATCAATTCTGGCAAATCTTGCTACAAAGTATGCAGTAATTCACTACTTTACGTGCGGCAATAGCCGAAGAAAACTGGCTGCTTGATGACCGTCTGGCAGTCGGTGGCGCGTCGATGGGCGCGATGa                                                <  2:911358/139‑1 (MQ=255)                                                                                          cgcgGCGGTTAAATCAATTCTGGCAAATCTTGCTACAAAGTATGCAGTAATTCACTACTTTACGTGCGGCAATAGCCGAAGAAAACTGGCTGCTTGATGACCGTCTGGCAGTCGGTGGCGCGTCGATGGGCGCGATGACg                                              >  2:183731/1‑140 (MQ=255)                                                                                             ggcggTTAAATCAATTCTGGCAAATCTTGCTACAAAGTATGCAGTAATTCACTACTTTACGTGCGGCAATAGCCGAAGAAAACTGGCTGCTTGATGACCGTCTGGCAGTCGGTGGCGCGTCGATGGGCGCGATGACGGCa                                           <  2:449164/140‑1 (MQ=255)                                                                                             ggcggTTAAATCAATTCTGGCAAATCTTGCTACAAAGTATGCAGTAATTCACTACTTTACGTGCGGCAATAGCCGAAGAAAACTGGCTGCTTGATGACCGTCTGGCAGTCGGTGGCGCGTCGATGGGCGCGATGACGGCa                                           <  2:549712/140‑1 (MQ=255)                                                                                                           ttCTGGCAAATCTTGCTACAAAGTATGCAGTAATTCACTACTTTACGTGCGGCAATAGCCGAAGAAAACTGGCTGCTTGATGACCGTCTGGCAGTCGGTGGCGCGTCGATGGGCGCGATGACGGCACTGGGGATTACCGc                             >  2:606528/1‑140 (MQ=255)                                                                                                            tCTGGCAAATCTTGCTACAAAGTATGCAGTAATTCACTACTTTACGTGCGGCAATAGCCGAAGAAAACTGGCTGCTTGATGACCGTCTGGCAGTCGGTGGCGCGTCGATGGGCGCGATGACGGCACTGGGGATTACCGCt                            >  1:173557/1‑140 (MQ=255)                                                                                                             cTGGCAAATCTTGCTACAAAGTATGCAGTAATTCACTACTTTACGTGCGGCAATAGCCGAAGAAAACTGGCTGCTTGATGACCGTCTGGCAGTCGGTGGCGCGTCGATGGGCGCGATGACGGCACTGGGGATTACCGCTc                           <  2:928814/140‑1 (MQ=255)                                                                                                              tGGCAAATCTTGCTACAAAGTATGCAGTAATTCACTACTTTACGTGCGGCAATAGCCGAAGAAAACTGGCTGCTTGATGACCGTCTGGCAGTCGGTGGCGCGTCGATGGGCGCGATGACGGCACTGGGGATTACCGCTCg                          <  1:606528/140‑1 (MQ=255)                                                                                                               ggCAAATCTTGCTACAAAGTATGCAGTAATTCACTACTTTACGTGCGGCAATAGCCGAAGAAAACTGGCTGCTTGATGACCGTCTGGCAGTCGGTGGCGCGTCGATGGGCGCGATGACGGCACTGGGGATTACCGCTCGc                         >  1:151477/1‑140 (MQ=255)                                                                                                                  aaaTCTTGCTACAAAGTATGCAGTAATTCACTACTTTACGTGCGGCAATAGCCGAAGAAAACTGGCTGCTTGATGACCGTCTGGCAGTCGGTGGCGCGTCGATGGGCGCGATGACGGCACTGGGGATTACCGCTCGCCAc                      >  1:633320/1‑140 (MQ=255)                                                                                                                   aaTCTTGCTACAAAGTATGCAGTAATTCACTACTTTACGTGCGGCAATAGCCGAAGAAAACTGGCTGCTTGATGACCGTCTGGCAGTCGGTGGCGCGTCGATGGGCGCGATGACGGCACTGGGGATTACCGCTCGCCAcc                     <  1:293088/140‑1 (MQ=255)                                                                                                                    aTCTTGCTACAAAGTATGCAGTAATTCACTACTTTACGTGCGGCAATAGCCGAAGAAAACTGGCTGCTTGATGACCGTCTGGCAGTCGGTGGCGCGTCGATGGGCGCGATGACGGCACTGGGGATTACCGCTCGCCAccc                    >  2:780009/1‑140 (MQ=255)                                                                                                                     tCTTGCTACAAAGTATGCAGTAATTCACTACTTTACGTGCGGCAATAGCCGAAGAAAACTGGCTGCTTGATGACCGTCTGGCAGTCGGTGGCGCGTCGATGGGCGCGATGACGGCACTGGGGATTACCGCTCGCCAcccc                   >  2:371708/1‑140 (MQ=255)                                                                                                                      cTTGCTACAAAGTATGCAGTAATTCACTACTTTACGTGCGGCAATAGCCGAAGAAAACTGGCTGCTTGATGACCGTCTGGCAGTCGGTGGCGCGTCGATGGGCGCGATGACGGCACTGGGGATTACCGCTCGCCACCCCa                  <  1:1275375/140‑1 (MQ=255)                                                                                                                        tGCTACAAAGTATGCAGTAATTCACTACTTTACGTGCGGCAATAGCCGAAGAAAACTGGCTGCTTGATGACCGTCTGGCAGTCGGTGGCGCGTCGATGGGCGCGATGACGGCACTGGGGATTACCGCTCGCCACCCCACg                <  1:467378/140‑1 (MQ=255)                                                                                                                        tGCTACAAAGTATGCAGTAATTCACTACTTTACGTGCGGCAATAGCCGAAGAAAACTGGCTGCTTGATGACCGTCTGGCAGTCGGTGGCGCGTCGATGGGCGCGATGACGGCACTGGGGATTACCGCTCGCCACCCCACg                <  1:292245/140‑1 (MQ=255)                                                                                                                                   aTGCAGTAATTCACTACTTTACGTGCGGCAATAGCCGAAGAAAACTGGCTGCTTGATGACCGTCTGGCAGTCGGTGGCGCGTCGATGGGCGCGATGACGGCACTGGGGATTACCGCTCGCCACCCCACGGTGAGATGTAc     >  2:468301/1‑140 (MQ=255)                                                                                                                                      cAGTAATTCACTACTTTACGTGCGGCAATAGCCGAAGAAAACTGGCTGCTTGATGACCGTCTGGCAGTCGGTGGCGCGTCGATGGGCGCGATGACGGCACTGGGGATTACCGCTCGCCACCCCACGGTGAGATGTACCGc  <  1:183731/140‑1 (MQ=255)                                                                                                                                         |                                                                                                                                         CTATTTTGCCGTTGCGCTGGCGCAGGCTGGTTTGCGGGTGATCATGCCGGATGCGCCCGATCACGGTAGCCGTTTTAGTGGTGACGCAGCGCGGCGGTTAAATCAATTCTGGCAAATCTTGCTACAAAGTATGCAGGAATTCACTACTTTACGTGCGGCAATAGCCGAAGAAAACTGGCTGCTTGATGACCGTCTGGCAGTCGGTGGCGCGTCGATGGGCGCGATGACGGCACTGGGGATTACCGCTCGCCACCCCACGGTGAGATGTACCGC  >  NC_000913/4417083‑4417355

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