Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,692,114	C→T	S165N (AGC→AAC)	bcsB ←	cellulose synthase periplasmic subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,692,114	0	C	T	100.0%	59.5 / NA	22	S165N (AGC→AAC)	bcsB	cellulose synthase periplasmic subunit
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (8/14);  total (8/14)

CAGCCCAACATCCGGCGCACCCGCAAAGACCATCGGCAAGGTGTTGGTGCGGTTATCGCTCGGGTCAAAGAATGGCACCGGGAAGTGTGACAGGTCATTCTTCACATTCAGGGTCTGATAGGTCAGATCCAGTCCACTGCTCCGCCCAACATCCAGCCAAAGCGTGGTGCTGGCCGGTTTTTCGCACACGTCCTGATAATGGCCGACAAACTCCAGCCGTACACGGTTGAAGTCGCTAATAAACAGTGGGTTAATGGGCATTTGCGCCAGCGTTTTTTTACCCAACTGTTCTTTGGTCACTGGCAGCACGCCCATCAGTTCATCATTGAGATAAACCTTTAACTGCGACTGGACAGGCAGTAACGATGGCGATGGGGTGTATTCGAGGTTGAGCATCGCCTTCGTCACCACTTCATCGCTGCGCATACCAAACTCAATGCTGCCGTTCGGGTTAATGCCACGT  >  NC_000913/3691879‑3692341                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                    cAGCCCAACATCCGGCGCACCCGCAAAGACCATCGGCAAGGTGTTGGTGAGGTTATCGCTCGGGTCAACGCATGGCACCGGGAAGCGTGACAGGTCATTCTTCACATTCAGGGTCTGCTAGGTCAGATCCAGTCCACTGCTCCGCCCAACATCCAGCCAAAGCGTGGTGCTGGCCGGTTTTTCGCACACGTCCTGATAATGGCCGACAAACTCCAGCCGTACACGGTTGAAGTCGTTAATAAACAGTGGGt                                                                                                                                                                                                                      >  1:215439/1‑251 (MQ=255)                                                 gCGGTTATCGCGCGGGTCAAAGAATGGCACCGGGAAGTGTGACAGGTCATTCTTCACATTCAGGGTCTGATAGGTCAGATCCAGTCCACTGCTCCGCCCAACATCCAGCCAAAGCGTGGTGCTGGCCGGTTTGTCGCACACGTCCTGATAATGGCCGACAAACTCCAGCCGTGCACGGTTGAAGTCGTTAATAAACAGTGGGTTAATGGGCATTTGCGCCAGCGTTTGTTTACCCAACTGTTCTTTGGTc                                                                                                                                                                       <  1:354888/250‑1 (MQ=255)                                                   ggTTATCGCTCGGGTCAAAGAATGGCACCGGGAAGTGTGACAGGTCATTCTTCACATTCAGGGTCTGATAGGTCAGATCCAGTCCACTGCTCCGCCCAACATCCAGCCAAAGCGTGGTGCTGGCCGGTTTTTCGCACACGTCCTGATAATGGCCGACAAACTCCAGCCGTACACGGTTGAAGTCGTTAATAAACAGTGGGTTAATGGGCATTTGCGCCAGCGTTTTTTTACCCAACTGTTCTTTGGTc                                                                                                                                                                       >  1:165064/1‑248 (MQ=255)                                                                gTCAAAGAATGGCACCGGGAAGTGTGACAGGTCATTCTTCACATTCAGGGTCTGATAGGTCAGATCCAGTCCACTGCTCCGCCCAACATCCAGCCAAAGCGTGGTGCTGGCCGGTTTTTCGCACACGTCCTGATAATGGCCGACAAACTCCAGCCGTACACGGTTGAAGTCGTTAATAAACAGTGGGTTAATGGGCATTTGCGCCAGCGTTTTTTTAcccn                                                                                                                                                                                     <  1:268313/221‑2 (MQ=255)                                                                  cAAAGAATGGCACCGGGAAGTGTGACAGGTCATTCTTCACATTCAGGGTCTGATAGGTCAGATCCAGTCCACTGCTCCGCCCAACATCCAGCCAAAGCGTGGTGCTGGCCGGTTTTTCGCACACGTCCTGATAATGGCCGACAAACTCCAGCCGTACACGGTTGAAGTCGTTAATAAACAGTGGGTTAATGGGCATTTGCGCCAGCGTTTTTTTACCCAACTGTTCTTTGGTCACTGGCAGCACGCCCATc                                                                                                                                                     >  1:206025/1‑251 (MQ=255)                                                                        aTGGCACCGGGAAGTGTGACAGGTCATTCTTCACATTCAGGGTCTGATAGGTCAGATCCAGTCCACTGCTCCGCCCAACATCCAGCCAAAGCGTGGTGCTGGCCGGTTTTTCGCACACGTCCTGATAATGGCCGACAAACTCCAGCCGTACACGGTTGAAGTCGTTAATAAACAGTGGGTTAATGGGCATTTgcgc                                                                                                                                                                                                      <  1:157512/196‑1 (MQ=255)                                                                                                                ggTCTGATAGGTCAGATCCAGTCCACTGCTCCGCCCAACATCCAGCCAAAGCGTGGTGCTGGCCGGTTTTTCGCACACGTCCTGATAATGGCCGACAAACTCCAGCCGTACACGGTTGAAGTCGTTAATAAACAGTGGGTTAATGGGCATTTGCGCCAGCGTTTTTTTaac                                                                                                                                                                                       <  1:417996/171‑3 (MQ=255)                                                                                                                    tGATAGGTCAGATCCAGTCCACTGCTCCCCCCAACATCCAGCCAAAGCGAGGAGATGACCGGTTTTTCGCCCCCGGCCTGAAAATGGCCGACAAACCCCAGCCGTACACGGTTGAACTCGTTAATAAACCGTGGGTTAATGGGCATTTGCGCCAGCGTTTTTTTACACAACTGTCCTTTGGACACTGGCAGCACGCCCATCAGTTCATCCTTGAGATAAACCTTTAACTGCGACTGGACAGGCAGTAAcg                                                                                                    >  1:229018/1‑250 (MQ=255)                                                                                                                         ggTCAGATCCAGTTCACTGCTCCGCCCAACATCCAGCCAAAGCGTGGTGCTGGCCGGTTTTTCGCACACGTCCTGATAATGGCCGACAAACTCCAGCCGTACACGGTTGAAGTCGTTAATAAACAGTGGGTTAATGGGCATTTGCGCCAGCGTTTTTTTACCCAACTGTTCtttn                                                                                                                                                                          <  1:129243/175‑2 (MQ=255)                                                                                                                          gTCAGATCCAGTCCACTGCTCCGCCCAACATCCAGCCAAAGCGTGGTGCTGTCCGGTTTTTCGCACACGTCCTGATAATGGCCGACAAACGCCAGCCGTGCACGGTGGAAGTCGTTAATAAACAGTGGGTTAATGGGCATTTGCGCCAGCGTTTGTTTACCCAACTg                                                                                                                                                                                 <  1:64059/167‑1 (MQ=255)                                                                                                                              gATCCAGTCCACTGCTCCGCCCAACATCCAGCCAAAGCGTGGTGCTGGCCGGTTTTTCGCACACGTCCTGATAATGGCCGACAAACTCCAGCCGTACACGGTTGAAGTCGTTAATAAACAGTGGGTTAATGGGCATTTgc                                                                                                                                                                                                        <  1:193937/140‑1 (MQ=255)                                                                                                                                         cTGCTCCGCCCAACATCCAGCCAAAGCGTGGTGCTGGCCGGTTTTTCGCACACGTCCTGATAATGGCCGACAAACTCCAGCCGTACACGGTTGCAGTCGTTAATAAACAGTGGGTTAATGGGCATTTGCGCCAGCGttttt                                                                                                                                                                                            >  1:227947/1‑141 (MQ=255)                                                                                                                                                 cccAACATCCAGCCAAATCGTTGTGCTGGCCGGTTTTTCGCACACGTCCTGATAAGGGCCGACAAACTCCAGCCGTACACGGTTGAAGTCGTTAATAAACAGTGGGTTAATGGGCATTTgcgcn                                                                                                                                                                                                     <  1:87317/124‑2 (MQ=255)                                                                                                                                                  ccAACATCCAGCCAAAGCGTAGTGCTGGCCGGTTTTTCGCACACGTCCTGATAATGGCCGACAAACTCCAGCCGTACACGGTTGAAGTCGTTAATAAACAGTGGGTTAATGGGCATTTGCGCCAGCGTTTTTTTACCCAACTGTTCTTTGGTCACTGGCAGCACGCCCa                                                                                                                                                       <  1:443834/169‑1 (MQ=255)                                                                                                                                                   cAACATCCAGCCAAAGCGTGGTGCTGGCCGGTTTTTCGCACACGTCCTGATAATGGCCGACAAACTACAGCAGTACACGGTTGAAGTCGTTAATAAACAGTGGGTTAATGGGCATTTGCGCCAGCGTTTTTTTACCCACCTGTTCTTTGGTCACTGGCAGCACGCCCATCAGTTCCTCATTGAGATAAACCTTTAACTGCGACTGGACAGGCAGTAACGATGGCGGTGGGGTGTATTCGAGGTTGAGCATc                                                                    >  1:312942/1‑251 (MQ=255)                                                                                                                                                         ncAGCCAAAGCGTGGTGCTGGCCGGTTTTTCGCACACGTCCTGATAATGGCCGACAAACTCCAGCCGTACACGGTTGAAGTCGTTAATAAACAGTGGGTTAATGGGCATTTgcg                                                                                                                                                                                                       >  1:269074/2‑114 (MQ=255)                                                                                                                                                         ccAGCCAAAGCGTGGTGCTGGCCGGTTTTTCGCACACGTCCTGATAATGGCCGACAAACTCCAGCCGTACACGGTTGAAGTCGTTAATAAACAGTGGGTTAATGGGCATTTgcg                                                                                                                                                                                                       <  2:269074/114‑1 (MQ=255)                                                                                                                                                                       ntgCTGGCCGGTTTTTCGCACACGTCCTGATAATGGCCGACAAACTCCAGCCGTACACGGTTGAAGTCGTTAATAAACAGTGGGTTAATGGGCCTTTGCGCCAGCGTTTTTTTACCCAACTGTTCTTTGGTCACTGGCAGCACGCCCATCAGTTCATCATTGAGATAAACCTTTAACTGCGACTGGACAGGCAGTAACGATGGCGATGGGGTGTATTCGAGGTTGAGCATCGCCTTCGTCACCACTTCATc                                                >  1:39575/2‑251 (MQ=255)                                                                                                                                                                                               tCCTGATAATGGCCGACAAACTCCAGCCGTACACGGTTGAAGTCGTTAATAAACAGTGGGTTAATGGGCATTTGCGCCAGCGTTTTTTTACCCAACTGTTCTTTGGTc                                                                                                                                                                       <  1:377416/108‑1 (MQ=255)                                                                                                                                                                                                   gATAATGGCCGACAAACTCCAGCCGTACACGGTTGAAGTCGTTAATAAACAGTGGGTTAATGGGCATTTGCGCCAGCGTTTTTTTACCCAACTGTTCTTTGGTCACTGGCAGCACGCCCATCAGTTCATCATTGAGATAAACCTTTAACTGCGNc                                                                                                                    <  1:348191/155‑1 (MQ=255)                                                                                                                                                                                                                   cTCCAGCCGTACACGGTTGAAGTCGTTAATAAACAGTGGGTTAATGGGCATTTGCGCCAGCGTTTTTTTACCCAACTGTTCTTTGGGCACTGGCATCACGCCCATCAGTTCATCATTGAGATAAAcc                                                                                                                                <  1:285856/127‑1 (MQ=255)                                                                                                                                                                                                                     ccAGCCGTACACGGTTGAAGTCGTTAATAAACAGTGGGTTAATGGGCATTTGCGCCAGCGTTTTTTTACCCAACTGTTCTTTGGTCACTGGCAGCACGCCCATCAGTTCATCATTGAGATAAACCTTTAACTGCGACTGGACAGGCAGTAACGATGGCGATGGGGTGTATTCGAGGTTGATCATCGCCTTCGTCACCACTTCATCGCTGCGCATACCAAACTCAATGCTGCCGTTCGGGTTAATGCCACGt  <  1:429274/251‑1 (MQ=255)                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                    CAGCCCAACATCCGGCGCACCCGCAAAGACCATCGGCAAGGTGTTGGTGCGGTTATCGCTCGGGTCAAAGAATGGCACCGGGAAGTGTGACAGGTCATTCTTCACATTCAGGGTCTGATAGGTCAGATCCAGTCCACTGCTCCGCCCAACATCCAGCCAAAGCGTGGTGCTGGCCGGTTTTTCGCACACGTCCTGATAATGGCCGACAAACTCCAGCCGTACACGGTTGAAGTCGCTAATAAACAGTGGGTTAATGGGCATTTGCGCCAGCGTTTTTTTACCCAACTGTTCTTTGGTCACTGGCAGCACGCCCATCAGTTCATCATTGAGATAAACCTTTAACTGCGACTGGACAGGCAGTAACGATGGCGATGGGGTGTATTCGAGGTTGAGCATCGCCTTCGTCACCACTTCATCGCTGCGCATACCAAACTCAATGCTGCCGTTCGGGTTAATGCCACGT  >  NC_000913/3691879‑3692341

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