Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,567,517	G→A	H173Y (CAC→TAC)	glgA ←	glycogen synthase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,567,517	0	G	A	100.0%	138.1 / NA	38	H173Y (CAC→TAC)	glgA	glycogen synthase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (17/21);  total (17/21)

GGAAAGAGATTTGTCCGTTGAATTCCAGCCCATGAATATTAAAGAATGACCATGGCAATTGGATGTCATTCATGTGATGTGCATAAAACATGCCTTGATAGGCCAGGTTGTGCACAGTAAACACCGACTTCGCCGGACGCCCGCGCGCCGC  >  NC_000913/3567442‑3567592                                                                            |                                                                            ggAAAGAGATTTGTCCGTTGAATTCCAGCCCATGAATATTAAAGAATGACCATGGCAATTGGATGTCATTCatgtaatgtg                                                                        >  1:678550/1‑81 (MQ=255)     agagATTTGTCCGTTGAATTCCAGCCCATGAATATTAAAGAATGACCATGGCAATTGGATGTCATTCATGTAATGTGCAt                                                                     >  1:678564/1‑80 (MQ=255)      gagaTTTGTCCGTTGAATTCCAGCCCATGAATATTAAAGAATGACCATGGCAATTGGATGTCATTCATGTAATGTGCATa                                                                    >  1:678559/1‑80 (MQ=255)          tttGTCCGTTGAATTCCAGCCCATGAATATTAAAGAATGACCATGGCAATTGGATGTCATTCATGTAATGTGCATAAAACa                                                               <  1:678516/81‑1 (MQ=255)               ccGTTGAATTCCAGCCCATGAATATTAAAGAATGACCATGGCAATTGGATGTCATTCATGTAATGTGCATAAAACATGCCt                                                          >  1:678575/1‑81 (MQ=255)                cGTTGAATTCCAGCCCATGAATATTAAAGAATGACCATGGCAATTGGATGTCATTCATGTAATGTGCATAAAACATGCCt                                                          <  1:678517/80‑1 (MQ=255)                 gTTGAATTCCAGCCCATGAATATTAAAGAATGACCATGGCAATTGGATGTCATTCATGTAATGTGCATAAAACATGCCTTg                                                        <  2:678518/81‑1 (MQ=255)                  ttGAATTCCAGCCCATGAATATTAAAGAATGACCATGGCAATTGGATGTCATTCATGTAATGTGCATAAAACATGCCTTGa                                                       <  2:678520/81‑1 (MQ=255)                  ttGAATTCCAGCCCATGAATATTAAAGAATGACCATGGCAATTGGATGTCATTCATGTAATGTGCATAAAACATGCCTTGa                                                       <  1:678521/81‑1 (MQ=255)                  ttGAATTCCAGCCCATGAATATTAAAGAATGACCATGGCAATTGGATGTCATTCATGTAATGTGCATAAAACATGCCTTGa                                                       <  1:678519/81‑1 (MQ=255)                          cAGCCCATGAATATTAAAGAATGACCATGGCAATTGGATGTCATTCATGTAATGTGCATAAAACATGCCTTGATAGGCCa                                                <  1:678522/80‑1 (MQ=255)                          cAGCCCATGAATATTAAAGAATGACCATGGCAATTGGATGTCATTCATGTAATGTGCATAAAACATGCCTTGATAGGCCAg                                               <  1:678523/81‑1 (MQ=255)                           aGCCCATGAATATTAAAGAATGACCATGGCAATTGGATGTCATTCATGTAATGTGCATAAAACATGCCTTGATAGGCCAgg                                              >  2:678573/1‑81 (MQ=255)                            gCCCATGAATATTAAAGAATGACCATGGCAATTGGATGTCATTCATGTAATGTGCATAAAACATGCCTTGATAGGCCAGGt                                             <  2:678524/81‑1 (MQ=255)                               cATGAATATTAAAGAATGACCATGGCAATTGGATGTCATTCATGTAATGTGCATAAAACATGCCTTGATAGGCCAGGTtgt                                          <  1:678525/81‑1 (MQ=255)                                  gAATATTAAAGAATGACCATGGCAATTGGATGTCATTCATGTAATGTGCATAAAACATGCCTTGATAGGCCAGGTTGTGca                                       >  2:678567/1‑81 (MQ=255)                                   aaTATTAAAGAATGACCATGGCAATTGGATGTCATTCATGTAATGTGCATAAAACATGCCTTGATAGGCCAGGTTGTGca                                       <  2:678526/80‑1 (MQ=255)                                   aaTATTAAAGAATGACCATGGCAATTGGATGTCATTCATGTAATGTGCATAAAACATGCCTTGATAGGCCAGGTTGTGca                                       >  1:678571/1‑80 (MQ=255)                                   aaTATTAAAGAATGACCATGGCAATTGGATGTCATTCATGTAATGTGCATAAAACATGCCTTGATAGGCCAGGTTGTGca                                       <  1:678527/80‑1 (MQ=255)                                      atTAAAGAATGACCATGGCAATTGGATGTCATTCATGTAATGTGCATAAAACATGCCTTGATAGGCCAGGTTGTGCACAg                                    >  1:678566/1‑80 (MQ=255)                                         aaaGAATGACCATGGCAATTGGATGTCATTCATGTAATGTGCATAAAACATGCCTTGATAGGCCAGGTTGTGCACAGTaaa                                >  2:678565/1‑81 (MQ=255)                                             aaTGACCATGGCAATTGGATGTCATTCATGTAATGTGCATAAAACATGCCTTGATAGGCCAGGTTGTGCACAGTAAACAcc                            >  1:678558/1‑81 (MQ=255)                                             aaTGACCATGGCAATTGGATGTCATTCATGTAATGTGCATAAAACATGCCTTGATAGGCCAGGTTGTGCACAGTAAACAcc                            <  1:678528/81‑1 (MQ=255)                                               tGACCATGGCAATTGGATGTCATTCATGTAATGTGCATAAAACATGCCTTGATAGGCCAGGTTGTGCACAGTAAACACCGa                          <  1:678529/81‑1 (MQ=255)                                                gACCATGGCAATTGGATGTCATTCATGTAATGTGCATAAAACATGCCTTGATAGGCCAGGTTGTGCACAGTAAACACCGAc                         >  2:678562/1‑81 (MQ=255)                                                gACCATGGCAATTGGATGTCATTCATGTAATGTGCATAAAACATGCCTTGATAGGCCAGGTTGTGCACAGTAAACACCGAc                         >  1:678563/1‑81 (MQ=255)                                                  ccATGGCAATTGGATGTCATTCATGTAATGTGCATAAAACATGCCTTGATAGGCCAGGTTGTGCACAGTAAACACCGACtt                       >  1:678561/1‑81 (MQ=255)                                                   cATGGCAATTGGATGTCATTCATGTAATGTGCATAAAACATGCCTTGATAGGCCAGGTTGTGCACAGTAAACACCGACTTc                      <  2:678530/81‑1 (MQ=255)                                                     tGGCAATTGGATGTCATTCATGTAATGTGCATAAAACATGCCTTGATAGGCCAGGTTGTGCACAGTAAACACTGACTTCGc                    <  1:678531/81‑1 (MQ=255)                                                          aTTGGATGTCATTCATGTAATGTGCATAAAACATGCCTTGATAGGCCAGGTTGTGCATAGTAAACACCGATTTCGCCGTAc               >  1:678580/1‑81 (MQ=255)                                                          aTTGGATGTCATTCATGTAATGTGCATAAAACATGCCTTGATAGGCCAGGTTGTGCACAGTAAACACCGACTTCGCCGGAc               >  1:678579/1‑81 (MQ=255)                                                           ttGGATGTCATTCATGTAATGTGCATAAAACATGCCTTGATAGGCCAGGTTGTGCACAGTAAACACCGACTTCGCCGGACg              <  1:678532/81‑1 (MQ=255)                                                              gATGTCATTCATGTAATGTGCATAAAACATGCCTTGATAGGCCAGGTTGTGCACAGTAAACACCGACTTCGCCGGACGccc           <  1:678533/81‑1 (MQ=255)                                                               aTGTCATTCATGTAATGTGCATAAAACATGCCTTGATAGGCCAGGTTGTGCACAGTAAACACCGACTTCGCCGGACGCCcg          <  2:678534/81‑1 (MQ=255)                                                                 gTCATTCATGTAATGTGCATAAAACATGCCTTGATAGGCCAGGTTGTGCACAGTAAACACCGACTTCGCCGGACGCCcgcg        >  1:678586/1‑81 (MQ=255)                                                                 gTCATTCATGTAATGTGCATAAAACATGCCTTGATAGGCCAGGTTGTGCACAGTAAACACCGACTTCGCCGGACGCCcgcg        >  2:678577/1‑81 (MQ=255)                                                                  tcattcatGTAATGTGCATAAAACATGCCTTGATAGGCCAGGTTGTGCACAGTAAACACCGACTTCGCCGGACGCCcgcgc       <  2:678535/81‑1 (MQ=255)                                                                    attcatGTAATGTGCATAAAACATGCCTTGATAGGCCAGGTTGTGCACAGTAAACACCGACTTCGCCGGACGCCcgcgcg      <  2:678536/80‑1 (MQ=255)                                                                      tcatGTAATGTGCATAAAACATGCCTTGATAGGCCAGGTTGTGCACAGTAAACACCGACTTCGCCGGACGCCCGCGcgccg   >  2:678592/1‑81 (MQ=255)                                                                        atgtaatgtgCATAAAACATGCCTTGATAGGCCAGGTTGTGCACAGTAAACACCGACTTCGCCGGACGCCCGCGcgccgc  >  2:678568/1‑80 (MQ=255)                                                                            |                                                                            GGAAAGAGATTTGTCCGTTGAATTCCAGCCCATGAATATTAAAGAATGACCATGGCAATTGGATGTCATTCATGTGATGTGCATAAAACATGCCTTGATAGGCCAGGTTGTGCACAGTAAACACCGACTTCGCCGGACGCCCGCGCGCCGC  >  NC_000913/3567442‑3567592

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