Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NZ_CP009273	2,122,267	A→C	S35A (TCG→GCG)	wcaF ←	colanic acid biosynthesis acetyltransferase WcaF

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NZ_CP009273	2,122,267	0	A	C	100.0%	130.3 / NA	38	S35A (TCG→GCG)	wcaF	colanic acid biosynthesis acetyltransferase WcaF
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base C (17/21);  total (17/21)

CTGACGGACGAATAACTACGTTTTTTCCTATTTTTGCTCCGAATAAACGTAATAAAAAAGCCCGCCAGCGATACAATACTTGTGGCGACCAGGCAAATATTGTTGCCTGTACTGCCCACCATAATTGCACTTTAATAGCGTTGCCGCCCCGGAACCCTTTCGGCACCGAGAAA  >  NZ_CP009273/2122180‑2122352                                                                                        |                                                                                      cTGGCGGACGAATAACTACGTTTTTTCCTATTATTGCCCCGAAAAACCTTAATAAAAAAGCCCGCCAGCGATACAATACTTGTGGCGccc                                                                                     <  1:1009762/90‑1 (MQ=255) cTGACGGACGAATAACTACGTTTTTTCCTATTTTTGCTCCGAATAAACGTAATAAAAAAGCCCGCCAGCGATACAATACTTGTGGCGccc                                                                                     <  1:152508/90‑1 (MQ=255) cTGACGGACGAATAACTACGTTTTTTCCTATTTTTGCTCCGAATAAACGTAATAAAAAAGCCCGCCAGCGATACAATACTTGTGGCGccc                                                                                     <  1:527094/90‑1 (MQ=255)   gacggacgAATAACTACGTTTTTTCCTATTTTTGCTCCGAATAAACGTAATAAAAAAGCCCGCCAGCGATACAATACTTGTGGCGCCCAg                                                                                   <  2:583901/90‑1 (MQ=255)   gacggacgAATAACTACGTTTTTTCCTATTTTTGCTCCGAATAAACGTAATAAAAAAGCCCGCCAGCGATACAATACTTGTGGCGCCCAg                                                                                   >  1:510835/1‑90 (MQ=255)     cggacgAATAACTACGTTTTTTCCTATTTTTGCTCCGAATAAACGTAATAAAAAAGCCCGCCAGCGATACAATACTTGTGGCGCCCAGGc                                                                                 >  1:603887/1‑90 (MQ=255)      ggacgAATAACTACGTTTTTTCCTATTTTTGCTCCGAATAAACGTAATAAAAAAGCCCGCCAGCGATACAATACTTGTGGCGCCCAGGCa                                                                                <  2:1089259/90‑1 (MQ=255)             tAACTACGTTTTTTCCTATTTTTGCTCCGAATAAACGTAATAAAAAAGCCCGCCAGCGATACAATACTTGTGGCGCCCAGGCAAATAttg                                                                         <  2:603887/90‑1 (MQ=255)             tAACTACGTTTTTTCCTATTTTTGCTCCGAATAAACGTAATAAAAAAGCCCGCCAGCGATACAATACTTGTGGCGCCCAGGCAAATAttg                                                                         <  2:297924/90‑1 (MQ=255)             tAACTACGTTTTTTCCTATTTTTGCTCCGAATAAACGTAATAAAAAAGCCCGCCAGCGATACAATACTTGTGGCGCCCAGGCAAATAttg                                                                         <  2:775059/90‑1 (MQ=255)             tAACTACGTTTTTTCCTATTTTTGCTCCGAATAAACGTAATAAAAAAGCCCGCCAGCGATACAATACTTGTGGCGCCCAGGCAAATAttg                                                                         <  1:366220/90‑1 (MQ=255)             tAACTACGTTTTTTCCTATTTTTGCTCCGAATAAACGTAATAAAAAAGCCCGCCAGCGATACAATACTTGTGGCGCCCAGGCAAATAttg                                                                         <  2:866974/90‑1 (MQ=255)                       ttttCCTATTTTTGCTCCGAATAAACGTAATAAAAAAGCCCGCCAGCGATACAATACTTGTGGCGCCCAGGCAAATAttg                                                                         >  1:685286/1‑80 (MQ=255)                       ttttCCTATTTTTGCTCCGAATAAACGTAATAAAAAAGCCCGCCAGCGATACAATACTTGTGGCGCCCAGGCAAATAttg                                                                         <  2:685286/80‑1 (MQ=255)                       ttttCCTATTTTTGCTCCGAATAAACGTAATAAAAAAGCCCGCCAGCGATACAATACTTGTGGCGCCCAGGCAAATATTGTTGCCTGTAc                                                               >  2:448398/1‑90 (MQ=255)                          tCCTATTTTTGCTCCGAAAAAACGTAAAAAAAAAGCCCGCCAGCGATACAATACTTGTGGCGCCCAGGCAAATATTGTTGCCTGTACTGc                                                            <  1:708802/90‑1 (MQ=255)                              aTTTTTGCTCCGAATAAACGTAATAAAAAAGCCCGCCAGCGATACAATACTTGTGGCGCCCAGGCAAATATTGTTGcc                                                                    >  1:727963/1‑78 (MQ=255)                              aTTTTTGCTCCGAATAAACGTAATAAAAAAGCCCGCCAGCGATACAATACTTGTGGCGCCCAGGCAAATATTGTTGcc                                                                    <  2:727963/78‑1 (MQ=255)                               tttttGCTCCGAATAAACGTAATAAAAAAGCCCGCCAGCGATACAATACTTGTGGCGCCCAGGCAAATATTGTTGCCTGTACTGCccacc                                                       <  2:1134321/90‑1 (MQ=255)                                       ccGAATAAACGTAATAAAAAAGCCCGCCAGCGATACAATACTTGTGGCGCCCAGGCAAATATTGTTGCCTGTACTGCCCACCATAATTGc                                               <  2:360174/90‑1 (MQ=255)                                       ccGAATAAACGTAATAAAAAAGCCCGCCAGCGATACAATACTTGTGGCGCCCAGGCAAATATTGTTGCCTGTACTGCCCACCATAATTGc                                               <  2:1101861/90‑1 (MQ=255)                                       ccGAATAAACGTAATAAAAAAGCCCGCCAGCGATACAATACTTGTGGCGCCCAGGCAAATATTGTTGCCTGTACTGCCCACCATAATTGc                                               <  1:771014/90‑1 (MQ=255)                                         gAATAAACGTAATAAAAAAGCCCGCCAGCGATACAATACTTGTGGCGCCCAGGCAAATATTGTTGCCTGTACTGCCCACCATAATTGCAc                                             <  1:780897/90‑1 (MQ=255)                                         gAATAAACGTAATAAAAAAGCCCGCCAGCGATACAATACTTGTGGCGCCCAGGCAAATATTGTTGCCTGTACTGCCCACCATAATTGCAc                                             >  2:1104998/1‑90 (MQ=255)                                                   aataaAAAAGCCCGCCAGCGATACAATACTTGTGGCGCCCAGGCAAATATTGTTGCCTGTACTGCCCACCATAATTGCACTTTAATAGCg                                   >  1:622350/1‑90 (MQ=255)                                                   aataaAAAAGCCCGCCAGCGATACAATACTTGTGGCGCCCAGGCAAATATTGTTGCCTGTACTGCCCACCATAATTGCACTTTAATAGCg                                   >  2:258654/1‑90 (MQ=255)                                                   aataaAAAAGCCCGCCAGCGATACAATACTTGTGGCGCCCAGGCAAATATTGTTGCCTGTACTGCCCACCATAATTGCACTTTAATAGCg                                   >  1:538200/1‑90 (MQ=255)                                                   aataaAAAAGCCCGCCAGCGATACAATACTTGTGGCGCCCAGGCAAATATTGTTGCCTGTACTGCCCACCATAATTGCACTTTAATAGCg                                   >  1:1213435/1‑90 (MQ=255)                                                          aaGCCCGCCAGCGATACAATACTTGTGGCGCCCAGGCAAATATTGTTGCCTGTACTGCCCACCATAATTGCACTTTAATAGCGTTgccgc                            <  1:448398/90‑1 (MQ=255)                                                                   aGCGATACAATACTTGTGGCGCCCAGGCAAATATTGTTGCCTGTACTGCCCACCATAATTGCACTTTAATAGCGTTGCCGCCCCGGAAcc                   <  1:450687/90‑1 (MQ=255)                                                                        tACAATACTTGTGGCGCCCAGGCAAATATTGTTGCCTGTACTGCCCACCATAATTGCACTTTAATAGCGTTGCCGCCCCGGAACCCTTTc              <  2:1218520/90‑1 (MQ=255)                                                                           aaTACTTGTGGCGCCCAGGCAAATATTGTTGCCTGTACTGCCCACCATAATTGCACTTTAATAGCGTTGCCGCCCCGGAACCCTTTCGGc           <  2:1187551/90‑1 (MQ=255)                                                                               cTTGTGGCGCCCAGGCAAATATTGTTGCCTGTACTGCCCACCATAATTGCACTTTAATAGCGTTGCCGCCCCGGAACCCTTTCGGCACCg       >  2:285316/1‑90 (MQ=255)                                                                                 tgtgGCGCCCAGGCAAATATTGTTGCCTGTACTGCCCACCATAATTGCAc                                             <  2:1100483/50‑1 (MQ=255)                                                                                 tgtgGCGCCCAGGCAAATATTGTTGCCTGTACTGCCCACCATAATTGCAc                                             >  1:1100483/1‑50 (MQ=255)                                                                                 tgtgGCGCCCAGGCAAATATTGTTGCCTGTACTGCCCACCATAATTGCACTTTAATAGCGTTGCCGCCCCGGAACCCTTTCGGCACCgag     >  2:916917/1‑90 (MQ=255)                                                                                 tgtgGCGCCCAGGCAAATATTGTTGCCTGTACTGCCCACCATAATTGCACTTTAATAGCGTTGCCGCCCCGGAACCCTTTCGGCACCgag     >  2:1073722/1‑90 (MQ=255)                                                                                 tgtgGCGCCCAGGCAAATATTGTTGCCTGTACTGCCCACCATAATTGCACTTTAATAGCGTTGCCGCCCCGGAACCCTTTCGGCACCgag     >  1:451338/1‑90 (MQ=255)                                                                                 tgtgGCGCCCAGGCAAATATTGTTGCCTGTACTGCCCACCATAATTGCACTTTAATAGCGTTGCCGCCCCGGAACCCTTTCGGCACCgag     >  1:18854/1‑90 (MQ=255)                                                                                 tgtgGCGCCCAGGCAAATATTGTTGCCTGTACTGCCCACCATAATTGCACTTTAATAGCGTTGCCGCCCCGGAACCCTTTCGGCACCgag     >  1:159132/1‑90 (MQ=255)                                                                                    ggCGCCCAGGCAAATATTGTTGCCTGTACTGCCCACCATAATTGCACTTTAATAGCGTTGCCGCCCCGGAACCCTTTCGGCACCGAGaaa  <  1:916917/90‑1 (MQ=255)                                                                                        |                                                                                      CTGACGGACGAATAACTACGTTTTTTCCTATTTTTGCTCCGAATAAACGTAATAAAAAAGCCCGCCAGCGATACAATACTTGTGGCGACCAGGCAAATATTGTTGCCTGTACTGCCCACCATAATTGCACTTTAATAGCGTTGCCGCCCCGGAACCCTTTCGGCACCGAGAAA  >  NZ_CP009273/2122180‑2122352

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