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.3%	146.9 / ‑0.1	43	coding (500/732 nt)	pgeF	polyphenol oxidase YfiH
Reads supporting (aligned to +/- strand):  ref base . (0/2);  new base C (20/21);  total (20/23)
Fisher's exact test for biased strand distribution p-value = 4.91e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.91e-01

GCCGCCGAAAATTTGCTCAACACCCACGTTCGCCAGACGCTGCCGGGCAAGCTGATAAATATCCGCCAGATACTTATCACCATGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACC‑T‑‑CCCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAACAGGAAACCGTCTCTTCCAGCACGCCTGCGCACAGTCCACGCCAGCCAGCA  >  NC_000913/2734396‑2734663                                                                                                                                              |                                                                                                                                  gccgccGAAAATTTGCTCAACACCCACGTTCGCCAGACGCTGCCGGGCAAGCTGATAAATATCCGCCAGATACTTATCACCATGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACC‑t                                                                                                                                    <  1:417217/140‑1 (MQ=255)   cgccgAAAATTTGCTCAACACCCACGTTCGCCAGACGCTGCCGGGCAAGCTGATAAATATCCGCCAGATACTTATCACCATGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACC‑T‑‑cc                                                                                                                                <  1:512581/140‑1 (MQ=255)   cgccgAAAATTTGCTCAACACCCACGTTCGCCAGACGCTGCCGGGCAAGCTGATAAATATCCGCCAGATACTTATCACCATGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACC‑T‑‑cc                                                                                                                                <  1:733004/140‑1 (MQ=255)   cgccgAAAATTTGCTCAACACCCACGTTCGCCAGACGCTGCCGGGCAAGCTGATAAATATCCGCCAGATACTTATCACCATGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACC‑T‑‑cc                                                                                                                                <  1:781277/140‑1 (MQ=255)    gccgAAAATTTGCTCAACACCCACGTTCGCCAGACGCTGCCGGGCAAGCTGATAAATATCCGCCAGATACTTATCACCATGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACC‑T‑‑ccc                                                                                                                               <  1:800181/140‑1 (MQ=255)      cgAAAATTTGCTCAACACCCACGTTCGCCAGACGCTGCCGGGCAAGCTGATAAATATCCGCCAGATACTTATCACCATGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACC‑T‑‑ccccc                                                                                                                             >  2:704388/1‑140 (MQ=255)             ttGCTCAACACCCACGTTCGCCAGACGCTGCCGGGCAAGCTGATAAATATCCGCCAGATACTTATCACCATGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACC‑T‑‑cccccccccact                                                                                                                      <  1:498536/140‑5 (MQ=255)                  cAACACCCACGTTCGCCAGACGCTGCCGGGCAAGCTGATAAATATCCGCCAGATACTTATCACCATGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACC‑TC‑CCCCCCCCACTTCGaa                                                                                                                  <  1:209345/140‑1 (MQ=255)                                   aGACGCTGCCGGGCAAGCTGATAAATATCCGCCAGATACTTATCACCATGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACC‑TC‑CCCCCCCCACTTCGAACGCGCGTGGACCAATTg                                                                                                 <  1:229814/140‑1 (MQ=255)                                    gACGCTGCCGGGCAAGCTGATAAATATCCGCCAGATACTTATCACCATGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACC‑TC‑CCCCCCCCACTTCGAACGCGCGTGGACCAATTGc                                                                                                <  1:39875/140‑1 (MQ=255)                                    gACGCTGCCGGGCAAGCTGATAAATATCCGCCAGATACTTATCACCATGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACC‑TC‑CCCCCCCCACTTCGAACGCGCGTGGACCAATTGc                                                                                                <  1:239101/140‑1 (MQ=255)                                    gACGCTGCCGGGCAAGCTGATAAATATCCGCCAGATACTTATCACCATGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACC‑TC‑CCCCCCCCACTTCGAACGCGCGTGGACCAATTGc                                                                                                <  1:394394/140‑1 (MQ=255)                                              ggCAAGCTGATAAATATCCGCCAGATACTTATCACCATGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACC‑TC‑CCCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAAc                                                                                      <  1:838190/140‑1 (MQ=255)                                                  aGCTGATAAATATCCGCCAGATACTTATCACCATGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACC‑TC‑CCCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAgg                                                                                  <  1:1035078/140‑1 (MQ=255)                                                              tCCGCCAGATACTTATCACCATGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACC‑TC‑CCCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATAtttt                                                                      >  2:203725/1‑140 (MQ=255)                                                              tCCGCCAGATACTTATCACCATGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACC‑TC‑CCCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATAtttt                                                                      >  2:400745/1‑140 (MQ=255)                                                                 gCCAGATACTTATCACCATGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACC‑TC‑CCCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCg                                                                   >  2:1092952/1‑140 (MQ=255)                                                                      aTACTTATCACCATGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACC‑TC‑CCCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTa                                                              >  1:230834/1‑140 (MQ=255)                                                                                   tGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACC‑TC‑CCCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCTCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAACAgg                                                 >  2:792314/1‑140 (MQ=255)                                                                                   tGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACC‑TC‑CCCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAACAgg                                                 >  1:579338/1‑140 (MQ=255)                                                                                    gCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACC‑TC‑CCCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAACAGGa                                                >  1:507807/1‑140 (MQ=255)                                                                                    gCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACC‑TC‑CCCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAACAGGa                                                >  1:28011/1‑140 (MQ=255)                                                                                        aatgaaAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACC‑TC‑CCCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAACAGGaaa                                              <  2:485426/138‑1 (MQ=255)                                                                                        aatgaaAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACC‑TC‑CCCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAACAGGaaa                                              >  1:485426/1‑138 (MQ=255)                                                                                        aatgaaAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACC‑TC‑CCCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAACAGGAAAcc                                            >  2:761282/1‑140 (MQ=255)                                                                                        aatgaaAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACC‑TC‑CCCCCCCCACTTCGAACGCGCGTGGAACAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAACAGGAAAcc                                            >  1:515165/1‑140 (MQ=255)                                                                                           gaaAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACC‑TC‑CCCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAACAGGAAACCGtc                                         >  1:1110496/1‑140 (MQ=255)                                                                                             aaGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACC‑TC‑CCCCCCCCACTTCGAACGCGCGGGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAACAGGAAACCGcctc                                       >  1:479963/1‑140 (MQ=255)                                                                                              aGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACC‑TC‑CCCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAACAGGAAACCGTctct                                      >  2:736493/1‑140 (MQ=255)                                                                                                     cTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACC‑TC‑CCCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAACAGGAAACCGTCTCTTCCAGCa                               <  1:792314/140‑1 (MQ=255)                                                                                                     cTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACC‑TC‑CCCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAACAGGAAACCGTCTCTTCCAGCa                               <  1:700952/140‑1 (MQ=255)                                                                                                     cTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACC‑TC‑CCCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAACAGGAAACCGTCTCTTCCAGCa                               <  1:1092952/140‑1 (MQ=255)                                                                                                         cTTTAGCGTCTACTGCCATAAACGCCTCGCGAACC‑TCCCCCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAACAGGAAACCGTCTCTTCCAGCACGc                            <  1:523965/140‑1 (MQ=255)                                                                                                          tttAGCGTCTACTGCCATAAACGCCTCGCGAACC‑TC‑CCCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAACAGGAAACCGTCTCTTCCAGCACGCCt                          <  1:305590/140‑1 (MQ=255)                                                                                                           ttAGCGTCTACTGCCATAAACGCCTCGCGAACC‑TC‑CCCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAACAGGAAACCGTCTCTTCCAGCACGCCTg                         >  1:286089/1‑140 (MQ=255)                                                                                                              gCGTCTACTGCCATAAACGCCTCGCGAACCAC‑‑CCCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAACAGGAAACCGTCTCTTCCAGCACGCCTgcgc                      <  1:963764/140‑1 (MQ=255)                                                                                                              gCGTCTACTGCCATAAACGCCTCGCGAACC‑TC‑CCCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAACAGGAAACCGTCTCTTCCAGCACGCCTgcgc                      >  2:307887/1‑140 (MQ=255)                                                                                                              gCGTCTACTGCCATAAACGCCTCGCGAACC‑T‑‑CCCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAACAGGAAACCGTCTCTTCCAGCACGCCTGCGca                     <  1:981476/140‑1 (MQ=255)                                                                                                               cGTCTACTGCCATAAACGCCTCGCGACCC‑TC‑CCCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAACAGGAAACCGTCTCTTCCAGCACGCCTGCGca                     <  1:145797/140‑1 (MQ=255)                                                                                                               cGTCTACTGCCATAAACGCCTCGCGAACC‑TC‑CCCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAACAGGAAACCGTCTCTTCCAGCACGCCTGCGca                     <  2:507807/140‑1 (MQ=255)                                                                                                               cGTCTACTGCCATAAACGCCTCGCGAACC‑TC‑CCCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAACAGGAAACCGTCTCTTCCAGCACGCCTGCGca                     <  1:1066625/140‑1 (MQ=255)                                                                                                               cGTCTACTGCCATAAACGCCTCGCGAACC‑TC‑CCCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAACAGGAAACCGTCTCTTCCAGCACGCCTGCGca                     <  2:286089/140‑1 (MQ=255)                                                                                                               cGTCTACTGCCATAAACGCCTCGCGAACC‑TC‑CCCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAACAGGAAACCGTCTCTTCCAGCACGCCTGCGca                     <  2:28011/140‑1 (MQ=255)                                                                                                               cGTCTACTGCCATAAACGCCTCGCGAACC‑TC‑CCCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAACAGGAAACCGTCTCTTCCAGCACGCCTGCGca                     <  1:761282/140‑1 (MQ=255)                                                                                                                   tACTGCCATAAACGCCTCGCGAACC‑TC‑CCCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAACAGGAAACCGTCTCTTCCAGCACGCCTGCGCACAGt                 >  2:1827/1‑140 (MQ=255)                                                                                                                   tACTGCCATAAACGCCTCGCGAACC‑TC‑CCCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAACAGGAAACCGTCTCTTCCAGCACGCCTGCGCACAGt                 <  1:791184/140‑1 (MQ=255)                                                                                                                        ccATAAACGCCTCGCGAACC‑TC‑CCCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAACAGGAAACCGTCTCTTCCAGCACGCCTGCGCACAGTCCACg            <  1:27714/140‑1 (MQ=255)                                                                                                                                  cTCGCGAACC‑TC‑CCCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAACAGGAAACCGTCTCTTCCAGCACGCCTGCGCACAGTCCACGCCAGCCAGCa  >  1:255924/1‑140 (MQ=255)                                                                                                                                  cTCGCGAACC‑TC‑CCCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAACAGGAAACCGTCTCTTCCAGCACGCCTGCGCACAGTCCACGCCAGCCAGCa  >  2:87222/1‑140 (MQ=255)                                                                                                                                  cTCGCGAACC‑TC‑CCCCCCCCACTTCGAACCCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAACAGGAAACCGTCTCTTCCAGCACGCCTGCGCACAGTCCACGCCAGCCAGCa  >  2:752092/1‑140 (MQ=255)                                                                                                                                              |                                                                                                                                  GCCGCCGAAAATTTGCTCAACACCCACGTTCGCCAGACGCTGCCGGGCAAGCTGATAAATATCCGCCAGATACTTATCACCATGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACC‑T‑‑CCCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAACAGGAAACCGTCTCTTCCAGCACGCCTGCGCACAGTCCACGCCAGCCAGCA  >  NC_000913/2734396‑2734663

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