Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,571,253	A→C	F23V (TTC→GTC)	glgX ←	limit dextrin alpha‑1,6‑glucohydrolase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,571,253	0	A	C	100.0%	179.7 / NA	49	F23V (TTC→GTC)	glgX	limit dextrin alpha‑1,6‑glucohydrolase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base C (28/21);  total (28/21)

TCATAGCGATGTTCCTGGCCATTGGCGTCAAAGACACACAGTTCTACCCGCTCGGCATGAGCGGAGAAAAGTGTGAAGTTGACGCCCTGACCGTCGTAATGCGCGCCGAGGGGAGCGGGTTTGCCAATGGCGAGTTGTGTCATTCTGCCTCC  >  NC_000913/3571177‑3571328                                                                             |                                                                            tCATAGCGATGTTCCTGGCCATTGGCGTCAAAGACACACAGTTCTACCCGCTCGGCATGAGCGGAGAAAAGTGTGACGtt                                                                          >  1:633150/1‑80 (MQ=255) tCATAGCGATGTTCCTGGCCATTGGCGTCAAAGACACACAGTTCTACCCGCTCGGCATGAGCGGAGAAAAGTGTGACGTTg                                                                         <  1:633090/81‑1 (MQ=255)   aTAGCGATGTTCCTGGCCATTGGCGTCAAAGACACACAGTTCTACCCGCTCGGCATGAGCGGAGAAAAGTGTGACGTTGAc                                                                       <  1:633091/81‑1 (MQ=255)   aTAGCGATGTTCCTGGCCATTGGCGTCAAAGACACACAGTTCTACCCGCTCGGCATGAGCGGAGAAAAGTGTGACGTTGAc                                                                       <  2:633092/81‑1 (MQ=255)   aTAGCGATGTTCCTGGCCATTGGCGTCAAAGACACACAGTTCTACCCGCTCGGCATGAGCGGAGAAAAGTGTGACGTTGAc                                                                       >  2:633159/1‑81 (MQ=255)   aTAGCGATGTTCCTGGCCATTGGCGTCAAAGACACACAGTTCTACCCGCTCGGCATGAGCGGAGAAAAGTGTGACGTTGAc                                                                       >  1:633174/1‑81 (MQ=255)   aTAGCGATGTTCCTGGCCATTGGCGTCAAAGACACACAGTTCTACCCGCTCGGCATGAGCGGAGAAAAGTGTGACGTTGAc                                                                       >  2:633164/1‑81 (MQ=255)       cGATGTTCCTGGCCATTGGCGTCAAAGACACACAGTTCTACCCGCTCGGCATGAGCGGAGAAAAGTGTGACGTTGACGccc                                                                   <  1:633093/81‑1 (MQ=255)         aTGTTCCTGGCCATTGGCGTCAAAGACACACAGTTCTACCCGCTCGGCATGAGCGGAGAAAAGTGTGACGTTGACGCCCt                                                                  >  2:633145/1‑80 (MQ=255)            ttCCTGGCCATTGGCGTCAAAGACACACAGTTCTACCCGCTCGGCATGAGCGGAGAAAAGTGTGACGTTGACGCCCTGAcc                                                              <  2:633094/81‑1 (MQ=255)              ccTGGCCATTGGCGTCAAAGACACACAGTTCTACCCGCTCGGCATGAGCGGAGAAAAGTGTGACGTTGACGCCCTGACcgt                                                            >  2:633172/1‑81 (MQ=255)                 ggCCATTGGCGTCAAAGACACACAGTTCTACCCGCTCGGCATGAGCGGAGAAAAGTGTGACGTTGACGCCCTGACcgtcg                                                          >  1:633149/1‑80 (MQ=255)                  gCCATTGGCGTCAAAGACACACAGTTCTACCCGCTCGGCATGAGCGGAGAAAAGTGTGACGTTGACGCCCTGACcgtcgt                                                         <  2:633095/80‑1 (MQ=255)                    cATTGGCGTCAAAGACACACAGTTCTACCCGCTCGGCATGAGCGGAGAAAAGTGTGACGTTGACGCCCTGACCGTCGTAAt                                                      >  1:633165/1‑81 (MQ=255)                      ttGGCGTCAAAGACACACAGTTCTACCCGCTCGGCATGAGCGGAGAAAAGTGTGACGTTGACGCCCTGACCGTCGTAATgc                                                    <  2:633097/81‑1 (MQ=255)                      ttGGCGTCAAAGACACACAGTTCTACCCGCTCGGCATGAGCGGAGAAAAGTGTGACGTTGACGCCCTGACCGTCGTAATgc                                                    <  1:633096/81‑1 (MQ=255)                       tGGCGTCAAAGACACACAGTTCTACCCGCTCGGCATGAGCGGAGAAAAGTGTGACGTTGACGCCCTGACCGTCGTAATgcg                                                   >  2:633169/1‑81 (MQ=255)                       tGGCGTCAAAGACACACAGTTCTACCCGCTCGGCATGAGCGGAGAAAAGTGTGACGTTGACGCCCTGACCGTCGTAATgcg                                                   >  1:633155/1‑81 (MQ=255)                        ggCGTCAAAGACACACAGTTCTACCCGCTCGGCATGAGCGGAGAAAAGTGTGACGTTGACGCCCTGACCGTCGTAATgcgc                                                  >  1:633139/1‑81 (MQ=255)                         gCGTCAAAGACACACAGTTCTACCCGCTCGGCATGAGCGGAGAAAAGTGTGACGTTGACGCCCTGACCGTCGTAATgcgc                                                  <  2:633098/80‑1 (MQ=255)                           gTCAAAGACACACAGTTCTACCCGCTCGGCATGAGCGGAGAAAAGTGTGACGTTGACGCCCTGACCGTCGTAATGCGCGcc                                               >  1:633163/1‑81 (MQ=255)                           gTCAAAGACACACAGTTCTACCCGCTCGGCATGAGCGGAGAAAAGTGTGACGTTGACGCCCTGACCGTCGTAATGCGCGcc                                               >  1:633157/1‑81 (MQ=255)                            tCAAAGACACACAGTTCTACCCGCTCGGCATGAGCGGAGAAAAGTGTGACGTTGACGCCCTGACCGTCGTAATGCGCGCCg                                              <  2:633099/81‑1 (MQ=255)                              aaaGACACACAGTTCTACCCGCTCGGCATGAGCGGAGAAAAGTGTGACGTTGACGCCCTGACCGTCGTAATGCGCGCCGAg                                            >  1:633158/1‑81 (MQ=255)                               aaGACACACAGTTCTACCCGCTCGGCATGCGCGGAGAAAAGTGTGACGTTGACGCCCTGACCGTCGTAATGCGCGCCGAgg                                           <  2:633100/81‑1 (MQ=255)                                      acaGTTCTACCCGCTCGGCATGAGCGGAGAAAAGTGTGACGTTGACGCCCTGACCGTCGTAATGCGCGCCGAGGGGAGCgg                                    >  2:633162/1‑81 (MQ=255)                                      acaGTTCTACCCGCTCGGCATGAGCGGAGAAAAGTGTGACGTTGACGCCCTGACCGTCGTAATGCGCGCCGAGGGGAGCgg                                    >  1:633152/1‑81 (MQ=255)                                        aGTTCTACCCGCTCGGCATGAGCGGAGAAAAGTGTGACGTTGACGCCCTGACCGTCGTAATGCGCGCCGAGGGGAGCGGGt                                  <  2:633101/81‑1 (MQ=255)                                         gTTCTACCCGCTCGGCATGAGCGGAGAAAAGTGTGACGTTGACGCCCTGACCGTCGTAATGCGCGCCGAGGGGAGCGGGtt                                 <  1:633102/81‑1 (MQ=255)                                          ttCTACCCGCTCGGCATGAGCGGAGAAAAGTGTGACGTTGACGCCCTGACCGTCGTAATGCGCGCCGAGGGGAGCGGGttt                                >  1:633171/1‑81 (MQ=255)                                          ttCTACCCGCTCGGCATGAGCGGAGAAAAGTGTGACGTTGACGCCCTGACCGTCGTAATGCGCGCCGAGGGGAGCGGGttt                                >  1:633173/1‑81 (MQ=255)                                               cccGCTCGGCATGAGCGGAGAAAAGTGTGACGTTGACGCCCTGACCGTCGTAATGCGCGCCGAGGGGAGCGGGTTTGcc                             <  2:633103/79‑1 (MQ=255)                                                      ggCATGAGCGGAGAAAAGTGTGACGTTGACGCCCTGACCGTCGTAATGCGCGCCGAGGGGAGCGGGTTTGCCAATGGCg                      <  1:633106/79‑1 (MQ=255)                                                      ggCATGAGCGGAGAAAAGTGTGACGTTGACGCCCTGACCGTCGTAATGCGCGCCGAGGGGAGCGGGTTTGCCAATGGCGa                     <  1:633104/80‑1 (MQ=255)                                                      ggCATGAGCGGAGAAAAGTGTGACGTTGACGCCCTGACCGTCGTAATGCGCGCCGAGGGGAGCGGGTTTGCCAATGGCGAg                    <  2:633105/81‑1 (MQ=255)                                                        cATGAGCGGAGAAAAGTGTGACGTTGACGCCCTGACCGTCGTAATGCGCGCCGAGGGGAGCGGGTTTGCCAATGGCGAGtt                  <  2:633107/81‑1 (MQ=255)                                                         aTGAGCGGAGAAAAGTGTGACGTTGACGCCCTGACCGTCGTAATGCGCGCCGAGGGGAGCGGGTTTGCCAATGGCGAGTtg                 >  1:633179/1‑81 (MQ=255)                                                         aTGAGCGGAGAAAAGTGTGACGTTGACGCCCTGACCGCCGTAATGCGCGCCGAGGGGAGCGGGTTTGCCAATGGCGAGTtg                 >  1:633188/1‑81 (MQ=255)                                                          tGAGCGGAGAAAAGTGTGACGTTGACGCCCTGACCGTCGTAATGCGCGCCGAGGGGAGCGGGTTTGCGAATGGCGAGTtgt                <  2:633108/81‑1 (MQ=255)                                                          tGAGCGGAGAAAAGTGTGACGTTGACGCCCTGACCGTCGTAATGCGCGCCGAGGGGAGCGGGTTTGCCAATGGCGAGTtgt                >  1:633148/1‑81 (MQ=255)                                                              cGGAGAAAAGTGTGACGTTGACGCCCTGACCGTCGTAATGCGCGCCGAGGGGAGCGGGTTTGCCAATGGCGAGTTGTGTCa            >  1:633177/1‑81 (MQ=255)                                                               ggAGAAAAGTGTGACGTTGACGCCCTGACCGTCGTAATGCGCGCCGAGGGGAGCGGGTTTGCCAATGGCGAGTTGTGTCAt           >  1:633161/1‑81 (MQ=255)                                                                gagaAAAGTGTGACGTTGACGCCCTGACCGTCGTAATGCGCGCCGAGGGGAGCGGGTTTGCCAATGGCGAGTTGTGTCAtt          <  2:633109/81‑1 (MQ=255)                                                                  gaAAAGTGTGACGTTGACGCCCTGACCGTCGTAATGCGCGCCGAGGGGAGCGGGTTTGCCAATGGCGAGTTGTGTCATTCt        >  2:633167/1‑81 (MQ=255)                                                                  gaAAAGTGTGACGTTGACGCCCTGACCGTCGTAATGCGCGCCGAGGGGAGCGGGTTTGCCAATGGCGAGTTGTGTCATTCt        >  2:633178/1‑81 (MQ=255)                                                                     aaGTGTGACGTTGACGCCCTGACCGTCGTAATGCGCGCCGAGGGGAGCGGGTTTGCCAATGGCGAGTTGTGTCATTCTGcc     >  1:633168/1‑81 (MQ=255)                                                                     aaGTGTGACGTTGACGCCCTGACCGTCGTAATGCGCGCCGAGGGGAGCGGGTTTGCCAATGGCGAGTTGTGTCATTCTGcc     <  1:633110/81‑1 (MQ=255)                                                                      aGTGTGACGTTGACGCCCTGACCGTCGTAATGCGCGCCGAGGGGAGCGGGTTTGCCAATGGCGAGTTGTGTCATTCTGCCt    >  1:633160/1‑81 (MQ=255)                                                                        tgtgACGTTGACGCCCTGACCGTCGTAATGCGCGCCGAGGGGAGCGGGTTTGCCAATGGCGAGTTGTGTCATTCTGCCTcc  >  1:633189/1‑81 (MQ=255)                                                                             |                                                                            TCATAGCGATGTTCCTGGCCATTGGCGTCAAAGACACACAGTTCTACCCGCTCGGCATGAGCGGAGAAAAGTGTGAAGTTGACGCCCTGACCGTCGTAATGCGCGCCGAGGGGAGCGGGTTTGCCAATGGCGAGTTGTGTCATTCTGCCTCC  >  NC_000913/3571177‑3571328

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