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%	151.1 / NA	47	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/24);  total (23/24)

GACGGGTATTGGCGCAGTCGGCGCGGCCATAACCGGCATTGTGCTGCTCGGTGAGTCCGCTAACCCGATGCGCCTGGCGAGTCTGGCGTTAATCGTATTGGGGATTATTGGTCTGAAACTCAGCACTCACTAACTACCAGGCTGCTGTACCCAAATAAATTTACTGACATCAAACCCTTCCCGGGTCGCGACTGCCAGCATCTCCTGTTTCACTTCGTCAGAAATGGTTGGCGTGCGGGAGAGTATCCACAGGTAGTCGC  >  NC_000913/4377060‑4377319                                                                                                                                  |                                                                                                                                   gACGGGTATTGGCGCAGTCGGCGCGGCCATAACCGGCATTGTGCTGCTCGGTGAGTCCGCTAACCCGATGCGCCTGGCGAGTCTGGCGTTAATCGTATTGGGGATTATTGGTCTGAAACTCAGCACTCATTAACTACCa                                                                                                                           >  2:1225587/1‑139 (MQ=255)   cGGGTATTGGCGCAGTCGGCGCGGCCATAACCGGCATTGTGCTGCTCGGTGAGTCCGCTAACCCGATGCGCCTGGCGAGTCTGGCGTTAATCGTATTGGGGATTATTGGTCTGAAACTCAGCACTCATTAACTACCAgg                                                                                                                         <  2:529879/139‑1 (MQ=255)   cGGGTATTGGCGCAGTCGGCGCGGCCATAACCGGCATTGTGCTGCTCGGTGAGTCCGCTAACCCGATGCGCCTGGCGAGTCTGGCGTTAATCGTATTGGGGATTATTGGTCTGAAACTCAGCACTCATTAACTACCAgg                                                                                                                         <  2:1241315/139‑1 (MQ=255)                aGTCGGCGCGGCCATAACCGGCATTGTGCTGCTCGGTGAGTCCGCTAACCCGATGCGCCTGGCGAGTCTGGCGTTAATCGTATTGGGGATTATTGGTCTGAAACTCAGCACTCATTAACTACCAGGCTGCTGTACCCaa                                                                                                            >  2:751874/1‑139 (MQ=255)                           ccATAACCGGCATTGTGCTGCTCGGTGAGTCCGCTAACCCGATGCGCCTGGCGAGTCTGGCGTTAATCGTATTGGGGATTATTGGTCTGAAACTCAGCACTCATTAACTACCAGGCTGCTGTACCCAAATAAATTTACt                                                                                                 <  1:1369623/139‑1 (MQ=255)                            cATAACCGGCATTGTGCTGCTCGGTGAGTCCGCTAACCCGATGCGCCTGGCGAGTCTGGCGTTAATCGTATTGGGGATTATTGGTCTGAAACTCAGCACTCATTAACTACCAgg                                                                                                                         <  2:703629/114‑1 (MQ=255)                            cATAACCGGCATTGTGCTGCTCGGTGAGTCCGCTAACCCGATGCGCCTGGCGAGTCTGGCGTTAATCGTATTGGGGATTATTGGTCTGAAACTCAGCACTCATTAACTACCAgg                                                                                                                         >  1:703629/1‑114 (MQ=255)                            cATAACCGGCATTGTGCTGCTCGGTGAGTCCGCTAACCCGATGCGCCTGGCGAGTCTGGCGTTAATCGTATTGGGGATTATTGGTCTGAAACTCAGCACTCATTAACTACCAGGCTGCTGTACCCAAATAAATTTACt                                                                                                 <  1:701252/138‑1 (MQ=255)                            cATAACCGACATTGTGCTGCTCGGTGAGTCCGCTAACCCGATGCGCCTGGCGAGTCTGGCGTTAATCGTATTGGGGATTATTGGTCTGAAACTCAGCACTCATTAACTACCAGGCTGCTGTACCCAAATAAATTTACt                                                                                                 >  2:701252/1‑138 (MQ=255)                                  cGGCATTGTGCTGCTCGGTGAGTCCGCTAACCCGATGCGCCTGGCGAGTCTGGCGTTAATCGTATTGGGGATTATTGGTCTGAAACTCAGCACTCATTAACTACCAGGCTGCTGTACCCAAATAAATTTACTGACATCa                                                                                          <  2:64818/139‑1 (MQ=255)                                    gCATTGTGCTGCTCGGTGAGTCCGCTAACCCGATGCGCCTGGCGAGTCTGGCGTTAATCGTATTGGGGATTATTGGTCTGAAACTCAGCACTCATTAACTACCAGGCTGCTGTACCCAAATAAATTTACTGACATCaaa                                                                                        <  1:518285/139‑1 (MQ=255)                                     cATTGTGCTGCTCGGTGAGTCCGCTAACCCGATGCGCCTGGCGAGTCTGGCGTTAATCGTATTGGGGATTATTGGTCTGAAACTCAGCACTCATTAACTACCAGGCTGCTGTACCCAAATAAATTTACTGACATCAAAc                                                                                       <  2:577036/139‑1 (MQ=255)                                     cATTGTGCTGCTCGGTGAGTCCGCTAACCCGATGCGCCTGGCGAGTCTGGCGTTAATCGTATTGGGGATTATTGGTCTGAAACTCAGCACTCATTAACTACCAGGCTGCTGTACCCAAATAAATTTACTGACATCAAAc                                                                                       >  2:1344216/1‑139 (MQ=255)                                        tgtgCTGCTCGGTGAGTCCGCTAACCCGATGCGCCTGGCGAGTCTGGCGTTAATCGTATTGGGGATTATTGGTCTGAAACTCAGCACTCATTAACTACCAGGCTGCTGTACCCAAATAAATTTACTGACATCAAACCCt                                                                                    >  1:1261280/1‑139 (MQ=255)                                          tgctgcTCGGTGAGTCCGCTAACCCGATGCGCCTGGCGAGTCTGGCGTTAATCGTATTGGGGATTATTGGTCTGAAACTCAGCACTCATTAACTACCAGGCTGCTGTACCCaaa                                                                                                           >  2:1418493/1‑114 (MQ=255)                                          tgctgcTCGGTGAGTCCGCTAACCCGATGCGCCTGGCGAGTCTGGCGTTAATCGTATTGGGGATTATTGGTCTGAAACTCAGCACTCATTAACTACCAGGCTGCTGTACCCaaa                                                                                                           <  1:1418493/114‑1 (MQ=255)                                                tCGGTGAGTCCGCTAACCCGATGCGCCTGGCGAGTCTGGCGTTAATCGTATTGGGGATTATTGGTCTGAAACTCAGCACTCATTAACTACCAGGCTGCTGTACCCAAATAAATTTACTGACATCAAACCCTTCCCGGGt                                                                            <  1:714711/139‑1 (MQ=255)                                                    tGAGTCCGCTAACCCGATGCGCCTGGCGAGTCTGGCGTTAATCGTATTGGGGATTATTGGTCTGAAACTCAGCACTCATTAACTACCAGgctg                                                                                                                      <  1:958915/93‑1 (MQ=255)                                                    tGAGTCCGCTAACACGATGCGCCTGGCGAGTCTGGCGTTAATCGTATTGGGGATTATTGGTCTGAAACTCAGCACTCATTAACTACCAGgctg                                                                                                                      >  2:958915/1‑93 (MQ=255)                                                          cGCTAACCCGATGCGCCTGGCGAGTCTGGCGTTAATCGTATTGGGGATTATTGGTCTGAAACTCAGCACTCATTAACTACCAGGCTGCTGTACCCAAATAAATTTACTGACATCAAACCCTTCCCGGGTCGCGACTGcc                                                                  <  2:436974/139‑1 (MQ=255)                                                                cccGATGCGCCTGGCGAGTCTGGCGTTAATCGTATTGGGGATTATTGGTCTGAAACTCAGCACTCATTAACTACCAGGCTGCTGTACCCAAATAAATTTACTGACATCAAACCCTTCCCGGGTCGCGACTGCCAGCAtc                                                            >  1:120759/1‑139 (MQ=255)                                                                         ccTGGCGAGTCTGGCGTTAATCGTATTGGGGATTATTGGTCTGAAACTCAGCACTCATTAACTACCAGgctg                                                                                                                      >  1:1264474/1‑72 (MQ=255)                                                                         ccTGGCGAGTCTGGCGTTAATCGTATTGGGGATTATTGGTCTGAAACTCAGCACTCATTAACTACCAGgctg                                                                                                                      <  2:1264474/72‑1 (MQ=255)                                                                         ccTGGCGAGTCTGGCGTTAATCGTATTGGGGATTATTGGTCTGAAACTCAGCACTCATTAACTACCAGGCTGCTGTACCCAAATAAATTTACTGACATCAAACCCTTCCCGGGTCGCGACTGCCAGCAtctc                                                          <  2:1094526/132‑1 (MQ=255)                                                                         ccTGGCGAGTCTGGCGTTAATCGTATTGGGGATTATTGGTCTGAAACTCAGCACTCATTAACTACCAGGCTGCTGTACCCAAATAAATTTACTGACATCAAACCCTTCCCGGGTCGCGACTGCCAGCAtctc                                                          >  1:1094526/1‑132 (MQ=255)                                                                                     ggCGTTAATCGTATTGGGGATTA‑TGGTCTGAAACTCAGCACTCATTAACTACCAGGCTGCTGTACCCAAATAAATTTACTGACATCAAACCCTTCCCGGGTCGCGACTGCCAGCATCTCCTGTTTCACTTCGTCAGaaa                                      <  1:401876/139‑1 (MQ=255)                                                                                          tAATCGTATTGGGGATTATTGGTCTGAAACTCAGCACTCATTAACTACCAGGCTGCTGTACCCAAATAAATTTACTGACATCAAACCCTTCCCGGGTCGCGACTGCCAGCATCTCCTGTTTCACTTCGTCAGAAATGGt                                  <  2:1261280/139‑1 (MQ=255)                                                                                            aTCGTATTGGGGATTATTGGTCTGAAACTCAGCACTCATTAACTACCAGGCTGCTGTACCCAAGTAAATTTACTGACATCAAACCCTTCCCGGGTCGCGACTGCCAGCATCTCCTGTTTCACTTCGTCAGAAATGGTTg                                >  2:932066/1‑139 (MQ=255)                                                                                                tATTGGGGATTATTGGTCTGAAACTCAGCACTCATTAACTACCAGGCTGCTGTACCCAAATAAATTTACt                                                                                                 >  1:402696/1‑70 (MQ=255)                                                                                                tATTGGGGATTATTGGTCTGAAACTCAGCACTCATTAACTACCAGGCTGCTGTACCCAAATAAATTTACt                                                                                                 <  2:402696/70‑1 (MQ=255)                                                                                                 aTTGGGGATTATTGGTCTGAAACTCAGCACTCATTAACTACCAGGCTGCTGTACCCAAATAAATTTACTGACATCAAACCCTTCCCGGGTCGCGACTGCCAGCATCTCCTGTTTCACTTCGTCAGAAATGGTTGGCGTg                           >  2:553623/1‑139 (MQ=255)                                                                                                       gATTATTGGTCTGAAACTCAGCACTCATTAACTACCAGGCTGCTGTACCCAAATAAATTTACTGACATCAAACCCTTCCCGGGTCGCGACTGCCAGCATCTCCTGTTTCACTTCGTCAGAAATGGTTGGCGTGCGGgag                     <  2:120759/139‑1 (MQ=255)                                                                                                          tattGGTCTGAAACTCAGCACTCATTAACTACCAGGCTGCTGTACCCAAATAAATTTACt                                                                                                 <  2:222202/60‑1 (MQ=255)                                                                                                          tattGGTCTGAAACTCAGCACTCATTAACTACCAGGCTGCTGTACCCAAATAAATTTACt                                                                                                 >  1:222202/1‑60 (MQ=255)                                                                                                          tattGGTCTGAAACTCAGCACTCATTAACTACCAGGCTGCTGTACCCAAATAAATTTACTGACATCAAACCCTTCCCGGGTCGCGACTGCCAGCATCTCCTGTTTCACTTCGTCAGAAATGGTTGGCGTGCGGGAGAGt                  >  2:958117/1‑139 (MQ=255)                                                                                                              ggTCTGAAACTCAGCACTCATTAACTACCAGGCTGCTGTACCCAAATAAATTTACTGACATCAAACCCTTCCCGGGTCGCGACTGCCAGCATCTCCTGTTTCACTTCGTCAGAAATGGTTGGCGTGCGGGAGAGTATcc              <  1:751874/139‑1 (MQ=255)                                                                                                               gTCTGAAACTCAGCACTCATTAACTACCAGGCTGCTGTACCCAAATAAATTTACTGACATCAAACCCTTCCCGGGTCGCGACTGCCAGCATCTCCTGTTTCACTTCGTCAGa                                        >  1:1336427/1‑112 (MQ=255)                                                                                                               gTCTGAAACTCAGCACTCATTAACTACCAGGCTGCTGTACCCAAATAAATTTACTGACATCAAACCCTTCCCGGGTCGCGACTGCCAGCATCTCCTGTTTCACTTCGTCAGa                                        <  2:1336427/112‑1 (MQ=255)                                                                                                               gTCTGAAACTCAGCACTCATTAACTACCAGGCTGCTGTACCCAAATAAATTTACTGACATCAAACCCTTCCCGGGTCGCGACTGCCAGCATCTCCTGTTTCACTTCGTCAGAAATGGTTg                                >  1:635749/1‑120 (MQ=255)                                                                                                               gTCTGAAACTCAGCACTCATTAACTACCAGGCTGCTGTACCCAAATAAATTTACTGACATCAAACCCTTCCCGGGTCGCGACTGCCAGCATCTCCTGTTTCACTTCGTCAGAAATGGTTg                                <  2:635749/120‑1 (MQ=255)                                                                                                                 cTGAAACTCAGCACTCATTAACTACCAGGCTGCTGTACCCAAATAAATTTACTGACATCAAACCCTTcc                                                                                 >  2:261950/1‑69 (MQ=255)                                                                                                                 cTGAAACTCAGCACTCATTAACTACCAGGCTGCTGTACCCAAATAAATTTACTGACATCAAACCCTTcc                                                                                 <  1:261950/69‑1 (MQ=255)                                                                                                                    aaaCTCAGCACTCATTAACTACCAGGCTGCTGTACCCAAATAAATTTACt                                                                                                 <  1:1265322/50‑1 (MQ=255)                                                                                                                    aaaCTCAGCACTCATTAACTACCAGGCTGCTGTACCCAAATAAATTTACt                                                                                                 >  2:1265322/1‑50 (MQ=255)                                                                                                                          aGCACTCATTAACTACCAGGCTGCTGTACCCAAATAAATTTACTGACATCAAACCCTTCCCGGGTCGCGACTGCCAGCATCTCCTGTTTCACTTCGTCAGAAATGGTTGGCGTGCGGGAGAGTATCCACAGGTAGTcgc  >  2:318301/1‑139 (MQ=255)                                                                                                                          aGCACTCATTAACTACCAGGCTGCTGTACCCAAATAAATTTACTGACATCAAACCCTTCCCGGGTCGCGACTGCCAGCATCTCCTGTTTCACTTCGTCAGAAATGGTTGGCGTGCGGGAGAGTATCCACAGGTAGTcgc  >  2:476794/1‑139 (MQ=255)                                                                                                                          aGCACTCATTAACTACCAGGCTGCTGTACCCAAATAAATTTACTGACATCAAACCCTTCCCGGGTCGCGACTGCCAGCATCTCCTGTTTCACTTCGTCAGAAATGGTTGGCGTGCGGGAGAGTATCCACAGGTAGTcgc  >  2:888348/1‑139 (MQ=255)                                                                                                                                  |                                                                                                                                   GACGGGTATTGGCGCAGTCGGCGCGGCCATAACCGGCATTGTGCTGCTCGGTGAGTCCGCTAACCCGATGCGCCTGGCGAGTCTGGCGTTAATCGTATTGGGGATTATTGGTCTGAAACTCAGCACTCACTAACTACCAGGCTGCTGTACCCAAATAAATTTACTGACATCAAACCCTTCCCGGGTCGCGACTGCCAGCATCTCCTGTTTCACTTCGTCAGAAATGGTTGGCGTGCGGGAGAGTATCCACAGGTAGTCGC  >  NC_000913/4377060‑4377319

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