Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	756,886	T→A	L327Q (CTG→CAG)	sdhA →	succinate dehydrogenase, flavoprotein subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	756,886	0	T	A	96.2%	272.0 / 0.1	78	L327Q (CTG→CAG)	sdhA	succinate dehydrogenase, flavoprotein subunit
Reads supporting (aligned to +/- strand):  ref base T (2/1);  new base A (41/34);  total (43/35)
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.99e-01

GATCGAAATCCGTGAAGGTCGCGGCTGTGATGGTCCGTGGGGGCCACACGCGAAACTGAAACTCGATCACCTGGGTAAAGAAGTTCTCGAATCCCGTCTGCCGGGTATCCTGGAGCTTTCCCGTACCTTCGCTCACGTCGATCCGGTGAAAGAGCCGATTCCGGTTATCCCAACCTGTCACTACATGATGGGCG  >  NC_000913/756788‑756981                                                                                                   |                                                                                                gATCGAAATCCGTGAAGGTCGCGGCTGTGATGGTCCGTGGGGGCCACACGCGAAACTGAAACTCGATCACCTGGGTAAAGAAGTTCTCGAATCCCGTCTg                                                                                                >  2:1898302/1‑100 (MQ=255)   tCGAAATCCGTGAAGGTCGCGGCTGTGATGGTCCGTGGGGGCCACACGCGAAACTGAAACTCGATCACCTGGGTAAAGAAGTTCTCGAATCCCGTCAGcc                                                                                              >  1:1252927/1‑100 (MQ=255)     gAAATCCGTGAAGGTCGCGGCTGTGATGGTCCGTGGGGGCCACACGCGAAACTGAAACTCGATCACCTGGGTAAAGAAGTTCTCGAATCCCGTCAGCCgg                                                                                            >  1:1530725/1‑100 (MQ=255)     gAAATCCGTGAAGGTCGCGGCTGTGATGGTCCGTGGGGGCCACACGCGAAACTGAAACTCGATCACCTGGGTAAAGAAGTTCTCGAATCCCGTCAGCCgg                                                                                            >  1:510371/1‑100 (MQ=255)     gAAATCCGTGAAGGTCGCGGCTGTGATGGTCCGTGGGGGCCACACGCGAAACTGAAACTCGATCACCTGGGTAAAGAAGTTCTCGAATCCCGTCAGCCgg                                                                                            <  1:446389/100‑1 (MQ=255)      aaaTCCGTGAAGGTCGCGGCTGTGATGGTCCGTGGGGGCCACACGCGAAACTGAAACTCGATCACCTGGGTAAAGAAGTTCTCGAATCCCGTCAGCCggg                                                                                           >  1:1116228/1‑100 (MQ=255)      aaaTCCGTGAAGGTCGCGGCTGTGATGGTCCGTGGGGGCCACACGCGAAACTGAAACTCGATCACCTGGGTAAAGAAGTTCTCGAATCCCGTCAGCCggg                                                                                           <  1:1818185/100‑1 (MQ=255)        aTCCGTGAAGGTCGCGGCTGTGATGGTCCGTGGGGGCCACACGCGAAACTGAAACTCGATCACCTGGGTAAAGAAGTTCTCGAATCCCGTCAGCCGGGTa                                                                                         <  1:472887/100‑1 (MQ=255)           cGTGAAGGTCGCGGCTGTGATGGTCCGTGGGGGCCACACGCGAAACTGAAACTCGATCACCTGGGTAAAGAAGTTCTCGAATCCCGTCAGCCGGGTATcc                                                                                      >  1:1312410/1‑100 (MQ=255)           cGTGAAGGTCGCGGCTGTGATGGTCCGTGGGGGCCACACGCGAAACTGAAACTCGATCACCTCGGTAAAGAAGTTCTCGAATCCCGTCAGCCGGGTATcc                                                                                      <  2:1252924/100‑1 (MQ=255)             tGAAGGTCGCGGCTGTGATGGTCCGTGGGGGCCACACGCGAAACTGAAACTCGATCACCTGGGTAAAGAAGTTCTCGAATCCCGTCAGCCGGGTATCCTg                                                                                    <  1:1734027/100‑1 (MQ=255)              gAAGGTCGCGGCTGTGATGGTCCGTGGGGGCCACACGCGAAACTGAAACTCGATCACCTGGGTAAAGAAGTTCTCGAATCCCGTCAGCCGGGTATCCTgg                                                                                   >  1:2168989/1‑100 (MQ=255)               aaGGTCGCGGCTGTGATGGTCCGTGGGGGCCACACGCGAAACTGAAACTCGATCACCTGGGTAAAGAAGTTCTCGAATCCCGTCAGCCGGGTATCCTGGa                                                                                  <  2:2117416/100‑1 (MQ=255)                aGGTCGCGGCTGTGATGGTCCGTGGGGGCCACACGCGAAACTGAAACTCGATCACCTGGGTAAAGAAGTTCTCGAATCCCGTCAGCCGGGTATCCTGGAg                                                                                 <  1:2051579/100‑1 (MQ=255)                aGGTCGCGGCTGTGATGGTCCGTGGGGGCCACACGCGAAACTGAAACTCGATCACCTGGGTAAAGAAGTTCTCGAATCCCGTCAGCCGGGTATCCTGGAg                                                                                 <  1:137843/100‑1 (MQ=255)                aGGTCGCGGCTGTGATGGTCCGTGGGGGCCACACGCGAAACTGAAACTCGATCACCTGGGTAAAGAAGTTCTCGAATCCCGTCAGCCGGGTATCCTGGAg                                                                                 <  1:736350/100‑1 (MQ=255)                       ggCTGTGATGGTCCGTGGGGGCCACACGCGAAACTGAAACTCGATCACCTGGGTAAAGAAGTTCTCGAATCCCGTCAGCCGGGTATCCTGGAGCTTTccc                                                                          <  2:412804/100‑1 (MQ=255)                         cTGTGATGGTCCGTGGGGGCCACACGCGAAACTGAAACTCGATCACCTGGGTAAAGAAGTTCTCGAATCCCGTCAGCCGGGTATCCTGGAGCTTTCCCGt                                                                        <  1:830833/100‑1 (MQ=255)                            tgATGGTCCGTGGGGGCCACACGCGAAACTGAAACTCGATCACCTGGGTAAAGAAGTTCTCGAATCCCGTCAGCCGGGTATCCTGGAGCTTTCCCGTAcc                                                                     >  2:1143554/1‑100 (MQ=255)                               tGGTCCGTGGGGGCCACACGCGAAACTGAAACTCGATCACCTGGGTAAAGAAGTTCTCGAATCCCGTCAGCCGGGTATCCTGGAGCTTTCCCGTACCTTc                                                                  >  2:2076649/1‑100 (MQ=255)                                 gTCCGTGGGGGCCACACGCGAAACTGAAACTCGATCACCTGGGTGAAGAAGTTCTCGAATCCCGTCAGCCGGGTATCCTGGAGCTTTCCCGTACCTTCGc                                                                >  1:1048632/1‑100 (MQ=255)                                 gTCCGTGGGGGCCACACGCGAAACTGAAACTCGATCACCTGGGTAAAGAAGTTCTCGAATCCCGTCAGCCGGGTATCCTGGAGCTTTCCCGTACCTTCGc                                                                >  1:1207964/1‑100 (MQ=255)                                 gTCCGTGGGGGCCACACGCGAAACTGAAACTCGATCACCTGGGTAAAGAAGTTCTCGAATCCCGTCAGCCGGGTATCCTGGAGCTTTCCCGTACCTTCGc                                                                >  1:1189993/1‑100 (MQ=255)                                 gTCCGTGGGGGCCACACGCGAAACTGAAACTCGATCACCTGGGTAAAGAAGTTCTCGAATCCCGTCAGCCGGGTATCCTGGAGCTTTCCCGTACCTTCGc                                                                >  1:1048554/1‑100 (MQ=255)                                  tCCGTGGGGGCCACACGCGAAACTGAAACTCGATCACCTGGGTAAAGAAGTTCTCGAATCCCGTCAGCCGGGTATCCTGGAGCTTTCCCGTACCTTCGCt                                                               >  2:1297838/1‑100 (MQ=255)                                   ccGTGGGGGCCACACGCGAAACTGAAACTCGATCACCTGGGTAAAGAAGTTCTCGAATCCCGTCAGCCGGGTATCCTGGAGCTTTCCCGTACCTTCGCTc                                                              >  2:1956828/1‑100 (MQ=255)                                   ccGTGGGGGCCACACGCGAAACTGAAACTCGATCACCTGGGTAAAGAAGTTCTCGAATCCCGTCAGCCGGGTATCCTGGAGCTTTCCCGTACCTTCGCTc                                                              >  2:620332/1‑100 (MQ=255)                                     gTGGGGGCCACACGCGAAACTGAAACTCGATCACCTGGGTAAAGAAGTTCTCGAATCCCGTCAGCCGGGTATCCTGGAGCTTTCCCGTACCTTCGCTCAc                                                            >  1:273099/1‑100 (MQ=255)                                      tGGGGGCCACACGCGAAACTGAAACTCGATCACCTGGGTAAAGAAGTTCTCGAATCCCGTCAGCCGGGTATCCTGGAGCTTTCCCGTACCTTCGCTCACg                                                           >  1:1262575/1‑100 (MQ=255)                                          ggCCACACGCGAAACTGAAACTCGATCACCTGGGTAAAGAAGTTCTCGAATCCCGTCAGCCGGGTATCCTGGAGCTTTCCCGTACCTTCGCTCACGTCGa                                                       <  2:510369/100‑1 (MQ=255)                                                 cgcgAAACTGCAACTCGATCACCTGGGTAAAGAAGTTCTCGAATCCCGTCAGCCGGGTATCCTGGATCTTTCCCGTACCTTCGCTCACGTCGATCCGGTg                                                <  2:1530721/100‑1 (MQ=255)                                                 cgcgAAACTGAAACTCGATCACCTGGGTAAAGAAGTTCTCGAATCCCGTCAGCCGGGTATCCTGGAGCTTTCCCGTACCTTCGCTCACGTCGATCCGGTg                                                <  2:651937/100‑1 (MQ=255)                                                 cgcgAAACTGAAACTCGATCACCTGGGTAAAGAAGTTCTCGAATCCCGTCAGCCGGGTATCCTGGAGCTTTCCCGTACCTTCGCTCACGTCGATCCGGTg                                                <  2:1312407/100‑1 (MQ=255)                                                    gaaactgaaactCGATCACCTGGGTAAAGAAGTTCTCGAATCCCGTCAGCCGGGTATCCTGGAGCTTTCCCGTACCTTCGCTCACGTCGATCCGGTGaaa                                             <  1:1956835/100‑1 (MQ=255)                                                    gaaactgaaactCGATCACCTGGGTAAAGAAGTTCTCGAATCCCGTCAGCCGGGTATCCTGGAGCTTTCCCGTACCTTCGCTCACGTCGATCCGGTGaaa                                             <  1:620333/100‑1 (MQ=255)                                                    gaaactgaaactCGATCACCTGGGTAAAGAAGTTCTCGAATCCCGTCAGCCGGGTATCCTGGAGCTTTCCCGTACCTTCGCTCACGTCGATCCGGTGaaa                                             <  1:744700/100‑1 (MQ=255)                                                      aactgaaactCGATCACCTGGGTAAAGAAGTTCTCGAATCCCGTCAGCCGGGTATCCTGGAGCTTTCCCGTACCTTCGCTCACGTCGATCCGGTGAAaga                                           >  2:2137200/1‑100 (MQ=255)                                                      aactgaaactCGATCACCTGGGTAAAGAAGTTCTCGAATCCCGTCAGCCGGGTATCCTGGAGCTTTCCCGTACCTTCGCTCACGTCGATCCGGTGAAaga                                           >  2:101248/1‑100 (MQ=255)                                                       actgaaactCGATCACCTGGGTAAAGAAGTTCTCGAATCCCGTCAGCCGGGTATCCTGGAGCTTTCCCGTACCTTCGCTCACGTCGATCCGGTGAAagag                                          >  2:2329609/1‑100 (MQ=255)                                                        cgctcttccGATCTCCTGGGTAAAGAAGTTCTCGAATCCCGTCAGCCGGGTATCCTGGAGCTTTCCCGTACCTTCGCTCACGTCGATCCGGTGAAAGAGc                                         <  2:1473757/92‑1 (MQ=255)                                                         tgaaactCGATCACCTGGGTAAAGAAGTTCTCGAATCCCGTCAGCCGGGTATCCTGGAGCTTTCCCGTACCTTCGCTCACGTCGATCCGGTGAAAGAGcc                                        <  1:2329618/100‑1 (MQ=255)                                                         tgaaactCGATCACCTGGGTAAAGAAGTTCTCGAATCCCGTCAGCCGGGTATCCTGGAGCTTTCCCGTACCTTCGCTCACGTCGATCCGGTGAAAGAGcc                                        >  1:1818226/1‑100 (MQ=255)                                                          gaaactCGATCACCTGGGTAAAGAAGTTCTCGAATCCCGTCTGCCGGGTATCCTGGAGCTTTCCCGTACCTTCGCTCACGTCGATCCGGTGAAAGAGCCg                                       <  1:1454010/100‑1 (MQ=255)                                                          gaaactCGATCACCTGGGTAAAGAAGTTCTCGAATCCCGTCAGCCGGGTATCCTGGAGCTTTCCCGTACCTTCGCTCACGTCGATCCGGTGAAAGAGcc                                        <  1:101250/99‑1 (MQ=255)                                                          gaaactCGATCACCTGGGTAAAGAAGTTCTCGAATCCCGTCAGCCGGGTATCCTGGAGCTTTCCCGTACCTTCGCTCACGTCGATCCGGTGAAAGAGCCg                                       <  1:1296323/100‑1 (MQ=255)                                                          gaaactCGATCACCTGGGTAAAGAAGTTCTCGAATCCCGTCAGCCGGGTATCCTGGAGCTTTCCCGTACCTTCGCTCACGTCGATCCGGTGAAAGAGCCg                                       <  2:2168982/100‑1 (MQ=255)                                                            aactCGATCACCTGGGTAAAGAAGTTCTCGAATCCCGTCAGCCGGGTATCCTGGAGCTTTCCCGTACCTTCGCTCACGTCGATCCGGTGAAAGAGCCGAt                                     >  1:1503184/1‑100 (MQ=255)                                                              ctCGATCACCTGGGTAAAGAAGTTCTCGAATCCCGTCAGCCGGGTATCCTGGAGCTTTCCCGTACCTTCGCTCACGTCGATCCGGTGAAAGAGCCGATTc                                   <  1:669552/100‑1 (MQ=255)                                                                 gATCACCTGGGTAAAGAAGTTCTCGAATCCCGTCAGCCGGGTATCCTGGAGCTTTCCCGTACCTTCGCTCACGTCGATCCGGTGAAAGAGCCGATTCCgg                                <  1:1143556/100‑1 (MQ=255)                                                                 gATCACCTGGGTAAAGAAGTTCTCGAATCCCGTCAGCCGGGTATCCTGGAGCTTTCCCGTACCTTCGCTCACGTCGATCCGGTGAAAGAGCCGATTCCgg                                >  1:2066405/1‑100 (MQ=255)                                                                 gATCACCTGGGTAAAGAAGTTCTCGAATCCCGTCAGCCGGGTATCCTGGAGCTTTCCCGTACCTTCGCTCACGTCGATCCGGTGAAAGAGCCGATTCCgg                                >  1:1277294/1‑100 (MQ=255)                                                                  aTCACCTGGGTAAAGAAGTTCTCGAATCCCGTCAGCCGGGTATCCTGGAGCTTTCCCGTACCTTCGCTCACGTCGATCCGGTGAAAGAGCCGATTCCGGt                               <  2:982052/100‑1 (MQ=255)                                                                  aTCACCTGGGTAAAGAAGTTCTCGAATCCCGTCAGCCGGGTATCCTGGAGCTTTCCCGTACCTTCGCTCACGTCGATCCGGTGAAAGAGCCGATTCCGGt                               >  1:1010114/1‑100 (MQ=255)                                                                      ccTGGGTAAAGAAGTTCTCGAATCCCGTCAGCCGGGTATCCTGGAGCTTTCCCGTACCTTCGCTCACGTCGATCCGGTGAAAGAGCCGATTCCGGTTATc                           >  2:454896/1‑100 (MQ=255)                                                                      ccTGGGTAAAGAAGTTCTCGAATCCCGTCAGCCGGGTATCCTGGAGCTTTCCCGTACCTTCGCTCACGTCGATCCGGTGAAAGAGCAGAtcggaagagca                           >  1:1473761/1‑90 (MQ=255)                                                                       tctggTAAAGAAGTTCTCGAATCCCGTCAGCCGGGTATCCTGGAGCTTTCCCGTACCTTCGCTCACGTCGATCCGGTGAAAGAGCCGATTCCGGTTATcc                          <  2:519629/97‑1 (MQ=255)                                                                       cTGGGTAAAGAAGTTCTCGAATCCCGTCAGCCGGGTATCCTGGAGCTTTCCCGTACCTTCGCTCACGTCGATCCGGTGAAAGAGCCGATTCCGGTTATcc                          <  1:1837875/100‑1 (MQ=255)                                                                        gatctAAAGAAGTTCTCGAATCCCGTCAGCCGGGTATCCTGGAGCTTTCCCGTACCTTCGCTCACGTCGATCCGGTGAAAGAGCCGATTCCGGTTATccc                         <  2:1635237/96‑1 (MQ=255)                                                                          ggTAAAGAAGTTCTCGAATCCCGTCAGCCGGGTATCCTGGAGCTTTCCCGTACCTTCGCTCACGTCGATCCGGTGAAAGAGCCGATTCCGGTTATccaga                       >  1:519631/1‑97 (MQ=255)                                                                            tAAAGAAGTTCTCGAATCCCGTCAGCCGGGTATCCTGGAGCTTTCCCGTACCTTCGCTCACGTCGATCCGGTGAAAGAGCCGATTCCGGTTATCCCAAcc                     >  1:1090340/1‑100 (MQ=255)                                                                            tAAAGAAGTTCTCGAATCCCGTCAGCCGGGTATCCTGGAGCTTTCCCGTACCTTCGCTCACGTCGATCCGGTGAAAGAGCCGATTCCGGTTATCCCAAcc                     >  2:2126567/1‑100 (MQ=255)                                                                            tAAAGAAGTTCTCGAATCCCGTCAGCCGGGTATCCTGGAGCTTTCCCGTACCTTCGCTCACGTCGATCCGGTGAAAGAGCCGATTCCGGTTATCCCAAcc                     <  2:1818221/100‑1 (MQ=255)                                                                            tAAAGAAGTTCTCGAATCCCGTCAGCCGGGTATCCTGGAGCTTTCCCGTACCTTCGCTCACGTCGATCCGGTGAAAGAGCCGATTCCGGTTATCCCAAcc                     <  2:372971/100‑1 (MQ=255)                                                                             aaaGAAGTTCTCGAATCCCGTCAGCCGGGTATCCTGGAGCTTTCCCGTACCTTCGCTCACGTCGATCCGGTGAAAGAGCCGATTCCGGTTATCCCagatc                    >  1:1635242/1‑96 (MQ=255)                                                                              aagaagTTCTCGAATCCCGTCAGCCGGGTATCCTGGAGCTTTCCCGTACCTTCGCTCACGTCGATCCGGTGAAAGAGCCGATTCCGGTTATCCCAACCTg                   >  2:1585067/1‑100 (MQ=255)                                                                              aagaagTTCTCGAATCCCGTCAGCCGGGTATCCTGGAGCTTTCCCGTACCTTCGCTCACGTCGATCCGGTGAAAGAGCCGATTCCGGTTATCCCAACCTg                   >  2:1500311/1‑100 (MQ=255)                                                                               agaagTTCTCGAATCCCGTCAGCCGGGTATCCTGGAGCTTTCCCGTACCTTCGCTCACGTCGATCCGGTGAAAGAGCCGATTCCGGTTATCCCAACCTGt                  <  1:1500315/100‑1 (MQ=255)                                                                               agaagTTCTCGAATCCCGTCAGCCGGGTATCCTGGAGCTTTCCCGTACCTTCGCTCACGTCGATCCGGTGAAAGAGCCGATTCCGGTTATCCCAACCTGt                  >  1:1925682/1‑100 (MQ=255)                                                                               agaagTTCTCGAATCCCGTCAGCCGGGTATCCTGGAGCTTTCCCGTACCTTCGCTCACGTAGATCCGGTGAAAGAGCCGCTTCCGGTTATCCCCACCTGt                  >  1:96662/1‑100 (MQ=255)                                                                                     tctcGAATCCCGTCAGCCGGGTATCCTGGAGCTTTCCCGTACCTTCGCTCACGTCGATCCGGTGAAAGAGCCGATTCCGGTTATCCCAACCTGTCACTAc            >  1:1617859/1‑100 (MQ=255)                                                                                        cGAATCCCGTCAGCCGGGTATCCTGGAGCTTTCCCGTACCTTCGCTCACGTCGATCCGGTGAAAGAGCCGATTCCGGTTATCCCAACCTGTCACTACatg         >  2:1742378/1‑100 (MQ=255)                                                                                        cGAATCCCGTCAGCCGGGTATCCTGGAGATTTCCCGTACCTTCGCTAACGTCGATCCGGTGAAAGAGCCGATTCCGGTTATCCCAACATGTCACTACatg         >  2:2179969/1‑100 (MQ=255)                                                                                         gAATCCCGTCTGCCGGGTATCCTGGAGCTTTCCCGTACCTTCGCTCACGTCGATCCGGTGAAAGAGCCGATTCCGGTTATCCCAACCTGTCACTACatga        >  2:1888580/1‑100 (MQ=255)                                                                                         gAATCCCGTCAGCCGGGTATCCTGGAGCTTTCCCGTACCTTCGCTCACGTCGATCCGGTGAAAGAGCCGATTCCGGTTATCCCAACCTGTCACTACatga        >  2:66520/1‑100 (MQ=255)                                                                                           aTCCCGTCAGCCGGGTATCCTGGAGCTTTCCCGTACCTTCGCTCACGTCGATCCGGTGAAAGAGCCGATTCCGGTTATCCCAACCTGTCACTACatgatg      >  2:1984702/1‑100 (MQ=255)                                                                                            tCCCGTCAGCCGGGTATCCTGGAGCTTTCCCGTACCTTCGCTCACGTCGATCCGGTGAAAGAGCCGATTCCGGTTATCCCAACCTGTCACTACATGATgg     <  2:1116226/100‑1 (MQ=255)                                                                                               cGTCAGCCGGGTATCCTGGAGCTTTCCCGTACCTTCGCTCACGTCGATCCGGTGAAAGAGCCGATTCCGGTTATCCCAACCTGTCACTACATGATGGGCg  <  2:1491323/100‑1 (MQ=255)                                                                                               cGTCAGCCGGGTATCCTGGAGCTTTCCCGTACCTTCGCTCACGTCGATCCGGTGAAAGAGCCGATTCCGGTTATCCCAACCTGTCACTACATGATGGGCg  >  2:2096318/1‑100 (MQ=255)                                                                                                   |                                                                                                GATCGAAATCCGTGAAGGTCGCGGCTGTGATGGTCCGTGGGGGCCACACGCGAAACTGAAACTCGATCACCTGGGTAAAGAAGTTCTCGAATCCCGTCTGCCGGGTATCCTGGAGCTTTCCCGTACCTTCGCTCACGTCGATCCGGTGAAAGAGCCGATTCCGGTTATCCCAACCTGTCACTACATGATGGGCG  >  NC_000913/756788‑756981

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