Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_007793	2,142,339	C→A	L119I (CTA→ATA)	SAUSA300_RS10910 →	CAAX amino protease

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_007793	2,142,339	0	C	A	98.8%	240.4 / ‑5.7	86	L119I (CTA→ATA)	SAUSA300_RS10910	CAAX amino protease
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base A (84/0);  minor base G (1/0);  total (86/0)
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.04e-01

GCTTGGATTTGCATTGGTAATGATTTATCAAATGGTAGTGAGCATTGTATTAACGCAAATTTATGGTGGACAACAAGTAAGTCCTAATACAGAAAAGCTAATTATTATTGCTCGAAAAATACCTATATTTATCTTCTTTGTATCTATTATTGGTCCTTTATTAGAAGAATATGTATTCAGAAAAGTAA  >  NC_007793/2142242‑2142429                                                                                                  |                                                                                           gCTTGGATTTGCATTGGTAATGATTTATCAAATGGTAGTGAGCATTGTATTAACGCAAATTTATGGTGGACAACAAGTAAGTCCTAATACAGAAAAGATa                                                                                          >  1:1507886/1‑100 (MQ=255) gCTTGGATTTGCATTGGTAATGATTTATCAAATGGTAGGGAGCATTGTATTAACGCAAATTTATGGTGGACAACAAGTAAGTCCTAATACAGAAAAGGTa                                                                                          >  1:6296330/1‑100 (MQ=255)          tGCATTGGTAATGATTTATCAAATGGTAGTGAGCATTGTATTAACGCAAATTTATGGTGGACAACAAGTAAGTCCTAATACAGAAAAGATAattattatt                                                                                 >  1:6265982/1‑100 (MQ=255)           gCATTGGTAATGATTTATCAAATGGTAGTGAGCATTGTATTAACGCAAATTTATGGTGGACAACAAGTAAGTCCTAATACAGAAAAGATAATTATTATTg                                                                                >  1:2718888/1‑100 (MQ=255)              ttGGTAATGATTTATCAAATGGTAGTGAGCATTGTATTAACGCAAATTTATGGTGGACAACAAGTAAGTCCTAATACAGAAAAGATAATTATTATTGCTc                                                                             >  1:4487411/1‑100 (MQ=255)               tGGTAATGATTTATCAAATGGTAGTGAGCATTGTATTAACGCAAATTTATGGTGGACAACAAGTAAGTCCTAATACAGAAAAGATAATTATTATTGCTCg                                                                            >  1:1410534/1‑100 (MQ=255)                    actgcttATCAAATGGTAGTGAGCATTGTATTAACGCAAATTTATGGTGGACAACAAGTAAGTCCTAATACAGAAAAGATAATTATTATTGCTCGaaaaa                                                                       >  1:2381462/6‑100 (MQ=255)                      gATTTATCAAATGGTAGTGAGCATTGTATTAACGCAAATTTATGGTGGACAACAAGTAAGTCCTAATACAGAAAAGATAATTATTATTGCTCGAAAAATa                                                                     >  1:2770428/1‑100 (MQ=255)                          tATCAAATGGTAGTGAGCATTGTATTAACGCAAATTTATGGAGGACAACAAGTAAGTCCTAATACAGAAAAGATAATTATTATTGCTCGAAAAATACCta                                                                 >  1:3469281/1‑100 (MQ=255)                           aTCAAATGGTAGTGAGCATTGTATTAACGCAAATTTATGGTGGACAACAAGTAAGTCCTAATACAGAAAAGATAATTATTATTGCTCGAAAAATACCtat                                                                >  1:6442632/1‑100 (MQ=255)                           aTCAAATGGTAGTGAGCATTGTATTAACGCAAATTTATGGTGGACAACAAGTAAGTCCTAATACAGAAAAGATAATTATTATTGCTCGAAAAATACCtat                                                                >  1:2183629/1‑100 (MQ=255)                             cAAATGGTAGTGAGCATTGTATTAACGCAAATTTATGGTGGACAACAAGTAAGTCCTAATACAGAAAAGATAATTATTATTGCTCGAAAAATACCTatat                                                              >  1:64320/1‑100 (MQ=255)                             cAAATGGTAGTGAGCATTGTATTAACGCAAATTTATGGTGGACAACAAGTAAGTCCTAATACAGAAAAGATAATTATTATTGCTCGAAAAATACCTatat                                                              >  1:4766649/1‑100 (MQ=255)                             cAAATGGTAGTGAGCATTGTATTAACGCAAATTTATGGTGGACAACAAGTAAGTCCTAATACAGAAAAGATAATTATTATTGCTCGAAAAATACCTatat                                                              >  1:99134/1‑100 (MQ=255)                             cAAATGGTAGTGAGCATTGTATTAACGCAAATTTATGGTGGACAACAAGTAAGTCCTAATACAGAAAAGATAATTATTATTGCTCGAAAAATACCTatat                                                              >  1:4441335/1‑100 (MQ=255)                             cAAATGGTAGTGAGCATTGTATTAACGCAAATTTATGGTGGACAACAAGTAAGTCCTAATACAGAAAAGATAATTATTATTGCTCGAAAAATACCTatat                                                              >  1:4268785/1‑100 (MQ=255)                             cAAATGGTAGTGAGCATTGTATTAACGCAAATTTATGGTGGACAACAAGTAAGTCCTAATACAGAAAAGATAATTATTATTGCTCGAAAAATACCTatat                                                              >  1:3737949/1‑100 (MQ=255)                             cAAATGGTAGTGAGCATTGTATTAACGCAAATTTATGGTGGACAACAAGTAAGTCCTAATACAGAAAAGATAATTATTATTGCTCGAAAAATACCTatat                                                              >  1:1700017/1‑100 (MQ=255)                             cAAATGGTAGTGAGCATTGTATTAACGCAAATTTATGGTGGACAACAAGTAAGTCCTAATACAGAAAAGATAATTATTATTGCTCGAAAAATACCTatat                                                              >  1:2523757/1‑100 (MQ=255)                             atgaTGGTAGTGAGCATTGTATTAACGCAAATTTATGGTGGACAACAAGTAAGTCCTAATACAGAAAAGATAATTATTATTGCTCGAAAAATACCTatat                                                              >  1:6250220/4‑100 (MQ=255)                              taaTGGTAGTGAGCATTGTATTAACGCAAATTTATGGTGGACAACAAGTAAGTCCTAATACAGAAAAGATAATTATTATTGCTCGAAAAATACCTATAtt                                                             >  1:5483123/2‑100 (MQ=255)                              taaTGGTAGTGAGCATTGTATTAACGCAAATTTATGGTGGACAACAAGTAAGTCCTAATACAGAAAAGATAATTATTATTGCTCGAAAAATACCTATAtt                                                             >  1:329142/2‑100 (MQ=255)                              ctaTGGTAGTGAGCATTGTATTAACGCAAATTTATGGTGGACAACAAGTAAGTCCTAATACAGAAAAGATAATTATTATTGCTCGAAAAATACCTATAtt                                                             >  1:1249626/3‑100 (MQ=255)                              caaTGGTAGTGAGCATTGTATTAACGCAAATTTATGGTGGACAACAAGTAAGTCCTAATACAGAAAAGATAATTATTATTGCTCGAAAAATACCTATAtt                                                             >  1:2377883/2‑100 (MQ=255)                              aaaTGGTAGTGAGCATTGTATTAACGCAAATTTATGGTGGACAACAAGTAAGTCCTAATACAGAAAAGATAATTATTATTGCTCGAAAAATACCTATAtt                                                             >  1:6277774/1‑100 (MQ=255)                              aaaTGGTAGTGAGCATTGTATTAACGCAAATTTATGGTGGACAACAAGTAAGTCCTAATACAGAAAAGATAATTATTATTGCTCGAAAAATACCTATAtt                                                             >  1:6171950/1‑100 (MQ=255)                              aaaTGGTAGTGAGCATTGTATTAACGCAAATTTATGGTGGACAACAAGTAAGTCCTAATACAGAAAAGATAATTATTATTGCTCGAAAAATACCTATAtt                                                             >  1:1576269/1‑100 (MQ=255)                              aaaTGGTAGTGAGCATTGTATTAACGCAAATTTATGGTGGACAACAAGTAAGTCCTAATACAGAAAAGATAATTATTATTGCTCGAAAAATACCTATAtt                                                             >  1:3396508/1‑100 (MQ=255)                              aaaTGGTAGTGAGCATTGTATTAACGCAAATTTATGGTGGACAACAAGTAAGTCCTAATACAGAAAAGATAATTATTATTGCTCGAAAAATACCTATAtt                                                             >  1:3941594/1‑100 (MQ=255)                               taTGGTAGTGAGCATTGTATTAACGCAAATTTATGGTGGACAACAAGTAAGTCCTAATACAGAAAAGATAATTATTATTGCTCGAAAAATACCTATAttt                                                            >  1:3100837/2‑100 (MQ=255)                               gaTGGTAGTGAGCATTGTATTAACGCAAATTTATGGTGGACAACAAGTAAGTCCTAATACAGAAAAGATAATTATTATTGCTCGAAAAATACCTATAttt                                                            >  1:2317206/2‑100 (MQ=255)                               agtGGTAGTGAGCATTGTATTAACGCAAATTTATGGTGGACAACAAGTAAGTCCTAATACAGAAAAGATAATTATTATTGCTCGACAAATACCTATAttt                                                            >  1:4943030/3‑100 (MQ=255)                               aaTGGTAGTGAGCATTGTCTTAACGCAAATTTATGGTGGACAACAAGTAAGTCCTAATACAGAAAAGATAATTATTATTGCTCGAAAAATACCTATAttt                                                            >  1:3364093/1‑100 (MQ=255)                               aaTGGTAGTGAGCATTGTATTAACGCAAATTTATGGTGGACAACAAGTAAGTCCTAATGCAGAAAAGATAATTATTATTGCTCGAAAAATACCTATAttt                                                            >  1:2516145/1‑100 (MQ=255)                               aaTGGTAGTGAGCATTGTATTAACGCAAATTTATGGTGGACAACAAGTAAGTCCTAATACAGAAAAGATAATTATTATTGCTCGACAAATACCTATAttt                                                            >  1:4890566/1‑100 (MQ=255)                               aaTGGTAGTGAGCATTGTATTAACGCAAATTTATGGTGGACAACAAGTAAGTCCTAATACAGAAAAGATAATTATTATTGCTCGAAAAATACCTATAttt                                                            >  1:802583/1‑100 (MQ=255)                               aaTGGTAGTGAGCATTGTATTAACGCAAATTTATGGTGGACAACAAGTAAGTCCTAATACAGAAAAGATAATTATTATTGCTCGAAAAATACCTATAttt                                                            >  1:3428217/1‑100 (MQ=255)                               aaTGGTAGTGAGCATTGTATTAACGCAAATTTATGGTGGACAACAAGTAAGTCCTAATACAGAAAAGATAATTATTATTGCTCGAAAAATACCTATAttt                                                            >  1:1508932/1‑100 (MQ=255)                               aaTGGTAGTGAGCATTGTATTAACGCAAATTTATGGTGGACAACAAGTAAGTCCTAATACAGAAAAGATAATTATTATTGCTCGAAAAATACCTATAttt                                                            >  1:548496/1‑100 (MQ=255)                               aaTGGTAGTGAGCATTGTATTAACGCAAATTTATGGTGGACAACAAGTAAGTCCTAATACAGAAAAGATAATTATTATTGCTCGAAAAATACCTATAttt                                                            >  1:4684705/1‑100 (MQ=255)                               aaTGGTAGTGAGCATTGTATTAACGCAAATTTATGGTGGACAACAAGTAAGTCCTAATACAGAAAAGATAATTATTATTGCTCGAAAAATACCTATAttt                                                            >  1:4660107/1‑100 (MQ=255)                               aaTGGTAGTGAGCATTGTATTAACGCAAATTTATGGTGGACAACAAGTAAGTCCTAATACAGAAAAGATAATTATTATTGCTCGAAAAATACCTATAttt                                                            >  1:2473961/1‑100 (MQ=255)                               aaTGGTAGTGAGCATTGTATTAACGCAAATTTATGGTGGACAACAAGTAAGTCCTAATACAGAAAAGATAATTATTATTGCTCGAAAAATACCTATAttt                                                            >  1:269929/1‑100 (MQ=255)                               aaTGGTAGTGAGCATTGTATTAACGCAAATTTATGGTGGACAACAAGTAAGTCCTAATACAGAAAAGATAATTATTATTGCTCGAAAAATACCTATAttt                                                            >  1:1913428/1‑100 (MQ=255)                                aTGGTAGTGAGCATTGTATTAACGCAAATTTATGGTGGACAACAAGTAAGTCCTAATACAGAAAAGATAATTATTATTGCTCGAAAATTACCTATATTTa                                                           >  1:321223/1‑100 (MQ=255)                                aTGGTAGTGAGCATTGTATTAACGCAAATTTATGGTGGACAACAAGTAAGTCCTAATACAGAAAAGATAATTATTATTGCTCGAAAAATACCTATATTTa                                                           >  1:5621493/1‑100 (MQ=255)                                aTGGTAGTGAGCATTGTATTAACGCAAATTTATGGTGGACAACAAGTAAGTCCTAATACAGAAAAGATAATTATTATTGCTCGAAAAATACCTATATTTa                                                           >  1:6627862/1‑100 (MQ=255)                                  ggTAGTGAGCATTGTATTAACGCAAATTTATGGTGGACAACAAGTAGGTCCTAATACAGAAAAGATAATTATTATTACTCGAAAAATACCTATATTTAtc                                                         >  1:5805656/1‑100 (MQ=255)                                  ggTAGTGAGCATTGTATTAACGCAAATTTATGGTGGACAACAAGTAAGTCCTAATACAGAAAAGATAATTATTATTGCTCGAAAAGTACCTATATTTAtc                                                         >  1:4961484/1‑100 (MQ=255)                                  ggTAGTGAGCATTGTATTAACGCAAATTTATGGTGGACAACAAGTAAGTCCTAATACAGAAAAGATAATTATTATTGCTCGAAAAATACCTATATTTAtc                                                         >  1:5961055/1‑100 (MQ=255)                                   gTAGTGAGCATTGTATTAACGCAAATTTATGGTGGACAACAAGTAAGTCCTAATACAGAAAAGATAATTATTATTGCTCGAAAAATACCTATATTTAtct                                                        >  1:3370429/1‑100 (MQ=255)                                     aGTGAGCATTGTATTAACGCAAATTTATGGTGGACAACAAGTAAGTCCTAATACAGAAAAGATAATTATTATTGCTCGAAAAATACCTATATTTAtcttc                                                      >  1:4220457/1‑100 (MQ=255)                                          tcATTGTATTAACGCAAATTTATGGTGGACAACAAGTAAGTCCTAATACAGAAAAGATAATTATTATTGCTCGAAAAATACCTATATTTATCTTCTTTGt                                                 >  1:788001/2‑100 (MQ=255)                                            aTTGTATTAACGCAAATTTATGGTGGACAACAAGTAAGTCCTAATACAGAAAAGATAATTATTATTGCTCGAAAAATACCTATATTTATCTTCTTTGTAt                                               >  1:4909650/1‑100 (MQ=255)                                            aTTGTATTAACGCAAATTTATGGTGGACAACAAGTAAGTCCTAATACAGAAAAGATAATTATTATTGCTCGAAAAATACCTATATTTATCTTCTTTGTAt                                               >  1:2703362/1‑100 (MQ=255)                                            aTTGTATTAACGCAAATTTATGGTGGACAACAAGTAAGTCCTAATACAGAAAAGATAATTATTATTGCTCGAAAAATACCTATATTTATCTTCTTTGTAt                                               >  1:4310793/1‑100 (MQ=255)                                              tGTATTAACGCAAATTTATGGTGGACAACAAGTAAGTCCTAATACAGAAAAGATAATTATTATTGCTCGAAAAATACCTATATTTATCTTCTTTGTATCt                                             >  1:4331268/1‑100 (MQ=255)                                              tGTATTAACGCAAATTTATGGTGGACAACAAGTAAGTCCTAATACAGAAAAGATAATTATTATTGCTCGAAAAATACCTATATTTATCTTCTTTGTATCt                                             >  1:806019/1‑100 (MQ=255)                                              tGTATTAACGCAAATTTATGGTGGACAACAAGTAAGTCCTAATACAGAAAAGATAATTATTATTGCTCGAAAAATACCTATATTTATCTTCTTTGTATCt                                             >  1:894306/1‑100 (MQ=255)                                                tATTAACGCAAATTTATGGTGGACAACAAGTAAGTCCTAATACAGAAAAGATAATTATTATTGCTCGAAAAATACCTATATTTATCTTCTTTGTATCtat                                           >  1:4912575/1‑100 (MQ=255)                                                 aTTAACGCAAATTTATGGTGGACAACAAGTAAGTCCTAATACAGAAAAGTTAATTATTATTGCTCGAAAAATACCTATATTTATCTTCTTTGTATCtatt                                          >  1:1926356/1‑100 (MQ=255)                                                  tcaaCGCAAATTTATGGTGGACAACAAGTAAGTCCTAATACAGAAAAGATAATTATTATTGCTCGAAAAATACCTATATTTATCTTCTTTGTATCtatta                                         >  1:4186028/3‑100 (MQ=255)                                                   tAACGCAAATTTATGGTGGACAACAAGTAAGTCCTAATACAGAAAAGATAATTATTATTGCTCGAAAAATACCTATATTTATCTTCTTTGTATCtattat                                        >  1:867747/1‑100 (MQ=255)                                                   tAACGCAAATTTATGGTGGACAACAAGTAAGTCCTAATACAGAAAAGATAATTATTATTGCTCGAAAAATACCTATATTTATCTTCTTTGTATCtattat                                        >  1:3915788/1‑100 (MQ=255)                                                   caaCGCAAATTTATGGTGGACAACAAGTAAGTCCTAATACAGAAAAGATAATTATTATTGCTCGAAAAATACCTATATTTATCTTCTTTGTATCtattat                                        >  1:6455564/2‑100 (MQ=255)                                                    aaCGCAAATTTATGGTGGACAACAAGTAAGTCCTAATACAGAGAAGATAATTATTATTGCTCGAAAAATACCTATATTTATCTTCTTTGTATCTattatt                                       >  1:6414433/1‑100 (MQ=255)                                                    aaCGCAAATTTATGGTGGACAACAAGTAAGTCCTAATACAGAAAAGATAATTATTATTGCTCGAAAAATACCTATATTTATCTTCTTTGTATCTattatt                                       >  1:3645167/1‑100 (MQ=255)                                                    aaCGCAAATTTATGGTGGACAACAAGTAAGTCCTAATACAGAAAAGATAATTATTATTGCTCGAAAAATACCTATATTTATCTTCTTTGTATCTattatt                                       >  1:694778/1‑100 (MQ=255)                                                    aaCGCAAATTTATGGTGGACAACAAGTAAGTCCTAATACAGAAAAGATAATTATTATTGCTCGAAAAATACCTATATTTATCTTCTTTGTATCTattatt                                       >  1:2935318/1‑100 (MQ=255)                                                    aaCGCAAATTTATGGTGGACAACAAGTAAGTCCTAATACAGAAAAGATAATTATTATTGCTCGAAAAATACCTATATTTATCTTCTTTGTATCTattatt                                       >  1:305413/1‑100 (MQ=255)                                                    aaCGCAAATTTATGGTGGACAACAAGTAAGTCCTAATACAGAAAAGATAATTATTATTGCTCGAAAAATACCTATATTTATCTTCTTTGTATCTattatt                                       >  1:4127330/1‑100 (MQ=255)                                                    aaCGCAAATTTATGGTGGACAACAAGTAAGTCCTAATACAGAAAAGATAATTATTATTGCTCGAAAAATACCTATATTTATCTTCTTTGTATCTattatt                                       >  1:5795910/1‑100 (MQ=255)                                                    aaCGCAAATTTATGGTGGACAACAAGTAAGTCCTAATACAGAAAAGATAATTATTATTGCTCGAAAAATACCTATATTTATCTTCTTTGTATCTattatt                                       >  1:5623703/1‑100 (MQ=255)                                                    aaCGCAAATTTATGGTGGACAACAAGTAAGTCCTAATACAGAAAAGATAATTATTATTGCTCGAAAAATACCTATATTTATCTTCTTTGTATCTattatt                                       >  1:5320575/1‑100 (MQ=255)                                                    aaCGCAAATTTATGGTGGACAACAAGTAAGTCCTAATACAGAAAAGATAATTATTATTGCTCGAAAAATACCTATATTTATCTTCTTTGTATCTattatt                                       >  1:5025679/1‑100 (MQ=255)                                                     aCGCAAATTTATGGTGGACAACAAGTAAGTCCTAATACAGAAAAGATAATTATTATTGCTCGAAAAATACCTATATTTATCTTCTTTGTATCtattatcg                                      >  1:4752215/1‑98 (MQ=255)                                                     aCGCAAATTTATGGTGGACAACAAGTAAGTCCTAATACAGAAAAGATAATTATTATTGCTCGAAAAATACCTATATTTATCTTCTTTGTATCTATTATTg                                      >  1:2662911/1‑100 (MQ=255)                                                     aCGCAAATTTATGGTGGACAACAAGTAAGTCCTAATACAGAAAAGATAATTATTATTGCTCGAAAAATACCTATATTTATCTTCTTTGTATCTATTATTg                                      >  1:235617/1‑100 (MQ=255)                                                     aCGCAAATTTATGGTGGAAAACAAGTAAGTCCTAATACAGAAAAGATAATTATTATTGCTCGAAAAATACCTATATTTATCTTCTTTGTATCTATTATTg                                      >  1:846667/1‑100 (MQ=255)                                                      ggCAAATTTATGGTGGACAACAAGTAAGTCCTAATACAGAAAAGATAATTATTATTGCTCGAAAAATACCTATATTTATCTTCTTTGTATCTATTATTgg                                     >  1:4024859/2‑100 (MQ=255)                                                            cttATGGTGGACAACAAGTAAGTCCTAATACAGAAAAGATAATTATTATTGCTCGAAAAATACCTATATTTATCTTCTTTGTATCTATTATTGGTCCttt                               >  1:835165/2‑100 (MQ=255)                                                              tATGGTGGACAACAAGTAAGTCCTAATACAGAAAAGATAATTATTATTGCTCGAAAAATACCTATATTTATCTTCTTTGTATCTATTATTGGTCCTttat                             >  1:4607956/1‑100 (MQ=255)                                                              tATGGTGGACAACAAGTAAGTCCTAATACAGAAAAGATAATTATTATTGCTCGAAAAATACCTATATTTATCTTCTTTGTATCTATTATTGGTCCTttat                             >  1:2044675/1‑100 (MQ=255)                                                                      acaacaAGTAAGTCCTAATACAGAAAAGATAATTATTATTGCTCGAAAAATACCTATATTTATCTTCTTTGTATCTATTATTGGTCCTTTATTAgaagaa                     >  1:6061502/1‑100 (MQ=255)                                                                                 ttCCTAATACAGAAAAGATAATTATTATTGCTCGAAAAATACCTATATTTATCTTCTTTGTATCTATTATTGGTCCTTTATTAGAAGAATATGTATTCAg          >  1:258887/2‑100 (MQ=255)                                                                                         aCAGAAAAGATAATTATTATTGCTCGAAAAATACCTATATTTATCTTCTTTGTATCTATTATTGGTCCTTTATTAGAAGAATATGTATTCAGAAAAGTaa  >  1:5892588/1‑100 (MQ=255)                                                                                                  |                                                                                           GCTTGGATTTGCATTGGTAATGATTTATCAAATGGTAGTGAGCATTGTATTAACGCAAATTTATGGTGGACAACAAGTAAGTCCTAATACAGAAAAGCTAATTATTATTGCTCGAAAAATACCTATATTTATCTTCTTTGTATCTATTATTGGTCCTTTATTAGAAGAATATGTATTCAGAAAAGTAA  >  NC_007793/2142242‑2142429

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