Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,694,840	G→A	H133Y (CAT→TAT)	bcsA ←	cellulose synthase, catalytic subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,694,840	0	G	A	100.0%	131.5 / NA	41	H133Y (CAT→TAT)	bcsA	cellulose synthase, catalytic subunit
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (17/24);  total (17/24)

ATAAACTGCGCCAGCGGGTTAAACGGCTGAGTAACGCAGATTAACGCCAGAATCAGCGAGAAGGTGACGATGATACCGAGGATCAACCGCCGCGCCCCGGCACTCAAATGACCGAGCTCTTTTTTCTCATCCAGATGCTGTGTTTTATGGCTAACGCGCTCAGGCAGCTCGTTCATCCATTGATGGTAACGTCCACGAATATTTTGCAGACCTGAAAATGCCCGCCTGCGCGGTTTCGGCGTTTCTTTGCGCGATGCACC  >  NC_000913/3694704‑3694963                                                                                                                                         |                                                                                                                            aTAAACTGCGCCAGCGGGTTAAACGGCTGAGTAACGCAGATTAACGCCAGAATCAGCGAGAAGGTGACGATGATACCGAGGATCAACCGCCGCGCCCCGGCACTCAAATGACCGAGCTCTTTTTTCTCATCCAGATACt                                                                                                                           <  2:197455/139‑1 (MQ=255)   aaaCTGCGCCAGCGGGTTAAACGGCTGAGTAACGCAGATTAACGCCAGAATCAGCGAGAAGGTGACGATGATACCGAGGATCAACCGCCGCGCCCCGGCACTCAAATGACCGAGCTCTTTTTTCTCATCCAGATACtgt                                                                                                                         <  2:303070/139‑1 (MQ=255)              gCGGGTTAAACGGCTGAGTAACGCAGATTAACGCCAGAATCAGCGAGAAGGTGACGATGATACCGAGGATCAACCGCCGCGCCCCGGCACTCAAATGACCGAGCTCTTTTTTCTCATCCAGATACTGTGTTTTATGGCt                                                                                                              <  1:725013/139‑1 (MQ=255)                gggTTAAACGGCTGAGTAACGCAGATTAACGCCAGAATCAGCGAGAAGGTGACGATGATACCGAGGATCAACCGCCGCGCCCCGGCACTCAAATGACCGAGCTCTTTTTTCTCATCCAGATACTGTGTTTTATGGCTaa                                                                                                            <  2:805706/139‑1 (MQ=255)                   ttAAACGGCTGAGTAACGCAGATTAACGCCAGAATCAGCGAGAAGGTGACGATGATACCGAGGATCAACCGCCGCGCCCCGGCACTCAAATGACCGAGCTCTTTTTTCTCATCCAGATACTGTGTTTTATGGCTAAcgc                                                                                                         >  1:127754/1‑139 (MQ=255)                             gAGTAACGCAGTTTAACGCCAGAATCAGCGAGAAGGTGACGATGATACCGAGGATCAACCGCCGCGCCCCGGCACTCAAATGACCGAGCTCTTTTTTCTCATCCAGATACTCTGTTTTATGGCTAACGCGCTCAGGCAg                                                                                               <  1:543134/139‑1 (MQ=255)                                     cAGATTAACGCCAGAATCAGCGAGAAGGTGACGATGATACCGAGGATCAACCGCCGCGCCCCGGCACTCAAATGACCGAGCTCTTTTTTCTCATCCAGATACTGTGTTTTATGGCTAACGCGCTCAGGCAGCTCGTTCa                                                                                       >  1:219286/1‑139 (MQ=255)                                        aTTAACGCCAGAATCAGCGAGAAGGTGACGATGATACCGAGGATCAACCGCCGCGCCCCGGCACTCAAATGACCGAGCTCTTTTTTCTCATCCAGATACTGTGTTTTATGGCTAACGCGCTCAGGCAGCTCGTTCATcc                                                                                    >  2:688765/1‑139 (MQ=255)                                          tAACGCCAGAATCAGCGAGAAGGTGACGATGATACCGAGGATCAACCGCCGCGCCCCGGCACTCAAATGACCGAGCTCTTTTTTCTCATCCAGATACTGTGTTTTATGGCTAACGCGCTCAGGCAGCTCGTTCATCCAt                                                                                  <  2:475509/139‑1 (MQ=255)                                                   aaTCAGCGAGAAGGTGACGATGATACCGAGGATCAACCGCCGCGCCCCGGCACTCAAATGACCGAGCTCTTTTTTCTCATCCAGATACtg                                                                                                                          >  2:759971/1‑90 (MQ=255)                                                   aaTCAGCGAGAAGGTGACGATGATACCGAGGATCAACCGCCGCGCCCCGGCACTCAAATGACCGAGCTCTTTTTTCTCATCCAGATACtg                                                                                                                          <  1:759971/90‑1 (MQ=255)                                                                    cGATGATACCGAGGATCAACCGCCGCGCCCCGGCACTCAAATGACCGAGCTCTTTTTTCTCATCCAGATACTGTGTTTTATGGCTAACGCGCTCAGGCAGCTCGTTCATCCATTGATGGTAACGTCCACGAATATTTTg                                                        <  2:219286/139‑1 (MQ=255)                                                                             cGAGGATCAACCGCCGCGCCCCGGCACTCAAATGACCGAGCTCTTTTTTCTCATCCAGATACTGTGTTTTATGGCTAACGCGCTCAGGCAGCTCGTTCATCCATTGATGGTAACGTCCACGAATATTTTGCAGACCTGa                                               <  1:64600/139‑1 (MQ=255)                                                                               aGGATCAACCGCCGCGCCCCGGCACTCAAATGACCGAGCTCTTTTTTCTCATCCAGATACTGTGTTTTATGGCTAACGCGCTCAGGCAGCTCg                                                                                           >  1:122380/1‑93 (MQ=255)                                                                               aGGATCAACCGCCGCGCCCCGGCACTCAAATGACCGAGCTCTTTTTTCTCATCCAGATACTGTGTTTTATGGCTAACGCGCTCAGGCAGCTCg                                                                                           <  2:122380/93‑1 (MQ=255)                                                                               aGGATCAACCGCCGCGCCCCGGCACTCAAATGACCGAGCTCTTTTTTCTCATCCAGATACTGTGTTTTATGGCTAACGCGCTCAGGCAGCTCGTTCATCCATTGATGGTAACg                                                                       >  1:608915/1‑113 (MQ=255)                                                                               aGGATCAACCGCCGCGCCCCGGCACTCAAATGACCGAGCTCTTTTTTCTCATCCAGATACTGTGTTTTATGGCTAACGCGCTCAGGCAGCTCGTTCATCCATTGATGGTAACg                                                                       <  2:608915/113‑1 (MQ=255)                                                                                  aTCAACCGCCGCGCCCCGGCACTCAAATGACCGAGCTCTTTTTTCTCATCCAGATACTGTGTTTTATGGCTAACGCGCTCAGGCAGCTCGTTCATCCATTGATGGTAACGTCCACGAATATTTTGCAGACCTGAAAATg                                          <  1:790499/139‑1 (MQ=255)                                                                                   tCAACCGCCGCGCCCCGGCACTCAAATGACCGAGCTCTTTTTTCTCATCCAGATACTGTGTTTTATGGCTAACGCGCTCAGGCAGCTCGTTCATCCATTGATGGTAACGTCCACGAATATTTTGCAGACCTGa                                               <  1:604892/133‑1 (MQ=255)                                                                                   tCAACCGCCGCGCCCCGGCACTCAAATGACCGAGCTCTTTTTTCTCATCCAGATACTGTGTTTTATGGCTAACGCGCTCAGGCAGCTCGTTCATCCATTGATGGTAACGTCCACGAATATTTTGCAGACCTGa                                               >  2:604892/1‑133 (MQ=255)                                                                                    cAACCGCCGCGCCCCGGCACTCAAATGACCGAGCTCTTTTTTCTCATCCAGATACTGTGTTTTATGGCCAACGCGCTCAGGCAGCTCGTTCATCCAt                                                                                  >  1:829162/1‑97 (MQ=255)                                                                                    cAACCGCCGCGCCCCGGCACTCAAATGACCGAGCTCTTTTTTCTCATCCAGATACTGTGTTTTATGGCCAACGCGCTCAGGCAGCTCGTTCAGCCAt                                                                                  <  2:829162/97‑1 (MQ=255)                                                                                      aCCGCCGCGCCCCGGCACTCAAATGACCGAGCTCTTTTTTCTCATCCAGATACtgt                                                                                                                         >  1:668456/1‑56 (MQ=255)                                                                                      aCCGCCGCGCCCCGGCACTCAAATGACCGAGCTCTTTTTTCTCATCCAGATACtgt                                                                                                                         <  2:668456/56‑1 (MQ=255)                                                                                      aCCGCCGCGCCCCGGCACTCAAATGACCGAGCTCTTTTTTCTCATCCAGATACTGTGTTTTATGGCTAACGCGCTCAGGCAGCTCGTTCATCCATTGATGGTAACGTCCACGAATATTTTGCAGACCTGAAAATGCCCg                                      >  1:741037/1‑139 (MQ=255)                                                                                      aCCGCCGCGCCCCGGCACTCAAATGACCGAGCTCTTTTTTCTCATCCAGATACTGTGTTTTATGGCTAACGCGCTCAGGCAGCTCGTTCATCCATTGATGGTAACGTCCACGAATATTTTGCAGACCTGAAAATGCCCg                                      >  1:357596/1‑139 (MQ=255)                                                                                          ccgcGCCCCGGCACTCAAATGACCGAGCTCTTTTTTCTCATCCAGATACTGTGTTTTATGGCTAACGCGCTCAGGCAGCTCGTTCATCCATTGATGGTAACGTCCACGAATATTTTGCAGACCTGAAAATGCCCGCCTg                                  >  2:231872/1‑139 (MQ=255)                                                                                           cgcgCCCCGGCACTCAAATGACCGAGCTCTTTTTTCTCATCCAGATACTGTGTTTTATGGCTAACGCGCTCAGGCAGCTCGTTCATCCATTGATGGTAACGTCCACGAATATTTTGCAGACCTGAAAATGCCCGCCTgc                                 >  2:418856/1‑139 (MQ=255)                                                                                            gcgcCCCGGCACTCAAATGACCGAGCTCTTTTTTCTCATCCAGATACTGTGTTTTATGGCTaa                                                                                                            >  1:108710/1‑63 (MQ=255)                                                                                            gcgcCCCGGCACTCAAATGACCGAGCTCTTTTTTCTCATCCAGATACTGTGTTTTATGGCTaa                                                                                                            <  2:108710/63‑1 (MQ=255)                                                                                                  cGGCACTCAAATGACCGAGCTCTTTTTTCTCATCCAGATACTGTGTTTTATGGCTAACGCGCTCAGGCAGCTCGTTCATCCAt                                                                                  >  1:189738/1‑83 (MQ=255)                                                                                                  cGGCACTCAAATGACCGAGCTCTTTTTTCTCATCCAGATACTGTGTTTTATGGCTAACGCGCTCAGGCAGCTCGTTCATCCAt                                                                                  <  2:189738/83‑1 (MQ=255)                                                                                                    gCACTCAAATGACCGAGCTCTTTTTTCTCATCCAGATACTGTGTTTTATGGCTAACGCGCTCAGGCAGCTCGTTCATCCATTGATGGTAACGTCCAc                                                                  <  2:423982/97‑1 (MQ=255)                                                                                                    gCACTCAAATGACCGAGCTCTTTTTTCTCATCCAGATACTGTGTTTTATGGCTAACGCGCTCAGGCAGCTCGTTCATCCATTGATGGTAACGTCCAc                                                                  >  1:423982/1‑97 (MQ=255)                                                                                                               aCCGAGCTCTTTTTTCTCATCCAGATACTGTGTTTTATGGCTAACGCGCTCAGGCAg                                                                                               <  2:564100/57‑1 (MQ=255)                                                                                                               aCCGAGCTCTTTTTTCTCATCCAGATACTGTGTTTTATGGCTAACGCGCTCAGGCAg                                                                                               >  1:564100/1‑57 (MQ=255)                                                                                                               aCCGAGCTCTTTTTTCTCATCCAGATACTGTGTTTTATGGCTAACGCGCTCAGGCAGCTCGTTCATCCATTGATGGTAACGTCCACGAATATTTTGCAGACCTGAAAATGCCCGCCTGCGCGGTTTCGGCGTTTCTTTg             <  1:688765/139‑1 (MQ=255)                                                                                                                 cGAGCTCTTTTTTCTCATCCAGATACTGTGTTTTATGGCTAACGCGCTCAGGCAGCTCGTTCATCCATTGATGGTAACGTCCACGAATATTTTGCAGACCTGAAAATGCCCGCCTGCGCGGTTTCGGCGTTTCTTTgcg           <  2:127754/139‑1 (MQ=255)                                                                                                                 cGAGCTCTTTTTTCTCATCCAGATACTGTGTTTTATGGCTAACGCGCTCAGGCAGCTCGTTCATCCATTGATGGTAACGTCCACGAATATTTTGCAGACCTGAAAATGCCCGCCTGCGCGGTTTCGGCGTTTCTTTgcg           <  1:753710/139‑1 (MQ=255)                                                                                                                        ttttttCTCATCCAGATACTGTGTTTTATGGCTAACGCGCTCAGGCAGCTCGTTCATCCATTGATGGTAACGTCCACGAATATTTTGCAGACCTGAAAATGCCCGCCTGCGCGGTTTCGGCGTTTCTTTGCGCGATGCa    <  2:666029/139‑1 (MQ=255)                                                                                                                          ttttCTCATCCAGATACTGTGTTTTATGGCTAACGCGCTCAGGCAGCTCGTTCATCCATTGATGGTAACGTCCACGAATATTTTGCAGACCTGAAAATGCCCGCCTGCGCGGTTTCGGCGTTTCTTTGCGCGATGCAcc  <  1:231872/139‑1 (MQ=255)                                                                                                                                         |                                                                                                                            ATAAACTGCGCCAGCGGGTTAAACGGCTGAGTAACGCAGATTAACGCCAGAATCAGCGAGAAGGTGACGATGATACCGAGGATCAACCGCCGCGCCCCGGCACTCAAATGACCGAGCTCTTTTTTCTCATCCAGATGCTGTGTTTTATGGCTAACGCGCTCAGGCAGCTCGTTCATCCATTGATGGTAACGTCCACGAATATTTTGCAGACCTGAAAATGCCCGCCTGCGCGGTTTCGGCGTTTCTTTGCGCGATGCACC  >  NC_000913/3694704‑3694963

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