Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,686,929	G→T	Q749K (CAG→AAG)	bcsC ←	cellulose synthase subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,686,929	0	G	T	100.0%	162.6 / NA	46	Q749K (CAG→AAG)	bcsC	cellulose synthase subunit
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base T (22/24);  total (22/24)

GGCGATAGAGGTCCGCCGCATCGCTGCGCACGCCACGTTTCAGCCAGTCATCTTTCTCGTCGTTACGGGTCAGTCGGGTAAAGGTGTCGTTATCCTGCGGACGCGTCGTAGTCACACCGGATGCGACCATGGCGTCTTTGTAGGTTTCCAGCGCCTGCGTTGGATCGCCCGCCTGCGCTTCAAACTTC  >  NC_000913/3686833‑3687020                                                                                                 |                                                                                            ggCGATAGAGGTCCGCCGCATCGCTGCGCACGCCACGTTTCAGCCAGTCATCTTTCTCGTCGTTACGGGTCAGTCGGGTAAAGGTGTCGTTATCCTTCgg                                                                                          >  2:813842/1‑100 (MQ=255) ggCGATAGAGGTCCGCCGCATCGCTGCGCACGCCACGTTTCAGCCAGTCATCTTTCTCGTCGTTACGGGTCAGTCGGGTAAAGGTGTCGTTATCCTTCgg                                                                                          >  2:122671/1‑100 (MQ=255)  gCGATAGAGGTCCGCCGCATCGCTGCGCACGCCACGTTTCAGCCAGTCATCTTTCTCGTCGTTACGGGTCAGTCGGGTAAAGGTGTCGTTATCCTTCGGa                                                                                         >  1:301708/1‑100 (MQ=255)    gATAGAGGTCCGCCGCATCGCTGCGCACGCCACGTTTCAGCCAGTCATCTTTCTCGTCGTTACGGGTCAGTCGGGTAAAGGTGTCGTTATCCTTCGGAcg                                                                                       <  2:611229/100‑1 (MQ=255)        gagGTCCGCCGCATCGCTGCGCACGCCACGTTTCAGCCAGTCATCTTTCTCGTCGTTACGGGTCAGTCGGGTAAAGGTGTCGTTATCCTTCGGACGcgtc                                                                                   >  1:763141/1‑100 (MQ=255)          ggTCCGCCGCATCGCTGCGCACGCCACGTTTCAGCCAGTCATCTTTCTCGTCGTTACGGGTCAGTCGGGTAAAGGTGTCGTTATCCTTCGGACGcgtcgt                                                                                 <  1:565394/100‑1 (MQ=255)            tCCGCCGCATCGCTGCGCACGCCACGTTTCAGCCAGTCATCTTTCTCGTCGTTACGGGTCAGTCGGGTAAAGGTGTCGTTATCCTTCGGACGCGTCGTAg                                                                               <  1:715551/100‑1 (MQ=255)             ccgccgCATCGCTGCGCACGCCACGTTTCAGCCAGTCATCTTTCTCGTCGTTACGGGTCAGTCGGGTAAAGGTGTCGTTATCCTTCGGACGCGTCGTAGt                                                                              <  1:1247973/100‑1 (MQ=255)             ccgccgCATCGCTGCGCACGCCACGTTTCAGCCAGTCATCTTTCTCGTCGTTACGGGTCAGTCGGGTAAAGGTGTCGTTATCCTTCGGACGCGTCGTAGt                                                                              >  1:356136/1‑100 (MQ=255)              cgccgcATCGCTGCGCACGCCACGTTTCAGCCAGTCATCTTTCTCGTCGTTACGGGTCAGTCGGGTAAAGGTGTCGTTATCCTTCGGACGCGTCGTAGTc                                                                             <  1:1166757/100‑1 (MQ=255)              cgccgcATCGCTGCGCACGCCACGTTTCAGCCAGTCATCTTTCTCGTCGTTACGGGTCAGTCGGGTAAAGGTGTCGTTATCCTTCGGACGCGTCGTAGTc                                                                             >  1:17650/1‑100 (MQ=255)               gccgcATCGCTGCGCACGCCACGTTTCAGCCAGTCATCTTTCTCGTCGTTACGGGTCAGTCGGGTAAAGGTGTCGTTATCCTTCGGACGCGTCGTAGTca                                                                            >  2:153005/1‑100 (MQ=255)                  gcATCGCTGCGCACGCCACGTTTCAGCCAGTCATCTTTCTCGTCGTTACGGGTCAGTCGGGTAAAGGTGTCGTTATCCTTCGGACGCGTCGTAGTCacac                                                                         >  2:597928/1‑100 (MQ=255)                         tGCGCACGCCACGTTTCAGCCAGTCATCTTTCTCGTCGTTACGGGTCAGTCGGGTAAAGGTGTCGTTATCCTTCGGACGCGTCGTAGTCACACCGGATGc                                                                  <  1:142979/100‑1 (MQ=255)                          gcgcACGCCACGTTTCAGCCAGTCATCTTTCTCGTCGTTACGGGTCAGTCGGGTAAAGGTGTCGTTATCCTTCGGACGCGTCGTAGTCACACCGGATGc                                                                  <  1:151569/99‑1 (MQ=255)                                   aCGTTTCAGCCAGTCATCTTTCTCGTCGTTACGGGTCAGTCGGGTAAAGGTGTCGTTATCCTTCGGACGCGTCGTAGTCACACCGGATGCGACCATGGCg                                                        <  2:15220/100‑1 (MQ=255)                                     gTTTCAGCCAGTCATCTTTCTCGTCGTTACGGGTCAGTCGGGTAAAGGTGTCGTTATCCTTCGGACGCGTCGTAGTCACACCGGATGCGACCATGGCGTc                                                      >  2:713994/1‑100 (MQ=255)                                     gTTTCAGCCAGTCATCTTTCTCGTCGTTACGGGTCAGTCGGGTAAAGGTGTCGTTATCCTTCGGACGCGTCGTAGTCACACCGGATGCGACCATGGCGTc                                                      <  1:724932/100‑1 (MQ=255)                                     gTTTCAGCCAGTCATCTTTCTCGTCGTTACGGGTCAGTCGGGTAAAGGTGTCGTTATCCTTCGGACGCGTCGTAGTCACACCGGATGCGACCATGGCGTc                                                      <  1:756673/100‑1 (MQ=255)                                       tctaGCCAGTCATCTTTCTCGTCGTTACGGGTCAGTCGGGTAAAGGTGTCGTTATCCTTCGGACGCGTCGTAGTCACACCGGATGCGACCATGGCGTCtt                                                    <  2:983102/97‑1 (MQ=255)                                        tCAGCCAGTCATCTTTCTCGTCGTTACGGGTCAGTCGGGTAAAGGTGTCGTTATCCTTCGGACGCGTCGTAGTCACACCGGATGCGACCATGGCGTCttt                                                   >  1:821016/1‑100 (MQ=255)                                          aGCCAGTCATCTTTCTCGTCGTTACGGGTCAGTCGGGTAAAGGTGTCGTTATCCTTCGGACGCGTCGTAGTCACACCGGATGCGACCATGGCGTCttaga                                                 >  1:983110/1‑97 (MQ=255)                                          aGCCAGTCATCTTTCTCGTCGTTACGGGTCAGTCGGGTAAAGGTGTCGTTATCCTTCGGACGCGTCGTAGTCACACCGGATGCGACCATGGCGTCTTTGt                                                 >  2:599534/1‑100 (MQ=255)                                           gCCAGTCATCTTTCTCGTCGTTACGGGTCAGTCGGGTAAAGGTGTCGTTATCCTTCGGACGCGTCGTAGTCACACCGGATGCGACCATGGCGTCTTTGTa                                                <  1:122673/100‑1 (MQ=255)                                               gTCATCTTTCTCGTCGTTACGGGTCAGTCGGGTAAAGGTGTCGTTATCCTTCGGACGCGTCGTAGTCACACCGGATGCGACCATGGCGTCTTTGTAGGtt                                            >  1:862640/1‑100 (MQ=255)                                                  aTCTTTCTCGTCGTTACGGGTCAGTCGGGTAAAGGTGTCGTTATCCTTCGGACGCGTCGTAGTCACACCGGATGCGACCATGGCGTCTTTGTAGGTTTcc                                         <  2:1238973/100‑1 (MQ=255)                                                  aTCTTTCTCGTCGTTACGGGTCAGTCGGGTAAAGGTGTAGTTATCCTTCGGACGCGTCGTAGTCACACCGGATGCGACCATGGCGTCTTTGTAGGTTTcc                                         <  2:49046/100‑1 (MQ=255)                                                    cTTTCTCGTCGTTACGGGTCAGTCGGGTAAAGGTGTCGTTATCCTTCGGACGCGTCGTAGTCACACCGGATGCGACCATGGCGTCTTTGTAGGTTTCCAg                                       >  2:781784/1‑100 (MQ=255)                                                      ttCTCGTCGTTACGGGTCAGTCGGGTAAAGGTGTCGTTATCCTTCGGACGCGTCGTAGTCACACCGGATGCGACCATGGCGTCTTTGTAGGTTTCCAgcg                                     <  2:426426/100‑1 (MQ=255)                                                        ctcGTCGTTACGGGTCAGTCGGGTAAAGGTGTCGTTATCCTTCGGACGCGTCGTAGTCACACCGGATGCGACCATGGCGTCTTTGTAGGTTTCCAGCGcc                                   <  1:173411/100‑1 (MQ=255)                                                        ctcGTCGTTACGGGTCAGTCGGGTAAAGGTGTCGTTATCCTTCGGACGCGTCGTAGTCACACCGGATGCGACCATGGCGTCTTTGTAGGTTTCCAGCGcc                                   >  1:63656/1‑100 (MQ=255)                                                          cgtcgtTACGGGTCAGTCGGGTAAAGGTGTCGTTATCCTTCGGACGCGTCGTAGTCACACCGGATGCGACCATGGCGTCTTTGTAGGTTTCCAGCGCCTg                                 >  2:1262647/1‑100 (MQ=255)                                                            tcgtTACGGGTCAGTCGGGTAAAGGTGTCGTTATCCTTCGGACGCGTCGTAGTCACACCGGATGCGACCATGGCGTCTTTGTAGGTTTCCAGCGCCTGCg                               <  1:742433/100‑1 (MQ=255)                                                               ttACGGGTCAGTCGGGTAAAGGTGTCGTTATCCTTCGGACGCGTCGTAGTCACACCGGATGCGACCATGGCGTCTTTGTAGGTTTCCAGCGCCTGCGTTg                            <  1:402365/100‑1 (MQ=255)                                                                tACGGGTCAGTCGGGTAAAGGTGTCGTTATCCTTCGGACGCGTCGTAGTCACACCGGATGCGACCATGGCGTCTTTGTAGGTTTCCAGCGCCTGCGTTgg                           >  1:1122365/1‑100 (MQ=255)                                                                  cGGGTCAGTCGGGTAAAGGTGTCGTTATCCTTCGGACGCGTCGTAGTCACACCGGATGCGACCATGGCGTCTTTGTAGGTTTCCAGCGCCTGCGTTGGAt                         >  1:461656/1‑100 (MQ=255)                                                                   gggTCAGTCGGGTAAAGGTGTCGTTATCCTTCGGACGCGTCGTAGTCACACCGGATGCGACCATGGCGTCTTTGTAGGTTTCCAGCGCCTGCGTTGGATc                        <  2:17650/100‑1 (MQ=255)                                                                    ggTCAGTCGGGTAAAGGTGTCGTTATCCTTCGGACGCGTCGTAGTCACACCGGATGCGACCATGGCGTCTTTGTAGGTTTCCAGCGCCTGCGTTGGATcg                       <  1:597932/100‑1 (MQ=255)                                                                    ggTCAGTCGGGTAAAGGTGTAGTTATCATTCGGACGCGTCGTAGTCACACCGGATGCGACCATAGCGACTTTGTAGGTTTCCAGCGCCTGCGTTGGATcg                       >  2:589013/1‑100 (MQ=255)                                                                         gTCGGGTAAAGGTGTCGTTATCCTTCGGACGCGTCGTAGTCACACCGGATGCGACCATGGCGTCTTTGTAGGTTTCCAGCGCCTGCGTTGGATcgcccgc                  >  1:135568/1‑100 (MQ=255)                                                                          tCGGGTAAAGGTGTCGTTATCCTTCGGACGCGTCGTAGTCACACCGGATGCGACCATGGCGTCTTTGTAGGTTTCCAGCGCCTGCGTTGGATcgcccgcc                 <  1:778260/100‑1 (MQ=255)                                                                              gTAAAGGTGTCGTTATCCTTCGGACGCGTCGTAGTCACACCGGATGCGACCATGGCGTCTTTGTAGGTTTCCAGCGCCTGCGTTGGATCGCCCGCCTgcg             <  1:970857/100‑1 (MQ=255)                                                                                  aGGTGTCGTTATCCTTCGGACGCGTCGTAGTCACACCGGATGCGACCATGGCGTCTTTGTAGGTTTCCAGCGCCTGCGTTGGATCGCCCGCCTGCGCTTc         >  2:1230373/1‑100 (MQ=255)                                                                                    gtgtCGTTATCCTTCGGACGCGTCGTAGTCACACCGGATGCGACCATGGCGTCTTTGTAGGTTTCCAGCGCCTGCGTTGGATCGCCCGCCTGCGCTTCaa       <  2:465764/100‑1 (MQ=255)                                                                                     tgtCGTTATCCTTCGGACGCGTCGTAGTCACACCGGATGCGACCATGGCGTCTTTGTAGGTTTCCAGCGCCTGCGTTGGATCGCCCGCCTGCGCTTCaaa      >  2:196569/1‑100 (MQ=255)                                                                                         gTTATCCTTCGGACGCGTCGTAGTCACACCGGATGCGACCATGGCGTCTTTGTAGGTTTCCAGCGCCTGCGTTGGATCGGCCGCCTGCGCTTCAAACTTc  <  2:796955/100‑1 (MQ=255)                                                                                                 |                                                                                            GGCGATAGAGGTCCGCCGCATCGCTGCGCACGCCACGTTTCAGCCAGTCATCTTTCTCGTCGTTACGGGTCAGTCGGGTAAAGGTGTCGTTATCCTGCGGACGCGTCGTAGTCACACCGGATGCGACCATGGCGTCTTTGTAGGTTTCCAGCGCCTGCGTTGGATCGCCCGCCTGCGCTTCAAACTTC  >  NC_000913/3686833‑3687020

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑