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%	166.6 / NA	47	Q144P (CAG→CCG)	wecC	UDP‑N‑acetyl‑D‑mannosaminuronic acid dehydrogenase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base C (24/23);  total (24/23)

CGCTGGTGATCCTTGAATCCACCTCGCCGGTGGGGTCAACCGAGAAGATGGCAGAATGGTTAGCAGAGATGCGTCCGGATCTCACTTTCCCGCAGCAGGTGGGCGAGCAGGCGGACGTCAACATTGCTTACTGCCCGGAACGCGTGTTACCAGGACAGGTAATGGTCGAGCTGATTAAAAACGATCGCGT  >  NC_000913/3971597‑3971786                                                                                              |                                                                                                 cgcTGGTGATCCTTGAATCCACCTCGCCGGTGGGGTCAACCGAGAAGATGGCAGAATGGTTAGCAGAGATGCGTCCGGATCTCACTTTCCCGCCGCAGGt                                                                                            >  1:526072/1‑100 (MQ=255)  gcTGGTGATCCTTGAATCCACCTCGCCGGTGGGGTCAACCGAGAAGATGGCAGAATGGTTAGCAGAGATGCGTCCGGATCTCACTTTCCCGCCGCAGGTg                                                                                           <  1:810737/100‑1 (MQ=255)    tGGTGATCCTTGAATCCACCTCGCCGGTGGGGTCAACCGAGAAGATGGCAGAATGGTTAGCAGAGATGCGTCCGGATCTCACTTTCCCGCCGCAGGTggg                                                                                         <  1:57286/100‑1 (MQ=255)       tGATCCTTGAATCCACCTCGCCGGTGGGGTCAACCGAGAAGATGGCAGAATGGTTAGCAGAGATGCGTCCGGATCTCACTTTCCCGCCGCAGGTGGGCGa                                                                                      <  1:888138/100‑1 (MQ=255)       tGATCCTTGAATCCACCTCGCCGGTGGGGTCAACCGAGAAGATGGCAGAATGGTTAGCAGAGATGCGTCCGGATCTCACTTTCCCGCCGCAGGTGGGCGa                                                                                      >  1:751861/1‑100 (MQ=255)          tCCTTGAATCCACCTCGCCGGTGGGGTCAACCGAGAAGATGGCAGAATGGTTAGCAGAGATGCGTCCGGATCTCACTTTCCCGCCGCAGGTGGGCGAGCa                                                                                   >  1:621951/1‑100 (MQ=255)          tCCTTGAATCCACCTCGCCGGTGGGGTCAACCGAGAAGATGGCAGAATGGTTAGCAGAGATGCGTCCGGATCTCACTTTCCCGCCGCAGGTGGGCGAGCa                                                                                   >  1:715946/1‑100 (MQ=255)          tCCTTGAATCCACCTCGCCGGTGGGGTCAACCGAGAAGATGGCAGAATGGTTAGCAGAGATGCCTCCGGATCTCACTTTCCCGCCGCAGGTGGGCGAGCa                                                                                   <  1:918067/100‑1 (MQ=255)              tGAATCCACCTCGCCGGTGGGGTCAACCGAGAAGATGGCAGAATGGTTAGCAGAGATGCGTCCGGATCTCACTTTCCCGCCGCAGGTGGGCGAGCAGGCg                                                                               <  2:621947/100‑1 (MQ=255)              tGAATCCACCTCGCCGGTGGGGTCAACCGAGAAGATGGCAGAATGGTTAGCAGAGATGCGTCCGGATCTCACTTTCCCGCCGCAGGTGGGCGAGCAGGCg                                                                               >  2:882534/1‑100 (MQ=255)               gAATCCCCATCGCCGGTGGGGTCAACCGAGAAGATGGCAGAATGGTTAGCAGAGATGCGTCCGGATCTCACTTTCCCGCCGCAGGTGGGCGAGCAGGCgg                                                                              <  1:795899/100‑1 (MQ=255)                    caccggcgaggTGGGGTCAACCGAGAAGATGGCAGCATGGTTCGCAGAGGTGCGTCCGGATCCCACGTTCCCGCCGCAGGTGGGCGAGCAGGCGGACg                                                                           >  2:170430/10‑98 (MQ=255)                    cACCTCGCCGGTGGGGTCAACCGAGAAGATGGCAGAATGGTTAGCAGAGATGCGTCCGGATCTCACTTTCCCGCCGCAGGTGGGCGAGCAGGCGGACGTc                                                                         <  2:872156/100‑1 (MQ=255)                       cTCGCCGGTGGGGTCAACCGAGAAGATGGCAGAATGGTTCGCAGGGATGCGTCCGGATCTCACTTTCCCGCCGCAGGTGGGCGCGCAGGCGGACGTCaa                                                                       >  2:876313/1‑99 (MQ=255)                       cTCGCCGGTGGGGTCAACCGAGAAGATGGCAGAATGGTTAGCAGAGATGCGTCCGGATCTCACTTTCCCGCCGCAGGTGGGCGAGCAGGCGGACGTCAAc                                                                      >  1:1276905/1‑100 (MQ=255)                            cGGTGGGGTCAACCGAGAAGATGACAGAATGGTTAGCAGAGATGCGTCCGGATCTCACTTTCCCGCCGCAGGTGGGCGAGCAGGCGGACGTCAACATTGc                                                                 >  1:1131690/1‑100 (MQ=255)                             ggTGGGGTCAACCGAGAAGATGGCAGAATGGTTAGCAGAGATGCGTCCGGATCTCACTTTCCCGCCGCAGGTGGGCGAGCAGGCGGACGTCAACATTGCt                                                                >  1:999007/1‑100 (MQ=255)                             ggTGGGGTCAACCGAGAAGATGGCAGAATGGTTAGCAGAGATGCGTCCGGATCTCACTTTCCCGCCGCAGGTGGGCGAGCAGGCGGACGTCAACATTGCt                                                                <  1:898698/100‑1 (MQ=255)                                 gggTCAACCGAGAAGATGGCAGAATGGTTAGCAGAGATGCGTCCGGATCTCACTTTCCCGCCGCAGGTGGGCGAGCAGGCGGACGTCAACATTGCTTACt                                                            >  1:63559/1‑100 (MQ=255)                                    tCAACCGAGAAGATGGCAGAATGGTTAGCAGAGATGCGTCCGGATCTCACTTTCCCGCCGCAGGTGGGCGAGCAGGCGGACGTCAACATTGCTTACTGcc                                                         <  2:265895/100‑1 (MQ=255)                                         cGAGAAGATGGCAGAATGGTTAGCAGAGATGCGTCCGGATCTCACTTTCCCGCCGCAGGTGGGCGAGCAGGCGGACGTCAACATTGCTTACTGCCCGGaa                                                    >  1:61899/1‑100 (MQ=255)                                             aagaTGGCAGAATGGTTAGCAGAGATGCGTCCGGATCTCACTTTCCCGCCGCAGGTGGGCGAGCAGGCGGACGTCAACATTGCTTACTGCCCGGAAcgcg                                                <  2:1147472/100‑1 (MQ=255)                                                 tGGCAGAATGGTTAGCAGAGATGCGTCCGGATCTCACTTTCCCGCCGCAGGTGGGCGAGCAGGCGGACGTCAACATTGCTTACTGCCCGGAACGCGTGtt                                            <  1:751836/100‑1 (MQ=255)                                                 tGGCAGAATGGTTAGCAGAGATGCGTCCGGATCTCACTTTCCCGCCGCAGGTGGGCGAGCAGGCGGACGTCAACATTGCTTACTGCCCGGAACGCGTGtt                                            <  2:670896/100‑1 (MQ=255)                                                    cAGAATGGTTAGCAGAGATGCGTCCGGATCTCACTTTCCCGCCGCAGGTGGGCGAGCAGGCGGACGTCAACATTGCTTACTGCCCGGAACGCGTGTTAcc                                         >  1:894248/1‑100 (MQ=255)                                                    cAGAATGGTTAGCAGAGATGCGTCCGGATCTCACTTTCCCGCCGCAGGTGGGCGAGCAGGCGGACGTCAACATTGCTTACTGCCCGGAACGCGTGTTAcc                                         <  1:729703/100‑1 (MQ=255)                                                     aGAATGGTTAGCAGAGATGCGTCCGGATCTCACTTTCCTGCCGCAGGTGGGCGAGCAGGCGGACGTCAACATTGCTTACTGCCCGGAACGCGTGTTACCa                                        >  2:805066/1‑100 (MQ=255)                                                     aGAATGGTTAGCAGAGATGCGTCCGGATCTCACTTTCCCGCCGCAGGTGGGCGAGCAGGCGGACGTCAACATTGCTTACTGCCCGGAACGCGTGTTACCa                                        <  2:1276894/100‑1 (MQ=255)                                                          ggTTAGCAGAGATGCGTCCGGATCTCACTTTCCCGCCGCAGGTGGGCGAGCAGGCGGACGTCAACATTGCTTACTGCCCGGAACGCGTGTTACCAGGACa                                   >  2:245184/1‑100 (MQ=255)                                                              aGCAGAGATGCGTCCGGATCTCACTTTCCCGCCGCAGGTGGGCGAGCAGGCGGACGTCAACATTGCTTACTGCCCGGAACGCGTGTTACCAGGACAGGt                                <  1:24648/99‑1 (MQ=255)                                                              aGCAGAGATGCGTCCGGATCTCACTTTCCCGCCGCAGGTGGGCGAGCAGGCGGACGTCAACATTGCTTACTGCCCGGAACGCGTGTTACCAGGACAGGTa                               <  1:1288307/100‑1 (MQ=255)                                                              aGCAGAGATGCGTCCGGATCTCACTTTCCCGCCGCAGGTGGGCGAGCAGGCGGACGTCAACATTGCTTACTGCCCGGAACGCGTGTTACCAGGACAGGTa                               <  1:190679/100‑1 (MQ=255)                                                                cAGAGATGCGTCCGGATCTCACTTTCCCGCCGCAGGTGGGCGAGCAGGCGGACGTCAACATTGCTTACTGCCCGGAACGCGTGTTACCAGGACAGGTAAt                             >  2:900264/1‑100 (MQ=255)                                                                 agagATGCGTCCGGATCTCACTTTCCCGCCGCAGGTGGGCGAGCAGGCGGACGTCAACATTGCTTACTGCCCGGAACGCGTGTTACCAGGACAGGTAATg                            >  2:627490/1‑100 (MQ=255)                                                                  gagaTGCGTCGGGATCTCACTTTCCCGCCGCAGGTGGGCGAGCAGGCGGACGTCAACATTGCTTACTGCCCGGAACGCGTGTTACCAGGACAGGTAATgg                           >  2:917856/1‑100 (MQ=255)                                                                   agaTGCGTCCGGATCTCACTTTCCCGCCGCAGGTGGGCGAGCAGGCGGACGTCAACATTGCTTACTGCCCGGAACGCGTGTTACCAGGACAGGTAATGGt                          <  1:1189101/100‑1 (MQ=255)                                                                       gCGTCCGGATCTCACTTTCCCGCCGCAGGTGGGCGAGCAGGCGGACGTCAACATTGCTTACTGCCCGGAACGCGTGTTACCAGGACAGGTAATGGTCGAg                      <  2:729732/100‑1 (MQ=255)                                                                         gTCCGGATCTCACTTTCCCGCCGCAGGTGGGCGAGCAGGCGGACGTCAACATTGCTTACTGCCCGGAACGCGTGTTACCAGGACAGGTAATGGTCGAGCt                    >  2:328507/1‑100 (MQ=255)                                                                          tCCGGATCTCACTTTCCCGCCGCAGGTGGGCGAGCAGGCGGACGTCAACATTGCTTACTGCCCGGAACGCGTGTTACCAGGACAGGTAATGGTCGAGCTg                   >  2:656220/1‑100 (MQ=255)                                                                             ggATCTCACTTTCCCGCCGCAGGTGGGCGAGCAGGCGGACGTCAACATTGCTTACTGCCCGGAACGCGTGTTACCAGGACAGGTAATGGTCGAGCTGAtt                <  2:531558/100‑1 (MQ=255)                                                                               aTCTCACTTTCCCGCCGCAGGTGGGCGAGCAGGCGGACGTCAACATTGCTTACTGCCCGGAACGCGTGTTACCAGGACAGGTAATGGTCGAGCTGATTaa              >  1:19750/1‑100 (MQ=255)                                                                                tctcACTTTCCCGCCGCAGGTGGGCGAGCAGGCGGACGTCAACATTGCTTACTGCCCGGAACGCGTGTTACCAGGACAGGTAATGGTCGAGCTGATTaaa             >  1:436807/1‑100 (MQ=255)                                                                                   cACTTTCCCGCCGCAGGTGGGCGAGCAGGCGGACGTCAACATTGCTTACTGCCCGGAACGCGTGTTACCAGGACAGGTAATGGTCGAGCTGATTAAAAAc          >  2:67923/1‑100 (MQ=255)                                                                                    aCTTTCCCGCCGCAGGTGGGCGAGCAGGCGGACGTCAACATTGCTTACTGCCCGGAACGCGTGTTACCAGGACAGGTAATGGTCGAGCTGATTAAAAACg         <  2:332707/100‑1 (MQ=255)                                                                                      tttCCCGCCGCAGGTGGGCGAGCAGGCGGACGTCAACATTGCTTACTGCCCGGAACGCTTGTTACCAGGACAGGTAATGGTCGAGCTGATTAAAAACGAt       <  2:712342/100‑1 (MQ=255)                                                                                      tttCCCGCCGCAGGTGGGCGAGCAGGCGGACGTCAACATTGCTTACTGCCCGGAACGCGTGTTACCAGGACAGGTAATGGTCGAGCTGATTAAAAACGAt       >  2:825343/1‑100 (MQ=255)                                                                                       ttCCCGCCGCAGGTGGGCGAGCAGGCGGACGTCAACATTGCTTACTGCCCGGAACGCGTGTTACCAGGACAGGTAATGGTCGAGCTGATTAAAAACGATc      <  1:794825/100‑1 (MQ=255)                                                                                          ccgccgCAGGTGGGCGAGCAGGCGGACGTCAACATTGCTTACTGCCCGGAACGCGTGTTACCAGGACAGGTAATGGTCGAGCTGATTAAAAACGATcgcg   >  2:290958/1‑100 (MQ=255)                                                                                           cgccgcAGGTGGGCGAGCAGGCGGACGTCAACATTGCTTACTGCCCGGAACGCGTGTTACCAGGACAGGTAATGGTCGAGCTGATTAAAAACGATCGCGt  <  2:514411/100‑1 (MQ=255)                                                                                              |                                                                                                 CGCTGGTGATCCTTGAATCCACCTCGCCGGTGGGGTCAACCGAGAAGATGGCAGAATGGTTAGCAGAGATGCGTCCGGATCTCACTTTCCCGCAGCAGGTGGGCGAGCAGGCGGACGTCAACATTGCTTACTGCCCGGAACGCGTGTTACCAGGACAGGTAATGGTCGAGCTGATTAAAAACGATCGCGT  >  NC_000913/3971597‑3971786

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