Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,148,681	C→T	G529S (GGC→AGC)	eptC ←	LPS heptose I phosphoethanolamine transferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,148,681	0	C	T	100.0%	178.5 / NA	56	G529S (GGC→AGC)	eptC	LPS heptose I phosphoethanolamine transferase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (25/31);  total (25/31)

ATAGGGCAGTGTGTCGTAATCGATCAGTGCGTTTTTCTTATACGGGTTACCAATCCAGCGGGTAGTTTCTTTGAACTGCGGATTCACCACTGAACGGGTTGGATCGTAACCGTCGTAAGATAAGCCCGCCAAATCTGACCAGGTGTGGATCAGTTCCGCCAGGCTGTATTTACGATCAACATCCTGCGAGAAATCACGGGGATGAGTCGCTTGCCATTTTTCTGACGTCCACAGCAGGAACGGAATGGT  >  NC_000913/4148556‑4148804                                                                                                                              |                                                                                                                            aTAGGGCAGTGTGTCGTAATCGATCAGTGCGTTTTTCTTATACGGGTTACCAATCCAGCGGGTAGTTTCTTTGAACTGCGGATTCACCACTGAACGGGTTGGATCGTAACCGTCGTAAGATAAGCTCGCCAAATCTGAc                                                                                                                <  2:763359/139‑1 (MQ=255)          tgtgtCGTAATCGATCAGTGCGTTTTTCTTATACGGGTTACCAATCCAGCGGGTAGTTTCTTTGAACTGCGGATTCACCACTGAACGGGTTGGATCGTAACCGTCGTAAGATAAGCTCGCCAAATCTGACCAGGTGTgg                                                                                                       >  2:46862/1‑139 (MQ=255)          tgtgtCGTAATCGATCAGTGCGTTTTTCTTATACGGGTTACCAATCCAGCGGGTAGTTTCTTTGAACTGCGGATTCACCACTGAACGGGTTGGATCGTAACCGTCGTAAGATAAGCTCGCCAAATCTGACCAGGTGTgg                                                                                                       >  1:705382/1‑139 (MQ=255)           gtgtCGTAATCGATCAGTGCGTTTTTCTTATACGGGTTACCAATCCAGCGGGTAGTTTCTTTGAACTGCGGATTCACCACTGAACGGGTTGGATCGTAACCGTCGTAAGATAAGCTCGCCAAATCTGACCAGGTGTGGa                                                                                                      <  1:543453/139‑1 (MQ=255)            tgtCGTAATCGATCAGTGCGTTTTTCTTATACGGGTTACCAATCCAGCGGGTAGTTTCTTTGAACTGCGGATTCACCACTGAACGGGTTGGATCGTAACCGTCGTAAGATAAGCTCGCCAAATCTGACCAGGTGTGGAt                                                                                                     <  1:926879/139‑1 (MQ=255)            tgtCGTAATCGATCAGTGCGTTTTTCTTATACGGGTTACCAATCCAGCGGGTAGTTTCTTTGAACTGCGGATTCACCACTGAACGGGTTGGATCGTAACCGTCGTAAGATAAGCTCGCCAAATCTGACCAGGTGTGGAt                                                                                                     >  1:752320/1‑139 (MQ=255)               cGTAATCGATCAGTGCGTTTTTCTTATACGGGTTACCAATCCAGCGGGTAGTTTCTTTGAACTGCGGATTCACCACTGAACGGGTTGGATCGTAACCGTCGTAAGATAAGCTCGCCAAATCTGACCAGGTGTGGATCAg                                                                                                  <  2:705382/139‑1 (MQ=255)               cGTAATCGATCAGTGCGTTTTTCTTATACGGGTTACCAATCCAGCGGGTAGTTTCTTTGAACTGCGGATTCACCACTGAACGGGTTGGATCGTAACCGTCGTAAGATAAGCTCGCCAAATCTGACCAGGTGTGGATCAg                                                                                                  >  2:63720/1‑139 (MQ=255)                gTAATCGATCAGTGCGTTTTTCTTATACGGGTTACCAATCCAGCGGGTAGTTTCTTTGAACTGCGGATTCACCACTGAACGGGTTGGATCGTAACCGTCGTAAGATAAGCTCGCCAAATCTGACCAGGTGTGGATCAGt                                                                                                 <  2:752320/139‑1 (MQ=255)                    tCGATCAGTGCGTTTTTCTTATACGGGTTACCAATCCAGCGGGTAGTTTCTTTGAACTGCGGATTCACCACTGAACGGGTTGGATCGTAACCGTCGTAAGATAAGCTCGCCAAATCTGACCAGGTGTGGATCAGTTCCg                                                                                             >  2:151162/1‑139 (MQ=255)                       aTCAGTGCGTTTTTCTTATACGGGTTACCAATCCAGCGGGTAGTTTCTTTGAACTGCGGATTCACCACTGAACGGGTTGGATCGTAACCGTCGTAAGATAAGCTCGCCAAATCTGACCAGGTGTGGATCAGTTCCGCCa                                                                                          <  1:347977/139‑1 (MQ=255)                       aTCAGTGCGTTTTTCTTATACGGGTTACCAATCCAGCGGGTAGTTTCTTTGAACTGCGGATTCACCACTGAACGGGTTGGATCGTAACCGTCGTAAGATAAGCTCGCCAAATCTGACCAGGTGTGGATCAGTTCCGCCa                                                                                          <  2:561367/139‑1 (MQ=255)                        tCAGTGCGTTTTTCTTATACGGGTTACCAATCCAGCGGGTAGTTTCTTTGAACTGCGGATTCACCACTGAACGGGTTGGATCGTAACCGTCGTAAGATAAGCTCGCCAAATCTGAcc                                                                                                               >  1:229424/1‑117 (MQ=255)                        tCAGTGCGTTTTTCTTATACGGGTTACCAATCCAGCGGGTAGTTTCTTTGAACTGCGGATTCACCACTGAACGGGTTGGATCGTAACCGTCGTAAGATAAGCTCGCCAAATCTGAcc                                                                                                               <  2:229424/117‑1 (MQ=255)                        tCAGTGCGTTTTTCTTATACGGGTTACCAATCCAGCGGGTAGTTTCTTTGAACTGCGGATTCACCACTGAACGGGTTGGATCGTAACCGTCGTAAGATAAGCTCGCCAAATCTGACCAGGTGTGGATCAGTTCCGCCAg                                                                                         >  1:178296/1‑139 (MQ=255)                         cAGTGCGTTTTTCTTATACGGGTTACCAATCCAGCGGGTAGTTTCTTTGAACTGCGGATTCACCACTGAACGGGTTGGATCGTAACCGTCGTAAGATAAGCTCGCCAAATCTGACCAGGTGTGGATCAGTTCCGCCAgg                                                                                        <  2:178296/139‑1 (MQ=255)                          aGTGCGTTTTTCTTATACGGGTTACCAATCCAGCGGGTAGTTTCTTTGAACTGCGGATTCACCACTGAACGGGTTGGATCGTAACCGTCGTAAGATAAGCTCGCCAAATCTGACCAGGTGTGGATCAGTTCCGCCAGGc                                                                                       <  2:325711/139‑1 (MQ=255)                            tGCGTTTTTCTTATACGGGTTACCAATCCAGCGGGTAGTTTCTTTGAACTGCGGATTCACCACTGAACGGGTTGGATCGTAACCGTCGTAAGATAAGCTCGCCAAATCTGACCAGGTGTGGATCAGTTCCGCCAGGCTg                                                                                     <  1:63720/139‑1 (MQ=255)                                      ttATACGGGTTACCAATCCAGCGGGTAGTTTCTTTGAACTGCGGATTCACCACTGAACGGGTTGGATCGTAACCGTCGTAAGATAAGCTCGCCAAATCTGACCAGGTGTGGATCAGTTCCGCCAGGCTGTATTTACGAt                                                                           >  2:203033/1‑139 (MQ=255)                                                        cAGCGGGTAGTTTCTTTGAACTGCGGATTCACCACTGAACGGGTTGGATCGTAACCGTCGTAAGATAAGCTCGCCAAATCTGAcc                                                                                                               <  2:319443/85‑1 (MQ=255)                                                        cAGCGGGTAGTTTCTTTGAACTGCGGATTCACCACTGAACGGGTTGGATCGTAACCGTCGTAAGATAAGCTCGCCAAATCTGAcc                                                                                                               >  1:319443/1‑85 (MQ=255)                                                           cGGGTAGTTTCTTTGAACTGCGGATTCACCACTGAACGGGTTGGATCGTAACCGTCGTAAGATAAGCTCGCCAAATCTGACCAGGTGTGGATCAGTTCCGCCAGGCTGTATTTACGATCAACATCCTGCGAGAAATCAc                                                      >  1:491432/1‑139 (MQ=255)                                                           cGGGTAGTTTCTTTGAACTGCGGATTCACCACTGAACGGGTTGGATCGTAACCGTCGTAAGATAAGCTCGCCAAATCTGACCAGGTGTGGATCAGTTCCGCCAGGCTGTATTTACGATCAACATCCTGCGAGAAATCAc                                                      >  2:592324/1‑139 (MQ=255)                                                            gggTAGTTTCTTTGAACTGCGGATTCACCACTGAACGGGTTGGATCGTAACCGTCGTAAGATAAGCTCGCCAAATCTGACCAGGTGTGGATCAGTTCCGCCAGGCTGTATTTACGATCAACATCCTGCGAGAAATCACg                                                     >  1:481496/1‑139 (MQ=255)                                                             ggTAGTTTCTTTGAACTGCGGATTCACCACTGAACGGGTTGGATCGTAACCGTCGTAAGATAAGCTCGCCAAATCTGACCAGGTGTGGATCAGTTCCGCCAGGCTGTATTTACGATCAACATCCTGCGAGAAATCACgg                                                    <  1:592324/139‑1 (MQ=255)                                                                    tCTTTGAACTGCGGATTCACCACTGAACGGGTTGGATCGTAACCGTCGTAAGATAAGCTCGCCAAATCTGACCAGGTGTGGATCAGTTCCGCCAGGCTGTATTTACGATCAACATCCTGCGAGAAATCACGGGGATGAg                                             >  1:937661/1‑139 (MQ=255)                                                                       ttGAACTGCGGATTCACCACTGAACGGGTTGGATCGTAACCGTCGTAAGATAAGCTCGCCAAATCTGACCAGGTGTGGATCAGTTCCGCCAGGCTGTATTTACGATCAACATCCTGCGAGAAATCACGGGGATGAGTCg                                          >  1:792727/1‑139 (MQ=255)                                                                            cTGCGGATTCACCACTGAACGGGTTGGATCGTAACCGTCGTAAGATAAGCTCGCCAAATCTGACCAGGTGTGGATCAGTTCCGCCAGGCTGTATTTACGATCAACATCCTGCGAGAAATCACGGGGATGAGTCGCTTGc                                     <  1:46862/139‑1 (MQ=255)                                                                               cGGATTCACCACTGAACGGGTTGGATCGTAACCGTCGTAAGATAAGCTCGCCAAATCTGACCAGGTGTGGATCAGTTCCGCCAGGCTGTAttt                                                                                >  1:909671/1‑93 (MQ=255)                                                                               cGGATTCACCACTGAACGGGTTGGATCGTAACCGTCGTAAGATAAGCTCGCCAAATCTGACCAGGTGTGGATCAGTTCCGCCAGGCTGTAttt                                                                                <  2:909671/93‑1 (MQ=255)                                                                               cGGATTCACCACTGAACGGGTTGGATCGTAACCGTCGTAAGATAAGCTCGCCAAATCTGACCAGGTGTGGATCAGTTCCGCCAGGCTGTATTTACGATCAACATCCTGCg                                                               <  1:140973/110‑1 (MQ=255)                                                                               cGGATTCACCACTGAACGGGTTGGATCGTAACCGTCGTAAGATAAGCTCGCCAAATCTGACCAGGTGTGGATCAGTTCCGCCAGGCTGTATTTACGATCAACATCCTGCg                                                               >  2:140973/1‑110 (MQ=255)                                                                               cGGATTCACCACTGAACGGGTTGGATCGTAACCGTCGTAAGATAAGCTCGCCAAATCTGACCAGGTGTGGATCAGTTCCGCCAGGCTGTATTTACGATCAACATCCTGCGAGAAATCACGGGGATGAGTCGCTTGCCAt                                  >  2:711837/1‑139 (MQ=255)                                                                               cGGATTCACCACTGAACGGGTTGGATCGTAACCGTCGTAAGATAAGCTCGCCAAATCTGACCAGGTGTGGATCAGTTCCGCCAGGCTGTATTTACGATCAACATCCTGCGAGAAATCACGGGGATGAGTCGCTTGCCAt                                  >  2:628343/1‑139 (MQ=255)                                                                               cGGATTCACCACTGAACGGGTTGGATCGTAACCGTCGTAAGATAAGCTCGCCAAATCTGACCAGGTGTGGATCAGTTCCGCCAGGCTGTATTTACGATCAACATCCTGCGAGAAATCACGGGGATGAGTCGCTTGCCAt                                  >  2:406669/1‑139 (MQ=255)                                                                                   ttCACCACTGAACGGGTTGGATCGTAACCGTCGTAAGATAAGCTCGCCAAATCTGACCAGGTGTGGATCAGTTCCGCCAGGCTGTATTTACGATCAACATCCTGCGAGAAATCACGGGGATGAGTCGCTTGCCAttttt                              <  2:937661/139‑1 (MQ=255)                                                                                      accacTGAACGGGTTGGATCGTAACCGTCGTAAGATAAGCTCGCCAAATCTGACCAGGTGTGGATCAGTTCCGCCAGGCTGTATTTACGATCAACATCCTGCGAGAAATCACGGGGATGAGTCGCTTGCCATTTTTCTg                           <  1:711837/139‑1 (MQ=255)                                                                                       ccacTGAACGGGTTGGATCGTAACCGTCGTAAGATAAGCTCGCCAAATCTGACCAGGTGTGGATCAGtt                                                                                                <  2:117783/69‑1 (MQ=255)                                                                                       ccacTGAACGGGTTGGATCGTAACCGTCGTAAGATAAGCTCGCCAAATCTGACCAGGTGTGGATCAGtt                                                                                                >  1:117783/1‑69 (MQ=255)                                                                                             aaCGGGTTGGATCGTAACCGTCGTAAGATAAGCTCGCCAAAt                                                                                                                     >  1:338223/1‑42 (MQ=255)                                                                                             aaCGGGTTGGATCGTAACCGTCGTAAGATAAGCTCGCCAAAt                                                                                                                     <  2:338223/42‑1 (MQ=255)                                                                                               cGGGTTGGGTCGTAACCGTCGTAAGATAAGCTCGCCAAATCTGACCAGGTGTGGATCAGTTCCGcc                                                                                           <  1:722037/66‑1 (MQ=255)                                                                                               cGGGTTGGATCGTAACCGTCGTAAGATAAGCTCGCCAAATCTGACCAGGTGTGGATCAGTTCCGcc                                                                                           >  2:722037/1‑66 (MQ=255)                                                                                                   ttGGATCGTAACCGTCGTAAGATAAGCTCGCCAAATCTGACCAGGTGTGGATCAGTTCCGCCAGGCTGTATTTACGATCAACATCCTGCGAGAAATCACGGGGATGAGTCGCTTGCCATTTTTCTGACGTCCAcagcag              <  2:792727/139‑1 (MQ=255)                                                                                                    tGGATCGTAACCGTCGTAAGATAAGCTCGCCAAATCTGACCAGGTGTGGATCAGTTCCGCCAGGCTGTCTTTACGATCAACATCCTGCGAGAAATCACGGGGATGAGTCGCTTGCCATTTTTCTGACGTCCACAGCAgg             >  2:628872/1‑139 (MQ=255)                                                                                                      gATCGTAACCGTCGTAAGATAAGCTCGCCAAATCTGACCAGGTGTGGATCAGTTCCGCCAGGCTGTATTTACGATCAACATCCTGCGAGAAATCACGGGGATGAGTCGCTTGCCATTTTTCTGACGTCCACAGCAGGaa           <  1:151162/139‑1 (MQ=255)                                                                                                          gTAACCGTCGTAAGATAAGCTCGCCAAATCTGACCAGGTGTGGATCAGTTCCGCCAGGCTGTATTTACGATCAACATCCTGCGAGAAATCACGGGGATGAGTCGCTTGCCATTTTTCTGACGTCCACAGCAGGAACGGa       <  1:628872/139‑1 (MQ=255)                                                                                                          gTAACCGTCGTAAGATAAGCTCGCCAAATCTGACCAGGTGTGGATCAGTTCCGCCAGGCTGTATTTACGATCAACATCCTGCGAGAAATCACGGGGATGAGTCGCTTGCCATTTTTCTGACGTCCACAGCAGGAACGGa       <  1:203033/139‑1 (MQ=255)                                                                                                           tAACCGTCGTAAGATAAGCTCGCCAAATCTGACCAGGTGTGGATCAGTTCCGCCAGGCTGTATTTACGATCAACATCCTGCGAGAAATCACGGGGATGAGTCGCTTGCCATTTTTCTGACGTCCACAGCAGGAACGGaa      <  2:343441/139‑1 (MQ=255)                                                                                                           tAACCGTCGTAAGATAAGCTCGCCAAATCTGACCAGGTGTGGATCAGTTCCGCCAGGCTGTATTTACGATCAACATCCTGCGAGAAATAACGGGGATGAGTCGCTTGCCATTTTTCTGACGTCCACAGCAGGAACGGaa      <  1:710845/139‑1 (MQ=255)                                                                                                            aaCCGTCGTAAGATAAGCTCGCCAAATCTGACCAGGTGTGGATCAGTTCCGCCAGGCTGTATTTACGATCAACATCCTGCGAGAAATCACGGGGATGAGTCGCTTGCCATTTTTCTGACGTCCACAGCAGGAACGGAAt     <  2:481496/139‑1 (MQ=255)                                                                                                               cgtcgtAAGATAAGCTCGCCAAATCTGACCAGGTGTGGATCAGTTCCGCCAGGCTGTATTTACGATCAACATCCTGCGAGAAATCACGGGGATGAGTCGCTTGCCATTTTTCTGACGTCCACAGCAGGAACGGAATGgt  <  2:542380/139‑1 (MQ=255)                                                                                                                   gtAAGATAAGCTCGCCAAATCTGACCAGGTGTGGATCAGTTCCGCCAgg                                                                                        >  1:468143/1‑49 (MQ=255)                                                                                                                   gtAAGATAAGCTCGCCAAATCTGACCAGGTGTGGATCAGTTCCGCCAgg                                                                                        <  2:468143/49‑1 (MQ=255)                                                                                                                        aTAAGCTCGCCAAATCTGACCAGGTGTGGATCAGTTCCGCCAGGCTGTATTTAc                                                                              >  1:308789/1‑54 (MQ=255)                                                                                                                        aTAAGCTCGCCAAATCTGACCAGGTGTGGATCAGTTCCGCCAGGCTGTATTTAc                                                                              <  2:308789/54‑1 (MQ=255)                                                                                                                              |                                                                                                                            ATAGGGCAGTGTGTCGTAATCGATCAGTGCGTTTTTCTTATACGGGTTACCAATCCAGCGGGTAGTTTCTTTGAACTGCGGATTCACCACTGAACGGGTTGGATCGTAACCGTCGTAAGATAAGCCCGCCAAATCTGACCAGGTGTGGATCAGTTCCGCCAGGCTGTATTTACGATCAACATCCTGCGAGAAATCACGGGGATGAGTCGCTTGCCATTTTTCTGACGTCCACAGCAGGAACGGAATGGT  >  NC_000913/4148556‑4148804

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