Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,693,326	C→G	R637R (CGG→CGC)	bcsA ←	cellulose synthase, catalytic subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,693,326	0	C	G	100.0%	106.4 / NA	31	R637R (CGG→CGC)	bcsA	cellulose synthase, catalytic subunit
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base G (14/17);  total (14/17)

CAAGAACAATCAGGTTGTAGAACACCCACACCATACTGACGACCACGGTGAGCATCTCGGTTGGCGGGCCATAGAAGTAGCGCCAGATGCCTACCGCAACGCCCACCAGGTTGAGCAGGACAAGGAAGATGTAGGGCCGCGAGATCACCCAGTCGACGTACTCTTCTTCCACCAGTCCACCTTTGGCGGTGACGTTAAATTTGCCTTTGTGCGGGTTAATCAGCGCCACCAGCGTCGGTGGTGCGATATACCACGCCAGCACCGTTTCGTAGATTT  >  NC_000913/3693190‑3693465                                                                                                                                         |                                                                                                                                            cAAGAACAATCAGGTTGTAGAACACCCACACCATACTGACGACCACGGTGAGCATCTCGGTTGGCGGGCCATAGAAGTAGCGCCAGATGCCTACCGCAACGCCCACCAGGTTGAGCAGGACAAGGAAGATGTAGGGGCGCGAGATCAcc                                                                                                                                 <  3:53479/149‑1 (MQ=255)   aGAACAATCAGGTTGTAGAACACCCACACCATACTGACGACCACGGTGAGCATCTCGGTTGGCGGGCCATAGAAGTAGCGCCAGATGCCTACCGCAACGCCCACCAGGTTGAGCAGGACAAGGAAGATGTAGGGGCGCGAGATCACCCa                                                                                                                               <  3:13349/149‑1 (MQ=255)             ggTTGTAGAACACCCACACCATACTGACGACCACGGGGAGCATCTCGGTGGGCGGGCCATAGAAGTAGCGCCAGATGCCTACCGCAACGCCCACCAGGTTGAGCAGGACAAGGAAGATGTAGGGGCGCGAGATCACCCAGTCGACGTAc                                                                                                                     <  2:152163/149‑1 (MQ=255)                   aGAACACCCACACCATACTGACGACCACGGTGAGCATCTCGGTTGGCGGGCCATAGAAGTAGCGCCAGATGCCTACCGCAACGCCCACCAGGTTGAGCAGGACAAGGAAGATGTAGGGGCGCGAGATCACCCAGTCGACGTACtcttct                                                                                                               >  2:30798/1‑149 (MQ=255)                    gAACACCCACACCATACTGACGACCACGGTGAGCATCTCGGTTGGCGGGCCATAGAAGTAGCGCCAGATGCCTACCGCAACGCCCACCAGGTTGAGCAGGACAAGGAAGATGTAGGGGCGCGAGATCACCCAGTCGACGTACtcttct                                                                                                               <  3:23865/148‑1 (MQ=255)                     aaCACCCACACCATACTGACGACCACGGTGAGCATCTCGGTTGGCGGGCCATAGAAGTAGCGCCAGATGCCTACCGCAACGCCCACCAGGTTGAGCAGGACAAGGAAGATGTAGGGGCGCGAGATCACCCAGTCGACGTACTCttcttc                                                                                                             <  6:27939/149‑1 (MQ=255)                          ccACACCATACTGACGACCACGGTGAGCATCTCGGTTGGCGGGCCATAGAAGTAGCGCCAGATGCCTACCGCAACGCCCACCAGGTTGAGCAGGACAAGGAAGATGTAGGGGCGCGAGATCACCCAGTCGACGTACTCTTCTTccac                                                                                                          >  5:149725/1‑147 (MQ=255)                           cacaCCATACTGACGACCACGGTGAGCATCTCGGTTGGCGGGCCATAGAAGTAGCGCCAGATGCCTACCGCAACGCCCACCAGGTTGAGCAGGACAAGGAAGATGTAGGGGCGCGAGATCACCCAGTCGACGTACTCTTCTTCCACCAg                                                                                                       >  6:54610/1‑149 (MQ=255)                            acacCATACTGACGACCACGGTGAGCATCTCGGTTGGCGGGCCATAGAAGTAGCGCCAGATGCCTACCGCAACGCCCACCAGGTTGAGCAGGACAAGGAAGATGTAGGGGCGCGAGATCACCCAGTCGACGTACTCTTCTTCCACCAGt                                                                                                      <  1:30798/149‑1 (MQ=255)                            acacCATACTGACGACCACGGTGAGCATCTCGGTTGGCGGGCCATAGAAGTAGCGCCAGATGCCTACCGCAACGCCCACCAGGTTGAGCAGGACAAGGAAGATGTAGGGGCGCGAGATCACCCAGTCGACGTACTCTTCTTCCACCAGt                                                                                                      >  1:44806/1‑149 (MQ=255)                               ccATACTGACGACCACGGTGAGCATCTCGGTTGGCGGGCCATAGAAGTAGCGCCAGATGCCTACCGCAACGCCCACCAGGTTGAGCAGGACAAGGAAGATGTAGGGGCGCGAGATCACCCAGTCGACGTACTCTTCTTCCACCAGTCCa                                                                                                   >  6:1713/1‑149 (MQ=255)                                     tGACGACCACGGTGAGCATCTCGGTTGGCGGGCCATAGAAGTAGCGCCAGATGCCTACCGCAACGCCCACCAGGTTGAGCAGGACAAGGAAGATGTAGGGGCGCGAGATCACCCAGTCGACATACTCTTCTTCCACCAGTCCACCTTTg                                                                                             >  3:97617/1‑149 (MQ=255)                                          acCACGGTGAGCATCTCGGTTGGCGGGCCATAGAAGTAGCGCCAGATGCCTACCGCAACGCCCACCAGGTTGAGCAGGACAAGGAAGATGTAGGGGCGCGAGATCACCCAGTCGACGTACTCTTCTTCCACCAGTCCACCTTTggtggt                                                                                        <  1:44764/149‑1 (MQ=255)                                            cACGGTGAGCATCTCGGTTGGCGGGCCATAGAAGTAGCGCCAGATGCCTACCGCAACGCCCACCAGGTTGAGCAGGACAAGGAAGATGTAGGGGCGCGAGATCACCCAGTCGACGTACTCTTCTTCCACCAGTCCACCTTTGGCGGTGa                                                                                      <  6:149725/149‑1 (MQ=255)                                                     cATCTCGGTTGGCGGGCCATAGAAGTAGCGCCAGATGCCTACCGCAACGCCCACCAGGTTGAGCAGGACAAGGAAGATGTAGGGGCGCGAGATCACCCAGTCGACGTACTCTTCTTCCACCAGTCCACCTTTGGCGGTGACGTTAAAtt                                                                             <  2:44806/149‑1 (MQ=255)                                                           ggTTGGCGGGCCATAGAAGTAGCGCCAGATGCCTACCGCAACGCCCACCAGGTTGAGCAGGACAAGGAAGATGTAGGGGCGCGAGATCACCCAGTCGACGTACTCTTCTTCCACCAGTCCACCTTTGGCGGTGACGTTAAATTTGc                                                                          <  5:1713/146‑1 (MQ=255)                                                                   ggCCATAGAAGTAGCGCCAGATGCCTACCGCAACGCCCACCAGGTTGAGCAGGACAAGGAAGATGTAGGGGCGCGAGATCACCCAGTCGACTTACTCTTCTTCCACCAGTCCACCTTTGGCGGTGACGTTAAATTTGCCTTTGTGCggg                                                               <  4:97617/149‑1 (MQ=255)                                                                            aGTAGCGCCAGATGCCTACCGCAACGCCCACCAGGTTGAGCAGGACAAGGAAGATGTAGGGGCGCGAGATCACCCAGTCGACGTACTCTTCTTCCACCAGTCCAc                                                                                                  <  2:130456/105‑1 (MQ=255)                                                                                         gCCTACCGCAACGCCCACCAGGTTGAGCAGGACAAGGAAGATGTAGGGGCGCGAGATCACCCAGTCGACGTACTCTTCTTCCACCAGTCCACCTTTGGCGGTGACGTTAAATTTGCCTTTGTGCGGGTTAATCAGCGccacca                                               >  6:29571/1‑143 (MQ=255)                                                                                         gCCTACCGCAACGCCCACCAGGTTGAGCAGGACAAGGAAGATGTAGGGGCGCGAGATCACCCAGTCGACGTACTCTTCTTCCACCAGTCCACCTTTGGCGGTGACGTTAAATTTGCCTTTGTGCGGGTTAATCAGCGccacca                                               <  5:29571/143‑1 (MQ=255)                                                                                            tACCGCAACGCCCACCAGGTTGAGCAGGACAAGGAAGATGTAGGGGCGCGAGATCACCCAGTCGACGTACTCTTCTTCCACCAGTCCACCTTTGGCGGTGACGTTAAATTTGCCTTTGTGCGGGTTAATCAGCGCCACCAGCGTCggtg                                      >  2:46795/1‑149 (MQ=255)                                                                                                gCAACGCCCACCAGGTTGAGCAGGAAAAGGAAGATGTAGGGGCTCGAGATCACCCAGTCGCCGAACTCTTCTTCCACCAGTCCACCTTTGGCGGTGACGTTAAATTTGACTTTGTGCGTGTTAATCAGCGCCACCAGCGACGGTGGTGc                                  >  4:13554/1‑149 (MQ=255)                                                                                                 cAACGCCCACCAGGTTGAGCAGGACAAGGAAGATGTAGGGGCGCGAGATCACCCAGTCGACGTACTCTTCTTCCACCAGTCCACCTTTGGCGGTGACGTTAAATTTGCCTTTGTGCGGGTTAATCAGCGCCACCAGCGTCggtg                                      >  4:20691/1‑144 (MQ=255)                                                                                                 cAACGCCCACCAGGTTGAGCAGGACAAGGAAGATGTAGGGGCGCGAGATCACCCAGTCGACGTACTCTTCTTCCACCAGTCCACCTTTGGCGGTGACGTTAAATTTGCCTTTGTGCGGGTTAATCAGCGCCACCAGCGTCggtg                                      <  3:20691/144‑1 (MQ=255)                                                                                                   aCGCCCACCAGGTTGAGCAGGACAAGGAAGATGTAGGGGCGCGAGATCACCCAGTCGACGTACTCTTCTTCCACCAGTCCACCTTTGGCGGTGACGTTAAATTTGCCTTTGTGCGGGTTAATCAGCGCCACCAGCGTCGGTGGTGCGa                                >  6:154871/1‑148 (MQ=255)                                                                                                     gCCCACCAGGTTGAGCAGGACAAGGAAGATGTAGGGGCGCGAGATCACCCAGTCGACGTACTCTTCTTCCACCAGTCCACCTTTGGCGGTGACGTTAAATTTGCCTTTGTGCGGGTTAATCAGCGCCACCAGCGTCGGTGGTGCGatat                             <  1:46795/149‑1 (MQ=255)                                                                                                       ccaccaGGTTGAGCAGGACAAGGAAGATGTAGGGGCGCGAGATCACCCAGTCGACGTACTCTTCTTCCACCAGTCCACCTTTGGc                                                                                           <  4:136849/85‑1 (MQ=255)                                                                                                       ccaccaGGTTGAGCAGGACAAGGAAGATGTAGGGGCGCGAGATCACCCAGTCGACGTACTCTTCTTCCACCAGTCCACCTTTGGc                                                                                           >  3:136849/1‑85 (MQ=255)                                                                                                                     aGGACAAGGAAGATGTAGGGGCGCGAGATCACCCAGTCGACGTACTCTTCTTCCACCAGTCCACCTTTGGCGGTGACGTTAAATTTGCCTTTGTGCGGGTTAATCAGCGCCACCAGCGTCGGTGGTGCGATATACCACGCCAGCACCGt             >  2:53473/1‑149 (MQ=255)                                                                                                                          aaGGAAGATGTAGGGGCGCGAGATCACCCAGTCGACGTACTCTTCTTCCACCAGTCCACCTTTGGCGGTGACGTTAAATTTGCCTTTGTGCGGGTTAATCAGCGCCACCAGCGTCGGTGGTGCGATATACCACGCCAGCACCGTTTCg         >  1:76584/1‑148 (MQ=255)                                                                                                                                gATGTAGGGGCGCGAGATCACCCAGTCGACGTACTCTTCTTCCACCAGTCCACCTTTGGCGGTGACGTTAAATTTGCCTTTGTGCGGGTTAATCAGCGCCACCAGCGTCGGTGGTGCGATATACCACGCCAGCACCGTTTCGTAGAttt  <  5:154871/149‑1 (MQ=255)                                                                                                                                         |                                                                                                                                            CAAGAACAATCAGGTTGTAGAACACCCACACCATACTGACGACCACGGTGAGCATCTCGGTTGGCGGGCCATAGAAGTAGCGCCAGATGCCTACCGCAACGCCCACCAGGTTGAGCAGGACAAGGAAGATGTAGGGCCGCGAGATCACCCAGTCGACGTACTCTTCTTCCACCAGTCCACCTTTGGCGGTGACGTTAAATTTGCCTTTGTGCGGGTTAATCAGCGCCACCAGCGTCGGTGGTGCGATATACCACGCCAGCACCGTTTCGTAGATTT  >  NC_000913/3693190‑3693465

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