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%	203.1 / NA	48	coding (34/1974 nt)	glgX	limit dextrin alpha‑1,6‑glucohydrolase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base . (29/19);  total (29/19)

AAAGACACACAGTTCTACCCGCTCGGCATGAGCGGAGAAAAGTGTGAAGTTGACGCCCTGACCGTCGTAATGCGCGCCGAGGGGAGCGGGTTTGCCAATGGCGAGTTGTGTCATTCTGCCTCCCGAACCAGCCAGATAGTGGCCAGCGGTGGTA  >  NC_000913/3571206‑3571359                                                                                 |                                                                          aaaGACACACAGTTCTACCCGCTCGGCATGAGCGGAGAAAAGTGTGAAGCTGACGCCCTGACCGTCGTAATGCGCGCCGAg                                                                           >  2:688450/1‑81 (MQ=255)   aGACACACAGTTCTACCCGCTCGGCATGAGCGGAGAAAAGTGTGAAGTTGACGCCCTGACCGTCGTAATGCGCGCCGAggg                                                                         <  2:688404/81‑1 (MQ=255)      cacacaGTTCTACCCGCTCGGCATGAGCGGAGAAAAGTGTGAAGTTGACGCCCTGACCGTCGTAATGCGCGCCGAgggagc                                                                      <  1:688406/81‑4 (MQ=255)            gTTCTACCCGCTCGGCATGAGCGGAGAAAAGTGTGAAGTTGACGCCCTGACCGTCGTAATGCGCGCCGA‑GGGAGCGGGttt                                                               <  2:688407/81‑1 (MQ=255)             ttCTACCCGCTCGGCATGAGCGGAGAAAAGTGTGAAGTTGACGCCCTGACCGTCGTAATGCGCGCCGA‑GGGAGCGGGTTTg                                                              <  2:688408/81‑1 (MQ=255)              tCTACCCGCTCGGCATGAGCGGAGAAAAGTGTGAAGTTGACGCCCTGACCGTCGTAATGCGCGCCGA‑GGGAGCGGGTTTGc                                                             >  2:688454/1‑81 (MQ=255)              tCTACCCGCTCGGCATGAGCGGAGAAAAGTGTGAAGTTGACGCCCTGACCGTCGTAATGCGCGCCGA‑GGGAGCGGGTTTGc                                                             >  2:688457/1‑81 (MQ=255)               cTACCCGCTCGGCATGAGCGGAGAAAAGTGTGAAGTTGACGCCCTGACCGTCGTAATGCGCGCCGA‑GGGAGCGGGTTTGc                                                             <  1:688409/80‑1 (MQ=255)                 aCCCGCTCGGCATGAGCGGAGAAAAGTGTGAAGTTGACGCCCTGACCGTCGTAATGCGCGCCGA‑GGGAGCGGGTTTGCCaa                                                          >  2:688458/1‑81 (MQ=255)                  cccGCTCGGCATGAGCGGAGAAAAGTGTGAAGTTGACGCCCTGACCGTCGTAATGCGCGCCGA‑GGGAGCGGGTTTGcc                                                            <  1:688410/78‑1 (MQ=255)                    cGCTCGGCATGAGCGGAGAAAAGTGTGAAGTTGACGCCCTGACCGTCGTAATGCGCGCCGA‑GGGAGCGGGTTTGCCAATgg                                                       <  1:688411/81‑1 (MQ=255)                     gCTCGGCATGAGCGGAGAAAAGTGTGAAGTTGACGCCCTGACCGTCGTAATGCGCGCCGA‑GGGAGCGGGTTTGCCAATGGc                                                      >  2:688459/1‑81 (MQ=255)                     gCTCGGCATGAGCGGAGAAAAGTGTGAAGTTGACGCCCTGACCGTCGTAATGCGCGCCGA‑GGGAGCGGGTTTGCCAATGGc                                                      >  1:688471/1‑81 (MQ=255)                         ggCATGAGCGGAGAAAAGTGTGAAGTTGACGCCCTGACCGTCGTAATGCGCGCCGA‑GGGAGCGGGTTTGCCAATGGCGAGt                                                  >  1:688476/1‑81 (MQ=255)                           cATGAGCGGAGAAAAGTGTGAAGTTGACGCCCTGACCGTCGTAATGCGCGCCGA‑GGGAGCGGGTTTGCCAATGGCGAGTtg                                                >  2:688486/1‑81 (MQ=255)                            aTGAGCGGAGAAAAGTGTGAAGTTGACGCCCTGACCGTCGTAATGCGCGCCGA‑GGGAGCGGGTTTGCCAATGGCGAGTtgt                                               <  1:688412/81‑1 (MQ=255)                               aGCGGAGAAAAGTGTGAAGTTGACGCCCTGACCGTCGTAATGCGCGCCGA‑GGGAGCGGGTTTGCCAATGGCGAGTTgtgt                                             <  1:688413/80‑1 (MQ=255)                                   gagaAAAGTGTGAAGTTGACGCCCTGACCGTCGTAATGCGCGCCGA‑GGGAGCGGGTTTGCCAATGGCGAGTTGTGTCATTc                                        >  2:688496/1‑81 (MQ=255)                                   gagaAAAGTGTGAAGTTGACGCCCTGACCGTCGTAATGCGCGCCGA‑GGGAGCGGGTTTGCCAATGGCGAGTTGTGTCATTc                                        <  1:688414/81‑1 (MQ=255)                                     gaAAAGTGTGAAGTTGACGCCCTGACCGTCGTAATGCGCGCCGA‑GGGAGCGGGTTTGCCAATGGCGAGTTGTGTCATTCTg                                      <  1:688415/81‑1 (MQ=255)                                       aaaGTGTGAAGTTGACGCCCTGACCGTCGTAATGCGCGCCGA‑GGGAGCGGGTTTGCCAATGGCGAGTTGTGTCATTCTGcc                                    >  2:688477/1‑81 (MQ=255)                                       aaaGTGTGAAGTTGACGCCCTGACCGTCGTAATGCGCGCCGA‑GGGAGCGGGTTTGCCAATGGCGAGTTGTGTCATTCTGcc                                    <  2:688416/81‑1 (MQ=255)                                       aaaGTGTGAAGTTGACGCCCTGACCGTCGTAATGCGCGCCGA‑GGGAGCGGGTTTGCCAATGGCGAGTTGTGTCATTCTGcc                                    >  1:688481/1‑81 (MQ=255)                                             tgAAGTTGACGCCCTGACCGTCGTAATGCGCGCCGA‑GGGAGCGGGTTTGCCAATGGCGAGTTGTGTCATTCTGCCTCCCg                               >  2:688515/1‑80 (MQ=255)                                              gAAGTTGACGCCCTGACCGTCGTAATGCGCGCCGA‑GGGAGCGGGTTTGCCAATGGCGAGTTGTGTCATTCTGCCTCCCGaa                             >  2:688465/1‑81 (MQ=255)                                              gAAGTTGACGCCCTGACCGTCGTAATGCGCGCCGA‑GGGAGCGGGTTTGCCAATGGCGAGTCGTGTCATCCTGCCTCCCGaa                             >  2:688489/1‑81 (MQ=255)                                               aaGTTGACGCCCTGACCGTCGTAATGCGCGCCGA‑GGGAGCGGGTTTGCCAATGGCGAGTTGTGTCATTCTGCCTCCCGAAc                            <  2:688417/81‑1 (MQ=255)                                                     aCGCCCTGACCGTCGTAATGCGCGCCGA‑GGGAGCGGGTTTGCCAATGGCGAGTTGTGTCATTCTGCCTCCCGAAccagcca                      <  2:688419/81‑1 (MQ=255)                                                     aCGCCCTGACCGTCGTAATGCGCGCCGA‑GGGAGCGGGTTTGCCAATGGCGAGTTGTGTCATTCTGCCTCCCGAAccagcca                      <  1:688418/81‑1 (MQ=255)                                                     aCGCCCTGACCGTCGTAATGCGCGCCGA‑GGGAGCGGGTTTGCCAATGGCGAGTTGTGTCATTCTGCCTCCCGAAccagc                        >  1:688482/1‑79 (MQ=255)                                                     aCGCCCTGACCGTCGTAATGCGCGCCGA‑GGGAGCGGGTTTGCCAATGGCGAGTTGTGTCATTCTGCCTCCCGAAccagc                        >  2:688470/1‑79 (MQ=255)                                                      cGCCCTGACCGTCGTAATGCGCGCCGA‑GGGAGCGGGTTTGCCAATGGCGAGTTGTGTCATTCTGCCTCCCGAAccagccag                     <  1:688420/81‑1 (MQ=255)                                                       gCCCTGACCGTCGTAATGCGCGCCGA‑GGGAGCGGGTTTGCCAATGGCGAGTTGTGTCATTCTGCCTCCCGAAccagc                        >  1:688485/1‑77 (MQ=255)                                                       gCCCTGACCGTCGTAATGCGCGCCGA‑GGGAGCGGGTTTGCCAATGGCGAGTTGTGTCATTCTGCCTCCCGAAccagc                        >  2:688503/1‑77 (MQ=255)                                                       gCCCTGACCGTCGTAATGCGCGCCGA‑GGGAGCGGGTTTGCCAATGGCGAGTTGTGTCATTCTGCCTCCCGAAccagc                        >  2:688494/1‑77 (MQ=255)                                                          cTGACCGTCGTAATGCGCGCCGA‑GGGAGCGGGTTTGCCAATGGCGAGTTGTGTCATTCTGCCTCCCGAACCAGCCAGATAg                 >  1:688508/1‑81 (MQ=255)                                                            gACCGTCGTAATGCGCGCCGA‑GGGAGCGGGTTTGCCAATGGCGAGTTGTGTCATTCTGCCTCCCGAACCAGCCAGATAGTg               >  2:688518/1‑81 (MQ=255)                                                            gACCGTCGTAATGCGCGCCGA‑GGGAGCGGGTTTGCCAATGGCGAGTTGTGTCATTCTGCCTCCCGAACCAGCCAGATAGTg               >  2:688492/1‑81 (MQ=255)                                                            gACCGTCGTAATGCGCGCCGA‑GGGAGCGGGTTTGCCAATGGCGAGTTGTGTCATTCTGCCTCCCGAACCAGCCAGATAGTg               >  2:688472/1‑81 (MQ=255)                                                             aCCGTCGTAATGCGCGCCGA‑GGGAGCGGGTTTGCCAATGGCGAGTTGTGTCATTCTGCCTCCCGAACCAGCCAGATAGTg               <  2:688421/80‑1 (MQ=255)                                                               cgtcgtAATGCGCGCCGA‑GGGAGCGGGTTTGCCAATGGCGAGTTGTGTCATTCTGCCTCCCGAACCAGCCAGATAGTGGc             >  2:688484/1‑80 (MQ=255)                                                               cgtcgtAATGCGCGCCGA‑GGGAGCGGGTTTGCCAATGGCGAGTTGTGTCATTCTGCCTCCCGAACCAGCCAGATAGTGGc             >  2:688510/1‑80 (MQ=255)                                                                gtcgtAATGCGCGCCGA‑GGGAGCGGGTTTGCCAATGGCGAGTTGTGTCATTCTGCCTCCCGAACCAGCCAGATAGTGGc             >  1:688474/1‑79 (MQ=255)                                                                gtcgtAATGCGCGCCGA‑GGGAGCGGGTTTGCCAATGGCGAGTTGTGTCATTCTGCCTCCCGAACCAGCCAGATAGTGGc             >  2:688495/1‑79 (MQ=255)                                                                gtcgtAATGCGCGCCGA‑GGGAGCGGGTTTGCCAATGGCGAGTTGTGTCATTCTGCCTCCCGAACCAGCCAGATAGTGGCCa           <  1:688422/81‑1 (MQ=255)                                                                  cgtAATGCGCGCCGA‑GGGAGCGGGTTTGCCAATGGCGAGTTGTGTCATTCTGCCTCCCGAACCAGCCAGATAGTGGCCAGc         >  1:688491/1‑81 (MQ=255)                                                                     aaTGCGCGCCGA‑GGGAGCGGGTTTGCCAATGGCGAGTTGTGTCATTCTGTCTCCCGAACCAGCCAGATAGTGGCCAGCggt      >  2:688502/1‑81 (MQ=255)                                                                      aTGCGCGCCGA‑GGGAGCGGGTTTGCCAATGGCGAGTTGTGTCATTCTGCCTCCCGAACCAGCCAGATAGTGGCCAGCggtg     <  1:688423/81‑1 (MQ=255)                                                                       tGCGCGCCGA‑GGGAGCGGGTTTGCCAATGGCGAGTTGTGTCATTCTGCCTCCCGAACCAGCCAGATAGTGGCCAGCggtgg    <  2:688424/81‑1 (MQ=255)                                                                        gcgcgcCGA‑GGGAGCGGGTTTGCCAATGGCGAGTTGTGTCATTCTGCCTCCCGAACCAGCCAGATAGTGGCCAGCggtgg    >  1:688487/1‑80 (MQ=255)                                                                         cgcgcCGA‑GGGAGCGGGTTTGCCAATGGCGAGTTGTGTCATTCTGCCTCCCGAACCAGCCAGATAGTGGCCAGCGGTGGTa  <  1:688425/81‑1 (MQ=255)                                                                                 |                                                                          AAAGACACACAGTTCTACCCGCTCGGCATGAGCGGAGAAAAGTGTGAAGTTGACGCCCTGACCGTCGTAATGCGCGCCGAGGGGAGCGGGTTTGCCAATGGCGAGTTGTGTCATTCTGCCTCCCGAACCAGCCAGATAGTGGCCAGCGGTGGTA  >  NC_000913/3571206‑3571359

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