Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	2,304,988	T→C	13.1%	intergenic (+595/+120)	eco → / ← mqo	ecotin, a serine protease inhibitor/malate dehydrogenase, FAD/NAD(P)‑binding domain

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,304,988	0	T	C	13.1%	138.7 / 12.6	56	intergenic (+595/+120)	eco/mqo	ecotin, a serine protease inhibitor/malate dehydrogenase, FAD/NAD(P)‑binding domain
Reads supporting (aligned to +/- strand):  ref base T (29/18);  new base C (2/5);  total (31/25)
Fisher's exact test for biased strand distribution p-value = 1.22e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.37e-01

CCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCAACATCTAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGCATACGCCACATCCGGCATACCATGCCGGATGTGGCGTATCATTACAACGCAATATCCGCCACTTCTTTT  >  NC_000913/2304866‑2305135                                                                                                                           |                                                                                                                                                    cccGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCaa                                                                                                                           >  7:221141/1‑149 (MQ=31)            gATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCAACATCTAACGcc                                                                                                                <  3:163857/149‑1 (MQ=35)                  gCGGTTACGCCGCATCCGACATCTAACGCCCCAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTGGGATAAGGCGGTTACGCCGCATCCAACATCTAACGCCCGAGcc                                                                                                          <  6:61266/149‑1 (MQ=255)                  gCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCAACATCTAACGCCCGAGcc                                                                                                          >  8:261584/1‑149 (MQ=37)                   cGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCAACATCTAACGCCCGAGCCg                                                                                                         <  7:261584/149‑1 (MQ=37)                               aTCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACCCCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCAACATCTAACGCCCGAGCCGGTTGCCTGATGc                                                                                             >  5:356874/1‑149 (MQ=37)                                  cGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCAACATCTAACGCCCGAGCCGGTTGCCTGATGCGAc                                                                                          <  7:47245/149‑1 (MQ=37)                                   gACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGAGTCTTATCAGGCATCCACCGATGAGAAGTCCTCAACCCCGTAGGTTCGATAAGGCGGTAGCGACGGAACCAACATCTATATCCCCAGCCGGATGACTGATGCGCCg                                                                                         >  6:263386/1‑149 (MQ=2)                                       tCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCAAc                                                                                                                          <  4:65585/112‑1 (MQ=32)                                       tCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCAAc                                                                                                                          >  3:65585/1‑112 (MQ=32)                                         tAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGTCGGTTACGCCGCATCCAACATCTAACGCCCGAGCCGGTTGCCTGAAGCGACGCTGgcg                                                                                   >  1:79590/1‑149 (MQ=255)                                          aaCGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCAACATCTAACGCCCGAGCCGGTTGCCTGATGCGACGCTGgcgc                                                                                  >  7:215184/1‑149 (MQ=37)                                          aaCGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCAACATCTAACGCCCGAGCCGGTTGCCTGATGCGACGCTGgcgc                                                                                  >  6:129941/1‑149 (MQ=37)                                           aCGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCAACATCTAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGcgcg                                                                                 >  4:362589/1‑149 (MQ=37)                                             gCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCAACATCTAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTc                                                                               >  6:30157/1‑149 (MQ=37)                                               ccAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTATGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCAACATCTAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCt                                                                              >  7:296571/1‑148 (MQ=37)                                                     cGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCAACATCTAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAg                                                                        >  1:316087/1‑148 (MQ=37)                                                          gCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCAACATCTAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGCATACg                                                                  >  8:15294/1‑149 (MQ=255)                                                           ccTGATGCGACGCTGGCGCGTCTTATCCGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCAACATCTAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGCATACGc                                                                 <  1:208457/149‑1 (MQ=255)                                                            cTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCAACATCTAACGCCCGAGCCGGTTGCCTGAt                                                                                               <  4:182069/118‑1 (MQ=39)                                                             tGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCAACATCTAACGCCCGAGCCGGTTGc                                                                                                    >  6:4328/1‑112 (MQ=39)                                                             tGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCAACATCTAACGCCCGAGCCGGTTGc                                                                                                    <  5:4328/112‑1 (MQ=39)                                                              gATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCAACATCTAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCGTATCAGCATCCGCca                                                               >  4:290476/1‑148 (MQ=255)                                                                tGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGACGCATCCAACATCAAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGCATACGCCACAt                                                            >  2:355280/1‑149 (MQ=255)                                                                          gcgcGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCAACATCTAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGCATACGCCACATCCGGCATAc                                                   >  7:228581/1‑148 (MQ=255)                                                                          gcgcGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCAACATCTAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGCATACGCCACATCCGGCATAc                                                   >  7:8096/1‑148 (MQ=255)                                                                               tCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCAACATCTAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGCATACGCCACATCCGGCATACCATGcc                                             >  7:75832/1‑149 (MQ=255)                                                                                    tCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCAACATCTAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGCATACGCCACATCCGGCATACCATGCCGGAt                                         <  5:129941/148‑1 (MQ=255)                                                                                     cAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCAACATCTAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGcgcg                                                                                 <  4:102274/107‑1 (MQ=17)                                                                                      aGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCAACATCTAACGCCCGAGCCGGTTGCCTGa                                                                                                <  2:222418/91‑1 (MQ=40)                                                                                      aGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCAACATCTAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGCATACGCCACATCCGGCATACCATGCCGGAtgtg                                      >  1:277213/1‑149 (MQ=255)                                                                                      aGGCCTACACCGCTGTGAAGGCCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCAACATCTAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGCATACGCCACATCCGGCATACCATGCCGGAtgtg                                      <  5:200294/149‑1 (MQ=255)                                                                                      aGGCCTACACCGCTGTGAACTGCGCCACCCCGTAGGATGGATAAGGCGGTTACGCCGCATCCAACATCTAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGCATACGCCACATCCGGCATACCATGCCGGAtgtg                                      <  5:267444/149‑1 (MQ=255)                                                                                        gCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCAACATCTAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGCATACGCCACATCCGGCATGGTATGCCGGATGTGGc                                    <  3:362589/149‑1 (MQ=255)                                                                                        gCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCAACATCTAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGCATACGCCACATCCGGCATACCATGCCGGATGTGGc                                    >  4:254595/1‑149 (MQ=255)                                                                                        gCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCAACATCTAACGCCCGAGCCGGTTGCCTGAt                                                                                               <  4:6670/90‑1 (MQ=18)                                                                                        gCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCAACATCTAACGCCCGAGCCGGTTGCCTGAt                                                                                               >  3:6670/1‑90 (MQ=18)                                                                                          cTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCAACATCTAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGCATACGCCACATCCGGCATACCATGCCGGATGTGGCg                                   >  2:338625/1‑148 (MQ=255)                                                                                          cTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCAACATCTAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGCATACGCCACATCCGGCATACCATGCCGGATGTGGCg                                   >  2:230235/1‑148 (MQ=255)                                                                                          cTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCAACATCTAACGCCCGAGCCGGTTGc                                                                                                    <  8:158971/83‑1 (MQ=18)                                                                                            acacCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCAACAt                                                                                                                        <  3:5221/61‑1 (MQ=35)                                                                                            acacCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCAACAt                                                                                                                        >  4:5221/1‑61 (MQ=35)                                                                                             cacCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCAACATCTAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGCATACGCCACATCCGGCATACCATGCCGGATGTGGCGTATc                               <  8:75832/149‑1 (MQ=255)                                                                                             cacCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCAACATCTAACGCCCGAGCCGGTTCCCTGATGCGACTCTGGCGCGTCTTATCAGCATACGCCACATCCGGCATACCATGCCGGATGTGGCGTAt                                >  4:179781/1‑148 (MQ=255)                                                                                                  cTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCAACATCTAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGt                                                                                <  2:156442/95‑1 (MQ=40)                                                                                                  cTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCAACATCTAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGt                                                                                >  1:156442/1‑95 (MQ=40)                                                                                                        agcgCTCCACCCCGTATCGAGGATAAGGCGCTTACGCCGCATCCCACATCTAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGCATACGCCACATCCGGCATACCATGCCGGATGTGGCGTATCATTACAACGCa                    <  2:390894/146‑1 (MQ=21)                                                                                                         gTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCAACATCTAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGCATACGCCACATCCGGCATACCATGCCGGATGTGGCGTATCATTACAACGCaa                   >  6:135962/1‑149 (MQ=255)                                                                                                         gTGCTCCACCCCGTAGGGTGGATAAGGCGGTTACGCCGCATCCAACATCTAACGCCCGAGCCGGTTGCCTGAAGCGACGCTGGCGCGTCTTATCAGCATACGCCACATCCGGCATACCATGCCGGATGTGGCGTATCATTACAACGCaa                   <  3:290476/149‑1 (MQ=255)                                                                                                            cTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCAACATCTAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTc                                                                               <  6:65790/86‑1 (MQ=18)                                                                                                            cTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCAACATCTAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTa                                                                            <  6:228214/89‑1 (MQ=18)                                                                                                            cTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCAACATCTAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTa                                                                            >  5:228214/1‑89 (MQ=18)                                                                                                              ccACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCAACATCTAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGCATACGCCACATCCGGCATACCATGCCGGATGTGGCGTATCATTACAACGCAATATcc              <  5:30157/149‑1 (MQ=255)                                                                                                               cACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCAACATCTAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGCATACGCCACATCCGGCATACCATGCCGGATGTGGCGTATCATTACAACGCAATATcc              <  7:15294/148‑1 (MQ=255)                                                                                                                  cccGTAGGTTGGATAAGGCGGTTACGCCGCATCCAACATCTAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGCATACGCCACATCCGGCATACCATGCCGGATGTGGCGTATCATTACAACGCAATATCCGCCa          >  2:168879/1‑149 (MQ=255)                                                                                                                         gTTGGATAAGGCGGTTACGCCGCATCCAACATCTAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGCATACGCCACATCCGGCATACCATGCCGGATGTGGCGTATCATTACAACGCAATATCCGCCACTGCttt   >  6:73367/1‑149 (MQ=255)                                                                                                                          ttGGATAAGGCGGTTACGCCGCATGCAACATCTAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGCATACGCCACATCCGGCATACCATGCCGGATGTGGCGTATCATTACAACGCAATATCCGCCACTTCtttt  <  5:137307/149‑1 (MQ=255)                                                                                                                           |                                                                                                                                                    CCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCAACATCTAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGCATACGCCACATCCGGCATACCATGCCGGATGTGGCGTATCATTACAACGCAATATCCGCCACTTCTTTT  >  NC_000913/2304866‑2305135

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