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%	170.0 / NA	53	H133Y (CAT→TAT)	bcsA	cellulose synthase, catalytic subunit
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (22/31);  total (22/31)

ATAAACTGCGCCAGCGGGTTAAACGGCTGAGTAACGCAGATTAACGCCAGAATCAGCGAGAAGGTGACGATGATACCGAGGATCAACCGCCGCGCCCCGGCACTCAAATGACCGAGCTCTTTTTTCTCATCCAGATGCTGTGTTTTATGGCTAACGCGCTCAGGCAGCTCGTTCATCCATTGATGGTAACGTCCACGAATATTTTGCAGACCTGAAAATGCCCGCCTGCGCGGTTTCGGCGTTTCTTTGCGCGATGCACCGATCAATAAC  >  NC_000913/3694704‑3694973                                                                                                                                         |                                                                                                                                      aTAAACTGCGCCAGCGGGTTAAACGGCTGAGTAACGCAGATTAACGCCAGAATCAGCGAGAAGGTGACGATGATACCGAGGATCAACCGCCGCGCCCCGGCACTCAAATGACCGAGCTCTTTTTTCTCATCCAGATACt                                                                                                                                     >  1:302316/1‑139 (MQ=255)        gcgcCAGCGGGTTAAACGGCTGAGTAACGCAGATTAACGCCAGAATCAGCGAGAAGGTGACGATGATACCGAGGATCAACCGCCGCGCCCCGGCACTCAAATGACCGAGCTCTTTTTTCTCATCCAGATACTGTGtttt                                                                                                                              >  2:1206008/1‑139 (MQ=255)          gcCAGCGGGTTAAACGGCTGAGTAACGCAGATTAACGCCAGAATCAGCGAGAAGGTGACGATGATACCGAGGATCAACCGCCGCGCCCCGGCACTCAAATGACCGAGCTCTTTTTTCTCATCCAGATACTGTGTTTTAt                                                                                                                            >  1:535810/1‑139 (MQ=255)                gggTTAAACGGCTGAGTAACGCAGATTAACGCCAGAATCAGCGAGAAGGTGACGATGATACCGAGGATCAACCGCCGCGCCCCGGCACTCAAATGACCGAGCTCTTTTTTCTCATCCAGATACTGTGTTTTATGGCTaa                                                                                                                      <  2:535810/139‑1 (MQ=255)                             gAGTAACGCAGATTAACGCCAGAATCAGCGAGAAGGTGACGATGATACCGAGGATCAACCGCCGCGCCCCGGCACTCAAATGACCGAGCTCTTTTTTCTCATCCAGATACTGTGTTTTATGGCTAACGCGCTCAGGCAg                                                                                                         <  1:241739/139‑1 (MQ=255)                              aGTAACGCAGATTAACGCCAGAATCAGCGAGAAGGTGACGATGATACCGAGGATCAACCGCCGCGCCCCGGCACTCAAATGACCGAGCTCTTTTTTCTCATCCAGATACTGTGTTTTATGGCTAACGCGCTCAGGCAGc                                                                                                        <  2:806640/139‑1 (MQ=255)                               gTAACGCAGATTAACGCCAGAATCAGCGAGAAGGTGACGATGATACCGAGGATCAACCGCCGCGCCCCGGCACTCAAATGACCGAGCTCTTTTTTCTCATCCAGATACTGTGTTTTATGGCTAACGCGCTCAGGCAGCt                                                                                                       <  2:325847/139‑1 (MQ=255)                                      aGATTAACGCCAGAATCAGCGAGAAGGTGACGATGATACCGAGGATCAACCGCCGCGCCCCGGCACTCAAATGACCGAGCTCTTTTTTCTCATCCAGATACTGTGTTTTATGGCTAACGCGCTCAGGCAGCTCGTTCAt                                                                                                <  1:1206008/139‑1 (MQ=255)                                      aGATTAACGCCAGAATCAGCGAGAAGGTGACGATGATACCGAGGATCAACCGCCGCGCCCCGGCACTCAAATGACCGAGCTCTTTTTTCTCATCCAGATACTGTGTTTTATGGCTAACGCGCTCAGGCAGCTCGTTCAt                                                                                                <  2:5047/139‑1 (MQ=255)                                        aTTAACGCCAGAATCAGCGAGAAGGTGACGATGATACCGAGGATCAACCGCCGCGCCCCGGCACTCAAATGACCGAGCTCTTTTTTCTCATCCAGATACTGTGTTTTATGGCTAAcg                                                                                                                    <  2:569690/117‑1 (MQ=255)                                        aTTAACGCCAGAATCAGCGAGAAGGTGACGATGATACCGAGGATCAACCGCCGCGCCCCGGCACTCAAATGACCGAGCTCTTTTTTCTCATCCAGATACTGTGTTTTATGGCTAAcg                                                                                                                    >  1:569690/1‑117 (MQ=255)                                        aTTAACGCCAGAATCAGCGAGAAGGTGACGATGATACCGAGGATCAACCGCCGCGCCCCGGCACTCAAATGACCGAGCTCTTTTTTCTCATCCAGATACTGTGTTTTATGGCTAACGCGCTCAGGCAGCTCGTTCATcc                                                                                              >  2:661754/1‑139 (MQ=255)                                          tAACGCCAGAATCAGCGAGAAGGTGACGATGATACCGAGGATCAACCGCCGCGCCCCGGCACTCAAATGACCGAGCTCTTTTTTCTCATCCAGATACTGTGTTTTATGGCTAACGCGCTCAGGCAGCTCGTTCATCCAt                                                                                            <  1:225981/139‑1 (MQ=255)                                                  gAATCAGCGAGAAGGTGACGATGATACCGAGGATCAACCGCCGCGCCCCGGCACTCAAATGACCGAGCTCTTTTTTCTCATCCAGATACTGTGTTTTATGGCTAACGCGCTCAGGCAGCTCGTTCATCCATTGacgtta                                                                                    >  1:133483/1‑134 (MQ=255)                                                    aTCAGCGAGAAGGTGACGATGATACCGAGGATCAACCGCCGCGCCCCGGCACTCAAATGACCGAGCTCTTTTTTCTCATCCAGATACTGTGTTTTATGGCTAACGCGCTCAGGCAGCTCGTTCATCCATTGacgttacc                                                                                  <  2:133483/139‑8 (MQ=255)                                                    aTCAGCGAGAAGGTGACGATGATACCGAGGATCAACCGCCGCGCCCCGGCACTCAAATGACCGAGCTCTTTTTTCTCATCCAGATACTGTGTTTTATGGCTAACGCGCTCAGGCAGCTCGTTCATCCATTGATGGTAAc                                                                                  >  2:2964/1‑139 (MQ=255)                                                     tCAGCGAGAAGGTGACGATGATACCGAGGATCAACCGCCGCGCCCCGGCACTCAAATGACCGAGCTCTTTTTTCTCATCCAGATACTGTGTTTTATGGCTAACGCGCTCAGGCAGCTCGTTCATCCATTGATGGTAACg                                                                                 <  2:422691/139‑1 (MQ=255)                                                     tCAGCGAGAAGGTGACGATGATACCGAGGATCAACCGCCGCGCCCCGGCACTCAAATGACCGAGCTCTTTTTTCTCATCCAGATACTGTGTTTTATGGCTAACGCGCTCAGGCAGCTCGTTCATCCATTGATGGTAACg                                                                                 <  2:1352290/139‑1 (MQ=255)                                                       aGCGAGAAGGTGACGATGATACCGAGGATCAACCGCCGCGCCCCGGCACTCAAATGACCGAGCTCTTTTTTCTCATCCAGATACTGTGTTTTATGGCTAACGCGCTCAGGCAGCTCGTTCATCCATTGATGGTa                                                                                    <  1:206948/134‑1 (MQ=255)                                                       aGCGAGAAGGTGACGATGATACCGAGGATCAACCGCCGCGCCCCGGCACTCAAATGACCGAGCTCTTTTTTCTCATCCAGATACTGTGTTTTATGGCTAACGCGCTCAGGCAGCTCGTTCATCCATTGATGGTa                                                                                    >  2:206948/1‑134 (MQ=255)                                                       aGCGAGAAGGTGACGATGATACCGAGGATCAACCGCCGCGCCCCGGCACTCAAATGACCGAGCTCTTTTTTCTCATCCAGATACTGTGTTTTATGGCTAACGCGCTCAGGCAGCTCGTTCATCCATTGATGGTAACGTc                                                                               <  1:524376/139‑1 (MQ=255)                                                             aaGGTGACGATGATACCGAGGATCAACCGCCGCGCCCCGGCACTCAAATGACCGAGCTCTTTTTTCTCATCCAGATACTGTGTTTTATGGCTAACGCGCTCAGGCAGCTCg                                                                                                     <  1:243877/111‑1 (MQ=255)                                                             aaGGTGACGATGATACCGAGGATCAACCGCCGCGCCCCGGCACTCAAATGACCGAGCTCTTTTTTCTCATCCAGATACTGTGTTTTATGGCTAACGCGCTCAGGCAGCTCg                                                                                                     >  2:243877/1‑111 (MQ=255)                                                                         atACCGAGGATCAACCGCCGCGCCCCGGCACTCAAATGACCGAGCTCTTTTTTCTCATCCAGATACTGTGTTTTATGGCTAACGCGCTCAGGCAGCTCGTTCATCCATTGATGGTAACGTCCACGAATATTTTGCAGAc                                                             <  2:1140079/139‑1 (MQ=255)                                                                               aGGATCAACCGCCGCGCCCCGGCACTCAGATGACCGAGCTCTTTTTTCTCATCCAGATACTGTGTTTTATGGCTAACGCGCTCAGGCAGCTCGTTCAt                                                                                                >  1:571497/1‑98 (MQ=255)                                                                               aGGATCAACCGCCGCGCCCCGGCACTCAAATGACCGAGCTCTTTTTTCTCATCCAGATACTGTGTTTTATGGCTAACGCGCTCAGGCAGCTCGTTCAt                                                                                                <  2:571497/98‑1 (MQ=255)                                                                                      aCCGCCGCGCCCCGGCACTCAAATGACCGAGCTCTTTTTTCTCATCCAGATACTGTGTTTTATGGCTAACGCGCTCa                                                                                                              <  1:176011/77‑1 (MQ=255)                                                                                      aCCGCCGCGCCCCGGCACTCAAATGACCGAGCTCTTTTTTCTCATCCAGATACTGTGTTTTATGGCTAACGCGCTCa                                                                                                              >  2:176011/1‑77 (MQ=255)                                                                                      aCCGCCGCGCCCCGGCACTCAAATGACCGAGCTCTTTTTTCTCATCCAGATACTGTGTTTTATGGCTAACGCGCTCAGGCAg                                                                                                         <  1:224943/82‑1 (MQ=255)                                                                                      aCCGCCGCGCCCCGGCACTCAAATGACCGAGCTCTTTTTTCTCATCCAGATACTGTGTTTTATGGCTAACGCGCTCAGGCAg                                                                                                         >  2:224943/1‑82 (MQ=255)                                                                                      aCCGCCGCGCCCCGGCACTCAAATGACCGAGCTCTTTTTTCTCATCCAGATACTGTGTTTTATGGCTAACGCGCTCAGGCAGCTCGTTCATCCATTGATGGTAACGTCCACGAATATTTTGCAGACCTGAAAATGCCCg                                                >  1:1242012/1‑139 (MQ=255)                                                                                       ccgccgCGCCCCGGCACTCAAATGACCGAGCTCTTTTTTCTCATCCAGATACTGTGTTTTATGGCTAACGCGCTCAgg                                                                                                            <  1:1296489/78‑1 (MQ=255)                                                                                       ccgccgCGCCCCGGCACTCAAATGACCGAGCTCTTTTTTCTCATCCAGATACTGTGTTTTATGGCTAACGCGCTCAgg                                                                                                            >  2:1296489/1‑78 (MQ=255)                                                                                           cgcgCCCGGGCACTCAAATGACCGAGCTCTTTTTTCTCATCCAGATACTGTGTTTTATGGCTAACGCGCTCAGGCAGCTCGTTCATCCATTGATGGTAACGTCCACGAATATTTTGCAGACCTGAAAATGCCCGCCTgc                                           >  1:895355/1‑139 (MQ=255)                                                                                           cgcgCCCCGGCACTCAAATGACCGAGCTCTTTTTTCTCATCCAGATACTGTGTTTTATGGCTAACGCGCTCAGGCAGCt                                                                                                       <  1:281817/79‑1 (MQ=255)                                                                                           cgcgCCCCGGCACTCAAATGACCGAGCTCTTTTTTCTCATCCAGATACTGTGTTTTATGGCTAACGCGCTCAGGCAGCt                                                                                                       >  2:281817/1‑79 (MQ=255)                                                                                           cgcgCCCCGGCACTCAAATGACCGAGCTCTTTTTTCTCATCCAGATACTGTGTTTTATGGCTAACGCGCTCAGGCAGCTCGTTCAt                                                                                                >  2:1107905/1‑86 (MQ=255)                                                                                           cgcgCCCCGGCACTCAAATGACCGAGCTCTTTTTTCTCATCCAGATACTGTGTTTTATGGCTAACGCGCTCAGGCAGCTCGTTCAt                                                                                                <  1:1107905/86‑1 (MQ=255)                                                                                           cgcgCCCCGGCACTCAAATGACCGAGCTCTTTTTTCTCATCCAGATACTGTGTTTTATGGCTAACGCGCTCAGGCAGCTCGTTCATcc                                                                                              >  1:1303846/1‑88 (MQ=255)                                                                                           cgcgCCCCGGCACTCAAATGACCGAGCTCTTTTTTCTCATCCAGATACTGTGTTTTATGGCTAACGCGCTCAGGCAGCTCGTTCATcc                                                                                              <  2:1303846/88‑1 (MQ=255)                                                                                                               aCCGAGCTCTTTTTTCTCATCCAGATACTGTGTTTTATGGCTAACGCGCTCAGGCAGCTCGTTc                                                                                                  >  1:438077/1‑64 (MQ=255)                                                                                                               aCCGAGCTCTTTTTTCTCATCCAGATACTGTGTTTTATGGCTAACGCGCTCAGGCAGCTCGTTc                                                                                                  <  2:438077/64‑1 (MQ=255)                                                                                                                 cGAGCTCTTTTTTCTCATCCAGATACTGTGTTTTATGGCTAACGCGCTCAGGCAGCTCGTTCATCCATTGATGGTAACGTCCACGAATATTTTGCATACCTGAAAATGCCCGCCTGCGCGGTTTCGGCGTTTCTTTgcg                     <  2:895355/139‑1 (MQ=255)                                                                                                                 cGAGCTCTTTTTTCTCATCCAGATACTGTGTTTTATGGCTAACGCGCTCAGGCAGCTCGTTCATCCATTGATGGTAACGTCCACGAATATTTTGCAGACCTGAAAATGCCCGCCTGCGCGGTTTCGGCGTTTCTTTgcg                     <  1:827035/139‑1 (MQ=255)                                                                                                                  gAGCTCTTTTTTCTCATCCAGATACTGTGTTTTATGGCTAACGCGCTCAGGCa                                                                                                          <  1:6939/53‑1 (MQ=255)                                                                                                                  gAGCTCTTTTTTCTCATCCAGATACTGTGTTTTATGGCTAACGCGCTCAGGCa                                                                                                          >  2:6939/1‑53 (MQ=255)                                                                                                                  gAGCTCTTTTTTCTCATCCAGATACTGTGTTTTATGGCTAACGCGCTCAGGCAGCTCGTTCATCCATTGATGGTAACGTCCACGAATATTTTGCAGACCTGaaa                                                       >  2:1335888/1‑104 (MQ=255)                                                                                                                  gAGCTCTTTTTTCTCATCCAGATACTGTGTTTTATGGCTAACGCGCTCAGGCAGCTCGTTCATCCATTGATGGTAACGTCCACGAATATTTTGCAGACCTGaaa                                                       <  1:1335888/104‑1 (MQ=255)                                                                                                                      tctTTTTTCTCTTCCAGATACTGTGTTTTATGGCTAACGCGCTCAGGCAGCTCGTTCATCCATTGATGGTAACGTCCACGAATATTTTGCAGACCTGAAAATGCCCGCCTGCGCGGTTTCGGCGTTTCTTTGCGCGATg                <  2:887518/139‑1 (MQ=255)                                                                                                                              ctcATCCAGATACTGTGTTTTATGGCTAACGCGCTCAGGCAGCTCGTTCATCCATTGATGGTAACGTCCACGAATATTTTGCAGACCTGAAAATGCCCGCCTGCGCGGTTTCGGCGTTTCTTTGCGCGATGCACCGATc        <  2:302316/139‑1 (MQ=255)                                                                                                                                  tCCAGATACTGTGTTTTATGGCTAACGCGCTCAGGCAGCTCGTTCATCCATTGATGGTAACGTCCACGAATATTTTGCAGACCTGAAAATGCCCGCCTgcg                                          <  2:185654/101‑1 (MQ=255)                                                                                                                                  tCCAGATACTGTGTTTTATGGCTAACGCGCTCAGGCAGCTCGTTCATCCATTGATGGTAACGTCCACGAATATTTTGCAGACCTGAAAATGCCCGCCTgcg                                          >  1:185654/1‑101 (MQ=255)                                                                                                                                    cAGATACTGTGTTTTATGGCTAACGCGCTCAGGCAGCTCGTTCATCCATTGATGGTAACGTCCACGAATATTTTGCAGACCTGAAAATGCCCGCCTGCGCGGTTTCGGCGTTTCTTTGCGCGATGCACCGATCAATAAc  <  1:661754/139‑1 (MQ=255)                                                                                                                                         |                                                                                                                                      ATAAACTGCGCCAGCGGGTTAAACGGCTGAGTAACGCAGATTAACGCCAGAATCAGCGAGAAGGTGACGATGATACCGAGGATCAACCGCCGCGCCCCGGCACTCAAATGACCGAGCTCTTTTTTCTCATCCAGATGCTGTGTTTTATGGCTAACGCGCTCAGGCAGCTCGTTCATCCATTGATGGTAACGTCCACGAATATTTTGCAGACCTGAAAATGCCCGCCTGCGCGGTTTCGGCGTTTCTTTGCGCGATGCACCGATCAATAAC  >  NC_000913/3694704‑3694973

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