Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,829,334	G→T	7.3%	P134Q (CCG→CAG)	astD ←	aldehyde dehydrogenase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,829,334	0	G	T	7.3%	135.1 / 3.4	56	P134Q (CCG→CAG)	astD	aldehyde dehydrogenase
Reads supporting (aligned to +/- strand):  ref base G (29/23);  new base T (2/2);  total (31/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 = 4.52e-01

TTACCGCTTCGCCACTCCACGGTGTCAGTTCGCTGGGTTTAAAGATAATGGTGTTACCTGCCAGCAATGCCGGAACGATATGTCCGTTCGGCAAATGACCAGGGAAATTATACGGCCCAAACACCGCCAGCACGCCGTGCGGGCGATGTCGCAGGCTCGCCGCGCCGTCCGGCATTTCACTACGCTGCTCGCCGGTACGAACGTGATACGCCTTAATTGATATCGCGATTTTATTGATCATCGCCGTCACTTCGGTTGCCGCTTCCCAGCGCGGCTT  >  NC_000913/1829194‑1829470                                                                                                                                             |                                                                                                                                         tgaCCGCTTCGCCACTCCACGGTGTCAGTTCGCTGGGTTTAAAGATAATGGTGTTACCTGCCAGCAATGCCGGGACGACATGTCCGTTCGGCAAATGACCAGGGAAATTATACGGCCCAAACACCGCCAGCACGCCGTGCggg                                                                                                                                        <  1:127277/141‑1 (MQ=255)    ccGCTCCGCCACTCCACGCTGTCAGTTCGCTGGGTTTAAAGATATTGGCGTCACCTGCCAGCAATGCCGGAACGATATGTCCGTTCGGCAAACGACCAGGGAAATTATACGGCCCAAACACCGCCCGCACGGCGTGCTGGCgg                                                                                                                                     <  1:894137/143‑2 (MQ=255)     cGCTTCGCCACTCCACGGTGTCAGTTCGCTGGGTTTAAAGATAATGGTGTTACCTGCCAGCAATGCCGGAACGATATGTCCGTTCGGCAAATGACCAGGGAAATTATACGGCCCAAACACCGCCAGCACGCCGTGCGGGCGac                                                                                                                                    >  2:677329/1‑142 (MQ=255)       cTTCGCCACTCGCCGGTGTCAGTTCGCTGGGTTTAAAGATAATGGTGTTACCTGCCAGCAATGCCGGAACGATATGTCCGTTCGGCAAATGACCAGGGAAATTATACGGCCCAACACCCGATAGTACGCCGTGCGGGCGATGt                                                                                                                                  <  2:111286/143‑1 (MQ=255)       cTTCGCCACTCCACGGTGTCAGTTCGCTGGGTTTAAAGATAATGGTGTTACCTGCCAGCAATGCCGGAACGATATGTCCGTTCGGCAAATGACCAGGGAAATTATACGGCCCAAACACCGCCAGCACGCCGTGCGGGCGATGt                                                                                                                                  <  1:50162/143‑1 (MQ=255)          cGCCACTCCACGGTGTCAGTTCGCTGGGTTTAAAGATAATGGTGTTACCTGCCAGCAATGCCGGAACGATATGTCCGTTCGGCAAATGACCAGGGAAATTATACGGCCCAAACACCGCCAGGACGCCGTGCGGGCGATGTCGc                                                                                                                               >  1:163202/1‑143 (MQ=255)          cGCCACTCCACGGTGTCAGTTCGCTGGGTTTAAAGATAATGCCGTTCCCTGCCAGCAATGCCGGAACGATATGTCCGTTCGGCAAATGACCAGGGAAATTATACGGCCCAAACATCGCCCGCACGGCGTGCTGGCGgtgtttg                                                                                                                               <  1:205019/143‑4 (MQ=255)            cctcTCCAGGGTGTCAGTGTGCTGGTGTTAAAGATAATGGTGTTGCCTGCCAGCAATGCCGGAACGATATGTCCGTCCGGCAAATGACCAGGGAAATTATACGGCCCAAACACCGCCAGCACGCCGTGCGGGCGATGTCGCAg                                                                                                                             <  2:859193/140‑1 (MQ=255)                      gtgtCAGTTCGCTTGGTTTCTATATATTGGTGTCACCTGCCAGCAATGCCGGAACGATATGTCCGGGCGGCAAATGACCAGGGAAATTATACGGCCCAAACACCGCCAGCACGCCGTGCGGGCGATGTCGCAGGCTcgccgcg                                                                                                                   <  1:6555/143‑1 (MQ=255)                         tCAGTTCGCTGGGTTTAAAGATAATGGTGTTACCTGCCAGCAATGCCGGAACGATATGTCCGTTCGGCAAATGACCAGGGAAATTATACGGCCCAAACACCGCCAGCACGCCGTGCGGCCGATGTCCCAGGCTcgccccaccc                                                                                                                >  2:774905/1‑137 (MQ=255)                             ttCGCTGGGTTTAAAGATAATGGTGTTACCTGCCAGCAATGCCGGAACGATATGTCCGTTCGGCAAATGACCCGGGAAATTATACGGCCCAAACACCGCCAGCACGCCGTGCGGGCGATGTCGCAGGCTCGCCGCGCCGCCCg                                                                                                            >  2:790173/1‑143 (MQ=255)                             ttCGCTGGGTTTAAAGATAATGGTGTTACCTGCCAGCAATGCCGGAACGATATGTCCGTTCGGCAAATGACCAGGGAAATTATACGGCCCAAACACCGCCAGCACGCCGTGCGGGCGATGTCGCAGGCTCGCCGCGCCGTCCg                                                                                                            >  2:379518/1‑143 (MQ=255)                             ttCGCTGGGTTTAAAGATAATGGTGTTACCTGCCAGCAATGCCGGAACGATATGTCCGTTCGGCAAATGACCAGGGAAATTATACGGCCCAAACACCGCCAGCACGCCGTGCGGGCGATGTCGCAGGCTCGCCGCGCCGTCCg                                                                                                            >  2:124746/1‑143 (MQ=255)                                gCTGGGTTTAAAGATAATGGTGTTACCTGCCAGCAATGCCGGAACGATATGTCCGTTCGGCAACTGACCCGGGAAATTATACGGCCCAAACACCGCCAGCACGCCGTGCGGGCGATGTCGGAGGCTcgccgcggcgccaccca                                                                                                         >  2:628764/1‑136 (MQ=255)                                     gTTGAAAGAGAATGGTGTTACCGGCCAGCAATGCCGGAACGATATGTCCGGTCGGCAAATGACCAGGGAAATTATACGGCCCAAACACCGCCAGCACGCCGTGCGGGCGATGTCGCAGGCTCGCCGCGCCGTCCGGCATTTCa                                                                                                    <  1:893407/143‑1 (MQ=255)                                          aaGATAATGGTGTTACCTGCCAGCAATGCCGGAACGATAGGTCCGTTCGGCAAATGACCAGGGAAATTATACGGACGAAACACCGCCAGCACGCCGTGCGGGCGAGGTCGCAGGCTCGCCGCCCCGTCCGTCATTTCACTACg                                                                                               >  1:430802/1‑143 (MQ=255)                                           aGACAATGGTGTGATTTGCCAGCAATGCCGTAACTATATGTCCGTTCGGCAAATGACCAGGGAAATTATACGGCCCAAACACCGCCAGCACGCCGTGCGGGCGATGTCGCAGGCTCGCCGCGCCGTCCGGCATTTCACTACgc                                                                                              <  2:40560/143‑1 (MQ=255)                                            gATAATGGTGTTACCTGCCAGCAATGCCGGAACGATATGTCCGTTCGGCAAATGACCAGGGAAATTATACGGCCCAAACACCGCCAGCACGCCGTGCGGGCGATGTCCCAGGCTCGCCCCGCCATCCGCCCTTCCACTACgcc                                                                                             >  2:421149/1‑142 (MQ=255)                                             aTAATGGTGTTACATGACAGCAATGCCGGAAAAATAGGTCCGTTCGGCAAATGACCAGGGCAATTATACGGCCCAAACACCGCCAGCACGCCGTGCGGGCGATGACGCAGGCCCGCCGCGCCCGCCGGCACTTCAcgaggcgc                                                                                            >  2:109518/1‑136 (MQ=255)                                                 tGGTGTTACCTGCCAGCAATGCCGGAACGATATGTCCGTTCGGCAAATGCCCCGGGATATTATACGTCCCAAACATCGCCAGAACGCAGTGCGGGCGATGTCGCAGGCTCGCCGCGCCGTCCGGCATTTCACTACGCTGCTCg                                                                                        >  1:131159/1‑143 (MQ=255)                                                     gtgaCCGGCCAGCAATGCCGGAACGATATGTCCGTTCGGCAAATGACCAGGGAAATTATACGGCCCAAACACCGCCAGCACGCCGTGCGGGCGATGTCGCAGGCTCGCCGCGCCGTCCGGCATTTCACTACGCTGCTCGCCgg                                                                                    <  1:685885/140‑1 (MQ=255)                                                     gtTCCCTGCCAGCCATGCCGGAACGATCTGGCCGTTCGGCAAATGACCAGGGAAATTATACGGCCCAAACACCGCCAGCACGCCGTGCGGGCGATGTCGCAGGCTCGCCGCGCCGTCCGGCATTTCACTACGCTGCTCGCCgg                                                                                    <  1:40668/143‑1 (MQ=255)                                                      ttACCTGCCAGCAATGCCGGAACGATATGTCCGTTCGGCAAATGACCAGGGAAATTATACGGCCTAAAAACCACCAGCACGCCGTGCGGGCGATGTCGCAGGCTCGCCGCGCCGTCCGGCATTTCCCTACGCTGCTCGcccgc                                                                                   >  1:126947/1‑140 (MQ=255)                                                      ttACCTGCCAGCAATGCCGGAACGATATGTCCGTTCGGCAAATGACCAGGGAAATTATACGGCCCACACACCGCCAGCACGCCGTGCGGGCGATGTCGCAGGCTCGCCGCGCCGTCCGGCATTTCACTACGCTGCTCGCCGGt                                                                                   <  1:866662/143‑1 (MQ=255)                                                        aCCTGCCAGCAATGCCGGAACGATATGTCCGTTCGGCAAATGACCAGGGAAATTATACGGCCCAAACACCGCCAGCACGCCGTGCGGGCGATGTCGCAGGCTCGCCGCGCCGTCCGGCCTTTCCCTACCCCGCTCGCCGGTAc                                                                                 >  1:228728/1‑143 (MQ=255)                                                           tGCCAGCAATGCCCGACCGACATGTCCGTTCGGCAAATGACCAGGGAAATTATACGGCCCAAACACCGCCAGCACGCCGTGCGGGCGATGTCGCAGGCTCGCCGCGCCGTCCGGCATTTCACTACGCTGCTCGCCGGTACGaa                                                                              <  2:162167/143‑1 (MQ=255)                                                            gCCAGCAATGCCGGAACGATATGTCCGTTCGGCAAATGACCAGGGAAATTATACGGACCAAACACCGCCAGCACGCCGTGCGGGCGATGTCGCAGGCTCGCCGCGCCGTCCGGCATTTGACCACACTACTCGCCGGTACGacc                                                                             >  2:806088/1‑141 (MQ=255)                                                                      ccGGAACGATATGTCCGTTCGGCAAATGACCAGGGAAATTATACGGCCCAAACACCGCCAGCACGCCGTGCGGGCGATGTCGCAGGCTCGCCGCGCCGTCCGGCATTTCACTACGCTGCTCGCCGGTACGAACGTGATCCGcc                                                                   >  1:794458/1‑143 (MQ=255)                                                                       cGGAACGATATGTCCGTTCGGCAAATGACCAGGGAAATTATACGGCCCAAACACCGCCAGCACGCAGTGCGGGCGATGTCGCAGGCTCGCCGCGCCGTCCGGCATTTCACTACGCTGCTCGCCGGTGCGAACGGGATAGGCCt                                                                  >  2:129387/1‑143 (MQ=255)                                                                        ggAACGATATGTCCGTTCGGCAAATGACCAGGGAAATTATACGGCCCAAACACCGCCAGCACGCCGTGCGGGCGATGTCGCAGGCTCGCCGCGCCGTCCGGCATTTCACTACGCTGCTCGCCGGTGCGAACGTGGTACGCCtt                                                                 >  1:425043/1‑143 (MQ=255)                                                                         gAACGATATGTCCGTTCGGCAAATGACCAGGGAAATTATACGGCCCAAACACCGCCAGCACGCCGTGCGGGCGATGTCGCAGGCTCGCCGCGCCGTCCGGCATTTCACTACGAGGCGCGCAGGTACGAACGTGGTACGCCtga                                                                >  1:632225/1‑141 (MQ=255)                                                                          aaCGATATGTCCGTTCGGCAAATGACCAGGGAAATTATACGGCCCAAACACCGCCAGCACGCCGTGCGGGCGATGTCGCAGGCTCGCCGCGCCGTCCGGCATTTCACTACGCTGCTCGCCGGTACGAACGTGATACGCCTTaa                                                               <  1:80794/143‑1 (MQ=255)                                                                            cGATATGTCCGCTCGGCAAATGACCAAAAGCAGACGTCTGCCCAAGCACCGCCAGCACGCCGTGCGGGCGATGCCGCAGTCTCGCCGCGCCGTCCGGCATTTCACTACGCCGCTCGTCGGTACGAACGTGATACGCCTTAAtt                                                             <  1:164880/143‑1 (MQ=255)                                                                            cGATATGCTCGTGCGGCAAATGACCAGGGAAATTATACGGCCCAAACACCGCCAGCACGCGGGGCGGGCGAGGTCGCAGGCGCGCCGCGCCGGTCGGCATTTCACTACGCTGCTCGTCGGTACGAACGTGATACGCCTTAAtt                                                             <  1:843132/143‑1 (MQ=255)                                                                               tatGTCCGTTCGGCAAATGACCAGGGAAATTATACGGCCCAAACACCGCCAGCACGCCGTGCGGGCGATGTCGCAGGCTCGCCGCGCCGTCCGGCATTTCACTACGCTGCTCGCCGGTACGAACGTGATACGCCTTAATTGat                                                          >  2:856934/1‑143 (MQ=255)                                                                                   tCCGTTCGGCAAATGACCAGGGAAATTATACGGCCCAAACACCGCCAGCACGCCGTGCGGGCGATGTCGCAGGCTCGCCGCGCCGTCCGGCATTTCACTACGCTGCTCGCCGGTACGAACGTGATACGCCTTAATTGATATcg                                                      <  1:629239/143‑1 (MQ=255)                                                                                     cGTTCGGCAAATGACCAGGGAAATTATACGGCCCAAACACCGCCAGCACGCCGTGCGGGCGATGTCGCAGGCTCGCCGCGCCGTCCGGCATTTCACTACGCTGCTCGCCGGTACGAACGTGATACGCCTTCCTTGCTATcgcg                                                    >  1:289976/1‑143 (MQ=255)                                                                                      gTTCGGCAAATGACCAGGGAAATTATACGGCCCAAACACCGCCAGCACGCCGTGCGGGCGTTGTCGCAGGCTCGCCGCGCCGTTCAGCATTTCATTACTCTGCTCGCCGGTACGAACGTGATACGCCTTCATTGATATAGCGa                                                   >  1:146394/1‑143 (MQ=255)                                                                                       ttCGGCAAATGACCAGGGAAATTATACGGCCCAAACACCGCCAGCACGCCGTGCGGGCGATGTCGCAGGCTCGCCGCGCCGTCCGGCATTTCACTACGCTGCTCGCCGGTACGAACGTGATACGCCTTAATTGATATCGCGAt                                                  >  1:731778/1‑143 (MQ=255)                                                                                          ggCAAATGACCAGGGAAATTATACGGCCCAAACACCGCCAGCACGCCGTGCGGGCGATGTCGCAGGCTCGCCGCGCCGTCCGGCATTTCACTACGCTGCTCGCCGGTACGAACGTGATACGCCTTAATTGATATCGCGAtttt                                               >  1:634387/1‑143 (MQ=255)                                                                                             aaaTGACCGGGGAAATTATACGGCCCAAACACCGCCAGCACGCCGTGCGGGCGATGTCGCAGGCGCGCCGCGCCGTCCGGCATTTCACTACGCTGCTCGCCGGTACGAACGTGATACGCCTTAATTGATATCGCGATGTTAtt                                            <  2:48510/143‑1 (MQ=255)                                                                                             aaaTGACCGGGAACATGACACGGCCCACACACCGCCAGCACGCCGTGCGGGCGATGCCGCAGGCTCGCCGCGCCGTCCGGCATTTCACTACGCTGCTCGCCTGTACGAACGTGATACGCCTTAATTGATATCGCGATTTTAtt                                            <  2:348740/143‑1 (MQ=255)                                                                                             aaaTGACCAGGTAATTTATCCGGCCCAATCACCGCCAGCACGCCGTGCGGGCCAGTTCGCAGGCTCGCCGCGCCGTCCGTCATTTCTCTACGCTGCTCGCCGGTACGAACTTGATACGCCTTAATTGATATCGCGATTTTAtt                                            <  2:751411/143‑1 (MQ=255)                                                                                                     aGGGAAATTATACGGCCCAAACACCGCCAGCACGCCGTGCGGGCGATGTCGCAGGCTCGCCGCGCCGTCCGGCATTTCACTACGCTGCTCGCCGGTACGAACGTGATACGCCTTAATTGATATCGCGATTTTCTTGATCATCg                                    >  2:558685/1‑143 (MQ=255)                                                                                                                aCGGCCCAAACACCGCCAGCACGCCGTGCGGGCGATGTCGCAGGCTCGCCGCGCCGTCCGGCATTTCACTACGCTGCTCGCCGGTACGAACGTGATACGCCTTAATTGATATCGCGAGTTTATTGATCATCGCCGTCACGTCg                         >  1:496722/1‑143 (MQ=255)                                                                                                                 cGGCCCAAACACCGCCAGCACGCCGTGCGGGCGATGTCGCAGGCTCGCCGCGCCGTCCGGCATTTCACTACGCTGCTCGCCGGTACGAACGTGATACGCCTTAATTGATATCGCGATTTTATTGATCATCGCCGTCACTTCgg                        <  1:718497/143‑1 (MQ=255)                                                                                                                    cccAGACCCCGCCAGCACGCCGTGCGGGCGATGTCGCAGGCTCGCCGCGCCGTCCGGCATTTCACTACGCTGCTCGCCGGTACGAACGTGATACGCCTTAATTGATATCGCGATTTTATTGATCATCGCCGTCACTTCGGTTg                     <  1:687967/143‑1 (MQ=255)                                                                                                                     tcatacacCGCCAGCACGCTGTTGGGGCGATGTCGCAGGCTCGCCGCGCCGTCCGGCATTTCACTACGCTGCTCGCCGGTATGAACGTGATACGCCTTAATTGATATCGCGATTTTATTGATCATCGCCGTCACTTCGGTTGc                    <  2:563501/139‑1 (MQ=255)                                                                                                                        agcacCGCCAGCACGCCGGGCGGGCGAGGTCGCAGGCTCGCCGCGCCGTCCGGCATTTCACTACGCTGCTCGCCGGTACGAACGTGATACGCCTTAATTGATATCGCGATTTTATTGATCATCGCCGTCACTTCGGTTGCCGc                 <  1:348956/141‑1 (MQ=255)                                                                                                                        aaCACCGCCAGCACGCCGTGCGGGCGATGTCGCAGGCTCGCCGCGCCGTCCGGCATTTCACTACGCTGCTCGCCGGTACGAACTGGATACGCCTTACTTGATATCGCGATTCTCTTGATCATCGCCGTCACTACGGTAGccac                 >  2:448082/1‑141 (MQ=255)                                                                                                                            tcGCCCGCACGGCGTGCTGGCGGTGTCGCAGGCTCGCCGCGCCGTCCGGCATTTCACTACGCTGCTCGCCGGTACGAACGTGATACGCCTTAATTGATATCGCGATTTTATTGATCATCGCCGTCACGTCGGTTGCCGCTTcc             >  2:671138/2‑143 (MQ=255)                                                                                                                            tcGCCCGCACGGCGTGCTGGCGATGTCGCAGGCTCGCCGCGCCGTCCGGCATTTCACTACGCCGCTCGCCGGTACGAACGCGATACGCCTTAATTGATATCGCGACTTTATTGCTCATCGCCGTCACCTCGGTTGCCGCTTcc             >  2:345286/2‑143 (MQ=255)                                                                                                                            ccGCCAGCACGCCGTGCGGGCGATGTCGCAGGCTCGCCGCGCCGTCCGGCATTTCACTACGCTGCTCGCCGGTACGAACGTGATACGCCTTAATTGATATCGCGATTTTATTGATCATCGCCGTCACTTCGGTTGCCGCTTcc             <  2:699971/143‑1 (MQ=255)                                                                                                                                      gCCGTGCGGGCGATGTCGCAGGCTCGCCGCGCCGTCCGGCATTTCACTACGCTGCTCGCCGGTACGAACGTGATACGCCTTAATTGATATCGCGATTTTATTGATCATCGCCGTCACTTCGGTTGCCGCCTCCCAGCGCGGCt   >  2:777971/1‑143 (MQ=255)                                                                                                                                      gCCGTGCGGGCGATGTCGCAGGCTCGCCGCGCCGTCCGGCATTTCACTACGCTGCTCGCCGGGACGAACGTGATACGCCTTAATTGATATCGCGATTTTATTGCTCATCACCGTCACCTCGGTTGCCGCTTCCCAGTGCGGCt   >  2:305978/1‑143 (MQ=255)                                                                                                                                       ccGTGCGGGCGATGTCGCAGGCTCGCCGCGCCGTCCGGCATTTCACTACGCTGCTCGCCGGTACGAACGTGATACGCCTTAATTGATATCGCGATTTTATTGATCATCGCCGTCACTTCGGTTGCCGCTTCCCAGCGCGGCtt  <  1:124746/143‑1 (MQ=255)                                                                                                                                       ccGTGCGGGCGATGTCGCAGGCTCGCCGCGCCGTCCGGCATTTCACTACGCTCCTCGCCGGCACGCCCGCGATACGCCTTCCTTGATATCGCGCTCTTCTTGACCATCGCCGCCACTTCGGCTGCCGCTTCCTCGCGCCGCtt  >  1:354959/1‑143 (MQ=255)                                                                                                                                             |                                                                                                                                         TTACCGCTTCGCCACTCCACGGTGTCAGTTCGCTGGGTTTAAAGATAATGGTGTTACCTGCCAGCAATGCCGGAACGATATGTCCGTTCGGCAAATGACCAGGGAAATTATACGGCCCAAACACCGCCAGCACGCCGTGCGGGCGATGTCGCAGGCTCGCCGCGCCGTCCGGCATTTCACTACGCTGCTCGCCGGTACGAACGTGATACGCCTTAATTGATATCGCGATTTTATTGATCATCGCCGTCACTTCGGTTGCCGCTTCCCAGCGCGGCTT  >  NC_000913/1829194‑1829470

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