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%	134.2 / NA	42	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 (20/22);  total (20/22)

TCTTTTAACTCTTTAGCAGAATGGACTAATGCGATAAAACTGACATTTTCTTTATTATTAGTTGCTGCCAAAAATTTTTTCAAAATGGTAAATGGACCGCCAGTGGTAAAATTGACCGCAGAAACAACGACTATGCTTTTTCCCATAATTTGGTCTCATGATTGTATTTCTTTATGATTTTTGCTGGTACTCCCGCAATGACAGTATTTTCGGGAATAGAACCTCTAACAACAGAATTGGCGCCGACTACGAC  >  NC_000913/2104364‑2104616                                                                                                                                   |                                                                                                                           tctttTAACTCTTTAGCAGAATGGACTAATGCGATAAAACTGACATTTTCTTTATTATTAGTTGCTGCCAAAAATTTTTTCAAAATGGTAAATGGACCGCCACTGGTAAAATTGACCGCAGAAACAACGATTATGCttt                                                                                                                    >  2:935196/1‑139 (MQ=255)            tttAGCAGAATGGACTAATGCGATAAAACTGACATTTTCTTTATTATTAGTTGCTGCCAAAAATTTTTTCAAAATGGTAAATGGACCGCCAGTGGTAAAATTGACCGCAGAAACAACGATTATGCTTTTTCCCATAAtt                                                                                                         <  2:555499/139‑1 (MQ=255)             ttAGCAGAATGGACTAATGCGATAAAACTGACATTTTCTTTATTATTAGTTGCTGCCAAAAATTTTTTCAAAATGGTAAATGGACCGCCAGTGGTAAAATTGACCGCAGAAACAACGATTATGCTTTTTCCCATAAttt                                                                                                        <  1:337015/139‑1 (MQ=255)                       ggACTAATGCGATAAAACTGACATTTTCTTTATTATTAGTTGCTGCCAAAAATTTTTTCAAAATGGTAAATGGACCGCCAGTGGTAAAATTGACCGCAGAAACAACGATTATGCTTTTTCCCATAATTTGGTCTCATGa                                                                                              <  2:1417713/139‑1 (MQ=255)                       ggACTAATGCGATAAAACTGACATTTTCTTTATTATTAGTTGCTGCCAAAAATTTTTTCAAAATGGTAAATGGACCGCCAGTGGTAAAATTGACCGCAGAAACAACGATTATGCTTTTTCCCATAATTTGGTCTCATGa                                                                                              <  1:885814/139‑1 (MQ=255)                              tGCGATAAAACTGACATTTTCTTTATTATTAGTTGCTGCCAAAAATTTTTTCAAAATGGTAAATGGACCGCCAGTGGTAAAATTGACCGCAGAAACAACGATTATGCTTTTTCCCATAATTTGGtct                                                                                                   <  1:1096226/127‑1 (MQ=255)                              tGCGATAAAACTGACATTTTCTTTATTATTAGTTGCTGCCAAAAATTTTTTCAAAATGGTAAATGGACCGCCAGTGGTAAAATTGACCGCAGAAACAACGATTATGCTTTTTCCCATAATTTGGtct                                                                                                   >  2:1096226/1‑127 (MQ=255)                              tGCGATAAAACTGACATTTTCTTTATTATTAGTTGCTGCCAAAAATTTTTTCAAAATGGTAAATGGACCGCCAGTGGTAAAATTGACCGCAGAAACAACGATTATGCTTTTTCCCATAATTTGGTCTCATGATTGTAtt                                                                                       >  1:1402352/1‑139 (MQ=255)                                   tAAAACTGACATTTTCTTTATTATTAGTTGCTGCCAAAAATTTTTTCAAAATGGTAAATGGACCGCCAGTGGTAAAATTGACCGCAGAAACAACGATTATGCTTTTTccc                                                                                                               <  2:371956/110‑1 (MQ=255)                                   tAAAACTGACATTTTCTTTATTATTAGTTGCTGCCAAAAATTTTTTCAAAATGGTAAATGGACCGCCAGTGGTAAAATTGACCGCAGAAACAACGATTATGCTTTTTccc                                                                                                               >  1:371956/1‑110 (MQ=255)                                        cTGACATTTTCTTTATTATTAGTTGCTGCCAAAAATTTTTTCAAAATGGTAAATGGACCGCCAGTGGTAAAATTGACCGCAGAAACAACGATTATGCTTTTTCCCATAATTTGGTCTCATGATTGTATTTCTTTATGAt                                                                             <  2:1364319/139‑1 (MQ=255)                                            cATTTTCTTTATTATTAGTTGCTGCCAAAAATTTTTTCAAAATGGTAAATGGACCGCCAGTGGTAAAATTGACCGCAGAAACAACGATTATGCtt                                                                                                                     >  1:358489/1‑95 (MQ=255)                                            cATTTTCTTTATTATTAGTTGCTGCCAAAAATTTTTTCAAAATGGTAAATGGACCGCCAGTGGTAAAATTGACCGCAGAAACAACGATTATGCtt                                                                                                                     <  2:358489/95‑1 (MQ=255)                                            cATTTTCTTTATTATTAGTTGCTGCCAAAAATTTTTTCAAAATGGTAAATGGACCGCCAGTGGTAAAATTGACCGCAGAAACAACGATTATGCTTTTTCCCATAATTTGGTCTCATGATTGTATTTCTTTATGAttttt                                                                         <  2:1299505/139‑1 (MQ=255)                                                  cTTTATTATTAGTTGCTGCCAAAAATTTTTTCAAAATGGTAAATGGACCGCCAGTGGTAAAATTGACCGCAGAAACAACGATTATGCTTTTTCCCATAATTTGGTCTCATGATTGTATTTCTTTATGATTTTTGCTGGt                                                                   >  1:547987/1‑139 (MQ=255)                                                           taGTTGCTGCCAAAAATTTTTTCAAAATGGTAAATGGACCGCCAGTGGTAAAATTGACCGCAGAAACAACGATTATGCttt                                                                                                                    >  1:1178664/1‑81 (MQ=255)                                                           taGTTGCTGCCAAAAATTTTTTCAAAATGGTAAATGGACCGCCAGTGGTAAAATTGACCGCAGAAACAACGATTATGCttt                                                                                                                    <  2:1178664/81‑1 (MQ=255)                                                                gctgcCAAAAATTTTTTCAAAATGGTAAATGGACCGCCAGTGGTAAAATTGACCGCAGAAACAACGATTATGCTTTTTCCCATAATTTGGTCTCATg                                                                                               >  2:758856/1‑97 (MQ=255)                                                                gctgcCAAAAATTTTTTCAAAATGGTAAATGGACCGCCAGTGGTAAAATTGACCGCAGAAACAACGATTATGCTTTTTCCCATAATTTGGTCTCATg                                                                                               <  1:758856/97‑1 (MQ=255)                                                                gctgcCAAAAATTTTTTCAAAATGGTAAATGGACCGCCAGTGGTAAAATTGACCGCAGAAACAACGATTATGCTTTTTCCCATAATTTGGTCTCATGATTGTATTTCTTTATGATTTTTGCTGGTACTCCCGCAATGAc                                                     >  2:94010/1‑139 (MQ=255)                                                                 ctgcCAAAAATTTTTTCAAAATGGTAAATGGACCGCCAGTGGTAAAATTGACCGCAGAAACAATGATTATGCTTTTTCCCATAATTTGGt                                                                                                     >  1:957761/1‑90 (MQ=255)                                                                 ctgcCAAAAATTTTTTCAAAATGGTAAATGGACCGCCAGTGGTAAAATTGACCGCAGAAACAACGATTATGCTTTTTCCCATAATTTGGt                                                                                                     <  2:957761/90‑1 (MQ=255)                                                                 ctgcCAAAAATTTTTTCAAAATGGTAAATGGACCGCCAGTGGTAAAATTGACCGCAGAAACAACGATTATGCTTTTTCCCATAATTTGGTCTCATGATTGTATTTCTTTATGAtt                                                                            <  1:1002785/115‑1 (MQ=255)                                                                 ctgcCAAAAATTTTTTCAAAATGGTAAATGGACCGCCAGTGGTAAAATTGACCGCAGAAACAACGATTATGCTTTTTCCCATAATTTGGTCTCATGATTGTATTTCTTTATGAtt                                                                            >  2:1002785/1‑115 (MQ=255)                                                                    ccAAAAATTTTTTCAAAATGGTAAATGGACCGCCAGTGGTAAAATTGACCGCAGAAACAACGATTATGCTTTTTCCCATAATTTGGTCTCATGATTGTATTTCTTTATGATTTTTGCTGGTACTCCCGCAATGACAGTa                                                 >  2:1233923/1‑139 (MQ=255)                                                                                 cAAAATGGTAAATGGACCGCCAGTGGTAAAATTGACCGCAGAAACAACGATTATGCTTTTTCCCATAATTTGGTCTCATGATTGTATTTCTTTATGATTTTTGCTGGTACTCCCGCAATGACAGTATTTTCGGGAATAg                                    <  1:1407450/139‑1 (MQ=255)                                                                                                aCCGCCAGTGGTAAAATTGACCGCAGAAACAACGATTATGCTTTTTCCCATAATTTGGTCTCATGATTGTATTTCTTTATGATTTTTGCTGGTACTCCCGCAATGACAGTATTTTCGGGAATAGAACCTCTAACAACAg                     >  2:688282/1‑139 (MQ=255)                                                                                                 ccGCCAGTGGTAAAATTGACCGCAGAAACAACGATTATGCTTTTTCCCATAATTTGGTCTCATGATTGTATTTCTTTATGATTTTTGCTGGTACTCCCGCAATGACAGTATTTTCGGGAATAGAACCTCTAACAACAGa                    >  1:1289826/1‑139 (MQ=255)                                                                                                  cGCCAGTGGTAAAATTGACCGCAGAAACAACGATTATGCTTTTTCCCATAATTTGGTCTCATGATTGTATTTCTTTATGATTTTTGCTGGTACTCCCGCaa                                                         <  1:1248712/101‑1 (MQ=255)                                                                                                  cGCCAGTGGTAAAATTGACCGCAGAAACAACGATTATGCTTTTTCCCATAATTTGGTCTCATGATTGTATTTCTTTATGATTTTTGCTGGTACTCCCGCaa                                                         >  2:1248712/1‑101 (MQ=255)                                                                                                          gTAAAATTGACCGCAGAAACAACGATTATGCTTTTTCCCATAATTTGGTCTCATGATTGTATTTCTTTATGATTTTTGCTGGTACTCCCGCAATGACAGTATTTTCGGGAATAGAACCTCTAACAACAGAATTGGCGcc           <  2:547987/139‑1 (MQ=255)                                                                                                                   aCCGCAGAAACAACGATTATGCTTTTTCCCATAATTTGGtct                                                                                                   <  1:1125266/42‑1 (MQ=255)                                                                                                                   aCCGCAGAAACAACGATTATGCTTTTTCCCATAATTTGGtct                                                                                                   >  2:1125266/1‑42 (MQ=255)                                                                                                                   aCCGCAGAAACAACGATTATGCTTTTTCCCATAATTTGGTCTCATGa                                                                                              >  2:460068/1‑47 (MQ=255)                                                                                                                   aCCGCAGAAACAACGATTATGCTTTTTCCCATAATTTGGTCTCATGa                                                                                              <  1:460068/47‑1 (MQ=255)                                                                                                                   aCCGCAGAAACAACGATTATGCTTTTTCCCATAATTTGGTCTCATGATTGTATTTCTTTATGATTTTTGCTGGTACTCCCGCAATGACAGTATTTTCGGGAATAGAACCTCTAACAACAGAATTGGCGCCGACTACGAc  <  1:1233923/139‑1 (MQ=255)                                                                                                                    ccGCAGAAACAACGATTATGCTTTTTCCCATAATTTGGTCTCAt                                                                                                >  1:930811/1‑44 (MQ=255)                                                                                                                    ccGCAGAAACAACGATTATGCTTTTTCCCATAATTTGGTCTCAt                                                                                                <  2:930811/44‑1 (MQ=255)                                                                                                                    ccGCAGAAACAACGATTATGCTTTTTCCCATAATTTGGTCTCATGATTGTATTTCTTTATGATTTTTGCTGGTACTCCCGCAATGACa                                                    >  1:857098/1‑88 (MQ=255)                                                                                                                    ccGCAGAAACAACGATTATGCTTTTTCCCATAATTTGGTCTCATGATTGTATTTCTTTATGATTTTTGCTGGTACTCCCGCAATGACa                                                    <  2:857098/88‑1 (MQ=255)                                                                                                                       cAGAAACAACGATTATGCTTTTTCCCATAATTTGGTCTCAt                                                                                                <  2:536886/41‑1 (MQ=255)                                                                                                                       cAGAAACAACGATTATGCTTTTTCCCATAATTTGGTCTCAt                                                                                                >  1:536886/1‑41 (MQ=255)                                                                                                                                   |                                                                                                                           TCTTTTAACTCTTTAGCAGAATGGACTAATGCGATAAAACTGACATTTTCTTTATTATTAGTTGCTGCCAAAAATTTTTTCAAAATGGTAAATGGACCGCCAGTGGTAAAATTGACCGCAGAAACAACGACTATGCTTTTTCCCATAATTTGGTCTCATGATTGTATTTCTTTATGATTTTTGCTGGTACTCCCGCAATGACAGTATTTTCGGGAATAGAACCTCTAACAACAGAATTGGCGCCGACTACGAC  >  NC_000913/2104364‑2104616

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