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%	202.2 / NA	62	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 (27/35);  total (27/35)

TGTGGACGGGTATTGGCGCAGTCGGCGCGGCCATAACCGGCATTGTGCTGCTCGGTGAGTCCGCTAACCCGATGCGCCTGGCGAGTCTGGCGTTAATCGTATTGGGGATTATTGGTCTGAAACTCAGCACTCACTA‑‑ACTACCAGGCTGCTGTACCCAAATAAATTTACTGACATCAAACCCTTCCCGGGTCGCGACTGCCAGCATCTCCTGTTTCACTTCGTCAGAAATGGTTGGCGTGCGGGAGAGTATCCACAGGTAGTCGC  >  NC_000913/4377056‑4377319                                                                                                                                      |                                                                                                                                     tgtgGACGGGTATTGGCGCAGTCGGCGCGGCCATAACCGGCATTGTGCTGCTCGGTGAGTCCGCTAACCCGATGCGCCTGGCGAGTCTGGCGTTAATCGTATTGGGGATTATTGGTCTGAAACTCAGCACTCatta‑‑act                                                                                                                               >  1:479505/1‑139 (MQ=255)   tgGACGGGTATTGGCGCAGTCGGCGCGGCCATAACCGGCATTGTGCTGCTCGGTGAGTCCGCTAACCCGATGCGCCTGGCGAGTCTGGCGTTAATCGTATTGGGGATTATTGGTCTGAAACTCAGCACTCATTA‑‑ACTAc                                                                                                                             <  1:113540/139‑1 (MQ=255)     gACGGGTATTGGCGCAGTCGGCGCGGCCATAACCGGTATTGTGCTGCTCGGTGAGTCCGCTAACCCGATGCGCCTGGCGAGTCTGGCGTTAATCGTATTGGGGATTATTGGTCTGAAACTCAGCACTCATTA‑‑ACTACCa                                                                                                                           <  2:1181236/139‑1 (MQ=255)     gACGGGTATTGGCGCAGTCGGCGCGGCCATAACCGGCATTGTGCTGCTCGGTGAGTCCGCTAACCCGATGCGCCTGGCGAGTCTGGCGTTAATCGTATTGGGGATTATTGGTCTGAAACTCAGCACTCATTA‑‑ACTACCa                                                                                                                           <  2:513749/139‑1 (MQ=255)      aCGGGTATTGGCGCAGTCGGCGCGGCCATAACCGGCATTGTGCTGCTCGGTGAGTCCGCTAACCCGATGCGCCTGGCGAGTCTGGCGTTAATCGTATTGGGGATTATTGGTCTGAAACTCAGCACTCATTA‑‑ACTACCAg                                                                                                                          <  2:407021/139‑1 (MQ=255)       cGGGTATTGGCGCAGTCGGCGCGGCCATAACCGGCATTGTGCTGCTCGGTGAGTCCGCTAACCCGATGCGCCTGGCGAGTCTGGCGTTAATCGTATTGGGGATTATTGGTCTGAAACTCAGCACTCATTA‑‑ACTACCAgg                                                                                                                         <  1:397241/139‑1 (MQ=255)             ttGGCGCAGTCGGCGCGGCCATAACCGGCATTGTGCTGCTCGGTGAGTCCGCTAACCCGATGCGCCTGGCGAGTCTGGCGTTAATCGTATTGGGGATTATTGGTCTGAAACTCAGCACTCATTA‑‑ACTACCAGGctgctg                                                                                                                   >  1:310014/1‑139 (MQ=255)               ggCGCAGTCGGCGCGGCCATAACCGGCATTGTGCTGCTCGGTGAGTCCGCTAACCCGATGCGCCTGGCGAGTCTGGCGTTAATCGTATTGGGGATTATTGGTCTGAAACTCAGCACTCATTA‑‑ACTACCAGGCTGCTGTa                                                                                                                 <  2:310014/139‑1 (MQ=255)                    aGTCGGCGCGGCCATAACCGGCATTGTGCTGCTCGGTGAGTCCGCTAACCCGATGCGCCTGGCGAGTCTGGCGTTAATCGTATTGGGGATTATTGGTCTGAAACTCAGCACTCATTA‑‑ACTACCAg                                                                                                                          >  2:1516191/1‑125 (MQ=255)                    aGTCGGCGCGGCCATAACCGGCATTGTGCTGCTCGGTGAGTCCGCTAACCCGATGCGCCTGGCGAGTCTGGCGTTAATCGTATTGGGGATTATTGGTCTGAAACTCAGCACTCATTA‑‑ACTACCAg                                                                                                                          <  1:1516191/125‑1 (MQ=255)                     gTCGGCGCGGCCATAACCGGCATTGTGCTGCTCGGTGAGTCCGCTAACCCGATGCGCCTGGCGAGTCTGGCGTTAATCGTATTGGGGATTATTGGTCTGAAACTCAGCACTCATTA‑‑ACTACCAgg                                                                                                                         >  1:924653/1‑125 (MQ=255)                     gTCGGCGCGGCCATAACCGGCATTGTGCTGCTCGGTGAGTCCGCTAACCCGATGCGCCTGGCGAGTCTGGCGTTAATCGTATTGGGGATTATTGGTCTGAAACTCAGCACTCATTA‑‑ACTACCAgg                                                                                                                         <  2:924653/125‑1 (MQ=255)                                cATAACCGGCATTGTGCTGCTCGGTGAGTCCGCTAACCCGATGCGCCTGGCGAGTCTGGCGTTAATCGTATTGGGGATTATTGGTCTGAAACTCAGCACTCATTA‑‑ACTACCAGGCTGCTGTACCCAAATAAATTTACTg                                                                                                <  1:108068/139‑1 (MQ=255)                                         cATTGTGCTGCACGGTGAGTCCGCTAACCCGATGCGCCTGGCGAGTCTGGCGTTAATCGTATTGGGGATTATTGGTCTGAAACTCAGCACTCATTA‑‑ACTACCAGGCTGCTGTACCCAAATAAATTTACt                                                                                                 <  2:1284324/129‑1 (MQ=255)                                         cATTGTGCTGCACGGTGAGTCCGCTAACCCGATGCGCCTGGCGAGTCTGGCGTTAATCGTATTGGGGATTATTGGTCTGAAACTCAGCACTCATTA‑‑ACTACCAGGCTGCTGTACCCAAATAAATTTACt                                                                                                 >  1:1284324/1‑129 (MQ=255)                                             gtgCTGCTCGGTGAGTCCGCTAACCCGATGCGCCTGGCGAGTCTGGCGTTAATCGTATTGGGGATTATTGGTCTGAAACTCAGCACTCATTA‑‑ACTACCAGGCTGCTGTACCCAAATAAATTTACTGACATCAAACCCtt                                                                                   >  1:13439/1‑139 (MQ=255)                                                            tCCGCTAACCCGATGCGCCTGGCGAGTCTGGCGTTAATCGTATTGGGGATTATTGGTCTGAAACTCAGCACTCATTA‑‑ACTACCAGGCTGCTGTACCCAAATAAATTTACTGACATCAAACCCTTCCCGGGTCGCGACTg                                                                    <  1:1513227/139‑1 (MQ=255)                                                             ccGCTAACCCGATGCGCCTGGCGAGTCTGGCGTTAATCGTATTGGGGATTATTGGTCTGAAACTCAGCACTCATTA‑‑ACTACCAGGCTGCTGTACCCAAATAAATTTACt                                                                                                 >  1:1399011/1‑109 (MQ=255)                                                             ccGCTAACCCGATGCGCCTGGCGAGTCTGGCGTTAATCGTATTGGGGATTATTGGTCTGAAACTCAGCACTCATTA‑‑ACTACCAGGCTGCTGTACCCAAATAAATTTACt                                                                                                 <  2:1399011/109‑1 (MQ=255)                                                                 tAACCCGATGCGCCTGGCGAGTCTGGCGTTAATCGTATTGGGGATTATTGGTCTGAAACTCAGCACTCATTA‑‑ACTACCAGGCTGCTGTACCCAAATAAATTTACTGACATCAAACCCTTCCCGGGTCGCGACTGCCAGc                                                               <  1:1595684/139‑1 (MQ=255)                                                                        aTGCGCCTGGCGAGTCTGGCGTTAATCGTATTGGGGATTATTGGTCTGAAACTCAGCACTCATTA‑‑ACTACCAGGCTGCTGTACCCAAATAAATTTACTGACATCAAACCCTTCCCGGGTCGCGACTGCCAGCAtctc                                                          >  1:257924/1‑137 (MQ=255)                                                                        aTGCGCCTGGCGAGTCTGGCGTTAATCGTATTGGGGATTATTGGTCTGAAACTCAGCACTCATTA‑‑ACTACCAGGCTGCTGTACCCAAATAAATTTACTGACATCAAACCCTTCCCGGGTCGCGACTGCCAGCAtctc                                                          <  2:257924/137‑1 (MQ=255)                                                                         tGCGCCTGGCGAGTCTGGCGTTAATCGTATTGGGGATTATTGGTCTGAAACTCAGCACTCATTAACACTACCAGGCTGCTGTACCCAAATAAATTTACTGACATCAAACCCTTCCCGGGTCGCGACTGCCAGCATCTcc                                                         >  1:115438/1‑139 (MQ=255)                                                                             ccTGGCGAGTCTGGCGTTAATCGTATTGGGGATTATTGGTCTGAAACTCAGCACTCATTA‑‑ACTACCAGGCTGCTGTACCCAAATAAATTTACTGACATCa                                                                                          >  2:127030/1‑100 (MQ=255)                                                                             ccTGGCGAGTCTGGCGTTAATCGTATTGGGGATTATTGGTCTGAAACTCAGCACTCATTA‑‑ACTACCAGGCTGCTGTACCCAAATAAATTTACTGACATCa                                                                                          <  1:127030/100‑1 (MQ=255)                                                                             ccTGGCGAGTCTGGCGTTAATCGTATTGGGGATTATTGGTCTGAAACTCAGCACTCATTA‑‑ACTACCAGGCTGCTGTACCCAAATAAATTTACTGACATCAAACCCTTCCCGGGTCGCGACTGCCAGCATCTCCTGTTTc                                                   >  1:695977/1‑139 (MQ=255)                                                                               tGGCGAGTCTGGCGTTAATCGTATTGGGGATTATTGGTCTGAAACTCAGCACTCATTA‑‑ACTACCAGGCTGCTGTACCCAAATAAATTTACTGACATCAAACCCTTCCCGGGTCGCGACTGCCAGCATCTCCTGTTTCAc                                                 <  2:13439/139‑1 (MQ=255)                                                                                 gCGAGTCTGGCGTTAATCGTATTGGGGATTATTGGTCTGAAACTCAGCACTCATTAACACTACCAGGCTGCTGTACCCAAATAAATTTACTGACATCAAACCCTTCCCGGGTCGCGACTGCCAGCATCTCCTGTTTCAc                                                 <  2:115438/139‑1 (MQ=255)                                                                                    aGTCTGGCGTTAATCGTATTGGGGATTATTGGTCTGAAACTCAGCACTCATTA‑‑ACTACCAGGCTGCTGTACCCAAATAAATTTACTGACATCAAACCCTTCCCGGGTCGCGACTGCCAGCATCTCCTGTTTCACTTCGt                                            <  1:162350/139‑1 (MQ=255)                                                                                      tCTGGCGTTAATCGTATTGGGGATTATTGGTCTGAAACTCAGCACTCATTA‑‑ACTACCAGGCTGCTGTACCCAAATAAATTTACTGACATCAAACCCTTCCCGGGTCGCGACTGCCAGCAtctc                                                          <  2:1705986/123‑1 (MQ=255)                                                                                      tCTGGCGTTAATCGTATTGGGGATTATTGGTCTGAAACTCAGCACTCATTA‑‑ACTACCAGGCTGCTGTACCCAAATAAATTTACTGACATCAAACCCTTCCCGGGTCGCGACTGCCAGCAtctc                                                          >  1:1705986/1‑123 (MQ=255)                                                                                      tCTGGCGTTAATCGTATTGGGGATTATTGGTCTGAAACTCAGCACTCATTA‑‑ACTACCAGGCTGCTGTACCCAAATAAATTTACTGACATCAAACCCTTCCCGGGTCGCGACTGCCAGCATCTCCTGTTTCACTTCGTCa                                          >  1:441780/1‑139 (MQ=255)                                                                                        tGGCGTTAATCGTATTGGGGATTATTGGTCTGAAACTCAGCACTCATTA‑‑ACTACCAGGCTGCTGTACCCAAATAAATTTACTGACATCAAACCCTTCCCGGGTCGCGACTGCCAGc                                                               <  1:424081/116‑1 (MQ=255)                                                                                        tGGCGTTAATCGTATTGGGGATTATTGGTCTGAAACTCAGCACTCATTA‑‑ACTACCAGGCTGCTGTACCCAAATAAATTTACTGACATCAAACCCTTCCCGGGTCGCGACTGCCAGc                                                               >  2:424081/1‑116 (MQ=255)                                                                                               aaTCGTATTGGGGATTATTGGTCTGAAACTCAGCACTCATTA‑‑ACTACCAGGCTGCTGTACCCAAATAAATTTACTGACATCAAACCCTTCCCGGGTCGCGACTGc                                                                   <  1:1182395/105‑1 (MQ=255)                                                                                               aaTCGTATTGGGGATTATTGGTCTGAAACTCAGCACTCATTA‑‑ACTACCAGGCTGCTGTACCCAAATAAATTTACTGACATCAAACCCTTCCCGGGTCGCGACTGc                                                                   >  2:1182395/1‑105 (MQ=255)                                                                                               aaTCGTATTGGGGATTATTGGTCTGAAACTCAGCACTCATTA‑‑ACTACCAGGCTGCTGTACCCAAATAAATTTACTGACATCAAACCCTTCCCGGGTCGCGACTGCCAGCATCTCCTGTTTCACTTCGTCAGAAATGGtt                                 >  1:719539/1‑139 (MQ=255)                                                                                                aTCGTATTGGGGATTATTGGTCTGAAACTCAGCACTCATTA‑‑ACTACCAGGCTGCTGTACCCAAATAAATTTACTGACATCAAACCCTTCCCGGGTCGCGACTGCCAGCATCTCCTGTTTCACTTCGTCAGAAATGGTTg                                >  2:1702476/1‑139 (MQ=255)                                                                                                aTCGTATTGGGGATTATTGGTCTGAAACTCAGCACTCATTA‑‑ACTACCAGGCTGCTGTACCCAAATAAATTTACTGACATCAAACCCTTCCCGGGTCGCGACTGCCAGCATCTCCTGTTTCACTTCGTCAGAAATGGTTg                                <  2:695977/139‑1 (MQ=255)                                                                                                  cGTATTGGGGATTATTGGTCTGAAACTCAGCACTCATTA‑‑ACTACCAGGCTGCTGTACCCAAATAAATTTACTGACATCAAACCCTTCCCGGGTCGCGACTGCCAGCATCTCCTGTTTCACTTCGTCAGAAATGGTTGGc                              <  1:166995/139‑1 (MQ=255)                                                                                                  cGTATTGGGGATTATTGGTCTGAAACTCAGCACTCATTA‑‑ACTACCAGGCTGCTGTACCCAAATAAATTTACTGACATCAAACCCTTCCCGGGTCGCGACTGCCAGCATCTCCTGTTTCACTTCGTCAGAAATGGTTGGc                              <  1:1561854/139‑1 (MQ=255)                                                                                                   gTATTGGGGATTATTGGTCTGAAACTCAGCACTCATTA‑‑ACTACCAGGCTGCTGTACCCAAATAAATTTACTGACATCAAACCCTTCCCGGGTCGCGACTGCCAGCATCTCCTGTTTCACTTCGTCAGAAATGGTTGGCg                             <  1:1702476/139‑1 (MQ=255)                                                                                                          ggATTATTGGTCTGAAACTCAGCACTCATTA‑‑ACTACCAGGCTGCTGTACCCAAATAAATTTACTGACATCAAACCCTTCCCGGGTCGCGa                                                                       >  2:358988/1‑90 (MQ=255)                                                                                                          ggATTATTGGTCTGAAACTCAGCACTCATTA‑‑ACTACCAGGCTGCTGTACCCAAATAAATTTACTGACATCAAACCCTTCCCGGGTCGCGa                                                                       <  1:358988/90‑1 (MQ=255)                                                                                                            attattGGTCTGAAACTCAGGACTCATTA‑‑ACTACCAgg                                                                                                                         >  1:1103530/1‑38 (MQ=255)                                                                                                            attattGGTCTGAAACTCAGCACTCATTA‑‑ACTACCAgg                                                                                                                         <  2:1103530/38‑1 (MQ=255)                                                                                                             ttattGGTCTGAAACTCAGCACTCATTA‑‑ACTACCAGGCTGCTGTACCCAAATAAATTTACTGACATCAAACCCTTCCCGGGTCGCGACTGCCAGCATCTCCTGTTTCACTTCGTCAGAAATGGTTGGCGTGCGGGagag                   <  2:441780/139‑1 (MQ=255)                                                                                                              tattGGTCTGAAACTCAGCACTCATTA‑‑ACTACCAGGCTGCTGTACCCAAATAAATTTACTGACATc                                                                                           >  1:1609882/1‑66 (MQ=255)                                                                                                              tattGGTCTGAAACTCAGCACTCATTA‑‑ACTACCAGGCTGCTGTACCCAAATAAATTTACTGACATc                                                                                           <  2:1609882/66‑1 (MQ=255)                                                                                                              tattGGTCTGAAACTCAGCACTCATTA‑‑ACTACCAGGCTGCTGTACCCAAATAAATTTACTGACATCAAACCCTTCCCGGGTCGCGACTGCCAGCATCTCCTGTTTCACTTCGTCAGa                                        <  2:622872/117‑1 (MQ=255)                                                                                                              tattGGTCTGAAACTCAGCACTCATTA‑‑ACTACCAGGCTGCTGTACCCAAATAAATTTACTGACATCAAACCCTTCCCGGGTCGCGACTGCCAGCATCTCCTGTTTCACTTCGTCAGa                                        >  1:622872/1‑117 (MQ=255)                                                                                                                   gTCTGAAACTCAGCACTCATTA‑‑ACTACCAGGCTGCTGTACCCAAATAAATTTACTGACATc                                                                                           <  2:1129718/61‑1 (MQ=255)                                                                                                                   gTCTGAAACTCAGCACTCATTA‑‑ACTACCAGGCTGCTGTACCCAAATAAATTTACTGACATc                                                                                           >  1:1129718/1‑61 (MQ=255)                                                                                                                   gTCTGAAACTCAGCACTCATTA‑‑ACTACCAGGCTGCTGTACCCAAATAAATTTACTGACATCAAACCCTTCCCGGGTCGCGACTGCCAGCATCTCCTGTTTCACTTCGTCAGAAATGGt                                  >  1:930388/1‑118 (MQ=255)                                                                                                                   gTCTGAAACTCAGCACTCATTA‑‑ACTACCAGGCTGCTGTACCCAAATAAATTTACTGACATCAAACCCTTCCCGGGTCGCGACTGCCAGCATCTCCTGTTTCACTTCGTCAGAAATGGt                                  <  2:930388/118‑1 (MQ=255)                                                                                                                   gTCTGAAACTCAGCACTCATTA‑‑ACTACCAGGCTGCTGTACCCAAATAAATTTACTGACATCAAACCCTTCCCGGGTCGCGACTGCCAGCATCTCCTGTTTCACTTCGTCAGAAATGGTTGGCGTGCGGGAGAGTATCca             >  1:1341639/1‑139 (MQ=255)                                                                                                                     cTGAAACTCAGCACTCATTA‑‑ACTACCAGGCTGCTGTACCCAAATAAATTTACTGACATCAAACCCTTCCCGGGTCGCGAc                                                                      >  2:351243/1‑80 (MQ=255)                                                                                                                     cTGAAACTCAGCACTCATTA‑‑ACTACCAGGCTGCTGTACCCAAATAAATTTACTGACATCAAACCCTTCCCGGGTCGCGAc                                                                      <  1:351243/80‑1 (MQ=255)                                                                                                                        aaaCTCAGCACTCATTA‑‑ACTACCAGGCTGCTGTACCCAAATAAATTTACTGACATCAAACCCTTCCCGGGTCGCGACTGCCAGCAtc                                                            >  1:530050/1‑87 (MQ=255)                                                                                                                        aaaCTCAGCACTCATTA‑‑ACTACCAGGCTGCTGTACCCAAATAAATTTACTGACATCAAACCCTTCCCGGGTCGCGACTGCCAGCAtc                                                            <  2:530050/87‑1 (MQ=255)                                                                                                                              aGCACTCATTA‑‑ACTACCAGGCTGCTGTACCCAAATAAATTTACTGACAt                                                                                            >  2:410507/1‑49 (MQ=255)                                                                                                                              aGCACTCATTA‑‑ACTACCAGGCTGCTGTACCCAAATAAATTTACTGACAt                                                                                            <  1:410507/49‑1 (MQ=255)                                                                                                                              aGCACTCATTA‑‑ACTACCAGGCTGCTGTACCCAAATAAATTTACTGACATCAAACCCTTCCCGGGTCGCGACTGCCAGCATCTCCTGTTTCACTTCGTCAGAAATGGTTGGCGTGCGGGAGAGTATCCACAGGTAGTcgc  >  1:527895/1‑139 (MQ=255)                                                                                                                                cactcattA‑‑ACTACCAGGCTGCTGTACCCAAATAAATTTACTGACATCAAACCCTTCCCGGGTCGCGACTGCCAGCATCTCCTGTTTCACTTCGTCAGAAATGGTTGGCGTGCGGGagag                   <  2:1448227/120‑1 (MQ=255)                                                                                                                                cactcattA‑‑ACTACCAGGCTGCTGTACCCAAATAAATTTACTGACATCAAACCCTTCCCGGGTCGCGACTGCCAGCATCTCCTGTTTCACTTCGTCAGAAATGGTTGGCGTGCGGGagag                   >  1:1448227/1‑120 (MQ=255)                                                                                                                                      |                                                                                                                                     TGTGGACGGGTATTGGCGCAGTCGGCGCGGCCATAACCGGCATTGTGCTGCTCGGTGAGTCCGCTAACCCGATGCGCCTGGCGAGTCTGGCGTTAATCGTATTGGGGATTATTGGTCTGAAACTCAGCACTCACTA‑‑ACTACCAGGCTGCTGTACCCAAATAAATTTACTGACATCAAACCCTTCCCGGGTCGCGACTGCCAGCATCTCCTGTTTCACTTCGTCAGAAATGGTTGGCGTGCGGGAGAGTATCCACAGGTAGTCGC  >  NC_000913/4377056‑4377319

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: ‑