Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	404,123	(C)5→4	coding (421/1116 nt)	yaiC →	diguanylate cyclase, cellulose regualtor

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	404,119	0	C	.	100.0%	165.1 / NA	39	coding (417/1116 nt)	yaiC	diguanylate cyclase, cellulose regualtor
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base . (18/21);  total (18/21)

TTAGCCGGGAAATTTACAACTTAAAAACCGATGCAGTATTAGCGGGAATGTGGGTAGGCGTAATGGGCGTAAACGTGCTGCCTTCCACCGCGATGTTGATGATTATGTGTCTGAATTTGATGGGGGCAGGCGGCCCCCGTCTGTTTGTCGCGGGTCTGGTGTTGATGGTGGTTTCCTGCCTTGTCACCCTCGAGCTGACGGGCATTACCGTGTCGTTCAATAGTGCGCCGCTGGAATGGTGGCTCTCCCTTCCCATTATTGTCA  >  NC_000913/403986‑404249                                                                                                                                      |                                                                                                                                   ttAGCCGGGAAATTTACAACTTAAAAACCGATGCAGTATTAGCGGGAATGTGGGTAGGCGTAATGGGCGTAAACGTGCTGCCTTCCACCGCGATGTTGATGATTATGTGTCTGAATTTGATGGGGGCAGGCGGccccgt                                                                                                                               <  2:556067/139‑3 (MQ=255)           aaTTTACAACTTAAAAACCGATGCAGTATTAGCGGGAATGTGGGTAGGCGTAATGGGCGTAAACGTGCTGCCTTCCACCGCGATGTTGATGATTATGTGTCTGAATTTGATGGGGGCAGGCGG‑CCCCGTCTGTTTGTcg                                                                                                                    <  1:560904/139‑1 (MQ=255)           aaTTTACAACTTAAAAACCGATGCAGTATTAGCGGGAATGTGGGTAGGCGTAATGGGCGTAAACGTGCTGCCTTCCACCGCGATGTTGATGATTATGTGTCTGAATTTGATGGGGGCAGGCGG‑CCCCGTCTGTTTGTcg                                                                                                                    <  2:701253/139‑1 (MQ=255)              ttACAACTTAAAAACCGATGCAGTATTAGCGGGAATGTGGGTAGGCGTAATGCGCGTAAACGTGCTGCCTTCCACCGCGATGTTGATGATTATGTGTCTGAATTTGATGGGGGCAGGCGG‑CCCCGTCTGTTTGTCGCgg                                                                                                                 <  2:734902/139‑1 (MQ=255)                          aaCCGATTCAGTATTAGCGGGAATGTGGGTAGGCGTAATGGGCGTAAACGTGCTGCCTTCCACCGCGATGTTGATGATTATGTGTCTGAATTTGATGGGGGCAGGCGG‑CCCCGTCTGTTTGTCGCGGGTCTGGTGTTGa                                                                                                     <  1:868008/139‑1 (MQ=255)                                    gTATTAGCGGGAATGTGGGTAGGCGTAATGGGCGTAAACGTGCTGCCTTCCACCGCGATGTTGATGATTATGTGTCTGAATTTGATGGGGGCAGGCGG‑CCCCGTCTGTTTGTCGCGGGTCTGGTGTTGATGGTGGTTTc                                                                                           <  1:539044/139‑1 (MQ=255)                                               aaTGTGGGTAGGCGTAATGGGCGTAAACGTGCTGCCTTCCACCGCGATGTTGATGATTATGTGTCTGAATTTGATGGGGGCAGGCGG‑CCCCGTCTGTTTGTCGCGGGTCTGGTGTTGATGGTGGTTTc                                                                                           <  1:1037175/128‑1 (MQ=255)                                               aaTGTGGGTAGGCGTAATGGGCGTAAACGTGCTGCCTTCCACCGCGATGTTGATGATTATGTGTCTGAATTTGATGGGGGCAGGCGG‑CCCCGTCTGTTTGTCGCGGGTCTGGTGTTGATGGTGGTTTc                                                                                           >  2:1037175/1‑128 (MQ=255)                                               aaTGTGGGTAGGCGTAATGGGCGTAAACGTGCTGCCTTCCACCGCGATGTTGATGATTATGTGTCTGAATTTGATGGGGGCAGGCGG‑CCCCGTCTGTTTGTCGCGGGTCTGGTGTTGATGGTGGTTTCCTGCCTTGTCa                                                                                >  2:354649/1‑139 (MQ=255)                                               aaTGTGGGTAGGCGTAATGGGCGTAAACGTGCTGCCTTCCACCGCGATGTTGATGATTATGTGTCTGAATTTGATGGGGGCAGGCGG‑CCCCGTCTGTTTGTCGCGGGTCTGGTGTTGATGGTGGTTTCCTGCCTTGTCa                                                                                >  2:240111/1‑139 (MQ=255)                                                 tgtgGGTAGGCGTAATGGGCGTAAACGTGCTGCCTTCCACCGCGATGTTGATGATTATGTGTCTGAATTTGATGGGGGCAGGCGG‑CCCCGTCTGTTTGTCGCGGGTCTGGTGTTGATGGTGGTTTCCTGCCTTGTCAcc                                                                              >  1:950224/1‑139 (MQ=255)                                                            gTAATGGGCGTAAACGTGCTGCCTTCCACCGCGATGTTGATGATTATGTGTCTGAATTTGATGGGGGCAGGCGG‑CCCCGTCTGTTTGTCGCGGGTCTGGTGTTGATGGTGGTTTCCTGCCTTGTCACCCTCGAGCTGAc                                                                   <  2:345748/139‑1 (MQ=255)                                                            gTAATGGGCGTAAACGTGCTGCCTTCCACCGCGATGTTGATGATTATGTGTCTGAATTTGATGGGGGCAGGCGG‑CCCCGTCTGTTTGTCGCGGGTCTGGTGTTGATGGTGGTTTCCTGCCTTGTCACCCTCGAGCTGAc                                                                   <  2:950224/139‑1 (MQ=255)                                                               aTGGGCGTAAACGTGCTGCCTTCCACCGCGATGTTGATGATTATGTGTCTGAATTTGATGGGGGCAGGCGG‑CCCCGTCTGTTTGTCGCGGGTCt                                                                                                             >  1:267368/1‑94 (MQ=255)                                                               aTGGGCGTAAACGTGCTGCCTTCCACCGCGATGTTGATGATTATGTGTCTGAATTTGATGGGGGCAGGCGG‑CCCCGTCTGTTTGTCGCGGGTCt                                                                                                             <  2:267368/94‑1 (MQ=255)                                                               aTGGGCGTAAACGTGCTGCCTTCCACCGCGATGTTGATGATTATGTGTCTGAATTTGATGGGGGCAGGCGG‑CCCCGTCTGTTTGTCGCGGGTCTGGTGTTGATGGTGGTTTCCTGCCTTGTCACCCTCGAGCTGACggg                                                                <  1:354649/139‑1 (MQ=255)                                                                    cGTAAACGTGCTGCCTTCCACCGCGATGTTGATGATTATGTGTCTGAATTTGATGGGGGCAGGCGG‑CCCCGTCTGTTTGTCGCGGGTCt                                                                                                             >  1:254660/1‑89 (MQ=255)                                                                    cGTAAACGTGCTGCCTTCCACCGCGATGTTGATGATTATGTGTCTGAATTTGATGGGGGCAGGCGG‑CCCCGTCTGTTTGTCGCGGGTCt                                                                                                             <  2:254660/89‑1 (MQ=255)                                                                    cGTAAACGTGCTGCCTTCCACCGCGATGTTGATGATTATGTGTCTGAATTTGATGGGGGCAGGCGG‑CCCCGTCTGTTTGTCGCGGGTCTGGTGTTGATGGTGGTTTCCTGCCTTGTCACCCTCGAGCTg                                                                     >  1:226022/1‑129 (MQ=255)                                                                    cGTAAACGTGCTGCCTTCCACCGCGATGTTGATGATTATGTGTCTGAATTTGATGGGGGCAGGCGG‑CCCCGTCTGTTTGTCGCGGGTCTGGTGTTGATGGTGGTTTCCTGCCTTGTCACCCTCGAGCTg                                                                     <  2:226022/129‑1 (MQ=255)                                                                    cGTAAACGTGCTGCCTTCCACCGCGATGTTGATGATTATGTGTCTGAATTTGATGGGGGCAGGCGG‑CCCCGTCTGTTTGTCGCGGGTCTGGTGTTGATGGTGGCTTCCTGCCTTGTCACCCTCGAGCTGACGGGCATTa                                                           >  1:777067/1‑139 (MQ=255)                                                                         aCGTGCTGCCTTCCACCGCGATGTTGATGATTATGTGTCTGAATTTGATGGGGGCAGGCGG‑CCCCGTCTGTTTGTCGCGGGTCTGGTGTTGATGGTGGTTTCCTGCCTTGTCACCCTCGAGTTGACGGGCATTACCgtg                                                      <  2:693990/139‑1 (MQ=255)                                                                         aCGTGCTGCCTTCCACCGCGATGTTGATGATTATGTGTCTGAATTTGATGGGGGCAGGCGG‑CCCCGTCTGTTTGTCGCGGGTCTGGTGTTGATGGTGGTTTCCTGCCTTGTCACCCTCGAGCTGACGGGCATTACCgtg                                                      <  1:1096442/139‑1 (MQ=255)                                                                            tgctgcCTTCCACCGCGATGTTGATGATTATGTGTCTGAATTTGATGGGGGCAGGCGG‑CCCCGTCTGTTTGTCGCGGGTCTGGTGTTGATGGTGGTTTCCTGCCTTGTCACCCTCGAGCTGACGGGCATTACCGTGTCg                                                   >  1:586502/1‑139 (MQ=255)                                                                                 ccTTCCACCGCGATGTTGATGATTATGTGTCTGAATTTGATGGGGGCAggcggc                                                                                                                                    >  1:518867/1‑54 (MQ=255)                                                                                 ccTTCCACCGCGATGTTGATGATTATGTGTCTGAATTTGATGGGGGCAggcggc                                                                                                                                    <  2:518867/54‑1 (MQ=255)                                                                                 ccTTCCACCGCGATGTTGATGATTATGTGTCTGAATTTGATGGGGGCAGGCGG‑CCCCGTCTGTTTGTCGCGGGTCTGGTGTTGATGGTGGTTTCCTGCCTTGTCACCCTCGAGCTGACGGGCATTACCGTGTCGTTCaa                                              <  2:586502/139‑1 (MQ=255)                                                                                                ttGATGATTATGTGTCTGAATTTGATGGGGGCAGGCGG‑CCCCGTCTGTTTGTCGCGGGTCTGGTGTTGATGGTGGTTTCCTGCCTTGTCACCCTCGAGCTGACGGGCATTACCGTGTCGTTCAATAGTGCGCCGCTGGa                               >  1:445172/1‑139 (MQ=255)                                                                                                        tATGTGTCTGAATTTGATGGGGGCAGGCGG‑CCCCGTCTGTTTGTCGCGGGTCTGGTGTTGATGGTGGTTTCCTGCCTTGTCACCCTCGAGCTGACGGGCATTACCGTGTCGTTCAATAGTGCGCCGCTGGAATGGTGGc                       >  1:185395/1‑139 (MQ=255)                                                                                                        tATGTGTCTGAATTTGATGGGGGCAGGCGG‑CCCCGTCTGTTTGTCGCGGGTCTGGTGTTGATGGTGGTTTCCTGCCTTGTCACCCTCGAGCTGACGGGCATTACCGTGTCGTTCAATAGTGCGCCGCTGGAATGGTGGc                       <  2:445172/139‑1 (MQ=255)                                                                                                          tgtgTCTGAATTTGATGGGGGCAGGCGG‑CCCCGTCTGTTTGTcg                                                                                                                    >  1:570162/1‑44 (MQ=255)                                                                                                          tgtgTCTGAATTTGATGGGGGCAGGCGG‑CCCCGTCTGTTTGTcg                                                                                                                    <  2:570162/44‑1 (MQ=255)                                                                                                          tgtgTCTGAATTTGATGGGGGCAGGCGG‑CCCCGTCTGTTTGTCGCGGGTCTGGTGTTGATGGTGGTTTCCTGCCTTGTCACCCTCGAGCTGACGGGCATTACCGTGTCGTTCAATAGTGCGCCGCTGGAATGGTGGctc                     >  1:560680/1‑139 (MQ=255)                                                                                                          tgtgTCTGAATTTGATGGGGGCAGGCGG‑CCCCGTCTGTTTGTCGCGGGTCTGGTGTTGATGGTGGTTTCCTGCCTTGTCACCCTCGAGCTGACGGGCATTACCGTGTCGTTCAATAGTGCGCCGCTGGAATGGTGGctc                     >  1:510398/1‑139 (MQ=255)                                                                                                             gtCTGAATTTGATGGGGGCAGGCGG‑CCCCGTCTGTTTGTCGCGGGTCTGGTGTTGATGGTGGTTTCCTGCCTTGTCACCCTCGAGCTGACggg                                                                >  1:201242/1‑93 (MQ=255)                                                                                                             gtCTGAATTTGATGGGGGCAGGCGG‑CCCCGTCTGTTTGTCGCGGGTCTGGTGTTGATGGTGGTTTCCTGCCTTGTCACCCTCGAGCTGACggg                                                                <  2:201242/93‑1 (MQ=255)                                                                                                               cTGAATTTGATGGGGGCAGGCGG‑CCCCGTCTGTTTGTCGCGGGTCTGGTGTTGATGGTGGTTTCCTGCCTTGTCACCCTCGAGCTGACGGGCATTACCGTGTCGTTCAATAGTGCGCCGCTGGAATGGTGGCTCTCCCt                >  2:574788/1‑139 (MQ=255)                                                                                                                        aTGGGGGCAGGCGG‑CCCCGTCTGTTTGTCGCGGGTCTGGTGTTGATGGTGGTTTCCTGCCTTGTCACCCTCGAGCTGACGGGCATTACCGTGTCGTTCAATAGTGCGCCGCTGGAATGGTGGctc                     <  1:685865/125‑1 (MQ=255)                                                                                                                        aTGGGGGCAGGCGG‑CCCCGTCTGTTTGTCGCGGGTCTGGTGTTGATGGTGGTTTCCTGCCTTGTCACCCTCGAGCTGACGGGCATTACCGTGTCGTTCAATAGTGCGCCGCTGGAATGGTGGctc                     >  2:685865/1‑125 (MQ=255)                                                                                                                         tGGGGGCAGGCGG‑CCCCGTCTGTTTGTCGCGGGTCTGGTGTTGATGGTGGTTTCCTGCCTTGTCACCCTCGAGCTGACGGGCATTACCGTGTCGTTCAATAGTGCGCCGCTGGAATGGTGGCTCTCCCTTCCCattatt      <  2:777067/139‑1 (MQ=255)                                                                                                                           ggggCAGGCGG‑CCCCGTCTGTTTGTCGCGGGTCTGGTGTTGATGGTGGTTTCCTGCCTTGTCACCCTCGAGCTGACGGGCATTACCGTGTCGTTCAATAGTGCGCCGCTGGAATGGTGGCTCTCCCTTCCCATTATTGt    <  2:185395/139‑1 (MQ=255)                                                                                                                             ggCAGGCGG‑CCCCGTCTGTTTGTCGCGGGTCTGGTGTTGATGGTGGTTTCCTGCCTTGTCACCCTCGAGCTGACGGGCATTACCGTGTCGTTCAATAGTGCGCCGCTGGAATGGTGGCTCTCCCTTCCCATTATTGTCa  >  2:249737/1‑139 (MQ=255)                                                                                                                                      |                                                                                                                                   TTAGCCGGGAAATTTACAACTTAAAAACCGATGCAGTATTAGCGGGAATGTGGGTAGGCGTAATGGGCGTAAACGTGCTGCCTTCCACCGCGATGTTGATGATTATGTGTCTGAATTTGATGGGGGCAGGCGGCCCCCGTCTGTTTGTCGCGGGTCTGGTGTTGATGGTGGTTTCCTGCCTTGTCACCCTCGAGCTGACGGGCATTACCGTGTCGTTCAATAGTGCGCCGCTGGAATGGTGGCTCTCCCTTCCCATTATTGTCA  >  NC_000913/403986‑404249

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