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	95.7%	161.3 / ‑0.1	47	coding (500/732 nt)	pgeF	polyphenol oxidase YfiH
Reads supporting (aligned to +/- strand):  ref base . (1/1);  new base C (26/19);  total (27/20)
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.70e-01

GCCGCCGAAAATTTGCTCAACACCCACGTTCGCCAGACGCTGCCGGGCAAGCTGATAAATATCCGCCAGATACTTATCACCATGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAA‑CCT‑‑CCCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAACAGGAAACCGTCTCTTCCAGCACGCCTGCGCACAGTCCACGCCAGCCAGCATGAGCG  >  NC_000913/2734396‑2734669                                                                                                                                              |                                                                                                                                        gccgccGAAAATTTGCTCAACACCCACGTTCGCCAGACGCTGCCGGGCAAGCTGATAAATATCCGCCAGATACTTATCACCATGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAA‑CCt                                                                                                                                          >  2:379558/1‑140 (MQ=255)              tGCTCAACACCCACGTTCGCCAGACGCTGCCGGGCAAGCTGATAAATATCCGCCAGATACTTATCACCATGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAA‑CCT‑‑cccccccccactt                                                                                                                           <  2:1066926/140‑6 (MQ=255)                         cACGTTCGCCAGACGCTGCCGGGCAAGCTGATAAATATCCGCCAGATACTTATCACCATGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAA‑CCTC‑CCCCCCCCACTTCGAACGCGCGt                                                                                                                 <  1:379558/140‑1 (MQ=255)                          aCGTTCGCCAGACGCTGCCGGGCAAGCTGATAAATATCCGCCAGATACTTATCACCATGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAA‑CCTC‑CCCCCCCCACTTCGaa                                                                                                                        >  2:152372/1‑132 (MQ=255)                          aCGTTCGCCAGACGCTGCCGGGCAAGCTGATAAATATCCGCCAGATACTTATCACCATGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAA‑CCTC‑CCCCCCCCACTTCGaa                                                                                                                        <  1:152372/132‑1 (MQ=255)                          aCGTTCGCCAGACGCTGCCGGGCAAGCTGATAAATATCCGCCAGATACTTATCACCATGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAA‑CCTC‑CCCCCCCCACTTCGAACGCGCGTg                                                                                                                >  2:51913/1‑140 (MQ=255)                               cGCCAGACGCTGCCGGGCAAGCTGATAAATATCCGCCAGATACTTATCACCATGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAA‑CCTC‑CCCCCCCCACTTCGAACGCGCGTGGACCa                                                                                                           >  2:1167418/1‑140 (MQ=255)                                    gACGCTGCCGGGCAAGCTGATAAATATCCGCCAGATACTTATCACCATGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAA‑CCTC‑CCCCCCCCACTTCGAACGCGCGTGGACCAATTGc                                                                                                      <  2:786083/140‑1 (MQ=255)                                    gACGCTGCCGGGCAAGCTGATAAATATCCGCCAGATACTTATCACCATGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAA‑CCTC‑CCCCCCCCACTTCGAACGCGCGTGGACCAATTGc                                                                                                      <  1:1167418/140‑1 (MQ=255)                                          gCCGGGCAAGCTGATAAATATCCGCCAGATACTTATCACCATGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAA‑CCTC‑CCCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGccc                                                                                                >  2:517930/1‑140 (MQ=255)                                           ccGGGCAAGCTGATAAATATCCGCCAGATACTTATCACCATGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAA‑CCTC‑CCCCCCCCACTTCGAACCCGCGGGGACCAATTGCCGGCCCt                                                                                               >  2:1157593/1‑140 (MQ=255)                                               gCAAGCTGATAAATATCCGCCAGATACTTATCACCATGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAA‑CCTC‑CCCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAAcc                                                                                           >  2:1182593/1‑140 (MQ=255)                                                  aGCTGATAAATATCCGCCAGATACTTATCACCATGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAA‑CCTC‑CCCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAgg                                                                                        <  1:1436102/140‑1 (MQ=255)                                                        tAAATATCCGCCAGATACTTATCACCATGCTGAATGAAAGCTGCGCTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAA‑CCTC‑CCCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGaa                                                                                  <  1:271417/140‑1 (MQ=255)                                                                 gCCAGATACTTATCACCATGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAA‑CCT‑‑CCCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCgg                                                                        >  1:355922/1‑140 (MQ=255)                                                                 gCCAGATACTTATCACCATGCTGAATGAAAGCTGCAATTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAA‑CCTC‑CCCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCg                                                                         >  2:1100144/1‑140 (MQ=255)                                                                                 cATGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAA‑CCTC‑CCCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAACa                                                         >  2:186595/1‑140 (MQ=255)                                                                                  aTGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAA‑CCTC‑CCCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAACAg                                                        >  2:1423462/1‑140 (MQ=255)                                                                                   tGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAA‑CCTC‑CCCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAACAgg                                                       >  1:586905/1‑140 (MQ=255)                                                                                   tGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAA‑CCTC‑CCCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAACAgg                                                       >  1:672253/1‑140 (MQ=255)                                                                                    gCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAA‑CCTC‑CCCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAACAGGa                                                      >  1:1266474/1‑140 (MQ=255)                                                                                    gCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAA‑CCTC‑CCCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAACAGGa                                                      >  2:1007129/1‑140 (MQ=255)                                                                                    gCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAA‑CCTC‑CCCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAACAGGa                                                      >  1:1324766/1‑140 (MQ=255)                                                                                    gCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAA‑CCTC‑CCCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAACAGGa                                                      >  1:1416318/1‑140 (MQ=255)                                                                                        aatgaaAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAA‑CCTC‑CCCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAACAGGAAAcc                                                  >  1:449614/1‑140 (MQ=255)                                                                                               gCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAA‑CCTCCCCCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAACAGGAAACCGTctct                                            <  1:947388/140‑1 (MQ=255)                                                                                                     cTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAA‑CCTC‑CCCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAACAGGAAACCGTCTCTTCCAGCa                                     <  1:517930/140‑1 (MQ=255)                                                                                                           ttAGCGTCTACTGCCATAAACGCCTCGCGAA‑CCTC‑CCCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAACAGGAAACCGTCTCTTCCAGCACGCCTg                               >  2:603928/1‑140 (MQ=255)                                                                                                             aGCGTCTACTGCCATAAACGCCTCGCGAA‑CCTC‑CCCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAACAGGAAACCGTCTCTTCCAGCACGCCTgcg                             >  2:853670/1‑140 (MQ=255)                                                                                                             aGCGTCTACTGCCATAAACGCCTCGCGAA‑CCTC‑CCCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAACAGGAAACCGTCTCTTCCAGCACGCCTgcg                             >  1:1281793/1‑140 (MQ=255)                                                                                                              gCGTCTACTGCCATAAACGCCTCGCGAA‑CCTC‑CCCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAACAGGAAACCGTCTCTTCCAGCACGCCTgcgc                            <  2:1380112/140‑1 (MQ=255)                                                                                                              gCGTCTACTGCCATAAACGCCTCGCGAA‑CCT‑‑CCCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAACAGGAAACCGTCTCTTCCAGCACGCCTGCGca                           <  2:355922/140‑1 (MQ=255)                                                                                                               cctttACTGCCATAAACGCCTCGCGAACCCC‑‑CCCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAACAGGAAACCGTCTCTTCCAGCACGCCTGCGca                           <  2:1281793/136‑1 (MQ=255)                                                                                                               cGTCTCCTGCCATAAACGCCTCGCGAA‑CCTC‑CCCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAACAGGAAACCGTCTCTTCCAGCACGCCTGCGca                           <  1:872981/140‑1 (MQ=255)                                                                                                               cGTCTACTGCCATAAACGCCTCGCGAA‑CCTC‑CCCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAACAGGAAACCGTCTCTTCCAGCACGCCTGCGca                           <  1:825326/140‑1 (MQ=255)                                                                                                               cGTCTACTGCCATAAACGCCTCGCGAA‑CCTC‑CCCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAACAGGAAACCGTCTCTTCCAGCACGCCTGCGca                           <  1:1060555/140‑1 (MQ=255)                                                                                                               cGTCTACTGCCATAAACGCCTCGCGAA‑CCTC‑CCCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAACAGGAAACCGTCTCTTCCAGCACGCCTGCGca                           <  2:1416318/140‑1 (MQ=255)                                                                                                               cGTCTACTGCCATAAACGCCTCGCGAA‑CCTC‑CCCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAACAGGAAACCGTCTCTTCCAGCACGCCTGCGca                           <  2:586905/140‑1 (MQ=255)                                                                                                               cGTCTACTGCCATAAACGCCTCGCGAA‑CCTC‑CCCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAACAGGAAACCGTCTCTTCCAGCACGCCTGCGca                           <  2:493202/140‑1 (MQ=255)                                                                                                                gTCTACTGCCATAAACGCCTCGCGAA‑CCTC‑CCCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAACAAGAAACCGTCTCTTCCAGCACGCCTGCGcac                          >  2:1379547/1‑140 (MQ=255)                                                                                                                   tACTGCCATAAACGCCTCGCGAA‑CCTC‑CCCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAACAGGAAACCGTCTCTTCCAGCACGcc                                 >  2:1244119/1‑130 (MQ=255)                                                                                                                   tACTGCCATAAACGCCTCGCGAA‑CCTC‑CCCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAACAGGAAACCGTCTCTTCCAGCACGcc                                 <  1:1244119/130‑1 (MQ=255)                                                                                                                   tACTGCCATAAACGCCTCGCGAA‑CCTC‑CCCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAACAGGAAACCGTCTCTTCCAGCACGCCTGCGCACAGt                       >  2:1055068/1‑140 (MQ=255)                                                                                                                    aCTGCCATAAACGCCTCGCGAA‑CCTC‑CCCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAACAGGAAACCGTCTCTTCCAGCACGCCTGCGCACAGt                       >  1:414880/1‑139 (MQ=255)                                                                                                                    aCTGCCATAAACGCCTCGCGAA‑CCTC‑CCCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAACAGGAAACCGTCTCTTCCAGCACGCCTGCGCACAGt                       <  2:414880/139‑1 (MQ=255)                                                                                                                         cATAAACGCCTCGCGAA‑CCTC‑CCCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAACAGGAAACCGTCTCTTCCAGCACGCCTGCGCACAGTCCACgc                 >  1:1362953/1‑140 (MQ=255)                                                                                                                           tAAACGCCTCGCGAA‑CCTC‑CCCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAACAGGAAACCGTCTCTTCCAGCACGCCTGCGCACAGTCCACgcca               >  2:1202450/1‑140 (MQ=255)                                                                                                                             aaCGCCTCGCGAA‑CCTC‑CCCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAACAGGAAACCGTCTCTTCCAGCACGCCTGCGCACAGTCCACgccagc             <  1:603928/140‑1 (MQ=255)                                                                                                                                  cTCGCGAA‑CCTC‑CCCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAACAGGAAACCGTCTCTTCCAGCACGCCTGCGCACAGTCCACGCCAGCCAGCa        >  1:386995/1‑140 (MQ=255)                                                                                                                                  cTCGCGAA‑CCTC‑CCCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAACAGGAAACCGTCTCTTCCAGCACGCCTGCGCACAGTCCACGCCAGCCAGCa        >  1:556508/1‑140 (MQ=255)                                                                                                                                        aa‑CCTC‑CCCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAACAGGAAACCGTCTCTTCCAGCACGCCTGCGCACAGTCCACGCCAGCCAGCATGAgcg  <  1:1007129/140‑1 (MQ=255)                                                                                                                                              |                                                                                                                                        GCCGCCGAAAATTTGCTCAACACCCACGTTCGCCAGACGCTGCCGGGCAAGCTGATAAATATCCGCCAGATACTTATCACCATGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAA‑CCT‑‑CCCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAACAGGAAACCGTCTCTTCCAGCACGCCTGCGCACAGTCCACGCCAGCCAGCATGAGCG  >  NC_000913/2734396‑2734669

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