Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,699,320	T→C	A369A (GCT→GCC)	bcsG →	DUF3260 family cellulose production inner membrane protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,699,320	0	T	C	100.0%	122.6 / NA	35	A369A (GCT→GCC)	bcsG	DUF3260 family cellulose production inner membrane protein
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (23/12);  total (23/12)

GATGGGACATAACGGCCAGTTCGGCGGTTTTTTGAAAGAAGTTCGCGAAAATGGCGGCATGCAGAGCGAATTGATGGATCAAACAAATCTGCCGGTTATTTTGCTGGGCTTTGATGGTTCGCCGGTTTATGACGATACCGCTGTGCTTAACCGCTGGCTGGACGTTACCGAAAAAGATAAAAACAGCCGTAGTGCCACGTTCTACAACACGCTTCCACTGCATGACGGCAACCATTATCCGGGGGTCAGCAAAACAGCGGATTACAAAGCGCGGGCGCAGAAATTCT  >  NC_000913/3699179‑3699465                                                                                                                                              |                                                                                                                                                  gatgGGACATAACGGCCAGTTCGGCGGTTTTTTGAAAGAAGTTCGCGAAAATGGCGGCATGCAGAGCGAATTGATGGATCAAACAAATCTGCCGGTTATTTTGCTGGGCTTTGATGGTTCGCCGGTTTATGACGATACCGCCGTGCTTaa                                                                                                                                           <  2:535916/150‑1 (MQ=255)   tgGGACATAACGGCCAGTCCGGCGGTTTTTTGAAAGAAGTTCGCGAAAATGGCGGCATGCAGAGCGAATTGATGGATCAACCAAATCTGCCGGTTATTTTGCTGGGCTTTGATGGTTCCCCGGTTTATGACGATACCGCCGTGCTTAAc                                                                                                                                          <  1:44728/149‑1 (MQ=255)            aCGGCCAGTTCGGCGGTTTTTTGAAAGAAGTTCGCGAAAATGGCGGCACGCAGAGCGGATTGATGGATCAAACAAATCTGCCGGTTATTTTGCTGGGCTTTGATGGTTCGCCGGTTTATGACGATACCGCCGTGCTTAACCGctggctgg                                                                                                                                >  1:36609/1‑150 (MQ=255)              ggCCAGTTCGGCGGTTTTTTGAAAGAAGTTCGCGAAAATGGCGGCATGCAGAGCGAATTGATGGATCAAACAAATCTGCCGGTTATTTTGCTGGGCTTTGATGGTTCGCCGGTTTATGACGATACCGCCGTGCTTAACCGCTGGCTGGAc                                                                                                                              <  1:602615/150‑1 (MQ=255)                     tCGGCGGTTTTTTGAAAGAAGTTCGCGAAAATGGCGGCATGCAGAGCGAATTGATGGATCAAACAAATCTGCCGGTTATTTTGCTGGGCTTTGATGGTTCGCCGGTTTATGACGATACCGCCGTGCTTAACCGCTGGCTGGACGTTACCg                                                                                                                       >  1:450565/1‑150 (MQ=255)                                  gAAAGAAGTTCGCGAAAATGGCGGCATGCAGAGCGAATTGATGGATCAAACAAATCTGCCGGTTATTTTGCTGGGCTTTGATGGTTCGCCGGTTTATGACGATACCGCCGTGCTTAACCGCTGGCTGGACGTTACCGAAAAAGATaaaaa                                                                                                          <  2:456271/150‑1 (MQ=255)                                   aaaGAAGTTCGCGAAAATGGCGGCATGCAGAGCGAATTGATGGATCAAACAAATCTGCCGGTTATTTTGCTGGGCTTTGATGGTTCGCCGGTTTATGACGATACCGCCGTGCTTAACCGCTGGCTGGACGTTACCGAAAAAGATAAAAAc                                                                                                         >  2:279094/1‑150 (MQ=255)                                   aaaGAAGTTCGCGAAAATGGCGGCATGCAGAGCGAATTGATGGATCAAACAAATCTGCCGGTTATTTTGCTGGGCTTTGATGGTTCGCCGGTTTATGACGATACCGCCGTGCTTAACCGCTGGCTGGACGTTACCGAAAAAGATAAAAACa                                                                                                        <  2:275506/151‑1 (MQ=255)                                     agaagTTCGCGAAAATGGCGGCATGCAGAGCGAATTGATGGATCAAACAAATCTGCCGGTTATTTTGCTGGGCTTTGATGGTTCGCCGGTTTATGACGATACCGCCGTGCTTAACCGCTGGCTGGACGTTACCGAAAAAGATAAAAACAGc                                                                                                      >  1:694938/1‑151 (MQ=255)                                        agTTCGCGAAAATGGCGGCATGCAGAGCGAATTGATGGATCAAACAAATCTGCCGGTTATTTTGCTGGGCTTTGATGGTTCGCCGGTTTATGACGATACCGCCGTGCTTAACCGCTGGCTGGACGTTACCGAAAAAGATAAAAACAGCCGt                                                                                                   >  2:12168/1‑151 (MQ=255)                                             gcgAAAATGGCGGCATGCAGAGCGAATTGATGGATCAAACAAATCTGCCGGTTATTTTGCTGGGCTTTGATGGTTCGCCGGTTTATGACGATACCGCCGTGCTTAACCGCTGGCTGGACGTTACCGAAAAAGATAAAAACAGCCGTAGTGc                                                                                              >  1:633385/1‑151 (MQ=255)                                              cgAAAATGGCGGCATGCAGAGCGAATTGATGGATCAAACAAATCTGCCGGTTATTTTGCTGGGCTTTGATGGTTCGCCGGTTTATGACGATACCGCCGTGCTTAACCGCTGGCTGGACGTTACCGAAAAAGATAAAAACAGCCGTAGTGcc                                                                                             >  1:714069/1‑151 (MQ=255)                                                     ggcggcATGCAGAGCGAATTGATGGATCAAACAAATCTGCCGGTTATTTTGCTGGGCTTTGATGGTTCGCCGGTTTATGACGATACCGCCGTGCTTAACCGCTGGCTGGACGTTACCGAAAAAGATAAAAACAGCCGTAGTGCCACGTTCt                                                                                      >  1:122320/1‑151 (MQ=255)                                                         gcATGCAGAGCGAATTGATGGATCAAACAAATCTGCCGGTTATTTTGCTGGGCTTTGATGGTTCGCCGGTTTATGACGATACCGCCGTGCTTAACCGCTGGCTGGACGTTACCGAAAAAGATAAAAACAGCCGTAGTGCCACGTTCTaca                                                                                   >  2:604454/1‑150 (MQ=255)                                                           aTGCAGAGCGAATTGATGGATCAAACAAATCTGCCGGTTATTTTGCTGGGCTTTGATGGTTCGCCGGTTTATGACGATACCGCCGTGCTTAACCGCTGGCTGGACGTTACCGAAAAAGATAAAAACAGCCGTAGTGCCACGTTCTacaac                                                                                 >  2:484048/1‑150 (MQ=255)                                                                gagCGAATTGATGGATCAAACAAATCTGCCGGTTATTTTGCTGGGCTTTGATGGTTCGCCGGTTTATGACGATACCGCCGTGCTTAACCGCTGGCTGGACGTTACCGAAAAAGATAAAAACAGCCGTAGTGCCACGTTCTACAACACGCtt                                                                           >  1:635583/1‑151 (MQ=255)                                                                  gCGAATTGATGGATCAAACAAATCTGCCGGTTATTTTGCTGGGCTTTGATGGTTCGCCGGTTTATGACGATACCGCCGTGCTTAACCGCTGGCTGGACGTTACCGAAAAAGATAAAAACAGCCGTAGTGCCACGTTCTACAACACGCTTcc                                                                         >  1:243782/1‑151 (MQ=255)                                                                      aTTGATGGATCAAACAAATCTGCCGGTTATTTTGCTGGGCTTTGATGGTTCGCCGGTTTATGACGATACCGCCGTGCTTAACCGCTGGCTGGACGTTACCGAAAAAGATAAAAACAGCCGTAGTGCCACGTTCTACAACACGCTTCCACTg                                                                     <  1:524926/151‑1 (MQ=255)                                                                         gATGGATCAAACAAATCTGCCGGTTATTTTGCTGGGCTTTGATGGTTCGCCGGTTTATGACGATACCGCCGTGCTTAACCGCTGGCTGGACGTTACCGAAAAAGATAAAAACAGCCGTAGTGCCACGTTCTACAACACGCTTCCACTGCAt                                                                  <  2:328867/151‑1 (MQ=255)                                                                              aTCAAACAAATCTGCCGGTTATTTTGCTGGGCTTTGATGGTTCGCCGGTTTATGACGATACCGCCGTGCTTAACCGCTGGCTGGACGTTACCGAAAAAGATAAAAACAGCCGTAGTGCCACGTTCTACAACACGCTTCCACTGCATGACgg                                                             <  1:426192/151‑1 (MQ=255)                                                                                     aaaTCTGCCGGTTATTTTGCTGGGCTTTGATGGTTCGCCGGTTTATGACGATACCGCCGTGCTTAACCGCTGGCTGGACGTTACCGAAAAAGATAAAAACAGCCGTAGTGCCACGTTCTACAACACGCTTCCACTGCATGACGGCAACCAt                                                      <  2:571716/151‑1 (MQ=255)                                                                                     aaaTCTGCCGGTTATTTTGCTGGGCTTTGATGGTTCGCCGGTTTATGACGATACCGCCGTGCTTAACCGCTGGCTGGACGTTACCGAAAAAGATAAAAACAGCCGTAGTGCCACGTTCTACAACACGCTTCCACTGCATGACGGCAACCAt                                                      >  2:235506/1‑151 (MQ=255)                                                                                     aaaTCTGCCGGTTATTTTGCTGGGCTTTGATGGTTCGCCGGTTTATGACGATACCGCCGTGCTTAACCGCTGGCTGGACGTTACCGAAAAAGATAAAAACAGCCGTAGTGCCACGTTCTACAACACGCTTCCACTGCATGACGGCAACCAt                                                      <  2:412758/151‑1 (MQ=255)                                                                                         cTGCCGGTTATTTTGCTGGGCTTTGATGGTTCGCCGGTTTATGACGATACCGCCGTGCTTAACCGCTGGCTGGACGTTACCGAAAAAGATAAAAACAGCCGTAGTGCCACGTTCTACAACACGCTTCCACTGCATGACGGCAACCATTATc                                                  >  2:766565/1‑151 (MQ=255)                                                                                           gCCGGTTATTTTGCTGGGCTTTGATGGTTCGCCGGTTTATGACGATACCGCCGTGCTTAACCGCTGGCTGGACGTTACCGAAAAAGATAAAAACAGCCGTAGTGCCACGTTCTACAACACGCTTCCACTGCATGACGGCAACCATTATc                                                  <  2:46661/149‑1 (MQ=255)                                                                                              ggTTATTTTGCTGGGCTTTGATGGTTCGCCGGTTTATGACGATACCGCCGTGCTTAACCGCTGGCTGGACGTTACCGAAAAAGATAAAAACAGCCGTAGTGCCACGTTCTACAACACGCTTCCACTGCATGACGGCAACCATTATCCgggg                                             >  2:163769/1‑151 (MQ=255)                                                                                                 tATTTTGCTGGGCTTTGATGGTTCGCCGGTTTATGACGATACCGCCGTGCTTAACCGCTGGCTGGACGTTACCGAAAAAGATAAAAACAGCCGTAGTGCCACGTTCTACAACACGCTTCCACTGCATGACGGCAACCATTATCCGGGGGTc                                          >  2:152995/1‑151 (MQ=255)                                                                                                 tATTTTGCTGGGCTTTGATGGTTCGCCGGTTTATGACGATACCGCCGTGCTTAACCGCTGGCTGGACGTTACCGAAAAAGATAAAAACAGCCGTAGTGCCACGTTCTACAACACGCTTCCACTGCATGACGGCAACCATTATCCGGGGGTc                                          >  2:699103/1‑151 (MQ=255)                                                                                                             cTTTGATGGTTCGCCGGTTTATGACGATACCGCCGTGCTTAACCGCTGGCTGGACGTTACCGAAAAAGATAAAAACAGCCGTAGTGCCACGTTCTACAACACGCTTCCACTGCATGACGGCAACCATTATCCGGGGGTCAGCAAAACAGCg                              >  1:489020/1‑151 (MQ=255)                                                                                                              tttGATGGTTCGCCGGTTTATGACGATACCGCCGTGCTTAACCGCTGGCTGGACGTTACCGAAAAAGATAAAAACAGCCGTAGTGCCACGTTCTACAACACGCTTCCACTGCATGACGGCAACCATTATCCGGGGGTCAGCAAAACAGCgg                             >  2:34985/1‑151 (MQ=255)                                                                                                                  aTGGTTCGCCGGTTTATGACGATACCGCCGTGCTTAACCGCTGGCTGGACGTTACCGAAAAAGATAAAAACAGCCGTAGTGCCACGTTCTACAACACGCTTCCACTGCATGACGGCAACCATTATCCGGGGGTCAGCAAAACAGCGGAtt                          >  2:453165/1‑150 (MQ=255)                                                                                                                           cGGTTTATGACGATACCGCCGTGCTTAACCGCTGGCTGGACGTTACCGAAAAAGATAAAAACAGCCGTAGTGCCACGTTCTACAACACGCTTCCACTGCATGACGGCAACCATTATCCGGGGGTCAGCAAAACAGCGGATTACAAAgcgc                 >  1:346460/1‑150 (MQ=255)                                                                                                                           cGGTTTATGACGATACCGCCGTGCTTAACCGCTGGCTGGACGTTACCGAAAAAGATAAAAACAGCCGTAGTGCCACGTTCTACAACACGCTTCCACTGCATGACGGCAACCATTATCCGGGGGTCAGCAAAACAGCGGATTACAAAGcgcg                >  2:449582/1‑151 (MQ=255)                                                                                                                                tATGACGATACCGCCGTGCTTAACCGCTGGCTGGACGTTACCGAAAAAGATAAAAACAGCCGTAGTGCCACGTTCTACAACACGCTTCCACTGCATGACGGCAACCATTATCCGGGGGTCAGCAAAACAGCGGATTACAAAGCGCGGgcg            <  2:95477/150‑1 (MQ=255)                                                                                                                                         aCCGCCGTGCTTAACCGCTGGCTGGACGTTACCGAAAAAGATAAAAACAGCCGTAGTGCCACGTTCTACAACACGCTTCCACTGCATGACGGCAACCATTATCCGGGGGTCAGCAAAACAGCGGATTACAAAGCGCGGGCGCAGAAATTCt  >  2:407942/1‑151 (MQ=255)                                                                                                                                              |                                                                                                                                                  GATGGGACATAACGGCCAGTTCGGCGGTTTTTTGAAAGAAGTTCGCGAAAATGGCGGCATGCAGAGCGAATTGATGGATCAAACAAATCTGCCGGTTATTTTGCTGGGCTTTGATGGTTCGCCGGTTTATGACGATACCGCTGTGCTTAACCGCTGGCTGGACGTTACCGAAAAAGATAAAAACAGCCGTAGTGCCACGTTCTACAACACGCTTCCACTGCATGACGGCAACCATTATCCGGGGGTCAGCAAAACAGCGGATTACAAAGCGCGGGCGCAGAAATTCT  >  NC_000913/3699179‑3699465

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