Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,734,543	(C)8→7	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,536	0	C	.	96.1%	203.7 / 2.1	51	coding (499/732 nt)	pgeF	polyphenol oxidase YfiH
Reads supporting (aligned to +/- strand):  ref base C (1/1);  new base . (19/30);  total (20/31)
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.86e-01

CCGAAAATTTGCTCAACACCCACGTTCGCCAGACGCTGCCGGGCAAGCTGATAAATATCCGCCAGATACTTATCACCATGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCTCCCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAACAGGAAACCGTCTCTTCCAGCACGCCTGCGCACAGTCCACGCCAGCCAGCATG  >  NC_000913/2734400‑2734665                                                                                                                                         |                                                                                                                                  ccgAAAATTTGCTCAACACCCACGTTCGCCAGACGCTGCCGGGCAAGCTGATAAATATCCGCCAGATACTTATCACCATGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCTcccc                                                                                                                                <  1:1112981/140‑1 (MQ=255)       aTTTGCTCAACACCCACGTTCGCCAGACGCTGCCGGGCAAGCTGATAAATATCCGCCAGATACTTATCACCATGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCTcccccccact                                                                                                                          <  2:1148360/140‑4 (MQ=255)            cTCAACACCCACGTTCGCCAGACGCTGCCGGGCAAGCTGATAAATATCCGCCAGATACTTATCACCATGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCT‑CCCCCCCACTTCGaa                                                                                                                    <  1:911513/140‑1 (MQ=255)            cTCAACACCCACGTTCGCCAGACGCTGCCGGGCAAGCTGATAAATATCCGCCAGATACTTATCACCATGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCT‑CCCCCCCACTTCGaa                                                                                                                    <  1:779207/140‑1 (MQ=255)                        gTTCGCCAGACGCTGCCGGGCAAGCTGATAAATATCCGCCAGATACTTATCACCATGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCT‑CCCCCCCACTTCGAACGCGCGTGGAcc                                                                                                        <  1:87115/140‑1 (MQ=255)                        gTTCGCCAGACGCTGCCGGGCAAGCTGATAAATATCCGCCAGATACTTATCACCATGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCT‑CCCCCCCACTTCGAACGCGCGTGGAcc                                                                                                        <  1:753048/140‑1 (MQ=255)                              cAGACGCTGCCGGGCAAGCTGATAAATATCCGCCAGATACTTATCACCATGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCT‑CCCCCCCACTTCGAACGCGCGTGGACCAATTGc                                                                                                  <  1:775452/140‑1 (MQ=255)                              cAGACGCTGCCGGGCAAGCTGATAAATATCCGCCAGATACTTATCACCATGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCT‑CCCCCCCACTTCGAACGCGCGTGGACCAATTGc                                                                                                  <  1:912299/140‑1 (MQ=255)                              cAGACGCTGCCGGGCAAGCTGATAAATATCCGCCAGATACTTATCACCATGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCT‑CCCCCCCACTTCGAACGCGCGTGGACCAATTGc                                                                                                  <  1:747489/140‑1 (MQ=255)                              cAGACGCTGCCGGGCAAGCTGATAAATATCCGCCAGATACTTATCACCATGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCT‑CCCCCCCACTTCGAACGCGCGTGGACCAATTGc                                                                                                  <  1:735433/140‑1 (MQ=255)                                   gCTGCCGGGCAAGCTGATAAATATCCGCCAGATACTTATCACCATGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCT‑CCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGcc                                                                                             <  1:1264608/140‑1 (MQ=255)                                   gCTGCCGGGCAAGCTGATAAATATCCGCCAGATACTTATCACCATGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCT‑CCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGcc                                                                                             <  1:1296725/140‑1 (MQ=255)                                            cAAGCTGATAAATATCCGCCAGATACTTATCACCATGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCT‑CCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAgg                                                                                    <  1:882951/140‑1 (MQ=255)                                                   aTAAATATCCGCCAGATACTTATCACCATGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCT‑CCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAat                                                                             <  1:249547/140‑1 (MQ=255)                                                             gCCAGATACTTATCACCATGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCT‑CCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGa                                                                   >  2:252751/1‑140 (MQ=255)                                                             gCCAGATACTTATCACCATGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCT‑CCCCCCCACTTCAAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGa                                                                   >  2:916612/1‑140 (MQ=255)                                                                    aCTTATCACCATGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCT‑CCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAg                                                            <  1:347262/140‑1 (MQ=255)                                                                    aCTTATCACCATGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCT‑CCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAg                                                            <  1:359964/140‑1 (MQ=255)                                                                    aCTTATCACCATGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCT‑CCCCCCCACTTCAAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAg                                                            <  1:916612/140‑1 (MQ=255)                                                                         tCACCATGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCT‑CCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCaaaa                                                       >  1:1328022/1‑140 (MQ=255)                                                                           aCCATGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCT‑CCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAACa                                                     >  1:613336/1‑140 (MQ=255)                                                                             cATGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCT‑CCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAACAgg                                                   >  2:415538/1‑140 (MQ=255)                                                                             cATGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCT‑CCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAACAgg                                                   >  1:824934/1‑140 (MQ=255)                                                                               tGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCTCCCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAACAGGa                                                  >  1:210889/1‑140 (MQ=255)                                                                               tGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCT‑CCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAACAGGaa                                                 >  2:506830/1‑140 (MQ=255)                                                                               tGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCT‑CCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAACAGGaa                                                 >  1:287299/1‑140 (MQ=255)                                                                               tGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCT‑CCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAACAGGaa                                                 >  1:1277326/1‑140 (MQ=255)                                                                                gCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCT‑CCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAACAGGaaa                                                >  1:68024/1‑140 (MQ=255)                                                                                gCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCT‑CCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAACAGGaaa                                                <  1:415538/140‑1 (MQ=255)                                                                                 cTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCT‑CCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAACAGGAAAc                                               >  1:320147/1‑140 (MQ=255)                                                                                    aatgaaAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCT‑CCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAACAGGAAACCGt                                            >  2:1294030/1‑140 (MQ=255)                                                                                    aatgaaAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCT‑CCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAACAGGAAACCGt                                            >  1:431405/1‑140 (MQ=255)                                                                                    aatgaaAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCT‑CCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAACAAAACAGGAAACCGt                                            >  2:387267/1‑140 (MQ=255)                                                                                     atgaaAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCT‑CCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAACAGGAAACCGtc                                           >  1:325720/1‑140 (MQ=255)                                                                                        aaaGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCT‑CCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAACAGGAAACCGTctct                                        <  2:325720/140‑1 (MQ=255)                                                                                               cACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCT‑CCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAACAGGAAACCGTCTCTTCCAGCa                                 <  2:320147/140‑1 (MQ=255)                                                                                               cACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCT‑CCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAACAGGAAACCGTCTCTTCCAGCa                                 <  1:1153670/140‑1 (MQ=255)                                                                                                   tGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCT‑CCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAACAGGAAACCGTCTCTTCCAGCACGcc                             <  2:287299/140‑1 (MQ=255)                                                                                                       ttAGCGTCTACTGCCATAAACGCCTCGCGAACCT‑CCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAACAGGAAACCGTCTCTTCCAGCACGCCTgcg                         <  1:252751/140‑1 (MQ=255)                                                                                                         aGCGTCTACTGCCATAAACGCCTCGCGAACCT‑CCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAACAGGAAACCGTCTCTTCCAGCACGCCTGCGca                       <  2:1328022/140‑1 (MQ=255)                                                                                                         aGCGTCTACTGCCATAAACGCCTCGCGAACCT‑CCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAACAGGAAACCGTCTCTTCCAGCACGCCTGCGca                       <  1:99108/140‑1 (MQ=255)                                                                                                         aGCGTCTACTGCCATAAACGCCTCGCGAACCT‑CCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAACAGGAAACCGTCTCTTCCAGCACGCCTGCGca                       <  1:1294030/140‑1 (MQ=255)                                                                                                         aGCGTCTACTGCCATAAACGCCTCGCGAACCT‑CCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAACAGGAAACCGTCTCTTCCAGCACGCCTGCGca                       <  1:1041151/140‑1 (MQ=255)                                                                                                         aGCGTCTACTGCCATAAACGCCTCGCGAACCT‑CCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAACAGGAAACCGTCTCTTCCAGCACGCCTGCGca                       <  2:824934/140‑1 (MQ=255)                                                                                                         aGCGTCTACTGCCATAAACGCCTCGCGAACCT‑CCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAACAGGAAACCGTCTCTTCCAGCACGCCTGCGca                       <  1:438571/140‑1 (MQ=255)                                                                                                         aGCGTCTACTGCCATAAACGCCTCGCGAACCT‑CCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAACAGGAAACCGTCTCTTCCAGCACGCCTGCGca                       >  2:967196/1‑140 (MQ=255)                                                                                                          gCGTCTACTGCCATAAACGCCTCGCGAACCTCCCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAACAGGAAACCGTCTCTTCCAGCACGCCTGCGca                       <  2:210889/140‑1 (MQ=255)                                                                                                                  tGCCATAAACGCCTCGCGAACCT‑CCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAACAGGAAACCGTCTCTTCCAGCACGCCTGCGCACAGTCCACg              >  1:1165895/1‑140 (MQ=255)                                                                                                                     cATAAACGCCTCGCGAACCT‑CCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAACAGGAAACCGTCTCTTCCAGCACGCCTGCGCACAGTCCACgcca           <  1:967196/140‑1 (MQ=255)                                                                                                                           cGCCTCGCGAACCT‑CCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAACAGGAAACCGTCTCTTCCAGCACGCCTGCGCACAGTCCACGCcagccagc     <  2:68024/140‑1 (MQ=255)                                                                                                                              cTCGCGAACCT‑CCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAACAGGAAACCGTCTCTTCCAGCACGCCTGCGCACAGTCCACgc             >  1:838684/1‑129 (MQ=255)                                                                                                                              cTCGCGAACCT‑CCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAACAGGAAACCGTCTCTTCCAGCACGCCTGCGCACAGTCCACgc             <  2:838684/129‑1 (MQ=255)                                                                                                                              cTCGCGAACCT‑CCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAACAGGAAACCGTCTCTTCCAGCACGCCTGCGCACAGTCCACGCCAGCCAGCATg  >  1:471692/1‑140 (MQ=255)                                                                                                                                         |                                                                                                                                  CCGAAAATTTGCTCAACACCCACGTTCGCCAGACGCTGCCGGGCAAGCTGATAAATATCCGCCAGATACTTATCACCATGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCTCCCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAACAGGAAACCGTCTCTTCCAGCACGCCTGCGCACAGTCCACGCCAGCCAGCATG  >  NC_000913/2734400‑2734665

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: ‑