Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	423,467	A→G	D317G (GAC→GGC)	malZ →	maltodextrin glucosidase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	423,467	0	A	G	97.0%	106.0 / ‑4.0	33	D317G (GAC→GGC)	malZ	maltodextrin glucosidase
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base G (14/18);  total (15/18)
Fisher's exact test for biased strand distribution p-value = 4.55e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.03e-01

CCCCGAATCGCCCTGGCGCGACTGGTACTCGTTTAGTGATGATGGCACGGCGCTCGACTGGCTTGGCTATGCCAGCTTGCCGAAGCTGGATTATCAGTCGGAAAGTCTGGTGAATGAAATTTATCGCGGGGAAGACAGTATTGTCCGCCACTGGCTGAAAGCGCCGTGGAATATGGACGGCTGGCGGCTGGATGTGGTGCATATGCTGGGGGAGGCGGGTGGGGCGCGCAATAATATGCAGCACGTTGCCGGGATCACCGAAGCGGCGAAA  >  NC_000913/423333‑423603                                                                                                                                       |                                                                                                                                         ccccGAATCGCCCTGGCGCGACTGGTACTCGTTTAGTGATGATGGCACGGCGCTCGACTGGCTTGGCTATGCCAGCTTGCCGAAGCTGGATTATCAGTCGGAAAGTCTGGTGAATGAAATTTATCGCGGGGAAGGCAGTa                                                                                                                                     <  2:53968/140‑1 (MQ=255) ccccGAATCGCCCTGGCGCGACTGGTACTCGTTTAGTGATGATGGCACGGCGCTCGACTGGCTTGGCTATGCCAGCTTGCCGAAGCTGGATTATCAGTCGGAAAGTCTGGTGAATGAAATTTATCGCGGGGAAGGCAGTa                                                                                                                                     <  1:82919/140‑1 (MQ=255)    cGAATCGCCCTGGCGCGACTGGTACTCGTTTAGTGATGATGGCACGGCGCTCGACTGGCTTGGCTATGCCAGCTTGCCGAAGCTGGATTATCAGTCGGAAAGTCTGGTGAATGAAATTTATCGCGGGGAAGGCAGTATTg                                                                                                                                  >  1:485763/1‑140 (MQ=255)           cccTGGCGCGACTGGTACTCGTTTAGTGATGATGGCACGGCGCTCGACTGGCTTGGCTATGCCAGCTTGCCGAAGCTGGATTATCAGTCGGAAAGTCTGGTGAATGAAATTTATCGCGGGGAAGGCAGTATTGTCCGCCa                                                                                                                           <  1:346570/140‑1 (MQ=255)                    gACTGGTACTCGTTTAGTGATGATGGCACGGCGCTCGACTGGCTTGGCTATGCCAGCTTGCCGAAGCTGGATTATCAGTCGGAAAGTCTGGTGAATGAAATTTATCGCGGGGAAGGCAGTATTGTCCGCCACTGGCTGaa                                                                                                                  >  1:1692/1‑140 (MQ=255)                                    gTGATGATGGCACGGCGCTCGACTGGCTTGGCTATGCCAGCTTGCCGAAGCTGGATTATCAGTCGGAAAGTCTGGTGAATGAAATTTATCGCGGGGAAGGCAGTATTGTCCGCCACTGGCTGAAAGCGCCGTGGAATATg                                                                                                  <  1:88502/140‑1 (MQ=255)                                         gatgGCACGGCGCTCGACTGGCTTGGCTATGCCAGCTTGCCGAAGCTGGATTATCAGTCGGAAAGTCTGGTGAATGAAATTTATCGCGGGGAAGGCAGTATTGTCCGCCACTGGCTGAAAGCGCCGTGGAATATGGAcgg                                                                                             >  2:184709/1‑140 (MQ=255)                                          atgGCACGGCGCTCGACTGGCTTGGCTATGCCAGCTTGCCGAAGCTGGATTATCAGTCGGAAAGTCTGGTGAATGAAATTTATCGCGGGGAAGGCAGTATTGTCCGCCACTGGCTGAAAGCGCCGTGGAATATGGAcggc                                                                                            >  1:397875/1‑140 (MQ=255)                                             gCACGGCGCTCGACTGGCTTGGCTATGCCAGCTTGCCGAAGCTGGATTATCAGTCGGAAAGTCTGGTGAATGAAATTTATCGCGGGGAAGGCAGTATTGTCCGCCACTGGCTGAAAGCGCCGTGGAATATGGAcggctgg                                                                                         <  2:65626/140‑1 (MQ=255)                                              cACGGCGCTCGACTGGCTTGGCTATGCCAGCTTGCCGAAGCTGGATTATCAGTCGGAAAGTCTGGTGAATGAAATTTATCGCGGGGAAGGCAGTATTGTCCGCCACTGGCTGAAAGCGCCGTGGAATATGGAcggctggc                                                                                        >  1:97159/1‑140 (MQ=255)                                                 ggCGCTCGACTGGCTTGGCTATGCCAGCTTGCCGAAGCTGGATTATCAGTCGGAAAGTCTGGTGAATGAAATTTATCGCGGGGAAGGCAGTATTGTCCGCCACTGGCTGAAAGCGCCGTGGAATATGGAcggctggcggc                                                                                     <  1:95199/140‑1 (MQ=255)                                                   cgcTCGACTGGCTTGGCTATGCCAGCTTGCCGAAGCTGGATTATCAGTCGGAAAGTCTGGTGAATGAAATTTATCGCGGGGAAGGCAGTATTGTCCGCCACTGGCTGAAAGCGCCGTGGAATATGGAcggctggcggctg                                                                                   >  1:15117/1‑140 (MQ=255)                                                    gcTCGACTGGCTTGGCTATGCCAGCTTGCCGAAGCTGGATTATCAGTCGGAAAGTCTGGTGAATGAAATTTATCGCGGGGAAGGCAGTATTGTCCGCCACTGGCTGAAAGCGCCGTGGAATATGGAcggctggcggctgg                                                                                  <  2:544281/140‑1 (MQ=255)                                                         aCTGGCTTGGCTATGCCAGCTTGCCGAAGCTGGATTATCAGTCGGAAAGTCTGGTGAATGAAATTTATCGCGGGGAAGGCAGTATTGTCCGCCACTGGCTGAAAGCGCCGTGGAATATGGACGGCTGGCGGCTGGAtgtg                                                                             <  2:70531/140‑1 (MQ=255)                                                               ttggctATGCCAGCTTGCCGAAGCTGGATTATCAGTCGGAAAGTCTGGTGAATGAAATTTATCGCGGGGAAGGCAGTATTGTCCGCCACTGGCTGAAAGCGCCGTGGAATATGGACGGCTGGCGGCTGGATGTGGTGCat                                                                       <  1:243756/140‑1 (MQ=255)                                                                  gctATGCCAGCTTGCCGAAGCTGGATTATCAGTCGGAAAGTCTGGTGAATGAAATTTATCGCGGGGAAGGCAGTATTGTCCGCCACTGGCTGAAAg                                                                                                                >  2:533972/1‑96 (MQ=255)                                                                  gctATGCCAGCTTGCCGAAGCTGGATTATCAGTCGGAAAGTCTGGTGAATGAAATTTATCGCGGGGAAGGCAGTATTGTCCGCCACTGGCTGAAAg                                                                                                                <  1:533972/96‑1 (MQ=255)                                                                                ccGAAGCTGGATTATCAGTCGGAAAGTCTGGTGAATGAAATTTATCGCGGGGAAGGCAGTATTGTCCGCCACTGGCTGAAAGCGCCGTGGAATATGGACGGCTGGCGGCTGGATGTGGTGCATATGCTGGGGGAGGCggg                                                      >  2:375406/1‑140 (MQ=255)                                                                                ccGAAGCTGGATTATCAGTCGGAAAGTCTGGTGAATGAAATTTATCGCGGGGAAGGCAGTATTGTCCGCCACTGGCTGAAAGCGCCGTGGAATATGGACGGCTGGCGGCTGGATGTGGTGCATATGCTGGGGGAGGCggg                                                      >  1:113926/1‑140 (MQ=255)                                                                                      cTGGATTATCAGTCGGAAAGTCTGGTGAATGAAATTTATCGCGGGGAAGGCAGTATTGTCCGCCACTGGCTGAAAGCGCCGTGGAATATGGACGGCTGGCGGCTGGATGTGGTGCATATGCTGGGGGAGGCGGGGGGGgc                                                >  1:215263/1‑140 (MQ=255)                                                                                         gATTATCAGTCGGAAAGTCTGGTGAATGAAATTTATCGCGGGGAAGGCAGTATTGTCCGCCACTGGCTGAAAGCGCCGTGGAATATGGACGGCTGGCGGCTGGATGTGGTGCATATGCTGGGGGAGGCGGGTGGGgcgcg                                             <  1:326638/140‑1 (MQ=255)                                                                                         gATTATCAGTCGGAAAGTCTGGTGAATGAAATTTATCGCGGGGAAGGCAGTATTGTCCGCCACTGGCTGAAAGCGCCGTGGAATATGGACGGCTGGCGGCTGGATGTGGTGCATATGCTGGGGGAGGCGGGTGGGgcgcg                                             <  1:416908/140‑1 (MQ=255)                                                                                          aTTATCAGTCGGAAAGTCTGGTGAATGAAATTTATCGCGGGGAAGGCAGTATTGTCCGCCACTGGCTGAAAGCGCCGTGGAATATGGACGGCTGGCGGCTGGATGTGGTGCATATGCTGGGGGAGGCGGGTGGGgcgcgc                                            >  2:448574/1‑140 (MQ=255)                                                                                           ttATCAGTCGGAAAGTCTGGTGAATGAAATTTATCGCGGGGAAGGCAGTATTGTCCGCCACTGGCTGAAAGCGCCGTGGAATATGGACGGCTGGCGGCTGGATGTGGTGCATATGCTGGGGGAGGCGGGTGGGGCGCGCa                                           <  1:448574/140‑1 (MQ=255)                                                                                                   cGGAAAGTCTGGTGAATGAAATTTATCGCGGGGAAGGCAGTATTGTCCGCCACTGGCTGAAAGCGCCGTGGAATATGGACGGCTGGCGGCTGGATGTGGTGCATATGCTGGGGGAGGCGGGTgg                                                   >  2:181756/1‑124 (MQ=255)                                                                                                   cGGAAAGTCTGGTGAATGAAATTTATCGCGGGGAAGGCAGTATTGTCCGCCACTGGCTGAAAGCGCCGTGGAATATGGACGGCTGGCGGCTGGATGTGGTGCATATGCTGGGGGAGGCGGGTgg                                                   <  1:181756/124‑1 (MQ=255)                                                                                                    ggAAAGTCTGGTGAATGAAATTTATCGCGGGGAAGGCAGTATTGTCCGCCACTGGCTGAAAGCGCCGTGGAATATGGACGGCTGGCGGCTGGATGTGGTGCATATGCTGGGGGAGGCGGGTGGGGCGCGCAATAATATgc                                  <  1:375406/140‑1 (MQ=255)                                                                                                     gAAAGTCTGGTGAATGAAATTTATCGCGGGGAAGGCAGTATTGTCCGCCACTGGCTGAAAGCGCCGTGGAATATGGACGGCTGGCGGCTGGATGTGGTGCATATGCTGGGGGAGGCGGGTGGGGCGCGCAATAATATgca                                 >  2:402281/1‑140 (MQ=255)                                                                                                          tCTGGTGAATGAAATTTATCGCGGGGAAGGCAGTATTGTCCGCCACTGGCTGAAAGCGCCGTGGAATATGGACGGCTGGCGGCTGGATGTGGTGCATATGCTGGGGGAGGCGGGTGGGGCGCGCAATAATATGCAGCACg                            <  2:515392/140‑1 (MQ=255)                                                                                                                    gaaATTTATCGCGGGGAAGGCAGTATTGTCCGCCACTGGCTGAAAGCGCCGTGGAATATGGACGGCTGGCGGCTGGATGTGGTGCATATGCTGGGGGAGGCGGGTGGGGCGCGCAATAATATGCAGCACGTTGCCGGGAt                  >  2:424721/1‑140 (MQ=255)                                                                                                                               cgGGGAAGGCAGTATTGTCCGCCACTGGCTGAAAGCGCCGTGGAATATGGACGGCTGGCGGCTGGATGTGGTGCATATGCTGGGGGAGGCGGGTGGGGCGCGCAATAATATGCAGCACGTTGCCGGGATCACCGAAgcgg       <  2:1692/140‑1 (MQ=255)                                                                                                                                ggggAAGGCAGTATTGTCCGCCACTGGCTGAAAGCGCCGTGGAATATGGACGGCTGGCGGCTGGATGTGGTGCATATGCTGGGGGAGGCGGGTGGGGCGCGCAATAATATGCAGCACGTTGCCGGGATCACCGAAgcggc      <  1:424721/140‑1 (MQ=255)                                                                                                                                    aaGACAGTATTGTCCGCCACTGGCTGAAAGCGCCGTGGAATATGGACGGCTGGCGGCTGGATGTGGTGCATATGCTGGGGGAGGCGGGTGGGGCGCGCAATAATATGCAGCACGTTGCCGGGATCACCGAAGCGGCgaaa  >  1:84177/1‑140 (MQ=255)                                                                                                                                       |                                                                                                                                         CCCCGAATCGCCCTGGCGCGACTGGTACTCGTTTAGTGATGATGGCACGGCGCTCGACTGGCTTGGCTATGCCAGCTTGCCGAAGCTGGATTATCAGTCGGAAAGTCTGGTGAATGAAATTTATCGCGGGGAAGACAGTATTGTCCGCCACTGGCTGAAAGCGCCGTGGAATATGGACGGCTGGCGGCTGGATGTGGTGCATATGCTGGGGGAGGCGGGTGGGGCGCGCAATAATATGCAGCACGTTGCCGGGATCACCGAAGCGGCGAAA  >  NC_000913/423333‑423603

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