Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_007793	700,471	A→T	E217D (GAA→GAT)	SAUSA300_RS03360 →	glycosyl transferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_007793	700,471	0	A	T	98.7%	204.9 / ‑5.6	76	E217D (GAA→GAT)	SAUSA300_RS03360	glycosyl transferase
Reads supporting (aligned to +/- strand):  ref base A (1/0);  major base T (71/3);  minor base C (1/0);  total (73/3)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.66e-01

TGTCAAGGCGATGCAAAAAGCTACGGATATACATTTAGTTTCACAGATGATTGTCGGAGATAACGATATAGATACACATGCTACAAGTAACGAAGAAGATTTTAATAGATATATCACAGAAATTATGAAAATAAGACAACGAGTCATGGAAATGTTACTATTACCTGAACAAAGGCTATTATATAGTG  >  NC_007793/700378‑700565                                                                                              |                                                                                               tGTCAAGGCGATGCAAAAAGCTACGGTTACACATTTAGTTTCACAGATGATTGTCGGAGATAACGATATAGATACACATGCCACAAGTAACAATGAAGAt                                                                                          >  1:2879821/1‑100 (MQ=255) tGTCAAGGCGATGCAAAAAGCTACGGATATACATTTAGTTTCACAGATGATTGTCGGAGATAACGATATAGATACACATGCTACAAGTAACGGTGAAGAt                                                                                          >  1:3563812/1‑100 (MQ=255) tGTCAAGGCGATGCAAAAAGCTACGGATATACATTTAGTTTCACAGATGATTGTCGGAGATAACGATATAGATACACATGCTACAAGTAACGATGAAGAt                                                                                          >  1:6536211/1‑100 (MQ=255)   acAAGGCGATGCAAAAAGCTACGGATATACATTTAGTTTCACAGATGATTGTCGGAGATAACGATATAGATACACATGCTACAAGTAATGATGAAGAtcg                                                                                        >  1:2859125/2‑98 (MQ=255)   acAAGGCGATGCAAAAAGCTACGGATATACATTTAGTTTCACAGATGATTGTCGGAGATAACGATATAGATACACATGCTACAAGTAACGCTGAAGAtcg                                                                                        >  1:4647412/2‑98 (MQ=255)    cAAGGCGATGCAAAAAGCTACGGGGATACATTTAGTTTCACAGATGATTGTCGGAGATAACGATATAGATACACATGCTACAAGTAATGATGAAGAtcgg                                                                                       >  1:3588862/1‑97 (MQ=255)    cAAGGCGATGCAAAAAGCTACGGATATACATTTAGTTTCACAGATGATTGTCGGAGATAACGATATAGATACACATGCTACAAGTGACGATGAAGAtcgg                                                                                       >  1:4786206/1‑97 (MQ=255)    cAAGGCGATGCAAAAAGCTACGGATATACATTTAGTTTCACAGATGATTGTCGGAGATAACGATATAGATACACATGCTACAAGTAACGATGAAGAtttt                                                                                       >  1:6605081/1‑100 (MQ=255)    cAAGGCGATGCAAAAAGCTACGGATATACATTTAGTTTCACAGATGATTGTCGGAGATAACGATATAGATACACATGCTACAAGTAACGATGAAGAtcgg                                                                                       >  1:2167017/1‑97 (MQ=255)     aaGGCGATGCAAAAAGCTACGGATATACATTTAGTGTCACAGATGATTGTCGGAGATAACGATATAGATACACATGCTACAAGTAACGATGAAGATTTTa                                                                                      >  1:5868663/1‑100 (MQ=255)       ggCGATGCAAAAAGCTACGGATATACATTTAGTTTCACAGATGATTGTCGGAGATAACGATATAGATACACATGCTACAAGTAACGATGAAGCTTTTAAt                                                                                    >  1:6404277/1‑100 (MQ=255)         cGATGCAAAAAGCTACGGATATACATTTAGTTTCACAGATGATTGTCGGAGATAACGATATAGATACACATGCTACAAGTAACGGCGAAGAtcggaagag                                                                                  >  1:3222995/1‑92 (MQ=255)          taTGCAAAAAGCTACGGATATACATTTAGTTTCACAGATGATTGTCGGAGATAACGATATAGATACACATGCTACAAGTAAAGATGAAGAtcggaagagc                                                                                 >  1:5664898/2‑91 (MQ=255)           ctGCAAAAAGCTACGGATATACATTTAGTTTCACAGATGATTGTCGGAGATAACGATATAGATACACATGCTACAAGTAACGATGAAGATTTTAATAGat                                                                                >  1:4609450/2‑100 (MQ=255)           aTGCAAAAAGCTACGGATATCCATTTAGTTTCACAGATGATTGTCGGAGATAACGATATAGATACACATGCTACAAGTAACAATGAAGAtcggaagagca                                                                                >  1:5839571/1‑90 (MQ=255)           aTGCAAAAAGCTACGGATATACATTTAGTTTCACAGATGATTGTCGGAGATAACGATATAGATACACATGCTACAAGTAACGGTGAAGAtcggaagagca                                                                                >  1:2162842/1‑90 (MQ=255)           aTGCAAAAAGCTACGGATATACATTTAGTTTCACAGATGATTGTCGGAGATAACGATATAGATACACATGCTACAAGTAACGAAGAAGAtcggaagagca                                                                                >  1:6192125/1‑90 (MQ=255)               aaaaaGCTACGGATATACATTTAGTTTCACAGATGATTGTCGGAGATAACGATATAGATACACATGCTACAAGTAACGATGAAGATTTTAATAGatatat                                                                            >  1:4958626/1‑100 (MQ=255)                aaaaGCTACGGATATACATTTAGTTTCACAGATGATTGTCGGAGATAACGATATAGATACACATGCTACAAGTAACGATGAAGATTTTAATAGATATATc                                                                           >  1:2215775/1‑100 (MQ=255)                  taGCTACGGATATACATTTAGTTTCACAGATGATTGTCGGAGATAACGATATAGATACACATGCTACAAGTAACGATGAAGATTTTAATAGATATATcac                                                                         >  1:3640029/2‑100 (MQ=255)                  aaGCTACGGATATACATTTAGTTTCACAGATGATTGTCGGAGATAACGATATAGATACACATGCTACAAGTAACGATGAAGATTTTAATAGATATATcac                                                                         <  1:3926346/100‑1 (MQ=255)                  aaGCTACGGATATACAATTAGTTTCACAGATGATTGTCGGAGATAACGATATAGATTCACATGTTACAAGTAACGATGAAGATTTTAATAGATATAGcac                                                                         <  1:6566398/100‑1 (MQ=255)                   aGCTACGGATATACATTTAGTTTCACAGATGATTGTCGGAGATAACGATATAGATACACATGCTACAAGTAACGATGAAGATTTTAATAGATATATcaca                                                                        >  1:1780816/1‑100 (MQ=255)                    gCTACGGATATACATTTAGTTTCACAGATGATTGTCGGAGATAACGATATAGATACACATGCTACAAGTAACGATGAAGATTTTAATAGATATATCACAg                                                                       >  1:1633000/1‑100 (MQ=255)                    acTACGGATATACATTTAGTTTCACAGATGATTGTCGGAGATAACGATATAGATACACATGCTACAAGTAACGATGAAGATTTTAATAGATATATCACAg                                                                       >  1:5182502/2‑100 (MQ=255)                       aCGGATATACATTTAGTTTCACAGATGATTGTCGGAGATAACGATATAGATACACATGCTACGAGTAACGATGAAGATTTTAATAGATATATCACAGaaa                                                                    >  1:1330360/1‑100 (MQ=255)                       aCGGATATACATTTAGTTTCACAGATGATTGTCGGAGATAACGATATAGATACACATGCTACAAGTAACGATGAAGATTTTAATAGATATATCACAGaaa                                                                    >  1:2218208/1‑100 (MQ=255)                            tataCATTTAGTTTCACAGATGATTGTCGGAGATAACGATATAGATACACATGCTACAAGTAACGATGAAGATTTTAATAGATATATCACAGAAATTATg                                                               >  1:3882762/1‑100 (MQ=255)                             ataCATTTAGTTTCACAGATGATTGTCGGAGATAACGATATAGATACACATGCTACAAGTAACGATGAAGATTTTAATAGATATATCACAGAAATTATGa                                                              >  1:5785715/1‑100 (MQ=255)                             ataCATTTAGTTTCACAGATGATTGTCGGAGATAACGATATAGATACACATGCTACAAGTAACGATGAAGATTTTAATAGATATATCACAGAAATTATGa                                                              >  1:3813482/1‑100 (MQ=255)                               aCATTTAGTTTCACAGATGATTGTCGGAGATAACGATATAGATACACATGCTACAAGTAACGATGAAGATTTTAATAGATATATCACAGAAATTATGaaa                                                            <  1:4239198/100‑1 (MQ=255)                                    tAGTTTCACAGATGATTGTCGGAGATAACGATATAGATACACATGCTACAAGTAACGATGAAGATTTTAATAGATATATCACAGAAATTATGAAAATAAg                                                       >  1:2597514/1‑100 (MQ=255)                                        ttCACAGATGATTGTCGGAGATAACGATATAGATACACATGCTACAAGTAACGATGAAGATTTTAATAGATATATCACAGAAATTATGAAAATAAGACaa                                                   >  1:882975/1‑100 (MQ=255)                                          cacaGATGATTGTCGGAGATAACGATATAGATACACATGCTACAAGTAACGATGAAGATTTTAATAGATATATCACAGAAATTATGAAAATAAGACAACg                                                 >  1:84888/1‑100 (MQ=255)                                          cacaGATGATTGTCGGAGATAACGATATAGATACACATGCTACAAGTAACGATGAAGATTTTAATAGATATATCACAGAAATTATGAAAATAAGACAACg                                                 >  1:455749/1‑100 (MQ=255)                                          cacaGATGATTGTCGGAGATAACGATATAGATACACATGCTACAAGTAACGATGAAGATTTTAATAGATATATCACAGAAATTATAAAAATGAGACAACg                                                 >  1:4935379/1‑100 (MQ=255)                                           acaGATGATTGTCGGAGATAACGATATAGATACACATGCTACAAGTAACGATGAAGATTTTAATAGATATATCACAGAAATTATGAAAATAAGACAACGa                                                >  1:1187603/1‑100 (MQ=255)                                            caGATGATGGTCGGAGATAACGATATAGATACACATGCTACAAGTAACGATGAAGATTTTAATAGATATATCACAGAAATTATGAAAATAAGACAACGAg                                               >  1:6055246/1‑100 (MQ=255)                                             tgatgatTGTCGGAGATAACGATATAGATACACATGCTACAAGTAACGATGAAGATTTTAATAGATATATCACAGAAATTATGAAAATAAGACAACGAGt                                              >  1:578182/2‑100 (MQ=255)                                             aGATGATTGTCGGAGATAACGATATAGATACACATGCTACAAGTAACGATGAAGATTTTAATAGATATATCACAGAAATTATGAAAATAAGACAACGAGt                                              >  1:1164820/1‑100 (MQ=255)                                               ctgatTGTCGGAGATAACGATATAGATACACATGCTACAAGTAACGATGAAGATTTTAATAGATATATCACAGAAATTATGAAAATAAGACAACGAGTCa                                            >  1:5153289/2‑100 (MQ=255)                                                 gatTGTCGGAGATAACGATATAGATACACATGCTACAAGTAACGATGAAGATTTTAATAGATATATCACAGAAATTATGAAAATAAGACAACGAGTCATg                                          >  1:5427503/1‑100 (MQ=255)                                                 gatTGTCGGAGATAACGATATAGATACACATGCTACAAGTAACGATGAAGATTTTAATAGATATATCACAGAAATTATGAAAATAAGACAACGAGTCATg                                          >  1:1416783/1‑100 (MQ=255)                                                    tGTCGGAGATAACGATATAGATACACATGCTACAAGTAACGATGAAGATTTTAATAGATATATCACAGAAATTATGAAAATAAGACAACGAgatcggaag                                       >  1:1131893/1‑92 (MQ=255)                                                           gaaaccGATATAGATACACATGCTACAAGTAACGATGAAGATTTTAATAGATATATCACAGAAATTATGAAAATAAGACAACGAGTCATGGAAATGTTAc                                >  1:2308462/6‑100 (MQ=255)                                                           gaTAACGATATAGATACACATGCTACAAGTAACGATGAAGATTTTAATAGATATATCACAGAAATTATGAAAATAAGACAACGAGTCATGGAAATGTTAc                                >  1:1081318/1‑100 (MQ=255)                                                             caaCGATATAGATACACATGCTACAAGTAACGATGAAGATTTTAATAGATATATCACAGAAATTATGAAAATAAGACAACGAGTCATGGAAATGTTACTa                              >  1:6486509/2‑100 (MQ=255)                                                              aaCGATATAGATACACATGCTACAAGTAACGATGAAGATTTTAATAGATATATCACAGAAATTATGAAAATAAGACAACGAGTCATGGAAATGTTACTAt                             >  1:3657474/1‑100 (MQ=255)                                                              aaCGATATAGATACACATGCTACAAGTAAAAATGAAGATTTTAATAGATATATCACAGAAATTATGAAAATAAGACAACGAGTCATGGAAATGTTACTAt                             >  1:3358881/1‑100 (MQ=255)                                                                   tataGATACACATGCTACAAGTAACGATGAAGATTTTAATAGATATATCACAGAAATTATGAAAATAAGACAACGAGTCATGGAAATGTTACTATTACCt                        >  1:552942/1‑100 (MQ=255)                                                                    ataGATACACATGCTACAAGTAACGATGAAGATTTTAATAGATATATCACAGAAATTATGAAAATAAGACAACGAGTCATGGAAATGTTACTATTACCTg                       >  1:2961265/1‑100 (MQ=255)                                                                      aGATACACATGCTACAAGTAACGATGAAGATTTTAATAGATATATCACAGAAATTATGAAAATAAGACAACGAGTCATGGAAATGTTACTATTACCTGaa                     >  1:6172550/1‑100 (MQ=255)                                                                         tACACATGCTACAAGTAACGATGAAGATTTTAATAGATATATCACAGAAATTATGAAAATAAGACAACGAGTCATGGAAATGTTACTATTACCTGAACaa                  >  1:2233987/1‑100 (MQ=255)                                                                          acacATGCTACAAGTAACGATGAAGATTTTAATAGATATATCACAGAAATTATGAAAATAAGACAACGAGTCATGGAAATGTTACTATTACCTGAACaaa                 >  1:1687932/1‑100 (MQ=255)                                                                          acacATGCTACAAGTAACGATGAAGATTTTAATAGATATATCACAGAAATTATGAAAATAAGACAACGAGTCATGGAAATGTTACTATTACCTGAACaaa                 >  1:3212489/1‑100 (MQ=255)                                                                           tacaTGCTACAAGTAACGATGAAGATTTTAATAGATATATCACAGAAATTATGAAAATAAGACAACGAGTCATGGAAATGTTACTATTACCTGGACAAAg                >  1:5914261/2‑100 (MQ=255)                                                                           cgcaTGCTACAAGTAACGATGAAGATTTTAATAGATATATCACAGAAATTATGAAAATAAGACAACGAGTCATGGAAATGTTACTATTACCTGAACAAAg                >  1:5127421/3‑100 (MQ=255)                                                                           cacaTGCTACAAGTAACGATGAAGATTTTAATAGATATATCACAGAAATTATGAAAATAAGACAACGAGTCATGGAAATGTTACTATTACCTGAACAAAg                >  1:5522343/1‑100 (MQ=255)                                                                            acaTGCTACAAGTAACGATGAAGATTTTAATAGATATATCACAGAAATTATGAAAATAAGACAACGAGTCATGGAAATGTTACTATTACCTGAACAAAgg               >  1:2120041/1‑100 (MQ=255)                                                                              gtGCTACAAGTAACGATGAAGATTTTAATAGATATATCACAGAAATTATGAAAATAAGACAACGAGTCATGGAAATGTTACTATTACCTGAACAAAGGCt             >  1:6430521/2‑100 (MQ=255)                                                                              aTGCTACAAGTAACGATGAAGATTTTAATAGATATATCACAGAAATTATGAAAATAAGACAACGAGTCATGGAAATGTTACTATTACCTGAACAAAGGCt             >  1:1481235/1‑100 (MQ=255)                                                                               tGCTACAAGTAACGATGAAGATTTTAATAGATATATCACAGAAATTATGAAAATAAGACAACGAGTCATGGAAATGTTACTATTACCTGAACAAAGGCta            >  1:6603574/1‑100 (MQ=255)                                                                               tGCTACAAGTAACGATGAAGATTTTAATAGATATATCACAGAAATTATGAAAATAAGACAACGAGTCATGGAAATGTTACTATTACCTGAACAAAGGCta            >  1:1511802/1‑100 (MQ=255)                                                                                gCTACAAGTAACGATGAAGATTTTAATAGATATATCACAGAAATTATGAAAATAAGACAACGAGTCATGGAAATGTTACTATTACCTGAACAAAGGCtat           >  1:919356/1‑100 (MQ=255)                                                                                gCTACAAGTAACGATGAAGATTTTAATAGATATATCACAGAAATTATGAAAATAAGACAACGAGTCATGGAAATGTTACTATTACCTGAACAAAGGCtat           >  1:4748105/1‑100 (MQ=255)                                                                                gCTACAAGTAACGATGAAGATTTTAATAGATATATCACAGAAATTATGAAAATAAGACAACGAGTCATGGAAATGTTACTATTACCTGAACAAAGGCtat           >  1:420508/1‑100 (MQ=255)                                                                                gCTACAAGTAACGATGAAGATTTTAATAGATATATCACAGAAATTATGAAAATAAGACAACGAGTCATGGAAATGTTACTATTACCTGAACAAAGGCtat           >  1:1904672/1‑100 (MQ=255)                                                                                   aCAAGTAACGATGAAGCTTTTAATAGCTATATCACAGAAATTATGAAAATAAGACAACGAGTCATGGAAATGTTACTATTACCTGAACAAAGGCtattat        >  1:2854401/1‑100 (MQ=255)                                                                                   aCAAGTAACGATGAAGATTTTAATAGATATATCACAGAAATTATGAAAATAAGACAACGAGTCATGGAAATGTTACTATTACCTGAACAAAGGCtattat        >  1:6463770/1‑100 (MQ=255)                                                                                     taGTAACGATGAAGATTTTAATAGATATATCACAGAAATTATGAAAATAAGACAACGAGTCATGGAAATGTTACTATTACCTGAACAAAGGCTATtatat      >  1:611922/2‑100 (MQ=255)                                                                                      aGTAACGATGAAGATTTTAATAGATATATCACAGAAATTATGAAAATAAGACAACGAGTCATGGAAATGTTACTATTACCTGAACCAAGGCTATtatata     >  1:4933014/1‑100 (MQ=255)                                                                                      aGTAACGATGAAGATTTTAATAGATATATCACAGAAATTATGAAAATAAGACAACGAGTCATGGAAATGTTACTATTACCTGAACAAAGGcaagatcgga     >  1:3307212/1‑91 (MQ=255)                                                                                      aGTAACGATGAAGATTTTAATAGATATATCACAGAAATTATGAAAATAAGACAACGAGTCATGGAAATGTTACTATTACCTGAACAAAGGCTATtatata     >  1:1318107/1‑100 (MQ=255)                                                                                        tAACGATGAAGATTTTAATAGATATATCACAGAACTTTTGAAAATAAGACAACGAGTCATGGAAATGTTACTATTACCTGAACAAAGGCTATTATATAGt   >  1:4675762/1‑100 (MQ=255)                                                                                        caaCGATGAAGATTTTAATAGATATATCACAGAAATTATGAAAATAAGACAACGAGTCATGGAAATGTTACTATTACCTGAACAAAGGCTATTATATAGt   >  1:31626/2‑100 (MQ=255)                                                                                         aaCGATGAAGATTTTAATAGATATATCACAGAAATTATGAAAATAAGACAACGAGTCATGGAAATGTTACTATTACCTGAACAAAGGCTATTATATAGTg  >  1:132069/1‑100 (MQ=255)                                                                                              |                                                                                               TGTCAAGGCGATGCAAAAAGCTACGGATATACATTTAGTTTCACAGATGATTGTCGGAGATAACGATATAGATACACATGCTACAAGTAACGAAGAAGATTTTAATAGATATATCACAGAAATTATGAAAATAAGACAACGAGTCATGGAAATGTTACTATTACCTGAACAAAGGCTATTATATAGTG  >  NC_007793/700378‑700565

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