Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,571,289	Δ1 bp	coding (31/1974 nt)	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,286	0	G	.	100.0%	165.0 / NA	39	coding (34/1974 nt)	glgX	limit dextrin alpha‑1,6‑glucohydrolase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base . (18/21);  total (18/21)

GTTCTACCCGCTCGGCATGAGCGGAGAAAAGTGTGAAGTTGACGCCCTGACCGTCGTAATGCGCGCCGAGGGGAGCGGGTTTGCCAATGGCGAGTTGTGTCATTCTGCCTCCCGAACCAGCCAGATAGTGGCCAGCGGTGGTAGCG  >  NC_000913/3571217‑3571362                                                                      |                                                                             gTTCTACCCGCTCGGCATGAGCGGAGAAAAGTGTGAAGTTGACGCCCTGACCGTCGTAATGCGCGCCGA‑GGGAGCGGGttt                                                                  <  1:622143/81‑1 (MQ=255) gTTCTACCCGCTCGGCATGAGCGGAGAAAAGTGTGAAGTTGACGCCCTGACCGTCGTAATGCGCGCCGA‑GGGAGCGGGttt                                                                  <  2:622143/81‑1 (MQ=255)  ttCTACCCGCTCGGCATGAGCGGAGAAAAGTGTGAAGTTGACGCCCTGACCGTCGTAATGCGCGCCGA‑GGGAGCGGGTTTg                                                                 >  2:622214/1‑81 (MQ=255)    cTACCCGCTCGGCATGAGCGGAGAAAAGTGTGAAGTTGACGCCCTGACCGTCGTAATGCGCGCCGA‑GGGAGCGGGTTTGc                                                                <  1:622145/80‑1 (MQ=255)    cTACCCGCTCGGCATGAGCGGAGAAAAGTGTGAAGTTGACGCCCTGACCGTCGTAATGCGCGCCGA‑GGGAGCGGGTTTGc                                                                <  2:622145/80‑1 (MQ=255)      aCCCGCTCGGCATGAGCGGAGAAAAGTGTGAAGTTGACGCCCTGACCGTCGTAATGCGCGCCGA‑GGGAGCGGGTTTGCCaa                                                             >  1:622205/1‑81 (MQ=255)       cccGCTCGGCATGAGCGGAGAAAAGTGTGAAGTTGACGCCCTGACCGTCGTAATGCGCGCCGA‑GGGAGCGGGTTTGcc                                                               <  2:622146/78‑1 (MQ=255)       cccGCTCGGCATGAGCGGAGAAAAGTGTGAAGTTGACGCCCTGACCGTCGTAATGCGCGCCGA‑GGGAGCGGGTTTGCCAAt                                                            <  2:622147/81‑1 (MQ=255)          gCTCGGCATGAGCGGAGAAAAGTGTGAAGTTGACGCCCTGACCGTCGTAATGCGCGCCGA‑GGGAGCGGGTTTGCCAATGGc                                                         <  1:622149/81‑1 (MQ=255)             cGGCATGAGCGGAGAAAAGTGTGAAGTTGACGCCCTGACCGTCGTAATGCGCGCCGA‑GGGAGCGGGTTTGCCAATGGCGAg                                                      >  2:622222/1‑81 (MQ=255)             cGGCATGAGCGGAGAAAAGTGTGAAGTTGACGCCCTGACCGTCGTAATGCGCGCCGA‑GGGAGCGGGTTTGCCAATGGCGAg                                                      >  1:622213/1‑81 (MQ=255)               gCATGAGCGGAGAAAAGTGTGAAGTTGACGCCCTGACCGTCGTAATGCGCGCCGA‑GGGAGCGGGTTTGCCAATGGCGAGtt                                                    >  1:622204/1‑81 (MQ=255)               gCATGAGCGGAGAAAAGTGTGAAGTTGACGCCCTGACCGTCGTAATGCGCGCCGA‑GGGAGCGGGTTTGCCAATGGCGAGtt                                                    <  2:622149/81‑1 (MQ=255)                cATGAGCGGAGAAAAGTGTGAAGTTGACGCCCTGACCGTCGTAATGCGCGCCGA‑GGGAGCGGGTTTGCCAATGGCGAGTtg                                                   <  1:622151/81‑1 (MQ=255)                    aGCGGAGAAAAGTGTGAAGTTGACGCCCTGACCGTCGTAATGCGCGCCGA‑GGGAGCGGGTTTGCCAATGGCGAGTTgtgt                                                >  1:622208/1‑80 (MQ=255)                    aGCGGAGAAAAGTGTGAAGTTGACGCCCTGACCGTCGTAATGCGCGCCGA‑GGGAGCGGGTTTGCCAATGGCGAGTTGTGTc                                               <  1:622152/81‑1 (MQ=255)                      cGGAGAAAAGTGTGAAGTTGACGCCCTGACCGTCGTAATGCGCGCCGA‑GGGAGCGGGTTTGCCAATGGCGAGTTGTGTCAt                                             <  2:622152/81‑1 (MQ=255)                         agaAAAGTGTGAAGTTGACGCCCTGACCGTCGTAATGCGCGCCGA‑GGGAGCGGGTTTGCCAATGGCGAGTTGTGTCATTCt                                          >  2:622202/1‑81 (MQ=255)                         agaAAAGTGTGAAGTTGACGCCCTGACCGTCGTAATGCGCGCCGA‑GGGAGCGGGTTTGCCAATGGCGAGTTGTGTCATTCt                                          <  2:622153/81‑1 (MQ=255)                             aaGTGTGAAGTTGACGCCCTGACCGTCGTAATGCGCGCCGA‑GGGAGCGGGTTTGCCAATGGCGAGTTGTGTCATTCTGCCt                                      >  2:622223/1‑81 (MQ=255)                               gtgtGAAGTTGACGCCCTGACCGTCGTAATGCGCGCCGA‑GGGAGCGGGTTTGCCAATGGCGAGTTGTGTCATTCTGCCTcc                                    <  1:622155/81‑1 (MQ=255)                                tgtgAAGTTGACGCCCTGACCGTCGTAATGCGCGCCGA‑GGGAGCGGGTTTGCCAATGGCGAGTTGTGTCATTCTGCCTccc                                   <  2:622155/81‑1 (MQ=255)                                  tgAAGTTGACGCCCTGACCGTCGTAATGCGCGCCGA‑GGGAGCGGGTTTGCCAATGGCGAGTTGTGTCATTCTGCCTCCCg                                  <  2:622156/80‑1 (MQ=255)                                   gAAGTTGACGCCCTGACCGTCGTAATGCGCGCCGA‑GGGAGCGGGTTTGCCAATGGCGAGTTGTGTCATTCTGCCTCCCGaa                                >  2:622208/1‑81 (MQ=255)                                   gAAGTTGACGCCCTGACCGTCGTAATGCGCGCCGA‑GGGAGCGGGTTTGCCAATGGCGAGTTGTGTCATTCTGCCTCCCGaa                                >  2:622216/1‑81 (MQ=255)                                         gACGCCCTGACCGTCGTAATGCGCGCCGA‑GGGAGCGGGTTTGCCAATGGCGAGTTGTGTCATTCTGCCTCCCGAAccagc                           >  2:622210/1‑80 (MQ=255)                                         gACGCCCTGACCGTCGTAATGCGCGCCGA‑GGGAGCGGGTTTGCCAATGGCGAGTTGTGTCATTCTGCCTCCCGAAccagc                           <  2:622157/80‑1 (MQ=255)                                          aCGCCCTGACCGTCGTAATGCGCGCCGA‑GGGAGCGGGTTTGCCAATGGCGAGTTGTGTCATTCTGCCTCCCGAAccagcc                          <  2:622158/80‑1 (MQ=255)                                          aCGCCCTGACCGTCGTAATGCGCGCCGA‑GGGAGCGGGTTTGCCAATGGCGAGTTGTGTCATTCTGCCTCCCGAAccagc                           >  2:622213/1‑79 (MQ=255)                                            gCCCTGACCGTCGTAATGCGCGCCGA‑GGGAGCGGGTTCGCCAATGGCGAGTTGTGTCATTCTGCCTCCTGAAccagc                           >  2:622219/1‑77 (MQ=255)                                              ccTGACCGTCGTAATGCGCGCCGA‑GGGAGCGGGTTTGCCAATGGCGAGTTGTGTCATTCTGCCTCCCGAACCAGCCAGa                       >  2:622235/1‑79 (MQ=255)                                              ccTGACCGTCGTAATGCGCGCCGA‑GGGAGCGGGTTTGCCAATGGCGAGTTGTGTCATTCTGCCTCCCGAACCAGCCAGATa                     <  2:622159/81‑1 (MQ=255)                                               cTGACCGTCGTAATGCGCGCCGA‑GGGAGCGGGTTTGCCAATGGCGAGTTGTGTCATTATGCCCCCCGAACCAGCCAGATAg                    >  2:622228/1‑81 (MQ=255)                                                  aCCGTCGTAATGCGCGCCGA‑GGGAGCGGGTTTGCCAATGGCGAGTTGTGTCATTCTGCCTCCCGAACCAGCCAGATAGt                   <  1:622161/79‑1 (MQ=255)                                                     gtcgtAATGCGCGCCGA‑GGGAGCGGGTTTGCCAATGGCGAGTTGTGTCATTCTGCCTCCCGAACCAGCCAGATAGTGGc                >  2:622215/1‑79 (MQ=255)                                                      tcgtAATGCGCGCCGA‑GGGAGCGGGTTTGCCAATGGCGAGTTGTGTCATTCTGCCTCCCGAACCAGCCAGATAGTGGCCAg             >  2:622236/1‑81 (MQ=255)                                                        gtAATGCGCGCCGA‑GGGAGCGGGTTTGCCAATGGCGAGTTGTGTCATTCTGCCTCCCGAACCAGCCAGATAGTGGCCAGc            <  2:622161/80‑1 (MQ=255)                                                                cgcCGA‑GGGAGCGGGTTTGCCAATGGCGAGTTGTGTCATTCTGCCTCCCGAACCAGCCAGATAGTGGCCAGCGGTGGTAGc   <  2:622162/81‑1 (MQ=255)                                                                 gcCGA‑GGGAGCGGGTTTGCCAATGGCGAGTTGTGTCATTCTGCCTCCCGAACCAGCCAGATAGTGGCCAGCGGTGGTAGCg  >  2:622226/1‑81 (MQ=255)                                                                      |                                                                             GTTCTACCCGCTCGGCATGAGCGGAGAAAAGTGTGAAGTTGACGCCCTGACCGTCGTAATGCGCGCCGAGGGGAGCGGGTTTGCCAATGGCGAGTTGTGTCATTCTGCCTCCCGAACCAGCCAGATAGTGGCCAGCGGTGGTAGCG  >  NC_000913/3571217‑3571362

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: ‑