Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,902,329	(C)5→4	coding (1373/1413 nt)	bglB ←	6‑phospho‑beta‑glucosidase B

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,902,329	0	C	.	100.0%	97.8 / NA	24	coding (1373/1413 nt)	bglB	6‑phospho‑beta‑glucosidase B
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base . (12/12);  total (12/12)

GCTTTGCGTTTCTCTTAACGCCTTTTCTAATAACTCAAGGCGTTGTTCCACCGTTAATGATTCAGCCAATGATTGTGCAGAAAAGGCTAACGGTGCCATTAATAAGATGGCAGAGGTAATAAGATTTCGTCTAAACATAATATCCCTTTATGGTGCAAAGAAAGAATTAACGCATCGCATCCAGACTGTTCTGAATGCGACGATAATTAAGGTGCTTTAATGGTTATTTTTTTTAATGACAGCCCCCGCGTCTTGATCACCTCTGCATACCATCCGAAGCTTTTCTTGCGTGTTCTTGTGAGGCTTCCTTCGCCATTATCATCAC‑GATCCACATAAATAAAGCCGTAGCGCTTAGACATTTGTGAATGAGATGCACTGACTAAATCAATTGGCCCCCAACTGGTGTACCCCATAATATCCACACCATCGGCAATCGCTTCATTTACCTGTACCAGGTGATCGTTTAAATAGGCAATTCGATAATCGTCCTGTA  >  NC_000913/3902083‑3902575                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                        gCTTTGCGTTTCTCTTAACGCCTTTTCTAATAACTCAAGGCGTTGTTCCACCGTTAATGATTCAGCCAATGATTGTGCAGAAAAGGCTAACGGTGCCATTAATAAGATGGCAGAGGTAATAAGATTTCGTCTAAACATAATATCCCTTTATGGTGCAAAGAAAGAATTAACGCATCGCATCCAGACTGTTCTGAATGCGACGATAATTAAGGTGCTTTAATGGTTATTTTTTTTAATGACAGCCCC‑GCGTc                                                                                                                                                                                                                                                    >  2:8362/1‑251 (MQ=255)          ttCTCTTAACGCCTTTTCTAATAACTCAAGGCGTTGTTCCACCGTTAATGATTCAGCCAATGATTGTGCAGAAAAGGCTAACGGTGCCATTAATAAGATGGCAGAGGTAATAAGATTTCGTCTAAACATAATATCCCTTTATGGTGCAAAGAAAGAATTAACGCATCGCATCCAGACTGTTCTGAATGCGACGATAATTAAGGTGCTTTAATGGTTATTTTTTTTAATGACAGCCCC‑GCGTCTTGATCAc                                                                                                                                                                                                                                            <  1:256960/250‑1 (MQ=255)                     ccTTTTCTAATAACTCAAGGCGTTGTTCCACCGTTAATGATTCAGCCAATGATTGTGCAGAAAAGGCTAACGGTGCCATTAATAAGATGGCAGAGGTAATAAGATTTCGTCTAAACATAATATCCCTTTATGGTGCAAAGAAAGAATTAACGCATCGCATCCAGACTGTTCTGAATGCGACGATAATTAAGGTGCTTTAATGGTTATTTTTTTTAATGACAGCCCC‑GCGTCTTGATCACCTCTGCATACCa                                                                                                                                                                                                                                <  1:163235/251‑1 (MQ=255)                         ttCTAATAACTCAAGGCGGTGTTCCACCGTTAATGATTCAGCCAATGATTGTGCAGAAAAGGCTACCGATGCCATTAATAAGATGGCAGAGGTAATAAGATTTCGTCTAAACATAATATCCCTTTATTGTGCAAAGAAAGAATTAACGCATCGCATCCAGACTGTTCTGAATGCGACGATAATTAAGGTGCTTTAATGGTTATTTTTTTTAATGACAGCCCC‑GCGTCTTGATCACCTCTGCATACCATcc                                                                                                                                                                                                                             <  1:175260/250‑1 (MQ=255)                                 aCTCAAGGCGTTGTTACACCGTTAATGATTCAGCCAATGATTGTGCAGAAAAGGCTAACGGTGCCATTAATAAGATGGCAGAGGTAATAAGATTTCGTCTACACATAATATCCCTTTATGGTGCAAAGAAAGAATTAACGCATCGCATCCAGACTGTTCTGAATGCGACGATAATTAAGGTGCTTTAATGGTTATTTTTTTTAATGACAGCCCA‑GCGTCTTGATCACCTCTGCATACCATCCGAAGCtttt                                                                                                                                                                                                                    >  2:150665/1‑251 (MQ=255)                                  cTCAAGGCGTTGTTCCACCGTTAATGATTCAGCCAATGATTGTGCAGAAAAGGCTAACGGTGCCATTAATAAGATGGCAGAGGTAATAAGATTTCGTCTAAACATAATATCCCTTTATGGTGCAAAGAAAGAATTAACGCATCGCATCCAGACTGTTCTGAATGCGACGATAATTAAGGTGCTTTAATGGTTATTTTTTTTAATGACAGCCCC‑GCGTCTTGATCACCTCTGCATACCATCCGAAGCTTTTc                                                                                                                                                                                                                   >  1:315706/1‑251 (MQ=255)                                              ttCCACCGTTAATGATTCAGCCAATGATTGTGCAGAAACGCCTAACGGTGCCATTAATAAGATGGCAGAGGTAATAAGATTTCGTCTAAACATAATATCCCTTTATGGTGCAAAGAAAGAATTAACGCATCGCATCCAGACTGTTCTGAATGCGACGATAATTAAGGTGCTTTAATGGTTATTTTTTTTAATGACAGCCCC‑GCGTCTTGATCACCTCTGCATACCATCCGAAGCTTTTCTTGCGTGTTCtt                                                                                                                                                                                                       >  1:34064/1‑251 (MQ=255)                                                        aaTGATTCAGCCAATGATTGTGCAGAAAAGGCTAACGGTGCCATTAATAAGATGGCAGAGGTAATAAGATTTCGTCTAAACATAATATCCCTTTATGGTGCAAAGAAAGAATTAACGCATCGCATCCAGACTGTTCTGAATGCGACGATAATTAAGGTGCTTTAATGGTTATTTTTTTTAATGACAGCCCC‑GCGTCTTGATCACCTCTGCATACCATCCGAAGCTTTTCTTGCGTGTTCTTGTGAGGcttc                                                                                                                                                                                             <  1:168567/251‑1 (MQ=255)                                                              tCAGCCAATGATTGTGCAGAAAAGGCTAACGGTGCCATTAATAAGATGGCAGAGGTAATAAGATTTCGTCTAAACATAATATCCCTTTATGGTGCAAAGAAAGAATTAACGCATCGCATCCAGACTGTTCTGAATGCGACGATAATTAAGGTGCTTTAATGGTTATTTTTTTTAATGACAGCCCC‑GCGTCTTGATCACCTCTGCATACCATCCGAAGCTTTTCTTGCGTGTTCTTGTGAGGCTTCCTTCGc                                                                                                                                                                                       >  1:128068/1‑251 (MQ=255)                                                                                         aaCGGTGCCATTAATAAGATGGCAGAGGGCATAAGATTTCCTCTAAACATAATATCCCTTTATGGTGCAAAGAAAGCATAAACGCATCGCATCCAGACTGTTCTGAATGCGACGATAATTAAGGTGCTTTAATGGTTATTTTTTTTAATGAAAGCCCC‑GAGTCTTGATCACCTCTGCATACCATCCGAAGCTTTTCTTGCGTGTTCTTGTGAGGCTTCCTTCGCCATTATCATCAC‑GATCCACataaata                                                                                                                                                            >  2:18131/1‑250 (MQ=255)                                                                                                    taataaGATGGCAGAGGTAATAAGATTTCGTCTAAACATAATATCCCTTTATGGTGCAAAGAAAGAATTAACGCATCGCATCCAGACTGTTCTGAATGCGACGATAATTAAGGAGCTTTAATGGTTATTTTTTTTAATGACAGCCCC‑GCGTCTTGATCACCTCTGCATACCATCCGAAGCTTTTCTTGCGTGTTCTTGTGAGGCTTCCTTCGCCATTATCATCAC‑GATCCACATAAATAAAGCCGTAgcgc                                                                                                                                                >  1:89904/1‑251 (MQ=255)                                                                                                      ataaGATGGCAGAGGTAATAAGATTTCGTCTAAACATAATATCCCTTTATGGTGCAAAGAAAGAATTAACGCATCGCATCCAGACTGTTCTGAATGCGACGATAATTAAGGTGCTTTAATGGTTATTTTTTTTAATGACAGCCCC‑GCGTCTTGATCACCTCTGCATACCATCCGAAGCTTTTCTTGCGTGTTCTTGTGAGGCTTCCTTCGCCATTATCATCAC‑GATCCACATAAATAAAGCCGTAGCGCtt                                                                                                                                              <  1:269452/251‑1 (MQ=255)                                                                                                             ggCAGAGGTAATAAGATTTCGTCTAAACATAATATCCCTTTATGGTGCAAAGAAAGAATTAACGCATCGCATCCAGACTGTTCTGAATGCGACGATAATTAAGGTGCTTTAATGGTTATTTTTTTTAATGACAGCCCC‑GCGTCTTGATCACCTCTGCATACCATCCGAAGCGTTTCTTGCGTGTTCTTGTGAGGCGTCCTTCGCCATTATCATCAC‑TATCCACATAAATAAAGCCGTAGCGCTTAGACAt                                                                                                                                        <  1:259176/250‑1 (MQ=255)                                                                                                                            aTTTCGTATAAACATAATATCCCTTTATGGTGCAAAGAAAGAATTACCGCATCGCATCCAGACTGTTCTGAATGCGACGATAATTAAGGTGCTTTAATGGTTATTTTTTTTAATGACAGCCCC‑GCGTCTTGATCAACTCTGCATACAATCCGAAGCTTTTCTTGCGTGTTCTTGTGAGGATTCCTTCGCCATTATCAT‑ACAGATCCACATAAATAAAGCCGTAGCGATTAGACATTTGTGAATGAGATGc                                                                                                                         >  1:185875/1‑250 (MQ=255)                                                                                                                                                      aTGGTGCAAAGAAAGAAGTAACGCATCGCATCCAGACTTTTCTGAATGCGACGATAATTAAGGTGCTTTAATGGTTATTTTTTTTAATGACAGCCCC‑GCGTCTTGATCACCTCTGCATACCATCCGAAGCTTTTCTTGCGTGTTCTTGTGAGGCTTCCTTCGCCATTATCATCAC‑GATCCACATAAATAAAGCCGTAGCTCTTAGACATTTGTGAATGAGATGCACTGACTAAATCAATTGGCCCCCAAc                                                                                               <  2:122981/250‑1 (MQ=255)                                                                                                                                                                    gaaTTAACGCATCGCATCCAGACTGTTCTGAATGCGACGATAATTAAGGTGCTTTAATGGTTATTTTTTTTAATGACAGCCCC‑GCGTCTTGATCACCTCTGCATACCATCCGAAGCTTTTCTTGCGTGTTCTTGTGAGGCTTCCTTCGCCATTATCATCAC‑GATCCACATAAATAAAGCCGTAGCGCTTAGACATTTGTGAATGAGATGCACTGACTAAATCAATTGGCCCCCAACTGGTGTACCCCataa                                                                                <  1:8362/251‑1 (MQ=255)                                                                                                                                                                                   atcCAGACTGTTCTGAATGCGACGATAATTAAGGAGCTTTAATGGTTATTTTTTTTAATGACAGCCCC‑GCGTCTTGATCACCTCTGCATACCATCCGAAGCTTTTCTTGCGTGTTCTTGTGAGGCTTCCTTCGCCATTATCATCAC‑GATCCACATAAATAAAGCCGTAGCGCTTAGACATTTGTGAATGAGATGCACTGACTAAATCAATTGGCCCCCAACTGGTGTACCCCATAATATCCACACCATCgg                                                                 >  1:212555/1‑251 (MQ=255)                                                                                                                                                                                   atcCAGACTGCTATGAATGCGCCGATAATTAAGGTGCTTTAATGGTTATTTTTTTTAATGACAGCCCC‑GCGTCTTGATCACCTCTGCATACCATCCGAAGCTTTTCTTGCGTGTTCTTGTGAGGCTTCCTTCGCCATTATCATCAC‑GATCCACATAAATAAAGCCGTAGCGCTTAGACATTTGTGAATGAGATGCACTGACTAAATCAATTGGCCCCCAACTGTTGTACCCCATAATATCCACACCATCgg                                                                 >  2:60205/1‑251 (MQ=255)                                                                                                                                                                                         aCTGTTCTGAATGCGACGATAATTAAGGTGCTTTAATGGTTATTTTTTTTAATGACAGCCCC‑GCGTCTTGATCACCTCTGCATACCATCCGAAGCTTTTCTTGCGTGTTCTTGTGAGGCTTCCTTCGCCATTATCATCAC‑GAGCCACATAAATAAAGCCGTAGCGCTTAGACATTTGTGAATGAGATGCACTGACTAAATCAATTGGCCCCCAACTGGTGTACCCCATAATATCCACACCATCGGCAATCg                                                           <  2:112190/251‑1 (MQ=255)                                                                                                                                                                                                                         ttAATGGTTATTTTTTTTAATGACAGCCCC‑GCGTCTTGATCACCTCTGCATACCATCCGAAGCTTTTCTTGCGTGTTCTTGTGAGGCTTCCTTCGCCATTATCATCAC‑GATCCACATAAATAAAGCCGTAGCGCTTAGACATTTGTGAATGAGATGCACTGACTAAATCAATTGGCCCCCAACTGGTGTACCCCATAATATCCACACCATCGGCAATCGCTTCATTTACCTGTACCAGGTGATCGTTTa                             <  1:150665/249‑1 (MQ=255)                                                                                                                                                                                                                           gaTGGTTATTTTTTTTAATGACAGCCCC‑GCGTCTTGATCACCTCTGCATACCATCCGAAGCTTTTCTTGCGTGTTCTTGTGAGGCTTCCTTCGCCATTATCATCAC‑GATCCACATAAATAAAGCCGTAGCGCTTAGACATTTGTGAATGAGATGCACTGACTAAATCAATTGGCCCCCAACTGGTGTACCCCATAATATCCACACCATCGGCAATCGCTTCATTTACCTGTACAAGGTGATCGTTTAAAt                          <  1:246640/249‑1 (MQ=255)                                                                                                                                                                                                                           aaTGGTTCTTTTTTTTAATGACCGCCCC‑GCGTCTTGATCACCTCTGCATACCATCCGAAGCTTTTCTTGCGTGTTCTTGTGAGGCTTCCTTCGCCATTATCCTCAC‑GATCCACATAACTAAAGCCGTAGCGCTTAGACATTTGTGACTGAGATGCACTGACTAAATCAATTGGCCCCCAACTGGTGTACCCCATAATATCCACACCATCGGCAATCGCTTCATTTACCTGTACCAGGTGATAGTTTaact                          >  2:198291/1‑248 (MQ=255)                                                                                                                                                                                                                                     ttttttAATGACAGCCCC‑GCGTCTTGATCACCTCTGCATACCATCCGAAGCTTTTCTTGCGTGTTCTTGTGAGGCTTCCTTCGCCATTATCATCAC‑GATCCACATAAATAAAGCCGTAGCGCTTAGTCATTTGTGAATGAGATGCACTGACTAAATCAATTGGCCCCCAACTGGTGTACCCCATAATATCCACACCATCGGCAATCGCTTCATTTACCCGTACCAGGTGATCGTTTAAATAGGCAATTCg                <  1:312584/250‑1 (MQ=255)                                                                                                                                                                                                                                      tttttAATGACAGCCCC‑GCGTCTTGATCACCTCTGCATACCATCCGAAGCTTTTCTTGCGTGTTCTTGTGAGGCTTCCTTCGCCATTATCATCAC‑GATCCACATAAATAAAGCCGTAGCGCTTAGACATTTGTGAATGAGATGCACTGACTAAATCAATTGGCCCCCAACTGGTGTACCCCATAATATCCACACCATCGGCAATCGCTTCATTTACCTGTACCAGGTGATCGTTTAAATAGGCAATTCGAt              >  1:293798/1‑251 (MQ=255)                                                                                                                                                                                                                                                   cccccGCGTCTTGATCACCTCTGCATACCATCCGAAGCTTTTCTTGCGTGTTCTTGTGAGGCTTCCTTCGCCATTATCATCAC‑GATCCACATAAATAAAGCCGTAGCGCTTAGACATTGGTGAATTAGATGCACTGACTAAATCAATTGTCCCCCAACTGGTGTACCCCATAATATCCTCACCATCGGCAATCGCTTCATTTACCTGTACCAGGTGATCGTTTAAATAGGCAATTCGATAATCGTCCTGTa  <  2:108208/251‑1 (MQ=255)                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                        GCTTTGCGTTTCTCTTAACGCCTTTTCTAATAACTCAAGGCGTTGTTCCACCGTTAATGATTCAGCCAATGATTGTGCAGAAAAGGCTAACGGTGCCATTAATAAGATGGCAGAGGTAATAAGATTTCGTCTAAACATAATATCCCTTTATGGTGCAAAGAAAGAATTAACGCATCGCATCCAGACTGTTCTGAATGCGACGATAATTAAGGTGCTTTAATGGTTATTTTTTTTAATGACAGCCCCCGCGTCTTGATCACCTCTGCATACCATCCGAAGCTTTTCTTGCGTGTTCTTGTGAGGCTTCCTTCGCCATTATCATCAC‑GATCCACATAAATAAAGCCGTAGCGCTTAGACATTTGTGAATGAGATGCACTGACTAAATCAATTGGCCCCCAACTGGTGTACCCCATAATATCCACACCATCGGCAATCGCTTCATTTACCTGTACCAGGTGATCGTTTAAATAGGCAATTCGATAATCGTCCTGTA  >  NC_000913/3902083‑3902575

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