Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,571,824	A→C	F560C (TTC→TGC)	glgB ←	1,4‑alpha‑glucan branching enzyme

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,571,824	0	A	C	100.0%	56.2 / NA	17	F560C (TTC→TGC)	glgB	1,4‑alpha‑glucan branching enzyme
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base C (7/10);  total (7/10)

CTGGACACCGTGGTGCCAGTTATCGCCGCCTTCCAACAGATGCCAGTCGAGGCTGGCGTCATGGTTCCACTCGCGGCCCTGGGCAAATTCGTTACCCATGAACAGTAGTTTCTTGCCCGGGAATGCCCACATCCAGCCATAGTAGGCGCGCAGGTTCGCGAATTTCTGCCATGCGTCGCCCGGCATGCGGTCGAGAATCGATTTTTTACCGTGGACCACTTCATCATGCGACAACGGCAGGACGAA  >  NC_000913/3571703‑3571948                                                                                                                          |                                                                                                                             cTGGACACCGTGGTGCCAGTTATCGCCGCCTTCCAACAGATGCCAGTCGAGGCTGGCGTCATGGTTCCACTCGCGGCCCTGGGCAAATTCGTTACCCATGAACAGTAGTTTCTTGCCCGGGCATGCCCACATc                                                                                                                   <  1:536202/133‑1 (MQ=255) cTGGACACCGTGGTGACATTTAGCGCCGCCCTCCAACAGATGCCAGTCGAGGCTGGCGTCATGGTTCCACTCGCGGCCCTGGGCAAATTCGTTACCCATGAACAGTAGTTTCTTGCCCGGGCATGCCCACATc                                                                                                                   <  1:270789/133‑1 (MQ=255)   ggACACCGTGCTGCCTGTTATCGCCGCCTTCCTACAGATGCCAGTCGAGGCTGGCGTCATGGTTCCACTCGCGGCCCTGGGCAAATTCGTTACCCATGAACAGTAGTTTCTTGCCCGGGCATGCCCACATCCa                                                                                                                 <  2:174288/133‑1 (MQ=255)         cGTGGTGCCAGTTATCGCCGCCTTCCAACAGATGCCAGTCGAGGCTGGCGTCATGGTTCCACTCGCGGCCCTGGGCAAATTCGTTACCCATGAACAGTAGTTTCTTGCCCGGGCATGCCCACATCCAGCCAta                                                                                                           >  2:604207/1‑133 (MQ=255)                   gTTATCGCCGCCTTCCAACAGATGCCAGTCGAGGCTGGCGTCATGGTTCCACACGCGGCCCTGGGCAAATTCGTTACCCATGAACAGTAGTTTCTTGCCCGGGCATGCCCACATCCAGCCATAGTAGgcgcgc                                                                                                 >  2:5740/1‑133 (MQ=255)                            gccTTCCAACAGATGCCAGTCGAGGCTGGCGTCATGGTTCCACTCGCGGCCCTGGGCAAATTCGTTACCCATGAACAGTAGTTTCTTGCCCGGGCATGCCCACATCCAGCCATAGTAGGCGCGCAGGTTcgcg                                                                                        <  2:272739/133‑1 (MQ=255)                              cTTCCAACAGATGCCAGTCGAGGCTGGCGTCATGGTTCCACTCGCGGCCCTGGGCAAATTCGTTACCCATGAACAGTAGTTTCTTGCCCGGGCATGCCCACATCCAGCCATAGTAGGCGCGCAGGTTCGCGaa                                                                                      >  2:20865/1‑133 (MQ=255)                                        aTGCCAGTCGAGGCTGGCGTCATGGTTCCACTCGCGGCCCTGGGCAAATTCGTTACCCATGAACAGTAGTTTCTTGCCCGGGCATGCCCACATCCAGCCATAGTAGGCGCGCAGGTTCGCGAATTTCTGCCAt                                                                            >  1:623200/1‑133 (MQ=255)                                                  aGGCTGGCGTCATGGTTCCACTCGCGGCCCTGGGCAAATTCGTTACCCATGAACAGTAGTTTCTTGCCCGGGCATGCCCACATCCAGCCATAGTAGGCGCGCAGGTTCGCGAATTTCTGCCATGCGTCACCCg                                                                  <  1:550157/133‑1 (MQ=255)                                                       ggCGTCATGGTTCCACTCGCGGCCCTGGGCAAATTCGTTACCCATGAACAGTAGTTTCTTGCCCGGGCATGCCCACATCCAGCCATAGTAGGCGCGCAGGTTCGCGAATTTCTGCCATGCGTCGCCCGGCATg                                                             <  1:586368/133‑1 (MQ=255)                                                             aTGGTTCCACTCGCGGCCCTGGGCAAATTCGTTACCCATGAACAGTAGTTTCTTGCCCGGGCATGCCCACATCCAGCCACAGTAGGCGCGCAGGTTCGCGAATTTCTGCCATGCGTCGCCCGGCATg                                                             >  1:284405/1‑127 (MQ=255)                                                             aTGGTTCCACTCGCGGCCCTGGGCAAATTCGTTACCCATGAACAGTAGTTTCTTGCCCGGGCATGCCCACATCCAGCCACAGTAGGCGCGCAGGTTCGCGAATTTCTGCCATGCGTCGCCCGGCATg                                                             <  2:284405/127‑1 (MQ=255)                                                                    cACTCGCGGCCCTGGGCAAATTCGTTACCCATGAACAGTAGTTTCTTGCCCGGGCATGCCCACATCCAGCCATAGTAGGCGCGCAGGTTCGCGAATTTCTGCCATGCGTCGCCCGGCATGCGGTCGAGAATCg                                                <  2:623200/133‑1 (MQ=255)                                                                            gCCCTGGGCAAATTCGTTACCCATGAACAGTAGTTTCTTGCCCGGGCATGCCCACATCCAGCCATAGTAGGCGCGCAGGTTCGCGAATTTCTGCCATGCGTCGCCCGGCATGCGGTCGAGAATCGATTTTTTa                                        <  1:5740/133‑1 (MQ=255)                                                                                  ggCAAATTCGTTACCCATGAACAGTAGTTTCTTGCCCGGGCATGCCCACATCCAGCCATAGTAGGCGCGCAGGTTCGCGAATTTCTGCCATGCGTCGCCCGGCATGCGGTCGAGAATCGATTTTTTACCGTgg                                  >  2:504456/1‑133 (MQ=255)                                                                                   gCAAATTCGTTACCCATGAACAGTAGTTTCTTGCCCGGGCATGCCCACATCCAGCCATAGTAGGCGCGCAGGTTCGCGAATTTCTGCCATGCGTCGCCCGGCATGCGGTCGAGAATCGATTTTTTACCGTGGa                                 <  1:560236/133‑1 (MQ=255)                                                                                                                  tGCCCGGGCATGCCCACATCCAGCCATAGTAGGCGCGCAGGTTCGCGAATTTCTACCATGCGTCGCCCGGCATGCGGTCGAGAATCGATTTTTTACCGTGGACCACTTCATCATGCGACAACGGCAGGACGaa  >  2:601073/1‑133 (MQ=255)                                                                                                                          |                                                                                                                             CTGGACACCGTGGTGCCAGTTATCGCCGCCTTCCAACAGATGCCAGTCGAGGCTGGCGTCATGGTTCCACTCGCGGCCCTGGGCAAATTCGTTACCCATGAACAGTAGTTTCTTGCCCGGGAATGCCCACATCCAGCCATAGTAGGCGCGCAGGTTCGCGAATTTCTGCCATGCGTCGCCCGGCATGCGGTCGAGAATCGATTTTTTACCGTGGACCACTTCATCATGCGACAACGGCAGGACGAA  >  NC_000913/3571703‑3571948

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