Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,976,605	T→C	F110L (TTT→CTT)	wzxE →	O‑antigen translocase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,976,605	0	T	C	100.0%	177.8 / NA	49	F110L (TTT→CTT)	wzxE	O‑antigen translocase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (24/25);  total (24/25)

TACGTTGCCCAGTACCATGATAATCCGCAACAGCTGCGCCGCGTGGTCGGCACTTCATCAGCGATGGTACTTGGTTTCTCTACGCTGATGGCGCTGGTTTTTGTGCTGGCAGCTGCGCCAATCAGCCAGGGATTGTTTGGTAATACCGACTATCAGGGGCTGGTGCGTTTAGTGGCGCTGGTGCAAATGGGGATCGCCTGGGGCAACCTGTTACTGGCGCTGATGAAAGGCTTTCGCGATGCCGCAGGTAA  >  NC_000913/3976470‑3976720                                                                                                                                        |                                                                                                                    tACGTTGCCCAGTACCATGATAATCCGCAACAGCTGCGCCGCGTGGTCGGCACTTCATCAGCGATGGTACTTGGTTTCTCTACGCTGATGGCGCTGGTTTTTGTGCTGGCAGCTGCGCCAATCAGCCAGGGATTGCTTg                                                                                                                  >  2:340195/1‑139 (MQ=255)  aCGTTGCCCAGTACCATGATAATCCGCAACAGCTGCGCCGCGTGGTCGGCACTTCATCAGCGATGGTACTTGGTTTCTCTACGCTGATGGCGCTGGTTTTTGTGCTGGCAGCTGCGCCAATCAGCCAGGGATTGCTTgg                                                                                                                 <  1:340195/139‑1 (MQ=255)      tGCCCAGTACCATGATAATCCGCAACAGCTGCGCCGCGTGGTCGGCACTTCATCAGCGATGGTACTTGGTTTCTCTACGCTGATGGCGCTGGTTTTTGTGCTGGCAGCTGCGCCAATCAGCCAGGGATTGCTTGGTAAt                                                                                                             >  2:726115/1‑139 (MQ=255)                  tGATAATCCGCAACAGCTGCGCCGCGTGGTCGGCACTTCATCAGCGATGGTACTTGGTTTCTCTACGCTGATGGCGCTGGTTTTTGTGCTGGCAGCTGCGCCAATCAGCCAGGGATTGCTTGGTAATACCGACTATCAg                                                                                                 >  2:1043814/1‑139 (MQ=255)                  tGATAATCCGCAACAGCTGCGCCGCGTGGTCGGCACTTCATCAGCGATGGTACTTGGTTTCTCTACGCTGATGGCGCTGGTTTTTGTGCTGGCAGCTGCGCCAATCAGCCAGGGATTGCTTGGTAATACCGACTATCAg                                                                                                 >  1:771299/1‑139 (MQ=255)                       aTCCGCAACAGCTGCGCCGCGTGGTCGGCACTTCATCAGCGATGGTACTTGGTTTCTCTACGCTGATGGCGCTGGTTTTTGTGCTGGCAGCTGCGCCAATCAGCCAGGGATTGCTTGGTa                                                                                                               >  2:740880/1‑120 (MQ=255)                       aTCCGCAACAGCTGCGCCGCGTGGTCGGCACTTCATCAGCGATGGTACTTGGTTTCTCTACGCTGATGGCGCTGGTTTTTGTGCTGGCAGCTGCGCCAATCAGCCAGGGATTGCTTGGTa                                                                                                               <  1:740880/120‑1 (MQ=255)                       aTCCGCAACAGCTGCGCCGCGTGGTCGGCACTTCATCAGCGATGGTACTTGGTTTCTCTACGCTGATGGCGCTGGTTTTTGTGCTGGCAGCTGCGCCAATCAGCCAGGGATTGCTTGGTAATACCGACTATCAGGGGCt                                                                                            >  1:334130/1‑139 (MQ=255)                       aTCCGCAACAGCTGCGCCGCGTGGTCGGCACTTCATCAGCGATGGTACTTGGTTTCTCTACGCTGATGGCGCTGGTTTTTGTGCTGGCAGCTGCGCCAATCAGCCAGGGATTGCTTGGTAATACCGACTATCAGGGGCt                                                                                            >  2:540073/1‑139 (MQ=255)                          cGCAACAGCTGCGCCGCGTGGTCGGCACTTCATCAGCGATGGTACTTGGTTTCTCTACGCTGATGGCGCTGGTTTTTGTGCTGGCAGCTGCGCCAATCAGCCAGGGATTGCTTGGTAATACCGACTATCAGGGGCTGGt                                                                                         >  1:856032/1‑139 (MQ=255)                            cAACAGCTGCGCCGCGTGGTCGGCACTTCATCAGCGATGGTACTTGGTTTCTCTACGCTGATGGCGCTGGTTTTTGTGCTGGCAGCTGCGCCAATCAGCCAGGGATTGCTTGGTAATACCGa                                                                                                        <  2:309454/122‑1 (MQ=255)                            cAACAGCTGCGCCGCGTGGTCGGCACTTCATCAGCGATGGTACTTGGTTTCTCTACGCTGATGGCGCTGGTTTTTGTGCTGGCAGCTGCGCCAATCAGCCAGGGATTGCTTGGTAATACCGa                                                                                                        >  1:309454/1‑122 (MQ=255)                            cAACAGCTGCGCCGCGTGGTCGGCACTTCATCAGCGATGGTACTTGGTTTCTCTACGCTGATGGCGCTGGTTTTTGTGCTGGCAGCTGCGCCAATCAGCCAGGGATTGCTTGGTAATACCGACTATCAGGGGCTGGTGc                                                                                       >  2:46482/1‑139 (MQ=255)                                       ccgcGTGGTCGGCACTTCATCAGCGATGGTACTTGGTTTCTCTACGCTGATGGCGCTGGTTTTTGTGCTGGCAGCTGCGCCAATCAGCCAGGGATTGCTTGGTAATACCGACTATCAGGGGCTGGTGCGTTTAGTGgcg                                                                            <  1:1103143/139‑1 (MQ=255)                                        cgcgTGGTCGGCACTTCATCAGCGATGGTACTTGGTTTCTCTACGCTGATGGCGCTGGTTTTTGTGCTGGCAGCTGCGCCAATCAGCCAGGGATTGCTTGGTAATACCGACTATCAGGGGCTGGTGCGTTTAGTGgcgc                                                                           <  1:46482/139‑1 (MQ=255)                                             ggTCGGCACTTCATCAGCGATGGTACTTGGTTTCTCTACGCTGATGGCGCTGGTTTTTGTGCTGGCAGCTGCGCCAATCAGCCAGGGATTGCTTGGTAATACCGACTATCAGGGGCTGGTGCGTTTAGTGGCGCTGGTg                                                                      <  2:856032/139‑1 (MQ=255)                                                 ggCACTTCATCAGCGATGGTACTTGGTTTCTCTACGCTGATGGCGCTGGTTTTTGTGCTGGCAGCTGCGCCAATCAGCCAGGGATTGCTTGGTAATACCGACTATCAGGGGCTGGTGCGTTTAGTGGCGCTGGTGCaaa                                                                  <  1:868792/139‑1 (MQ=255)                                                  gCACTTCATCAGCGATGGTACTTGGTTTCTCTACGCTGATGGCGCTGGTTTTTGTGCTGGCAGCTGCGCCAATCAGCCAGGGATTGCTTGGTAATACCGACTATCAgg                                                                                                <  2:205735/108‑1 (MQ=255)                                                  gCACTTCATCAGCGATGGTACTTGGTTTCTCTACGCTGATGGCGCTGGTTTTTGTGCTGGCAGCTGCGCCAATCAGCCAGGGATTGCTTGGTAATACCGACTATCAgg                                                                                                >  1:205735/1‑108 (MQ=255)                                                       tcatcaGCGATGGTACTTGGTTTCTCTACGCTGATGGCGCTGGTTTTTGTGCTGGCAGCTGCGCCAATCAGCCAGGGATTGCTTGGTAATACCGACTATCAGGGGCTGGTGCGTTTAGTGGCGCTGGTGCAAATGGGGa                                                            <  2:1260703/139‑1 (MQ=255)                                                         atcaGCGATGGTACTTGGTTTCTCTACGCTGATGGCGCTGGTTTTTGTGCTGGCAGCTGCGCCAATCAGCCAGGGATTGCTTGGTAATACCGACTATCAGGGGCTGGt                                                                                         >  1:94202/1‑108 (MQ=255)                                                         atcaGCGATGGTACTTGGTTTCTCTACGCTGATGGCGCTGGTTTTTGTGCTGGCAGCTGCGCCAATCAGCCAGGGATTGCTTGGTAATACCGACTATCAGGGGCTGGt                                                                                         <  2:94202/108‑1 (MQ=255)                                                                 tGGTACTTGGTTTCTCTACGCTGATGGCGCTGGTTTTTGTGCTGGCAGCTGCGCCAATCAGCCAGGGATTGCTTGGTa                                                                                                               <  1:1121029/78‑1 (MQ=255)                                                                 tGGTACTTGGTTTCTCTACGCTGATGGCGCTGGTTTTTGTGCTGGCAGCTGCGCCAATCAGCCAGGGATTGCTTGGTa                                                                                                               >  2:1121029/1‑78 (MQ=255)                                                                  ggTACTTGGTTTCTCTACGCTGATGGCGCTGGTTTTTGTGCTGGCAGCTGCGCCAATCAGCCAGGGATTGCTTGGTAATACCGACTATCAGGGGCTGGTGCGTTTAGTGGCGCTGGTGCAAATGGGGATCGCCTGGGGc                                                 >  1:883908/1‑139 (MQ=255)                                                                            ttCTCTACGCTGATGGCGCTGGTTTTTGTGCTGGCAGCTGCGCCAATCAGCCAGGGATTGCTTGGTAATACCGACTATCAGGGGCTGGTGCGTTTAGTGGCGCTGGTGCAAATGGGGATCGCCTGGGGCAACCTGTTAc                                       <  2:771299/139‑1 (MQ=255)                                                                             tctcTACGCTGATGGCGCTGGTTTTTGTGCTGGCAGCTGCGCCAATCAGCCAGGGATTGCTTGGTAATACCGACTATCAGGGGCTGGt                                                                                         <  2:342032/88‑1 (MQ=255)                                                                             tctcTACGCTGATGGCGCTGGTTTTTGTGCTGGCAGCTGCGCCAATCAGCCAGGGATTGCTTGGTAATACCGACTATCAGGGGCTGGt                                                                                         >  1:342032/1‑88 (MQ=255)                                                                              ctctACGCTGATGGCGCTGGTTTTTGTGCTGGCAGCTGCGCCAATCAGCCAGGGATTGCTTGGTAATACCGACTATCAGGGGCTGGTGCGTTTAGTGGCGCTGGTGCAAATGGGGATCGCCTGGGGCAACCTGTTACTg                                     >  1:1030876/1‑139 (MQ=255)                                                                                        aTGGCGCTGGTTTTTGTGCTGGCAGCTGCGCCAATCAGCCAGGGATTGCTTGGTAATACCGACTATCAGGGGCTGGTGCGTTTAGTGGCGCTGGTGCAAATGGGGATCGCCTGGGGCAACCTGTTACTGGCGCtgatga                           >  1:701243/1‑139 (MQ=255)                                                                                        aTGGCGCTGGTTTTTGTGCTGGCAGCTGCGCCAATCAGCCAGGGATTGCTTGGTAATACCGACTATCAGGGGCTGGTGCGTTTAGTGGCGCTGGTGCAAATGGGGATCGCCTGGGGCAACCTGTTACTGGCGCtgatga                           >  2:212196/1‑139 (MQ=255)                                                                                             gcTGGTTTTTGTGCTGGCAGCTGCGCCAATCAGCCAGGGATTGCTTGGTAATACCGACTATCAGGGGCTGGTGCGTTTAGTGGCGCTGGTGCAAATGGGGATCGCCTGGGGCaa                                               >  1:28056/1‑114 (MQ=255)                                                                                             gcTGGTTTTTGTGCTGGCAGCTGCGCCAATCAGCCAGGGATTGCTTGGTAATACCGACTATCAGGGGCTGGTGCGTTTAGTGGCGCTGGTGCAAATGGGGATCGCCTGGGGCaa                                               <  2:28056/114‑1 (MQ=255)                                                                                                  tttttGTGCTGGCAGCTGCGCCAATCAGCCAGGGATTGCTTGGTAATACCGACTATCAGGGGCTGGTGCGTTTAGTGGCGCTGGTGCAAATGGGGATCGCCTGGGGCAACCTGTTACTGgcgc                                 >  2:316486/1‑123 (MQ=255)                                                                                                  tttttGTGCTGGCAGCTGCGCCAATCAGCCAGGGAGTGCTTGGTAATACCGACTATCAGGGGCTGGTGCGTTTAGTGGCGCTGGTGCAAATGGGGATCGCCTGGGGCAACCTGTTACTGgcgc                                 <  1:316486/123‑1 (MQ=255)                                                                                                    tttGTGCTGGCAGCTGCGCCAATCAGCCAGGGATTGCTTGGTAATACCGACTATCAGGGGCTGGTGCGTTTAGTGGCGCTGGTGCAAATGGGGATCGCCTGGGGCAACCTGTTACTGGCGCTGATGAAAGGCTTTcgcg               <  2:334130/139‑1 (MQ=255)                                                                                                    tttGTGCTGGCAGCTGCGCCAATCAGCCAGGGATTGCTTGGTAATACCGACTATCAGGGGCTGGTGCGTTTAGTGGCGCTGGTGCAAATGGGGATCGCCTGGGGCAACCTGTTACTGGCGCTGATGAAAGGCTTTcgcg               <  1:1043814/139‑1 (MQ=255)                                                                                                    tttGTGCTGGCAGCTGCGCCAATCAGCCAGGGATTGCTTGGTAATACCGACTATCAGGGGCTGGTGCGTTTAGTGGCGCTGGTGCAAATGGGGATCGCCTGGGGCAACCTGTTACTGGCGCTGATGAAAGGCTTTcgcg               <  2:701243/139‑1 (MQ=255)                                                                                                        tgCTGGCAGCTGCGCCAATCAGCCAGGGATTGCTTGGTAATACCGACTATCAGGGGCTGGTGCGTTTAGTGGCGCTGGTGCAAATGGGGATCGCCTGGGGCAACCTGTTACTGGCGCTGATGAAAGGCTTTCGCGATGc           <  2:883908/139‑1 (MQ=255)                                                                                                          cTGGCAGCTGCGCCAATCAGCCAGGGATTGCTTGGTAATACCGACTATCAGGGGCTGGTGCGTTTAGTGGCGCTGGTGCAAATGGGGATCGCCTGGGGCAACCTGTTACTGGCGCTGATGAAAGGCTTTCGCGATGCCg         <  1:726115/139‑1 (MQ=255)                                                                                                                 cTGCGCCAATCAGCCAGGGATTGCTTGGTAATACCGACTATCAGGGGCTGGTGCGTTTAGTGGCGCTGGTGCAAATGGGGATCGCCTGGGGCAACCTGTTACTGGCGCTGATGAAAGGCTTTCGCGATGCCGCAGGTaa  <  1:834087/139‑1 (MQ=255)                                                                                                                  tGCGCCAATCAGCCAGGGATTGCTTGGTAATACCGACTATCAGGGGCTGGTGCGTTTAGTgg                                                                              <  1:946855/62‑1 (MQ=255)                                                                                                                  tGCGCCAATCAGCCAGGGATTGCTTGGTAATACCGACTATCAGGGGCTGGTGCGTTTAGTgg                                                                              >  2:946855/1‑62 (MQ=255)                                                                                                                  tGCGCCAATCAGCCAGGGATTGCTTGGTAATACCGACTATCAGGGGCTGGTGCGTTTAGTgg                                                                              >  2:124242/1‑62 (MQ=255)                                                                                                                  tGCGCCAATCAGCCAGGGATTGCTTGGTAATACCGACTATCAGGGGCTGGTGCGTTTAGTgg                                                                              <  1:124242/62‑1 (MQ=255)                                                                                                                         aTCAGCCAGGGATTGCTTGGTAATACCGACTATCAGGGGCTGGTGCGTTTAGTGGCGCTGGTGCAAATGGGGa                                                            >  1:1008305/1‑73 (MQ=255)                                                                                                                         aTCAGCCAGAGATTGCTTGGTAATACCGACTATCAGGGGCTGGTGCGTTTAGTGGCGCTGGTGCAAATGGGGa                                                            <  2:1008305/73‑1 (MQ=255)                                                                                                                                  ggATTGCTTGGTAATACCGACTATCAGGGGCTGGTGCGTTTAGTGGCGCt                                                                          <  2:1023450/50‑1 (MQ=255)                                                                                                                                  ggATTGCTTGGTAATACCGACTATCAGGGGCTGGTGCGTTTAGTGGCGCt                                                                          >  1:1023450/1‑50 (MQ=255)                                                                                                                                        |                                                                                                                    TACGTTGCCCAGTACCATGATAATCCGCAACAGCTGCGCCGCGTGGTCGGCACTTCATCAGCGATGGTACTTGGTTTCTCTACGCTGATGGCGCTGGTTTTTGTGCTGGCAGCTGCGCCAATCAGCCAGGGATTGTTTGGTAATACCGACTATCAGGGGCTGGTGCGTTTAGTGGCGCTGGTGCAAATGGGGATCGCCTGGGGCAACCTGTTACTGGCGCTGATGAAAGGCTTTCGCGATGCCGCAGGTAA  >  NC_000913/3976470‑3976720

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