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%	134.7 / NA	42	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 (18/24);  total (18/24)

TGTGGACGGGTATTGGCGCAGTCGGCGCGGCCATAACCGGCATTGTGCTGCTCGGTGAGTCCGCTAACCCGATGCGCCTGGCGAGTCTGGCGTTAATCGTATTGGGGATTATTGGTCTGAAACTCAGCACTCACTAACTACCAGGCTGCTGTACCCAAATAAATTTACTGACATCAAACCCTTCCCGGGTCGCGACTGCCAGCATCTCCTGTTT‑‑CACTTCGTCAGAAATGGTTGGCGTGCGGGAGAGTATCCACAGGTAGTCG  >  NC_000913/4377056‑4377318                                                                                                                                      |                                                                                                                                    tgtgGACGGGTATTGGCGCAGTCGGCGCGGCCATAACCGGCATTGTGCTGCTCGGTGAGTCCGCTAACCCGATGCGCCTGGCGAGTCTGGCGTTAATCGTATTGGGGATTATTGGTCTGAAACTCAGCACTCattaact                                                                                                                                >  2:1400813/1‑139 (MQ=255)  gtgGACGGGTATTGGCGCAGTCGGCGCGGCCATAACCGGCATTGTGCTGCTCGGTGAGTCCGCTAACCCGATGCGCCTGGCGAGTCTGGCGTTAATCGTATTGGGGATTATTGGTCTGAAACTCAGCACTCattaacta                                                                                                                               <  1:229174/139‑1 (MQ=255)   tgGACGGGTATTGGCGCAGTCGGCGCGGCCATAACCGGCATTGTGCTGCTCGGTGAGTCCGCTAACCCGATGCGCCTGGCGAGTCTGGCGTTAATCGTATTGGGGATTATTGGTCTGAAACTCAGCACTCATTAACTAc                                                                                                                              >  1:597456/1‑139 (MQ=255)     gACGGGTATTGGCGCAGTCGGCGCGGCCATAACCGGCATTGTGCTGCTCGGTGAGTCCGCTAACCCGATGCGCCTGGCGAGTCTGGCGTTAATCGTATTGGGGATTATTGGTCTGAAACTCAGCACTCATTAACTACCa                                                                                                                            <  1:1400813/139‑1 (MQ=255)       cGGGTATTGGCGCAGTCGGCGCGGCCATAACCGGCATTGTGCTGCTCGGTGAGTCCGCTAACCCGATGCGCCTGGCGAGTCTGGCGTTAATCGTATTGGGGATTATTGGTCTGAAACTCAGCACTCATTAACTACCAgg                                                                                                                          <  2:1003642/139‑1 (MQ=255)       cGGGTATTGGCGCAGTCGGCGCGGCCATAACCGGCATTGTGCTGCTCGGTGAGTCCGCTAACCCGATGCGCCTGGCGAGTCTGGCGTTAATCGTATTGGGGATTATTGGTCTGAAACTCAGCACTCATTAACTACCAgg                                                                                                                          <  2:1457187/139‑1 (MQ=255)            aTTGGCGCAGTCGGCGCGGCCATAACCGGCATTGTGCTGCTCGGTGAGTCCGCTAACCCGATGCGCCTGGCGAGTCTGGCGTTAATCGTATTGGGGATTATTGGTCTGAAACTCAGCACTCATTAACTACCAGgctgct                                                                                                                     <  1:627039/139‑1 (MQ=255)                   cAGTCGGCGCGGCCATAACCGGCATTGTGCTGCTCGGTGAGTCCGCTAACCCGATGCGCCTGGCGAGTCTGGCGTTAATCGTATTGGGGATTATTGGTCTGAAACTCAGCACTCATTAACTACCAGgctgc                                                                                                                      <  1:950574/131‑1 (MQ=255)                   cAGTCGGCGCGGCCATAACCGGCATTGTGCTGCTCGGTGAGTCCGCTAACCCGATGCGCCTGGCGAGTCTGGCGTTAATCGTATTGGGGATTATTGGTCTGAAACTCAGCACTCATTAACTACCAGgctgc                                                                                                                      >  2:950574/1‑131 (MQ=255)                     gTCGGCGCGGCCATAACCGGCATTGTGCTGCTCGGTGAGTCCGCTAACCCGATGCGCCTGGCGAGTCTGGCGTTAATCGTATTGGGGATTATTGGTCTGAAACTCAGCACTCATTAACTACCAGgctgct                                                                                                                     <  1:1033382/130‑1 (MQ=255)                     gTCGGCGCGGCCATAACCGGCATTGTGCTGCTCGGTGAGTCCGCTAACCCGATGCGCCTGGCGAGTCTGGCGTTAATCGTATTGGGGATTATTGGTCTGAAACTCAGCACTCATTAACTACCAGgctgct                                                                                                                     >  2:1033382/1‑130 (MQ=255)                                        gCATTGTGCTGCTCGGTGAGTCCGCTAACCCGATGCGCCTGGCGAGTCTGGCGTTAATCGTATTGGGGATTATTGGTCTGAAACTCAGCACTCATTAACTACCAGGCTGCTGTACCCAAATAAATTTACTGACATCaaa                                                                                         <  2:597456/139‑1 (MQ=255)                                              tgctgcTCGGTGAGTCCGCTAACCCGATGCGCCTGGCGAGTCTGGCGTTAATCGTATTGGGGATTATTGGTCTGAAACTCAGCACTCATTAACTACCAGGCTGCTGTACCCAAATAAATTTACTGACATCAAACCCTTc                                                                                   >  1:77939/1‑139 (MQ=255)                                                  gctCGGTGAGTCCGCTAACCCGATGCGCCTGGCGAGTCTGGCGTTAATCGTATTGGGGATTATTGGTCTG‑AACTCAGCACTCATTAACTACCAGGCTGCTGTACCCAAATAAATTTACTGACATCAAACCCTTCCCggg                                                                              <  1:1082367/139‑1 (MQ=255)                                                         gAGTCCGCTAACCCGATGCGCCTGGCGAGTCTGGCGTTAATCGTATTGGGGATTATTGGTCTGAAACTCAGCACTCATTAACTACCAGGCTGCTGTACCCAAATAAATTTACTGACATCAAACCCTTCCCGGGTCGCGa                                                                        <  2:135404/139‑1 (MQ=255)                                                             ccGCTAACCCGATGCGCCTGGCGAGTCTGGCGTTAATCGTATTGGGGATTATTGGTCTGAAACTCAGCACTCATTAACTACCAGGCTGCTGTACCCAAATAAATTTACTGACATCAAACCCTTCCCGGGTCGCGACTGc                                                                    >  1:577715/1‑139 (MQ=255)                                                                      cGATGCGCCTGGCGAGTCTGGCGTTAATCGTATTGGGGATTATTGGTCTGAAACTCAGCACTCATTAACTACCAGGCTGCTGTACCCAAATAAATTTACTGACATCAAACCCTTCCCGGGTCGCGACTGCCAGCAtctc                                                           <  2:248890/139‑1 (MQ=255)                                                                           cgcCTGGCGAGTCTGGCGTTAATCGTATTGGGGATTATTGGTCTGAAACTCAGCACTCATTAACTACCAGGCTGCTGTACCCAAATAAATTTACTGACATCAAACCCTTCCCGGGTCGCGACTGCCAGCATCTCCTGtt                                                      >  2:283722/1‑139 (MQ=255)                                                                             ccTGGCGAGTCTGGCGTTAATCGTATTGGGGATTATTGGTCTGAAACTCAGCACTCattaacta                                                                                                                               <  1:954215/64‑1 (MQ=255)                                                                             ccTGGCGAGTCTGGCGTTAATCGTATTGGGGATTATTGGTCTGAAACTCAGCACTCattaacta                                                                                                                               >  2:954215/1‑64 (MQ=255)                                                                                        tGGCGTTAATCGTATTGGGGATTATTGGTCTGAAACTCAGCACTCATTAACTACCAGGCTGCTGTACCCAAATAAATTTACTGACATCAAACCCTTCCCGGGTCGCGACTGCCAGCATCTCCTGTTT‑‑CACTTCGTCAGa                                       <  2:77939/139‑1 (MQ=255)                                                                                        tGGCGTTAATCGTATTGGGGATTATTGGTCTGAAACTCAGCACTCATTAACTACCAGGCTGCTGTACCCAAATAAATTTACTGACATCAAACCCTTCCCGGGTCGCGACTGCCAGCATCTCCTGTTT‑‑CACTTCGTCAGa                                       >  1:124629/1‑139 (MQ=255)                                                                                               aaTCGTATTGGGGATTATTGGTCTGAAACTCAGCACTCATTAACTACCAGGCTGCTGTACCCAAATAAATTTACTGACATCAAACCCTTCCCGGGTCGCGACTGCCAGCATCTCCTGt                                                       >  1:1125415/1‑118 (MQ=255)                                                                                               aaTCGTATTGGGGATTATTGGTCTGAAACTCAGCACTCATTAACTACCAGGCTGCTGTACCCAAATAAATTTACTGACATCAAACCCTTCCCGGGTCGCGACTGCCAGCATCTCCTGt                                                       <  2:1125415/118‑1 (MQ=255)                                                                                               aaTCGTATTGGGGATTATTGGTCTGAAACTCAGCACTCATTAACTACCAGGCTGCTGTACCCAAATAAATTTACTGACATCAAACCCTTCCCGGGTCGCGACTGCCAGCATCTCCTGTTT‑‑CACTTCGTCAGAAATGGtt                                >  2:11362/1‑139 (MQ=255)                                                                                                  cGTATTGGGGATTATTGGTCTGAAACTCAGCACTCATTAACTACCAGGCTGCTGTACCCAAATAAATTTACTGACATCAAACCCTTCCCGGGTCGCGACTGCCAGCATCTCCTGTTT‑‑CACTTCGTCAGAAATGGTTGGc                             <  2:601921/139‑1 (MQ=255)                                                                                                    tATTGGGGATTATTGGTCTGAAACTCAGCACTCATTAACTACCAGGCTGCTGTACCCAAATAAATTTACTGACATCAAACCCTTCCCGGGTCGCGa                                                                        >  1:206412/1‑96 (MQ=255)                                                                                                    tATTGGGGATTATTGGTCTGAAACTCAGCACTCATTAACTACCAGGCTGCTGTACCCAAATAAATTTACTGACATCAAACCCTTCCCGGGTCGCGa                                                                        <  2:206412/96‑1 (MQ=255)                                                                                                     aTTGGGGATTATTGGTCTGAAACTCAGCACTCATTAACTACCAGGCTGCTGTACCCAAATAAATTTACTGACATCAAACCCTTCCCGGGTCGCGACTGCCAGCATCTCCTGTTT‑‑CACTTCGTCAGAAATGGTTGGCGTg                          <  2:124629/139‑1 (MQ=255)                                                                                                     aTTGGGGATTATTGGTCTGAAACTCAGCACTCATTAACTACCAGGCTGCTGTACCCAAATAAATTTACTGACATCAAACCCTTCCCGGGTCGCGACTGCCAGCATCTCCTGTTT‑‑CACTTCGTCAGAAATGGTTGGCGTg                          <  2:577715/139‑1 (MQ=255)                                                                                                         gggATTATTGGTCTGAAACTCAGCACTCATTAACTACCAGGCTGCTGTACCCAAATAAATTTACTGACATCAAACCCTTCCCGGGTCGCGACTGCCAGCATCTCCTGTTTCACACTTCGTCAGAAATGGTTGGCGTGCg                        >  1:641158/1‑139 (MQ=255)                                                                                                         gggATTATTGGTCTGAAACTCAGCACTCATTAACTACCAGGCTGCTGTACCCAAATAAATTTACTGACATCAAACCCTTCCCGGGTCGCGACTGCCAGCATCTCCTGTTT‑‑CACTTCGTCAGAAATGGTTGGCGTGCggg                      <  1:1383418/139‑1 (MQ=255)                                                                                                          ggATTATTGGTCTGAAACTCAGCACTCATTAACTACCAGGCTGCTGTACCCAAATAAATTTACTGACATCAAACCCTTCCCGGGTCGCGACTGCCAGCATCTCCTGTTT‑‑CACTTCGTCAGAAATGGTTGGCGTGCGGga                     >  2:1462706/1‑139 (MQ=255)                                                                                                           gATTATTGGTCTGAAACTCAGCACTCATTAACTACCa                                                                                                                            <  1:1186384/37‑1 (MQ=255)                                                                                                           gATTATTGGTCTGAAACTCAGCACTCATTAACTACCa                                                                                                                            >  2:1186384/1‑37 (MQ=255)                                                                                                              tattGGTCTGAAACTCAGCACTCATTAACTACCAGGCTGCTGTACCCAAATAAATTTACTGACATCAAACCCTTCCCGGGTCGCGAc                                                                       >  2:155731/1‑87 (MQ=255)                                                                                                              tattGGTCTGAAACTCAGCACTCATTAACTACCAGGCTGCTGTACCCAAATAAATTTACTGACATCAAACCCTTCCCGGGTCGCGAc                                                                       <  1:155731/87‑1 (MQ=255)                                                                                                              tattGGTCTGAAACTCAGCACTCATTAACTACCAGGCTGCTGTACCCAAATAAATTTACTGACATCAAACCCTTCCCGGGTCGCGACTGCCAGCATCTCCTGTTTCACACTTCGTCAGAAATGGTTGGCGTGCGGgaga                   <  2:641158/139‑1 (MQ=255)                                                                                                              tattGGTCTGAAACTCAGCACTCATTAACTACCAGGCTGCTGTACCCAAATAAATTTACTGACATCAAACCCTTCCCGGGTCGCGACTGCCAGCATCTCCTGTTT‑‑CACTTCGTCAGAAATGGTTGGCg                            <  2:682287/128‑1 (MQ=255)                                                                                                              tattGGTCTGAAACTCAGCACTCATTAACTACCAGGCTGCTGTACCCAAATAAATTTACTGACATCAAACCCTTCCCGGGTCGCGACTGCCAGCATCTCCTGTTT‑‑CACTTCGTCAGAAATGGTTGGCg                            >  1:682287/1‑128 (MQ=255)                                                                                                                   gTCTGAAACTCAGCACTCATTAACTACCAGGCTGCTGTACCCAAATAAATTTACTGACATCAAACCCTTCCCGGGTCGCGACTGCCAg                                                                 <  1:11726/88‑1 (MQ=255)                                                                                                                   gTCTGAAACTCAGCACTCATTAACTACCAGGCTGCTGTACCCAAATAAATTTACTAACATCAAACCCTTCCCGGGTCGCGACTGCCAg                                                                 >  2:11726/1‑88 (MQ=255)                                                                                                                        aaaCTCAGCACTCATTAACTACCAGGCTGCTGTACCCAAATAAATTTACTGACATCAAACCCTTCCCGGGTCGCGACTGCCAGCATCTCCTGTTT‑‑CACTTCGTCAGAAATGGTTGGCGTGCGGGAGAGTATCCACAGGt       >  2:1561362/1‑139 (MQ=255)                                                                                                                             cAGCACTCATTAACTACCAGGCTGCTGTACCCAAATAAATTTACTGACATCAAACCCTTCCCGGGTCGCGACTGCCAGCATCTCCTGTTT‑‑CACTTCGTCAGAAATGGTTGGCGTGCGGGAGAGTATCCACAGGTAGTcg  <  1:1294178/139‑1 (MQ=255)                                                                                                                              aGCACTCATTAACTACCAGGCTGCTGTACCCAAATAAATTTACTGACATCAAACCCTTCCCGGGTCGCGACTGCCAGCATCTCCTGTTT‑‑CACTTCGTCAGAAATGGTTGGCGTGCggg                      <  2:1485480/118‑1 (MQ=255)                                                                                                                              aGCACTCATTAACTACCAGGCTGCTGTACCCAAATAAATTTACTGACATCAAACCCTTCCCGGGTCGCGACTGCCAGCATCTCCTGTTT‑‑CACTTCGTCAGAAATGGTTGGCGTGCggg                      >  1:1485480/1‑118 (MQ=255)                                                                                                                                      |                                                                                                                                    TGTGGACGGGTATTGGCGCAGTCGGCGCGGCCATAACCGGCATTGTGCTGCTCGGTGAGTCCGCTAACCCGATGCGCCTGGCGAGTCTGGCGTTAATCGTATTGGGGATTATTGGTCTGAAACTCAGCACTCACTAACTACCAGGCTGCTGTACCCAAATAAATTTACTGACATCAAACCCTTCCCGGGTCGCGACTGCCAGCATCTCCTGTTT‑‑CACTTCGTCAGAAATGGTTGGCGTGCGGGAGAGTATCCACAGGTAGTCG  >  NC_000913/4377056‑4377318

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