Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,131,907	A→G	V55A (GTC→GCC)	wcaC ←	putative colanic acid biosynthesis glycosyl transferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,131,907	0	A	G	100.0%	126.6 / NA	35	V55A (GTC→GCC)	wcaC	putative colanic acid biosynthesis glycosyl transferase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (18/17);  total (18/17)

AAAATGCAGGACCACCGGACCCGCTGTGCGAGTAATGGTGCGATATAACTCATTGAAATTGCCAAACAGATCGCGATTAAACAGACGAAACAGAGCAATATTCGCCATCGCGGTCATCCGCGGCGTATGTTTGATGACCTGCGGATAGTTCTGATGGCTGACGCTCTCTTTGCCGCCTTTGCCGTAACCGTAGACAAAATGTGACGCCAGCCCCTGTTGCAGCGCACGCTGGTGGAGATCTAACGCCACACCTGCTGCCCCGCCTTC  >  NC_000913/2131771‑2132037                                                                                                                                         |                                                                                                                                   aaaaTGCAGGACCACCGGACCCGCTGTGCGAGTAATGGTGCGATATAACTCATTGAAATTGCCAAACAGATCGCGATTAAACAGACGAAACAGAGCAATATTCGCCATCGCGGTCATCCGCGGCGTATGTTTGATGGCCt                                                                                                                                 <  2:302309/140‑1 (MQ=255)          gACCACCGGACCCGCTGTGCGAGTAATGGTGCGATATAACTCATTGAAATTGCCAAACAGATCGCGATTAAACAGACGAAACAGAGCAATATTCGCCATCGCGGTCATCCGCGGCGTATGTTTGATGGCCTGCGGATAGt                                                                                                                        >  2:139193/1‑140 (MQ=255)               ccGGACCCGCTGTGCGAGTAATGGTGCGATATAACTCATTGAAATTGCCAAACAGATCGCGATTAAACAGACGAAACAGAGCAATATTCGCCATCGCGGTCATCCGCGGCGTATGTTTGATGGCCTGCGGATAGTTCTg                                                                                                                    <  1:222926/139‑1 (MQ=255)               ccGGACCCGCTGTGCGAGTAATGGTGCGATATAACTCATTGAAATTGCCAAACAGATCGCGATTAAACAGACGAAACAGAGCAATATTCGCCATCGCGGTCATCCGCGGCGTATGTTTGATGGCCTGCGGATAGTTCTg                                                                                                                    >  2:222926/1‑139 (MQ=255)                       gCTGTGCGAGTAATGGTGCGATATAACTCATTGAAATTGCCAAACAGATCGCGATTAAACAGACGAAACAGAGCAATATTCGCCATCGCGGTCATCCGCGGCGTATGTTTGATGGCCTGCGGATAGTTCTGATGGCTGAc                                                                                                           >  2:105411/1‑140 (MQ=255)                        cTGTGCGAGTAATGGTGCGATATAACTCATTGAAATTGCCAAACAGATCGCGATTAAACAGACGAAACAGAGCAATATTCGCCATCGCGGTCATCCGCGGCGTATGTTTGATGGCCTGCGGATAGTTCTGATGGCTGACg                                                                                                          <  1:105411/140‑1 (MQ=255)                        cTGTGCGAGTAATGGTGCGATATAACTCATTGAAATTGCCAAACAGATCGCGATTAAACAGACGAAACAGAGCAATATTCGCCATCGCGGTCATCCGCGGCGTATGTTTGATGGCCTGCGGATAGTTCTGATGGCTGACg                                                                                                          <  1:294164/140‑1 (MQ=255)                                   aTGGTGCGATATAACTCATTGAAATTGCCAAACAGATCGCGATTAAACAGACGAAACAGAGCAATATTCGCCATCGCGGTCATCCGCGGCGTATGTTTGATGGCCTGCGGATAGTTCTGATGGCTGACGCTCTCTTTgcc                                                                                               >  2:508420/1‑140 (MQ=255)                                     ggTGCGATATAACTCATTGAAATTGCCAAACAGATCGCGATTAAACAGACGAAACAGAGCAATATTCGCCATCGCGGTCATCCGCGGCGTATGTTTGATGGCCTGCGGATAGTTCTGATGGCTGACGCTCTCTTTgccgc                                                                                             >  2:468480/1‑140 (MQ=255)                                        gCGATATAACTCATTGAAATTGCCAAACAGATCGCGATTAAACAGACGAAACAGAGCAATATTCGCCATCGCGGTCATCCGCGGCGTATGTTTGATGGCCTGCGGATAGTTCTGATGGCTGACGCTCTCTTTGCCGCCtt                                                                                          <  2:74880/140‑1 (MQ=255)                                            tataACTCATTGAAATTGCCAAACAGATCGCGATTAAACAGACGAAACAGAGCAATATTCGCCATCGCGGTCATCCGCGGCGTATGTTTGATGGCCTGCGGATAGTTCTGATGGCTGACGCTCTCTTTGCCGCCTTTGcc                                                                                      >  1:444136/1‑140 (MQ=255)                                             ataACTCATTGAAATTGCCAAACAGATCGCGATTAAACAGACGAAACAGAGCAATATTCGCCATCGCGGTCATCCGCGGCGTATGTTTGATGGCCTGCGGATAGTTCTGATGGCTGACGCTCTCTTTGCCGCCTTTGCCg                                                                                     <  1:148923/140‑1 (MQ=255)                                             ataACTCATTGAAATTGCCAAACAGATCGCGATTAAACAGACGAAACAGAGCAATATTCGCCATCGCGGTCATCCGCGGCGTATGTTTGATGGCCTGCGGATAGTTCTGATGGCTGACGCTCTCTTTGCCGCCTTTGCCg                                                                                     <  2:444136/140‑1 (MQ=255)                                                  tCATTGAAATTGCCAAACAGATCGCGATTAAACAGACGAAACAGAGCAATATTCGCCATCGCGGTCATCCGCGGCGTATGTTTGATGGCCTGCGGATAGTTCTGATGGCTGACGCTCTCTTTGCCGCCt                                                                                           >  1:107505/1‑129 (MQ=255)                                                  tCATTGAAATTGCCAAACAGATCGCGATTAAACAGACGAAACAGAGCAATATTCGCCATCGCGGTCATCCGCGGCGTATGTTTGATGGCCTGCGGATAGTTCTGATGGCTGACGCTCTCTTTGCCGCCt                                                                                           <  2:107505/129‑1 (MQ=255)                                                   cATTGAAATTGCCAAACAGATCGCGATTAAACAGACGAAACAGAGCAATATTCGCCATCGCGGTCATCCGCGGCGTATGTTTGATGGCCTGCGGATGGTTCTGATGGCTGACGCTCTCTTTGCCGCCTTTGCCGTAACCg                                                                               <  2:298853/140‑1 (MQ=255)                                                     ttGAAATTGCCAAACAGATCGCGATTAAACAGACGAAACAGAGCAATATTCGCCATCGCGGTCATCCGCGGCGTATGTTTGATGGCCTGCGGATAGTTCTGATGGCTGACGCTCTCTTTGCCGCCTTTGCCGTAACCGTa                                                                             >  1:20436/1‑140 (MQ=255)                                                      tGAAATTGCCAAACAGATCGCGATTAAACAGACGAAACAGAGCAATATTCGCCATCGCGGTCATCCGCGGCGTATGTTTGATGGCCTGCGGATAGTTCTGATGGCTGACGCTCTCTTTGCCGCCTTTGCCGTAACCGTAg                                                                            >  1:128894/1‑140 (MQ=255)                                                      tGAAATTGCCAAACAGATCGCGATTAAACAGACGAAACAGAGCAATATTCGCCATCGCGGTCATCCGCGGCGTATGTTTGATGGCCTGCGGATAGTTCTGATGGCTGACGCTCTCTTTGCCGCCTTTGCCGTAACCGTAg                                                                            >  2:507133/1‑140 (MQ=255)                                                      tGAAATTGCCAAACAGATCGCGATTAAACAGACGAAACAGAGCAATATTCGCCATCGCGGTCATCCGCGGCGTATGTTTGATGGCCTGCGGATAGTTCTGATGGCTGACGCTCTCTTTGCCGCCTTTGCCGTAACCGTAg                                                                            >  1:412484/1‑140 (MQ=255)                                                      tGAAATTGCCAAACAGATCGCGATTAAACAGACGAAACAGAGCAATATTCGCCATCGCGGTCATCCGCGGCGTATGTTTGATGGCCTGCGGATAGTTCTGATGGCTGACGCTCTCTTTGCCGCCTTTGCCGTAACCGTAg                                                                            >  2:112983/1‑140 (MQ=255)                                                             gacAAACAGATCGCGATTAAACAGACGAAACAGAGCAATATTCGCCATCGCGGTCATCCGCGGCGTATGTTTGATGGCCTGCGGATAGTTCTGATGGCTGACGCTCTCTTTGCCGCCTTTGCCGTa                                                                                   <  1:50617/124‑1 (MQ=255)                                                             gCCAAACAGATCGCGATTAAACAGACGAAACAGAGCAATATTCGCCATCGCGGTCATCCGCGGCGTATGTTTGATGGCCTGCGGATAGTTCTGATGGCTGACGCTCTCTTTGCCGCCTTTGCCGTa                                                                                   >  2:50617/1‑126 (MQ=255)                                                                   cAGATCGCGATTAAACAGACGAAACAGAGCAATATTCGCCATCGCGGTCATCCGCGGCGTATGTTTGATGGCCTGCGGATAGTTCTGATGGCTGACGCTCTCTTTGCCGCCTTTGCCGTAACCGTAGACAAAATGTGACg                                                               <  1:507133/140‑1 (MQ=255)                                                                   cAGATCGCGATTAAACAGACGAAACAGAGCAATATTCGCCATCGCGGTCATCCGCGGCGTATGTTTGATGGCCTGCGGATAGTTCTGATGGCTGACGCTCTCTTTGCCGCCTTTGCCGTAACCGTAGACAAAATGTGACg                                                               <  1:198529/140‑1 (MQ=255)                                                                      aTCGCGATTAAACAGACGAAACAGAGCAATATTCGCCATCGCGGTCATCCGCGGCGTATGTTTGATGGCCTGCGGATAGTTCTGATGGCTGACGCTCTCTTTGCCGCCTTTGCCGTAACCGTAGACAAAATGTGACGCCa                                                            <  1:112983/140‑1 (MQ=255)                                                                            aTTAAACAGACGAAACAGAGCAATATTCGCCATCGCGGTCATCCGCGGCGTATGTTTGATGGCCTGCGGATAGTTCTGATGGCTGACGCTCTCTTTGCCGCCTTTGCCGTAACCGTAGACAAAATGTGACGCCAGCCCCt                                                      >  1:407415/1‑140 (MQ=255)                                                                                                      tCGCCATCGCGGTCATCCGCGGCGTATGTTTGATGGCCTGCGGATAGTTCTGATGGCTGACGCTCTCTTTGCCGCCTTTGCCGTAACCGTAGACAAAATGTGACGCCAGCCCCTGTTGCAGCGCACGCTGGTGGAGATCt                            <  2:20436/140‑1 (MQ=255)                                                                                                             cgcgGTCATCCGCGGCGTATGTTTGATGGCCTGCGGATAGTTCTGATGGCTGACGCTCTCTTTGCCGCCTTTGCCGTAACCGTAGACAAAATGTGACGCCAGCCCCTGTTGCAGCGCACGCtggtg                                   >  2:365848/1‑126 (MQ=255)                                                                                                             cgcgGTCATCCGCGGCGTATGTTTGATGGCCTGCGGATAGTTCTGATGGCTGACGCTCTCTTTGCCGCCTTTGCCGTAACCGTAGACAAAATGTGACGCCAGCCCCTGTTGCAGCGCACGCtggtg                                   <  1:365848/126‑1 (MQ=255)                                                                                                                        gcggcgTATGTTTGATGGCCTGCGGATAGTTCTGATGGCTGACGCTCTCTTTGCCGCCTTTGCCGTAACCGTAGACAAAATGTGACGCCAGCCCCTGTTGCAGCGCACGCTGGTGGAGATCTAACGCCACACCTGCTGcc          <  1:508420/140‑1 (MQ=255)                                                                                                                         cggcgTATGTTTGATGGCCTGCGGATAGTTCTGATGGCTGACGCTCTCTTTGCCGCCTTTGCCGTAACCGTAGACAAAATGTGACGCCAGCCCCTGTTGCAGCGCACGCTGGTGGAGATCTAACGCCACACCTGCTGccc         <  1:468480/140‑1 (MQ=255)                                                                                                                            cgTATGTTTGATGGCCTGCGGATAGTTCTGATGGCTGACGCTCTCTTTGCCGCCTTTGCCGTAACCGTAGACAAAATGTGACGCCAGCCCCTGTTGCAGCGCACGCTGGTGGAGATCTAACGCCACACCTGCTGCCCCGc      >  1:93961/1‑140 (MQ=255)                                                                                                                              tATGTTTGATGGCCTGCGGATAGTTCTGATGGCTGACGCTCTCTTTGCCGCCTTTGCCGTAACCGTAGACAAAATGTGACGCCAGCCCCTGTTGCAGCGCACGCTGGTGGAGATCTAACGCCACACCTGCTGCCCCGCCt    >  1:17721/1‑140 (MQ=255)                                                                                                                                tGTTTGATGGCCTGCGGATAGTTCTGATGGCTGACGCTCTCTTTGCCGCCTTTGCCGTAACCGTAGACAAAATGTGACGCCAGCCCCTGTTGCAGCGCACGCTGGTGGAGATCTAACGCCACACCTGCTGCCCCGCCTTc  >  1:340184/1‑140 (MQ=255)                                                                                                                                         |                                                                                                                                   AAAATGCAGGACCACCGGACCCGCTGTGCGAGTAATGGTGCGATATAACTCATTGAAATTGCCAAACAGATCGCGATTAAACAGACGAAACAGAGCAATATTCGCCATCGCGGTCATCCGCGGCGTATGTTTGATGACCTGCGGATAGTTCTGATGGCTGACGCTCTCTTTGCCGCCTTTGCCGTAACCGTAGACAAAATGTGACGCCAGCCCCTGTTGCAGCGCACGCTGGTGGAGATCTAACGCCACACCTGCTGCCCCGCCTTC  >  NC_000913/2131771‑2132037

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