Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	2,906,609	G→C	13.0%	intergenic (+26/+34)	ygcG → / ← eno	TPM domain protein, putative phosphatase/enolase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,906,609	0	G	C	13.0%	156.5 / 13.2	54	intergenic (+26/+34)	ygcG/eno	TPM domain protein, putative phosphatase/enolase
Reads supporting (aligned to +/- strand):  ref base G (24/23);  new base C (4/3);  total (28/26)
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.52e-01

AGCATCTGGTAGAAGTTATCATTCAGACAACAGCGGTTCCGCTGGAGGCTCAGATAGCGGCGGATTTAGCGGCGGGGGCGGTTCTTCTGGAGGTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCC‑CGGAGGCTGGCATTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGATCTCTTTACGACCGTTGTACGGTGCTTTTTCGCCCAGAGCTTCTTCGATACGAATCAGCTGGTTGTATTTAGCAACACGGTCAGAACGGCTCATAGAA  >  NC_000913/2906466‑2906757                                                                                                                                                 |                                                                                                                                                     agCATCTGGTAGAAGTTATCATTCAGACAACAGCGGTTCCGCTGGAGGCTCAGATAGCGGCGGATTTAGCGGCGGGGGGGGTTTTTTTGGCGGGGGGGGCGCGTCCGGCCGCTGGTAACTCAACCATAAAAAATGGCAGCC‑CGGAGGCt                                                                                                                                                 >  4:693158/1‑149 (MQ=255)  gCATCTGGTAGAAGTTATCATTCAGACAACAGCGGTTCCGCTGGAGGCTCAGATAGCGGCGGATTTAGCGGCGGGGGCGGTTCTTCTGGAGGTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCC‑CGGAGGCTg                                                                                                                                                >  3:56952/1‑149 (MQ=255)    aTCTGGTAGAAGTTATCATTCAGACAACAGCGGTTCCGCTGGAGGCTCAGATAGCGGCGGATTTAGCGGCGGGGGCGGTTCTTCTGGAGGTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCC‑CGGAGGCTGGc                                                                                                                                              >  4:440290/1‑149 (MQ=255)       tGGTAGAATTTATCCTTCAGAGAACAGCGGTTCCGCTGGAGGCTCAGATAGCGGCGGATTTAGCGGCGGGCGCGTTTCTTCTGGAGTTGGCGTCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCC‑TCCGGGCTGGCAtt                                                                                                                                           <  3:468829/149‑1 (MQ=255)         gTAGAAGTTATCATTCAGACAACAGCGGTTCCGCTGGAGGCTCAGATAGCGGCGGATTTAGCGGCGGGGGCGGTGCTTCTGGCGGTGGCAGCGCATCCCGGCGCTGGGAACTCCACCATAAAAAATGCGGGGC‑CGGAGGGTTGCCtttt                                                                                                                                         >  2:36874/1‑149 (MQ=255)          tAGAAGTTATCATTCAGACAACAGCGGTTCCGCTGGAGGCTCAGATAGCGGCGGATTTAGCGGCGGGGGCGGTTCTTCTGGAGGTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCC‑CGGAGGCTGGCAttttt                                                                                                                                        >  3:445396/1‑149 (MQ=255)          tAGAAGTTATCATTCAGACAACAGCGGTTCCGCTGGAGGCTCAGATAGCGGCGGATTTAGCGGCGGGGGCGGCGTTTCTGGAGGTGGCGGCGCATCCGGGCGCTGGTAAATCAACCATAAAAAATACCAGCC‑CCGAGGCTGGCAttttt                                                                                                                                        >  3:701870/1‑149 (MQ=255)           aGAAGTTATCATTCAGACAACAGCGGTTCCGCTGGAGGCTCAGATAGCGGCGGATTTAGCGGCGGGGGCGGTTCTTCTGGAGGCGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATAGCAGCC‑CGGAGGCTGGCATTTTTa                                                                                                                                       >  3:580766/1‑149 (MQ=255)              aGTTATCATTCAGACAACAGCGGTTCCGCTGGAGGCTCAGATAGCGGCGGATTTAGCGGCGGGGGCGGTTCTTCTGGAGGTGGCGGTGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCC‑CGGAGGCTGGCATTTTTAAAt                                                                                                                                    <  1:586117/149‑1 (MQ=255)                  aTCATTCAGACAACAGCGGTTCCGCTGGAGGCTCAGATAGCGGCGGATTTAGCGGCGGGGGCGGTTCTTCTGGAGGTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCC‑CGGAGGCTGGCATTTTTAAATCAGa                                                                                                                                <  3:440290/149‑1 (MQ=255)                      ttCAGACAACAGCGGTTCCGCTGGAGGCTCAGATAGCGGCGGATTTAGCGGCGGGGGCGGTTCTTCTGGAGGTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCC‑CGGGGGGTGGCATTTTTAAATCAGATaaa                                                                                                                            >  3:481237/1‑149 (MQ=255)                      ttCAGACAACAGCGGTTCCGCTGGAGGCTCAGATAGCGGCGGATTTAGCGGCGGGGGCGGCTCTTCTGGGGGTGGCGGCGCATCCGGCCGCTGGTCAATCAACCATAAACAATACCAGGC‑CCGAGGCTGGCATTTTTAAATCAGATaaa                                                                                                                            >  4:565663/1‑149 (MQ=255)                           acaacaGCGGTTCCGCTGGAGGCTCAGATAGCGGCGGATTTAGCGGCGGGGGCGGTTCTTCTGGAGGTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCC‑CGGAGGCTGGCATTTTTAAATCAGATAAagtcag                                                                                                                       <  1:61025/149‑1 (MQ=255)                             aacaGCGGTTCCGCTGGAGGCTCAGATAGCGGCGGATTTAGCGGCGGGGGCGGTTTTTCTGGAGGTGGAGGCGCATCCGGCCGCTGGTTACTCAACCATAAAAAATGCCAGCC‑CGGAGGCTGGCATTTTTAAATCAGATAAagtcagtc                                                                                                                     >  4:470981/1‑149 (MQ=255)                              acaGCGGTTCCGCTGGAGGCTCAGATAGCGGCGGATTTAGCGGCGGGGGCGGTTCTTCTGGAGGTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCC‑CGGAGGCTGGCATTTTTAAATCAGATAAAGTCAGTCt                                                                                                                    >  1:669558/1‑149 (MQ=255)                                   ggTTCCGCTGGAGGCTCAGATAGCGGCGGATTTAGCGGCGGGGGCGGTTCTTCTGGAGGTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCC‑CGGAGGCTGGCATTTTTAAATCAGATAAAGTCAGTCTTATGc                                                                                                               <  1:228173/149‑1 (MQ=255)                                   ggTTCCGCTGGAGGCTCAGATAGCGGCGGATTTAGCGGCGGGGGCGGTTCTGCTGGAGGTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCC‑TCCGGGCTGGCAtttt                                                                                                                                         <  1:668937/123‑1 (MQ=255)                                   ggTTCCGCTGGAGGCTCAGATAGCGGCGGATTTAGCGGCGGGGGCGGTTCCTCTGGAGGTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCC‑TCCGGGCTGGCAtttt                                                                                                                                         >  2:668937/1‑123 (MQ=255)                                      tCCGCTGGAGGCTCAGATAGCGGCGGATTTAGCGGCGGGGGCGGTTCTTCTGGAGGTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCC‑CGGAGGCTGGCATTTTTAAATCAGATAAAGTCAGTCTTATGCCTg                                                                                                            <  1:487335/149‑1 (MQ=255)                                           tGGAGGCTCAGATAGCGGCGGATTTAGCGGCGGGGGCGGTTCTTCTGGAGGTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCC‑CGGAGGCTGGCATTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCtt                                                                                                       <  1:221444/149‑1 (MQ=255)                                            ggAGGCTCAGGTAGCGGCGGATTTAGCGGCGGGGGTGGTTCTTCTGGAGGTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCC‑CGGAGGCTGGCATTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCttt                                                                                                      <  1:36874/149‑1 (MQ=255)                                               ggCTCAGATAGCGGCGGATTTAGCGGCGGGGGCGGTTCTTCTGGAGGTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCC‑CGGAGGCTGGCATTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGa                                                                                                    <  2:137899/148‑1 (MQ=255)                                               ggCTCAGATAGCGGCGGATTTAGCGGCGGGGGCGGTTCTTCTGGAGGTGGCGGCGCATCCGGCCGCTGGGAACTCAACCATAAAAAATGCCAGCC‑CGGAGGCTGGCATTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGa                                                                                                    >  1:137899/1‑148 (MQ=255)                                                  tCATATAGCGGCGGATTTAGCGGCGGGGGCGGTTCTTCTGGAGGTGGCGGGGCGTCCGGGCGCTGGGAACGCAACCATAAAAAATGCCAGCC‑CGGGAGCCGGCGTTTTTAAATCAGATAAAGTCAGTCCTATGCCGGGGGTTTGAtctc                                                                                                >  4:645909/1‑149 (MQ=255)                                                  tCAGATAGCGGCGGATTTAGCGGCGGGGGCGGTTCTTCTGGAGGTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCC‑CGGAGGCTGGCATTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTTAtctc                                                                                                >  1:376767/1‑149 (MQ=255)                                                  tCAGATAGCGGCGGATTTAGCGGCGGGGGCGGTTCTTCTGGAGGTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCC‑CGGAGGCTGGCATTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGAtctc                                                                                                <  1:745523/149‑1 (MQ=255)                                                              ggATTTAGCGGCGGGGGCGGTTCTTCTGGAGGTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCC‑CGGAGGCTGGCATTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGATCTCTTTACGACCGtt                                                                                    <  1:576430/149‑1 (MQ=255)                                                                     gcggcgGGGGCGGTTCTTCTGGAGGTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCC‑CGGAGGCTGGCATTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGATCTCTTTACGACCGTTGTACGGt                                                                             <  2:669558/149‑1 (MQ=255)                                                                             ggCGGTTCTTCTGGAGGTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCC‑CGGAGGCTGGCATTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGATCTCTTTACGACCGTTGTACGGTGCTTTTTc                                                                     <  1:60980/149‑1 (MQ=255)                                                                                ggTTCTTCTGGAGGTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCC‑TCCGGGCTGGCAtttt                                                                                                                                         >  4:192173/1‑78 (MQ=25)                                                                                ggTTCTTCTGGAGGTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCC‑TCCGGGCTGGCAtttt                                                                                                                                         <  3:192173/78‑1 (MQ=37)                                                                                ggTTCTTCTGGAGGTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCC‑CGGAGGCTGGCATTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGATCTCTTTACGACCGTTGTACGGTGCttttt                                                                      <  1:305856/145‑1 (MQ=255)                                                                                ggTTCTTCTGGAGGTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCC‑CGGAGGCTGGCATTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGATCTCTTTACGACCGTTGTACGGTGCttttt                                                                      >  2:305856/1‑145 (MQ=255)                                                                                  ttcttcTGGAGGTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCC‑CGGAGGCTGGCATTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGATCTCTTTACGACCGTTGTACGGTGCTTTTTCGCCCa                                                                >  2:442099/1‑149 (MQ=255)                                                                                           aGGTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCC‑CGGAGGCTGGCATTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGATCTCTTTACGACCGTTGTACGGTGCTTTTTCGCCCAGAGcttctt                                                       <  3:211099/149‑1 (MQ=255)                                                                                            ggTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATATAAAATGCCAGCC‑CGGAGGCTGGCATTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGATCTCTTTACGACCGTTGTACGGTGCTTTTTCGCCCAGAGCttcttc                                                      <  3:715387/149‑1 (MQ=255)                                                                                                 cggcgCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCC‑CGGAGGCTGGCATTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGATCTCTTTACGACCGTTGTACGGTGCTTTTTCGCCCAGAGCTTCTTCGATAc                                                 <  4:580766/149‑1 (MQ=255)                                                                                                   gcgcATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCC‑CGGAGGCTGGCATTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGATCTCTTTACGACCGTTGTACGGTGCTTTTTCGCCCaga                                                              <  1:126668/134‑1 (MQ=255)                                                                                                   gcgcATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCC‑CGGAGGCTGGCATTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGATCTCTTTACGACCGTTGTACGGTGCTTTTTCGCCCaga                                                              >  2:126668/1‑134 (MQ=255)                                                                                                              cGCTGGTAACTCAACCATAAAAAATGCCAGCC‑CGGAGGCTGGCATTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGATCTCTTTACGACCGTTGTACGGTGCTTTTTCGCCCAGAGCTTCTTCGATACGAATCAGCTGGtt                                    <  2:845/149‑1 (MQ=255)                                                                                                                   gTAACTCAACCATAAAAAATGCCAGCC‑CGGAGGCTGGCATTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGATCTCTTTACGACCGTTGTACGGTGCTTTTTCGCCCAGAGCTTCTTCGATACGAATCAGCTGGTTGTAtt                               <  2:259992/149‑1 (MQ=255)                                                                                                                     aaCTCAACCATAAAAAATGCCAGCC‑CGGAGGCTGGCATTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGATCTCTTTACGACCGTTGTACGGTGCTTTTTCGCCCAGAGCTTCTTCGATACGAATCAGCTGGTTGTATTTa                             <  2:398484/149‑1 (MQ=255)                                                                                                                        tCAACCATAAAAAATGCCAGCC‑CGGAGGCTGGCATTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGATCTCTTTACGACCGTTGTACGGTGCTTTTTCGCCCAGAGCTTCTTCGATACGAATCAGCTGGTTGTATTTAGCa                          >  3:136491/1‑149 (MQ=255)                                                                                                                        tCAACCATAAAAAATGCCAGCC‑CGGAGGCTGGCATTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGATCTCTTTACGACCGTTGTACGGTGCTTTTTCGCCCAGAGCTTCTTCGATACGAATCAGCTGGTTGTATTTAGCa                          >  4:84859/1‑149 (MQ=255)                                                                                                                            gatttaaaaaTGCCAGCCTCCG‑GGCTGGCATTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGATCTCTTTACGACCGTTGTACGGTGCTTTTTCGCCCAGAGcttc                                                         >  3:254347/6‑114 (MQ=37)                                                                                                                            gatttaaaaaTGCCAGCCTCCG‑GGCTGGCATTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGATCTCTTTACGACCGTTGTACGGTGCTTTTTCGCCCAGAGcttc                                                         <  4:254347/109‑1 (MQ=37)                                                                                                                            gatttaaaaaTGCCAGCC‑TCCGGGCTGGCATTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGATCTCTTTACGACCGTTGTACGGTGCTTTTTCGCCCAGAGCTTCTTCGATACGAATCAGCTGGTTGTATTTAGCAacac                      >  1:489629/6‑149 (MQ=39)                                                                                                                             cATAAAAAATGCCAGCC‑CGGAGGCTGGCATTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGATCTCTTTACGACCGTTGTACGGTGCTTTTTCGCCCAGAGCTTCTTCGATACGAATCAGCTGGTTGTATTTAGCAACACg                     >  1:522298/1‑149 (MQ=255)                                                                                                                                   aaaTGCCAGCC‑CGGAGGCTGGCATTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGATCTCTTTACGACCGTTGTACGGTGCTTTTTCGCCCAGAGCTTCTTCGATACGAATCAGCTGGTTGTATTTAGCAACACGGTCAGa               >  2:135475/1‑149 (MQ=255)                                                                                                                                    aaTGCCAGCC‑TCCGGGCTGGCATTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGATCTCTTTACGACCGTTGTACGGTGCTTTTTCGCCCAGAGCTTCTTCGATACGAATCAGCTGGTTGTATTTAGCAACACGGTCAGaa              >  1:548271/1‑149 (MQ=255)                                                                                                                                         cAGCC‑CGGAGGCTGGCATTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGATCTCTTTACGACCGTTGTACGGTGCTTTTTCGCCCAGAGCTTCTTCGATACGAATCAGCTGGTTGTATTTAGCAACACGGTCAGAACGGCt         <  3:470981/149‑1 (MQ=255)                                                                                                                                          aGCC‑CGGAGGCTGGCATTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGATCTCTTTACGACCGTTGTACGGTGCTTTTTCGCCCAGAGCTTCTTCGATACGAATCAGCTGGTTGTATTTAGCAACACGGTCAGAACGGCTc        >  3:712487/1‑149 (MQ=255)                                                                                                                                          aGCC‑CGGAGGCTGGCATTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGATCTCTTTACGACCGTTGTACGGTGCTTTTTCGCCCAGAGCTTCTTCGATACGAATCAGCTGGTTGTATTTAGCAACACGGTCAGAACGGCTc        >  4:422415/1‑149 (MQ=255)                                                                                                                                               cGGAGGCTGGCATTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGATCTCTTTACGACCGTTGTACGGTGCTTTTTCGCCCAGAGCTTCTTCGATACGAATCAGCTGGCTGTATTTAGCAACACGGTCAGAACGGCTCATAg    <  3:645909/149‑1 (MQ=255)                                                                                                                                                 gAGGCTGGCATTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGATCTCTTTACGACCGTTGTACGGTGCTTTTTCGCCCAGAGCTTCTTCGATACGAATCAGCTGGTTGTATTTAGCAACACGGTCAGAACGGCTCATAGaa  <  3:636898/149‑1 (MQ=255)                                                                                                                                                 gAGGCTGGCATTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGATCTCTTTACGACCGTTGTACGGTGCTTTTTCGCCCAGAGCTTCTTCGATACGAATCAGCTGGTTGTATTTAGCAACACGGTCAGAACGGCTCATAGaa  <  3:445168/149‑1 (MQ=255)                                                                                                                                                 |                                                                                                                                                     AGCATCTGGTAGAAGTTATCATTCAGACAACAGCGGTTCCGCTGGAGGCTCAGATAGCGGCGGATTTAGCGGCGGGGGCGGTTCTTCTGGAGGTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCC‑CGGAGGCTGGCATTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGATCTCTTTACGACCGTTGTACGGTGCTTTTTCGCCCAGAGCTTCTTCGATACGAATCAGCTGGTTGTATTTAGCAACACGGTCAGAACGGCTCATAGAA  >  NC_000913/2906466‑2906757

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