Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,693,321	T→C	Y639C (TAC→TGC)	bcsA ←	cellulose synthase catalytic subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,693,321	0	T	C	100.0%	92.7 / NA	27	Y639C (TAC→TGC)	bcsA	cellulose synthase catalytic subunit
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (13/14);  total (13/14)

AGAACAATCAGGTTGTAGAACACCCACACCATACTGACGACCACGGTGAGCATCTCGGTTGGCGGGCCATAGAAGTAGCGCCAGATGCCTACCGCAACGCCCACCAGGTTGAGCAGGACAAGGAAGATGTAGGGCCGCGAGATCACCCAGTCGACGTACTCTTCTTCCACCAGTCCACCTTTGGCGGTGACGTTAAATTTGCCTTTGTGCGGGTTAATCAGCGCCACCAGCGTCGGTGGTGCGA  >  NC_000913/3693192‑3693435                                                                                                                                  |                                                                                                                   aGAACAATCAGGTTGGAGCACACCCACACCTTACTGACGACCACGGTGAGCATCTCGGTTGGCGGGCCATAGAAGTAGCGCCAGATGCCTACCGCAACGCCCACCAGGTTGAGCAGGACAAGGAAGATGCAgg                                                                                                                 <  1:138557/133‑1 (MQ=255)    aCAATCAGTTTGTAGCACACCCCCACCATACTGACGACCACGGTGAGCATCTCGGTTGGCGGGCCATAGAAG‑AGCGCCAGATGCCTACCGCAACGCCCACCAGGTTGAGCAGGACAAGGAAGATGCAGGGCcg                                                                                                             <  1:250538/133‑1 (MQ=255)           ggTTGTAGAAGACCCCCACCATACTGACGACCGCGGTGAGCATCTCGGTTGGCGGGCCATAGAAGTAGCGCCAGATGCCTACCGCAACGCCCACCAGGTTGAGCAGGACAAGGAAGATGCAGGGCCGCGAGAt                                                                                                       <  2:220197/133‑1 (MQ=255)             ttGTAGAACACCCACACCATACTGACGACCACGGTGAGCATCTCGGTTGGCGGGCCATAGAAGTAGCGCCAGATGCTTACAGCAACGCCCACCAGGTTGAGCAGGACAAGGAAGATGCAGGGCCGCGAGACCa                                                                                                     >  2:85251/1‑133 (MQ=255)             ttGTAGAACACCCACACCATACTGACGACCACGGTGAGCATCTCGGTTGGCGGGCCATAGAAGTAGCGCCAGATGCCTACCGCAACGCCCACCAGGTTGAGCAGGACAAGGAAGATGCAGGGCCGCGAGATCa                                                                                                     >  2:34613/1‑133 (MQ=255)                      acCCACACCATACTGACGACCACGGTGAGCATCTCGGTTGGCGGGCCATTGAAGTACCGCCAGATACCTACCGCAACGCCCACGAGGTTGCAGACGACGAGGAAGAGCAAGGGCCGCAAGATCACCCAGTcaa                                                                                            >  2:230637/1‑131 (MQ=255)                          acacCATACTGACGACCACGGTGAGCATCTCGGTTGGCGGGCCATAGAAGTAGCGCCAGATGCCTACCGCAACGCCCACCAGGTTGAGCAGGACAAGGAAGATGCAGGGCCGCGAGATCACCCAGTCGACGTa                                                                                        <  1:88840/133‑1 (MQ=255)                            acCATACTGATGACCTCGGTGAGCATCTCGGTTGGCGGGCCCTAGAAGTAGCGCCAGATGCCTACCGCAACGCCCACCAGGTTGAGCAGGACAAGGAAGATGCAGGGCCGCGAGATCACCCAGTCGACGTAct                                                                                      <  2:207063/133‑1 (MQ=255)                            acCATACTGACGACCACGGTGAGCATCTCGGTTGGCGGGCCATAGAAGTAGCGCCAGATGCCTACCGCAACGCCCACCAGGTTGAGCAGGACAAGGAAGATGCAGGGCCGCGAGATCACCCAGTCGACGTAct                                                                                      <  1:70368/133‑1 (MQ=255)                                tACTGACGACCACGGTGAGCATCTCGGTTGGCGGGCCATAGAAGTAGCGCCAGATGCCTACCGCAACGCCCACCAGGTTGAGCAGGACAAGGAAGATGCAGGGCCGCGAGATCACCCAGTCGACGTACtcttc                                                                                  >  1:265779/1‑133 (MQ=255)                                    gacgacCACGGTGAGCATCTCGGTTGGCGGGCCATAGAAGTAGCGCCAGATGCCTACCGCAACGCCCACCAGGTTGAGCAGGACAAGGAAGATGCAGGGCCGCGAGATCACCCAGTCGACGTACTCTTCTTcc                                                                              >  1:36753/1‑133 (MQ=255)                                     acgacCACGGTGAGCATCTCGGTTGGCGGGCCATAGAAGTAGCGCCAGATGCCTACCGCAACGCCCACCAGGTTGAGCAGGACAAGGAAGATGCAGGGCCGCGAGATCACCCAGTCGACGTACTCTTCTTcca                                                                             >  2:120144/1‑133 (MQ=255)                                     acgacCACGGTGAGCATCTCGGTTGGCGGGCCATAGAAGTAGCGCCAGATGCCTACCGCAACGCCCACCAGGTTGAGCAGGACAAGGAAGATGCAGGGCCGCGAGATCACCCAGTCGACGTACTCTTCTTcca                                                                             >  1:286945/1‑133 (MQ=255)                                        acCACGGTGAGCATCTCGGTTGGCGGGCCATAGAAGTAGCGCCAGATGCCTACCGCAACGCCCACCAGGTTGAGCAGGACAAGGAAGATGCAGGGCCGCGAGATCACCCAGTCGACGTACTCTTCTTccacca                                                                          >  2:434992/1‑133 (MQ=255)                                        acCACGGTGAGCATCTCGGTTGGCGGGCCATAGAAGTAGCGCCAGATGCCTACCGCAACGCCCACCAGGTTGAGCAGGACAAGGAAGATGCAGGGCCGCGAGATCACCCAGTCGACGTACTCTTCTTccacca                                                                          >  1:312353/1‑133 (MQ=255)                                                 aGCATCTCGGTTGGCGGGCCATAGAAGTAGCGCCAGATGCCTACCGCAACGCCCACCAGGTTGAGCAGGACAAGGAAGATGCAGGGCCGCGAGATCACCCAGTCGACGTACTCTTCTTCCACCAGTCCACCtt                                                                 >  2:273887/1‑133 (MQ=255)                                                       tcGGTTGGCGGGCCATAGAAGTAGCGGCAGATGCCTGCCGCAACGCCCACCTGGTTGAGCAAGACGATGAAGATGCAGGGCCGCGAGATCACCCAGTCGACGTACTCTTCTTCCACCAGTCCACCTTTGGCgg                                                           <  2:12762/133‑1 (MQ=255)                                                             ggCGGGCCATAGAAGTAGCGCCAGATGCCTACCGCAACGCCCACCAGGGTGAGCAGGACAAGGAAGATGCAGGGCCGCGAGATCACCCAGTCGACGTACTATTCATCCGCCAGTCCAACTTTGGCGGTGACGt                                                     >  2:192923/1‑133 (MQ=255)                                                                        gAAGTAGCTCCAGATTCCCACCTCAACGCCAACCAGGTTGAGCAGGACAAGGAAGATGCAGGGCCGCGAGATCACCCAGTCGACGTACTCTTCTTCCACCAGTCCACCTTTGGCGGTGACGTTAAATTTGCCt                                          <  1:192923/133‑1 (MQ=255)                                                                            tAGCGCCAGATGCCTACCGCAACGCCCACCAGGTTGAGCAGGACAAGGAAGATGCAGGGCCGCGAGATCACCCAGTCGACGTACTCTTCTTCCACCAGTCCACCTTTGGCGGTGACGTTAAATTTGCCTTtgt                                      <  1:483060/133‑1 (MQ=255)                                                                            tAGCGCCAGATGCCTACCGCAACGCCAACCAGGTTGAGCAGGACAAGGAAGATGCAGGGCCGCGAGATCACCCAGTCGACGTACTCTTCTTCCACCAGTCTACCATTGGCGGTGATGGTAAACTTGCCTTtgt                                      >  2:179075/1‑133 (MQ=255)                                                                                 ccAGATGCCTACCGCAACGCCCACCAGGTTGAGCAGGACAAGGAAGATGCAGGGCCGCGAGATCACCCAGTCGACGTACTCTTCTTCCACCAGTCCACCTTTGGAGGTGACGTTAAACTTGCCTTTGTGCggg                                 >  2:188679/1‑133 (MQ=255)                                                                                   aGATGCCTACCGCGACGCCCACCAGGTGCGGCAGGACAAGGAAGATGCAGTTCCGCGAGACCACCCAGTCGACGTACTCTTCTTCCACCAGTCCCCCTTTGGCGGAGCCGGTACATTTGCCTTTGTGCGGGtt                               >  2:362788/1‑133 (MQ=255)                                                                                         cTACCGCAACGCCCACCAGGTTGAGCAGGACAAGGAAGATGCAGGGCCGCGAGATCACCCAGTCGACGTACTCTTCTTCCACCAGTCCACCTTTGGCGGTGACGTTAAATTTGCCTTTGTGCGGGTTAATCAg                         <  1:34613/133‑1 (MQ=255)                                                                                              gCAACGCCCACCAGGTTGAGCAGGACAAGGAAGATGCAGGGCCGCGAGATCACCCAGTCGACGTACTCTTCTTCCACCAGTCCACCTTTGGCGGTGACGTTAAATTTGCCTTTGTGCGGGTTAATCAGCGcca                    <  2:36753/133‑1 (MQ=255)                                                                                                aaCGCCCACCAGGTTGAGCAGGACAAGGAAGATGCAGGGCCGCGAGATCACCCAGTCGACGTACTCTTCTTCCACCAGTCCACCTTTGGCGGTGACGTTAAATTTGCCTTTGTGCGGGTTAATCAGCGccacc                  >  2:505151/1‑133 (MQ=255)                                                                                                          aGGTTGAGCAGGACAAGGAAGATGCAGGGCCGCGAGATCACCCAGTCGACGTACTCTTCTTCCACCAGTCCACCTTTGGCGGTGACGTTAAATTTGCCTTTGTGCGGGTTAATCAGCGCCACCAGCGTCggtg        <  1:434992/133‑1 (MQ=255)                                                                                                          aGGTTGAGCAGGACAAGGAAGATGCAGGGCCGCGAGATCACCCAGTCGACGTACTCTTCTTCCACCAGTCCACCTTTGGCGGTGACGTTAAATTTGCCTTTGTGCGGGTTAATCAGCGCCACCAGCGTCggtg        <  2:265779/133‑1 (MQ=255)                                                                                                          aGGTTGAGCAGGACAAGGAAGATGCAGGGCCGCGAGATCACCCAGTCGACGTACTCTTCTTCCACCAGTCCACCTTTGGCGGTGACGTTAAATTTGCCTTTGTGCGGGTTAATCAGCGCCACCAGCGTCggtg        <  1:505151/133‑1 (MQ=255)                                                                                                           ggTTGAGCAGGACAAGGAAGATGCAGGGCCGCGAGATCACCCAGTCGACGTACTCTTCTTGCACCAGTCCACCTTTGGCGGTGACGTTAAATTTGCCGTTGTGCGGGTTAATCAGCGCCACCAGCGTCggtgg       <  1:98430/133‑1 (MQ=255)                                                                                                                aGCAGGACAAGGAAGATGCAGGGCCGCGAGATCACCCAGTCGACGTACTCTTCTTCCACCAGTCCACCTTTGGCGGTGACGTTAAATTTGCCTTTGTGCGGGTTAATCAGCGCCACCAGCGTCGGTGGTGCGa  >  1:436805/1‑133 (MQ=255)                                                                                                                                  |                                                                                                                   AGAACAATCAGGTTGTAGAACACCCACACCATACTGACGACCACGGTGAGCATCTCGGTTGGCGGGCCATAGAAGTAGCGCCAGATGCCTACCGCAACGCCCACCAGGTTGAGCAGGACAAGGAAGATGTAGGGCCGCGAGATCACCCAGTCGACGTACTCTTCTTCCACCAGTCCACCTTTGGCGGTGACGTTAAATTTGCCTTTGTGCGGGTTAATCAGCGCCACCAGCGTCGGTGGTGCGA  >  NC_000913/3693192‑3693435

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