Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,377,189	C→T	H105H (CAC→CAT)  *178* (TAG→TAA)	sugE → blc ←	multidrug efflux system protein outer membrane lipoprotein cell division and growth lipocalin

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,377,189	0	C	T	100.0%	180.0 / NA	56	H105H (CAC→CAT)  *178* (TAG→TAA)	sugE blc	multidrug efflux system protein outer membrane lipoprotein cell division and growth lipocalin
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (29/27);  total (29/27)

GTGGACGGGTATTGGCGCAGTCGGCGCGGCCATAACCGGCATTGTGCTGCTCGGTGAGTCCGCTAACCCGATGCGCCTGGCGAGTCTGGCGTTAATCGTATTGGGGATTATTGGTCTGAAACTCAGCACTCACTAACTACCAGGCTGCTGTACCCAAATAAATTTACTGACATCAAACCCTTCCCGGGTCGCGACTGCCAGCATCTCCTGTTTCACTTCGTCAGAAATGGTTGGCGTGCGGGAGAGTATCCACAGGTAGTCGCGG  >  NC_000913/4377057‑4377321                                                                                                                                     |                                                                                                                                     gtgGACGGGTATTGGCGCAGTCGGCGCGGCCATAACCGGCATTGTGCTGCTCGGTGAGTCCGCTAACCCGATGCGCCTGGCGAGTCTGGCGTTAATCGTATTGGGGATTATTGGTCTGAAACTCAGCACTCattaacta                                                                                                                                >  1:1190429/1‑139 (MQ=255)    gACGGGTATTGGCGCAGTCGGCGCGGCCATAACCGGCATTGTGCTGCTCGGTGAGTCCGCTAACCCGATGCGCCTGGCGAGTCTGGCGTTAATCGTATTGGGGATTATTGGTCTGAAACTCAGCACTCATTAACTACCa                                                                                                                             <  1:1044102/139‑1 (MQ=255)              ggCGCAGTCGGCGCGGCCATAACCGGCATTGTGCTGCTCGGTGAGTCCGCTAACCCGATGCGCCTGGCGAGTCTGGCGTTAATCGTATTGGGGATTATTGGTCTGAAACTCAGCACTCATTAACTACCAGGCTGCTGTa                                                                                                                   >  2:100292/1‑139 (MQ=255)              ggCGCAGTCGGCGCGGCCATAACCGGCATTGTGCTGCTCGGTGAGTCCGCTAACCCGATGCGCCTGGCGAGTCTGGCGTTAATCGTATTGGGGATTATTGGTCTGAAACTCAGCACTCATTAACTACCAGGCTGCTGTa                                                                                                                   <  2:676789/139‑1 (MQ=255)                    gTCGGCGCGGCCATAACCGGCATTGTGCTGCTCGGTGAGTCCGCTAACCCGATGCGCCTGGCGAGTCTGGCGTTAATCGTATTGGGGATTATTGGTCTGAAACTCAGCACTCATTAACTACCAGGCTACTGTACCCaaa                                                                                                             >  2:330814/1‑139 (MQ=255)                              ccATAACCGGCATTGTGCTGCTCGGTGAGTCCGCTAACCCGATGCGCCTGGCGAGTCTGGCGTTAATCGTATTGGGGATTATTGGTCTGAAACTCAGCACTCATTAACTACCAGGCTGCTGTACCCAAATAAATTTACt                                                                                                   <  2:191388/139‑1 (MQ=255)                               cATAACCGGCATTGTGCTGCTCGGTGAGTCCGCTAACCCGATGCGCCTGGCGAGTCTGGCGTTAATCGTATTGGGGATTATTGGTCTGAAACTCAGCACTCATTAACTACCAgg                                                                                                                           >  1:1029173/1‑114 (MQ=255)                               cATAACCGGCATTGTGCTGCTCGGTGAGTCCGCTAACCCGATGCGCCTGGCGAGTCTGGCGTTAATCGTATTGGGGATTATTGGTCTGAAACTCAGCACTCATTAACTACCAgg                                                                                                                           <  2:1029173/114‑1 (MQ=255)                               cATAACCGGCATTGTGCTGCTCGGTGAGTCCGCTAACCCGATGCGCCTGGCGAGTCTGGCGTTAATCGTATTGGGGATTATTGGTCTGAAACTCAGCACTCATTAACTACCAGGCTGCTGTACCCAAATAAATTTACTg                                                                                                  <  2:863458/139‑1 (MQ=255)                                    ccGGCATTGTGCTGCTCGGTGAGTCCGCTAACCCGATGCGCCTGGCGAGTCTGGCGTTAATCGTATTGGGGATTATTGGTCTGAAACTCAGCACTCATTAACTACCAGGCTGCTGTACCCAAATAAATTTACTGACATc                                                                                             <  2:1277337/139‑1 (MQ=255)                                        cATTGTGCTGCTCGGTGAGTCCGCTAACCCGATGCGCCTGGCGAGTCTGGCGTTAATCGTATTGGGGATTATTGGTCTGAAACTCAGCACTCATTAACTACCAGGCTGCTGTACCCAAATAAATTTACTGACATCAAAc                                                                                         >  2:1358765/1‑139 (MQ=255)                                        cATTGTGCTGCTCGGTGAGTCCGCTAACCCGATGCGCCTGGCGAGTCTGGCGTTAATCGTATTGGGGATTATTGGTCTGAAACTCAGCACTCATTAACTACCAGGCTGCTGTACCCAAATAAATTTACTGACATCAAAc                                                                                         >  2:93894/1‑139 (MQ=255)                                            gtgCTGCTCGGTGAGTCCGCTAACCCGATGCGCCTGGCGAGTCTGGCGTTAATCGTATTGGGGATTATTGGTCTGAAACTCAGCACTCATTAACTACCAGGCTGCTGTa                                                                                                                   <  2:644993/109‑1 (MQ=255)                                            gtgCTGCTCGGTGAGTCCGCTAACCCGATGCGCCTGGCGAGTCTGGCGTTAATCGTATTGGGGATTATTGGTCTGAAACTCAGCACTCATTAACTACCAGGCTGCTGTa                                                                                                                   >  1:644993/1‑109 (MQ=255)                                             tgctgcTCGGTGAGTCCGCTAACCCGATGCGCCTGGCGAGTCTGGCGTTAATCGTATTGGGGATTATTGGTCTGAAACTCAGCACTCATTAACTACCAGGCTGCTGTACCCAAATAAATTTACTGACATCAAACCCTTc                                                                                    >  1:1365791/1‑139 (MQ=255)                                                       tGAGTCCGCTAACCCGATGCGCCTGGCGAGTCTGGCGTTAATCGTATTGGGGATTATTGGTCTGAAACTCAGCACTCATTAACTACCAgg                                                                                                                           <  1:302600/90‑1 (MQ=255)                                                       tGAGTCCGCTAACCCGATGCGCCTGGCGAGTCTGGCGTTAATCGTATTGGGGATTATTGGTCTGAAACTCAGCACTCATTAACTACCAgg                                                                                                                           >  2:302600/1‑90 (MQ=255)                                                       tGAGTCCGCTAACCCGATGCGCCTGGCGAGTCTGGCGTTAATCGTATTGGGGATTATTGGTCTGAAACTCAGCACTCATTAACTACCAGGCTGCTGTACCCAAATAAATTTACTGACATCAAACCCTTCCCGGGTcgcg                                                                          >  2:112061/1‑139 (MQ=255)                                                           tCCGCTAACCCGATGCGCCTGGCGAGTCTGGCGTTAATCGTATTGGGGATTATTGGTCTGAAACTCAGCACTCATTAACTACCAGGCTGCTGTACCCAAATAAATTTACTGACATCAAACCCTTCCCGGGTCGCGACTg                                                                      <  1:100292/139‑1 (MQ=255)                                                            ccGCTAACCCGATGCGCCTGGCGAGTCTGGCGTTAATCGTATTGGGGATTATTGGTCTGAAACTCAGCACTCATTAACTACCAGGCTGCTGTACCCAAATAAATTTACTGACATCAAACCCTTCCCGGGTCGCGACTGc                                                                     >  2:28827/1‑139 (MQ=255)                                                            ccGCTAACCCGATGCGCCTGGCGAGTCTGGCGTTAATCGTATTGGGGATTATTGGTCTGAAACTCAGCACTCATTAACTACCAGGCTGCTGTACCCAAATAAATTTACTGACATCAAACCCTTCCCGGGTCGCGACTGc                                                                     >  1:164448/1‑139 (MQ=255)                                                            ccGCTAACCCGATGCGCCTGGCGAGTCTGGCGTTAATCGTATTGGGGATTATTGGTCTGAAACTCAGCACTCATTAACTACCAGGCTGCTGTACCCAAATAAATTTACTGACATCAAACCCTTCCCGGGTCGCGACTGc                                                                     >  1:826919/1‑139 (MQ=255)                                                               cTAACCCGATGCGCCTGGCGAGTCTGGCGTTAATCGTATTGGGGATTATTGGTCTGAAACTCAGCACTCATTAACTACCAGGCTGCTGTACCCAAATAAATTTACTGACATCAAACCCTTCCCGGGTCGCGACTGCCAg                                                                  >  2:734902/1‑139 (MQ=255)                                                               cTAACCCGATGCGCCTGGCGAGTCTGGCGTTAATCGTATTGGGGATTATTGGTCTGAAACTCAGCACTCATTAACTACCAGGCTGCTGTACCCAAATAAATTTACTGACATCAAACCCTTCCCGGGTCGCGACTGCCAg                                                                  >  2:1019914/1‑139 (MQ=255)                                                                 aaCCCGATGCGCCTGGCGAGTCTGGCGTTAATCGTATTGGGGATTATTGGTCTGAAACTCAGCACTCATTAACTACCAGGCTGCTGTACCCAAATAAATTTACTGACATCAAACCCTTCCCGGGTCGCGACTGCCAGCa                                                                >  1:293718/1‑139 (MQ=255)                                                                   cccGATGCGCCTGGCGAGTCTGGCGTTAATCGTATTGGGGATTATTGGTCTGAAACTCAGCACTCATTAACTACCAGGCTGCTGTACCCAAATAAATTTACTGACATCAAACCCTTCCCGGGTCGCGACTGCCAGCAtc                                                              <  1:1019914/139‑1 (MQ=255)                                                                   cccGATGCGCCTGGCGAGTCTGGCGTTAATCGTATTGGGGATTATTGGTCTGAAACTCAGCACTCATTAACTACCAGGCTGCTGTACCCAAATAAATTTACTGACATCAAACCCTTCCCGGGTCGCGACTGCCAGCAtc                                                              <  1:734902/139‑1 (MQ=255)                                                                     cGATGCGCCTGGCGAGTCTGGCGTTAATCGTATTGGGGATTATTGGTCTGAAACTCAGCACTCATTAACTACCAGGCTGCTGTACCCAAATAAATTTACTGACATCAAACCCTTCCCGGGTCGCGACTGCCAGCAtctc                                                            >  1:742/1‑139 (MQ=255)                                                                       aTGCGCCTGGCGAGTCTGGCGTTAATCGTATTGGGGATTATTGGTCTGAAACTCAGCACTCATTAACTACCAGGCTGCTGTACCCAAATAAATTTACTGACATCAAACCCTTc                                                                                    <  1:211990/113‑1 (MQ=255)                                                                       aTGCGCCTGGCGAGTCTGGCGTTAATCGTATTGGGGATTATTGGTCTGAAACTCAGCACTCATTAACTACCAGGCTGCTGTACCCAAATAAATTTACTGACATCAAACCCTTc                                                                                    >  2:211990/1‑113 (MQ=255)                                                                              tGGCGAGTCTGGCGTTAATCGTATTGGGGATTATTGGTCTGAAACTCAGCACTCATTAACTACCAGGCTGCTGTACCCAAATAAATTTACTGACATCAAACCCTTCCCGGGTCGCGACTGCCAGCATCTCCTGTTTCAc                                                   <  2:293718/139‑1 (MQ=255)                                                                                       tGGCGTTAATCGTATTGGGGATTATTGGTCTGAAACTCAGCACTCATTAACTACCAGGCTGCTGTACCCAAATAAATTTACTGACATCAAACCCTTCCCGGGTCGCGACTGCCAGCATCTCCTGTTTCACTTCGTCAGa                                          <  1:112061/139‑1 (MQ=255)                                                                                       tGGCGTTAATCGTATTGGGGATTATTGGTCTGAAACTCAGCACTCATTAACTACCAGGCTGCTGTACCCAAATAAATTTACTGACATCAAACCCTTCCCGGGTCGCGACTGCCAGCATCTCCTGTTTCACTTCGTCAGa                                          <  2:742/139‑1 (MQ=255)                                                                                         gCGTTAATCGTATTGGGGATTATTGGTCTGAAACTCAGCACTCATTAACTACCAGGCTGCTGTACCCAAATAAATTTACTGACATCAAACCCTTCCCGGGTcg                                                                            <  1:1239622/103‑1 (MQ=255)                                                                                         gCGTTAATCGTATTGGGGATTATTGGTCTGAAACTCAGCACTCATTAACTACCAGGCTGCTGTACCCAAATAAATTTACTGACATCAAACCCTTCCCGGGTcg                                                                            >  2:1239622/1‑103 (MQ=255)                                                                                              aaTCGTATTGGGGATTATTGGTCTGAAACTCAGCACTCATTAACTACCAGgctg                                                                                                                        <  2:970258/54‑1 (MQ=255)                                                                                              aaTCGTATTGGGGATTATTGGTCTGAAACTCAGCACTCATTAACTACCAGGCTGCTGTACCCa                                                                                                               <  2:774817/63‑1 (MQ=255)                                                                                              aaTCGTATTGGGGATTATTGGTCTGAAACTCAGCACTCATTAACTACCAGGCTGCTGTACCCa                                                                                                               >  1:774817/1‑63 (MQ=255)                                                                                              aaTCGTATTGGGGATTATTGGTCTGAAACTCAGCACTCATTAACTACCAGGCTGCTGTACCCAAATAAATTTACTGACATCAAACCCTTCCCGGGTCGCGACTGCCAGCATCTCCTGTTTCACTTCGTCAGAAATGGtt                                   >  1:993310/1‑139 (MQ=255)                                                                                              aaTCGTATTGGGGATTATTGGTCTGAAACTCAGCACTCATTAACAACCAGgctg                                                                                                                        >  1:970258/1‑54 (MQ=255)                                                                                                        gggATTATTGGTCTGAAACTCAGCACTCATTAACTACCAGGCTGCTGTACCCAAATAAATTTACTGACATCAAACCCTTCCCGGGTCGCGACTGCAAGCATCTCCTGTTTCACTTCGTCAGAAATGGTTGGCGTGCggg                         >  1:846677/1‑139 (MQ=255)                                                                                                         ggATTATTGGTCTGAAACTCAGCACTCATTAACTACCAGGCTGCTGTACCCAAATAAATTTACTGACATCAAACCCTTCCCGGGTCGCGACTGCCAGCATCTCCTGTTTCACTTCGTCAGAAATGGTTGGCGTGCGGga                        <  2:1181329/139‑1 (MQ=255)                                                                                                         ggATTATTGGTCTGAAACTCAGCACTCATTAACTACCAGGCTGCTGTACCCAAATAAATTTACTGACATCAAACCCTTCCCGGGTCGCGACTGCCAGCATCTCCTGTTTCACTTCGTCAGAAATGGTTGGCGTGCGGga                        >  1:1181329/1‑139 (MQ=255)                                                                                                           attattGGTCTGAAACTCAGCACTCATTAACTACCAGGCTGCTGTACCCAAATAAATTTACTGACATCAAACCCTTCCCGGGTCGCGACTGCCAGCATCTCCTGTTTCACTTCGTCAGAAATGGTTGGCGTGCGGgaga                      <  2:164448/139‑1 (MQ=255)                                                                                                            ttattGGTCTGAAACTCAGCACTCATTAACTACCAGGCTGCTGTACCCAAATAAATTTACTGACATCAAACCCTTCCCGGGTCGCGACTGCCAGCATCTCCTGTTTCACTTCGTCAGAAATGGTTGGCGTGCGGGagag                     <  2:667074/139‑1 (MQ=255)                                                                                                             tattGGTCTGAAACTCCGCACTCATTAACTACCAGGCTGCTGTACCCAAATAAATTTACTGACATCAAACCCTTCCCGGGTCGCGACTGCCAGCATCTCCTGTTTCCCTTCGTCAGAAATGGTTGGCGTGCGGGAGAGt                    >  1:997700/1‑139 (MQ=255)                                                                                                                tGGTCTGAAACTCAGCACTCATTAACTACCAGGCTGCTGTACCCAAATAAATTTACTGACATCAAACCCTTCCCGGGTCGCGACTGCCAGCATCTCCTGTTTCACTTCGTCAGAAATGGTTGGCGTGCGGGAGAGTATc                 <  2:427498/139‑1 (MQ=255)                                                                                                                  gTCTGAAACTCAGCACTCATTAACTACCAGGCTGCTGTACCCAAATAAATTTACTGACATCAAACCCTTCCCGGGTCGCGACTGCCAGCATCTCCTGTTTCACTTCGTCAGAAATGGTTGGCGTGCGGGAGAGTATCca               >  1:486850/1‑139 (MQ=255)                                                                                                                  gTCTGAAACTCAGCACTCATTAACTACCAGGCTGCTGTACCCAAATAAATTTACTGACATCAAACCCTTCCCGGGTCGCGACTGCCAGCATATCCTGTTTCACTTCGTCAGAAATGGTTg                                  >  1:549039/1‑120 (MQ=255)                                                                                                                  gTCTGAAACTCAGCACTCATTAACTACCAGGCTGCTGTACCCAAATAAATTTACTGACATCAAACCCTTCCCGGGTCGCGACTGCCAGCATATCCTGTTTCACTTCGTCAGAAATGGTTg                                  <  2:549039/120‑1 (MQ=255)                                                                                                                      gAAACTCAGCACTCATTAACTACCAGGCTGCTGTACCCAAATAAATTTACTGACATCAAACCCTTCCCGGGTCGCGACTGCCAGCATCTCCTGTTTCACTTCGTCAGAAATGGTTGGCGTGCGGGAGAGTATCCACAgg           <  1:93894/139‑1 (MQ=255)                                                                                                                      gAAACTCAGCACTCATTAACTACCAGGCTGCTGTACCCAAATAAATTTACTGACATCAAACCCTTCCCGGGTCGCGACTGCCAGCATCTCCTGTTTCACTTCGTCAGAAATGGTTGGCGTGCGGGAGAGTATCCACAgg           <  2:1190429/139‑1 (MQ=255)                                                                                                                       aaaCTCAGCACTCATTAACTACCAGGCTGCTGTACCCAAATAAATTTACTGACATCAAACCCTTCCCGGGTCGCGACTGCCAGCATCTCCTGTTTCACTTCGTCAGAAATGGTTGGCGTGCGGgaga                      >  1:608437/1‑127 (MQ=255)                                                                                                                       aaaCTCAGCACTCATTAACTACCAGGCTGCTGTACCCAAATAAATTTACTGACATCAAACCCTTCCCGGGTCGCGACTGCCAGCATCTCCTGTTTCACTTCGTCAGAAATGGTTGGCGTGCGGgaga                      <  2:608437/127‑1 (MQ=255)                                                                                                                             aGCACTTATTAACTACCAGGCTGCTGTACCCAAATAAATTTACTg                                                                                                  >  2:102875/1‑45 (MQ=38)                                                                                                                             aGCACTTATTAACTACCAGGCTGCTGTACCCAAATAAATTTACTg                                                                                                  <  1:102875/45‑1 (MQ=38)                                                                                                                             aGCACTCATTAACTACCAGGCTGCTGTACCCAAATAAATTTACTGACATCAAACCCTTCCCGGGTCGCGACTGCCAGCATCTCCTGTTTCACTTCGTCAGAAATGGTTGGCGTGCGGGAGAGTATCCACAGGTAGTcgc    >  1:933539/1‑139 (MQ=255)                                                                                                                               cactcattAACTACCAGGCTGCTGTACCCAAATAAATTTACTGACATCAAACCCTTCCCGGGTCGCGACTGCCAGCATCTCCTGTTTCACTTCGTCAGAAATGGTTGGCGTGCgg                          <  2:62290/115‑1 (MQ=255)                                                                                                                               cactcattAACTACCAGGCTGCTGTACCCAAATAAATTTACTGACATCAAACCCTTCCCGGGTCGCGACTGCCAGCATCTCCTGTTTCACTTCGTCAGAAATGGTTGGCGTGCgg                          >  1:62290/1‑115 (MQ=255)                                                                                                                               cactcattAACTACCAGGCTGCTGTACCCAAATAAATTTACTGACATCAAACCCTTCCCGGGTCGCGACTGCCAGCATCTCCTGTTTCACTTCGTCAGAAATGGTTGGCGTGCGGGAGAGTATCCACAGGTAGTCGCgg  >  1:1138710/1‑139 (MQ=255)                                                                                                                                     |                                                                                                                                     GTGGACGGGTATTGGCGCAGTCGGCGCGGCCATAACCGGCATTGTGCTGCTCGGTGAGTCCGCTAACCCGATGCGCCTGGCGAGTCTGGCGTTAATCGTATTGGGGATTATTGGTCTGAAACTCAGCACTCACTAACTACCAGGCTGCTGTACCCAAATAAATTTACTGACATCAAACCCTTCCCGGGTCGCGACTGCCAGCATCTCCTGTTTCACTTCGTCAGAAATGGTTGGCGTGCGGGAGAGTATCCACAGGTAGTCGCGG  >  NC_000913/4377057‑4377321

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