Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,696,621	A→G	H55R (CAC→CGC)	bcsE →	cellulose production protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,696,621	0	A	G	100.0%	145.4 / NA	44	H55R (CAC→CGC)	bcsE	cellulose production protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (20/24);  total (20/24)

CCTGTATTCTCTATCGGTATCTCATCATTATGGGATGAGCTGCGACATATGCCAGCAGGCGGCGTCTGGTGGTTTAACGTCGATCGCCATGAAGATGCTATCAGTCTGGCGAATCAAACAATTGCATCCCAGGCTGAAACCGCACACGTCGCGGTCATTAGCATGGACAGCGATCCGGCGAAAATCTTTCAATTAGATGATTCTCAAGGGCCGGAAAAAATAAAATTATTTTCAATGCTAAATCATGAAAAAGGTCTATACTATTTGACCCGTGATTTGCAGTGTTCTATTGA  >  NC_000913/3696476‑3696768                                                                                                                                                  |                                                                                                                                                    ccTGTATTCTCTATCGGTATCTCATCATTATGGGATGAGCTGCGACATATGCCAGCAGGCGGCGTCTGGTGGTTTAACGTCGATCGCCATGAAGATGCTATCAGTCTGGCGAATCAAACAATTGCATCCCAGGCTGAAACCGCACGCGTcg                                                                                                                                                >  2:375393/1‑151 (MQ=255)       ttCTCTATCGGTATCTCATCATTATGGGATGAGCTGCGACATATGCCAGCAGGCGGCGTCTGGTGGTTTAACGTCGATCGCCATGAAGATGCTATCAGTCTGGCGAATCAAACAATTGCATCCCAGGCTGAAACCGCACGCGTCGCGGTc                                                                                                                                           >  2:483249/1‑150 (MQ=255)                            ttATGGGATGAGCTGCGACATATGCCAGCAGGCGGCGTCTGGTGGTTTAACGTCGATCGCCATGAAGATGCTATCAGTCTGGCGAATCAAACAATTGCATCCCAGGCTGAAACCGCACGCGTCGCGGTCATTAGCATGGACAGCGATCCgg                                                                                                                     >  1:353176/1‑151 (MQ=255)                              aTGGGATGAGCTGCGACATATGCCAGCAGGCGGCGTCTGGTGGTTTAACGTCGATCGCCATGAAGATGCTATCAGTCTGGCGAATCAAACAATTGCATCCCAGGCTGAAACCGCACGCGTCGCGGTCATTAGCATGGACAGCGATCCGGCg                                                                                                                   >  2:123718/1‑151 (MQ=255)                                  gATGAGCTGCGACATATGCCAGCAGGCGGCGTCTGGTGGTTTAACGTCGATCGCCATGAAGATGCTATCAGTCTGGCGAATCAAACAATTGCATCCCAGGCTGAAACCGCACGCGTCGCGGTCATTAGCATGGACAGCGATCCGGCGaaaa                                                                                                               <  2:138259/151‑1 (MQ=255)                                   aTGAGCTGCGACATATGCCAGCAGGCGGCGTCTGGTGGTTTAACGTCGATCGCCATGAAGATGCTATCAGTCTGGCGAATCAAACAATTGCATCCCAGGCTGAAACCGCACGCGTCGCGGTCATTAGCATGGACAGCGATCCGGCGAAAAt                                                                                                              >  1:228543/1‑151 (MQ=255)                                          gCGACATATGCCAGCAGGCGGCGTCTGGTGGTTTAACGTCGATCGCCATGAAGATGCTATCAGTCTGGCGAATCAAACAATTGCATCCCAGGCTGAAACCGCACGCGTCGCGGTCATTAGCATGGACAGCGATCCGGCGAAAATCttt                                                                                                          <  2:473844/148‑1 (MQ=255)                                          gCGACATATGCCAGCAGGCGGCGTCTGGTGGTTTAACGTCGATCGCCATGAAGATGCTATCAGTCTGGCGAATCAAACAATTGCATCCCAGGCTGAAACCGCACGCGTCGCGGTCATTAGCATGGACAGCGATCCGGCGAAAATCTTTc                                                                                                         <  2:228538/149‑1 (MQ=255)                                            gACATATGCCAGCAGGCGGCGTCTGGTGGTTTAACGTCGATCGCCATGAAGATGCTATCAGTCTGGCGAATCAAACAATTGCATCCCAGGCTGAAACCGCACGCGTCGCGGTCATTAGCATGGACAGCGATCCGGCGAAAATCTTTCAAtt                                                                                                     <  1:188128/151‑1 (MQ=255)                                               atatGCCAGCAGGCGGCGTCTGGTGGTTTAACGTCGATCGCCATGAAGATGCTATCAGTCTGGCGAATCAAACAATTGCATCCCAGGCTGAAACCGCACGCGTCGCGGTCATTAGCATGGAAAGCGATCCGGCGAAAATCTTTCAATTa                                                                                                    <  2:558873/149‑1 (MQ=255)                                                tatGCCAGCAGGCGGCGTCTGGTGGTTTAACGTCGATCGCCATGAAGATGCTATCAGTCTGGCGAATCAAACAATTGCATCCCAGGCTGAAACCGCACGCGTCGCGGTCATTAGCATGGACAGCGATCCGGCGAAAATCTTTCAAtt                                                                                                     >  2:30276/1‑147 (MQ=255)                                                 atGCCAGCAGGCGGCGTCTGGTGGTTTAACGTCGATCGCCATGAAGATGCTATCAGTCTGGCGAATCAAACAATTGCATCCCAGGCTGAAACCGCACGCGTCGCGGTCATTAGCATGGACAGCGATCCGGCGAAAATCTTTCAATTAgatg                                                                                                <  2:621056/151‑1 (MQ=255)                                                      agcagGCGGCGTCTGGTGGTTTAACGTCGATCGCCATGAAGATGCTATCAGTCTGGCGAATCAAACAATTGCATCCCAGGCTGAAACCGCACGCGTCGCGGTCATTAGCATGGACAGCGATCCGGCGAAAATCTTTCAATTAGATGATtct                                                                                           >  2:256751/1‑151 (MQ=255)                                                       gcagGCGGCGTCTGGTGGTTTAACGTCGATCGCCATGAAGATGCTATCAGTCTGGCGAATCAAACAATTGCATCCCAGGCTGAAACCGCACGCGTCGCGGTCATTAGCATGGACAGCGATCCGGCGAAAATCTTTCAATTAGATGATtctc                                                                                          <  1:414384/151‑1 (MQ=255)                                                         agGCGGCGTCTGGTGGTTTAACGTCGATCGCCATGAAGATGCTATCAGTCTGGCGAATCAAACAATTGCATCCCAGGCTGAAACCGCACGCGTCGCGGTCATTAGCATGGACAGCGATCCGGCGAAAATCTTTCAATTAGATGATTCTCaa                                                                                        <  2:518241/151‑1 (MQ=255)                                                                gTCTGGTGGTTTAACGTCGATCGCCATGAAGATGCTATCAGTCTGGCGAATCAAACAATTGCATCCCAGGCTGAAACCGCACGCGTCGCGGTCATTAGCATGGACAGCGATCCGGCGAAAATCTTTCAATTAGATGATTCTCAAGGGCCgg                                                                                 >  1:604609/1‑151 (MQ=255)                                                                  cTGGTGGTTTAACGTCGATCGCCATGAAGATGCTATCAGTCTGGCGAATCAAACAATTGCATCCCAGGCTGAAACCGCACGCGTCGCGGTCATTAGCATGGACAGCGATCCGGCGAAAATCTTTCAATTAGATGATTCTCAAGGGCCGGa                                                                                >  2:247425/1‑150 (MQ=255)                                                                   tggtggTTTAACGTCGATCGCCATGAAGATGCTATCAGTCTGGCGAATCAAACAATTGCATCCCAGGCTGAAACCGCACGCGTCGCGGTCATTAGCATGGACAGCGATCCGGCGAAAATCTTTCAATTAGATGATTCTCAAGGGCCGGaa                                                                               >  2:22525/1‑150 (MQ=255)                                                                    ggtggtTTAACGTCGATCGCCATGAAGATGCTATCAGTCTGGCGAATCAAACAATTGCATCCCAGGCTGAAACCGCACGCGTCGCGGTCATTAGCATGGACAGCGATCCGGCGAAAATCTTTCAATTAGATGATTCTCAAGGGCCGGaaa                                                                              <  1:144396/150‑1 (MQ=255)                                                                    ggtggtTTAACGTCGATCGCCATGAAGATGCTATCAGTCTGGCGAATCAAACAATTGCATCCCAGGCTGAAACCGCACGCGTCGCGGTCATTAGCATGGACAGCGATCCGGCGAAAATCTTTCAATTAGATGATTCTCAAGGGCCGGa                                                                                <  2:511631/148‑1 (MQ=255)                                                                     gtgctgtAACGTCGATCGCCATGAAGATGCTATCAGTCTGGCGAAGCAAACAATTGCATCCCAGGCTGAAACCGCACGCGTCGCGGTCATTAGCATGGACAGCGATCCGGCGAAAATCTTTCAATTAGATGATTCTCAAAGGCCGGaaaaa                                                                            <  2:126870/145‑1 (MQ=255)                                                                              cGTCGATCGCCATGAAGATGCTATCAGTCTGGCGAATCAAACAATTGCATCCCAGGCTGAAACCGCACGCGTCGCGGTCATTAGCATGGACAGCGATCCGGCGAAAATCTTTCAATTAGATGATTCTCAAGGGCCGGAAAAAATAAAatta                                                                   <  2:395900/151‑1 (MQ=255)                                                                                tCGATCGCCATGAAGATGCTATCAGTCTGGCGAATCAAACAATTGCATCCCAGGCTGAAACCGCACGCGTCGCGGTCATTAGCATGGACAGCGATCCGGCGAAAATCTTTCAATTAGATGATTCTCAAGGGCCGGAAAAAATAAAattatt                                                                 >  2:306350/1‑151 (MQ=255)                                                                                   aTCGCCATGAAGATGCTATCAGTCTGGCGAATCAAACAATTGCATCCCAGGCTGAAACCGCACGCGTCGCGGTCATTAGCATGGACAGCGATCCGGCGAAAATCTTTCAATTAGATGATTCTCAAGGGCCGGAAAAAATAAAATTATTTTc                                                              <  1:61098/151‑1 (MQ=255)                                                                                         aTGAAGATGCTATCAGTCTGGCGAATCAAACAATTGCATCCCAGGCTGAAACCGCACGCGTCGCGGTCATTAGCATGGACAGCGATCCGGCGAAAATCTTTCAATTAGATGATTCTCAAGGGCCGGAAAAAATAAAATTATTTTCAATGCt                                                        <  1:179239/151‑1 (MQ=255)                                                                                            aaGATGCTATCAGTCTGGCGAATCAAACAATTGCATCCCAGGCTGAAACCGCACGCGTCGCGGTCATTAGCATGGACAGCGATCCGGCGAAAATCTTTCAATTAGATGATTCTCAAGGGCCGGAAAAAATAAAATTATTTTCAATGCTaa                                                      >  1:600482/1‑150 (MQ=255)                                                                                            aaGATGCTATCAGTCTGGCGAATCAAACAATTGCATCCCAGGCTGAAACCGCACGCGTCGCGGTCATTAGCATGGACAGCGATCCGGCGAAAATCTTTCAATTAGATGATTCTCAAGGGCCGGAAAAAATAAAATTATTTTCAATGCTaa                                                      >  1:599962/1‑150 (MQ=255)                                                                                               aTGCTATCAGTCTGGCGAATCAAACAATTGCATCCCAGGCTGAAACCGCACGCGTCGCGGTCATTAGCATGGACAGCGATCCGGCGAAAATCTTTCAATTAGATGATTCTCAAGGGCCGGAAAAAATAAAATTATTTTCAATGCTAAATCa                                                  <  1:11612/151‑1 (MQ=255)                                                                                                    aTCAGTCTGGCGAATCAAACAATTGCATCCCAGGCTGAAACCGCACGCGTCGCGGTCATTAGCATGGACAGCGATCCGGCGAAAATCTTTCAATTAGATGATTCTCAAGGGCCGGAAAAAATAAAATTATTTTCAATGCTAAATCATGaaa                                             <  1:505479/151‑1 (MQ=255)                                                                                                      cAGTCTGGCGAATCAAACAATTGCATCCCAGGCTGAAACCGCACGCGTCGCGGTCATTAGCATGGACAGCGATCCGGCGAAAATCTTTCAATTAGATGATTCTCAAGGGCCGGAAAAAATAAAATTATTTTCAATGCTAAATCATGaaaa                                            >  2:249657/1‑150 (MQ=255)                                                                                                       aGTCTGGCGAATCAAACAATTGCATCCCAGGCTGAAACCGCACGCGTCGCGGTCATTAGCATGGACAGCGATCCGGCGAAAATCTTTCAATTAGATGATTCTCAAGGGCCGGAAAAAATAAAATTATTTTCAATGCTAAATCATGaaaa                                            >  2:164036/1‑149 (MQ=255)                                                                                                        gtatGGCGAATCAAACAATTGCACCCCAGGCTGAAACCGCACGCGTCGCGGTCATGAGCACGGACAGCGATCCAGCGAAAATCTTTCTATTAGATGATTCTCAAGGGCCGGAACAAATAAAATCATTTTCACTGCTAAATCATGAAAAAgg                                         >  2:435553/4‑151 (MQ=255)                                                                                                          cTGGCGAATCAAACAATTGCATCCCAGGCTGAAACCGCACGCGTCGCGGTCATTAGCATGGACAGCGATCCGGCGAAAATCTTTCAATTAGATGATTCTCAAGGGCCGGAAAAAATAAAATTATTTTCAATGCTAAATCATGAAAAAGGTc                                       >  2:572276/1‑151 (MQ=255)                                                                                                              cGAATCAAACAATTGCATCCCAGGCTGAAACCGCACGCGTCGCGGTCATTAGCATGGACAGCGATCCGGCGAAAATCTTTCAATTAGATGATTCTCAAGGGCCGGAAAAAATAAAATTATTTTCAATGCTAAATCATGAAAAAGGTCtata                                   <  2:529931/151‑1 (MQ=255)                                                                                                               gAATCAAACAATTGCATCCCAGGCTGAAACCGCACGCGTCGCGGTCATTAGCATGGACAGCGATCCGGCTAAAATCTTTCAATTAGATGATTCTCAAGGGCCGGAAAAAATAAAATTATTTTCAATGCTAAATCATGAAAAAGGTCTATAc                                  <  1:265657/151‑1 (MQ=255)                                                                                                                            gCATCCCAGGCTGAAACCGCACGCGTCGCGGTCATTAGCATGGACAGCGATCCGGCGAAAATCTTTCAATTAGATGATTCTCAAGGGCCGGAAAAAATAAAATTATTTTCAATGCTAAATCATGAAAAAGGTCTATACTATTTGACCCGTg                     <  1:123721/151‑1 (MQ=255)                                                                                                                            gCATCCCAGGCTGAAACCGCACGCGTCGCGGTCATTAGCATGGACAGCGATCCGGCGAAAATCTTTCAATTAGATGATTCTCAAGGGCCGGAAAAAATAAAATTATTTTCAATGCTAAATCATGAAAAAGGTCTATACTATTTGACCCGTg                     >  2:367445/1‑151 (MQ=255)                                                                                                                              aTCCCAGGCTGAAACCGCACGCGTCGCGGTCATTAGCATGGACAGCGATCCGGCGAAAATCTTTCAATTAGATGATTCTCAAGGGCCGGAAAAAATAAAATTATTTTCAATGCTAAATCATGAAAAAGGTCTATACTATTTGACCCGTGAt                   >  1:217702/1‑151 (MQ=255)                                                                                                                                cccAGGCTGAAACCGCACGCGTCGCGGTCATTAGCATGGACAGCGATCCGGCGAAAATCTTTCAATTAGATGATTCTCAAGGGCCGGAAAAAATAAAATTATTTTCAATGCTAAATCATGAAAAAGGTCTATACTATTTGACCCGTGAttt                 <  1:191403/151‑1 (MQ=255)                                                                                                                                cccAGGCTGAAACCGCACGCGTCGCGGTCATTAGCATGGACAGCGATCCGGCGAAAATCTTTCAATTAGATGATTCTCAAGGGCCGGAAAAAATAAAATTATTTTCAATGCTAAATCATGAAAAAGGTCTATACTATTTGACCCGTGAtt                  <  1:143107/150‑1 (MQ=255)                                                                                                                                          aaCCGCACGCGTCGCGGTCATTAGCATGGACAGCGATCCGGCGAAAATCTTTCAATTAGATGATTCTCAAGGGCCGGAAAAAATAAAATTATTTTCAATGCTAAATCATGAAAAAGGTCTATACTATTTGACCCGTGATTTGCAGTGTTCt       <  2:649537/151‑1 (MQ=255)                                                                                                                                           aCCGCACGCGTCGCGGTCATTAGCATGGACAGCGATCCGGCGAAAATCTTTCAATTAGATGATTCTCAAGGGCCGGAAAAAATAAAATTATTTTCAATGCTAAATCATGAAAAAGGTCTATACTATTTGACCCGTGATTTGCAGTGTTCTa      <  1:343461/151‑1 (MQ=255)                                                                                                                                            ccGCACGCGTCGCGGTCATTAGCATGGACAGCGATCCGGCGAAAATCTTTCAATTAGATGATTCTCAAGGGCCGGAAAAAATAAAATTATTTTCAATGCTAAATCATGAAAAAGGTCTATACTATTTGACCCGTGATTTGCAGTGTTCTAt     <  2:30667/151‑1 (MQ=255)                                                                                                                                              gCACGCGTCGCGGTCATTAGCATGGACAGCGATCCGGCGAAAATCTTTCAATTAGATGATTCTCAAGGGCCGGAAAAAATAAAATTATTTTCAATGCTAAATCATGAAAAAGGTCTATACTATTTGACCCGTGATTTGCAGTGTTCTATTg   >  1:225046/1‑151 (MQ=255)                                                                                                                                               cacgCGTCGCGGTCATTAGCATGGACAGCGATCCGGCGAAAATCTTTCAATTAGATGATTCTCAAGGGCCGGAAAAAATAAAATTATTTTCAATGCTAAATCATGAAAAAGGTCTATACTATTTGACCCGTGATTTGCAGTGTTCTATTg   >  2:386894/1‑150 (MQ=255)                                                                                                                                               cacgCGTCGCGGTCATTAGCATGGACAGCGATCCGGCGAAAATCTTTCAATTAGATGATTCTCAAGGGCCGGAAAAAATAAAATTATTTTCAATGCTAAATCATGAAAAAGGTCTATACTATTTGACCCGTGATTTGCAGTGTTCTATTGa  <  1:133459/151‑1 (MQ=255)                                                                                                                                                  |                                                                                                                                                    CCTGTATTCTCTATCGGTATCTCATCATTATGGGATGAGCTGCGACATATGCCAGCAGGCGGCGTCTGGTGGTTTAACGTCGATCGCCATGAAGATGCTATCAGTCTGGCGAATCAAACAATTGCATCCCAGGCTGAAACCGCACACGTCGCGGTCATTAGCATGGACAGCGATCCGGCGAAAATCTTTCAATTAGATGATTCTCAAGGGCCGGAAAAAATAAAATTATTTTCAATGCTAAATCATGAAAAAGGTCTATACTATTTGACCCGTGATTTGCAGTGTTCTATTGA  >  NC_000913/3696476‑3696768

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