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%	154.1 / NA	48	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 (23/25);  total (23/25)

TGTGGACGGGTATTGGCGCAGTCGGCGCGGCCATAACCGGCATTGTGCTGCTCGGTGAGTCCGCTAACCCGATGCGCCTGGCGAGTCTGGCGTTAATCGTATTGGGGATTATTGGTCTGAAACTCAGCACTCACTAACTACCAGGCTGCTGTACCCAAATAAATTTACTGACATCAAACCCTTCCCGGGTCGCGACTGCCAGCATCTCCTGTTTCACTTCGTCAGAAATGGTTGGCGTGCGGGAGAGTATCCACAGGTAGTCGCGG  >  NC_000913/4377056‑4377321                                                                                                                                      |                                                                                                                                     tgtgGACGGGTATTGGCGCAGTCGGCGCGGCCATAACCGGCATTGTGCTGCTCGGTGAGTCCGCTAACCCGATGCGCCTGGCGAGTCTGGCGTTAATCGTATTGGGGATTATTGGTCTGAAACTCAGCACTCattaact                                                                                                                                 >  2:1464278/1‑139 (MQ=255)  gtgGACGGGTATTGGCGCAGTCGGCGCGGCCATAACCGGCATTGTGCTGCTCGGTGAGTCCGCTAACCCGATGCGCCTGGCGAGTCTGGCGTTAATCGTATTGGGGATTATTGGTCTGAAACTCAGCACTCattaacta                                                                                                                                <  2:1187487/139‑1 (MQ=255)     gACGGGTATTGGCGCAGTCGGCGCGGCCATAACCGGCATTGTGCTGCTCGGTGAGTCCGCTAACCCGATGCGCCTGGCGAGTCTGGCGTTAATCGTATTGGGGATTATTGGTCTGAAACTCAGCACTCATTAACTACCa                                                                                                                             <  2:1651190/139‑1 (MQ=255)       cGGGTATTGGCGCAGTCGGCGCGGCCATAACCGGCATTGTGCTGCTCGGTGAGTCCGCTAACCCGATGCGCCTGGCGAGTCTGGCGTTAATCGTATTGGGGATTATTGGTCTGAAACTCAGCACTCATTAACTACCAgg                                                                                                                           <  2:646451/139‑1 (MQ=255)                               ccATAACCGGCATTGTGCTGCTCGGTGAGTCCGCTAACCCGATGCGCCTGGCGAGTCTGGCGTTAATCGTATTGGGGATTATTGGTCTGAAACTCAGCACTCATTAACTACCAGGCTGCTGTACCCAAATAAATTTACt                                                                                                   <  2:1385113/139‑1 (MQ=255)                                         cATTGTGCTGCTCGGTGAGTCCGCTAACCCGATGCGCCTGGCGAGTCTGGCGTTAATCGTATTGGGGATTATTGGTCTGAAACTCAGCACTCATTAACTACCAGGCTGCTGTACCCAAATAAATTTACTGACATCAAAc                                                                                         >  2:1343337/1‑139 (MQ=255)                                            tgtgCTGCTCGGTGAGTCCGCTAACCCGATGCGCCTGGCGAGTCTGGCGTTAATCGTATTGGGGATTATTGGTCTGAAACTCAGCACTCATTAACTACCAGGCTGCTGTACCCAAATAAATTTACTGACATCAAACCCt                                                                                      >  1:495613/1‑139 (MQ=255)                                              tgctgcTCGGTGAGTCCGCTAACCCGATGCGCCTGGCGAGTCTGGCGTTAATCGTATTGGGGATTATTGGTCTGAAACTCAGCACTCATTAACTACCAGGCTGCTGTACCCAAATAAATTTACTGACATCAAACCCTTc                                                                                    >  2:407586/1‑139 (MQ=255)                                               gctgctCGGTGAGTCCGCTAACCCGATGCGCCTGGCGAGTCTGGCGTTAATCGTATTGGGGATTATTGGTCTGAAACTCAGCACTCATTAACTACCAGGCTGCTGTACCCAAATAAATTTACTGACATCAAACCCTTcc                                                                                   <  1:1464278/139‑1 (MQ=255)                                                        tGAGTCCGCTAACCCGATGCGCCTGGCGAGTCTGGCGTTAATCGTATTGGGGATTATTGGTCTGAAACTCAGCACTCATTAACTACCAGGCTGCTGt                                                                                                                    <  2:1544553/97‑1 (MQ=255)                                                        tGAGTCCGCTAACCCGATGCGCCTGGCGAGTCTGGCGTTAATCGTATTGGGGATTATTGGTCTGAAACTCAGCACTCATTAACTACCAGGCTGCTGt                                                                                                                    >  1:1544553/1‑97 (MQ=255)                                                            tCCGCTAACCCGATGCGCCTGGCGAGTCTGGCGTTAATCGTATTGGGGATTATTGGTCTGAAACTCAGCACTCATTAACTACCAGGCTGCTGTACCCAAATAAATTTACTGACATCAAACCCTTCCCGGGTCGCGACTg                                                                      <  2:946346/139‑1 (MQ=255)                                                             ccGCTAACCCGATGCGCCTGGCGAGTCTGGCGTTAATCGTATTGGGGATTATTGGTCTGAAACTCAGCACTCATTAACTACCAGGCTGCTGTACCCAAATAAATTTACTGACATCAAACCCTTCCCGGGTCGCGACTGc                                                                     >  2:930215/1‑139 (MQ=255)                                                                        aTGCGCCTGGCGAGTCTGGCGTTAATCGTATTGGGGATTATTGGTCTGAAACTCAGCACTCATTAACTACCAGGCTGCTGTACCCAAATAAATTTACTGACATCAAACCCTTCCCGGGTCGCGACTg                                                                      <  1:826767/127‑1 (MQ=255)                                                                        aTGCGCCTGGCGAGTCTGGCGTTAATCGTATTGGGGATTATTGGTCTGAAACTCAGCACTCATTAACTACCAGGCTGCTGTACCCAAATAAATTTACTGACATCAAACCCTTCCCGGGTCGCGACTg                                                                      >  2:826767/1‑127 (MQ=255)                                                                         tGCGCCTGGCGAGTCTGGCGTTAATCGTATTGGGGATTATTGGTCTGAAACTCAGCACTCATTAACTACCAGGCTGCTGTACCCAAATAAATTTACTGACATc                                                                                             >  2:1639854/1‑103 (MQ=255)                                                                         tGCGCCTGGCGAGTCTGGCGTTAATCGTATTGGGGATTATTGGTCTGAAACTCAGCACTCATTAACTACCAGGCTGCTGTACCCAAATAAATTTACTGACATc                                                                                             <  1:1639854/103‑1 (MQ=255)                                                                             ccTGGCGAGTCTGGCGTTAATCGTATTGGGGATTATTGGTCTGAAACTCAGCACTCATTAACTACCAGGCTGCTGTACCCAAATAAATTTACTGACATCa                                                                                            >  1:691431/1‑100 (MQ=255)                                                                             ccTGGCGAGTCTGGCGTTAATCGTATTGGGGATTATTGGTCTGAAACTCAGCACTCATTAACTACCAGGCTGCTGTACCCAAATAAATTTACTGACATCa                                                                                            <  2:691431/100‑1 (MQ=255)                                                                             ccTGGCGAGTCTGGCGTTAATCGTATTGGGGATTATTGGTCTGAAACTCAGCACTCATTAACTACCAGGCTGCTGTACCCAAATAAATTTACTGACATCAAACCCTTCCCGGGTCGCGACTGCCAGCATCTCCTGTTTc                                                     >  2:663884/1‑139 (MQ=255)                                                                               tGGCGAGTCTGGCGTTAATCGTATTGGGGATTATTGGTCTGAAACTCAGCACTCATTAACTACCAGGCTGCTGTACCCAAATAAATTTACTGACATCAAACCCTTCCCGGGTCGCGACTGCCAGCATCTCCTGTTTCAc                                                   <  1:407586/139‑1 (MQ=255)                                                                                      tCTGGCGTTAATCGTATTGGGGATTATTGGTCTGAAACTCAGCACTCATTAACTACCAgg                                                                                                                           <  2:1678125/60‑1 (MQ=255)                                                                                      tCTGGCGTTAATCGTATTGGGGATTATTGGTCTGAAACTCAGCACTCATTAACTACCAgg                                                                                                                           >  1:1678125/1‑60 (MQ=255)                                                                                      tCTGGCGTTAATCGTATTGGGGATTATTGGTCTGAAACTCAGCACTCATTAACTACCAGGCTGCTGTACCCAAATAAATTTACTGACATCAAACCCTTCCCGGGTCGCGACTGCCAGCATCTCCTGTTTCACTTCGt                                              >  1:1482918/1‑137 (MQ=255)                                                                                      tCTGGCGTTAATCGTATTGGGGATTATTGGTCTGAAACTCAGCACTCATTAACTACCAGGCTGCTGTACCCAAATAAATTTACTGACATCAAACCCTTCCCGGGTCGCGACTGCCAGCATCTCCTGTTTCACTTCGt                                              <  2:1482918/137‑1 (MQ=255)                                                                                      tCTGGCGTTAATCGTATTGGGGATTATTGGTCTGAAACTCAGCACTCATTAACTACCAGGCTGCTGTACCCAAATAAATTTACTGACATCAAACCCTTCCCGGGTCGCGACTGCCAGCATCTCCTGTTTCACTTCGTCa                                            >  2:777274/1‑139 (MQ=255)                                                                                      tCTGGCGTTAATCGTATTGGGGATTATTGGTCTGAAACTCAGCACTCATTAACTACCAGGCTGCTGTACCCAAATAAATTTACTGACATCAAACCCTTCCCGGGTCGCGACTGCCAGCATCTCCTGTTTCACTTCGTCa                                            >  1:826929/1‑139 (MQ=255)                                                                                      tCTGGCGTTAATCGTATTGGGGATTATTGGTCTGAAACTCAGCACTCATTAACTACCAGGCTGCTGTACCCAAATAAATTTACTGACATCAAACCCTTCCCGGGTCGCGACTGCCAGCATCTCCTGTTTCACTTCGTCa                                            <  1:663884/139‑1 (MQ=255)                                                                                        tGGCGTTAATCGTATTGGGGATTATTGGTCTGAAACTCAGCACTCATTAACTACCAg                                                                                                                            <  1:151056/57‑1 (MQ=255)                                                                                        tGGCGTTAATCGTATTGGGGATTATTGGTCTGAAACTCAGCACTCATTAACTACCAg                                                                                                                            >  2:151056/1‑57 (MQ=255)                                                                                        tGGCGTTAATCGTATTGGGGATTATTGGTCTGAAACTCAGCACTCATTAACTACCAGGCTGCTGTACCCAAATAAATTTACTGACATCAAACCCTACCCGGGTCGCGACTGCCAGCATCTCCTGTTTCACTTCGTCAGa                                          <  2:826929/139‑1 (MQ=255)                                                                                           cGTTAATCGTATTGGGGATTATTGGTCTGAAACTCAGCACTCATTAACTACCAGGCTGCTGTACCCAAATAAATTTACTGACATCAAACCCTTCCCGGGTCGCGACTGCCAGCAtctc                                                            <  2:1198549/118‑1 (MQ=255)                                                                                           cGTTAATCGTATTGGGGATTATTGGTCTGAAACTCAGCACTCATTAACTACCAGGCTGCTGTACCCAAATAAATTTACTGACATCAAACCCTTCCCGGGTCGCGACTGCCAGCAtctc                                                            >  1:1198549/1‑118 (MQ=255)                                                                                                aTCGTATTGGGGATTATTGGTCTGAAACTCAGCACTCATTAACTACCAGGCTGCTGTACCCAAATAAATTTACTGACATCAAACCCTTCCCGGGTCGCGACTGCCAGCATCTCCTGTTTCACTTCGTCAGAAATGGTTg                                  <  1:930215/139‑1 (MQ=255)                                                                                                  cGTATTGGGGATTATTGGTCTGAAACTCAGCACTCATTAACTACCAGGCTGCTGTACCCAAATAAATTTACTGACATCAAACCCTTCCCGGGTCGCGACTGCCAGCATCTCCTGTTTCACTTCGTCAGAAATGGTTGGc                                <  1:777274/139‑1 (MQ=255)                                                                                                              tattGGTCTGAAACTCAGCACTCATTAACTACCAGGCTGCTGTACCCAAATAAATTTACTGACATCAAACCCTTCCCGGGTCGCGACTGCCa                                                                   <  2:775210/92‑1 (MQ=255)                                                                                                              tattGGTCTGAAACTCAGCACTCATTAACTACCAGGCTGCTGTACCCAAATAAATTTACTGACATCAAACCCTTCCCGGGTCGCGACTGCCa                                                                   >  1:775210/1‑92 (MQ=255)                                                                                                              tattGGTCTGAAACTCAGCACTCATTAACTACCAGGCTGCTGTACCCAAATAAATTTACTGACATCAAACCCTTCCCGGGTCGCGACTGCCAGCATCTCCTGTTTCACTTCGTCAGAAATGGTTGGCGTGCGGGAGAGt                    >  1:845473/1‑139 (MQ=255)                                                                                                              tattGGTCTGAAACTCAGCACTCATTAACTACCAGGCTGCTGTACCCAAATAAATTTACTGACATCAAACCCTTCCCGGGTCGCGACTGCCAGCATCTCCTGTTTCACTTCGTCAGAAATGGTTGGCGTGCGGGAGAGt                    >  2:1300748/1‑139 (MQ=255)                                                                                                                   gTCTGAAACTCAGCACTCATTAACTACCAGGCTGCTGTACCCAAATAAATTTACTGACATCaa                                                                                           <  1:88030/63‑1 (MQ=255)                                                                                                                   gTCTGAAACTCAGCACTCATTAACTACCAGGCTGCTGTACCCAAATAAATTTACTGACATCaa                                                                                           >  2:88030/1‑63 (MQ=255)                                                                                                                   gTCTGAAACTCAGCACTCATTAACTACCAGGCTGCTGTACCCAAATAAATTTACTGACATCAAACCCTTCCCGGGTCGCGACTGCCAGCATCTCCTGTTTCACTTCGTCAGAAATGGt                                    >  2:1164556/1‑118 (MQ=255)                                                                                                                   gTCTGAAACTCAGCACTCATTAACTACCAGGCTGCTGTACCCAAATAAATTTACTGACATCAAACCCTTCCCGGGTCGCGACTGCCAGCATCTCCTGTTTCACTTCGTCAGAAATGGt                                    <  1:1164556/118‑1 (MQ=255)                                                                                                                   gTCTGAAACTCAGCACTCATTAACTACCAGGCTGCTGTACCCAAATAAATTTACTGACATCAAACCCTTCCCGGGTCGCGACTGCCAGCATCTCCTGTTTCACTTCGTCAGAAATGGTTg                                  >  2:1248148/1‑120 (MQ=255)                                                                                                                   gTCTGAAACTCAGCACTCATTAACTACCAGGCTGCTGTACCCAAATAAATTTACTGACATCAAACCCTTCCCGGGTCGCGACTGCCAGCATCTCCTGTTTCACTTCGTCAGAAATGGTTg                                  <  1:1248148/120‑1 (MQ=255)                                                                                                                        aaaCTCAGCACTCATTAACTACCAGGCTGCTGTACCCAAATAAATTTACTGACATCAAACCCTTCCCGGGTCGCGACTGCCAGCATCTCCTGTTTCACt                                                  >  2:942240/1‑99 (MQ=255)                                                                                                                        aaaCTCAGCACTCATTAACTACCAGGCTGCTGTACCCAAATAAATTTACTGACATCAAACCCTTCCCGGGTCGCGACTGCCAGCATCTCCTGTTTCACt                                                  <  1:942240/99‑1 (MQ=255)                                                                                                                             cAGCACTCATTAACTACCAGGCTGCTGTACCCAAATAAATTTACTGACATCAAACCCTTCCCGGGTCGCGACTGCCAGCATCTCCTGTTTCACTTCGTCAGAAATGGTTGGCGTGCGGGAGAGTATCCACAGGTAGTcg     <  1:1343337/139‑1 (MQ=255)                                                                                                                              aGCACTCATTAACTACCAGGCTGCTGTACCCAAATAAATTTACTGACATCAAACCCTTCCCGGGTCGCGACTGCCAGCAtc                                                              <  1:873393/81‑1 (MQ=255)                                                                                                                              aGCACTCATTAACTACCAGGCTGCTGTACCCAAATAAATTTACTGACATCAAACCCTTCCCGGGTCGCGACTGCCAGCAtc                                                              >  2:873393/1‑81 (MQ=255)                                                                                                                                cactcattAACTACCAGGCTGCTGTACCCAAATAAATTTACTGACATCAAACCCTTCCCGGGTCGCGACTGCCAGCATCTCCTGTTTCACTTCGTCAGAAATGGTTGGCGTg                             <  1:439864/112‑1 (MQ=255)                                                                                                                                cactcattAACTACCAGGCTGCTGTACCCAAATAAATTTACTGACATCAAACCCTTCCCGGGTCGCGACTGCCAGCATCTCCTGTTTCACTTCGTCAGAAATGGTTGGCGTg                             >  2:439864/1‑112 (MQ=255)                                                                                                                                cactcattAACTACCAGGCTGCTGTACCCAAATAAATTTACTGACATCAAACCCTTCCCGGGTCGCGACTGCCAGCATCTCCTGTTTCACTTCGTCAGAAATGGTTGGCGTGCGGGAGAGTATCCACAGGTAGTCGCgg  >  2:1603535/1‑139 (MQ=255)                                                                                                                                  ctcattAACTACCAGGCTGCTGTACCCAAATAAATTTACt                                                                                                   <  2:1223529/40‑1 (MQ=255)                                                                                                                                  ctcattAACTACCAGGCTGCTGTACCCAAATAAATTTACt                                                                                                   >  1:1223529/1‑40 (MQ=255)                                                                                                                                   tcattAACTACCAGGCTGCTGTACCCAAATAAATTTACTGACATCAAACCCTTCCCGGGTCGCGACTGCCAGCATCTCCTGTTTCACTTCGTCAGAAATGGTTGGCGTGCGGGAGAGTATCcaca             >  1:937174/1‑125 (MQ=255)                                                                                                                                   tcattAACTACCAGGCTGCTGTACCCAAATAAATTTACTGACATCAAACCCTTCCCGGGTCGCGACTGCCAGCATCTCCTGTTTCACTTCGTCAGAAATGGTTGGCGTGCGGGAGAGTATCcaca             <  2:937174/125‑1 (MQ=255)                                                                                                                                      |                                                                                                                                     TGTGGACGGGTATTGGCGCAGTCGGCGCGGCCATAACCGGCATTGTGCTGCTCGGTGAGTCCGCTAACCCGATGCGCCTGGCGAGTCTGGCGTTAATCGTATTGGGGATTATTGGTCTGAAACTCAGCACTCACTAACTACCAGGCTGCTGTACCCAAATAAATTTACTGACATCAAACCCTTCCCGGGTCGCGACTGCCAGCATCTCCTGTTTCACTTCGTCAGAAATGGTTGGCGTGCGGGAGAGTATCCACAGGTAGTCGCGG  >  NC_000913/4377056‑4377321

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