Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,829,322	A→C	7.5%	L138R (CTG→CGG)	astD ←	aldehyde dehydrogenase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,829,322	0	A	C	7.5%	121.5 / 2.3	55	L138R (CTG→CGG)	astD	aldehyde dehydrogenase
Reads supporting (aligned to +/- strand):  ref base A (28/22);  new base C (3/2);  total (31/24)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.08e-01

CCATAAACGCATTACCGCTTCGCCACTCCACGGTGTCAGTTCGCTGGGTTTAAAGATAATGGTGTTACCTGCCAGCAATGCCGGAACGATATGTCCGTTCGGCAAATGACCAGGGAAATTATACGGCCCAAACACCGCCAGCACGCCGTGCGGGCGATGTCGCAGGCTCGCCGCGCCGTCCGGCATTTCACTACGCTGCTCGCCGGTACGAACGTGATACGCCTTAATTGATATCGCGATTTTATTGATCATCGCCGTCACTTCGGTTGCCGCTTCCCAGC  >  NC_000913/1829183‑1829463                                                                                                                                            |                                                                                                                                              ccATAAACGCATTACCGCTTCGCCACTCCACGGTGTCAGTTCGCTGGGTATAAAGATAAGGGTGTTACCTGCCAGCAATGCCGGAACGATATGGCCGTTCGGCAAGTGACCAGGGAAATTATACGGCCCCACCACCGCCAGcg                                                                                                                                            >  1:370034/1‑142 (MQ=255)       aCGCATTACCGCTTCGCCACTCCACGGTGTCAGTTCGCTGGGTTTAAAGATAATGGTGTTACCTGCCAGCAATGCCGGAACGATATGTCCGTTCGGCAAATGACCAGGGAAATTATACTGCCCACACACCGCCAGCACGCCGt                                                                                                                                      >  1:1001380/1‑143 (MQ=255)       aCGCATTACCGCTTCGCCACTCCACGGTGTCAGTTCGCTGGGTTTAAACATAATGGTGTTACCTGCCAGCAATGCCGGAACGATATGTCCGTTCGGCAACGGACCAGGGAAATTATACGGCCCAAACAACGCCAGCACGCCgc                                                                                                                                      >  2:724427/1‑142 (MQ=255)        cGCATTACCGCTTCGCCACTCCACGGTGTCAGTTCGCTAGGTTTAAAGCCAATGGTGTTCCCTGCCAGCAATGCCGGAACGATATGTCCGTCCGGCAAATGACCAGGTAAATTATACGGCCCAAACACCGCCAGCACGCCGTg                                                                                                                                     <  1:37332/143‑1 (MQ=255)         gCATTACCGCTTCGCCACTCCACGGTGTCAGTGCGCTGGGTTTAAAGATAATGGTGTTACCTGCCAGCAATGCCGGAACGATATGTCCGTTCGGCAAATGACCAGGGAAATTATACGGCCCAAACACCGCCAGCACGCCGTGc                                                                                                                                    <  1:550636/143‑1 (MQ=255)             tACCGCTTCGCCACTCCACGGTGTCAGTTCGCTGGGTTTAAAGATAATGGTGTTACCTGCCAGCAATGCCGGAACGATATGTCCGTTCGGCAAATGACCAGGGAAATTATACGGCCCAAACACCGCCAGCACGCCGTGCGGGc                                                                                                                                >  1:788415/1‑143 (MQ=255)             gcccGCTGCGCCACTCGACGGTGTCAGTTCGCTGGGTGTAAAGATAGTGGTCTTACCTGCCAGCAATGCCGGAACGATATGTCCGGTCGGCAAATGACCAGGGAAATTATACGGCCCAAACACCGCCAGCACGCCGTGCGGGc                                                                                                                                <  2:1013727/141‑1 (MQ=255)              aCCGCCTCCGCACTGCACGGTCTCAGTTCGCGGGGTTTAAAGATAATGGGGGTACCTGCCAGCAATGCCGGAACGATATGTCCGTTCGGCAAATGACCAGGGAAATTATACGGCCCAAACACCGCCCGCACGGCGTGCTGGCg                                                                                                                               <  2:204370/143‑1 (MQ=255)                 gCTTCGCCACTCCACGGTGTCAGTTCGCTGGGTTTAAAGATAATGGTGTTACCTGCCAGCAATGCCGGAACGATATGTCCGTTCGGCAAATGACCAGGGCAATTATACGGCCCAAACACCGCCAGCACGCCGCGCGGGCCATg                                                                                                                            >  2:408162/1‑143 (MQ=255)                   ttCGCCACTCCACGGTGTCAGTTCGCTGGGTTTAAAGATAATGGTGTTACCTGCCAGCAATGCCGGAACGATATGTGCGTTCGGCAAATGACCAGGGAAATTATACGGCCCAACCACCCGCCGCACGCCGTGCGGGCGATGTc                                                                                                                          >  2:318075/1‑143 (MQ=255)                   ttCGCCACTCCACGGTGTCAGTTCGCTGGGTTTAAAGATAATGGTGTTACCTGCCAGCAATGCCGGAACGATATGTCCGTTCGGCAAATGACCAGGGAAATTATACGGCCCAAACATCGCCCGCACGGCGTGCTGGCGgtgtt                                                                                                                          <  2:834482/143‑2 (MQ=255)                      gCCACTCCACGGTGTCAGTTCGCTGGGTTTAAAGATAATGGTGTTACCTGCCAGCAATGCCGGAACGATATGTCCGTTCGGCAAATGACCAGGGAAATTATACGGCCCAAACACCGCCAGCACGCCGTGCGGGCGATGTCGCa                                                                                                                       >  2:221706/1‑143 (MQ=255)                       ccACTCCACGGTGTCAGTTCGCTGGGTTTAAAGATAATGGTGTTACCTGCCAGCAATGCCGGAACGATATGTCCGTTCGGCAAATGACCAGGGAAATTATACGGCCCAAACACCGCCAGCACGCCGTGCGGGCGATGTCGCAg                                                                                                                      >  1:154915/1‑143 (MQ=255)                         aCTCCACGGTGTCAGTTCGCTGGGTTTAAAGATAATGGTGTTACCTGCCAGCAATGCCGGAACGATATGTCCGTTCGGCAAATGACCAGGGAAATTATACGGCCCAAACACCGCAAGCACGCCGTGCGGGCCATGTCGCAGGc                                                                                                                    >  2:310367/1‑143 (MQ=255)                          cTCCACGGTGTCAGTTCGCTGGGTTTAACGATAATGGTGTTACCTGCCAGCAATGCCGGAACGATATGTCCGTTCGGCAAATGACCAGGGAAATTATACGGCCCAAACACCGCCAGCACGCCGTGCGGGCGATGTCGCAGGCt                                                                                                                   <  1:588476/143‑1 (MQ=255)                                ggTGTCAGTTCGCTGGCTTTAACGATAATGGTGTTACCTGCCAGCAATGCCGGAACGATATGTCCGTTCGGCAAATGACCAGGGAAATTATACGGCCCAAACACCGCCAGCACGCCGTGCGGGCGATGTCGCAGGCTcgccgc                                                                                                             <  1:681416/143‑1 (MQ=255)                                     cAGTTCGCTGGGTTTAAAGATAATGGTGTTACCTGCCAGCAATGCCGGAACGATATGTCCGTTCGGCAAATGAGCAGGGAAATTATACGGCCCAAACACCGCCAGCACGCCGTGCGGGCGGTGTCGCAGGGTGGCCGCGGCGt                                                                                                        >  1:92960/1‑143 (MQ=255)                                      aGTTCGCTGGGTTTAAAGATAATGGTGTTACCTGCCAGCAATGCCGGAACGATATGTCCGTTCGGCAAATGACCAGGGAAATTATACGGCCCAAACACCGCCAGCACGCCGTGCGGGCGATGTCGCAGGCTCGCCGCGCCGTc                                                                                                       <  1:837160/143‑1 (MQ=255)                                        ttCGCTGGGTTTAAAGATAATGGTGTTACCTGCCAGCAATGCCGGAACGATATGTCCGTTCGGCAAATGACCAGGGAAATTATACGGCCCAAACACCGCCAGCGCGCCGTGCGGGCGATGTCGCAGCCTCGCCGCGCCGTCCg                                                                                                     >  2:1013617/1‑143 (MQ=255)                                        ttCGCTGGGTTTAAAGATAATGGTGTTACCTGCCAGCAATGCCGGAACGATATGTCCGTTCGGCAAATGACCAGGGAAATTATACGGCCCAAACACCGCCAGCACGCCGTGCGGGCGATGTCGCAGGCTCGCCGCgccggccg                                                                                                     >  1:70150/1‑143 (MQ=255)                                              gggTTTAAAGATAATGGTGTTACCTGCCAGCAATGCCGGAACGATATGTCCGTTCGGCAAATGACCAGGGAAATTATACGGCCCAAACACCGCCAGCACGCCGTGCGGGCGATGTCGCAGGCTCGCCGCGCCGTCCGGCAttt                                                                                               <  1:720647/143‑1 (MQ=255)                                              gggTTTAAAGATAATGGTGTTACCTGCCAGCAATGCCGGAACGATATGTCCGTTCGGCAAATGACCAGGGAAATTATACGGCCCAAACACCGCCAGCACGCCGTGCGGGCGATGTCGCAGGCTCGCCGCGCCGTCCGGCAttt                                                                                               >  1:925351/1‑143 (MQ=255)                                              ggcgttAAAGATAATGGTGTTACCTGCCAGCAATGCCGGAACGATATGTCCGTTCGGCAAATGACCAGGGAAATTATACGGCCCAAACACCGCCAGCACGCCGTGCGGGCGATGTCGCAGGCTCGCCGCGCCGTCCGGCAttt                                                                                               <  1:1064881/139‑1 (MQ=255)                                                    aaaGATAATGGTGTTACCTGCCAGCAATGCCGGAACGATATGTCCGTTCGGCAAATGACCAGGGAAATTATACGGCCCAAACACCGCCAGCACGCCGTGCGGGCGATGTCGCAGGCTCGCCGCGCCGTCCGGCATTTCACtcc                                                                                         >  1:686285/1‑141 (MQ=255)                                                    aaaGATAATGGTGTTACCTGCCAGCAATGCCGGAACGATATGTCCGTTCGGCAAATGACCAGGGAAATTATACGGCCCAAACACCGCCAGCACGCCGTGCGGGCGATGTCGCAGGCTCGCCGCGCCGTCCGGCATTTCACTAc                                                                                         >  1:200311/1‑143 (MQ=255)                                                      cgATAATGGGGTTACCTGCCAGCAATGCCGGAACGATATGTCCGTTCGGCAAATGACCAGGGAAATTATACGGCCCAAACACCGCCAGCACGCCGTGCGGGCGATGTCGCAGGCTCGCCGCGCCGTCCGGCATTTCACTACgc                                                                                       <  1:873435/142‑1 (MQ=255)                                                          aaTGGTGTTACCTGCCAGCAATGCCGGAACGACATGTCCGTTCGGCAAATGACCAGGGAAATTATACGGCCCAAACACCGCCAGCACGCCGTGCGGGCGATGTCGCAGGCTCGCCGCGCCGTCCGGCATTTCACTACgccgct                                                                                   >  2:826566/1‑143 (MQ=255)                                                                     cTGCCAGCAATGCCGGAACGATATGTCCGTTCGGCAAATGACCAGGGAAATTATATGCCCCAAACACCGCCAGCACGCCGTGCGGGCGATGTCGCAGGCTCGCCGCGCCGTCAGGCATTGCACTACGCCGCTCGCCGGGACGa                                                                        >  2:176295/1‑143 (MQ=255)                                                                      tGCCAGCAATGCCGGAACGATATGTCCGTTCGGCGAATGACCAGGGAAATTATACGGCCCAAACACCGCCAGCACGCCGTGCGGGCGATGTCGCAGGCGCGCCGCGGCGTCCGGCATTTCACTACGGCGCGCTCCGGTGCgga                                                                       >  1:566307/1‑141 (MQ=255)                                                                       gCCAGCAATGCCGGAACGATATGTCCGTTCGGCAAATGACCAGGGCAATTATACGGCCCAAACACCGCCAGCACGCCGTGCGGGCTATGTCGCAGGCTGGCCGGGCCGTCCGGCGTTTCACTATCCTGTACGCCGGTAcaacc                                                                      >  2:200232/1‑139 (MQ=255)                                                                        ccAGCAATGCCGGAACGATATGTCCGTTCGGCAAATGACCAGGGAAATTATACGGCCCAAACACCGCCAGCACGCCGTGCGGGCGATGTCGCAGGCTCGCCGCGCCGTCCGGCATTTCACTACGCTGCTCGCCGGTACgcccg                                                                     >  2:1043719/1‑139 (MQ=255)                                                                         cAGCAATGCCGGAACGATATGTCCCTTCGGCAAATGACCAGGGAAATTATACGGCCCAACCACCGCCAGCACGCCGTGCGGGCGATGTCGCAGGCTCGCCGCGCCGTCCGGCATTTCACTACGCTGCTCGCCGGTACGAACGt                                                                    <  1:351995/143‑1 (MQ=255)                                                                                 ccGGAACGATATGTCCGTTCGGCAAATGACCAGGGAAATTATACGGCCCAAACACCGCCAGCACGCCGTGCGGGCGATGTCGCAGGCTCGCCGCGCCGTCCGGCATTTCACTACGCTGCTCGCCGGTACGCACGTGATACGcc                                                            >  1:977964/1‑143 (MQ=255)                                                                                    gAACGATATGTCCGTTCGGCAAATGACCAGGGAAATTATACGGCCCAAACACCGCCAGCACGCCGTGCGGGCGATGTCGCAGGCTCGCCGCGCCGTCCGGCATTTCACTACGCTGCTCGCCGGTACGAACGTGATACGCCttc                                                         >  2:859218/1‑142 (MQ=255)                                                                                    gAACGATATGTCCGTTCGGCAAATGACCAGGGAAATCATACGGCCCAAACACCGCCAGCACGCCGTGCGGGCGATCTCGCCGGCTCGCCGCGCCGTCCGGCATTGCACTACGCTGCGCGCCGGCACGACCGTGCTCCCCCttc                                                         >  2:422757/1‑142 (MQ=255)                                                                                       cgacacgTCCGTTCGGCAAATGACCAGGGAAATTATACGGCCCAAACACCGCCAGCACGCCGTGCGGGCGATGTCGCAGGCTCGCCGCGCCGTCCGGCATTTCACTACGCTGCTCGCCGGTACGAACGTGATACGCCCTAAtt                                                      <  2:885163/137‑1 (MQ=255)                                                                                       cGATATGTCCGTTCGGCAAATGACCAGGGAAATTATACGGCCCAAACACCGCCAGCACGCCGTGCGGGCGATGTCGCAGGCTCGCCGCGCCGTCCGGCATTTCACTACGCTGCTCGCCGGTACGCACGTGATACGCCTTAAtt                                                      >  1:793303/1‑143 (MQ=255)                                                                                       cGATATGTCCGTTCGGCAAATGACCAGGGAAATTATACGGCCCAAACACCGCCAGCACGCCGTGCGGGCGATGTCGCAGGCTCGCCGCGCCGTCCGGCATTTCACTACGCTGCTAGCCGGTACGAACGTGATACGCCTTAAtt                                                      >  2:903338/1‑143 (MQ=255)                                                                                          tatGGCCGTTCGGCAAATGACCAGGGAAATTATACGGCCCAAACACCGCCAGCACGCCGTGCGGGCGATGTCGCAGGCTCGCCGCGCCGTCCGGCATTTCACTACGCTGCTCGCCGGTACGAACGTGATACGCCTTAATTGat                                                   <  2:547147/143‑1 (MQ=255)                                                                                             gTCCGTTCGGCAAATGACCAGGGAAATTATACGGCCCAAACACCGCCAGCACGCCGTGCGGGCGATGTCGCAGGCTCGCCGCGCCGTCCGGCATTTCACTACGCTGCTCGCCGGTACGAACGTGCTACGCCTTAGTTGATATc                                                >  1:964234/1‑143 (MQ=255)                                                                                                cGTTCGGCAAATGACCAGGGAAATTATACGGCCCAAACACCGCCAGCACGCCGTGCGGGCGATGTCGCAGGCTCGCCGCGCCGTCCGGCATTTCACTACGCTGCTCGCCGGTACGAACGCGATACGCCTTCATTCATATcgcg                                             >  2:238861/1‑143 (MQ=255)                                                                                                   tCGGCAAATAACCAGGGAAATTATACGGCCCAAACACCGCCCGCACGCCGTGCGCGCGATGTCGCAGGCGCGCCGCGCCGTCAGGCATTTCACTACGCTGCTCGCCGGCACGAACCTGCTACGCCTTCACTGAGATcgcgctc                                          >  2:236303/1‑140 (MQ=255)                                                                                                     ggCAAATGACCAGGGAAATTATACGGCCCAAACAACAACAACAGGCCTTTTGGGCGATGTCGCAGGCTCGCCGCGCCGTCCGGCATTTCACTACGCTGCTCGCCGGTACGAACGTGATACGCCTTAATTGATATCGCGAtctt                                        >  1:729443/1‑143 (MQ=255)                                                                                                      gCAAATGACCAGGGAAATTATACGGCCCAAACACCGCCAGCACGCCGTGCGGGCGATGTCGCAGGCTCGCCGCGCCGTCCGGCATTTCACTACGCTGCTCGCCGGTACGAACGTGATACGCCCTAATTGATACCGCGAtttaa                                       >  2:888332/1‑141 (MQ=255)                                                                                                       cAAATGACCAGGGAAATTATACGGCCCAAACACCGCCAGCACGCCGTGCGGGCGATGTCGCAGGCTCGCCGCGCCGTCCGGCATTTCACTACGCTGCTCGCCGGTACGAACGTGATACGCCTTAATTGATATCGCGATTTTAt                                      <  1:256404/143‑1 (MQ=255)                                                                                                        aaaTGACCAGGGAAATTATACGGCCCAAACCCCGCCAGCACGCCGTGCGGGCGATGTCGCAGGCTCGCCGCGCCGTCCGGCATTTCACTACGCTGCTCGCCGGTACGAACGTGATACGCCTTAATTGATATCGCGATTTTAtt                                     <  2:437421/143‑1 (MQ=255)                                                                                                            gACCCGGGAAATTATACGGCCCAAACACCGCCAGCACGCCGTGCGGGCGATGTCGCAGGCTCGCCGCGCCGTCCGGCATTTCACTACGCTGCTCGCCGGTACGAACGTGATACGCCTTAATTGATATCGCGATTTTATTGatc                                 <  2:772042/143‑1 (MQ=255)                                                                                                                aGGGAAGTTAGACGGCCCAAACACCGCCAGCACGCCGTGCGGGCGATGTCGCAGGCTCGCCGCGCCGTCCGGCATTTCACTACGCTGCTCGCCGGTACGAACGTGATACGCCTTAATTGATATCGCGATTTTATTGATCATCg                             <  2:100666/143‑1 (MQ=255)                                                                                                                  ggAACTTATACGGCCCAATCACCGCCAGAACGCGGTGCGGGCGATGTCGCAGGCTCGCCGCGCCGTCCGGCATTTCACTACGCTGCTCGCCGGTACGAACGTGATACGCCTTAATTGATATCGCGATTTTATTGATCATCGcc                           <  2:260376/143‑1 (MQ=255)                                                                                                                              gCCCAACCTCCGCCAGCACGCCGTGCGGGCGTTGTCGCAGGCTCGCCGCGCCGTCCGGCATTTCACTACGCTGCTCGCCGGTACGAACGTGATACGCCTTAATTGATATCGCGATTTTATTGATCATCGCCGTCACTTCGGtt               <  1:599195/143‑1 (MQ=255)                                                                                                                               cccAAACACCGCCAGCACGCCGCGCGGGCGATGTCGCAGGCTCGCCGCGCCGTCCGGCATTTCACTACGCTGCTCGCCGGTACGGACGTGATACGCCTTAATTGATATCGCGATTTTATTGATCATCGCCGTCACTTCGGTTg              <  2:526779/143‑1 (MQ=255)                                                                                                                                   gacatcGCCCGCACGGCGTGCTGGCGATGTCGCAGGCTCGCCGCGCCGTCCGGCATTTCACTACGCTGCTCGCCGGTACGAACGTGATACGCCTTAATTCATATCGCGATTCTATTGATCATCGCCGGCACTTCGGTTGCCGc          >  2:664831/6‑143 (MQ=255)                                                                                                                                     ccccGCCAGCACGCCGCGCGGGCGATGTCGCAGGCTCGCCGAGGTGGGCAGCATTGCACTACGCTGCTCGCCGGTACGAACGTGATACGCCTTAATTGATATCGCGATTTTATTGATCATCGCCGTCACTTCGGTTGCCGCtt        <  1:724427/141‑1 (MQ=255)                                                                                                                                     cacCGCCAGCACGCCGTGCGGGCGGTGGCGCAGGCTCGCCGCGCCGTCCGGCATTTCACTACGCTGCTCGCCGGTACGAACGTGATACGCCTTAATTGATATCGCGATTTTATTGATCATCGCCGTCACTTCGGTTGCCGCtt        <  2:855896/143‑1 (MQ=255)                                                                                                                                           cAGCACGCCGTGCGGGCGATGTCGCAGGCTCGCCGCGCCGCCCGGCATTTCACTACGCTGCTCGCCGGTACGAACGTGATACGCCTTAATTGATATCGCGATTTTATTGATCATCGCCGTCACTTCGGTTGCCGCTTCCCAgc  <  1:881237/143‑1 (MQ=255)                                                                                                                                            |                                                                                                                                              CCATAAACGCATTACCGCTTCGCCACTCCACGGTGTCAGTTCGCTGGGTTTAAAGATAATGGTGTTACCTGCCAGCAATGCCGGAACGATATGTCCGTTCGGCAAATGACCAGGGAAATTATACGGCCCAAACACCGCCAGCACGCCGTGCGGGCGATGTCGCAGGCTCGCCGCGCCGTCCGGCATTTCACTACGCTGCTCGCCGGTACGAACGTGATACGCCTTAATTGATATCGCGATTTTATTGATCATCGCCGTCACTTCGGTTGCCGCTTCCCAGC  >  NC_000913/1829183‑1829463

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