Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	REL606	3,940,843	T→A	I148N (ATC→AAC)	wecF →	putative enterobacterial common antigen polymerase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	REL606	3,940,843	0	T	A	100.0%	167.0 / NA	46	I148N (ATC→AAC)	wecF	putative enterobacterial common antigen polymerase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base A (22/24);  total (22/24)

CGTTGGTAAGCGTCGGCATCTTCTTCATGCACAACGGTTTTTTGCTGTTCCGTCTTAACTCCTACAGTCAGATCTTCTCCAGTGAAGTCTCCGGCGTGGCGTTAAAACGTTTCTTTTACTTTTTCATCCCGGCGATGCTGGTGGT  >  REL606/3940771‑3940915                                                                         |                                                                         cGTTGGTAAGCGTCGGCATCTTCTTCATGCACAACGGTTTTTTGCTGTTCCGTCTTAACTCCTACAGTCAGAactt                                                                       >  1:789348/1‑76 (MQ=255) cGTTGGTAAGCGTCGGCATCTTCTTCATGCACAACGGTTTTTTGCTGTTCCGTCTTAACTCCTACAGTCAGAactt                                                                       <  1:596833/76‑1 (MQ=255) cGTTGGTAAGCGTCGGCATCTTCTTCATGCACAACGGTTTATTGCTGTTCCGTCTTAACTCCTACAGTCAGAactt                                                                       >  1:1140457/1‑76 (MQ=255)    tGGTAAGCGTCGGCATCTTCTTCATGCACAACGGTTTTTTGCTGTTCCGTCTTAACTCCTACAGTCAGAACTtctc                                                                    <  1:603377/76‑1 (MQ=255)     ggTAAGCGTCGGCATCTTCTTCATGCACAACGGTTTTTTGCTGTTCCGTCTTAACTCCTACAGTCAGAACTTCTcc                                                                   >  2:647185/1‑76 (MQ=255)       tAAGCGTCGGCATCTTCTTCATGCACAACGGTTTTTTGCTGTTCCGTCTTAACTCCTACAGTCAGAACTTCTCCAg                                                                 <  1:1183138/76‑1 (MQ=255)        aaGCGTCGGCATCTTCTTCATGCACAACGGTTTTTTGCTGTTCCGTCTTAACTCCTACAGTCAGAACTTCTCCAGt                                                                >  2:327712/1‑76 (MQ=255)         aGCGTCGGCATCTTCTTCATGCACAACGGTTTTTTGCTGTTCCGTCTTAACTCCTACAGTCAGAACTTCTCCAGTg                                                               <  1:664391/76‑1 (MQ=255)          gCGTCGGCATCTTCTTCATGCACAACGGTTTTTTGCTGTTCCGTCTTAACTCCTACAGTCAGAACTTCTCCAGTGa                                                              <  1:1328800/76‑1 (MQ=255)           cGTCGGCATCTTCTTCATGCACAACGGTTTTTTGCTGTTCCGTCTTAACTCCTACAGTCAGAACTTCTCCAGTGaa                                                             <  2:1653666/76‑1 (MQ=255)            gTCGGCATCTTCTTCATGCACAACGGTTTTTTGCTGTTCCGTCTTAACTCCTACAGTCAGAACTTCTCCAGTGAAg                                                            <  1:687857/76‑1 (MQ=255)             tCGGCATCTTCTTCATGCACAACGGTTTTTTGCTGTTCCGTCTTAACTCCTACAGTCAGAACTTCTCCAGTGAAGt                                                           >  2:476617/1‑76 (MQ=255)              cGGCATCTTCTTCATGCACAACGGTTTTTTGCTGTTCCGTCTTAACTCCTACAGTCAGAACTTCTCCAGTGAAGtc                                                          >  1:1008937/1‑76 (MQ=255)                gCATCTTCGTCATGCACAACGGTTTTTTGCTGTTCCGTCTTAACTCCTACAGTCAGAACTTCTCCAGTGAAGtctc                                                        <  2:562675/76‑1 (MQ=255)                  aTCTTCTTCATGCACAACGGTTTTTTGCTGTTCCGTCTTAACTCCTACAGTCAGAACTTCTCCAGTGAAGTCTCCg                                                      >  1:1295718/1‑76 (MQ=255)                  aTCTTCTTCATGCACAACGGTTTTTTGCTGTTCCGTCTTAACTCCTACAGTCAGAACTTCTCCAGTGAAGTCTCCg                                                      <  1:1138407/76‑1 (MQ=255)                    cttcttCATGCACAACGGTTTTTTGCTGTTCCGTCTTAACTCCTACAGTCAGAACTTCTCCAGTGAAGTCTCCggc                                                    <  2:928327/76‑1 (MQ=255)                         tcATGCACAACGGTTTTTTGCTGTTCCGTCTTAACTCCTACAGTCAGAACTTCTCCAGTGAAGTCTCCggcgtggc                                               >  2:142311/1‑76 (MQ=255)                            tGCACAACGGTTTTTTGCTGTTCCGTCTTAACTCCTACAGTCAGAACTTCTCCAGTGAAGTCTCCGGCGTGGCGtt                                            <  2:1008936/76‑1 (MQ=255)                               acaACGGTTTTTTGCTGTTCCGTCTTAACTCCTACAGTCAGAACTTCTCCAGTGAAGTCTCCGGCGTGGCGTTaaa                                         >  1:447117/1‑76 (MQ=255)                               acaACGGTTTTTTGCTGTTCCGTCTTAACTCCTACAGTCAGAACTTCTCCAGTGAAGTCTCCGGCGTGGCGTTaaa                                         <  1:1432608/76‑1 (MQ=255)                               acaACGGTTTTTTGCTGTTCCGTCTTAACTCCTACAGTCAGAACTTCTCCAGTGAAGTCTCCGGCGTGGCGTTaaa                                         <  1:1365997/76‑1 (MQ=255)                               acaACGGTTTTTTGCTGTTCCGTCTTAACTCCTACAGTCAGAACTTCTCCAGTGAAGTCTCCGGCGTGGCGTTaaa                                         <  1:1286833/76‑1 (MQ=255)                                 aaCGGTTTTTTGCTGTTCCGTCTTAACTCCTACAGTCAGAACTTCTCCAGTGAAGTCTCCGGCGTGGCGTTAAAAc                                       >  1:973013/1‑76 (MQ=255)                                  aCGGTTTTTTGCTGTTCCGTCTTAACTCCTACAGTCAGAACTTCTCCAGTGAAGTCTCCGGCGTGGCGTTAAAACg                                      >  1:204108/1‑76 (MQ=255)                                  aCGGTTTTTTGCTGTTCCGTCTTAACTCCTACAGTCAGAACTTCTCCAGTGAAGTCTCCGGCGTGGCGTTAAAACg                                      <  1:360301/76‑1 (MQ=255)                                        ttttGCTGTTCCGTCTTAACTCCTACAGTCAGAACTTCTCCAGTGAAGTCTCCGGCGTGGCGTTAAAACGTTTCtt                                <  2:1176647/76‑1 (MQ=255)                                         tttGCTGTTCCGTCTTAACTCCTACAGTCAGAACTTCTCCAGTGAAGTCTCCGGCGTGGCGTTAAAACGTTTCttt                               >  2:943467/1‑76 (MQ=255)                                          ttGCTGTTCCGTCTTAACTCCTACAGTCAGAACTTCTCCAGTGAAGTCTCCGGCGTGGCGTTAAAACGTTTCtttt                              <  1:1009914/76‑1 (MQ=255)                                          ttGCTGTTCCGTCTTAACTCCTACAGTCAGAACTTCTCCAGTGAAGTCTCCGGCGTGGCGTTAAAACGTTTCtttt                              <  1:694276/76‑1 (MQ=255)                                              tGTTCCGTCTTAACTCCTACAGTCAGAACTTCTCCAGTGAAGTCTCCGGCGTGGCGTTAAAACGTTTCTTTTACtt                          >  1:1234764/1‑76 (MQ=255)                                               gTTCCGTCTTAACTCCTACAGTCAGAACTTCTCCAGTGAAGTCTCCGGCGTGGCGTTAAAACGTTTCTTTTACttt                         >  2:527544/1‑76 (MQ=255)                                               gTTCCGTCTTAACTCCTACAGTCAGAACTTCTCCAGTGAAGTCTCCGGCGTGGCGTTAAAACGTTTCTTTTACctt                         >  2:1141325/1‑76 (MQ=255)                                                     tCTTAACTCCTACAGTCAGAACTTCTCCAGTGAAGTCTCCGGCGTGGCGTTAAAACGTTTCTTTTACTTTTTCATc                   >  1:378740/1‑76 (MQ=255)                                                      cTTAACTCCTACAGTCAGAACTTCTCCAGTGAAGTCTCCGGCGTGGCGTTAAAACGTTTCTTTTACTTTTTCATcc                  >  2:619998/1‑76 (MQ=255)                                                      cTTAACTCCTACAGTCAGAACTTCTCCAGTGAAGTCTCCGGCGTGGCGTTAAAACGTTTCTTTTACTTTTTCATcc                  >  1:728493/1‑76 (MQ=255)                                                       ttAACTCCTACAGTCAGAACTTCTCCAGTGAAGTCTCCGGCGTGGCGTTAAAACGTTTCTTTTATTTTTTCATccc                 >  1:1052393/1‑76 (MQ=255)                                                       ttAACTCCTACAGTCAGAACTTCTCCAGTGAAGTCTCCGGCGTGGCGTTAAAACGTTTCTTTTACTTTTTCATccc                 <  1:1313632/76‑1 (MQ=255)                                                        tAACTCCTACAGTCAGAACTTCTCCAGTGAAGTCTCCGGCGTGGCGTTAAAACGTTTCTTTTACTTTTTCATCCCg                >  2:1663373/1‑76 (MQ=255)                                                        tAACTCCTACAGTCAGAACTTCTCCAGTGAAGTCTCCGGCGTGGCGTTAAAACGTTTCTTTTACTTTTTCATCCCg                >  2:134031/1‑76 (MQ=255)                                                         aaCTCCTACAGTCAGAACTTCTCCAGTGAAGTCTCCGGCGTGGCGTTAAAACGTTTCTTTTACTTTTTCATCCCgg               <  1:966043/76‑1 (MQ=255)                                                            tCCTACAGTCAGAACTTCTCCAGTGAAGTCTCCGGCGTNGCGTTAAAACGTTTCTTTTACTTTTTCATCCCGGCGa            <  1:1776691/76‑1 (MQ=255)                                                             ccTACAGTCAGAACTTCTCCAGTGAAGTCTCCGGCGTGGCGTTAAAACGTTTCTTTTACTTTTTCATCCCGGCGAt           >  1:1395348/1‑76 (MQ=255)                                                                aCAGTCAGAACTTCTCCAGTGAAGTCTCCGGCGTGGCGTTAAAACGTTTCTTTTACTTTTTCATCCCGGCGATGCt        <  2:1309352/76‑1 (MQ=255)                                                                aCAGTCAGAACTTCTCCAGTGAAGTCTCCGGCGTGGCGTTAAAACGTTTCTTTTACTTTTTCATCCCGGCGATGCt        <  1:998306/76‑1 (MQ=255)                                                                  aGTCAGAACTTCTCCAGTGAAGTCTCCGGCGTGGCGTTAAAACGTTTCTTTTACTTTTTCATCCCGGCGATGCtgg      >  1:212242/1‑76 (MQ=255)                                                                  aGTCAGAACTTCTCCAGTGAAGTCTCCGGCGTGGCGTTAAAACGTTTCTTTTACTTTTTCATCCCGGCGATGCtgg      >  1:299030/1‑76 (MQ=255)                                                                    tCAGAACTTCTCCAGTGAAGTCTCCGGCGTGGCGTTAAAACGTTTCTTTTACTTTTTCATCCCGGCGATGCtggtg    >  2:696504/1‑76 (MQ=255)                                                                      aGAACTTCTCCAGTGAAGTCTCCGGCGTGGCGTTAAAACGTTTCTTTTACTTTTTCATCCCGGCGATGCTggtggt  <  2:555396/76‑1 (MQ=255)                                                                      aGAACTTCTCCAGTGAAGTCTCCGGCGTGGCGTTAAAACGTTTCTTTTACTTTTTCATCCCGGCGATGCTggtggt  <  1:476618/76‑1 (MQ=255)                                                                         |                                                                         CGTTGGTAAGCGTCGGCATCTTCTTCATGCACAACGGTTTTTTGCTGTTCCGTCTTAACTCCTACAGTCAGATCTTCTCCAGTGAAGTCTCCGGCGTGGCGTTAAAACGTTTCTTTTACTTTTTCATCCCGGCGATGCTGGTGGT  >  REL606/3940771‑3940915

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