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%	202.7 / NA	64	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 (33/31);  total (33/31)

TCTTTTAACTCTTTAGCAGAATGGACTAATGCGATAAAACTGACATTTTCTTTATTATTAGTTGCTGCCAAAAATTTTTTCAAAATGGTAAATGGACCGCCAGTGGTAAAATTGACCGCAGAAACAACGACTATGCTTTTTCCCATAATTTGGTCTCATGATTGTATTTCTTTATGATTTTTGCTGGTACTCCCGCAATGACAGTATTTTCGGGAATAGAACCTCTAACAACAGAATTGGCGCCGACTACGA  >  NC_000913/2104364‑2104615                                                                                                                                   |                                                                                                                          tctttTAACTCTTTAGCAGAATGGACTAATGCGATAAAACTGACATTTTCTTTATTATTAGTTGCTGCCAAAAATTTTTTCAAAATGGTAAATGGACCGCCAGTGGTAAAATTGACCGCAGAAACAACGATTATGCttt                                                                                                                   >  2:1247951/1‑139 (MQ=255) tctttTAACTCTTTAGCAGAATGGACTAATGCGAT‑AAACTGACATTTTCTTTATTATTAGTTGCTGCCAAAAATTTTTTCAAAATGGTAAATGGACCGCCAGTGGTAAAATTGACCGCAGAAACAACGATTATGCtttt                                                                                                                  >  1:419619/1‑139 (MQ=255)            tttAGCAGAATGGACTGATGCGATAAAACTGACATTTTCTTTATTATTAGTTGCTGCCAAAAATTTTTTCAAAATGGTAAATGGACCGCCAGTGGTAAAATTGACCGCAGAAACAACGCTTATGCTTTTTCCCattatt                                                                                                        <  1:913570/139‑1 (MQ=255)            tttAGCAGAATGGACTAATGCGATAAAACTGACATTTTCTTTATTATTAGTTGCTGCCAAAAATTTTTTCAAAATGGTAAATGGACCGCCAGTGGTAAAATTGACCGCAGAAACAACGATTATGCTTTTTCCCATAAtt                                                                                                        <  2:1195892/139‑1 (MQ=255)             ttAGCAGAATGGACTAATGCGATAAAACTGACATTTTCTTTATTATTAGTTGCTGCCAAAAATTTTTTCAAAATGGTAAATGGACCGCCAGTGGTAAAATTGACCGCAGAAACAACGATTATGCTTTTTcc                                                                                                               >  2:281616/1‑131 (MQ=255)             ttAGCAGAATGGACTAATGCGATAAAACTGACATTTTCTTTATTATTAGTTGCTGCCAAAAATTTTTTCAAAATGGTAAATGGACCGCCAGTGGTAAAATTGACCGCAGAAACAACGATTATGCTTTTTcc                                                                                                               <  1:281616/131‑1 (MQ=255)                      tGGACTAATGCGATAAAACTGACATTTTCTTTATTATTAGTTGCTGCCAAAAATTTTTTCAAAATGGTAAATGGACCGCCAGTGGTAAAATTGACCGCAGAAACAACGATTATGCTTTTTCCCATAATTTGGTCTCATg                                                                                              >  1:711896/1‑139 (MQ=255)                       ggACTAATGCGATAAAACTGACATTTTCTTTATTATTAGTTGCTGCCAAAAATTTTTTCAAAATGGTAAATGGACCGCCAGTGGTAAAATTGACCGCAGAAACAACGATTATGCTTTTTCCCATGATTTGGTCTCATGa                                                                                             <  1:950242/139‑1 (MQ=255)                           tAATGCGATAAAACTGACATTTTCTTTATTATTAGTTGCTGCCAAAAATTTTTTCAAAATGGTAAATGGACCGCCAGTGGTAAAATTGACCGCAGAAACAACGATTATGC‑TTTTCCCATAATTTGGTCTCATGATTGTa                                                                                        >  2:1067609/1‑139 (MQ=255)                                cGATAAAACTGACATTTTCTTTATTATTAGTTGCTGCCAAAAATTTTTTCAAAATGGTAAATGGACCGCCAGTGGTAAAATTGACCGCAGAAACAACGATTATGCTTTTTccc                                                                                                              <  2:351487/113‑1 (MQ=255)                                cGATAAAACTGACATTTTCTTTATTATTAGTTGCTGCCAAAAATTTTTTCAAAATGGTAAATGGACCGCCAGTGGTAAAATTGACCGCAGAAACAACGATTATGCTTTTTccc                                                                                                              >  1:351487/1‑113 (MQ=255)                                cGATAAAACTGACATTTTCTTTATTATTAGTTGCTGCCAAAAATTTTTTCAAAATGGTAAATGGACCGCCAGTGGTAAAATTGACCGCAGAAACAACGATTATGCTTTTTCCCATAATTTGGt                                                                                                    <  2:388525/123‑1 (MQ=255)                                cGATAAAACTGACATTTTCTTTATTATTAGTTGCTGCCAAAAATTTTTTCAAAATGGTAAATGGACCGCCAGTGGTAAAATTGACCGCAGAAACAACGATTATGCTTTTTCCCATAATTTGGt                                                                                                    >  1:388525/1‑123 (MQ=255)                                  aTAAAACTGACATTTTCTTTATTATTAGTTGCTGCCAAAAATTTTTTCAAAATGGTAAATGGACCGCCAGTGGTAAAATTGACCGCAGAAACAACGATTATGCTTTTTCCCATAATTTg                                                                                                      >  1:223851/1‑119 (MQ=255)                                  aTAAAACTGACATTTTCTTTATTATTAGTTGCTGCCAAAAATTTTTTCAAAATGGTAAATGGACCGCCAGTGGTAAAATTGACCGCAGAAACAACGATTATGCTTTTTCCCATAATTTg                                                                                                      <  2:223851/119‑1 (MQ=255)                                                     tattattaGTTGCTGCCAAAAATTTTTTCAAAATGGTAAATGGACCGCCAGTGGTAAAATTGACCGCAGAAACAACGATTATGCTTTTTCCCATAATTTGGTCTCATGATTGTATTTCTTTATGATTTTTGCTGGTACt                                                               <  1:311115/139‑1 (MQ=255)                                                     tattattaGTTGCTGCCAAAAATTTTTTCAAAATGGTAAATGGACCGCCAGTGGTAAAATTGACCGCAGAAACAACGATTATGCTTTTTCCCATAATTTGGTCTCATGATTGTATTTCTTTATGATTTTTGCTGGTACt                                                               <  2:711896/139‑1 (MQ=255)                                                            aGTTGCTGCCAAAAATTTTTTCAAAATGGTAAATGGACCGCCAGTGGTAAAATTGACCGCAGAAACAACGATTATGCTTTTTCCCATAATTTGGTCTCATGATTGTa                                                                                        >  2:688014/1‑107 (MQ=255)                                                            aGTTGCTGCCAAAAATTTTTTCAAAATGGTAAATGGACCGCCAGTGGTAAAATTGACCGCAGAAACAACGATTATGCTTTTTCCCATAATTTGGTCTCATGATTGTa                                                                                        <  1:688014/107‑1 (MQ=255)                                                            aGTTGCTGCCAAAAATTTTTTCAAAATGGTAAATGGACCGACAGTGGTAAAATTGACCGCAGAAACAACGATTATGCTTTTTCCCATAATTTGGTCTCATGATTGTATTTCTTTATGATTTTTGCTGGTACTCCCGCaa                                                        <  1:1092464/139‑1 (MQ=255)                                                             gTTGCTGCCAAAAATTTTTTCAAAATGGTAAATGGACCGCCAGTGGTAAAATTGACCGCAGAAACAACGATTATGCTTTTTCCCATAATTTGGTCTCATGATTGTATTTCTTTATGATTTTTGCTGGTa                                                                 <  2:915026/129‑1 (MQ=255)                                                             gTTGCTGCCAAAAATTTTTTCAAAATGGTAAATGGACCGCCAGTGGTAAAATTGACCGCAGAAACAACGATTATGCTTTTTCCCATAATTTGGTCTCATGATTGTATTTCTTTATGATTTTTGCTGGTa                                                                 >  1:915026/1‑129 (MQ=255)                                                                gctgcCAAAAATTTTTTCAAAATGGTAAATGGACCGCCAGTGGTAAAATTGACCGCAGAAACAACGATTATGCTTTTTCCCATAATTTGGTCTCAt                                                                                               <  1:1137567/96‑1 (MQ=255)                                                                gctgcCAAAAATTTTTTCAAAATGGTAAATGGACCGCCAGTGGTAAAATTGACCGCAGAAACAACGATTATGCTTTTTCCCATAATTTGGTCTCAt                                                                                               >  2:1137567/1‑96 (MQ=255)                                                                gctgcCAAAAATTTTTTCAAAATGGTAAATGGACCGCCAGTGGTAAAATTGACCGCAGAAACAACGATTATGC‑TTTTCCCATAATTTGGTCTCATGATTGTATTTCTTTATGATTTTTGCTGGTACTCCCGCAATGACa                                                   <  1:1067609/139‑1 (MQ=255)                                                                 ctgcCAAAAATTTTTTCAAAATGGTAAATGGACCGCCAGTGGTAAAATTGACCGCAGAAACAACGATTATGCTTTTTCCCATAATTTGGTCTCATGATTGTATTTCTTTATGATTTTTGCTGGTACTCCCGCAATGACa                                                   >  2:525001/1‑139 (MQ=255)                                                                  tgcCAAAAATTTTTTCAAAATGGTAAATGGACCGCCAGTGGTAAAATTGACCGCAGAAACAACGATTATGCTTTTTCCCATAATTTGGTCTCATGATTGTATTTCTTTATGATTTTTGCTGGTACTCCCGCAATGACa                                                   <  2:143157/138‑1 (MQ=255)                                                                  tgcCAAAAATTTTTTCAAAATGGTAAATGGACCGCCAGTGGTAAAATTGACCGCAGAAACAACGATTATGCTTTTTCCCATAATTTGGTCTCATGATTGTATTTCTTTATGATTTTTGCTGGTACTCCCGCAATGACa                                                   >  1:143157/1‑138 (MQ=255)                                                                  tgcCAAAAATTTTTTCAAAATGGTAAATGGACCGCCAGTGGTAAAATTGACCGCAGAAACAACGATTATGCTTTTTCCCATAATTTGGTCTCATGATTGTATTTCTTTATGATTTTTGCTGGTACTCCCGCAATGACAg                                                  >  1:906958/1‑139 (MQ=255)                                                                    ccAAAAATTTTTTCAAAATGGTAAATGGACCGCCAGTGGTAAAATTGACCGCAGAAACAACGATTATGCTTTTTCCCATAATTTGGTCTCATGATTGTATTTCTTTATGATTTTTGCTGGTACTCCCGCAATGACAGTa                                                >  1:500028/1‑139 (MQ=255)                                                                         aaTTTTTTCAAAATGGTAAATGGACCGCCAGTGGTAAAATTGACCGCAGAAACAACGATTATGCTTTTTCCCATAATTTGGTCTCATGATTGTATTTCTTTATGATTTTTGCTGGTACTCCCGCAATGACAGTATTTTc                                           <  2:500028/139‑1 (MQ=255)                                                                                                aCCGCCAGTGGTAAAATTGACCGCAGAAACAACGATTATGCTTTTTCCCATAATTTGGTCTCATGATTGTATTTCTTTATGATTTTTGCTGGTACTCCCGCAATGACAGTATTTTc                                           <  2:934826/116‑1 (MQ=255)                                                                                                aCCGCCAGTGGTAAAATTGACCGCAGAAACAACGATTATGCTTTTTCCCATAATTTGGTCTCATGATTGTATTTCTTTATGATTTTTGCTGGTACTCCCGCAATGACAGTATTTTc                                           >  1:934826/1‑116 (MQ=255)                                                                                                aCCGCCAGTGGTAAAATTGACCGCAGAAACAACGATTATGCTTTTTCCCATAATTTGGTCTCATGATTGTATTTCTTTATGATTTTTGCTGGTACTCCCGCAATGACAGTATTTTCGGGAATAGAACCTCTAACAACAg                    >  1:595150/1‑139 (MQ=255)                                                                                                 ccGCCAGTGGTAAAATTGACCGCAGAAACAACGATTATGCTTTTTCCCATAATTTGGTCTCATGATTGTATTTCTTTATGATTTTTGCTGGTa                                                                 >  1:874079/1‑93 (MQ=255)                                                                                                 ccGCCAGTGGTAAAATTGACCGCAGAAACAACGATTATGCTTTTTCCCATAATTTGGTCTCATGATTGTATTTCTTTATGATTTTTGCTGGTa                                                                 <  2:874079/93‑1 (MQ=255)                                                                                                  cGCCAGTGGTAAAATTGACCGCAGAAACAACGATTATGCTTTTTCCCATAAtt                                                                                                        <  1:355303/53‑1 (MQ=255)                                                                                                  cGCCAGTGGTAAAATTGACCGCAGAAACAACGATTATGCTTTTTCCCATAAtt                                                                                                        >  2:355303/1‑53 (MQ=255)                                                                                                  cGCCAGTGGTAAAATTGACCGCAGAAACAACGATTATGCTTTTTCCCATAATTTGGTCTCAt                                                                                               >  2:1157115/1‑62 (MQ=255)                                                                                                  cGCCAGTGGTAAAATTGACCGCAGAAACAACGATTATGCTTTTTCCCATAATTTGGTCTCAt                                                                                               <  1:1157115/62‑1 (MQ=255)                                                                                                  cGCCAGTGGTAAAATTGACCGCAGAAACAACGATTATGCTTTTTCCCATAATTTGGTCTCATGAt                                                                                            <  2:666749/65‑1 (MQ=255)                                                                                                  cGCCAGTGGTAAAATTGACCGCAGAAACAACGATTATGCTTTTTCCCATAATTTGGTCTCATGAt                                                                                            >  1:666749/1‑65 (MQ=255)                                                                                                  cGCCAGTGGTAAAATTGACCGCAGAAACAACGATTATGCTTTTTCCCATAATTTGGTCTCATGATTGTATTTCTTTATGAtt                                                                           >  1:1212675/1‑82 (MQ=255)                                                                                                  cGCCAGTGGTAAAATTGACCGCAGAAACAACGATTATGCTTTTTCCCATAATTTGGTCTCATGATTGTATTTCTTTATGAtt                                                                           <  2:1212675/82‑1 (MQ=255)                                                                                                  cGCCAGTGGTAAAATTGACCGCAGAAACAACGATTATGCTTTTTCCCATAATTTGGTCTCATGATTGTATTTCTTTATGATTTTTGCTGGTACTCCCGCa                                                         >  1:296981/1‑100 (MQ=255)                                                                                                  cGCCAGTGGTAAAATTGACCGCAGAAACAACGATTATGCTTTTTCCCATAATTTGGTCTCATGATTGTATTTCTTTATGATTTTTGCTGGTACTCCCGCa                                                         <  2:296981/100‑1 (MQ=255)                                                                                                  cGCCAGTGGTAAAATTGACCGCAGAAACAACGATTATGCTTTTTCCCATAATTTGGTCTCATGATTGTATTTCTTTATGATTTTTGCTGGTACTCCCGCAATGACAGTATTTTCGGGAATAGAACCTCTAACAACAGaa                  >  2:941621/1‑139 (MQ=255)                                                                                                       gTGGTAAAATTGACCGCAGAAACAACGATTATGCTTTTTCCCATAATTTGGTCTCATGATTGTATTTCTTTATGATTTTTGCTGGTACTCCCGCAATGACAGTATTTTCGGGAATAGAACCTCTAACAACAGAATTGgc             >  2:999721/1‑139 (MQ=255)                                                                                                          gTAAAATTGACCGCAGAAACAACGATTATGCTTTTTCCCATAATTTGGTCTCATGATTGTATTTCTTTATGATTTTTGCTGGTACTCCCGCAATGACAGTATTTTCGGGAATAGAACCTCTAACAACAGAATTGGCGcc          <  2:595150/139‑1 (MQ=255)                                                                                                              aaTTGACCGCAGAAACAACGATTATGCTTTTTCCCATAATTTGGTCTCATGATTGTATTTCTTTATGATTTTTGCTGGTACTCCCGCAATGACAGTATTTTCGGGAATAGAACCTCTAACAACAGAATTGGCGCCGACt      <  1:941621/139‑1 (MQ=255)                                                                                                               aTTGACCGCAGAAACAACGATTATGCTTTTTCCCATAATTTGGTCTCATGATTGTATTTCTTTATGATTTTTGCTGGTACTCCCGCAATGACAGTATTTTCGGGAATAGAACCTCTAACAACAGAATTGGCGCCGACTa     >  1:419375/1‑139 (MQ=255)                                                                                                                  gACCGCAGAAACAACGATTATGCTTTTTCCCATAATTTGGt                                                                                                    >  2:1146928/1‑41 (MQ=255)                                                                                                                  gACCGCAGAAACAACGATTATGCTTTTTCCCATAATTTGGt                                                                                                    <  1:1146928/41‑1 (MQ=255)                                                                                                                  gACCGCAGAAACAACGATTATGCTTTTTCCCATAATTTGGTCTCa                                                                                                <  1:775177/45‑1 (MQ=255)                                                                                                                  gACCGCAGAAACAACGATTATGCTTTTTCCCATAATTTGGTCTCa                                                                                                >  2:775177/1‑45 (MQ=255)                                                                                                                  gACCGCAGAAACAACGATTATGCTTTTTCCCATAATTTGGTCTCATGATTGTATTTCTTTATGATTTTTGCTGGTACTCCCGCAATGACAGTATTTTCGGGAATAGAACCTCTAACAACAGAATTGGCGCCGACTACGa  >  1:1260998/1‑139 (MQ=255)                                                                                                                   aCCGCAGAAACAACGATTATGCTTTTTCCCATAATTTGGTCTCATGa                                                                                             <  2:427744/47‑1 (MQ=255)                                                                                                                   aCCGCAGAAACAACGATTATGCTTTTTCCCATAATTTGGTCTCATGa                                                                                             >  1:427744/1‑47 (MQ=255)                                                                                                                   aCCGCAGAAACAACGATTATGCTTTTTCCCATAATTTGGTCTCATGATTGTATTTCTTTATGAt                                                                            <  2:472926/64‑1 (MQ=255)                                                                                                                   aCCGCAGAAACAACGATTATGCTTTTTCCCATAATTTGGTCTCATGATTGTATTTCTTTATGAt                                                                            >  1:472926/1‑64 (MQ=255)                                                                                                                     cGCAGAAACAACGATTATGCTTTTTCCCATAATTTGGTCTCATGATTGTATTTCTTTATGATTTTTGCTGGTACTCCCGCAATGACAGTATTTTCGGGAATAGaa                                 >  1:920502/1‑105 (MQ=255)                                                                                                                     cGCAGAAACAACGATTATGCTTTTTCCCATAATTTGGTCTCATGATTGTATTTCTTTATGATTTTTGCTGGTACTCCCGCAATGACAGTATTTTCGGGAATAGaa                                 <  2:920502/105‑1 (MQ=255)                                                                                                                       cAGAAACAACGATTATGCTTTTTCCCATAATTTGGTCTCATGATTGTATTTCTTTATGATTTTTGCTGGTACTCCCGCAATGACAGTATTTTCGGGAATAGAACCTCTaacaac                      >  2:1122769/1‑114 (MQ=255)                                                                                                                       cAGAAACAACGATTATGCTTTTTCCCATAATTTGGTCTCATGATTGTATTTCTTTATGATTTTTGCTGGTACTCCCGCAATGACAGTATTTTCGGGAATAGAACCTCTaacaac                      <  1:1122769/114‑1 (MQ=255)                                                                                                                                   |                                                                                                                          TCTTTTAACTCTTTAGCAGAATGGACTAATGCGATAAAACTGACATTTTCTTTATTATTAGTTGCTGCCAAAAATTTTTTCAAAATGGTAAATGGACCGCCAGTGGTAAAATTGACCGCAGAAACAACGACTATGCTTTTTCCCATAATTTGGTCTCATGATTGTATTTCTTTATGATTTTTGCTGGTACTCCCGCAATGACAGTATTTTCGGGAATAGAACCTCTAACAACAGAATTGGCGCCGACTACGA  >  NC_000913/2104364‑2104615

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