Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,104,494	C→T	V6I (GTC→ATC)  *197* (TAG→TAA)	wbbK ← wbbJ ←	lipopolysaccharide biosynthesis protein putative lipopolysaccharide biosynthesis O‑acetyl transferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,104,494	0	C	T	100.0%	142.3 / NA	45	V6I (GTC→ATC)  *197* (TAG→TAA)	wbbK wbbJ	lipopolysaccharide biosynthesis protein putative lipopolysaccharide biosynthesis O‑acetyl transferase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (21/24);  total (21/24)

AATGGACTAATGCGATAAAACTGACATTTTCTTTATTATTAGTTGCTGCCAAAAATTTTTTCAAAATGGTAAATGGACCGCCAGTGGTAAAATTGACCGCAGAAACAACGACTATGCTTTTTCCCATAATTTGGTCTCATGATTGTATTTCTTTATGATTTTTGCTGGTACTCCCGCAATGACAGTATTTTCGGGAATAGAACCTCTAACAACAGAATTGGCGCC  >  NC_000913/2104383‑2104607                                                                                                                |                                                                                                                  aaTGGACTAATGCGATAAAACTGACATTTTCTTTATTATTAGTTGCTGCCAAAAATTTTTTCAAAATGGTAAATGGACCGCCAGTGGTAAAATTGACCGCAGAAACAACGATTATGCTTTTTCCCATAATTTGGTCTCa                                                                                        <  1:794906/139‑1 (MQ=255)  aTGGACTTATGCGATAAAACTGAAATTTTCTTTATTATTAGTTGCTGCCAAAAATTTTTTCAAAATGGTAAATGGACCGCCAGTGGTAAAATTTACCGCAGAAACAACGATTATGCTTTTTCCCATAATTTGGTCTCAt                                                                                       <  2:169497/139‑1 (MQ=255)   tGGACTAATGCGATAAAACTGACATTTTCTTTATTATTAGTTGCTGCCAAAAATTTTTTCAAAATGGTAAATGGACCGCCAGTGGTAAAATTGACCGCAGAAACAACGATTATGCTTTTTCCCATAATTTGGTCTCAt                                                                                       >  1:149571/1‑138 (MQ=255)   tGGACTAATGCGATAAAACTGACATTTTCTTTATTATTAGTTGCTGCCAAAAATTTTTTCAAAATGGTAAATGGACCGCCAGTGGTAAAATTGACCGCAGAAACAACGATTATGCTTTTTCCCATAATTTGGTCTCAt                                                                                       <  2:149571/138‑1 (MQ=255)         aaTGCGATAAAACTGACATTTTCTTTATTATTAGTTGCTGCCAAAAATTTTTTCAAAATGGTAAATGGACCGCCAGTGGTAAAATTGACCGCAGAAACAACGATTATGCTTTTTCCCATAATTTGGTCTCATGATTGTa                                                                                >  2:411958/1‑139 (MQ=255)           tGCGATAAAACTGACATTTTCTTTATTATTAGTTGCTGCCAAAAATTTTTTCAAAATGGTAAATGGACCGCCAGTGGTAAAATTGACCGCAGAAACAACGATTATg                                                                                                               <  1:69672/106‑1 (MQ=255)           tGCGATAAAACTGACATTTTCTTTATTATTAGTTGCTGCCAAAAATTTTTTCAAAATGGTAAATGGACCGCCAGTGGTAAAATTGACCGCAGAAACAACGATTATg                                                                                                               >  2:69672/1‑106 (MQ=255)           tGCGATAAAACTGACATTTTCTTTATTATTAGTTGCTGCCAAAAATTTTTTCAAAATGGTAAATGGACCGCCAGTGGTAAAATTGACCGCAGAAACAACGATTATGCTTTTTCCCATAATTTGGTCTCAt                                                                                       >  1:346220/1‑130 (MQ=255)           tGCGATAAAACTGACATTTTCTTTATTATTAGTTGCTGCCAAAAATTTTTTCAAAATGGTAAATGGACCGCCAGTGGTAAAATTGACCGCAGAAACAACGATTATGCTTTTTCCCATAATTTGGTCTCAt                                                                                       <  2:346220/130‑1 (MQ=255)               aTAAAACTGACATTTTCTTTATTATTAGTTGCTGCCAAAAATTTTTTCAAAATGGTAAATGGACCGCCAGTGGTAAAATTGACCGCAGAAACAACGATTATGCTTTTTCCCATAATTTGGTCTCAt                                                                                       >  2:486026/1‑126 (MQ=255)               aTAAAACTGACATTTTCTTTATTATTAGTTGCTGCCAAAAATTTTTTCAAAATGGTAAATGGACCGCCAGTGGTAAAATTGACCGCAGAAACAACGATTATGCTTTTTCCCATAATTTGGTCTCAt                                                                                       <  1:486026/126‑1 (MQ=255)                tAAAACTGACATTTTCTTTATTATTAGTTGCTGCCAAAAATTTTTTCAAAATGGTAAATGGACCGCCAGTGGTAAAATTGACCGCAGAAACAACGATTATGCTTTTTCCCATAATTTGGTCTCATGATTGTATTTCttt                                                                         >  2:150469/1‑139 (MQ=255)                                tttATTATTAGTTGCTGCCAAAAATTTTTTCAAAATGGTAAATGGACCGCCAGTGGTAAAATTGACCGCAGAAACAACGATTATGCTTTTTCCCATAATTTGGTCTCATGATTGTATTTCTTTATGATTTTTGCTGGTa                                                         <  1:150469/139‑1 (MQ=255)                                  tattatttGTTGCTGCCAAAAATTTTTTCAAAATGGTAAATGGACCGCCAGTGGTAAAAATGACCGCAGAAACAACGATTATGCTTTTTCCCATAAATTGGTCTCATGATTGTATTTCTTTATGATTTTTGCTGGTACt                                                       <  2:410007/139‑1 (MQ=255)                                  tattattaGTTGCTGCCAAAAATTTTTTCAAAATGGTAAATGGACCGCCAGTGGTAAAATTGACCGCAGAAACAACGATTATGCTTTTTCCCATAATTTGGTCTCATGATTGTATTTCTTTATGATTTTTGCTGGTACt                                                       >  2:836170/1‑139 (MQ=255)                                       ttaGTTGCTGCCAAAAATTTTTTCAAAATGGTAAATGGACCGCCAGTGGTAAAATTGACCGCAGAAACAACGATTATGCTTTTTCCCATAATTTGGTCTGATGATTGTATTTCTTTATGATTTTTGCTGGTACt                                                       >  1:654116/1‑134 (MQ=255)                                       ttaGTTGCTGCCAAAAATTTTTTCAAAATGGTAAATGGACCGCCAGTGGTAAAATTGACCGCAGAAACAACGATTATGCTTTTTCCCATAATTTGGTCTCATGATTGTATTTCTTTATGATTTTTGCTGGTACt                                                       <  2:654116/134‑1 (MQ=255)                                              ctgcCAAAAATTTTTTCAAAATGGTAAATGGACCGCCAGTGGTAAAATTGACCGCAGAAACAACGATTATGCTTTTTCCCATAATTTGGtct                                                                                          >  1:998601/1‑92 (MQ=255)                                              ctgcCAAAAATTTTTTCAAAATGGTAAATGGACCGCCAGTGGTAAAATTGACCGCAGAAACAACGATTATGCTTTTTCCCATAATTTGGtct                                                                                          <  2:998601/92‑1 (MQ=255)                                              ctgcCAAAAATTTTTTCAAAATGGTAAATGGACCGCCAGTGGTAAAATTGACCGCAGAAACAACGATTATGCTTTTTCCCATAATTTGGTCTCATGATTGTATTTCTTTATGATTTTTGCTGGTACTCCCGCAATGACa                                           >  2:935198/1‑139 (MQ=255)                                              ctgcCAAAAATTTTTTCAAAATGGTAAATGGACCGCCAGTGGTAAAATTGACCGCAGAAACAACGATTATGCTTTTTCCCATAATTTGGTCTCATGATTGTATTTCTTTATGATTTTTGCTGGTACTCCCGCAATGACa                                           >  2:987356/1‑139 (MQ=255)                                              ctgcCAAAAATTTTTTCAAAATGGTAAATGGACCGCCAGTGGTAAAATTGACCGCAGAAACAACGATTATGCTTTTTCCCATAATTTGGTCTCATGATTGTATTTCTTTATGATTTTTGCTGGTACTCCCGCAATGACa                                           >  1:124653/1‑139 (MQ=255)                                               tgcCAAAAATTTTTTCAAAATGGTAAATGGACCGCCAGTGGTAAAATTGACCGCAGAAACAACGATTATGCtt                                                                                                            <  1:397762/73‑1 (MQ=255)                                               tgcCAAAAATTTTTTCAAAATGGTAAATGGACCGCCAGTGGTAAAATTGACCGCAGAAACAACGATTATGCtt                                                                                                            >  2:397762/1‑73 (MQ=255)                                               tgcCAAAAATTTTTTCAAAATGGTAAATGGACCGCCAGTGGTAAAATTGACCGCAGAAACAACGATTATGCTTTTTCCCATAATTTGGt                                                                                            >  1:1279/1‑89 (MQ=255)                                               tgcCAAAAATTTTTTCAAAATGGTAAATGGACCGCCAGTGGTAAAATTGACCGCAGAAACAACGATTATGCTTTTTCCCATAATTTGGt                                                                                            <  2:1279/89‑1 (MQ=255)                                               tgcCAAAAATTTTTTCAAAATGGTAAATGGACCGCCAGTGGTAAAATTGACCGCAGAAACAACGATTATGCTTTTTCCCATAATTTGGTCTCATGATTGTATTTCTTTATGATTTTTGCTGGTACTCCCGCAATGACAg                                          >  1:498196/1‑139 (MQ=255)                                                 ccAAAAATTTTTTCAAAATGGTAAATGGACCGCCAGTGGTAAAATTGACCGCAGAAACAACGATTATGCTTTTTCCCATAATTTGGTCTCATGATTGTATTTCTTTATGATTTTTGCt                                                             <  1:553486/118‑1 (MQ=255)                                                 ccAAAAATTTTTTCAAAATGGTAAATGGACCGCCAGTGGTAAAATTGACCGCAGAAACAACGATTATGCTTTTTCCCATAATTTGGTCTCATGATTGTATTTCTTTATGATTTTTGCt                                                             >  2:553486/1‑118 (MQ=255)                                                 ccAAAAATTTTTTCAAAATGGTAAATGGACCGCCAGTGGTAAAATTGACCGCAGAAACAACGATTATGCTTTTTCCCATAATTTGGTCTCATGATTGTATTTCTTTATGATTTTTGCTGGTACTCCCGCAATGACAGTa                                        >  2:139865/1‑139 (MQ=255)                                                      aaTTTTTTCAAAATGGTAAATGGACCGCCAGTGGTAAAATTGACCGCAGAAACAACGATTATGCTTTTTCCCATAATTTGGTCTCATGATTGTATTTCTTTATGATTTTTGCTGGTACTCCCGCAATGACAGTATTTTc                                   <  1:363857/139‑1 (MQ=255)                                                             tCAAAATGGTAAATGGACCGCCAGTGGTAAAATTGACCGCAGAAACAACGATTATGCTTTTTCCCATAATTTGGTCTCATGATTGTATTTCTTTATGATTTTTGCTGGTACTCCCGCAATGACAGTATTTTCGGGAATa                            <  2:498196/139‑1 (MQ=255)                                                              cAAAATGGTAAATGGACCGCCAGTGGTAAAATTGACCGCAGAAACAACGATTATGCTTTTTCCCATAATTTGGTCTCATGATTGTATTTCTTTATGATTTTTGCTGGTACTCCCGCAATGACAGTATTTTCGGGAATAg                           <  2:355223/139‑1 (MQ=255)                                                                aaaTGGTAAATGGACCGCCAGTGGTAAAATTGACCGCAGAAACAACGATTATGCTTTTTCCCATAATTTGGTCTCATGATTGTATTTCTTTATGATTTTTGCTGGTACTCCCGCAATGACAGTATTTTCGGGAATAGaa                         <  2:660238/139‑1 (MQ=255)                                                                     gTAAATGGACCGCCAGTGGTAAAATTGACCGCAGAAACAACGATTATGCTTTTTCCCATAATTTGGTCTCATGATTGTATTTCTTTATGATTTTTGCTGGTACTCCCGCAATGACAGTATTTTCGGGAATAGAACctct                    <  1:836170/139‑1 (MQ=255)                                                                             aCCGCCAGTGGTAAAATTGACCGCAGAAACAACGATTATGCTTTTTCCCATAATTTGGTCTCATGa                                                                                     <  1:256256/66‑1 (MQ=255)                                                                             aCCGCCAGTGGTAAAATTGACCGCAGAAACAACGATTATGCTTTTTCCCATAATTTGGTCTCATGa                                                                                     >  2:256256/1‑66 (MQ=255)                                                                               cGCCAGTGGTAAAATTGACCGCAGAAACAACGATTATGCTTTTTCCCa                                                                                                     <  2:613683/48‑1 (MQ=255)                                                                               cGCCAGTGGTAAAATTGACCGCAGAAACAACGATTATGCTTTTTCCCa                                                                                                     >  1:613683/1‑48 (MQ=255)                                                                                     tGGTAAAATTGACCGCAGAAACAACGATTATGCTTTTTCCCATAATTTGGTCTCATGATTGTATTTCTTTATGATTTTTGCTGGTACTCCCGCAATGACAGTATTTTCGGGAATAGAACCTCTAACAACAGAATTGgcg    <  1:411958/139‑1 (MQ=255)                                                                                       gTAAAATTGACCGCAGAAACAACGATTATGCTTTTTCCCATAATTTg                                                                                              <  1:828613/47‑1 (MQ=255)                                                                                       gTAAAATTGACCGCAGAAACAACGATTATGCTTTTTCCCATAATTTg                                                                                              >  2:828613/1‑47 (MQ=255)                                                                                       gTAAAATTGACCGCAGAAACAACGATTATGCTTTTTCCCATAATTTGGTCTCATGATTGTATTTCTTTATGATTTTTGCTGGTACTCCCGCAATGACAGTATTTTCGGGAATAGAACCTCTAACAACAGAATTGGCGcc  <  1:962023/139‑1 (MQ=255)                                                                                                aCCGCAGAAACAACGATTATGCTTTTTCCCATAATTTGGTCTCATGATTGTATTTCTTTATGATTTTTGCTGGTACTCCCg                                                   >  2:26015/1‑81 (MQ=255)                                                                                                aCCGCAGAAACAACGATTATGCTTTTTCCCATAATTTGGTCTCATGATTGTATTTCTTTATGATTTTTGCTGGTACTCCCg                                                   <  1:26015/81‑1 (MQ=255)                                                                                                                |                                                                                                                  AATGGACTAATGCGATAAAACTGACATTTTCTTTATTATTAGTTGCTGCCAAAAATTTTTTCAAAATGGTAAATGGACCGCCAGTGGTAAAATTGACCGCAGAAACAACGACTATGCTTTTTCCCATAATTTGGTCTCATGATTGTATTTCTTTATGATTTTTGCTGGTACTCCCGCAATGACAGTATTTTCGGGAATAGAACCTCTAACAACAGAATTGGCGCC  >  NC_000913/2104383‑2104607

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