Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,686,902	T→C	T758A (ACC→GCC)	bcsC ←	cellulose biosynthesis protein BcsC

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,686,902	0	T	C	100.0%	134.6 / NA	37	T758A (ACC→GCC)	bcsC	cellulose biosynthesis protein BcsC
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (19/18);  total (19/18)

CCACCGGTGCCGCTCGAACCCCAGTAATCGTGCTCAAGGGTGACGTTAAGATCCTGCTGGCGATAGAGGTCCGCCGCATCGCTGCGCACGCCACGTTTCAGCCAGTCATCTTTCTCGTCGTTACGGGTCAGTCGGGTAAAGGTGTCGTTATCCTGCGGACGCGTCGTAGTCACACCGGATGCGACCATGGCGTCTTTGTAGGTTTCCAGCGCCTGCGTTGGATCGCCCGCCTGCGCTTCAAACTTCGCACCATCACGCAGCA  >  NC_000913/3686775‑3687036                                                                                                                                |                                                                                                                                       ccaccGGTGCCGCTCGAACCCCAGTAATCGTGCTCAAGGGTGACGTTAAGATCCTGCTGGCGATAGAGGTCCGCCGCATCGCTGCGCACGCCACGTTTCAGCCAGTCATCTTTCTCGTCGTTACGGGCCAg                                                                                                                                     <  2:660209/131‑1 (MQ=255) ccaccGGTGCCGCTCGAACCCCAGTAATCGTGCTCAAGGGTGACGTTAAGATCCTGCTGGCGATAGAGGTCCGCCGCATCGCTGCGCACGCCACGTTTCAGCCAGTCATCTTTCTCGTCGTTACGGGCCAg                                                                                                                                     >  1:660209/1‑131 (MQ=255)  caccGGTGCCGCTCGAACCCCAGTAATCGTGCTCAAGGGTGACGTTAAGATCCTGCTGGCGATAGAGGTCCGCCGCATCGCTGCGCACGCCACGTTTCAGCCAGTCATCTTTCTCGTCGTTACGGGCCAGTCGGGTAAAg                                                                                                                           <  2:262531/140‑1 (MQ=255)     cGGTGCCGCTCGAACCCCAGTAATCGTGCTCAAGGGTGACGTTAAGATCCTGCTGGCGATAGAGGTCCGCCGCATCGCTGCGCACGCCACGTTTCAGCCAGTCATCTTTCTCGTCGTTACGGGCCAGTCGGGTAAAGgtg                                                                                                                        >  2:160953/1‑140 (MQ=255)      ggTGCCGCTCGAACCCCAGTAATCGTGCTCAAGGGTGACGTTAAGATCCTGCTGGCGATAGAGGTCCGCCGCATCGCTGCGCACGCCACGTTTCAGCCAGTCATCTTTCTCGTCGTTACGGGCCAGTCGGGTAAAGgtgt                                                                                                                       >  2:770441/1‑140 (MQ=255)      ggTGCCGCTCGAACCCCAGTAATCGTGCTCAAGGGTGACGTTAAGATCCTGCTGGCGATAGAGGTCCGCCGCATCGCTGCGCACGCCACGTTTCAGCCAGTCATCTTTCTCGTCGTTACGGGCCAGTCGGGTAAAGgtgt                                                                                                                       >  2:486489/1‑140 (MQ=255)        tGCCGCTCGAACCCCAGTAATCGTGCTCAAGGGTGACGTTAAGATCCTGCTGGCGATAGAGGTCCGCCGCATCGCTGCGCACGCCACGTTTCAGCCAGTCATCTTTCTCGTCGTTACGGGCCAGTCGGGTAAAGGTGTCg                                                                                                                     <  2:215511/140‑1 (MQ=255)          ccGCTCGAACCCCAGTAATCGTGCTCAAGGGTGACGTTAAGATCCTGCTGGCGATAGAGGTCCGCCGCATCGCTGCGCACGCCACGTTTCAGCCAGTCATCTTTCTCGTCGTTACGGGCCAGTCGGGTAAAGGTGTCGtt                                                                                                                   >  1:263124/1‑140 (MQ=255)          ccGCTCGAACCCCAGTAATCGTGCTCAAGGGTGACGTTAAGATCCTGCTGGCGATAGAGGTCCGCCGCATCGCTGCGCACGCCACGTTTCAGCCAGTCATCTTTCTCGTCGTTACGGGCCAGTCGGGTAAAGGTGTCGtt                                                                                                                   >  2:362506/1‑140 (MQ=255)            gCTCGAACCCCAGTAATCGTGCTCAAGGGTGACGTTAAGATCCTGCTGGCGATAGAGGTCCGCCGCATCGCTGCGCACGCCACGTTTCAGCCAGTCATCTTTCTCGTCGTTACGGGCCAGTCGGGTAAAGGTGTCGTTAt                                                                                                                 >  1:144528/1‑140 (MQ=255)             cTCGAACCCCAGTAATCGTGCTCAAGGGTGACGTTAAGATCCTGCTGGCGATAGAGGTCCGCCGCATCGCTGCGCACGCCACGTTTCAGCCAGTCATCTTTCTCGTCGTTACGGGCCAGTCGGGTAAAGGTGTCGTTATc                                                                                                                >  2:484720/1‑140 (MQ=255)             cTCGAACCCCAGTAATCGTGCTCAAGGGTGACGTTAAGATCCTGCTGGCGATAGAGGTCCGCCGCATCGCTGCGCACGCCACGTTTCAGCCAGTCATCTTTCTCGTCGTTACGGGCCAGTCGGGTAAAGGTGTCGTTATc                                                                                                                >  1:93729/1‑140 (MQ=255)                               ggCTAAAGGGTAACGTTAAAATCCTGCTGGCGATAGAGGTCCGCCGCAGCGCTGCGCACGCCACGTTTCAGCCAGTCATCTTTCTCGTCGTTACGGGCCAGTCGGGTAAAGGTGTCGTTATCCTGCGGACGCGTCGTAGt                                                                                              <  2:616560/139‑1 (MQ=255)                                gCTCAAGGGTGACGTTAAGATCCTGCTGGCGATAGAGGTCCGCCGCATCGCTGCGCACGCCACGTTTCAGCCAGTCATCTTTCTCGTCGTTACGGGCCAGTCGGGTAAAGGTGTCGTTATCCTGCGGACGCGTCGTAGTc                                                                                             >  1:760657/1‑140 (MQ=255)                                    aaGGGTGACGTTAAGATCCTGCTGGCGATAGAGGTCCGCCGCATCGCTGCGCACGCCACGTTTCAGCCAGTCATCTTTCTCGTCGTTACGGGCCAGTCGGGTAAAGGTGTCGTTATCCTGCGGACGCGTCGTAGTCacac                                                                                         <  2:11604/140‑1 (MQ=255)                                              ttAAGATCCTGCTGGCGATAGAGGTCCGCCGCATCGCTGCGCACGCCACGTTTCAGCCAGTCATCTTTCTCGTCGTTACGGGCCAGTCGGGTAAAGGTGTCGTTATCCTGCGGACGCGTCGTAGTCACACCGGATGCGAc                                                                               >  1:703495/1‑140 (MQ=255)                                              ttAAGATCCTGCTGGCGATAGAGGTACGCCGCATCGCTGCGCACGCCACGTTTCAGCCAGTCATCTTTCTCGTCGTTACGGGCCAGTCGGGTAAAGGTGTCGTTATCCTGCGGACGCGTCGTAGTCACACCGGATGCGAc                                                                               >  2:356733/1‑140 (MQ=255)                                                aaGATCCTGCTGGCGATAGAGGTCCGCCGCATCGCTGCGCACGCCACGTTTCAGCCAGTCATCTTTCTCGTCGTTACGGGCCAGTCGGGTAAAGGTGTCGTTATCCTGCGGACGCGTCGTAGTCACACCGGATGCGACCa                                                                             <  1:160953/140‑1 (MQ=255)                                                     ccTGCTGGCGATAGAGGTCCGCCGCATCGCTGCGCACGCCACGTTTCAGCCAGTCATCTTTCTCGTCGTTACGGGCCAGTCGGGTAAAGGTGTCGTTATCCTGCGGACGCGTCGTAGTCACACCGGATGCGACCATGGCg                                                                        <  2:613306/140‑1 (MQ=255)                                                                      tCCGCCGCATCGCTGCGCACGCCACGTTTCAGCCAGTCATCTTTCTCGTCGTTACGGGCCAGTCGGGTAAAGGTGTCGTTATCCTGCGGACGCGTCGTAGTCACACCGGATGCGACCATGGCGTCTTTGTAGGTTTCCAg                                                       >  2:456487/1‑140 (MQ=255)                                                                      gccgccgCGTCGCTGCGCACGCCGCGTTTCAGCCAGTCATCTTTCTCGTCGTTACGGGCCAGTCGGGTAAAGGGGGCGTTATCCTGCGGGCGCGTCGTGGTCACGCCGGATGCGACCATGGCGTCTTTGTAGGTTTCCAg                                                       <  2:273886/139‑1 (MQ=255)                                                                         gccgcATCGCTGCGCACGCCACGTTTCAGCCAGTCATCTTTCTCGTCGTTACGGGCCAGTCGGGTAAAGGTGTCGTTATCCTGCGGACGCGTCGTAGTCACACCGGATGCGACCATGGCGTCTTTGTAGGTTTCCAgcgc                                                    <  1:193966/140‑1 (MQ=255)                                                                         gccgcATCGCTGCGCACGCCACGTTTCAGCCAGTCATCTTTCTCGTCGTTACGGGCCAGTCGGGTAAAGGTGTCGTTATCCTGCGGACGCGTCGTAGTCACACCGGATGCGACCATGGCGTCTTTGTAGGTTTCCAgcgc                                                    <  2:760657/140‑1 (MQ=255)                                                                         gccgcATCGCTGCGCACGCCACGTTTCAGCCAGTCATCTTTCTCGTCGTTACGGGCCAGTCGGGTAAAGGTGTCGTTATCCTGCGGACGCGTCGTAGTCACACCGGATGCGACCATGGCGTCTTTGTAGGTTTCCAgcgc                                                    <  2:745266/140‑1 (MQ=255)                                                                         gccgcATCGCTGCGCACGCCACGTTTCAGCCAGTCATCTTTCTCGTCGTTACGGGCCAGTCGGGTAAAGGTGTCGTTATCCTGCGGACGCGTCGTAGTCACACCGGATGCGACCATGGCGTCTTTGTAGGTTTCCAgcgc                                                    <  1:456487/140‑1 (MQ=255)                                                                         gccgcATCGCTGCGCACGCCACGTTTCAGCCAGTCATCTTTCTCGTCGTTACGGGCCAGTCGGGTAAAGGTGTCGTTATCCTGCGGACGCGTCGTAGTCACACCGGATGCGACCATGGCGTCTTTGTAGGTTTCCAgcgc                                                    <  1:356733/140‑1 (MQ=255)                                                                          ccgcATCGCTGCGCACGCCACGTTTCAGCCAGTCATCTTTCTCGTCGTTACGGGCCAGTTGGGTAAAGGTGTCGTTATCCTGCGGACGCGTCGTAGTCACACCGGATGCGACCATGGCGTCTTTGTAGGTTTCCAGCGcc                                                   >  2:790852/1‑140 (MQ=255)                                                                             cATCGCTGCGCACGCCACGTTTCAGCCAGTCATCTTTCTCGTCGTTACGGGCCAGTCGGGTAAAGGTGTCGTTATCCTGCGGACGCGTCGTAGTCACACCGGATGCGACCATGGCGTCTTTGTAGGTTTCCAGCGCCTGc                                                >  2:680065/1‑140 (MQ=255)                                                                                 gCTGCGCACGCCACGTTTCAGCCAGTCATCTTTCTCGTCGTTACGGGCCAGTCGGGTAAAGGTGTCGTTATCCTGCGGACGCGTCGTAGTCACACCGGATGCGACCATGGCGTCTTTGTAGGTTTCCAGCGCCTGCGTTg                                            <  1:341448/140‑1 (MQ=255)                                                                                   tGCGCACGCCACGTTTCAGCCAGTCATCTTTCTCGTCGTTACGGGCCAGTCGGGTAAAGGTGTCGTTATCCTGCGGACGCGTCGTAGTCACACCGGATGCGACCATGGCGTCTTTGTAGGTTTCCAGCGCCTGCGTTGGa                                          <  1:362506/140‑1 (MQ=255)                                                                                     cgcACGCCACGTTTCAGCCAGTCATCTTTCTCGTCGTTACGGGCCAGTAGGGTAAAGGTTTCGTTATCCTGCGGACGCGTCGTAGTCACACCGGATGCGACCATGGCGTCTTTGTAGGTTTCCAGCGCCTGCGTTGGATc                                        >  2:573631/1‑140 (MQ=255)                                                                                          gCCACGTTTCAGCCAGTCATCTTTCTCGTCGTTACGGGCCAGTCGGGTAAAGGTGTCGTTATCCTGCGGACGCGTCGTAGTCACACCGGATGCGACCATGGCGTCTTTGTAGGTTTCCAGCGCCTGCGTTGGATcgcccg                                   >  1:548000/1‑140 (MQ=255)                                                                                                    aGCCAGTCATCTTTCTCGTCGTTACGGGCCAGTCGGGTAAAGGTGTCGTTATCCTGCGGACGCGTCGTAGTCACACCGGATGCGACCATGGCGTCTTTGTAGGTTTCCAGCGCCTGCGTTGGATCGCCCGCCTGCGCTTc                         >  2:492747/1‑140 (MQ=255)                                                                                                            aTCTTTCTCGTCGTTACGGGCCAGTCGGGTAAAGGTGTCGTTATCCTGCGGACGCGTCGTAGTCACACCGGATGCGACCATGGCGTCTTTGTAGGTTTCCAGCGCCTGCGTTGGATCGCCCGCCTGCGCTTCAAACTTCg                 <  1:770441/140‑1 (MQ=255)                                                                                                               tttCTCGTCGTTACGGGCCAGTTGGGTAAAGGTGTCGTTATCCTGCGGACGCGTCGTAGTCACACCGGATGCGACCATGGCGTCTTTGTAGGTTTCCAGCGCCTGCGTTGGATCGCCCGCCTGCGCTTCAAACTTCGCAc              <  1:790852/140‑1 (MQ=255)                                                                                                                  ctcGTCGTTACGGGCCAGTCGGGTAAAGGTGTCGTTATCCTGCGGACGCGTCGTAGTCACACCGGATGCGACCATGGCGTCTTTGTAGGTTTCCAGCGCCTGCGTTGGATCGCCCGCCTGCGCTTCAAACTTCGCACCAt           <  1:492747/140‑1 (MQ=255)                                                                                                                    cgtcgtTACGGGCCAGTCGGGTAAAGGTGTCGTTATCCTGCGGACGCGTCGTAGTCACACCGGATGCGACCATGGCGTCTTTGTAGGTTTCCAGCGCCTGCGTTGGATCGCCCGCCTGCGCTTCAAACTTCGCACCATCa         <  1:483021/140‑1 (MQ=255)                                                                                                                           aCGGGCCAGTCGGGTAAAGGTGTCGTTATCCTGCGGACGCGTCGTAGTCACACCGGATGCGACCATGGCGTCTTTGTAGGTTTCCAGCGCCTGCGTTGGATCGCCCGCCTGCGCTTCAAACTTCGCACCATCACgcagca  >  2:668901/1‑140 (MQ=255)                                                                                                                                |                                                                                                                                       CCACCGGTGCCGCTCGAACCCCAGTAATCGTGCTCAAGGGTGACGTTAAGATCCTGCTGGCGATAGAGGTCCGCCGCATCGCTGCGCACGCCACGTTTCAGCCAGTCATCTTTCTCGTCGTTACGGGTCAGTCGGGTAAAGGTGTCGTTATCCTGCGGACGCGTCGTAGTCACACCGGATGCGACCATGGCGTCTTTGTAGGTTTCCAGCGCCTGCGTTGGATCGCCCGCCTGCGCTTCAAACTTCGCACCATCACGCAGCA  >  NC_000913/3686775‑3687036

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