Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_007793	2,580,941	C→T	T40M (ACG→ATG)	SAUSA300_RS13270 →	carboxylesterase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_007793	2,580,941	0	C	T	97.9%	132.1 / ‑5.0	48	T40M (ACG→ATG)	SAUSA300_RS13270	carboxylesterase
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base T (39/8);  total (40/8)
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.86e-01

GGTCAAATTCATGGTATTACACAAGATGGTTTAGATATCTTCTTAGGCATTCCTTATGCAGAACCACCAGTTCATGACAATCGCTTTAAACATTCTACGTTAAAAACACAATGGTCAGAGCCAATTGATGCAACTGAAATACAACCCATCCCACCGCAACCAGACAACAAATTAGAAGATTTTTTCTCCTCAC  >  NC_007793/2580844‑2581036                                                                                                  |                                                                                                ggtCAAATTCATGGTATTACACAAGATGGTTTAGATATCTTCTTAGGCATTCCTTATGCAGAACCACCAGTTCATGACAATCGCTTTAAACATTCTATGt                                                                                               >  1:1055072/1‑100 (MQ=255) ggtCAAATTCATGGTATTACACAAGATGGTTTAGATATCTTCTTAGGCATTCCTTATGCAGAACCACCAGTTCATGACAATCGCTTTAAACATTCTATGt                                                                                               >  1:994986/1‑100 (MQ=255) ggtCAAATTCATGGTATTACACAAGATGGTTTAGATATCTTCTTAGGCATTCCTTATGCAGAACCACCAGTTCATGACAATCGCTTTAAACATTCTATGt                                                                                               >  1:2692492/1‑100 (MQ=255) ggtCAAATTCATGGTATTACACAAGATGGTTTAGATATCTTCTTAGGCATTCCTTATGCAGAACCACCAGTTCATGACAATCGCTTTAAACATTCTATGt                                                                                               >  1:2258304/1‑100 (MQ=255) ggtCAAATTCATGGTATTACACAAGATGGTTTAGATATCTTCTTAGGCATTCCTTATGCAGAACCACCAGTTCATGACAATCGCTTTAAACATTCTATGt                                                                                               >  1:1460204/1‑100 (MQ=255)        ttCATGGTATTACACAAGATGGTTTAGATATCTTCTTAGGCATTCCTTATGCAGAACCACCAGTTCATGACAATCGCTTTAAACATTCTATGTTAAAAac                                                                                        >  1:54594/1‑100 (MQ=255)             ggTATTACACAAGATGGTTTAGATATCTTCTTAGGCATTCATTATGCAGAACCACAAGTTCATGACAATCGCTTTAAACATTCTATGTTAAAAACAAAAt                                                                                   >  1:2057765/1‑100 (MQ=255)                aTTACACAAGATGGTTTAGATATCTTCTTAGGCATTCCTTATGCAGAACCACCAGTTCATGACAATCGCTTTAAACATTCTATGTTAAAAACACAATGGt                                                                                >  1:4908320/1‑100 (MQ=255)                  tACACAAGATGGTTTAGATATCTTCTTAGGCATTCCTTATGCAGAACCACCAGTTCATGACAATCGCTTTAAACATTCTATGTTAAAAACACaagatcgg                                                                              >  1:3489899/1‑94 (MQ=255)                  tACACAAGATGGTTTAGATATCTTCTTAGGCATTCCTTATGCAGAACCACCAGTTCATGACAATCGCTTTAAACATTCTATGTTAAAAACACAATGGTCa                                                                              >  1:580338/1‑100 (MQ=255)                  tACACAAGATGGTTTAGATATCTTCTTAGGCATTCCTTATGCAGAACCACCAGTTCATGACAATCGCTTTAAACATTCTATGTTAAAAACACAATGGTCa                                                                              >  1:1355469/1‑100 (MQ=255)                   acacAAGATGGTTTAGATATCTTCTTAGGCATTCCTTATGCAGAACCACCAGTTCATGACAATCGCTTTAAACATTCTATGTTAAAAAcacagatcggaa                                                                             >  1:148431/1‑92 (MQ=255)                   acacAAGATGGTTTAGATATCTTCTTAGGCATTCCTTATGCAGAACCACCAGTTCATGACAATCGCTTTAAACATTCTATGTTAAAAACACAATGGTCag                                                                             >  1:855507/1‑100 (MQ=255)                   acacAAGATGGTTTAGATATCTTCTTAGGCATTCCTTATGCAGAACCACCAGTTCATGACAATCGCTTTAAACATTCTATGTTAAAAACACAATGGTCag                                                                             >  1:1378906/1‑100 (MQ=255)                   acacAAGATGGTTTAGATATCTTCTTAGGCATTCCTTATGCAGAACCACCAGTTCATGACAATCGCTTTAAACATTCTATGTTAAAAACACAATGGTCag                                                                             >  1:466659/1‑100 (MQ=255)                         caTGGTTTAGATATCTTCTTAGGCATTCCTTATGCAGAACCACCAGTTCATGACAATCGCTTTAAACATTCTATGTTAAAAACACAATGGTCAGAGCCaa                                                                       >  1:4122268/2‑100 (MQ=255)                                       cttcttAGGCATTCCTTATGCAGAACCACCAGTTCATGACAATCGCTTTAAACATTCTATGTTAAAAACACAATGGTCAGAGCCAATTGATGCAACTGaa                                                         >  1:2363366/1‑100 (MQ=255)                                         tcttAGGCATTCCTTATGCAGAACCACCAGTTCATGACAATCGCTTTAAACATTCTATGTTAAAAACACAATGGTCAGAGCCAATTGATGCAACTGAAAt                                                       >  1:3633042/1‑100 (MQ=255)                                         tcttAGGCATTCCTTATGCAGAACCACCAGTTCATGACAATCGCTTTAAACATTCTATGTTAAAAACACAATGGTCAGAGCCAATTGATGCAACTGAAAt                                                       >  1:986436/1‑100 (MQ=255)                                         tcttAGGCATTCCTTATGCAGAACCACCAGTTCATGACAATCGCTTTAAACATTCTATGTTAAAAACACAATGGTCAGAGCCAATTGATGCAACTGAAAt                                                       >  1:1388561/1‑100 (MQ=255)                                          cttAGGCATTCCTTATGCAGAACCACCAGTTCATGACAATCGCTTTAAACATTCTATGTTAAAAACACAATGGTCAGAGCCAATTGATGCAACTGAAATa                                                      >  1:2072481/1‑100 (MQ=255)                                          cttAGGCATTCCTTATGCAGAACCACCAGTTCATGACAATCGCTTTAAACATTCTATGTTAAAAACACAATGGTCAGAGCCAATTGATGCAACTGAAATa                                                      >  1:197703/1‑100 (MQ=255)                                          cttAGGCATTCCTTATGCAGAACCACCAGTTCATGACAATCGCTTTAAACATTCTACGTTAAAAACACAATGGTCAGAGCCAATTGATGCAACTGAAATa                                                      >  1:4905089/1‑100 (MQ=255)                                            tAGGCATTCCTTATGCAGAACCACCAGTTCATGACAATCGCTTTAAACATTCTATGTTAAAAACACAATGGTCAGAGCCAATTGATGCAACTGAAATACa                                                    >  1:3138441/1‑100 (MQ=255)                                             aGGCATTCCTTATGCAGAACCACCAGTTCATGACAATCGCTTTAAACATTCTATGTTAAAAACACAATGGTCAGAGCCAATTGATGCAACTGAAATACaa                                                   >  1:4687967/1‑100 (MQ=255)                                                    ccTTATGCAGAACCACCAGTTCATGACAATCGCTTTAAACATTCTATGTTAAAAACACAATGGTCAGAGCCAATTGATGCAACTGAAATACAACCCATcc                                            >  1:3328514/1‑100 (MQ=255)                                                       tATGCAGAACCACCAGTTCATGACAATCGCTTTAAACATTCTATGTTAAAAACACAATGGTCAGAGCCAATTGATGCAACTGAAATACAACCCATCCCat                                         <  1:4298057/100‑2 (MQ=255)                                                       tATGCAGAACCACCAGTTCATGACAATCGCTTTAAACATTCTATGTTAAAAACACAATGGTCAGAGCCAATTGATGCAACTGAAATACAACCCATCCCAc                                         >  1:651064/1‑100 (MQ=255)                                                       tATGCAGAACCACCAGTTCATGACAATCGCTTTAAACATTCTATGTTAAAAACACAATGGTCAGAGCCAATTGATGCAACTGAAATACAACCCATCCCAc                                         >  1:1954272/1‑100 (MQ=255)                                                       tATGCAGAACCACCAGTTCATGACAATCGCTTTAAACATTCTATGTTAAAAACACAATGGTCAGAGCCAATTGATGCAACTGAAATACAACCCATCCCAc                                         >  1:3265223/1‑100 (MQ=255)                                                       caTGCAGAACCACCAGTTCATGACAATCGCTTTAAACATTCTATGTTAAAAACACAATGGTCAGAGCCAATTGATGCAACTGAAATACAACCCATCCCAc                                         >  1:4126530/2‑100 (MQ=255)                                                           cAGAACCACCAGTTCATGACAATCGCTTTAAACATTCTATGTTAAAAACACAATGGTCAGAGCCAATTGATGCAACTGGAATACAACCCATCCCAcccag                                     <  1:3686712/100‑4 (MQ=255)                                                            tgAACCACCAGTTCATGACAATCGCTTTAAACATTCTATGTTAAAAACACAATGGTCAGAGCCAATTGATGCAACTGAAATACAACCCATCCCACCGCag                                    >  1:2056879/2‑99 (MQ=255)                                                             gAACCACCAGTTCATGACAATCGCTTTAAACATTCTATGTTAAAAACACAATGGTCAGAGCCAATTGATGCAACTGAAATACAACCCATCCCACCGCAAc                                   >  1:2810740/1‑100 (MQ=255)                                                             gAACCACCAGTTCATGACAATCGCTTTAAACATTCTATGTTAAAAACACAATGGTCAGAGCCAATTGATGCAACTGAAATACAACCCATCCCACCGCAAc                                   >  1:4367145/1‑100 (MQ=255)                                                             aaaCCACCAGTTCATGACAATCGCTTTAAACATTCTATGTTAAAAACACAATGGTCAGAGCCAATTGATGCAACTGAAATACAACCCATCCCACCGCAAc                                   >  1:5117168/2‑100 (MQ=255)                                                              aaCCACCAGTTCATGACAATCGCTTTAAACATTCTATGTTAAAAACACAATGGTCAGAGCCAATTGATGCAACTGAAATACAACCCATCCCACCGCAAcc                                  <  1:5199117/100‑1 (MQ=255)                                                              aaCCACCAGTTCATGACAATCGCTTTAAACATTCTATGTTAAAAACACAATGGTCAGAGCCAATTGATGCAACTGAAATACAACCCATCCCACCGCAAcc                                  >  1:1025180/1‑100 (MQ=255)                                                              aaCCACCAGTTCATGACAATCGCTTTAAACATTCTATGTTAAAAACACAATGGTCAGAGCCAATTGATGCAACTGAAATACAACCCATCCCACCGCAAcc                                  <  1:3085259/100‑1 (MQ=255)                                                              aaCCACCAGTTCATGACAATCGCTTTAAACATTCTATGTTAAAAACACAATGGTCAGAGCCAATTGATGCAACTGAAATACAACCCATCCCACCGCAAcc                                  <  1:3458521/100‑1 (MQ=255)                                                               tccaccaGTTCATGACAATCGCTTTAAACATTCTATGTTAAAAACACAATGGTCAGAGCCAATTGATGCAACTGAAATACAACCCATCCCACCGCAACCa                                 >  1:3494835/2‑100 (MQ=255)                                                               accaccAGTTCATGACAATCGCTTTAAACATTCTATGTTAAAAACACAATGGTCAGAGCCAATTGATGCAACTGAAATACAACCCATCCCACCGCAACCa                                 >  1:452623/1‑100 (MQ=255)                                                               accaccAGTTCATGACAATCGCTTTAAACATTCTATGTTAAAAACACAATGGTCAGAGCCAATTGATGCAACTGAAATACAACCCATCCCACCGCAACCa                                 >  1:4001797/1‑100 (MQ=255)                                                               accaccAGTTCATGACAATCGCTTTAAACATTCTATGTTAAAAACACAATGGTCAGAGCCAATTGATGCAACTGAAATACAACCCATCCCACCGCAACCa                                 <  1:2039311/100‑1 (MQ=255)                                                                ccaccaGTTCATGACAATCGCTTTAAACATTCTATGTTAAAAACACAATGGTCAGAGCCAATTGATGCAACTGAAATACAACCCATCCCACCGCAACCAg                                >  1:4212328/1‑100 (MQ=255)                                                                ccaccaGTTCATGACAATCGCTTTAAACATTCTATGTTAAAAACACAATGGTCAGAGCCAATTGATGCAACTGAAATACAACCCATCCCACCGCAACCAg                                <  1:3394544/100‑1 (MQ=255)                                                                          aTGACAATCGCTTTAAACATTCTATGTTAAAAACACAATGGTCAGAGCCAATTGATGCAACTGAAATACAACCCAGCCCACCGCAACCAGACAACAAAtt                      <  1:510383/100‑1 (MQ=255)                                                                                              tCTATGTTAAAAACACAATGGTCAGAGCCAATTGATGCAACTGAAATACAACCCATCCCACCGCAACCAGACAACAAATTAGAAGATTTTTTCTCCTcac  >  1:2755658/1‑100 (MQ=255)                                                                                                  |                                                                                                GGTCAAATTCATGGTATTACACAAGATGGTTTAGATATCTTCTTAGGCATTCCTTATGCAGAACCACCAGTTCATGACAATCGCTTTAAACATTCTACGTTAAAAACACAATGGTCAGAGCCAATTGATGCAACTGAAATACAACCCATCCCACCGCAACCAGACAACAAATTAGAAGATTTTTTCTCCTCAC  >  NC_007793/2580844‑2581036

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