Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,549,669	C→T	P134P (CCG→CCA)	malQ ←	4‑alpha‑glucanotransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,549,669	0	C	T	100.0%	117.0 / NA	32	P134P (CCG→CCA)	malQ	4‑alpha‑glucanotransferase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (19/13);  total (19/13)

CCCCAGTTTTTTTCCGATCGCAGCGTATAAAGCTGAACGCAGGCACCCCACAGCTTTTGTTTATTCAGCAACGCCTGCGGTTCGTAACAGCGTTTCGGGGCGACAATCACCCGGCAATGCGCGCGCTGGTCGTCCTGGGTGAGTGTCAGCGTGTG  >  NC_000913/3549592‑3549746                                                                              |                                                                              ccccAGTTTTTTTCCGATCGCAGCGTATAAAGCTGAACGCAGGCACCCCACAGCTTTTGTTTATTCAGCAACGCCTGTGGt                                                                            >  2:681095/1‑81 (MQ=255) ccccAGTTTTTTTCCGATCGCAGCGTATAAAGCTGAACGCAGGCACCCCACAGCTTTTGTTTATTCAGCAACGCCTGTGGt                                                                            >  2:681092/1‑81 (MQ=255)  cccAGTTTTTTTCCGATCGCAGCGTATAAAGCTGAACGCAGGCACCCCACAGCTTTTGTTTATTCAGCAACGCCTGTGGtt                                                                           >  1:681107/1‑81 (MQ=255)       tttttttCCGATCGCAGCGTATAAAGCTGAACGCAGGCACCCCACAGCTTTTGTTTATTCAGCAACGCCTGTGGTTCGTaa                                                                      >  2:681116/1‑81 (MQ=255)         tttttCCGATCGCAGCGTATAAAGCTGAACGCAGGCACCCCACAGCTTTTGTTTATTCAGCAACGCCTGTGGTTCGTAAc                                                                     <  1:681055/80‑1 (MQ=255)          ttttCCGATCGCAGCGTATAAAGCTGAACGCAGGCACCCCACAGCTTTTGTTTATTCAGCAACGCCTGTGGTTCGTAACAg                                                                   >  1:681104/1‑81 (MQ=255)            ttCCGATCGCAGCGTATAAAGCTGAACGCAGGCACCCCACAGCTTTTGTTTATTCAGCAACGCCTGTGGTTCGTAACAGCg                                                                 >  1:681086/1‑81 (MQ=255)            ttCCGATCGCAGCGTATAAAGCTGAACGCAGGCACCCCACAGCTTTTGTTTATTCAGCAACGCCTGTGGTTCGTAACAGCg                                                                 >  1:681087/1‑81 (MQ=255)             tCCGATCGCAGCGTATAAAGCTGAACGCAGGCACCCCACAGCTTTTGTTTATTCAGCAACGCCTGTGGTTCGTAACAGCGt                                                                >  1:681094/1‑81 (MQ=255)               cGATCGCAGCGTATAAAGCTGAACGCAGGCACCCCACAGCTTTTGTTTATACAGCAACGCCTGTGGTTCGTAACAGCGtt                                                               <  2:681056/80‑1 (MQ=255)                 aTCGCAGCGTATAAAGCTGAACGCAGGCACCCCACAGCTTTTGTTTATTCAGCAACGCCTGTGGTTCGTAACAGCGTTTCg                                                            >  1:681098/1‑81 (MQ=255)                     cAGCGTATAAAGCTGAACGCAGGCACCCCACAGCTTTTGTTTATTCAGCAACGCCTGTGGTTCGTAACAGCGTTTCGGGGc                                                        >  1:681109/1‑81 (MQ=255)                      aGCGTATAAAGCTGAACGCAGGCACCCCACAGCTTTTGTTTATTCAGCAACGCCTGTGGTTCGTAACAGCGTTTCGGGGc                                                        >  2:681101/1‑80 (MQ=255)                      aGCGTATAAAGCTGAACGCAGGCACCCCACAGCTTTTGTTTATTCAGCAACGCCTGTGGTTCGTAACAGCGTTTCGGGGCg                                                       <  1:681057/81‑1 (MQ=255)                        cGTATAAAGCTGAACGCAGGCACCCCACAGCTTTTGTTTATTCAGCAACGCCTGTGGTTCGTAACAGCGTTTCGGGGCGAc                                                     <  1:681058/81‑1 (MQ=255)                         gTATAAAGCTGAACGCAGGCACCCCACAGCTTTTGTTTATTCAGCAACGCCTGTGGTTCGTAACAGCGTTTCGGGGCGAc                                                     >  1:681108/1‑80 (MQ=255)                           ataAAGCTGAACGCAGGCACCCCACAGCTTTTGTTTATTCAGCAACGCCTGTGGTTCGTAACAGCGTTTCGGGGCGACAAt                                                  >  1:681078/1‑81 (MQ=255)                                    aaCGCAGGCACCCCACAGCTTTTGTTTATTCAGCAACGCCTGTGGTTCGTAACAGCGTTTCGGGGCGACAATCACCCGGc                                          <  1:681059/80‑1 (MQ=255)                                    aaCGCAGGCACCCCACAGCTTTTGTTTATTCAGCAACGCCTGTGGTTCGTAACAGCGTTTCGGGGCGACAATCACCCGGc                                          >  2:681102/1‑80 (MQ=255)                                          ggCACCCCACAGCTTTTGTTTATTCAGCAACGCCTGTGGTTCGTAACAGCGTTTCGGGGCGACAATCACCCGGCAATgcgc                                   <  1:681060/81‑1 (MQ=255)                                            cACCCCACAGCTTTTGTTTATTCAGCAACGCCTGTGGTTCGTAACAGCGTTTCGGGGCGACAATCACCCGGCAATgcgcgc                                 >  1:681123/1‑81 (MQ=255)                                               cccACAGCTTTTGTTTATTCAGCAACGCCTGTGGTTCGTAACAGCGTTTCGGGGCGACAATCACCCGGCAATGCGCGCGCt                              <  2:681061/81‑1 (MQ=255)                                                    aGCTTTTGTTTATTCAGCAACGCCTGTGGTTCGTAACAGCGTTTCGGGGCGACAATCACCCGGCAATGCGCGCGCTGgtcg                         >  2:681110/1‑81 (MQ=255)                                                     gCTTTTGTTTATTCAGCAACGCCTGTGGTTCGTAACAGCGTTTCGGGGCGACAATCACCCGGCAATGCGCGCGCTGgtcgt                        <  1:681062/81‑1 (MQ=255)                                                         ttGTTTATTCAGCAACGCCTGTGGTTCGTAACAGCGTTTCGGGGCGACAATCACCCGGCAATGCGCGCGCTGGTCGTcc                      <  2:681063/79‑1 (MQ=255)                                                                ttCAGCAACGCCTGTGGTTCGTAACAGCGTTTCGGGGCGACAATCACCCGGCAATGCGCGCGCTGGTCGTCCTGGGTGAgt             >  2:681111/1‑81 (MQ=255)                                                                ttCAGCAACGCCTGTGGTTCGTAACAGCGTTTCGGGGCGACAATCACCCGGCAATGCGCGCGCTGGTCGTCCTGGGTGAg              <  1:681064/80‑1 (MQ=255)                                                                    gCAACGCCTGTGGTTCGTAACAGCGTTTCGGGGCGACAATCACCCGGCAATGCGCGCGCTGGTCGTCCTGGGTGAgtgt           >  1:681115/1‑79 (MQ=255)                                                                         gCCTGTGGTTCGTAACAGCGTTTCGGGGCGACAATCACCCGGCAATGCGCGCGCTGGTCGTCCTGGGTGAGTGTCAGCgtg    <  2:681065/81‑1 (MQ=255)                                                                          ccTGTGGTTCGTAACAGCGTTTCGGGGCGACAATCACCCGGCAATGCGCGCGCTGGTCGTCCTGGGTGAGTGTCAGCgtgt   >  1:681121/1‑81 (MQ=255)                                                                           cTGTGGTTCGTAACAGCGTTTCGGGGCGACAATCACCCGGCAATGCGCGCGCTGGTCGTCCTGGGTGAGTGTCAGCGtgtg  <  2:681066/81‑1 (MQ=255)                                                                           cTGTGGTTCGTAACAGCGTTTCGGGGCGACAATCACCCGGCAATGCGCGCGCTGGTCGTCCTGGGTGAGTGTCAGCGtgtg  <  2:681067/81‑1 (MQ=255)                                                                              |                                                                              CCCCAGTTTTTTTCCGATCGCAGCGTATAAAGCTGAACGCAGGCACCCCACAGCTTTTGTTTATTCAGCAACGCCTGCGGTTCGTAACAGCGTTTCGGGGCGACAATCACCCGGCAATGCGCGCGCTGGTCGTCCTGGGTGAGTGTCAGCGTGTG  >  NC_000913/3549592‑3549746

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