Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,690,503	G→C	A702G (GCC→GGC)	bcsB ←	cellulose synthase periplasmic subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,690,503	0	G	C	100.0%	188.1 / NA	52	A702G (GCC→GGC)	bcsB	cellulose synthase periplasmic subunit
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base C (22/30);  total (22/30)

CCTACGTAATAAACGTCGCCAACACGTAGGCTGTTGATACCGGACTCGCGGATCACCGCGACCGAACCGAACATGGTGGCGCGTTTGCCGCTATCGTTCACCGCATCGTTAAGCATTTCATAACCGCGTGGGCTATCTGCCAACAGCGCAATCAC  >  NC_000913/3690425‑3690579                                                                               |                                                                             ccTACGTAATAAACGTCGCCAACACGTAGGCTGTTGATACCGGACTCGCGGATCACCGCGACCGAACCGAACATGGTGccg                                                                            <  1:708223/81‑1 (MQ=255)       taataaACGTCGCCAACACGTAGGCTGTTGATACCGGACTCGCGGATCACCGCGACCGAACCGAACATGGTGCCGCGTTTg                                                                      >  1:708306/1‑81 (MQ=255)        aataaACGTCGCCAACACGTAGGCTGTTGATACCGGACTCGCGGATCACCGCGACCGAACCGAACATGGTGCCGCGTTTGc                                                                     >  1:708310/1‑81 (MQ=255)            aaCGTCGCCAACACGTAGGCTGTTGATACCGGACTCGCGGATCACCGCGACCGAACCGAACATGGTGCCGCGTTTGCCGCt                                                                 <  2:708224/81‑1 (MQ=255)             aCGTCGCCAACACGTAGGCTGTTGATACCGGACTCGCGGATCACCGCGACCGAACCGAACATGGTGCCGCGTTTGCCGCTa                                                                <  2:708225/81‑1 (MQ=255)               gTCGCCAACACGTAGGCTGTTGATACCGGACTCGCGGATCACCGCGACCGAACCGAACATGGTGCCGCGTTTGCCGCTATc                                                              <  1:708226/81‑1 (MQ=255)                tCGCCAACACGTAGGCTGTTGATACCGGACTCGCGGATCACCGCGACCGAACCGAACATGGTGCCGCGTTTGCCGCTATCg                                                             >  2:708329/1‑81 (MQ=255)                 cGCCAACACGTAGGCTGTTGATACCGGACTCGCGGATCACCGCGACCGAACCGAACATGGTGCCGCGTTTGCCGCTATCGt                                                            >  2:708295/1‑81 (MQ=255)                  gCCAACACGTAGGCTGTTGATACCGGACTCGCGGATCACCGCGACCGAACCGAACATGGTGCCGCGTTTGCCGCTATCGtt                                                           <  1:708227/81‑1 (MQ=255)                   ccAACACGTAGGCTGTTGATACCGGACTCGCGGATCACCGCGACCGAACCGAACATGGTGCCGCGTTTGCCGCTATCGtt                                                           >  1:708313/1‑80 (MQ=255)                     aaCACGTAGGCTGTTGATACCGGACTCGCGGATCACCGCGACCGAACCGAACATGGTGCCGCGTTTGCCGCTATCGTTCAc                                                        <  1:708228/81‑1 (MQ=255)                      acacGTAGGCTGTTGATACCGGACTCGCGGATCACCGCGACCGAACCGAACATGGTGCCGCGTTTGCCGCTATCGTTCAcc                                                       <  1:708230/81‑1 (MQ=255)                      acacGTAGGCTGTTGATACCGGACTCGCGGATCACCGCGACCGAACCGAACATGGTGCCGCGTTTGCCGCTATCGTTCAcc                                                       <  2:708229/81‑1 (MQ=255)                          gTAGGCTGTTGATACCGGACTCGCGGATCACCGCGACCGAACCGAACATGGTGCCGCGTTTGCCGCTATCGTTCACCGCAt                                                   >  2:708309/1‑81 (MQ=255)                           tAGGCTGTTGATACCGGACTCGCGGATCACCGCGACCGAACCGAACATGGTGCCGCGTTTGCCGCTATCGTTCACCGCATc                                                  <  2:708231/81‑1 (MQ=255)                             ggCTGTTGATACCGGACTCGCGGATCACCGCGACCGAACCGAACATGGTGCCGCGTTTGCCGCTATCGTTCACCGCATCg                                                 <  2:708233/80‑1 (MQ=255)                             ggCTGTTGATACCGGACTCGCGGATCACCGCGACCGAACCGAACATGGTGCCGCGTTTGCCGCTATCGTTCACCGCATCg                                                 <  2:708232/80‑1 (MQ=255)                              gCTGTTGATACCGGACTCGCGGATCACCGCGACCGAACCGAACATGGTGCCGCGTTTGCCGCTATCGTTCACCGCATCGtt                                               <  2:708234/81‑1 (MQ=255)                                       aCCGGACTCGCGGATCACCGCGACCGAACCGAACATGGTGCCGCGTTTGCCGCTATCGTTCACCGCATCGTTAAGCAttt                                       >  2:708319/1‑80 (MQ=255)                                       aCCGGACTCGCGGATCACCGCGACCGAACCGAACATGGTGCCGCGTTTGCCGCTATCGTTCACCGCATCGTTAAGCAttt                                       >  1:708331/1‑80 (MQ=255)                                         cGGACTCGCGGATCACCGCGACCGAACCGAACATGGTGCCGCGTTTGCCGCTATCGTTCACCGCATCGTTAAGCATTTCa                                     >  2:708302/1‑80 (MQ=255)                                          ggACTCGCGGATCACCGCGACCGAACCGAACATGGTGCCGCGTTTGCCGCTATCGTTCACCGCATCGTTAAGCATTTCAt                                    >  2:708305/1‑80 (MQ=255)                                            aCTCGCGGATCACCGCGACCGAACCGAACATGGTGCCGCGTTTGCCGCTATCGTTCACCGCATCGTTAAGCATTTCATAAc                                 <  1:708235/81‑1 (MQ=255)                                            aCTCGCGGATCACCGCGACCGAACCGAACATGGTGCCGCGTTTGCCGCTATCGTTCACCGCATCGTTAAGCATTTCATAAc                                 <  1:708236/81‑1 (MQ=255)                                             cTCGCGGATCACCGCGACCGAACCGAACATGGTGCCGCGTTTGCCGCTATCGTTCACCGCATCGTTAAGCATTTCATAAcc                                <  2:708237/81‑1 (MQ=255)                                               cgcgGATCACCGCGACCGAACCGAACATGGTGCCGCGTTTGCCGCTATCGTTCACCGCATCGTTAAGCATTTCATAACcgc                              >  1:708340/1‑81 (MQ=255)                                               cgcgGATCACCGCGACCGAACCGAACATGGTGCCGCGTTTGCCGCTATCGTTCACCGCATCGTTAAGCATTTCATAACcgc                              >  2:708321/1‑81 (MQ=255)                                                gcgGATCACCGCGACCGAACCGAACATGGTGCCGCGTTTGCCGCTATCGTTCACCGCATCGTTAAGCATTTCATAACcgcg                             >  1:708325/1‑81 (MQ=255)                                                gcgGATCACCGCGACCGAACCGAACATGGTGCCGCGTTTGCCGCTATCGTTCACCGCATCGTTAAGCATTTCATAACcgc                              <  1:708238/80‑1 (MQ=255)                                                gcgGATCACCGCGACCGAACCGAACATGGTGCCGCGTTTGCCGCTATCGTTCACCGCATCGTTAAGCATTTCATAACcgc                              <  1:708239/80‑1 (MQ=255)                                                   gATCACCGCGACCGAACCGAACATGGTGCCGCGTTTGCCGCTATCGTTCACCGCATCGTTAAGCATTTCATAACCGCGTgg                          <  1:708240/81‑1 (MQ=255)                                                   gATCACCGCGACCGAACCGAACATGGTGCCGCGTTTGCCGCTATCGTTCACCGCATCGTTAAGCATTTCATAACCGCGTgg                          >  2:708320/1‑81 (MQ=255)                                                     tCACCGCGACCGAACCGAACATGGTGCCGCGTTTGCCGCTATCGTTCACCGCATCGTTAAGCATTTCATAACCGCGTGGGc                        <  2:708241/81‑1 (MQ=255)                                                       aCCGCGACCGAACCGAACATGGTGCCGCGTTTGCCGCTATCGTTCACCGCATCGTTAAGCATTTCATAACCGCGTGGGCt                       <  2:708244/80‑1 (MQ=255)                                                       aCCGCGACCGAACCGAACATGGTGCCGCGTTTGCCGCTATCGTTCACCGCATCGTTAAGCATTTCATAACCGCGTGGGCt                       <  1:708242/80‑1 (MQ=255)                                                       aCCGCGACCGAACCGAACATGGTGCCGCGTTTGCCGCTATCGTTCACCGCATCGTTAAGCATTTCATAACCGCGTGGGCt                       <  1:708243/80‑1 (MQ=255)                                                       aCCGCGACCGAACCGAACATGGTGCCGCGTTTGCCGCTATCGTTCACCGCATCGTTAAGCATTTCATAACCGCGTGGGCTa                      >  2:708312/1‑81 (MQ=255)                                                          gcgACCGAACCGAACATGGTGCCGCGTTTGCCGCTATCGTTCACCGCATCGTTAAGCATTTCATAACCGCGTGGGCTATCt                   >  1:708334/1‑81 (MQ=255)                                                           cgACCGAACCGAACATGGTGCCGCGTTTGCCGCTATCGTTCACCGCATCG‑TAAGCATTTCATAACCGCGTGGGCTATCTGc                 <  2:708245/81‑1 (MQ=255)                                                            gACCGAACCGAACATGGTGCCGCGTTTGCCGCTATCGTTCACCGCATCGTTAAGCATTTCATAACCGCGTGGGCTATCTGc                 >  1:708326/1‑81 (MQ=255)                                                            gACCGAACCGAACATGGTGCCGCGTTTGCCGCTATCGTTCACCGCATCGTTAAGCATTTCATAACCGCGTGGGCTATCTGc                 <  2:708246/81‑1 (MQ=255)                                                             accgaaccgaACATGGTGCCGCGTTTGCCGCTATCGTTCACCGCATCGTTAAGCATTTCATAACCGCGTGGGCTATCTGcc                <  2:708247/81‑1 (MQ=255)                                                             accgaaccgaACATGGTGCCGCGTTTGCCGCTATCGTTCACCGCATCGTTAAGCATTTCATAACCGCGTGGGCTATCTGc                 >  2:708315/1‑80 (MQ=255)                                                               cgaaccgaacATGGTGCCGCGTTTGCCGCTATCGTTCACCGCATCGTTAAGCATTTCATAACCGCGTGGGCTATCTGCCaa              >  1:708342/1‑81 (MQ=255)                                                                gaaccgaacATGGTGCCGCGTTTGCCGCTATCGTTCACCGCATCGTTAAGCATTTCATAACCGCGTGGGCTATCTGCCaa              >  2:708332/1‑80 (MQ=255)                                                                gaaccgaacATGGTGCCGCGTTTGCCGCTATCGTTCACCGCATCGTTAAGCATTTCATAACCGCGTGGGCTATCTGCCAAc             <  1:708248/81‑1 (MQ=255)                                                                  accgaacATGGTGCCGCGTTTGCCGCTATCGTTCACCGCATCGTTAAGCATTTCATAACCGCGTGGGCTATCTGCCAACAg           <  2:708250/81‑1 (MQ=255)                                                                  accgaacATGGTGCCGCGTTTGCCGCTATCGTTCACCGCATCGTTAAGCATTTCATAACCGCGTGGGCTATCTGCCAACAg           <  2:708249/81‑1 (MQ=255)                                                                    cgaacATGGTGCCGCGTTTGCCGCTATCGTTCACCGCATCGTTAAGCATTTCATAACCGCGTGGGCTATCTGCCAACAgcg         >  1:708347/1‑81 (MQ=255)                                                                       acATGGTGCCGCGTTTGCCGCTATCGTTCACCGCATCGTTAAGCATTTCATACCCGCGTGGGCTATCTGCCAACAGCGCaa      <  2:708251/81‑1 (MQ=255)                                                                       acATGGTGCCGCGTTTGCCGCTATCGTTCACCGCATCGTTAAGCATTTCATAACCGCGTGGGCTATCTGCCAACAGCGCaa      >  1:708335/1‑81 (MQ=255)                                                                         aTGGTGCCGCGTTTGCCGCTATCGTTCACCGCATCGTTAAGCATTTCATAACCGCGTGGGCTATCTGCCAACAGCGCAATc    <  2:708253/81‑1 (MQ=255)                                                                         aTGGTGCCGCGTTTGCCGCTATCGTTCACCGCATCGTTAAGCATTTCATAACCGCGTGGGCTATCTGCCAACAGCGCAATc    <  2:708252/81‑1 (MQ=255)                                                                          tggtgcCGCGTTTGCCGCTATCGTTCACCGCATCGTTAAGCATTTCATAACCGCGTGGGCTATCTGCCAACAGCGCAATCa   >  1:708323/1‑81 (MQ=255)                                                                          tggtgcCGCGTTTGCCGCTATCGTTCACCGCATCGTTAAGCATTTCATAACCGCGTGGGCTATCTGCCAACAGCGCAATCa   >  2:708330/1‑81 (MQ=255)                                                                           ggtgcCGCGTTTGCCGCTATCGTTCACCGCATCGTTAAGCATTTCATAACCGCGTGGGCTATCTGCCAACAGCGCAATCAc  >  1:708333/1‑81 (MQ=255)                                                                            gtgcCGCGTTTGCCGCTATCGTTCACCGCATCGTTAAGCATTTCATAACCGCGTGGGCTATCTGCCAACAGCGCAATCAc  <  1:708254/80‑1 (MQ=255)                                                                               |                                                                             CCTACGTAATAAACGTCGCCAACACGTAGGCTGTTGATACCGGACTCGCGGATCACCGCGACCGAACCGAACATGGTGGCGCGTTTGCCGCTATCGTTCACCGCATCGTTAAGCATTTCATAACCGCGTGGGCTATCTGCCAACAGCGCAATCAC  >  NC_000913/3690425‑3690579

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