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%	117.3 / NA	37	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 (17/20);  total (17/20)

TGTGGACGGGTATTGGCGCAGTCGGCGCGGCCATAACCGGCATTGTGCTGCTCGGTGAGTCCGCTAACCCGATGCGCCTGGCGAGTCTGGCGTTAATCGTATTGGGGATTATTGGTCTGAAACTCAGCACTCACTAACTACCAGGCTGCTGTACCCAAATAAATTTACTGACATCAAACCCTTCCCGGGTCGCGACTGCCAGCATCTCCTGTTTCACTTCGTCAGAAATGGTTGGCGTGCGGGAGAGTATCCACAGGTAGTCG  >  NC_000913/4377056‑4377318                                                                                                                                      |                                                                                                                                  tgtgGACGGGTATTGGCGCAGTCGGCGCGGCCATAACCGGCATTGTGCTGCTCGGTGAGTCCGCTAACCCGATGCGCCTGGCGAGTCTGGCGTTAATCGTATTGGGGATTATTGGTCTGAAACTCAGCACTCattaact                                                                                                                              >  2:437010/1‑139 (MQ=255) tgtgGACGGGTATTGGCGCAGTCGGCGCGGCCATAACCGGCATTGTGCTGCTCGGTGAGTCCGATAACCCGATGCGCCTGGCGAGTCTGGCGTTAATCGTATTGGGGATTATTGGTCTGAAACTCAGCACTCattaact                                                                                                                              >  2:395371/1‑139 (MQ=255)    ggACGGGTATTGGCGCAGTCGGCGCGGCCATAACCGGCATTGTGCTGCTCGGTGAGTCCGCTAACCCGATGCGCCTGGCGAGTCTGGCGTTAATCGTATT‑GGGATTATTGGTCTGAAACTCAGCACTCATTAACTACCa                                                                                                                          <  2:385518/139‑1 (MQ=255)       cGGGTATTGGCGCAGTCGGCGCGGCCATAACCGGCATTGTGCTGCTCGGTGAGTCCGCTAACCCGATGCGCCTGGCGAGTCTGGCGTTAATCGTATTGGGGATTATTGGTCTGAAACTCAGCACTCATTAACTACCAgg                                                                                                                        <  2:982319/139‑1 (MQ=255)           tATTGGCGCAGTCGGCGCGGCCATAACCGGCATTGTGCTGCTCGGTGAGTCCGCTAACCCGATGCGCCTGGCGAGTCTGGCGTTAATCGTATTGGGGATTATTGGTCTGAAACTCAGCACTCATTAACTACCAGgctgc                                                                                                                    <  2:51030/139‑1 (MQ=255)             ttGGCGCAGTCGGCGCGGCCATAACCGGCATTGTGCTGCTCGGTGAGTCCGCTAACCCGATGCGCCTGGCGAGTCTGGCGTTAATCGTATTGGGGATTATTGGTCTGAAACTCAGCACTCATTAACTACCAgg                                                                                                                        <  1:762485/133‑1 (MQ=255)             ttGGCGCAGTCGGCGCGGCCATAACCGGCATTGTGCTGCTCGGTGAGTCCGCTAACCCGATGCGCCTGGCGAGTCTGGCGTTAATCGTATTGGGGATTATTGGTCTGAAACTCAGCACTCATTAACTACCAgg                                                                                                                        >  2:762485/1‑133 (MQ=255)                    aGTCGGCGCGGCCATAACCGGCATTGTGCTGCTCGGTGAGTCCGCTAACCCGATGCGCCTGGCGAGTCTGGCGTTCATCGTATTGGGGATTATAGGTCTGAAACTCAGCACTCATTAACTACCAGGCTGCTGTACCCaa                                                                                                           >  2:864699/1‑139 (MQ=255)                    aGTCGGCGCGGCCATAACCGGCATTGTGCTGCTCGGTGAGTCCGCTAACCCGATGCGCCTGGCGAGTCTGGCGTTAATCGTATTGGGGATTATTGGTCTGAAACTCAGCACTCATTAACTACCAGGCTGCTGTACCCa                                                                                                            >  2:336207/1‑138 (MQ=255)                    aGTCGGCGCGGCCATAACCGGCATTGTGCTGCTCGGTGAGTCCGCTAACCCGATGCGCCTGGCGAGTCTGGCGTTAATCGTATTGGGGATTATTGGTCTGAAACTCAGCACTCATTAACTACCAGGCTGCTGTACCCa                                                                                                            <  1:336207/138‑1 (MQ=255)                     gTCGGCGCGGCCATAACCGGCATTGTGCTGCTCGGTGAGTCCGCTAACCCGATGCGCCTGGCGAGTCTGGCGTTAATCGTATTGGGGATTATTGGTCTGAAACTCAGCACTCATTAACTACCAGGCTGCTGTACCCaaa                                                                                                          >  1:1199596/1‑139 (MQ=255)                      tCGGCGCGGCCATAACCGGCATTGTGCTGCTCGGTGAGTCCGCTAACCCGATGCGCCTGGCGAGTCTGGCGTTAATCGTATTGGGGATTATTGGTCTGAAACTCAGCACTCATTAACTACCAGGCTGCTGTACCCaaat                                                                                                         <  1:1253316/139‑1 (MQ=255)                               ccATAACCGGCATTGTGCTGCTCGGTGAGTCCGCTAACCCGATGCGCCTGGCGAGTCTGGCGTTAATCGTATTGGGGATTATTGGTCTGAAACTCAGCACTCATTAACTACCAGGCTGCTGTACCCAAATAAATTTACt                                                                                                <  1:254642/139‑1 (MQ=255)                                         cATTGTGCTGCTCGGTGAGTCCGCTAACCCGATGCGCCTGGCGAGTCTGGCGTTAATCGTATTGGGGATTATTGGTCTGAAACTCAGCACTCATTAACTACCAGGCTGCTGTACCCAAATAAATTTACTg                                                                                               <  2:513065/130‑1 (MQ=255)                                         cATTGTGCTGCTCGGTGAGTCCGCTAACCCGATGCGCCTGGCGAGTCTGGCGTTAATCGTATTGGGGATTATTGGTCTGAAACTCAGCACTCATTAACTACCAGGCTGCTGTACCCAAATAAATTTACTg                                                                                               >  1:513065/1‑130 (MQ=255)                                         cATTGTGCTGCTCGGTGAGTCCGCTAACCCGATGCGCCTGGCGAGTCTGGCGTTAATCGTATTGGGGATTATTGGTCTGAAACTCAGCACTCATTAACTACCAGGCTGCTGTACCCAAATAAATTTACTGACATCAAAc                                                                                      >  1:163961/1‑139 (MQ=255)                                              tgctgcTCGGTGAGTCCGCTAACCCGATGCGCCTGGCGAGTCTGGCGTTAATCGTATTGGGGATTATTGGTCTGAAACTCAGCACTCATTAACTACCAGGCTGCTGTACCCAAATAAATTTACTGACATCAAACCCTTc                                                                                 >  1:227229/1‑139 (MQ=255)                                              tgctgcTCGGTGAGTCCGCTAACCCGATGCGCCTGGCGAGTCTGGCGTTAATCGTATTGGGGATTATTGGTCTGAAACTCAGCACTCATTAACTACCAGGCTGCTGTACCCAAATAAATTTACTGACATCAAACCCTTc                                                                                 >  2:240232/1‑139 (MQ=255)                                                              cGCTAACCCGATGCGCCTGGCGAGTCTGGCGTTAATCGTATTGGGGATTATTGGTCTGAAACTCAGCACTCATTAACTACCAGGCTGCTGTACCCAAATAAATTTACTGACATCAAACCCTTCCCGGGTCGCGACTGcc                                                                 <  1:240232/139‑1 (MQ=255)                                                                 tAACCCGATGCGCCTGGCGAGTCTGGCGTTAATCGTATTGGGGATTATTGGTCTGAAACTCAGCACTCATTAACTACCAGGCTGCTGTACCCAAATAAATTTACTGACATCAAACCCTTCCCGGGTCGCGACTGCCAGc                                                              <  2:227229/139‑1 (MQ=255)                                                                      cGATGCGCCTGGCGAGTCTGGCGTTAATCGTATTGGGGATTATTGGTCTGAAACTCAGCACTCATTAACTACCAGGCTGCTGTACCCAAATAAATTTACTGACATCAAACCCTTCCCGGGTCGCGACTGCCAGCAtctc                                                         >  1:13311/1‑139 (MQ=255)                                                                             ccTGGCGAGTCTGGCGTTAATCGTATTGGGGATTATTGGTCTGAAACTCAGCACTCATTAACTACCAGGCTGCTGTACCCAAATAAATTTACTGACATCAAACCCTTCCCGGGTCGCGa                                                                      >  2:223158/1‑119 (MQ=255)                                                                             ccTGGCGAGTCTGGCGTTAATCGTATTGGGGATTATTGGTCTGAAACTCAGCACTCATTAACTACCAGGCTGCTGTACCCAAATAAATTTACTGACATCAAACCCTTCCCGGGTCGCGa                                                                      <  1:223158/119‑1 (MQ=255)                                                                               tGGCGAGTCTGGCGTTAATCGTATTGGGGATTATTGGTCTGAAACTCAGCACTCATTAACTACCAGGCTGCTGTACCCAAATAAATTTACTGACATCAAACCCTTCCCGGGTCGCGACTGCCAGCATCTCCTGTTTCAc                                                >  1:1173623/1‑139 (MQ=255)                                                                                ggCGAGTCTGGCGTTAATCGTATTGGGGATTATTGGTCTGAAACTCAGCACTCATTAACTACCAGGCTGCTGTACCCAAATAAATTTACTGACATCAAACCCTTCCCGGGTCGCGACTGCCAGCATCTCCTGTTTCACt                                               <  1:1191160/139‑1 (MQ=255)                                                                                       cTGGCGTTAATCGTATTGGGGATTATTGGTCTGAAACTCAGCACTCATTAACTACCAGGCTGCTGTACCCAAATAAATTTACTGACATCAAACCCTTCCCGGGTCGCGACTGCCAGCa                                                             <  2:245828/118‑1 (MQ=255)                                                                                       cTGGCGTTAATCGTATTGGGGATTATTGGTCTGAAACTCAGCACTCATTAACTACCAGGCTGCTGTACCCAAATAAATTTACTGACATCAAACCCTTCCCGGGTCGCGACTGCCAGCa                                                             >  1:245828/1‑118 (MQ=255)                                                                                          gCGTTAATCGTATTGGGGATTATTGGTTTGAAACTCAGCACTCATTAACTACCAGGCTGCTGTACCCAAATAAATTTACTGACATCAAACCCTTCCCGGGTCGCGAc                                                                     <  2:486830/107‑1 (MQ=255)                                                                                          gCGTTAATCGTATTGGGGATTATTGGTTTGAAACTCAGCACTCATTAACTACCAGGCTGCTGTACCCAAATAAATTTACTGACATCAAACCCTTCCCGGGTCGCGAc                                                                     >  1:486830/1‑107 (MQ=255)                                                                                                    tATTGGGGATTATTGGTCTGAAACTCAGCACTCATTAACTACCAGGCTGCTGTACCCAAATAAATTTACTGACATCAAACCCTTCCCGGGTCGCGACTGCCAGCATCTCCTGTTTCACTTCGTCAGAAATGGTTGGCGt                           >  1:694984/1‑139 (MQ=255)                                                                                                       tGGGGATTATTGGTCTGAAACTCAGCACTCATTAACTACCAGGCTGCTGTACCCAAATAAATTTACTGACATCAAACCCTTCCCGGGTCGCGACTGCCAGCATCTCCTGTTTCACTTCGTCAGAAATGGTTGGCGTGCg                        <  2:1199596/139‑1 (MQ=255)                                                                                                        ggggATTATTGGTCTGAAACTCAGCACTCATTAACTACCAGGCTGCTGTACCCAAATAAATTTACTGACATCAAACCCTTCCCGGGTCGCGACTGCCAGCATCTCCTGTTTCACTTCGTCAGAAATGGTTGGCGTGCgg                       <  2:163961/139‑1 (MQ=255)                                                                                                             ttattGGTCTGAAACTCAGCACTCATTAACTACCAGGCTGCTGTACCCAAATAAATTTACTGACATCAAACCCTTCCCGGGTCGCGACTGCCAGCATCTCCTGTTTCACTTCGTCAGAAATGGTTGGCGTGCGGGagag                  <  2:13311/139‑1 (MQ=255)                                                                                                                       gAAACTCAGCACTCATTAACTACCAGGCTGCTGTACCCAAATAAATTTACTGACATCAAACCCTTCCCGGGTCGCGACTGCCAGCATCTCCTGTTTCACTTCGTCAGAAATGGTTGGCGTGCGGGAGAGTATCCACAgg        <  2:116405/139‑1 (MQ=255)                                                                                                                        aaaCTCAGCACTCATTAACTACCAGGCTGCTGTACCCAAATAAATTTACt                                                                                                <  2:453230/50‑1 (MQ=255)                                                                                                                        aaaCTCAGCACTCATTAACTACCAGGCTGCTGTACCCAAATAAATTTACt                                                                                                >  1:453230/1‑50 (MQ=255)                                                                                                                             cAGCACTCATTAACTACCAGGCTGCTGTACCCAAATAAATTTACTGACATCAAACCCTTCCCGGGTCGCGACt                                                                    >  2:1276587/1‑73 (MQ=255)                                                                                                                             cAGCACTCATTAACTACCAGGCTGCTGTACCCAAATAAATTTACTGACATCAAACCCTTCCCGGGTCGCGACt                                                                    <  1:1276587/73‑1 (MQ=255)                                                                                                                             cAGCACTCATTAACTACCAGGCTGCTGTACCCAAATAAATTTACTGACATCAAACCCTTCCCGGGTCGCGACTGCCAGCATCTCCTGTTTCACTTCGTCAGAAATGGTTGGCGTGCGGGAGAGTATCCACAGGTAGTcg  >  2:211089/1‑139 (MQ=255)                                                                                                                                      |                                                                                                                                  TGTGGACGGGTATTGGCGCAGTCGGCGCGGCCATAACCGGCATTGTGCTGCTCGGTGAGTCCGCTAACCCGATGCGCCTGGCGAGTCTGGCGTTAATCGTATTGGGGATTATTGGTCTGAAACTCAGCACTCACTAACTACCAGGCTGCTGTACCCAAATAAATTTACTGACATCAAACCCTTCCCGGGTCGCGACTGCCAGCATCTCCTGTTTCACTTCGTCAGAAATGGTTGGCGTGCGGGAGAGTATCCACAGGTAGTCG  >  NC_000913/4377056‑4377318

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