Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,132,787	A→C	I204S (ATC→AGC)	wcaA ←	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,132,787	0	A	C	100.0%	141.5 / NA	50	I204S (ATC→AGC)	wcaA	putative colanic acid biosynthesis glycosyl transferase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base C (26/24);  total (26/24)

CTGGTAGAGGGTAAACAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGATCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGCCATGCCCAGGTAAAGACCTGGTT  >  NC_000913/2132651‑2132926                                                                                                                                         |                                                                                                                                            cTGGTAGAGGGTAAACAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATGTGCAGGCTCTGCg                                                                                                                                       >  1:1046872/1‑143 (MQ=255)   ggTAGAGGGTAAACAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGCTTGATATGCAGGCTCTGCGTc                                                                                                                                     >  2:203362/1‑143 (MQ=255)    gTAGAGGGTAAACAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGCAGGACACGGCGAATTTTTTCGCTGGCTAGGTCATCTGCATTTCACCGTGATTGATACGAAGGCACTgcggcg                                                                                                                                    >  2:390118/1‑143 (MQ=255)    gTAGAGGGTAAACAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCGTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCg                                                                                                                                    >  1:639598/1‑143 (MQ=255)        agGGTAAACAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCctc                                                                                                                                >  1:906370/1‑143 (MQ=255)          ggTAAACAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCtctt                                                                                                                              >  2:594202/1‑143 (MQ=255)             aaaCAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTa                                                                                                                           <  2:874122/143‑1 (MQ=255)              aagcgCTCATATTTTTTGCTGGGGCGGTCGAATTTGTCTTTGTGCTGGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTAc                                                                                                                          <  2:969694/139‑1 (MQ=255)                cAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACtt                                                                                                                        >  1:5056/1‑143 (MQ=255)                 aGCTGATATTTTTTGCTGGCACGGTCGAATTTGTCATTGTGCTTGCGGCAAAAGTGGCAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCACTTCCCCGTGATCGATATCCAGGCTCAGCGCCCCCTCTTCCAcctc                                                                                                                       >  2:365907/1‑140 (MQ=255)                   cTGATATTTTTTGCTGGCGCGGTCGAATTTGTCGTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTcc                                                                                                                     <  1:469474/143‑1 (MQ=255)                     gATATTTTTGGCTGGCGCGGTCGAAGTTGGCTGTGGGCGTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCAc                                                                                                                   <  1:435267/143‑1 (MQ=255)                      atatTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTGTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCGTCTCCTTTCGACg                                                                                                                  >  2:571175/1‑143 (MQ=255)                        atTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACCTTCCACGGt                                                                                                                >  1:350458/1‑143 (MQ=255)                           ttttGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGCAGTACCCGGAGCATTTTTTCGGTGACGAGGTGATATGCATTTCACCGTGATTGTTATGCAGGCGCTGCGTCGCCTCTTCTACGTTCCGCGGTTCg                                                                                                             >  2:37762/1‑143 (MQ=255)                             ttGCGGGCGCGGTCGAATTTGTCTTTGTCTCTGCGGTCCAGGTGAAAGTACACGGAGAATTTTTTCGGTGACTTGGTGATGAGAAATACACTGCGATTGATATGGAGGCTCTGCGTTGCCTCTTCTTCTTTCAACggatagtc                                                                                                           <  2:1101554/143‑7 (MQ=255)                                        gTCGAGTGTGTCTTTGTGCTGGCGGTAACAGTGGAGGTACCTGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTccacca                                                                                                <  2:9724/143‑1 (MQ=255)                                        gTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGCTTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCAACGGTTCGCCGTACTccacca                                                                                                >  2:349029/1‑143 (MQ=255)                                        gTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTccacca                                                                                                <  2:1013373/143‑1 (MQ=255)                                        gTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGACTCTCCTACTTTCCACGGTTCGCCGTCCTccacca                                                                                                >  2:298216/1‑143 (MQ=255)                                          cGAATTAGACTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCATCGCCTCTTCTACTTTCCACGGTTCGCCGTACTccaccacc                                                                                              <  1:1100585/143‑1 (MQ=255)                                                tGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTCCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGc                                                                                        >  1:866614/1‑143 (MQ=255)                                                   cTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAgg                                                                                     <  1:980243/143‑1 (MQ=255)                                                          cctGCGGTAAAAGTGGAAGTACCCGTAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGt                                                                              <  1:489064/141‑1 (MQ=255)                                                           ttGCGGTAAAAGTGGAAGTACCCGGAGAATTGTTTCGGTGACGAGGTGAGCTGCATTGCACCGTGATTGATAGGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTc                                                                             <  1:83927/143‑1 (MQ=255)                                                             gCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGt                                                                           >  1:172260/1‑143 (MQ=255)                                                              cGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATCTCACCGTGATTGATATGCCGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTGCACCACCATCCGCAGGACGACGTCGTa                                                                          >  2:480633/1‑143 (MQ=255)                                                               ggTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCGCCACCATCCGCAGGAAGCTGTCGGaa                                                                         >  1:495876/1‑143 (MQ=255)                                                                 tAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATc                                                                       <  1:1011605/143‑1 (MQ=255)                                                                    aaGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGATGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCtg                                                                    >  1:720782/1‑143 (MQ=255)                                                                        ggAAGTACCCGGAGAATTGGTGCGGTGACGAGGTGAGCTGCACTGCACCGTGATTGATGTGCAGGCTCTGCGTCGCGTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCAGCCGCAGGAAGATGTCGTAATCCTGTGcc                                                                <  1:1088905/143‑1 (MQ=255)                                                                          taGTGCCCGGAGCATGTTTTCGGTGGCGAGGTGATCTGCATTTCACCGTGATTGAGATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGc                                                              <  2:703575/142‑1 (MQ=255)                                                                            ttcctcGGAGAATTTTGCCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCtt                                                            <  1:302816/138‑1 (MQ=255)                                                                                   gagaATGTGTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCGACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTc                                                     <  1:748206/143‑1 (MQ=255)                                                                                    agcaTTTTTTCGGGGACGTGGGGATCTGCTGTTCACCGTGATTTATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCg                                                    <  2:199950/140‑1 (MQ=255)                                                                                     gtcttttgtCGGTGCCGGGTTGGTCTGCGTTCCATCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACGTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCgg                                                   <  2:505103/140‑1 (MQ=255)                                                                                               gggggcGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGaacaa                                         <  2:747249/138‑1 (MQ=255)                                                                                               ggTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTACTACTTTCAACGGTTCGCAGTAGTCCACCACCATCCGCAGGAAGATGTCGTAATACTGTGCCGCGTACCGGTCGGTGTCGaacaa                                         >  2:1000245/1‑143 (MQ=255)                                                                                                 tGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGaacaaac                                       >  2:1027238/1‑143 (MQ=255)                                                                                                  gtcGAGGTGATCTCCGTTGCGCCGTGATTGATGTCCGGGCGCTGCGTCCCCTCTTCTACTTTCCACGGTGCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGaacaaaca                                      <  2:938279/141‑1 (MQ=255)                                                                                                                         gTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGCCATGCCCAgg               <  1:536264/143‑1 (MQ=255)                                                                                                                          tgattgatATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGCCATGCCCAGGt              <  1:273390/143‑1 (MQ=255)                                                                                                                            attgatATGCAGGCTCTGCGTCCCCTCTACTACTTTCCACGCTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACACACACTCTTTGAAACGCCATGCCCAGGTaa            >  2:976495/1‑143 (MQ=255)                                                                                                                             ttgatATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGCCATGCCCAGGTaaa           <  1:780386/143‑1 (MQ=255)                                                                                                                               ggcatGCAGGCTCTGCGTCGCCTCATATCCCTTCCACTGTTCGCCGGACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGCCATGCCCAGGTAAAGa         <  1:203362/140‑1 (MQ=255)                                                                                                                                atatGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCCGTGCCGCTTTCAGTTCCGTATCGAACCACCACTCTGTGAAGCGCCCTGCCCAGGtcagcgc        >  2:639767/1‑137 (MQ=255)                                                                                                                                atatGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACCGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAACCACTCTTTGAAACGCCATGCCCAGGTAAATAc        >  1:724057/1‑143 (MQ=255)                                                                                                                                    gCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAGACGCCATGCCCAGGTAAACACCTgg    >  1:796713/1‑143 (MQ=255)                                                                                                                                    gCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTCCGGGATCGCACAACCACTCTTTGAAACGCCATGCCAAGTTAAAGATCTgg    >  1:414398/1‑143 (MQ=255)                                                                                                                                      aGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGACGGTGTCGGAATACTGTGCCGCTCTCAGTTCGGTAGCGAACAGACACTCTGCGACACGCCCTGCCCAGTTAGAGAACGGGtt  >  2:80160/1‑143 (MQ=255)                                                                                                                                         |                                                                                                                                            CTGGTAGAGGGTAAACAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGATCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGCCATGCCCAGGTAAAGACCTGGTT  >  NC_000913/2132651‑2132926

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