Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,669,920	A→G	E110G (GAA→GGA)	treF →	cytoplasmic trehalase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,669,920	0	A	G	97.4%	273.4 / ‑2.8	79	E110G (GAA→GGA)	treF	cytoplasmic trehalase
Reads supporting (aligned to +/- strand):  ref base A (0/2);  new base G (34/42);  total (34/45)
Fisher's exact test for biased strand distribution p-value = 5.02e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.94e-01

TCCCGACTGCGCACCTAAAATGGACCCGCTGGATATCTTAATCCGCTACCGTAAAGTGCGCCGTCATCGTGATTTTGACTTGCGCAAGTTTGTTGAAAACCACTTCTGGCTGCCGGAGGTCTACTCCAGCGAGTATGTATCGGACCCGCAAAATTCCCTGAAAGAGCATATCGACCAGCTGTGGCCGGTGCT  >  NC_000913/3669825‑3670016                                                                                                |                                                                                                 tCCCGACTGCGCACCTAAAATGGACCCGCTGGATATCTTAATCCGCTACCGTAAAGTGCGCCGTCATCGTGATTTTGACTTGCGCAAGTTTGTTGGAAAc                                                                                              <  1:796544/100‑1 (MQ=255) tCCCGACTGCGCACCTAAAATGGACCCGCTGGATATCTTAATCCGCTACCGTAAAGTGCGCCGTCATCGTGATTTTGACTTGCGCAAGTTTGTTGGAAAc                                                                                              <  2:368929/100‑1 (MQ=255)  cccGACTGCGCACCTAAAATGGACCCGCTGGATATCTTAATCCGCTACCGTAAAGTGCGCCGTCATCGTGATTTTGACTTGCGCAAGTTTGTTGGAAAcc                                                                                             >  1:1666244/1‑100 (MQ=255)   ccGACTGCGCACCTAAAATGGACCCGCTGGATATCTTAATCCGCTACCGTAAAGTGCGCCGTCATCGTGATTTTGACTTGCGCAAGTTTGTTGGAAACCa                                                                                            >  1:1582742/1‑100 (MQ=255)          cgcACCTAGAATGGACCCGCTGGATATCTTAATCCGCTACCGTAAAGTGCGCCGTCATCGTGATTTTGACTTGCGCAAGTTTGTTGAAAACCACTTCTgg                                                                                     <  1:122311/100‑1 (MQ=255)          cgcACCTAAAATGGACCCGCTGGATATCTTAATCCGCTACCGTAAAGTGCGCCGTCATCGTGATTTTGACTTGCGCAAGTTTGTTGGAAACCACTTCTgg                                                                                     >  2:1440600/1‑100 (MQ=255)          cgcACCTAAAATGGACCCGCTGGATATCTTAATCCGCTACCGTAAAGTGCGCCGTCATCGTGATTTTGACTTGCGCAAGTTTGTTGGAAACCACTTCTgg                                                                                     >  1:167087/1‑100 (MQ=255)          cgcACCTAAAATGGACCCGCTGGATATCTTAATCCGCTACCGTAAAGTGCGCCGTCATCGTGATTTTGACTTGCGCAAGTTTGTTGGAAACCACTTCTgg                                                                                     <  2:493120/100‑1 (MQ=255)           gcACCTAAAATGGACCCGCTGGATATCTTAATCCGCTACCGTAAAGTGCGCCGTCATCGTGATTTTGACTTGCGCAAGTTTGTTGGAAACCACTTCTGGc                                                                                    <  2:258903/100‑1 (MQ=255)           gcACCTAAAATGGACCCGCTGGATATCTTAATCCGCTACCGTAAAGTGCGCCGTCATCGTGATTTTGACTTGCGCAAGTTTGTTGGAAACCACTTCTGGc                                                                                    >  1:1451112/1‑100 (MQ=255)            cACCTAAAATGGACCCGCTGGATATCTTAATCCGCTACCGTAAAGTGCGCCGTCATCGTGATTTTGACTTGCGCAAGTTTGTTGGAAACCACTTCTGGCt                                                                                   >  1:385251/1‑100 (MQ=255)            cACCTAAAATGGACCCGCTGGATATCTTAATCCGCTACCGTAAAGTGCGCCGTCATCGTGATTTTGACTTGCGCAAGTTTGTTGGAAACCACTTCTGGCt                                                                                   >  1:1835547/1‑100 (MQ=255)            cACCTAAAATGGACCCGCTGGATATCTTAATCCGCTACCGTAAAGTGCGCCGTCATCGTGATTTTGACTTGCGCAAGTTTGTTGGAAACCACTTCTGGCt                                                                                   >  1:1504449/1‑100 (MQ=255)             aCCTAAAATGGACCCGCTGGATATCTTAATCCGCTACCGTAAAGTGCGCCGTCATCGTGATTTTGACTTGCGCAAGTTTGTTGGAAACCACTTCTGGCTg                                                                                  >  2:437602/1‑100 (MQ=255)              ccTAAAATGGACCCGCTGGATATCTTAATCCGCTACCGTAAAGTGCGCCGTCATCGTGATTTTGACTTGCGCACGTTTGTTGGAAACCACTTCTGGCTGc                                                                                 >  2:1581932/1‑100 (MQ=255)                  aaaTGGACCCGCTGGATATCTTAATCCGCTACCGTAAAGTGCGCCGTCATCGTGATTTTGACTTGCGCAAGTTTGTTGGAAACCACTTCTGGCTGCCgg                                                                              >  1:1076187/1‑99 (MQ=255)                   aaTGGACCCGCTGGATATCTTAATCCGCTACCGTAAAGTGCGCCGTCATCGTGATTTTGACTTGCGCAAGTTTGTTGGAAACCACTTCTGGCTGCCGGAg                                                                            >  2:1821936/1‑100 (MQ=255)                   aaTGGACCCGCTGGATATCTTAATCCGCTACCGTAAAGTGCGCCGTCATCGTGATTTTGACTTGCGCAAGTTTGTTGGAAACCACTTCTGGCTGCCGGAg                                                                            >  2:651149/1‑100 (MQ=255)                            gCTGGATATCTTAATCCGCTACCGTAAAGTGCGCCGTCATCGTGATTTTGACTTGCGCAAGTTTGTTGGAAACCACTTCTGGCTGCCGGAGGTCTACTcc                                                                   <  1:1450089/100‑1 (MQ=255)                            gCTGGATATCTTAATCCGCTACCGTAAAGTGCGCCGTCATCGTGATTTTGACTTGCGCAAGTTTGTTGGAAACCACTTCTGGCTGCCGGAGGTCTACTcc                                                                   <  2:1858333/100‑1 (MQ=255)                            gCTGGATATCTTAATCCGCTACCGTAAAGTGCGCCGTCATCGTGATTTTGACTTGCGCAAGTTTGTTGGAAACCACTTCTGGCTGCCGGAGGTCTACTcc                                                                   <  1:363776/100‑1 (MQ=255)                              ttgATATCTTAATCCGCTACCGTAAAGTGCGCCGTCATCGTGATTTTTACTTGCTCAAGTTTGTTGGAAACCACTTCTGGCTGCCGGAGGTCTACTCCAg                                                                 <  1:1121772/98‑1 (MQ=255)                              tGGATATCTTAATCCGCTACCGTAAAGTGGGCCGTCATCGTGATTTTGACTTGCGCAAGTTTGTTGTAAACCACTTCTGGCTGCCGGAGGTCTACTCCAg                                                                 <  1:1168795/100‑1 (MQ=255)                                gATATCTTAATCCGCTACCGTAAAGTGCGCCGTCATCGTGATTTTGACTTGCGCAAGTTTGTTGGAAACCACTTCTGGCTGCCGGAGGTCTACTCCAGCg                                                               <  2:601088/100‑1 (MQ=255)                                atcttCTTAATCCGCTACCGTAAAGTGCGCCGTCATCGTGATTTTGACTTGCGCAAGTTTGTTGGAAACCACTTCTGGCTGCCGGAGGTCTACTCCAGCg                                                               <  2:264843/96‑1 (MQ=255)                                 atatCTTAATCCGCTACCGTAAAGTGCGCCGTCATCGTGATTTTGACTTGCGCAAGTTTGTTGGAAACCACTTCTGGCTGCCGGAGGTCTACTCCAGCGa                                                              <  2:1856873/100‑1 (MQ=255)                                    tCTTAATCCGCTACCGTAAAGTGCGCCGTCATCGTGATTTTGACTTGCGCAAGTTTGTTGGAAACCACTTCTGGCTGCCGGAGGTCTACTCCAGCGAgtn                                                           <  2:575864/100‑2 (MQ=255)                                    tCTTAATCCGCTACCGTAAAGTGCGCCGTCATCGTGATTTTGACTTGCGCAAGTTTGTTGGAAACCACTTCTGGCTGCCGGAGGTCTACTCCAGCGAgta                                                           >  2:396178/1‑100 (MQ=255)                                    tCTTAATCCGCTACCGTAAAGTGCGCCGTCATCGTGATTTTGACTTGCGCAAGTTTGTTGGAAACCACTTCTGGCTGCCGGAGGTCTACTCCAGCGAgta                                                           <  2:878179/100‑1 (MQ=255)                                    tCTTAATCCGCTACCGTAAAGTGCGCCGTCATCGTGATTTTGACTTGCGCAAGTTTGTTGGAAACCACTTCTGGCTGCCGGAGGTCTACTCCAGCGAgat                                                           >  1:264841/1‑98 (MQ=255)                                     cTTAATCCGCTACCGTAAAGTGCGCCGTCATCGTGATTTTGACTTGCGCAAGTTTGTTGGAAACCACTTCTGGCTGCCGGAGGTCTACTCCAGCGAgtat                                                          >  1:802682/1‑100 (MQ=255)                                     cTTAATCCGCTACCGTAAAGTGCGCCGTCATCGTGATTTTGACTTGCGCAAGTTTGTTGGAAACCACTTCTGGCTGCCGGAGGTCTACTCCAGCGAgtat                                                          >  1:1927049/1‑100 (MQ=255)                                      ttAATCCGCTACCGTAAAGTGCGCCGTCATCGTGATTTTGACTTGCGCAAGTTTGTTGGAAACCACTTCTGGCTGCCGGAGGTCTACTCCAGCGAgtatg                                                         <  1:1559470/100‑1 (MQ=255)                                      ttAATCCGCTACCGTAAAGTGCGCCGTCATCGTGATTTTGACTTGCGCAAGTTTGTTGGAAACCACTTCTGGCTGCCGGAGGTCTACTCCAGCGAgtatg                                                         <  2:1239202/100‑1 (MQ=255)                                      ttAATCCGCTACCGTAAAGTGCGCCGTCATCGTGATTTTGACTTGCGCAAGTTTGTTGGAAACCACTTCTGGCTGCCGGAGGTCTACTCCAGCGAgtatg                                                         <  2:808804/100‑1 (MQ=255)                                        aaTCCGCTACCGTAAAGTGCGCCGTCATCGTGATTTTGACTTGCGCAAGTTTGTTGGAAACCACTTCTGGCTGCCGGAGGTCTACTCCAGCGAgtatgta                                                       <  2:366746/100‑1 (MQ=255)                                        aaTCCGCTACCGTAAAGTGCGCCGTCATCGTGATTTTGACTTGCGCAAGTTTGTTGGAAACCACTTCTGGCTGCCGGAGGTCTACTCCAGCGAgtatgta                                                       <  1:1812239/100‑1 (MQ=255)                                          tCCGCTACCGTAAAGTGCGCCGTCATCGTGATTTTGACTTGCGCAAGTTTGTTGGAAACCACTTCTGGCTGCCGGAGGTCTACTCCAGCGAGTATGTATc                                                     <  1:135165/100‑1 (MQ=255)                                           gcGCTACCGTAAAGTGCGCCGTCATCGTTATTTTGACTTGCCCAAGTTTGTTGGAAACCACTTCTGGCTGCCCTAGGTGTACTCCAGCAAGTATGTATc                                                     <  1:1184937/98‑1 (MQ=255)                                           ccGCTACCGTAAAGTGCGCCGTCATCGTGATTTTGACTTGCGCAAGTTTGTTGGAAACCACTTCTGGCTGCCGGAGGTCTACTCCAGCGAGTATGTATCg                                                    >  1:20816/1‑100 (MQ=255)                                            cGCTACCGTAAAGTGCGCCGTCATCGTGATTTTGACTTGCGCAAGTTTGTTGGAAACCACTTCTGGCTGCCGGAGGTCTACTCCAGCGAGTATGTATCgg                                                   >  1:143421/1‑100 (MQ=255)                                              cTACCGTAAAGTGCGCCGTCATCGTGATTTTGACTTGCGCAAGTTTGTTGGAAACCACTTCTGGCTGCCGGAGGTCTACTCCAGCGAGTATGTATCGGAc                                                 <  1:493688/100‑1 (MQ=255)                                               tACCGTAAAGTGCGCCGTCATCGTGATTTTGACTTGCGCAAGTTTGTTGGAAACCACTTCTGGCTGCCGGAGGTCTACTCCAGCGAGTATGTATCGGAcc                                                >  1:665868/1‑100 (MQ=255)                                                 ccGTAAAGTGCGCCGTCATCGTGATTTTGACTTGCGCAAGTTTGTTGGAAACCACTTCTGGCTGCCGGAGGTCTACTCCAGCGAGTATGTATCGGACCCg                                              <  2:1131917/100‑1 (MQ=255)                                                       aGTGCGCCGTCATCGTGATTTTGACTTGCGCAAGTTTGTTGAAAACCACTTCTGGCTGCCGGAGGTCTACTCCAGCGAGTATGTATCGGACCCGCAAAAt                                        <  2:353896/100‑1 (MQ=255)                                                          gcgcCGTCATCGTGATTTTGACTTGCGCAAGTCTGTTGGAAACCACTTCTGGCTGCCGGAGGTCTACTCCAGCGAGTATGTATCGGACCCGCAAAATTcc                                     <  1:1634920/100‑1 (MQ=255)                                                           cgcCGTCATCGTGATTTTGACTTGCGCAAGTTTGTTGGAAACCACTTCTGGCTGCCGGAGGTCTACTCCAGCGAGTATGTATCGGACCCGCAAAATTccc                                    <  2:12000/100‑1 (MQ=255)                                                           cgcCGTCATCGTGATTTTGACTTGCGCAAGTTTGTTGGAAACCACTTCTGGCTGCCGGAGGTCTACTCCAGCGAGTATGTATCGGACCCGCAAAATTccc                                    <  1:1888436/100‑1 (MQ=255)                                                            gcCGTCATCGTGATTTTGACTTGCGCAAGTTTGTTGGAAACCACTTCTGGCTGCCGGAGGTCTACTCCAGCGAGTATGTATCGGACCCGCAAAATTCCCt                                   >  2:92646/1‑100 (MQ=255)                                                            gcCGTCATCGTGATTTTGACTTGCGCAAGTTTGTTGGAAACCACTTCTGGCTGCCGGAGGTCTACTCCAGCGAGTATGTATCGGACCCGCAAAATTCCCt                                   <  2:1510552/100‑1 (MQ=255)                                                            gcCGTCATCGTGATTTTGACTTGCGCAAGTTTGTTGGAAACCACTTCTGGCTGCCGGAGGTCTACTCCAGCGAGTATGTATCGGACCCGCAAAATTCCCt                                   >  1:940495/1‑100 (MQ=255)                                                            gcCGTCATCGTGATTTTGACTTGCGCAAGTTTGTTGGAAACCACTTCTGGCTGCCGGAGGTCTACTCCAGCGAGTATGTATCGGACCCGCAAAATTCCCt                                   <  2:1825071/100‑1 (MQ=255)                                                                  aTCTTGATTTTGACTTGCGCAAGTTTGTTGGAAACCACTTCTGGCTGCCGGAGGTCTACTCCAGCGAGGATGTATCGGACCCGCAAAATTCCCTGAAag                              <  1:1544175/99‑1 (MQ=255)                                                                    cGTGATTTTGACTTGCGCAAGTTTGTTGGAAACCACTTCTGGCTGCCGGAGGTCTACTCCAGCGAGTATGTATCGGACCCGCAAAATTCCCTGAAAGAGc                           <  1:176982/100‑1 (MQ=255)                                                                    cGTGATTTTGACTTGCGCAAGTTTGTTGGAAACCACTTCTGGCTGCCGGAGGTCTACTCCAGCGAGTATGTATCGGACCCGCAAAATTCCCTGAAAGAGc                           <  2:913550/100‑1 (MQ=255)                                                                    cGTGATTTTGACTTGCGCAAGTTTGTTGGAAACCACTTCTGGCTGCCGGAGGTCTACTCCAGCGAGTATGTATCGGACCCGCAAAATTCCCTGAAAGAGc                           <  2:347475/100‑1 (MQ=255)                                                                      tGATTTTGACTTGCTCAAGTTTGTTAGAAACTACTTCTGGCTGCCGGAGGTCTACTCCAGCGAGTATGTATCGTACCCGCAAAATTCCCTGAAAGAGc                           <  1:278508/98‑1 (MQ=255)                                                                        aTTTTGACTTGCGCAAGTTTGTTGGAAACCACTTCTGGCTGCCGGAGGTCTACTCCAGCGAGTATGTATCGGACCCGCAAAATTCCCTGAAAGAGCatat                       <  2:1451120/100‑1 (MQ=255)                                                                         ttcttaCTTGCGCAAGTTTGTTGGAAACCACTTCTGTCTGCCGGAGGTCTACTCCAGCGAGTATGTATCGGACCTGCAAAATTCCCTGAAAGAGCATATc                      <  2:1847944/95‑1 (MQ=255)                                                                          tttGACTTGCGCAAGTTTGTTGGAAACCACTTCTGGCTGCCGGAGGTCTACTCCAGCGAGTATGTATCGGACCCGCAAAATTCCCTGAAAGAGCATATCg                     <  2:1504457/100‑1 (MQ=255)                                                                             gACTTGCGCAAGTTTGTTGGAAACCACTTCTGGCTGCCGGAGGTCTACTCCAGCGAGTATGTATCGGACCCGCAAAATTCCCTGAAAGAGCATATCGAcc                  <  2:954145/100‑1 (MQ=255)                                                                             gACTTGCGCAAGTTTGTTGGAAACCACTTCTGGCTGCCGGAGGTCTACTCCAGCGAGTATGTATCGGACCCGCAAAATTCCCTGAAAGAGCATATCGAcc                  <  2:1625021/100‑1 (MQ=255)                                                                              aCTTGCGCAAGTTTGTTGGAAACCACTTCTGGCTGCCGGAGGTCTACTCCAGCGAGTATGTATCGGACCCGCAAAATTCCCTGAAAGAGCATATCGACCa                 >  1:1009130/1‑100 (MQ=255)                                                                               cTTGCGCAAGTTTGTTGGAAACCACTTCTGGCTGCCGGAGGTCTACTCCAGCGAGTATGTATCGGACCCGCAAAATTCCCTGAAAGAGCATATCGACCAg                >  2:1790716/1‑100 (MQ=255)                                                                               cTTGCGCAAGTTTGTTGGAAACCACTTCTGGCTGCCGGAGGTCTACTCCAGCGAGTATGTATCGGACCCGCAAAATTCCCTGAAAGAGCATATCGACCAg                >  2:44341/1‑100 (MQ=255)                                                                               cTTGCGCAAGTTTGTTGGAAACCACTTCTGGCTGCCGGAGGTCTACTCCAGCGAGTATGTATCGGACCCGCAAAATTCCCTGAAAGAGCATATCGACCAg                >  1:1077357/1‑100 (MQ=255)                                                                                 tGCGCAAGTTTGTTGGAAACCACTTCTGGCTGCCGGAGGTCTACTCCAGCGAGTATGTATCGGACCCGCAAAATTCCCTGAAAGAGCATATCGACCAGCt              >  2:1399243/1‑100 (MQ=255)                                                                                  gcgcAAGTTTGTTGGAAACCACTTCTGGCTGCCGGAGGTCTACTCCAGCGAGTATGTATCGGACCCGCAAAATTCCCTGAAAGAGCATATCGACCAGCtg             <  1:1758082/100‑1 (MQ=255)                                                                                    gcAAGTTTGTTGGAAACCACTTCTGGCTGCCGGAGGTCTACTCCAGCGAGTATGTATCGGACCCGCAAAATTCCCTGAAAGAGCATATCGACCAGCtgtg           <  1:1790705/100‑1 (MQ=255)                                                                                    gcAAGTTTGTTGGAAACCACTTCTGGCTGCCGGAGGTCTACTCCAGCGAGTATGTATCGGACCCGCAAAATTCCCTGAAAGAGCATATCGACCAGCtgtg           >  2:1619371/1‑100 (MQ=255)                                                                                     cAAGTTTGTTGGAAACCACTTCTGGCTGCCGGAGGTCTACTCCAGCGAGTATGTATCGGACCCGCAAAATTCCCTGAAAGAGCATATCGACCAGCTGTgg          >  1:143687/1‑100 (MQ=255)                                                                                      aaGTTTGTTGGAAACCACTTCTGGCTGCCGGAGGTCTACTCCAGCGAGTATGTATCGGACCCGCAAAATTCCCTGAAAGAGCATATCGACCAGCTGTGGc         <  1:701923/100‑1 (MQ=255)                                                                                      aaGTTTGTTGGAAACCACTTCTGGCTGCCGGAGGTCTACTCCAGCGAGTATGTATCGGACCCGCAAAATTCCCTGAAAGAGCATATCGACCAGCTGTGGc         <  1:1315936/100‑1 (MQ=255)                                                                                        gTTTGTTGGAAACCACTTCTGGCTGCCGGAGGTCTACTCCAGCGAGTATGTATCGGACCCGCAAAATTCCCTGAAAGAGCATATCGACCAGCTGTGGcc        <  1:1581921/99‑1 (MQ=255)                                                                                        gTTTGTTGGAAACCACTTCTGGCTGCCGGAGGTCTACTCCAGCGAGTATGTATCGGACCCGCAAAATTCCCTGAAAGAGCATATCGACCAGCTGTGGCCg       >  2:50454/1‑100 (MQ=255)                                                                                         tttGTTGGAAACCACTTCTGGCTGCCGGAGGTCTACTCCAGCGAGTATGTATCGGACCCGCAAAATTCCCTGAAAGAGCATATCGACCAGCTGTGGCCgg      >  2:1945937/1‑100 (MQ=255)                                                                                          ttgttgGAAACCACTTCTGGCTGCCGGAGGTCTACTCCAGCGAGTATGTATCGGACCCGCAAAATTCCCTGAAAGAGCATATCGACCAGCTGTGGCCGGt     >  2:33413/1‑100 (MQ=255)                                                                                           tgttgGAAACCACTTCTGGCTGCCGGAGGTCTACTCCAGCGAGTATGTATCGGACCCGCAAAATTCCCTGAAAGAGCATATCGACCAGCTGTGGCCGGTg    >  1:273859/1‑100 (MQ=255)                                                                                             ttgGAAACCACTTCTGGCTGCCGGAGGTCTACTCCAGCGAGTATGTATCGGACCCGCAAAATTCCCTGAAAGAGCATATCGACCAGCTGTGGCCGGTGCt  >  1:721139/1‑100 (MQ=255)                                                                                                |                                                                                                 TCCCGACTGCGCACCTAAAATGGACCCGCTGGATATCTTAATCCGCTACCGTAAAGTGCGCCGTCATCGTGATTTTGACTTGCGCAAGTTTGTTGAAAACCACTTCTGGCTGCCGGAGGTCTACTCCAGCGAGTATGTATCGGACCCGCAAAATTCCCTGAAAGAGCATATCGACCAGCTGTGGCCGGTGCT  >  NC_000913/3669825‑3670016

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