Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,971,690	A→C	Q144P (CAG→CCG)	wecC →	UDP‑N‑acetyl‑D‑mannosaminuronic acid dehydrogenase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,971,690	0	A	C	100.0%	216.1 / NA	61	Q144P (CAG→CCG)	wecC	UDP‑N‑acetyl‑D‑mannosaminuronic acid dehydrogenase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base C (26/35);  total (26/35)

CGCGCTGGTGATCCTTGAATCCACCTCGCCGGTGGGGTCAACCGAGAAGATGGCAGAATGGTTAGCAGAGATGCGTCCGGATCTCACTTTCCCGCAGCAGGTGGGCGAGCAGGCGGACGTCAACATTGCTTACTGCCCGGAACGCGTGTTACCAGGACAGGTAATGGTCGAGCTGATTAAAAACGATCGCG  >  NC_000913/3971595‑3971785                                                                                                |                                                                                                cgcgcTGGTGATCCTTGAATCCACCTCGCCGGTGGGGTCAACCGAGAAGATGGCAGAATGGTTAGCAGAGATGCGTCCGGATCTCACTTTCCCGccgcag                                                                                             <  1:1389416/100‑1 (MQ=255)        gTGATCCTTGAATCCACCTCGCCGGTGGGGTCAACCGAGAAGATGGCAGAATGGTTAGCAGAGATGCGTCCGGATCTCACTTTCCCGCCGCAGGTGGGCg                                                                                      <  2:384480/100‑1 (MQ=255)        gTGATCCTTGAATCCACCTCGCCGGTGGGGTCAACCGAGAAGATGGCAGAATGGTTAGCAGAGATGCGTCCGGATCTCACTTTCCCGCCGCAGGTGGGCg                                                                                      <  1:731224/100‑1 (MQ=255)            tCCTTGAATCCACCTCGCCGGTGGGGTCAACCGAGAAGATGGCAGAATGGTTAGCAGAGATGCGTCCGGATCTCACTTTCCCGCCGCAGGTGGGCGAGCa                                                                                  <  1:1182708/100‑1 (MQ=255)             ccTTGAATCCACCTCGCCGGTGGGGTCAACCGAGAAGATGGCAGAATGGTTAGCAGAGATGCGTCCGGCTCTCACTTTCCCGCCGCAGGTGGGCGAGCAg                                                                                 >  2:376905/1‑100 (MQ=255)             ccTTGAATCCACCTCGCCGGTGGGGTCAACCGAGAAGATGGCAGAATGGTTAGCAGAGATGCGTCCGGATCTCACTTTCCCGCCGCAGGTGGGCGAGCAg                                                                                 >  1:1373497/1‑100 (MQ=255)               ttGAATCCACCTCGCCGGTGGGGTCAACCGAGAAGATGGCAGAATGGTTAGCAGAGATGCGTCCGGATCTCACTTTCCCGCCGCAGGTGGGCGAGCAgg                                                                                <  2:509195/99‑1 (MQ=255)                 gAATCCACCTCGCCGGTGGGGTCAACCGAGAAGATGGCAGAATGGTTAGCAGAGATGCGTCCGGATCTCACTTTCCCGCCGCAGGTGGGCGAGCAGGCgg                                                                             <  1:933944/100‑1 (MQ=255)                  aaTCCACCTCGCCGGTGGGGTCAACCGAGAAGATGGCAGAATGGTTAGCAGAGATGCGTCCGGATCTCACTTTCCCGCCGCAGGTGGGCGAGCAGGCgg                                                                             <  1:1446253/99‑1 (MQ=255)                     ccACCTCGCCGGTGGGGTCAACCGAGAAGATGGCAGAAGGGTTAGCAGAGATGCGTCCGGATCTCACTTTCCCGCCGCAGGTGGGAGAGCAGGAGGACGt                                                                         >  2:1404276/1‑100 (MQ=255)                      cACCTCGCCGGTGGGGTCAACCGAGAAGATGGCAGAATGGTTAGCAGAGATGCGTCCGGATCTCACTTTCCCGCCGCAGGTGGGCGAGCAGGCGGACGTc                                                                        >  1:1544021/1‑100 (MQ=255)                        ccTCGCCGGTGGGGTCAACCGAGAAGATGGCAGAATGGTTAGCAGAGATGCGTCCGGATCTCACTTTCCCGCCGCAGGTGGGCGAGCAGGCGGACGTCaa                                                                      >  1:653089/1‑100 (MQ=255)                          tCGCCGGTGGGGTCAACCGAGAAGATGGCAGAATGGTTAGCAGAGATGCGTCCGGATCTCACTTTCCCGCCGCAGGTGGGCGAGCAGGCGGACGTCAACa                                                                    <  2:631505/100‑1 (MQ=255)                              cGGTGGGGTCAACCGAGAAGATGGCAGAATGGTTAGCAGAGATGCGTCCGGATCTCACTTTCCCGCCGCAGGTGGGCGAGCAGGTGGACGTCAACATTGc                                                                <  1:1070941/100‑1 (MQ=255)                               ggTGGGGTCAACCGAGAAGATGGCAGAATGGTTAGCAGAGATGCGTCCGGATCTCACTTTCCCGCCGCAGGTGGGCGAGCAGGCGGACGTCAACATTGCt                                                               >  1:1321788/1‑100 (MQ=255)                                gTGGGGTCAACCGAGAAGATGGCAGAATGGTTAGCAGAGATGCGTCCGGATCTCACTTTCCCGCCGCAGGTGGGCGAGCAGGCGGACGTCAACATTGCtt                                                              <  1:424572/100‑1 (MQ=255)                                    ggTCAACCGAGAAGATGGCAGAATGGTTAGCAGAGATGCGTCCGGATCTCACTTTCCCGCCGCAGGTGGGCGAGCAGGCGGACGTCAACATTGCTTACTg                                                          <  2:985148/100‑1 (MQ=255)                                     gTCAACCGAGAAGATGGCAGAATGGTTAGCAGAGATGCGTCCGGATCTCACTTTCCCGCCGCAGGTGGGCGAGCAGGCGGACGTCAACATTGCTTACTGc                                                         <  2:461870/100‑1 (MQ=255)                                          ccGAGAAGATGGCAGAATGGTTAGCAGAGATGCGTCCGGATCTCACTTTCCCGCCGCAGGTGGGCGAGCAGGCGGACGTCAACATTGCTTACTGCCCGGa                                                    <  2:751458/100‑1 (MQ=255)                                          ccGAGAAGATGGCAGAATGGTTAGCAGAGATGCGTCCGGATCTCACTTTCCCGCCGCAGGTGGGCGAGCAGGCGGACGTCAACATTGCTTACTGCCCGGa                                                    >  1:239071/1‑100 (MQ=255)                                              gaagaTGGCAGAATGGTTAGCAGAGATGCGTCCGGATCTCACTTTCCCGCCGCAGGTGGGCGAGCAGGCGGACGTCAACATTGCTTACTGCCCGGAAcgc                                                >  2:66170/1‑100 (MQ=255)                                                agaTGGCAGAATGGTTAGCAGAGATGCGTCCGGATCTCACTTTCCCGCCGCAGGTGGGCGAGCAGGCGGACGTCAACATTGCTTACTGCCCGGAACGCgt                                              >  1:1384172/1‑100 (MQ=255)                                                 atctgCAGAATGGTTAGCAGAGATGCGTCCGGATCTCACTTTCCCGCCGCAGGTGGGCGAGCAGGCGGACGTCAACATTGCTTACTGCCCGGAACGCgtg                                             <  1:799156/96‑1 (MQ=255)                                                   tGGCAGAATGGTTAGCAGAGATGCGTCCGGATCTCACTTTCCCGCCGCAGGTGGGCGAGCAGGCGGACGTCAACATTGCTTACTGCCCGGAACGCGTGtt                                           >  1:696572/1‑100 (MQ=255)                                                   tGGCAGAATGGCTAGCAGAGATGCGTCCGGCTCTCACTATCCCGCCGCCGGTGGGCGGGCAGGCGGACGTCAACATTGCTTACTGATCGGACCGCgtg                                             <  2:1202364/98‑1 (MQ=255)                                                     gCAGAATGGTTAGCAGAGATGCGTCCGGATCTCACTTTCCCGCCGCAGGTGGGCGAGCAGGCGGACGTCAACATTGCTTACTGCCCGGAACGCgtgagat                                         >  2:799151/1‑96 (MQ=255)                                                     gCAGAATGGTTAGCAGAGATGCGTCCGGATCTCACTTTCCCGCCGCAGGTGGGCGAGCAGGCGGACGTCAACATTGCTTACTGCCCGGAACGCGTGTTAc                                         <  1:106116/100‑1 (MQ=255)                                                      cAGAATGGTTAGCAGAGATGCGTCCGGATCTCACTTTCCCGCCGCAGGTGGGCGAGCAGGCGGACGTCAACATTGCTTACTGCCCGGAACGCGTGTTAcc                                        >  2:1502740/1‑100 (MQ=255)                                                       aGAATGGTTAGCAGAGATGCGTCCGGATCTCACTTTCCCGCCGCAGGTGGGCGAGCAGGCGGACGTCAACATTGCTTACTGCCCGGAACGCGTGTTACCa                                       <  1:50654/100‑1 (MQ=255)                                                        gAATGGTTAGCAGAGATGCGTCCGGATCTCACTTTCCCGCCGCAGGTGGGCGAGCAGGCGGACGTCAACATTGCTTACTGCCCGGAACGCGTGTTACCAg                                      >  2:442746/1‑100 (MQ=255)                                                          aTGGTTAGCAGAGATGCGTCCGGATCTCACTTTCCCGCCGCAGGTGGGCGAGCAGGCGGACGTCAACATTGCTTACTGCCCGGAACGCGTGTTACCAgg                                     <  2:757003/99‑1 (MQ=255)                                                           tGGTTAGCAGAGATGCGTCCGGATCTCACTTTCCCGCCGCAGGTGGGCGAGCAGGCGGACGTCAACATTGCTTACTGCCCGGAACGCGTGTTACCAGGAc                                   <  2:1544011/100‑1 (MQ=255)                                                             gTTAGCAGAGATGCGTCCGGATCTCACTTTCCCGCCGCAGGTGGGCGAGCAGGCGGACGTCAACATTGCTTACTGCCCGGAACGCGTGTTACCAGGACAg                                 <  2:1170814/100‑1 (MQ=255)                                                               tAGCAGAGATGCGTCCTGATCTCACTTTCCCGCCGCAGGTGGGCGAGCAGGCGGACGTCAACATTGCTTACTGCCCGGAACGCGTGTTACCAGGACAGGt                               <  2:920070/100‑1 (MQ=255)                                                               tAGCAGAGATGCGTCCCGATCTCACTTTCCCGCCGCAGGTGGGCGAGCAGGCGGACGTCAACATTGCTTACTGCCCGGAACGCGTGTTACCAGGACAGGt                               <  2:919960/100‑1 (MQ=255)                                                                 gCAGAGATGCGTCCGGATCTCACTTTCCCTCCGCAGGTGGGCGAGCAGGCGGACGTCAACATGGCTTACTGCCCGGAACGCGTGTTACCAGGACAgg                                <  2:121650/97‑1 (MQ=255)                                                                 gCAGAGATGCGTCCGGATCTCACTTTCCCGCCGCAGGTGGGCGAGCAGGCGGACGTCAACATTGCTTACTGCCCGGAACGCGTGTTACCAGGACAGGTaa                             >  2:584107/1‑100 (MQ=255)                                                                  cAGAGATGCGTCCGGATCTCACTTTCCCGCCGCAGGTGGGCGAGCAGGCGGACGTCAACATTGCTTACTGCCCGGAACGCGTGTTACCAGGACAGGTAAt                            <  2:2682/100‑1 (MQ=255)                                                                    gagaTGCGTCCGGATCTCACTTTCCCGCCGCAGGTGGGCGAGCAGGCGGACGTCAACATTGCTTACTGCCCGGAACGCGTGTTACCAGGACAGGTAATgg                          >  2:1145049/1‑100 (MQ=255)                                                                    gagaTGCGTCCGGATCTCACTTTCCCGCCGCAGGTGGGCGAGCAGGCGGACGTCAACATTGCTTACTGCCCGGAACGCGTGTTACCAGGACAGGTAATgg                          <  1:66170/100‑1 (MQ=255)                                                                     agaTGCGTCCGGATCTCACTTTCCCGCCGCAGGTGGGCGAGCAGGCGGACGTCAACATTGCTTACTGCCCGGAACGCGTGTTACCAGGACAGGTAATGGt                         >  2:469965/1‑100 (MQ=255)                                                                      gaTGCGTCCGGATCTCACTTTCCCGCCGCAGGTGGGCGAGCAGGCGGACGTCAACATTGCTTACTGCCCGGAACGCGTGTTACCAGGACAGGTAATGGTc                        >  1:841662/1‑100 (MQ=255)                                                                       aTGCGTCCGGATCTCACTTTCCCGCCGCAGGTGGGCGAGCAGGCGGACGTCAACATTGCTTACTGCCCGGAACGCGTGTTACCAGGACAGGTAATGGTCg                       <  2:422959/100‑1 (MQ=255)                                                                        tGCGTCGGGATCTCACTTTCCCGCCGCAGGTGGGCGGGCAGGCGGACGTCAACATTGGTTACTGCCCGGAACGCGTGTTACCAGGACAGGTAATGGTCg                       <  2:1486135/99‑1 (MQ=255)                                                                        tGCGTCCGGATCTCACTTTCCCGCCGCAGGTGGGCGAGCAGGCGGACGTCAACATTGCTTACTGCCCGGAACGCGTGTTACCAGGACAGGTAATGGTCg                       >  2:507366/1‑99 (MQ=255)                                                                         gCGTCCGGATCTCACTTTCCCGCCGCAGGTGGGCGAGCAGGCGGACGTCAACATTGCTTACTGCCCGGAACGCGTGTTACCAGGACAGGTAATGGTCGAg                     <  2:1321779/100‑1 (MQ=255)                                                                          cGTCCGGATCTCACTTTCCCGCCGCAGGTGGGCGAGCAGGCGGACGTCAACATTGCTTACTGCCCGGAACGCGTGTTACCAGGACAGGTAATGGTCGAGc                    >  1:370340/1‑100 (MQ=255)                                                                           gTCCGGATCTCACTTTCCCGCCGCAGGTGGGCGAGCAGGCGGACGTCAACATTGCTTACTGCCCGGAACGCGTGTTACCAGGACAGGTAATGGTCGAGCt                   <  1:1355379/100‑1 (MQ=255)                                                                           gTCCGGATCTCACTTTCCCGCCGCAGGTGGGCGAGCAGGCGGACGTCAACATTGCTTACTGCCCGGAACGCGTGTTACCAGGACAGGTAATGGTCGAGCt                   >  1:120274/1‑100 (MQ=255)                                                                             ccGGATCTCACTTTCCCGCCGCAGGTGGGCGAGCAGGCGGACGTCAACATTGCTTACTGCCCGGAACGCGTGTTACCAGGACAGGTAATGGTCGAGCTGa                 <  1:1238924/100‑1 (MQ=255)                                                                              cGGATCTCACTTTCCCGCCGCAGGTGGGCGAGCAGGCGGACGTCAACATTGCTTACTGCCCGGAACGCGTGTTACCAGGACAGGTAATGGTCGAGCTGAt                >  2:1421615/1‑100 (MQ=255)                                                                               ggATCTCACTTTCCCGCCGCAGGTGGGCGCGCAGGCGGACGTCAACATTGCTTACTGCCCGGGACGCGTGTTACCCGGACAGGTAATGGTCGAGCTGAtt               >  2:310564/1‑100 (MQ=255)                                                                               ggATCTCACTTTCCCGCCGCAGGTGGGCGAGCAGGCGGACGTCAACATTGCTTACTGCCCGGAACGCGTGTTACCAGGACAGGTAATGGTCGAGCTGAtt               >  2:251961/1‑100 (MQ=255)                                                                               ggATCTCACTTTCCCGCCGCAGGTGGGCGAGCAGGCGGACGTCAACATTGCTTACTGCCCGGAACGCGTGTTACCAGGACAGGTAATGGTCGAGCTGAtt               <  1:7480/100‑1 (MQ=255)                                                                               ggATCTCACTTTCCCGCCGCAGGTGGGCGAGCAGGCGGACGTCAACATTGCTTACTGCCCGGAACGCGTGTTACCAGGACAGGTAATGGTCGAGCTGAtt               <  1:237377/100‑1 (MQ=255)                                                                                gATCTCACTTTCCCGCCGCAGGTGGGCGAGCAGGCGGACGTCAACATTGCTTACTGCCCGGAACGCGTGTTACCAGGACAGGTAATGGTCGAGCTGATTa              <  1:1211417/100‑1 (MQ=255)                                                                                 aTCTCACTTTCCCGCCGCAGGTGGGCGAGCAGGCGGACGTCAACATTGCTTACTGCCCGGAACGCGTGTTACCAGGACAGGTACTGGTCGAGCTCATTa              >  2:757250/1‑99 (MQ=255)                                                                                 aTCTCACTTTCCCGCCGCAGGTGGGCGAGCAGGCGGACGTCAACATTGCTTACTGCCCGGAACGCGTGTTACCAGGACAGGTAATGGTCGAGCTGATTaa             >  2:384775/1‑100 (MQ=255)                                                                                       cTTTCCCGCCGCAGGTGGGCGAGCAGGCGGACGTCAACATTGCTTACTGCCCGGAACGCGTGTTACCAGGACAGGTAATGGTCGAGCTGATTAAAAACGa       >  1:279924/1‑100 (MQ=255)                                                                                        tttCCCGCCGCAGGTGGGCGAGCAGGCGGACGTCAACATTGATTACTGCCCGGAACGCGTGTTACCAGGACAGGTAATGGTCGAGCTGATTAAAAACGAt      <  2:702991/100‑1 (MQ=255)                                                                                         ttCCCGCCGCAGGTGGGCGAGCAGGCGGACGTCAACATTGCTTACTGCCCGGAACGCGTGTTACCAGGACAGGTAATGGTCGAGCTGATTAAAAACGATc     <  1:876620/100‑1 (MQ=255)                                                                                         ttCCCCCCGCAGGTGGGCGAGCAGGCGGACGTCAACATTGCTTACTGCCCGGAACGCGTGTTACCAGGACAGGTAATGGTCGAGCTGATTAAAAACGATc     <  1:1404286/100‑1 (MQ=255)                                                                                            ccgccgCAGGTGGGCGAGCAGGCGGACGTCAACATTGCTTACTGCCCGGAACGCGTGTTACCAGGACAGGTAATGGTCGAGCTGATTAAAAACGCTcgcg  >  2:167523/1‑100 (MQ=255)                                                                                                |                                                                                                CGCGCTGGTGATCCTTGAATCCACCTCGCCGGTGGGGTCAACCGAGAAGATGGCAGAATGGTTAGCAGAGATGCGTCCGGATCTCACTTTCCCGCAGCAGGTGGGCGAGCAGGCGGACGTCAACATTGCTTACTGCCCGGAACGCGTGTTACCAGGACAGGTAATGGTCGAGCTGATTAAAAACGATCGCG  >  NC_000913/3971595‑3971785

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