Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	756,180	G→T	G92W (GGG→TGG)	sdhA →	succinate dehydrogenase, flavoprotein subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	756,180	0	G	T	96.0%	245.4 / ‑0.8	74	G92W (GGG→TGG)	sdhA	succinate dehydrogenase, flavoprotein subunit
Reads supporting (aligned to +/- strand):  ref base G (0/3);  new base T (41/30);  total (41/33)
Fisher's exact test for biased strand distribution p-value = 8.42e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACCCATGAAGATAACTGGGAATGGCATATGTACGACACCGTGAAAGGGTCGGACTATATCGGTGACCAGGACGCGATTGAATATATGTGTAAAACCGGGCCGGAAGCGATTCTGGAACTCGAACACATGGGCCTGCCGTTCTCGCGTCTCGATGATGGTCGTATCTATCAACGTCCGTTTGGCGGTCAG  >  NC_000913/756084‑756272                                                                                                 |                                                                                             aCCCATGAAGATAACTGGGAATGGCATATGTACGACACCGTGAAAGGGTCGGACTATATCGGTGACCAGGACGCGATTGAATATATCTGTAAAACCTGGc                                                                                           <  2:1157251/100‑1 (MQ=255)  cccATGAAGATAACTGGGAATGGCATATGTACGACACCGTGAAAGGGTCGGACTATATCGGTGACCAGGACGCGATTGAATATATGTGTAAAACCTGGcc                                                                                          <  2:1510244/100‑1 (MQ=255)  cccATGAAGATAACTGGGAATGGCATATGTACGACACCGTGAAAGGGTCGGACTATATCGGTGACCAGGACGCGATTGAATATATGTGTAAAACCTGGcc                                                                                          >  2:309687/1‑100 (MQ=255)  cccATGAAGATAACTGGGAATGGCATATGTACGACACCGTGAAAGGGTCGGACTATATCGGGGACCAGGACGCGATTGAATATATGTGTAAAACCTGGcc                                                                                          >  2:320920/1‑100 (MQ=255)   ccATGAAGATAATTGGGAATGGCATATGTATGACACCGTTAAAGGGTCGGACTATATCGGTGACCAGGACGCTATTGAATATATGTGTAATACCTGGCCg                                                                                         <  1:469522/100‑1 (MQ=255)    cATGAAGATAACTGGGAATGGCATATGTACGACACCGTGAAAGGGTCGGACTATATCGGTGACCAGGACGCGATTGAATATATGTGTAAAACCTGGCCgg                                                                                        >  2:519137/1‑100 (MQ=255)    cATGAAGATAACTGGGAATGGCATATGTACGACACCGTGAAAGGGTCGGACTATATCGGTGACCAGGACGCGATTGAATATATGTGTAAAACCTGGCCgg                                                                                        >  2:244558/1‑100 (MQ=255)     aTGAAGATAACTGGGAATGGCATATGTACGACACCGTGAAAGGGTCGGACTATATCGGTGACCAGGACGCGATTGAATATATGTGTAAAACCTGGCCgg                                                                                        <  1:1618360/99‑1 (MQ=255)      tGAAGATAACTGGGAATGGCATATGTACGACACCGTGAAAGGGTCGGACTATATCGGTGACCAGGACGCGATTGAATATATGTGTAAAACCTGGCCGGaa                                                                                      >  2:1287408/1‑100 (MQ=255)       gAAGATAACTGGGAATGGCATATGTACGACACCGTGAAAGGGTCGGACTATATCGGTGACCAGGACGCGATTGAATATATGTGTAAAACCTGGCCGGAAg                                                                                     <  2:1635003/100‑1 (MQ=255)         aGATAACTGGGAATGGCATATGTACGACACCGTGAAAGGGTCGGACTATATCGGTGACCAGGACGCGATTGAATAAATGTGTAAAACATGGCCGGAAGCg                                                                                   >  1:1681588/1‑100 (MQ=255)            tAACTGGGAATGGCATATGTACGACACCGTGAAAGGGTCGGACTATATCGGTGACCAGGACGCGATTGAATATATGTGTAAAACCTGGCCGGAAGCGAtt                                                                                <  2:169616/100‑1 (MQ=255)              aCTGGGAATGGCATATGTACGACACCGTGAAAGGGTCGGACTATATCGGTGACCAGGACGCGATTGAATATATGTGTAAAACCTGGCCGGAAGCGATTCt                                                                              >  1:51257/1‑100 (MQ=255)               cTGGGAATGGCATATGTACGACACCGTGAAAGGGTCGGACTATATCGGTGACCAGGACGCGATTGAATATATGTGTAAAACCTGGCCGGAAGCGATTCTg                                                                             >  2:1665292/1‑100 (MQ=255)                tGGGAATGGCATATGTACGACACCGTGAAAGGGTCGGACTATATCGGTGACCAGGACGCGATTGAATATATGTGTAAAACCGGGCCGGAAGCGATTCt                                                                              <  1:1210446/98‑1 (MQ=255)                 gggAATGGCATATGTACGACACCGTGAAAGGGTCGGACTATATCGGTGACCAGGACGCGATTGAATATATGTGTAAAACCTGGCCGGAAGCGATTCTGGa                                                                           <  1:519136/100‑1 (MQ=255)                 gggAATGGCATATGTACGACACCGTGAAAGGGTCGGACTATATCGGTGACCAGGACGCGATTGAATATATGTGTAAAACCTGGCCGGAAGCGATTCTGGa                                                                           <  1:948319/100‑1 (MQ=255)                  ggAATGGCATATGTACGACACCGTGAAAGGGTCGGACTATATCGGTGACCAGGACGCGATTGAATATATGTGTAAAACCTGGCCGGAAGCGATTCTGGaa                                                                          >  1:162672/1‑100 (MQ=255)                    aaTGGCATATATACGACACCGTGAAAGGGTCGGACTATATCGGTGACCAGGACGCGATTGAATATATGTGTAAAACCTGGCCGGAAGCGATTCTGGAACt                                                                        <  1:1114983/100‑1 (MQ=255)                     aTGGCATATGTACGACACCGTGAAAGGGTCGGACTATATCGGTGACCAGCACGCGATTGAATATATGTGTAAAACCTGGCCGGCCGCGATTCTGGAACTc                                                                       >  1:1041598/1‑100 (MQ=255)                       ggCATATGTACGACACCGTGAAAGGGTCGGACTATATCGTTGACCAGGCCGCGATTGAATATATGTGTAAAACCTGGCGGGACGCGATTCTggca                                                                          >  1:1325206/1‑93 (MQ=255)                       ggCATATGTACGACACCGTGAAAGGGTCGGACTATATCGGTGACCAGGACGCGATTGAATATATGTGTAAAACCTGGCCGGAAGCGATTCTGGAACTCGa                                                                     <  2:51257/100‑1 (MQ=255)                         cATATGTACGACACCTGGAAATGGGCGGACTATATCTTTGACCAGGACTCGATTGAATATATTTGTAAAACCTGGCCGGAAGCGATTCTGGAACTCGAac                                                                   <  2:992733/100‑1 (MQ=255)                           tatGTACGACACCGTGAAAGGGTCGGACTATATCGGTGACCAGGACGCGATTGAATATATGTGTAAAACCTGGCCGGAAGCGATTCTGGAACTCGAacac                                                                 >  1:1680815/1‑100 (MQ=255)                             tGTACGACACCGTGAAAGGGTCGGACTATATCGGTGACCAGGACGCGATTGGATATATGTGTAAAACCTGGCCGGAAGCGATTCTGGAACTCGCac                                                                   >  1:712160/1‑96 (MQ=255)                              gTACGACACCGTGAAAGGGTCGGACTATATCGGTGACCAGGACGCGATTGAATATATGTGTAAAACCTGGCCGGAAGCGATTCTGGAACTCGAACACATg                                                              >  1:906327/1‑100 (MQ=255)                                aCGACACCGTGAAAGGGTCGGACTATATCGGTGACCAGGACGCGATTGAATATATGTGTAAAACCTGGCCGGAAGCGATTCTGGAACTCGAACACATggg                                                            >  1:411045/1‑100 (MQ=255)                                 cGACACCGTGAAAGGGTCGGACTATATCGGTGACCATGTAGCGCTTGAAGATATGCGCAAACCCTGGCCGGAAGCAATTCTGGAACTCGACcaca                                                                >  1:1309528/1‑95 (MQ=255)                                 cGACACCGTGAAAGGGTCGGACTATATCGGTGACCAGGACGCGATTGAATATATGTGTAAAACCTGGCCGGAAGCGATTCTGGAACTCGAACACATGGGc                                                           >  2:292883/1‑100 (MQ=255)                                 cGACACCGTGAAAGGGTCGGACTATATCGGTGACCAGGACGCGATTGAATATATGTGTAAAACCTGGCCGGAAGCGATTCTGGAACTCGAACACATGGGc                                                           >  2:1088817/1‑100 (MQ=255)                                    cacCGTGAAAGGGTCGGACTATATCGGTGACCAGGACGCGATTGAATATATGTGTAAAACCTGGCCGGAAGCGATTCTGGAACTCGAACACATGGGCCTg                                                        >  2:729869/1‑100 (MQ=255)                                    cacCGTGAAAGGGTCGGACTATATCGGTGACCAGGACGCGATTGAATATATGTGTAAAACCTGGCCGGAAGCGATTCTGGAACTCGAACACATGGGCCTg                                                        >  2:1177428/1‑100 (MQ=255)                                        gTGAAAGGGTCGGACTATATCGGTGACCAGGACGCGATTGAATATATGTGTAAAACCTGGCCGGAAGCGATTCTGGAACTCGAACACATGGGCCTGCCGt                                                    >  2:892970/1‑100 (MQ=255)                                        gTGAAAGGGTCGGACTATATCGGTGACCAGGACGCGATTGAATATATGTGTAAAACCTGGCCGGAAGCGATTCTGGAACTCGAACACATGGGCCTGCCGt                                                    >  1:288700/1‑100 (MQ=255)                                        gTGAAAGGGTCGGACTATATCGGTGACCAGGACGCGATTGAATATATGTGTAAAACCTGGCCGGAAGCGATTCTGGAACTCGAACACATGGGCCTGCCGt                                                    <  2:1069348/100‑1 (MQ=255)                                        gTGAAAGGGTCGGACTATATCGGTGACCAGGACGCGATTGAATATATGTGTAAAACCTGGCCGGAAGCGATTCTGGAACTCGAACACATAGGCCTGCCGt                                                    >  1:221106/1‑100 (MQ=255)                                               ggTCGGACTATATCGGTGACCAGGACGCGATTGAATATATGTGTAAAACCTGGCCGGAAGCGATTCTGGAACTCGAACACATGGGCCTGCCGTTCTcgcg                                             >  2:857591/1‑100 (MQ=255)                                                gTCGGACTATATCGGTGACCAGGACGCGATTGAATATATGTGTAAAACCTGGCCGGAAGCGATTCTGGAACTCGAACACATGGGCCTGCCGTTCTCGCGt                                            >  2:1355960/1‑100 (MQ=255)                                                    gACTATATCGGTGACCAGGACGCGATTGAATATATGTGTAAAACCTGGCCGGAAGCGATTCTGGAACTCGAACACATGGGCCTGCCGTTCTCGCGTCTCg                                        <  1:1074689/100‑1 (MQ=255)                                                      cTATATCGGTGACCAGGACGCGATTGAATATATGTGTAAAACCTGGCCGGAAGCGATTCTGGAACTCGAACACATGGGCCTGCCGTTCTCGCGTCTCgat                                      >  1:40583/1‑100 (MQ=255)                                                       tataTCGGTGACCAGGACGCGATTGAATATATGTGTAAAACCTGGCCGGAAGCGATTCTGGAACTCGAACACATGGGCCTGCCGTTCTCGCGTCTCgatg                                     <  1:676651/100‑1 (MQ=255)                                                       tataTCGGTGACCAGGACGCGATTGAATATATGTGTAAAACCTGGCCGGAAGCGATTCTGGAACTCGAACACATGGGCCTGCCGTTCTCGCGTCTCgatg                                     <  2:1681595/100‑1 (MQ=255)                                                        atatCGGTGACCAGGACGCGATTGAATATATGTGTAAAACCTGGCCGGAAGCGATTCTGGAACTCGAACTCATGGGCCTGCCGTTCTCGCGTCTCgatga                                    <  1:553984/100‑1 (MQ=255)                                                           tCGGTGACCAGGACGCGATTGAATATATGTGTAAAACCTGGCCGGAAGCGATTCTGGAACTCGAACACATGGGCCTGCCGTTCTCGCGTCTCGATGATgg                                 <  2:363698/100‑1 (MQ=255)                                                            cGGTGACCAGGACGCGATTGAAAAAAAAAGTAAAACCTGGCCGGAAGCGATTCTGGAACTAGAACACATGGGCCTGCCGTTCTCGCGTCTCGATGATGGt                                >  1:733327/1‑100 (MQ=255)                                                             ggTGACCAGGACGCGATTGAATATATGTGTAAAACCTGGCCGGAAGCGATTCTGGAACTCGAACACATGGGCCTGCCGTTCTCGCGTCTCGATGATGGTc                               >  1:1581456/1‑100 (MQ=255)                                                              gTGACCAGGACGCGATTGAATATATGTGTAAAACCTTGCCGGAAGCGATTCTGGAACTCGAACACATGGGCCTGCCGTTCTCGCGTCTCGATGATGGTCg                              <  2:1193812/100‑1 (MQ=255)                                                                gACCAGGACGCGATTGAATATATGTGTAAAACCTGGCCGGAAGCGATTCTGGAACTCGAACACATGGGCCTGCCGTTCTCGCGTCTCGATGATGGTCGta                            >  1:468119/1‑100 (MQ=255)                                                                  ccAGGACGCGATTGAATATATGTGTAAAACCTGGCCGGAAGCGATTCTGGAACTCGAACACATGGGCCTGCCGTTCTCGCGTCTCGATGATGGTCGtatc                          >  1:1293503/1‑100 (MQ=255)                                                                  ccAGGACGCGATTGAATATATGTGTAAAACCTGGCCGGAAGCGATTCTGGAACTCGAACACATGGGCCTGCCGTTCTCGCGTCTCGATGATGGTCGtatc                          >  1:623319/1‑100 (MQ=255)                                                                   cAGGACGCGATTGAATATATGTGTAAACCCGGGCCGGAAGCGATTCTGGAACTAGAACGCATGGGCCTGCCGTTCGCGCGTCTCGATGATGGTCGtatct                         <  1:169756/100‑1 (MQ=255)                                                                   cAGGACGCGATTGAATATATGTGTAAAACCGGGCCGGAAGCGATTCTGGAACTCGAACACATGGGCCTGCCGTTCTCGCGTCTCGATGATGGTCGtatct                         <  1:913879/100‑1 (MQ=255)                                                                    aGGACGCGATTGAATATATGTGTAAAACCTGGCCGGAAGCGATTCTGGAACTCGAACACATGGGCCTGCCGTTCTCGCGTCTCGATGATGGTCGtatcta                        <  2:843515/100‑1 (MQ=255)                                                                     ggACGCGATTGAATATATGTGTAAAACCTGGCCGGAAGCGATTCTGGAACTCGAACACATGGGCCTGCCGTTCTCGCGTCTCGATGATGGTCGtatctat                       <  2:655545/100‑1 (MQ=255)                                                                        cgtgATTGAATATATGTGTAAGGTCTGGCCGGAAGCGTTTCTGTAACTCGAACACTTGGGCCTGCCGTTCTCGCGTCTCGATGATGGTCGTATCTATCa                     <  1:296132/96‑1 (MQ=255)                                                                        cgcgATTGAATATATGTGTAAAACCTGGCCGGAAGCGATTCTGGAACTCGAACACATGGGCCTGCCGTTCTCGCGTCTCGATGATGGTCGTATCTATCaa                    >  1:61272/1‑100 (MQ=255)                                                                        cgcgATTGAATATATGTGTAAAACCTGGCCGGAAGCGATTCTGGAACTCGAACACATGGGCCTGCCGTTCTCGCGTCTCGATGATGGTCGTATCTATCaa                    <  2:954766/100‑1 (MQ=255)                                                                         gcgATTGAATATATGTGTAAAACCTGGCCGGAAGCGATTCTGGAACTCGAACACATGGGCCTGCCGTTCTCGCGTCTCGATGATGGTCGTATCTATCAAc                   <  2:65623/100‑1 (MQ=255)                                                                         gcgATTGAATATATGTGTAAAACCTGGCCGGAAGCGATTCTGGAACTCGAACACATGGGCCTGCCGTTCTCGCGTCTCGATGATGGTCGTATCTATCAAc                   <  1:548282/100‑1 (MQ=255)                                                                           gATTGAATATATGTGTAAAACCTGGCCGGAAGCGATTCTGGAACTCGAACACATGGGCCTGCCGTTCTCGCGTCTCGATGATGGTCGTATCTATCAACGt                 >  1:794329/1‑100 (MQ=255)                                                                            aTTGAATATATGTGTAAAACCTGGCCGGAAGCGATTCTGGAACTCGAACACATGGGCCTGCCGTTCTCGCGTCTCGATGATGGTCGTATCTATCAACGTc                >  2:558215/1‑100 (MQ=255)                                                                             ttGAATATATGTGTAAAACCTGGCCGGAAGCGATTCTGGAACTCGAACACATGGGCCTGCCGTTCTCGCGTCTCGATGATGGTCGTATCTATCAACGTcc               >  2:671503/1‑100 (MQ=255)                                                                              tGAATATATGTGTAAAACCTGGCCGGAAGCGATTCTGGAACTCGAACACATGGGCCTGCCGTTCTCGCGTCTCGATGATGGTCGTATCTATCAACGTCCg              <  2:162672/100‑1 (MQ=255)                                                                              tGAATATATGTGGAATACCTGGCCTGAAGCGATTCTGGACTTCTAACATATGGGCCTGCCGTTCTCGCGTATCGATGATGGTCGTATCTATCAACGTCCg              <  1:320919/100‑1 (MQ=255)                                                                                aaTATATGTGTAAAACCTGGCCGGAAGCGATTCTGGAACTCGAACACATGGGCCTGCCGTTCTCGCGTCTCGATGATGGTCGTATCTATCAACGTCCGtt            <  2:712162/100‑1 (MQ=255)                                                                                 atatatGTGTAAAACCTGGCCGGAAGCGATTCTGGAACTCGAACACATGGGCCTGCCGTTCTCGCGTCTCGATGATGGTCGTATCTATCAACGTCCGttt           >  2:963085/1‑100 (MQ=255)                                                                                   atatGTGTAAAACCTGGCCGGAAGCGATTCTGGAACTCGAACACATGGGCCTGCCGTTCTCGCGTCTCGATGATGGTCGTATCTATCAACGTCCGTTTgg         >  2:119789/1‑100 (MQ=255)                                                                                    tatGTGTAAAACCTGGCCGGAAGCGATTCTGGAACTCGAACACATGGGCCTGCCGTTCTCGCGTCTCGATGATGGTCGTATCTATCAACGTCCGTTTGGc        >  2:828927/1‑100 (MQ=255)                                                                                    tatGTGTAAAACCTGGCCGGAAGCGATTCTGGAACTCGAACACATGGGCCTGCCGTTCTCGCGTCTCGATGATGGTCGTATCTATCAACGTCCGTTTGGc        >  2:883712/1‑100 (MQ=255)                                                                                     atGTGTAAAACCTGGCCGGAAGCGATTCTGGAACTCGAACACATGGGCCTGCCGTTCTCGCGTCTCGATGATGGTCGTATCTATCAACGTCCGTTTGGCg       <  1:143625/100‑1 (MQ=255)                                                                                      tgtgTAAAACCTGGCCGGAAGCGATTCTGGAACTCGAACACATGGGCCTGCCGTTCTCGCGTCTCGATGATGGTCGTATCTATCAACGTCCGTTTGGCgg      <  2:592225/100‑1 (MQ=255)                                                                                       gtgtAAAACCTGGCCGGAAGCGATTCTGGAACTCGAACACATGGGCCTGCCGTTCTCGCGTCTCGATGATGGTCGTATCTATCAACGTCCGTTTGGCGGt     >  2:375973/1‑100 (MQ=255)                                                                                        tgtAAAACCTGGCCGGAAGCGATTCTGGAACTCGAACACATGGGCCTGCCGTTCTCGCGTCTCGATGATGGTCGTATCTATCAACGTCCGTTTGGCGGTc    <  2:1549517/100‑1 (MQ=255)                                                                                          tAAAACCTGGCCGGAAGCGATTCTGGAACTCGAACACATGGGCCTGCCGTTCTCGCGTCTCGATGATGGTCGTATCTATCAACGTCCGTTTGGCGGTCAg  >  1:425778/1‑100 (MQ=255)                                                                                                 |                                                                                             ACCCATGAAGATAACTGGGAATGGCATATGTACGACACCGTGAAAGGGTCGGACTATATCGGTGACCAGGACGCGATTGAATATATGTGTAAAACCGGGCCGGAAGCGATTCTGGAACTCGAACACATGGGCCTGCCGTTCTCGCGTCTCGATGATGGTCGTATCTATCAACGTCCGTTTGGCGGTCAG  >  NC_000913/756084‑756272

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