Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,689,231	T→G	D344A (GAT→GCT)	bcsZ ←	endo‑1,4‑D‑glucanase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,689,231	0	T	G	100.0%	95.5 / NA	27	D344A (GAT→GCT)	bcsZ	endo‑1,4‑D‑glucanase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base G (15/12);  total (15/12)

CGGTTGGTGCTGCCTCGACCATCGGCATGACGGCCAGACCGAGGGAAAGCGTGAATATGTTTAGTGTGAATTTGCGCATTCCTGGCCCCAGTCAGGTAATAACTCACCTTTTGTCGAGAAGCGGAAACGGTGTTGATCCCAGCCTTGTCCAAACAGGGTCAGCACATAGTTGTAATAGGCATCGCTGCCGGGAAAGTTATCGGCCACGCGCTGGCGCTGAACGGCCTGCGCATCGCGGTTTTGTAAAAAGGGCAGCATGGCGGCAGAAAAAC  >  NC_000913/3689095‑3689366                                                                                                                                         |                                                                                                                                        cGGTTGGTGCTGCCTCGACCATCGGCATGACGGCCAGACCGAGGGAAAGCGTGAATATGTTTAGTGTGAATTTGCGCATTCCTGGCCCCAGTCAGGTAATAACTCACCTTTTGTCGAGAAGCGGAAACGGTGTTGAGccc                                                                                                                                      >  2:693023/1‑140 (MQ=255) cGGTTGGTGCTGCCTCGACCATCGGCATGACGGCCAGACCGAGGGAAAGCGTGAATATGTTTAGTGTGAATTTGCGCATTCCTGGCCCCAGTCAGGTAATAACTCACCTTTTGTCGAGAAGCGGAAACGGTGTTGAGccc                                                                                                                                      >  2:149271/1‑140 (MQ=255)               tCGACCATCGGCATGACGGCCAGACCGAGGGAAAGCGTGAATATGTTTAGTGTGAATTTGCGCATTCCTGGCCCCAGTCAGGTAATAACTCACCTTTTGTCGAGAAGCGGAAACGGTGTTGAGCCCAGCCTTGTCCAAAc                                                                                                                        <  1:12350/140‑1 (MQ=255)                 gACCATCGGCATGACGGCCAGACCGAGGGAAAGCGTGAATATGTTTAGTGTGAATTTGCGCATTCCTGGCCCCAGTCAGGTAATAACTCACCTTTTGTCGAGAAGCGGAAACGGTGTTGAGCCCAGCCTTGTCCAAACAg                                                                                                                      <  2:388651/140‑1 (MQ=255)                  aCCATCGGCATGACGGCCAGACCGAGGGAAAGCGTGAATATGTTTAGTGTGAATTTGCGCATTCCTGGCCCCAGTCAGGTAATAACTCACCTTTTGTCGAGAAGCGGAAACGGTGTTGAGCCCAGCCTTGTCCAAACAgg                                                                                                                     >  1:13975/1‑140 (MQ=255)                  aCCATCGGCATGACGGCCAGACCGAGGGAAAGCGTGAATATGTTTAGTGTGAATTTGCGCATTCCTGGCCCCAGTCAGGTAATAACTCACCTTTTGTCGAGAAGCGGAAACGGTGTTGAGCCCAGCCTTGTCCAAACAgg                                                                                                                     >  2:246692/1‑140 (MQ=255)                         gCATGACGGCCAGACCGAGGGAAAGCGTGAATATGTTTAGTGTGAATTTGCGCATTCCTGGCCCCAGTCAGGTAATAACTCACCTTTTGTCGAGAAGCGGAAACGGTGTTGAGCCCAGCCTTGTCCAAACAGGGTCAGca                                                                                                              <  1:387822/140‑1 (MQ=255)                         gCATGACGGCCAGACCGAGGGAAAGCGTGAATATGTTTAGTGTGAATTTGCGCATTCCTGGCCCCAGTCAGGTAATAACTCACCTTTTGTCGAGAAGCGGAAACGGTGTTGAGCCCAGCCTTGTCCAAACAGGGTCAGca                                                                                                              <  2:538968/140‑1 (MQ=255)                           aTGACGGCCAGACCGAGGGAAAGCGTGAATATGTTTAGTGTGAATTTGCGCATTCCTGGCCCCAGTCAGGTAATAACTCACCTTTTGTCGAGAAGCGGAAACGGTGTTGAGCCCAGCCTTGTCCAAACAGGGTCAGcaca                                                                                                            <  1:633279/140‑1 (MQ=255)                            tGACGGCCAGACCGAGGGAAAGCGTGAATATGTTTAGTGTGAATTTGCGCATTCCTGGCCCCAGTCAGGTAATAACTCACCTTTTGTCGAGAAGCGGAAACGGTGTTGAGCCCAGCCTTGTCCAAACAGGGTCAGCACAt                                                                                                           <  1:142585/140‑1 (MQ=255)                              aCGGCCAGACCGAGGGAAAGCGTGAATATGTTTAGTGTGAATTTGCGCATTCCTGGCCCCAGTCAGGTAATAACTCACCTTTTGTCGAGAAGCGGAAACGGTGTTGAGCCCAGCCTTGTCCAAACAGGGTCAGCACATAg                                                                                                         <  1:161679/140‑1 (MQ=255)                                                  cGTGAATATGTTTAGTGTGAATTTGCGCATTCCTGGCCCCAGTCAGGTAATAACTCACCTTTTGTCGAGAAGCGGAAACGGTGTTGAGCCCAGCCTTGTCCAAACAGGGTCAGCACATAGTTGTAATAGGCATCGCTGcc                                                                                     <  1:149271/140‑1 (MQ=255)                                                                gtgtGAATTTGCGCATTCCTGGCCCCAGTCAGGTAATAACTCACCTTTTGTCGAGAAGCGGAAACGGTGTTGAGCCCAGCCTTGTCCAAACAGGGTCAGCACATAGTTGTAATAGGCATCGCTGCCGGGAAAGTTATCgg                                                                       >  1:250951/1‑140 (MQ=255)                                                                   tgAATTTGCGCATTCCTGGCCCCAGTCAGGTAATAACTCACCTTTTGTCGAGAAGCGGAAACGGTGTTGAGCCCAGCCTTGTCCAAACAGGGTCAGCACATAGTTGTAATAGGCATCGCTGCCGGGAAAGTTATCGGCCa                                                                    >  1:299795/1‑140 (MQ=255)                                                                        ttGCGCATTCCTGGCCCCAGTCAGGTAATAGCTCACCTTTTGTCGAGAAGCGGAAACGGTGTTGAGCCCAGCCTTGTCCAAACAGGGTCAGCACATAGTTGTAATAGGCATCGCTGCCGGGAAAGTTATCGGCCACgcgc                                                               <  2:330631/140‑1 (MQ=255)                                                                           cgcATTCCTGGCCCCAGTCAGGTAATAACTCACCTTTTGTCGAGAAGCGGAAACGGTGTTGAGCCCAGCCTTGTCCAAACAGGGTCAGCACATAGTTGTAATAGGCATCGCTGCCGGGAAAGTTATCGGCCACgcgctgg                                                            >  1:344430/1‑140 (MQ=255)                                                                               ttCCTGGCCCCAGTCAGGTAATAACTCACCTTTTGTCGAGAAGCGGAAACGGTGTTGAGCCCAGCCTTGTCCAAACAGGGTCAGCACATAGTTGTAATAGGCATCGCTGCCGGGAAAGTTATCGGCCACgcgctggcgct                                                        >  2:368341/1‑140 (MQ=255)                                                                                 ccTGGCCCCAGTCAGGTAATAACTCACCTTTTGTCGAGAAGCGGAAACGGTGTTGAGCCCAGCCTTGTCCAAACAGGGTCAGCACATAGTTGTAATAGGCATCGCTGCCGGGAAAGTTATCGGTCACGCGCTGGCGCTGa                                                      <  2:418381/140‑1 (MQ=255)                                                                                          aGTCAGGTAATAACTCACCTTTTGTCGAGAAGCGGAAACGGTGTTGAGCCCAGCCTTGTCCAAACAGGGTCAGCACATAGTTGTAATAGGCATCGCTGCCGGGAAAGTTATCGGCCACGCGCTGGCGCTGAACGGCCTgc                                             >  1:695388/1‑140 (MQ=255)                                                                                            tCAGGTAATAACTCACCTTTTGTCGAGAAGCGGAAACGGTGTTGAGCCCAGCCTTGTCCAAACAGGGTCAGCACATAGTTGTAATAGGCATCGCTGCCGGGAAAGTTATCGGCCACGCGCTGGCGCTGAACGGCCTgcgc                                           >  2:354261/1‑140 (MQ=255)                                                                                                          aCCTTTTGTCGAGAAGCGGAAACGGTGTTGAGCCCAGCCTTGTCCAAACAGGGTCAGCACATAGTTGTAATAGGCATCGCTGCCGGGAAAGTTATCGGCCACGCGCTGGCGCTGAACGGCCTGCGCATCGCGGTTTTGTa                             >  2:618903/1‑140 (MQ=255)                                                                                                              tttGTCGAGAAGCGGAAACGGTGTTGAGCCCAGCCTTGTCCAAACAGGGTCAGCACATAGTTGTAATAGGCATCGCTGCCGGGAAAGTTATCGGCCACGCGCTGGCGCTGAACGGCCTGCGCATCGCGGTTTTGTaaaaa                         >  2:569624/1‑140 (MQ=255)                                                                                                                      gaAGCGGAAACGGTGTTGAGCCCAGCCTTGTCCAAACAGGGTCAGCACATAGTTGTAATAGGCATCGCTGCCGGGAAAGTTATCGGCCACGCGCTGGCGCTGAACGGCCTGCGCATCGCGGTTTTGTAAAAAGGgcagca                 >  2:527434/1‑140 (MQ=255)                                                                                                                            gAAACGGTGTTGAGCCCAGCCTTGTCCAAACAGGGTCAGCACATAGTTGTAATAGGCATCGCTGCCGGGAAAGTTATCGGCCACGCGCTGGCGCTGAACGGCCTGCGCATCGCGGTTTTGTAAAAAGGGCAGCATggcgg           <  2:13975/140‑1 (MQ=255)                                                                                                                            gAAACGGTGTTGAGCCCAGCCTTGTCCAAACAGGGTCAGCACATAGTTGTAATAGGCATCGCTGCCGGGAAAGTTATCGGCCACGCGCTGGCGCTGAACGGCCTGCGCATCGCGGTTTTGTAAAAAGGGCAGCATggcgg           >  1:507569/1‑140 (MQ=255)                                                                                                                                    gtTGAGCCCAGCCTTGTCCAAACAGGGTCAGCACATAGTTGTAATAGGCATCGCTGCCGGGAAAGTTATCGGCCACGCGCTGGCGCTGAACGGCCTGCGCATCGCGGTTTTGTAAAAAGGGCAGCATGGCGGCAGaaaaa   <  1:415675/140‑1 (MQ=255)                                                                                                                                     ttGAGCCCAGCCTTGTCCAAACAGGGTCAGCACATAGTTGTAATAGGCATCGCTGCCGGGAAAGTTATCGGCCACGCGCTGGCGCTGAACGGCCTGCGCATCGCGGTTTTGTAAAAAGGGCAGCATGGCGGCAGAAAAAc  >  2:368226/1‑140 (MQ=255)                                                                                                                                         |                                                                                                                                        CGGTTGGTGCTGCCTCGACCATCGGCATGACGGCCAGACCGAGGGAAAGCGTGAATATGTTTAGTGTGAATTTGCGCATTCCTGGCCCCAGTCAGGTAATAACTCACCTTTTGTCGAGAAGCGGAAACGGTGTTGATCCCAGCCTTGTCCAAACAGGGTCAGCACATAGTTGTAATAGGCATCGCTGCCGGGAAAGTTATCGGCCACGCGCTGGCGCTGAACGGCCTGCGCATCGCGGTTTTGTAAAAAGGGCAGCATGGCGGCAGAAAAAC  >  NC_000913/3689095‑3689366

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