Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,979,952	A→G	*451* (TAA→TGA)	wzyE →	enterobacterial common antigen polymerase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,979,952	0	A	G	100.0%	146.4 / NA	42	*451* (TAA→TGA)	wzyE	enterobacterial common antigen polymerase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (15/27);  total (15/27)

GTTGTACTGGCTTTTTGAAAGCGCCGGACTCATTCATAAACGTACAAAATCATCGCTCCGGACGCAGGTTGAAGGATAACAATGAATAACAACACCACGGCACCAACCTATACGCTGCGTGGCTTACAGTTGATTGGTTGGCGTGATATGCAGC  >  NC_000913/3979875‑3980028                                                                              |                                                                             gttgtACTGGCTTTTTGAAAGCGCCGGACTCATTCATAAACGTACAAAATCATCGCTCCGGACGCAGGTTGAAGGATGAc                                                                            >  1:728929/1‑80 (MQ=255) gttgtACTGGCTTTTTGAAAGCGCCGGACTCATTCATAAACGTACAAAATCATCGCTCCGGACGCAGGTTGAAGGATGAc                                                                            >  1:728921/1‑80 (MQ=255)    gtACTGGCTTTTTGAAAGCGCCGGACTCATTCATAAACGTACAAAATCATCGCTCCGGACGCAGGTTGAAGGATGACAAt                                                                         >  1:728942/1‑80 (MQ=255)     tACTGGCTTTTTGAAAGCGCCGGACTCATTCATAAACGTACAAAATCATCGCTCCGGACGCAGGTTGAAGGATGACAATg                                                                        <  1:728874/80‑1 (MQ=255)       cTGGCTTTTTGAAAGCGCCGGACTCATTCATAAACGTACAAAATCATCGCTCCGGACGCAGGTTGAAGGATGACAATGaa                                                                      >  1:728940/1‑80 (MQ=255)                tGAAAGCGCCGGACTCATTCATAAACGTACAAAATCATCGCTCCGGACGCAGGTTGAAGGATGACAATGAATAACAacac                                                             <  2:728874/80‑1 (MQ=255)                  aaaGCGCCGGACTCATTCATAAACGTACAAAATCATCGCTCCGGACGCAGGTTGAAGGATGACAATGAATAACAAcaccac                                                          <  2:728875/81‑1 (MQ=255)                    aGCGCCGGACTCATTCATAAACGTACAAAATCATCGCTCCGGACGCAGGTTGAAGGATGACAATGAATAACAACACCACgg                                                        <  2:728876/81‑1 (MQ=255)                     gcgcCGGACTCATTCATAAACGTACAAAATCATCGCTCCGGACGCAGGTTGAAGGATGACAATGAATAACAACACCACGGc                                                       <  2:728877/81‑1 (MQ=255)                       gcCGGACTCATTCATAAACGTACAAAATCATCGCTCCGGACGCAGGTTGAAGGATGACAATGAATAACAACACCACGGCAc                                                     <  1:728879/81‑1 (MQ=255)                       gcCGGACTCATTCATAAACGTACAAAATCATCGCTCCGGACGCAGGTTGAAGGATGACAATGAATAACAACACCACGGCAc                                                     >  2:728936/1‑81 (MQ=255)                       gcCGGACTCATTCATAAACGTACAAAATCATCGCTCCGGACGCAGGTTGAAGGATGACAATGAATAACAACACCACGGCAc                                                     >  1:728945/1‑81 (MQ=255)                          ggACTCATTCATAAACGTACAAAATCATCGCTCCGGACGCAGGTTGAAGGATGACAATGAATAACAACACCACGGCACCaa                                                  <  1:728880/81‑1 (MQ=255)                           gACTCATTCATAAACGTACAAAATCATCGCTCCGGACGCAGGTTGAAGGATGACAATGAATAACAACACCACGGCACCAAc                                                 <  2:728880/81‑1 (MQ=255)                                attcatAAACGTACAAAATCATCGCTCCGGACGCAGGTTGAAGGATGACAATGAATAACAACACCACGGCACCAACCtata                                            <  1:728882/81‑1 (MQ=255)                                attcatAAACGTACAAAATCATCGCTCCGGACGCAGGTTGAAGGATGACAATGAATAACAACACCACGGCACCAACCtat                                             <  1:728883/80‑1 (MQ=255)                                    atAAACGTACAAAATCATCGCTCCGGACGCAGGTTGAAGGATGACAATGAATAACAACACCACGGCACCAACCTATACGCt                                        <  2:728883/81‑1 (MQ=255)                                      aaaCGTACAAAATCATCGCTCCGGACGCAGGTTGAAGGATGACAATGAATAACAACACCACGGCACCAACCTATACGCTGc                                      >  1:728950/1‑81 (MQ=255)                                         cGTACAAAATCATCGCTCCGGACGCAGGTTGAAGGATGACAATGAATAACAACACCACGGCACCAACCTATACGCTGCGTg                                   <  2:728884/81‑1 (MQ=255)                                             cAAAATCATCGCTCCGGACGCAGGTTGAAGGATGACAATGAATAACAACACCACGGCACCAACCTATACGCTGCGTGGCt                                >  1:728938/1‑80 (MQ=255)                                              aaaaTCATCGCTCCGGACGCAGGTTGAAGGATGACAATGAATAACAACACCACGGCACCAACCTATACGCTGCGTGGCTTa                              <  2:728887/81‑1 (MQ=255)                                              aaaaTCATCGCTCCGGACGCAGGTTGAAGGATGACAATGAATAACAACACCACGGCACCAACCTATACGCTGCGTGGCTTa                              <  1:728887/81‑1 (MQ=255)                                               aaaTCATCGCTCCGGACGCAGGTTGAAGGATGACAATGAATAACAACACCACGGCACCAACCTATACGCTGCGTGGCTTAc                             <  1:728889/81‑1 (MQ=255)                                                aaTCATCGCTCCGGACGCAGGTTGAAGGATGACAATGAATAACAACACCACGGCACCAACCTATACGCTGCGTGGCTTACa                            <  1:728890/81‑1 (MQ=255)                                                   catcGCTCCGGACGCAGGTTGAAGGATGACAATGAATAACAACACCACGGCACCAACCTATACGCTGCGTGGCTTACAGtt                         >  2:728942/1‑81 (MQ=255)                                                    atcGCTCCGGACGCAGGTTGAAGGATGACAATGAATAACAACACCACGGCACCAACCTATACGCTGCGTGGCTTACAGtt                         >  2:728935/1‑80 (MQ=255)                                                    atcGCTCCGGACGCAGGTTGAAGGATGACAATGAATAACAACACCACGGCACCAACCTATACGCTGCGTGGCTTACAGtt                         <  2:728890/80‑1 (MQ=255)                                                       gCTCCGGACGCAGGTTGAAGGATGACAATGAATAACAACACCACGGCACCAACCTATACGCTGCGTGGCTTACAGTTGAt                      <  2:728891/80‑1 (MQ=255)                                                           cGGACGCAGGTTGAAGGATGACAATGAATAACAACACCACGGCACCAACCTATACGCTGCGTGGCTTACAGTTGAttggtt                 >  2:728950/1‑81 (MQ=255)                                                           cGGACGCAGGTTGAAGGATGACAATGAATAACAACACCACGGCACCAACCTATACGCTGCGTGGCTTACAGTTGAttgg                   <  1:728893/79‑1 (MQ=255)                                                             gACGCAGGTTGAAGGATGACAATGAATAACAACACCACGGCACCAACCTATACGCTGCGTGGCTTACAGTTGAttggttgg               <  1:728894/81‑1 (MQ=255)                                                              aCGCAGGTTGAAGGATGACAATGAATAACAACACCACGGCACCAACCTATACGCTGCGTGGCTTACAGTTGATTGGTTGGc              <  2:728894/81‑1 (MQ=255)                                                               cGCAGGTTGAAGGATGACAATGAATAACAACACCACGGCACCAACCTATACGCTGCGTGGCTTACAGTTGATTGGTTGGCg             <  2:728895/81‑1 (MQ=255)                                                                 cAGGTTGAAGGATGACAATGAATAACAACACCACGGCACCAACCTATACGCTGCGTGGCTTACAGTTGATTGGTTGGCGTg           <  1:728897/81‑1 (MQ=255)                                                                  aGGTTGAAGGATGACAATGAATAACAACACCACGGCACCAACCTATACGCTGCGTGGCTTACAGTTGATTGGTTGGCGTGa          >  2:728943/1‑81 (MQ=255)                                                                    gTTGAAGGATGACAATGAATAACAACACCACGGCACCAACCTATACGCTGCGTGGCTTACAGTTGATTGGTTGGCGTGata        <  2:728897/81‑1 (MQ=255)                                                                     ttGAAGGATGACAATGAATAACAACACCACGGCACCAACCTATACGCTGCGTGGCTTACAGTTGATTGGTTGGCGTGatat       <  2:728898/81‑1 (MQ=255)                                                                     ttGAAGGATGACAATGAATAACAACACCACGGCACCAACCTATACGCTGCGTGGCTTACAGTTGATTGGTTGGCGTGatat       >  2:728940/1‑81 (MQ=255)                                                                     ttGAAGGATGACAATGAATAACAACACCACGGCACCAACCTATACGCTGCGTGGCTTACAGTTGATTGGTTGGCGTGatat       >  1:728939/1‑81 (MQ=255)                                                                          ggATGACAATGAATAACAACACCACGGCACCGACCTATACGCTGCGTGGCTTACAGTTGATTGGTTGGCGTGATATgcagc  <  2:728900/81‑1 (MQ=255)                                                                          ggATGACAATGAATAACAACACCACGGCACCAACCTATACGCTGCGTGGCTTACAGTTGATTGGTTGGCGTGATATgcagc  <  1:728900/81‑1 (MQ=255)                                                                          ggATGACAATGAATAACAACACCACGGCACCAACCTATACGCTGCGTGGCTTACAGTTGATTGGTTGGCGTGATATgcag   >  2:728962/1‑80 (MQ=255)                                                                              |                                                                             GTTGTACTGGCTTTTTGAAAGCGCCGGACTCATTCATAAACGTACAAAATCATCGCTCCGGACGCAGGTTGAAGGATAACAATGAATAACAACACCACGGCACCAACCTATACGCTGCGTGGCTTACAGTTGATTGGTTGGCGTGATATGCAGC  >  NC_000913/3979875‑3980028

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