Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	2,906,607	C→T	12.9%	intergenic (+24/+36)	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,607	0	C	T	12.9%	142.9 / 7.8	53	intergenic (+24/+36)	ygcG/eno	TPM domain protein, putative phosphatase/enolase
Reads supporting (aligned to +/- strand):  ref base C (24/22);  new base T (4/3);  total (28/25)
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.77e-01

AGAAAAGAGCATCTGGTAGAAGTTATCATTCAGACAACAGCGGTTCCGCTGGAGGCTCAGATAGCGGCGGATTTAGCGGCGGGGGCGGTTCTTCTGGAGGTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCC‑CGGAGGCTGGCATTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGATCTCTTTACGACCGTTGTACGGTGCTTTTTCGCCCAGAGCTTCTTCGATACGAATCAGCTGGTTGTATTTAGCAACACGGTCAGAACGGCTCATAG  >  NC_000913/2906459‑2906755                                                                                                                                                      |                                                                                                                                                     aGAAAAGAGCATCTGGTAGAAGTTATCATTCAGACAACAGCGGTTCCGCTGGAGGCTCAGATAGCGGCGGATTTAGCGGCGGGGGCGGTTCTTCTGGAGGTGGCGGCGCATCCGGCCGCTGGTAACTGAACCATATAAAATGCCAGcc‑c                                                                                                                                                      >  4:715387/1‑149 (MQ=255)        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      >  4:422415/1‑149 (MQ=255)                                                                                                                                                 aGCC‑CGGAGGCTGGCATTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGATCTCTTTACGACCGTTGTACGGTGCTTTTTCGCCCAGAGCTTCTTCGATACGAATCAGCTGGTTGTATTTAGCAACACGGTCAGAACGGCTc      >  3:712487/1‑149 (MQ=255)                                                                                                                                                      cGGAGGCTGGCATTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGATCTCTTTACGACCGTTGTACGGTGCTTTTTCGCCCAGAGCTTCTTCGATACGAATCAGCTGGCTGTATTTAGCAACACGGTCAGAACGGCTCATAg  <  3:645909/149‑1 (MQ=255)                                                                                                                                                      |                                                                                                                                                     AGAAAAGAGCATCTGGTAGAAGTTATCATTCAGACAACAGCGGTTCCGCTGGAGGCTCAGATAGCGGCGGATTTAGCGGCGGGGGCGGTTCTTCTGGAGGTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCC‑CGGAGGCTGGCATTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGATCTCTTTACGACCGTTGTACGGTGCTTTTTCGCCCAGAGCTTCTTCGATACGAATCAGCTGGTTGTATTTAGCAACACGGTCAGAACGGCTCATAG  >  NC_000913/2906459‑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: ‑