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	88.5%	76.5 / 3.6	26	coding (500/732 nt)	pgeF	polyphenol oxidase YfiH
Reads supporting (aligned to +/- strand):  ref base . (1/2);  new base C (13/10);  total (14/12)
Fisher's exact test for biased strand distribution p-value = 5.80e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CAACACCCACGTTCGCCAGACGCTGCCGGGCAAGCTGATAAATATCCGCCAGATACTTATCACCATGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCT‑CC‑CCCCC‑CACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAACAGGAAACCGTCTCTTCCAGCACGCCTGCGCACAGTCCACGCCAGCCAGCA  >  NC_000913/2734413‑2734663                                                                                                                            |                                                                                                                                   cAACACCCACGTTCGCCAGACGCTGCCGGGCAAGCTGATAAATATCCGCCAGATACTTATCACCATGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCTCCC‑CCCCC‑CACTTCGaa                                                                                                                  <  1:158986/140‑1 (MQ=255)              cGCCAGACGCTGCCGGGCAAGTTGATAAATATCCGCGAGATACTTATCACCATGCTGGATGAAAGCTGCACTTGCTTTAGCGTGTACTGCCCTAAACGCCGCGGCAACAT‑CC‑CCCCCTCACTCCGAACGCGCGTGGACCa                                                                                                     >  2:341703/1‑140 (MQ=255)               gCCAGACGCTGCCGGGCAAGCTGATAAATATCCGCCAGATACTTATCACCATGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCTCCC‑CCCCC‑CACTTCGAACGCGCGTGGACCaa                                                                                                    >  2:684453/1‑140 (MQ=255)                        tGCCGGGCAAGCTGATAAATATCCGCCAGATACTTATCACCATGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCTCCC‑CCCCC‑CACTTCGAACGCGCGTGGACCAATTGCCGGcc                                                                                           <  1:767866/140‑1 (MQ=255)                               cAAGCTGATAAATATCCGCCAGATATTTAACACCATGATGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCCCGCGAACCTCCC‑CCCCC‑CACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCa                                                                                    <  1:341703/140‑1 (MQ=255)                                 aGCTGATAAATATCCGCCAGATACTTATCACCATGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCTCCC‑CCCCC‑CACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAgg                                                                                  <  1:684453/140‑1 (MQ=255)                                            atCCGCCAGATACTTATCACCATGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCTCCC‑CCCCC‑CACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATAttt                                                                       >  1:1606/1‑140 (MQ=255)                                            atCCGCCAGATACTTATCACCAAGCTGAATGCAAGCTGCTCTTGCTTTTGCGTCTACTGCCATAAACGCCTCGCGAACCTCCC‑CCCCC‑CACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATAttt                                                                       <  2:288846/140‑1 (MQ=255)                                                gCCAGATACTTATCACCATGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCTCCC‑CCCCC‑CACTTCGAACGCGCGTGGGCCAATTGCCGGCCCCAACCAGGCGAGAATATTTTCCg                                                                   >  1:288846/1‑140 (MQ=255)                                                gCCAGATACTTATCACCATGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCTCCC‑CCCCC‑CACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCg                                                                   >  1:474540/1‑140 (MQ=255)                                                gCCAGATACTTATCACCATGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCT‑CC‑CCCCC‑CACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCgg                                                                  >  1:851792/1‑140 (MQ=255)                                                  cAGATACTTATCACCATGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCTCCC‑CCCCC‑CACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGa                                                                 <  2:1606/140‑1 (MQ=255)                                                       aCTTATCACCATGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCTCCC‑CCCCC‑CACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATc                                                            <  2:474540/140‑1 (MQ=255)                                                         ttATCACCATGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCTCCC‑CCCCC‑CACTTCGAACGCGCGTTGACCAATTGCCGACCCGAACCAGGCAAGAAAATTTTCCGGATTATCAg                                                          >  2:281070/1‑140 (MQ=255)                                                                  tGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCTCCC‑CCCCC‑CACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAACAgg                                                 >  1:35537/1‑140 (MQ=255)                                                                   gCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCTCCC‑CCCCC‑CACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATTAGCAAAACAGGa                                                >  2:184851/1‑140 (MQ=255)                                                                   gCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCTCCC‑CCCCC‑CACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAACAGGa                                                >  1:515361/1‑140 (MQ=255)                                                                             aGCTGCACTTGCTTTAGCGTCTTCTCCCACAAACGCCCCGCGAACCTCCC‑CCCCC‑CACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAACAGGAAACCGTctct                                      <  1:281070/140‑1 (MQ=255)                                                                                 gCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCTCCC‑CCCCC‑CCCTCCGAACGCCCGTGGACCAATTGCCGGCCCTAAACAGGCGAGAATATTTTCCGGATTATCAGCAAAACAGGAAACCGTCTCCTCCa                                  >  2:772694/1‑140 (MQ=255)                                                                                 gCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCTCCC‑CCCCC‑CACTTCCAACCCCCGTGGACCAATTGCCGGCGCTAACCAGGCGAGAATATTTTCCGCATTATTAGCAAAACATGAAACCCTCTCTTCCa                                  >  2:765047/1‑140 (MQ=255)                                                                                    cTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCTCCC‑CCCCC‑CACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAACAGGAAACCGTCTCTTCCAGCa                               <  1:184851/140‑1 (MQ=255)                                                                                         ggacgggTCTACTGCCATAAACACCCCGCGAGCCG‑CCTCCCCC‑CACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAACAGGAAACCGTCTCTTCCAGCTCGCCt                          <  1:323951/134‑1 (MQ=255)                                                                                            aGCGTCTACTGCCATAAACGCCTCGCGAACCTCCC‑CCCCC‑CACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAAAATTTTCCGGATTATCAGCAAAACAGGAAACCGTCTCTTCCAGCCCGCCTgcc                       >  1:426924/1‑139 (MQ=255)                                                                                             gCGTCTACTGCCATAAACGCCTCGCGAACCT‑CC‑CCCCC‑CACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAACAGGAAACCGTCTCTTCCAGCACGCCTGCGca                     <  2:851792/140‑1 (MQ=255)                                                                                              cGTCTACTGCCATAAGCGCCTCGCGAACCT‑CC‑CCCCCACACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAACGGGAAACCGTCTCTTCCAGCACGCCTGCGca                     <  1:331309/140‑1 (MQ=255)                                                                                              cGTCTACTGCCATAAACGCCTCGCGAACCTCCC‑CCCCC‑CACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAACAGGAAACCGTCTCTTCCAGCACGCCTGCGca                     <  1:45540/140‑1 (MQ=255)                                                                                                          tAAACGCCTCGCGAACCTCCC‑CCCCC‑CACTTCGAACGCGCGGGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCCGAAAAACAGGAAACAGTCCCTTCCCGCGCGCCACGCCACAGTccctgtct         >  2:323951/1‑134 (MQ=255)                                                                                                                 cTCGCGAACCTCCC‑CCCCC‑CACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAACAGGAAACCGTCTCTTCCAGCACGCCTGCGCACAGTCCACGCCAGCCAGCa  >  1:53777/1‑140 (MQ=255)                                                                                                                 cTCGCGAACCTCCC‑CCCCC‑CACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAACAGGAAACCGTCTCTTCCAGCACGCCTGCGCACAGTCCACGCCAGCCAGCa  >  2:73904/1‑140 (MQ=255)                                                                                                                 cTCGCGAACCTCCC‑CCCCC‑CACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAACAGGAAACCGTCTCTTCCAGCACGCCTGCGCACAGTCCACGCCAACCAGCa  >  2:441516/1‑140 (MQ=255)                                                                                                                            |                                                                                                                                   CAACACCCACGTTCGCCAGACGCTGCCGGGCAAGCTGATAAATATCCGCCAGATACTTATCACCATGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCT‑CC‑CCCCC‑CACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAACAGGAAACCGTCTCTTCCAGCACGCCTGCGCACAGTCCACGCCAGCCAGCA  >  NC_000913/2734413‑2734663

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