Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	2,906,608	G→C	20.2%	intergenic (+25/+35)	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,608	0	G	C	20.2%	124.3 / 24.8	54	intergenic (+25/+35)	ygcG/eno	TPM domain protein, putative phosphatase/enolase
Reads supporting (aligned to +/- strand):  ref base G (21/22);  new base C (7/4);  total (28/26)
Fisher's exact test for biased strand distribution p-value = 5.05e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.11e-01

AGCATCTGGTAGAAGTTATCATTCAGACAACAGCGGTTCCGCTGGAGGCTCAGATAGCGGCGGATTTAGCGGCGGGGGCGGTTCTTCTGGAGGTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCC‑CGGAGGCTGGCATTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGATCTCTTTACGACCGTTGTACGGTGCTTTTTCGCCCAGAGCTTCTTCGATACGAATCAGCTGGTTGTATTTAGCAACACGGTCAGAACGGCTCATAG  >  NC_000913/2906466‑2906755                                                                                                                                                |                                                                                                                                                    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                                                                                                                                       <  3:192173/78‑1 (MQ=37)                                                                                ggTTCTTCTGGAGGTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCC‑TCCGGGCTGGCAtttt                                                                                                                                       >  4:192173/1‑78 (MQ=25)                                                                                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                        >  4:84859/1‑149 (MQ=255)                                                                                                                        tCAACCATAAAAAATGCCAGCC‑CGGAGGCTGGCATTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGATCTCTTTACGACCGTTGTACGGTGCTTTTTCGCCCAGAGCTTCTTCGATACGAATCAGCTGGTTGTATTTAGCa                        >  3:136491/1‑149 (MQ=255)                                                                                                                            gatttaaaaaTGCCAGCCTCCG‑GGCTGGCATTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGATCTCTTTACGACCGTTGTACGGTGCTTTTTCGCCCAGAGcttc                                                       <  4:254347/109‑1 (MQ=37)                                                                                                                            gatttaaaaaTGCCAGCCTCCG‑GGCTGGCATTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGATCTCTTTACGACCGTTGTACGGTGCTTTTTCGCCCAGAGcttc                                                       >  3:254347/6‑114 (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)                                                                                                                                                |                                                                                                                                                    AGCATCTGGTAGAAGTTATCATTCAGACAACAGCGGTTCCGCTGGAGGCTCAGATAGCGGCGGATTTAGCGGCGGGGGCGGTTCTTCTGGAGGTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCC‑CGGAGGCTGGCATTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGATCTCTTTACGACCGTTGTACGGTGCTTTTTCGCCCAGAGCTTCTTCGATACGAATCAGCTGGTTGTATTTAGCAACACGGTCAGAACGGCTCATAG  >  NC_000913/2906466‑2906755

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