Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	16,941	A→G	V7A (GTG→GCG)	gef ←	Gef

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	16,941	0	A	G	97.4%	232.5 / ‑2.5	77	V7A (GTG→GCG)	gef	Gef
Reads supporting (aligned to +/- strand):  ref base A (1/1);  new base G (10/65);  total (11/66)
Fisher's exact test for biased strand distribution p-value = 2.67e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.20e-01

CGGTGATACAGATGACGATCAGGGCGACAATCATCGCCTTATGCTGCTTCATTGCTCTCTTCTCCTTGACCTTACGGTCAGTAAGAGGCACTCTACATGTGTTCAGCATATAGGAGGCCTCGGGTTGATGGTAAAATATCACTCGGGGCTTTTCTCTATCTGCCGTTCAGCTAATGCC  >  NC_000913/16852‑17029                                                                                          |                                                                                         cggTGATACAGATGACGATCAGGGCGACAATCATCGCCTTATGCTGCTTCATTGCTCTCTTCTCCTTGACCTTACGGTCAGTAAGAGGCGctc                                                                                       <  2:4005209/93‑1 (MQ=255)  ggTGATACAGATGACGATCAGGGCGACAATCATCGCCTTATGCTGCTTCATTGCTCTCTTCTCCTTGACCTTACGGTCAGTAAGAGGCGctct                                                                                      <  1:2479100/93‑1 (MQ=255)   gTGATACAGATGACGATCAGGGCGACAATCATCGCCTTATGCTGCTTCATTGCTCTCTTCTCCTTGACCTTACGGTCAGTAAGAGGCGCTCTa                                                                                     <  1:4285925/93‑1 (MQ=255)    tGATACAGATGACGATCAGGGCGACAATCATCGCCTTATGCTGCTTCATTGCTCTCTTCTCCTTGACCTTACGGTCAGTAAGAGGCGCTCTAc                                                                                    <  1:1814079/93‑1 (MQ=255)     gATACAGATGACGATCAGGGCGACAATCATCGCCTTATGCTGCTTCATTGCTCTCTTCTCCTTGACCTTACGGTCAGTAAGAGGCGCTCTACa                                                                                   <  1:1359199/93‑1 (MQ=255)     gATACAGATGACGATCAGGGCGACAATCATCGCCTTATGCTGCTTCATTGCTCTCTTCTCCTTGACCTTACGGTCAGTAAGAGGCGCTCTACa                                                                                   <  2:938096/93‑1 (MQ=255)     gATACAGATGACGATCAGGGCGACAATCATCGCCTTATGCTGCTTCATTGCTCTCTTCTCCTTGACCTTACGGTCAGTAAGAGGCGCTCTACa                                                                                   <  2:2543326/93‑1 (MQ=255)      aTACAGATGACGATCAGGGCGACAATCATCGCCTTATGCTGCTTCATTGCTCTCTTCTCCTTGACCTTACGGTCAGTAAGAGGCGCTCTACAt                                                                                  <  2:133821/93‑1 (MQ=255)       tACAGATGACGATCAGGGCGACAATCATCGCCTTATGCTGCTTCATTGCTCTCTTCTCCTTGACCTTACGGTCAGTAAGAGGCGCTCTACAtg                                                                                 <  1:1750068/93‑1 (MQ=255)         cAGATGACGATCAGGGCGACAATCATCGCCTTATGCTGCTTCATTGCTCTCTTCTCCTTGACCTTACGGTCAGTAAGAGGCGCTCTACAtgtg                                                                               <  1:2131579/93‑1 (MQ=255)           gATGACGATCAGGGCGACAATCATCGCCTTATGCTGCTTCATTGCTCTCTTCTCCTTGACCTTACGGTCAGTAAGAGGCGCTCTACATGTGtt                                                                             >  1:1366076/1‑93 (MQ=255)           gATGACGATCAGGGCGACAATCATCGCCTTATGCTGCTTCATTGCTCTCTTCTCCTTGACCTTACGGTCAGTAAGAGGCGCTCTACATGTGtt                                                                             >  2:1934836/1‑93 (MQ=255)           gATGACGATCAGGGCGACAATCATCGCCTTATGCTGCTTCATTGCTCTCTTCTCCTTGACCTTACGGTCAGTAAGAGGCGCTCTACATGTGtt                                                                             >  2:1934837/1‑93 (MQ=255)            aTGACGATCAGGGCGACAATCATCGCCTTATGCTGCTTCATTGCTCTCTTCTCCTTGACCTTACGGTCAGTAAGAGGCGCTCTACATGTGTTc                                                                            <  1:1429069/93‑1 (MQ=255)             tGACGATCAGGGCGACAATCATCGCCTTATGCTGCTTCATTGCTCTCTTCTCCTTGACCTTACGGTCAGTAAGAGGCGCTCTACATGTGTTCa                                                                           <  2:2959533/93‑1 (MQ=255)             tGACGATCAGGGCGACAATCATCGCCTTATGCTGCTTCATTGCTCTCTTCTCCTTGACCTTACGGTCAGTAAGAGGCGCTCTACATGTGTTCa                                                                           <  1:1082807/93‑1 (MQ=255)             tGACGATCAGGGCGACAATCATCGCCTTATGCTGCTTCATTGCTCTCTTCTCCTTGACCTTACGGTCAGTAAGAGGCGCTCTACATGTGTTCa                                                                           <  2:3371613/93‑1 (MQ=255)             tGACGATCAGGGCGACAATCATCGCCTTATGCTGCTTCATTGCTCTCTTCTCCTTGACCTTACGGTCAGTAAGAGGCGCTCTACATGTGTTCa                                                                           >  2:168521/1‑93 (MQ=255)              gACGATCAGGGCGACAATCATCGCCTTATGCTGCTTCATTGCTCTCTTCTCCTTGACCTTACGGTCAGTAAGAGGCGCTCTACATGTGTTCAg                                                                          <  2:2932444/93‑1 (MQ=255)              gACGATCAGGGCGACAATCATCGCCTTATGCTGCTTCATTGCTCTCTTCTCCTTGACCTTACGGTCAGTAAGAGGCGCTCTACATGTGTTCAg                                                                          <  2:1212171/93‑1 (MQ=255)                cGATCAGGGCGACAATCATCGCCTTATGCTGCTTCATTGCTCTCTTCTCCTTGACCTTACGGTCAGTAAGAGGCGCTCTACATGTGTTCAGCa                                                                        <  2:3662785/93‑1 (MQ=255)                cGATCAGGGCGACAATCATCGCCTTATGCTGCTTCATTGCTCTCTTCTCCTTGACCTTACGGTCAGTAAGAGGCGCTCTACATGTGTTCAGCa                                                                        <  2:3970615/93‑1 (MQ=255)                cGATCAGGGCGACAATCATCGCCTTATGCTGCTTCATTGCTCTCTTCTCCTTGACCTTACGGTCAGTAAGAGGCGCTCTACATGTGTTCAGCa                                                                        <  2:3970613/93‑1 (MQ=255)                cGATCAGGGCGACAATCATCGCCTTATGCTGCTTCATTGCTCTCTTCTCCTTGACCTTACGGTCAGTAAGAGGCGCTCTACATGTGTTCAGCa                                                                        <  2:3099725/93‑1 (MQ=255)                 gATCAGGGCGACAATCATCGCCTTATGCTGCTTCATTGCTCTCTTCTCCTTGACCTTACGGTCAGTAAGAGGCGCTCTACATGTGTTCAGCat                                                                       <  1:1870967/93‑1 (MQ=255)                     aGGGCGACAATCATCGCCTTATGCTGCTTCATTGCTCTCTTCTCCTTGACCTTACGGTCAGTAAGAGGCGCTCTACATGTGTTCAGCATATag                                                                   <  1:742414/93‑1 (MQ=255)                     aGGGCGACAATCATCGCCTTATGCTGCTTCATTGCTCTCTTCTCCTTGACCTTACGGTCAGTAAGAGGCGCTCTACATGTGTTCAGCATATag                                                                   <  1:1774587/93‑1 (MQ=255)                      gggCGACAATCATCGCCTTATGCTGCTTCATTGCTCTCTTCTCCTTGACCTTACGGTCAGTAAGAGGCGCTCTACATGTGTTCAGCATATagg                                                                  <  2:4712415/93‑1 (MQ=255)                      gggCGACAATCATCGCCTTATGCTGCTTCATTGCTCTCTTCTCCTTGACCTTACGGTCAGTAAGAGGCGCTCTACATGTGTTCAGCATATagg                                                                  <  1:2977519/93‑1 (MQ=255)                      gggCGACAATCATCGCCTTATGCTGCTTCATTGCTCTCTTCTCCTTGACCTTACGGTCAGTAAGAGGCGCTCTACATGTGTTCAGCATATagg                                                                  <  2:4672457/93‑1 (MQ=255)                        gCGACAATCATCGCCTTATGCTGCTTCATTGCTCTCTTCTCCTTGACCTTACGGTCAGTAAGAGGCGCTCTACATGTGTTCAGCATATaggag                                                                <  1:1835833/93‑1 (MQ=255)                          gACAATCATCGCCTTATGCTGCTTCATTGCTCTCTTCTCCTTGACCTTACGGTCAGTAAGAGGCGCTCTACATGTGTTCAGCATATAGGAGGc                                                              >  1:3468737/1‑93 (MQ=255)                          gACAATCATCGCCTTATGCTGCTTCATTGCTCTCTTCTCCTTGACCTTACGGTCAGTAAGAGGCGCTCTACATGTGTTCAGCATATAGGAGGc                                                              <  2:4550814/93‑1 (MQ=255)                           aCAATCATCGCCTTATGCTGCTTCATTGCTCTCTTCTCCTTGACCTTACGGTCAGTAAGAGGCGCTCTACATGTGTTCAGCATATAGGAGGcc                                                             <  1:1168534/93‑1 (MQ=255)                           aCAATCATCGCCTTATGCTGCTTCATTGCTCTCTTCTCCTTGACCTTACGGTCAGTAAGAGGCGCTCTACATGTGTTCAGCATATAGGAGGcc                                                             <  1:2037546/93‑1 (MQ=255)                           aCAATCATCGCCTTACGCTGCTTCATTGCTCTCTTCTCCTTGACCTTACGGTCAGTAAGAGGCGCTCTACATGTGTTCAGCATATAGGAGGcc                                                             <  2:4460492/93‑1 (MQ=255)                             aaTCATCGCCTTATGCTGCTTCATTGCTCTCTTCTCCTTGACCTTACGGTCAGTAAGAGGCGCTCTACATGTGTTCAGCATATAGGAGGCCTc                                                           <  2:334301/93‑1 (MQ=255)                             aaTCATCGCCTTATGCTGCTTCATTGCTCTCTTCTCCTTGACCTTACGGTCAGTAAGAGGCGCTCTACATGTGTTCAGCATATAGGAGGCCTc                                                           <  2:4519987/93‑1 (MQ=255)                             aaTCATCGCCTTATGCTGCTTCATTGCTCTCTTCTCCTTGACCTTACGGTCAGTAAGAGGCGCTCTACATGTGTTCAGCATATAGGAGGCCTc                                                           <  1:245600/93‑1 (MQ=255)                              atcatcGCCTTATGCTGCTTCATTGCTCTCTTCTCCTTGACCTTACGGTCAGTAAGAGGCGCTCTACATGTGTTCAGCATATAGGAGGCCTCg                                                          <  2:4590327/93‑1 (MQ=255)                               tcatcGCCTTATGCTGCTTCATTGCTCTCTTCTCCTTGACCTTACGGTCAGTAAGAGGCGCTCTACATGTGTTCAGCATATAGGAGGCCTCgg                                                         <  1:383157/93‑1 (MQ=255)                                 atcGCCTTATGCTGCTTCATTGCTCTCTTCTCCTTGACCTTACGGTCAGTAAGAGGCGCTCTACATGTGTTCAGCATATAGGAGGCCTCGGGt                                                       >  2:3723177/1‑93 (MQ=255)                                 atcGCCTTATGCTGCTGCATTGCTCTCTTCTCCTTAACCTTACGGTCAGTAAGAGTCGCTCTGCATGTGTTCAGCATATAGGAGGCCTCGGGt                                                       <  1:4475192/93‑1 (MQ=255)                                  tcGCCTTATGCTGCTTCATTGCTCTCTTCTCCTTGACCTTACGGTCAGTAAGAGGCGCTCTACATGTGTTCAGCATATAGGAGGCCTCGGGtt                                                      <  2:876571/93‑1 (MQ=255)                                  tcGCCTTATGCTGCTTCATTGCTCTCTTCTCCTTGACCTTACGGTCAGTAAGAGGCGCTCTACATGTGTTCAGCATATAGGAGGCCTCGGGtt                                                      <  2:1443682/93‑1 (MQ=255)                                  tcGCCTTATGCTGCTTCATTGCTCTCTTCTCCTTGACCTTACGGTCAGTAAGAGGCGCTCTACATGTGTTCAGCATATAGGAGACCTCGGGtt                                                      >  2:2153340/1‑93 (MQ=255)                                    gCCTTATGCTGCTTCATTGCTCTCTTCTCCTTGACCTTACGGTCAGTAAGAGGCGCTCTACATGTGTTCAGCATATAGGAGGCCTCGGGTTGa                                                    >  2:2227844/1‑93 (MQ=255)                                     ccTTATGCTGCTTCATTGCTCTCTTCTCCTTGACCTTACGGTCAGTAAGAGGCGCTCTACATGTGTTCAGCATATAGGAGGCCTCGGGTTGAt                                                   <  1:2990557/93‑1 (MQ=255)                                      cTTATGCTGCTTCATTGCTCTCTTCTCCTTGACCTTACGGTCAGTAAGAGGCACTCTACATGTGTTCAGCATATAGGAGGCCTCGGGTTGATg                                                  >  1:798994/1‑93 (MQ=255)                                       ttATGCTGCTTCATTGCTCTCTTCTCCTTGACCTTACGGTCAGTAAGAGGCGCTCTACATGTGTTCAGCATATAGGAGGCCTCGGGTTGATgg                                                 >  1:4730117/1‑93 (MQ=255)                                        tATGCTTCTTCATTTCGTTCTCCTCCATGACCTTACGGTCATTAAGAGGCACTCTACATGTGTTCAGCATATAGGAGGCCTCGGGTTGATCGt                                                <  2:695018/93‑1 (MQ=255)                                          tgctgcTTCATTGCTCTCTTCTCCTTGACCTTACGGTCAGTAAGAGGCGCTCTACATGTGTTCAGCATATAGGAGGCCTCGGGTTGATGGTaa                                              <  2:3113318/93‑1 (MQ=255)                                                ttCATTGCTCTCTTCTCCTTGACCTTACGGTCAGTAAGAGGCACTCTACATGTGTTCAGCATATAGGAGGCCTCGGGTTGATGGTAAAATATc                                        <  2:356879/93‑1 (MQ=255)                                                   aTTGCTCTCTTCTCCTTGACCTTACGGTCAGTAAGAGGCGCTCTACATGTGTTCAGCATATAGGAGGCCTCGGGTTGATGGTAAAATATCACt                                     <  1:4454363/93‑1 (MQ=255)                                                   aTTGCTCTCTTCTCCTTGACCTTACGGTCAGTAAGAGGCGCTCTACATGTGTTCAGCATATAGGAGGCCTCGGGTTGATGGTAAAATATCACt                                     <  1:458831/93‑1 (MQ=255)                                                    ttGCTCTCTTCTCCTTGACCTTACGGTCAGTAAGAGGCGCTCTACATGTGTTCAGCATATAGGAGGCCTCGGGTTGATGGTAAAATATCACTc                                    <  2:3226668/93‑1 (MQ=255)                                                      tctctctTCTCCTTGACCTTACGGTCAGTAAGAGGCGCTCTACATGTGTTCAGCATATAGGAGGCCTCGGGTTGATGGTAAAATATCACTCgg                                  <  1:3894743/92‑1 (MQ=255)                                                      gCTCTCTTCTCCTTGACCTTACGGTCAGTAAGAGGCGCTCTACATGTGTTCAGCATATAGGAGGCCTCGGGTTGATGGTAAAATATCACTCgg                                  <  1:4460115/93‑1 (MQ=255)                                                      gCTCTCTTCTCCTTGACCTTACGGTCAGTAAGAGGCGCTCTACATGTGTTCAGCATATAGGAGGCCTCGGGTTGATGGTAAAATATCACTCgg                                  <  2:4681353/93‑1 (MQ=255)                                                       ctctctTCTCCTTGACCTTACGGTCAGTAAGAGGCGCTCTACATGTGTTCAGCATATAGGAGGCCTCGGGTTGATGGTAAAATATCACTCggg                                 <  2:4578663/93‑1 (MQ=255)                                                       ctctctTCTCCTTGACCTTACGGTCAGTAAGAGGCGCTCTACATGTGTTCAGCATATAGGAGGCCTCGGGTTGATGGTAAAATATCACTCggg                                 <  2:3747500/93‑1 (MQ=255)                                                        tctgtTCTCCTTGACCTTACGGTCAGTAAGAGGCGCTCTACATGTGTTCAGCATATAGGAGGCCTCGGGTTGATGGTAAAATATCACTCgggg                                <  2:3825918/93‑1 (MQ=255)                                                        tctctTCTCCTTGACCTTACGGTCAGTAAGAGGCGCTCTACATGTGTTCAGCATATAGGAGGCCTCGGGTTGATGGTAAAATATCACTCgggg                                >  1:593312/1‑93 (MQ=255)                                                        tctctTCTCCTTGACCTTACGGTCAGTAAGAGGCGCTCTACATGTGTTCAGCATATAGGAGGCCTCGGGTTGATGGTAAAATATCACTCgggg                                <  2:1285262/93‑1 (MQ=255)                                                         ctctTCTCCTTGACCTTACGGTCAGTAAGAGGCGCTCTACATGTGTTCAGCATATAGGAGGCCTCGGGTTGATGGTAAAATATCACTCGGGGc                               <  2:1197815/93‑1 (MQ=255)                                                          tcttctCCTTGACCTTACGGTCAGTAAGAGGCGCTCTACATGTGTTCAGCATATAGGAGGCCTCGGGTTGATGGTAAAATATCACTCGGGGCt                              <  2:502595/93‑1 (MQ=255)                                                             tctcCTTGACCTTACGGTCAGTAAGAGGCGCTCTACATGTGTTCAGCATATAGGAGGCCTCGGGTTGATGGTAAAATATCACTCGGGGCtttt                           <  2:938823/93‑1 (MQ=255)                                                              ctcCTTGACCTTACGGTCAGTAAGAGGCGCTCTACATGTGTTCAGCATATAGGAGGCCTCGGGTTGATGGTAAAATATCACTCGGGGCTTTtc                          <  1:3951850/93‑1 (MQ=255)                                                              ctcCTTGACCTTACGGTCAGTAAGAGGCGCTCTACATGTGTTCAGCATATAGGAGGCCTCGGGTTGATGGTAAAATATCACTCGGGGCTTTtc                          <  1:2328539/93‑1 (MQ=255)                                                               tcCTTGACCTTACGGTCAGTAAGAGGCGCTCTACATGTGTTCAGCATATAGGAGGCCTCGGGTTGATGGTAAAATATCACTCGGGGCTTTtct                         <  2:1557093/93‑1 (MQ=255)                                                                ccTTGACCTTACGGTCAGTAAGAGGCGCTCTACATGTGTTCAGCATATAGGAGGCCTCGGGTTGATGGTAAAATATCACTCGGGGCTTTtctc                        <  1:768388/93‑1 (MQ=255)                                                                    gACCTTACGGTCAGTAAGAGGCGCTCTACATGTGTTCAGCATATAGGAGGCCTCGGGTTGATGGTAAAATATCACTCGGGGCTTTTCTCTATc                    <  1:2698048/93‑1 (MQ=255)                                                                    gACCTTACGGTCAGTAAGAGGCGCTCTACATGTGTTCAGCATATAGGAGGCCTCGGGTTGATGGTAAAATATCACTCGGGGCTTTTCTCTATc                    <  1:1524544/93‑1 (MQ=255)                                                                        ttACGGTCAGTAAGAGGCGCTCTACATGTGTTCAGCATATAGGAGGCCTCGGGTTGATGGTAAAATATCACTCGGGGCTTTTCTCTATCTGcc                <  2:1358539/93‑1 (MQ=255)                                                                            ggTCAGTAAGAGGCGCTCTACATGTGTTCAGCATATAGGAGGCCTCGGGTTGATGGTAAAATATCACTCGGGGCTTTTCTCTATCTGCCGTTc            <  1:2919034/93‑1 (MQ=255)                                                                             gTCAGTAAGAGGCGCTCTACATGTGTTCAGCATATAGGAGGCCTCGGGTTGATGGTAAAATATCACTCGGGGCTTTTCTCTATCTGCCGTTCa           <  2:903426/93‑1 (MQ=255)                                                                                aGTAAGAGGCGCTCTACATGTGTTCAGCATATAGGAGGCCTCGGGTTGATGGTAAAATATCACTCGGGGCTTTTCTCTATCTGCCGTTCAGCt        <  2:2378595/93‑1 (MQ=255)                                                                                   aaGAGGCGCTCTACATGTGTTCAGCATATAGGAGGCCTCGGGTTGATGGTAAAATATCACTCGGGGCTTTTCTCTATCTGCCGTTCAGCTAAt     <  2:2357653/93‑1 (MQ=255)                                                                                     gagGCGCTCTACATGTGTTCAGCATATAGGGGGCCTCGGGTTGATGGTAAAATATCACTCGGGGCTTTTCTCTATCTGCCGTTCAGCTAATGc   <  1:1491164/93‑1 (MQ=255)                                                                                      agGCGCTCTACATGTGTTCAGCATAGAGGAGTCCTCGGGTTGATGGTAAAATATCACTCGGGGCTTTTCTCTATCTGCCGTTCAGCTAATGcc  <  2:3333454/93‑1 (MQ=255)                                                                                          |                                                                                         CGGTGATACAGATGACGATCAGGGCGACAATCATCGCCTTATGCTGCTTCATTGCTCTCTTCTCCTTGACCTTACGGTCAGTAAGAGGCACTCTACATGTGTTCAGCATATAGGAGGCCTCGGGTTGATGGTAAAATATCACTCGGGGCTTTTCTCTATCTGCCGTTCAGCTAATGCC  >  NC_000913/16852‑17029

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