Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,734,543	(C)8→9	coding (492/732 nt)	pgeF ←	polyphenol oxidase YfiH

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,734,535	1	.	C	87.8%	126.7 / 10.6	41	coding (500/732 nt)	pgeF	polyphenol oxidase YfiH
Reads supporting (aligned to +/- strand):  ref base . (2/3);  new base C (17/19);  total (19/22)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.73e-01

GCCGCCGAAAATTTGCTCAACACCCACGTTCGCCAGACGCTGCCGGGCAAGCTGATAAATATCCGCCAGATACTTATCACCATGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCT‑CCCCC‑CC‑CACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAACAGGAAACCGTCTCTTCCAGCACGCCTGCGCACAGTCCACGCCAGCCAGCA  >  NC_000913/2734396‑2734663                                                                                                                                             |                                                                                                                                   gccgccGAAAATTTGCTCAACACCCACGTTCGCCAGACGCTGCCGGGCAAGCTGATAAATATCCGCCAGATACTTATCACCATGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCt                                                                                                                                     <  1:23070/140‑1 (MQ=255) gccgccGAAAATTTGCTCAACACCCACGTTCGCCAGACGCTGCCGGGCAAGCTGATAAATATCCGCCAGATACTTATCACCATGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCt                                                                                                                                     >  2:23070/1‑140 (MQ=255)      cgAAAATTTGCTCAACACCCACGTTCGCCAGACGCTGCCGGGCAAGCTGATAAATATCCGCCAGATACTTATCACCATGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCT‑ccccc                                                                                                                               >  1:923135/1‑140 (MQ=255)         aaaTTTGCTCAACACCCACGTTCGCCAGACGCTGCCGGGCAAGCTGATAAATATCCGCCAGATACTTATCACCATGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCT‑ccccc‑cc‑c                                                                                                                          <  1:128769/140‑1 (MQ=255)                   aaCACCCACGTTCGCCAGACGCTGCCGGGCAAGCTGATAAATATCCGCCAGATACTTATCACCATGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCTCCCCCC‑CC‑CACTTCGAAc                                                                                                                 <  1:30970/140‑1 (MQ=255)                            gTTCGCCAGACGCTGCCGGGCAAGCTGATAAATATCCGCCAGATACTTATCACCATGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCT‑CCCCC‑CC‑CACTTCGAACGCGCGTGGAc                                                                                                       <  2:923135/140‑1 (MQ=255)                             ttCGCCAGACGCTGCCGGGCAAGCTGATAAATATCCGCCAGATACTTATCACCATGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCTCCCCCC‑CC‑CACTTCGAACGCGCGTGGAc                                                                                                       <  1:438549/140‑1 (MQ=255)                                   aGACGCTGCCGGGCAAGCTGATAAATATCCGCCAGATACTTATCACCATGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCT‑CCCCC‑CC‑CACTTCGAACGCGCGTGGACCAATTGc                                                                                                <  1:838633/140‑1 (MQ=255)                                    gACGCTGCCGGGCAAGCTGATAAATATCCGCCAGATACTTATCACCATGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCTCCCCCC‑CC‑CACTTCGAACGCGCGTGGACCAATTGc                                                                                                <  1:759325/140‑1 (MQ=255)                                    gACGCTGCCGGGCAAGCTGATAAATATCCGCCAGATACTTATCACCATGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCTCCCCCC‑CC‑CACTTCGAACGCGCGTGGACCAATTGc                                                                                                <  1:506245/140‑1 (MQ=255)                                         tGCCGGGCAAGCTGATAAATATCCGCCAGATACTTATCACCATGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCTCCCCCC‑CC‑CACTTCGAACGCGCGTGGACCAATTGCCGGcc                                                                                           <  1:215723/140‑1 (MQ=255)                                                     tGATAAATATCCGCCAGATACTTATCACCATGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCTCCCCCC‑CC‑CACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCga                                                                               >  1:946707/1‑140 (MQ=255)                                                        tAAATATCCGCCAGATACTTATCACCATGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCTCCCCCC‑CC‑CACTTCGAACGAGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGaa                                                                            >  2:669967/1‑140 (MQ=255)                                                                            aTCACCATGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCTCCCCCC‑CC‑CACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCa                                                        <  1:669967/140‑1 (MQ=255)                                                                              cACCATGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCTCCCCCC‑CC‑CACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCaaa                                                      >  1:478111/1‑140 (MQ=255)                                                                                   tGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCTCCCCCC‑CC‑CACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAACAgg                                                 >  1:581980/1‑140 (MQ=255)                                                                                   tGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCTCCCCCC‑CC‑CACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAACAgg                                                 >  2:434590/1‑140 (MQ=255)                                                                                   tGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCTCCCCCC‑CC‑CACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAACAgg                                                 >  1:952278/1‑140 (MQ=255)                                                                                    gCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCTCCCCCC‑CC‑CACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAACAGGa                                                >  2:641418/1‑140 (MQ=255)                                                                                    gCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCTCCCCCC‑CC‑CACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAACAGGa                                                >  2:51880/1‑140 (MQ=255)                                                                                    gCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCTCCCCCC‑CC‑CACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAACAGGa                                                >  2:245727/1‑140 (MQ=255)                                                                                    gCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCT‑CCCCC‑CCACACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAACAGGa                                                >  2:414145/1‑140 (MQ=255)                                                                                     cTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCTCCCCCC‑CC‑CACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAACAGGaa                                               <  1:218135/140‑1 (MQ=255)                                                                                             aaGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCTCCCCCC‑CC‑CACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAACAGGAAACCGtctc                                       >  2:552803/1‑140 (MQ=255)                                                                                             aaGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCTCCCCCC‑CC‑CACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAACAGGAAACCGtctc                                       >  2:15851/1‑140 (MQ=255)                                                                                              aGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCTCCCCCC‑CC‑CACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAACAGGaaa                                              <  1:456200/132‑1 (MQ=255)                                                                                              aGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCTCCCCCC‑CC‑CACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAACAGGaaa                                              >  2:456200/1‑132 (MQ=255)                                                                                                  gCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCTCCCCCC‑CC‑CACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAACa                                                   <  2:64790/123‑1 (MQ=255)                                                                                                  gCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCTCCCCCC‑CC‑CACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAACa                                                   >  1:64790/1‑123 (MQ=255)                                                                                                  gCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCTCCCCCC‑CC‑CACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAACAGGAAACCGTCTCTTCCa                                  >  1:730097/1‑140 (MQ=255)                                                                                                  gCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCT‑CCCCCTCC‑CACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAACAGGAAACCGTCTCTTCCa                                  >  2:387576/1‑140 (MQ=255)                                                                                                     cTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCTCCCCCC‑CC‑CACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAACAGGAAACCGTCTCTTCCAGCa                               <  1:245727/140‑1 (MQ=255)                                                                                                       tGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCTCCCCCC‑CC‑CACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAACAGGAAACCGTCTCTTCCAGCACg                             <  2:946707/140‑1 (MQ=255)                                                                                                          tttAGCGTCTACTGCCATAAACGCCTCGCGAACCTCCCCCC‑CC‑CACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAACAGGAAACCGTCTCTTCCAGCACGCCt                          <  1:641418/140‑1 (MQ=255)                                                                                                              gCGTCTACTGCCATAAACGCCTCGCGAACCT‑CCCCC‑CC‑CACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAACAGGAAACCGTCTCTTCCAGCACGCCTGCGca                     <  1:414231/140‑1 (MQ=255)                                                                                                               cGTCTACTGCCATAAACGCCTCGCGAACCTCCCCCC‑CC‑CACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAACAGGAAACCGTCTCTTCCAGCACGCCTGCGca                     <  1:697196/140‑1 (MQ=255)                                                                                                               cGTCTACTGCCATAAACGCCTCGCGAACCTCCCCCC‑CC‑CACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAACAGGAAACCGTCTCTTCCAGCACGCCTGCGca                     <  1:298685/140‑1 (MQ=255)                                                                                                               cGTCTACTGCCATAAACGCCTCGCGAACCTCCCCCC‑CC‑CACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAACAGGAAACCGTCTCTTCCAGCACGCCTGCGca                     <  1:552803/140‑1 (MQ=255)                                                                                                               cGTCTACTGCCATAAACGCCTCGCGAACCTCCCCCC‑CC‑CACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAACAGGAAACCGTCTCTTCCAGCACGCCTGCGca                     <  1:232484/140‑1 (MQ=255)                                                                                                               cGTCTACTGCCATAAACGCCTCGCGAACCTCCCCCC‑CC‑CACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAACAGGAAACCGTCTCTTCCAGCACGCCTGCGca                     <  1:427954/140‑1 (MQ=255)                                                                                                                   tACTGCCATAAACGCCTCGCGAACCTCCCCCC‑CC‑CACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAACAGGAAACCGTCTCTTCCAGCACGCCTGCGCACAGt                 >  2:411601/1‑140 (MQ=255)                                                                                                                         cATAAACGCCTCGCGAACCTCCCCCC‑CC‑CACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAACAGGAAACCGTCTCTTCCAGCACGCCTGCGCACAGTCCTCgc           <  1:414145/140‑1 (MQ=255)                                                                                                                               cGCCTCGCGAACCTCCCCCC‑CC‑CACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAACAGGAAACCGTCTCTTCCAGCACGCCTGCGCACAGTCCa              >  2:922880/1‑131 (MQ=255)                                                                                                                               cGCCTCGCGAACCTCCCCCC‑CC‑CACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAACAGGAAACCGTCTCTTCCAGCACGCCTGCGCACAGTCCa              <  1:922880/131‑1 (MQ=255)                                                                                                                                  cTCGCGAACCTCCCCCC‑CC‑CACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAACAGGAAACCGTCTCTTCCAGCACGCCTGCGCACAGTCCACGCCAGCCAGCa  >  2:552756/1‑140 (MQ=255)                                                                                                                                             |                                                                                                                                   GCCGCCGAAAATTTGCTCAACACCCACGTTCGCCAGACGCTGCCGGGCAAGCTGATAAATATCCGCCAGATACTTATCACCATGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCT‑CCCCC‑CC‑CACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAACAGGAAACCGTCTCTTCCAGCACGCCTGCGCACAGTCCACGCCAGCCAGCA  >  NC_000913/2734396‑2734663

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