Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_007793	2,503,835	T→C	D111D (GAT→GAC)	SAUSA300_RS12860 →	general stress protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_007793	2,503,835	0	T	C	97.2%	86.9 / ‑2.9	36	coding (333/423 nt)	SAUSA300_RS12860	general stress protein
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base C (33/2);  minor base . (1/0);  total (34/2)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.33e-01

TTGAGATGGAAGCGACAATCGAAATCGTTACAGAACAAGAAGTGATTGATTGGCTATGGGAAACACAAGACAAAAGCTTTTTCAGTTCAAAAGAAGATCCAGAGTTATGTGTTTTAAGAGTAGTTCCGCAATCCATTAAGCTAATGAATGATAAATCATTAGATACACCTATCAAAATCGATTTATAACACAAA  >  NC_007793/2503738‑2503931                                                                                                  |                                                                                                 ctGAGATGGAAGCGACAATCGAAATCGTTACAGAACAAGAAGTGATTGATTGGCTATGGGAAACACAAGACAAAAGCTTTTTCAGTTCAAAAGAGGAccc                                                                                                >  1:2366716/2‑100 (MQ=255)   gagaTGGAAGCGACAATCGAAATCGTTACAGAACAAGAAGTGATTGATTGGCTATGGGAAACACAAGACAAAAGCTTTTTCAGTTCAAAAGAAGACCCag                                                                                              >  1:4112861/1‑100 (MQ=255)     aaTGGAAGCGACAATCGAAATCGTTACAGAACAAGAAGTGATTGATTGGCTATGGGAAACACAAGACAAAAGCTTTTTCAGTTCAAAAGAAGACCCagag                                                                                            >  1:4379640/2‑100 (MQ=255)      ctGGAAGCGACAATCGAAATCGTTACAGAACAAGAAGTGATTGATTGGCTATGGGAAACACAAGACAAAAGCTTTTTCAGTTCAAAAGAGGACCCagatc                                                                                           >  1:1961222/2‑98 (MQ=255)       gtgtggCGACAATCGAAATCGTTACAGAACAAGAAGTGATTGATTGGCTATGGGAAACACAAGACAAAAGCTTTTTCAGTTCAAAAGAAGA‑CCagagatc                                                                                         >  1:3150995/6‑97 (MQ=255)              gACAATCGAAATCGTTACAGAACAAGAAGTGATTGATTGGCTATGGGAAACACAAGACAAAAGCTTTTTCAGTTCAAAAGAAGACCCagatcggaagagc                                                                                   >  1:3373565/1‑90 (MQ=255)              gACAATCGAAATCGTTACAGAACAAGAAGTGATTGATTGGCTATGGGAAACACAAGACAAAAGCTTTTTCAGTTCAAAAGAAGACCCAGAGTTATgtgtg                                                                                   <  1:2807638/100‑2 (MQ=255)                 aaTCGAAATCGTTACAGAACAAGAAGTGATTGATTGGCTATGGGAAACACAAGACAAAAGCTTTTTCAGTTCAAAAGAAGACCCAGAGTTATGTGTTTTa                                                                                <  1:1073657/100‑1 (MQ=255)                  aTCGAAATCGTTACAGAACAAGAAGTGATTGATTGGCTATGGGAAACACAAGACAAAAGCTTTTTCAGTTCAAAAGAAGACCCAGAGTTATGTGTTTTaa                                                                               >  1:2007741/1‑100 (MQ=255)                             tACAGAACAAGAAGTGATTGATTGGCTATGGGAAACACAAGACAAAAGCTTTTTCAGTTCAAAAGAAGACCCAGAGTTATGTGTTTTAAGAGTAGTTCCg                                                                    >  1:2252211/1‑100 (MQ=255)                             tACAGAACAAGAAGTGATTGATTGGCTATGGGAAACACAAGACAAAAGCTTTTTCAATTCAAAAGAAGACCCAGAGTTATGTGTTTTAAGAGTAGTTCCg                                                                    >  1:5576361/1‑100 (MQ=255)                              ccAGAACAAGAAGTGATTGATTGGCTATGGGAAACACAAGACAAAAGCTTTTTCAGTTCAAAAGAAGACCCAGAGTTATGTGTTTTAAGAGTAGTTCCGc                                                                   >  1:1334435/2‑100 (MQ=255)                               ctgcaCAAGAAGTGATTGATTGGCTATGGGAAACACAAGACAAAAGCTTTTTCAGTTCAAAAGAAGACCCAGAGTTATGTGTTTTAAGAGTAGTTCCGCa                                                                  >  1:5779057/5‑100 (MQ=255)                               cAGAACAAGAAGTGATTGATTGGCTATGGGAAACACAAGACAAAAGCTTTTTCAGTTCAAAAGAAGACCCAGAGTTATGTGTTTTAAGAGTAGTTCCGCa                                                                  >  1:3196849/1‑100 (MQ=255)                               cAGAACAAGAAGAGATTGATTGGCTATGGGAAACACAAGACAAAAGCTTTTTCAGTTCAAAAGAAGACCCAGAGTTATGTGTTTTAAGAGTAGTTCCGCa                                                                  >  1:2132512/1‑100 (MQ=255)                                   aCAAGAAGTGATTGATTGGCTATGGGAAACACAAGACAAAAGCTTTTTCAGTTCAAAAGAAGACCCAGAGTTATGTGTTTTAAGAGTAGTTCCGCAATcc                                                              >  1:89878/1‑100 (MQ=255)                                   aCAAGAAGTGATTGATTGGCTATGGGAAACACAAGACAAAAGCTTTTTCAGTTCAAAAGAAGACCCAGAGTTATGTGTTTTAAGAGTAGTTCCGCAATcc                                                              >  1:950196/1‑100 (MQ=255)                                   aCAAGAAGTGATTGATTGGCTATGGGAAACACAAGACAAAAGCTTTTTCAGTTCAAAAGAAGACCCAGAGTTATGTGTTTTAAGAGTAGTTCCGCAATcc                                                              >  1:1456105/1‑100 (MQ=255)                                            gcttgattgGCTATGGGAAACACAAGACAAAAGCTTTTTCAGTTCAAAAGAAGACCCAGAGTTATGTGTTTTAAGAGTAGTTCCGCAATCCATTAAGCTa                                                     >  1:1070783/3‑100 (MQ=255)                                               tgattgGCTATGGGAAACACAAGACAAAAGCTTTTTCAGTTCAAAAGAAGACCCAGAGTTATGTGTTTTAAGAGTAGTTCCGCAATCCATTAAGCTaatg                                                  >  1:1006844/1‑100 (MQ=255)                                               tgattgGCTATGGGAAACACAAGACAAAAGCTTTTTCAGTTCAAAAGAAGACCCAGAGTTATGTGTTTTAAGAGTAGTTCCGCAATCCATTAAGCTaatg                                                  >  1:5590040/1‑100 (MQ=255)                                               tgattgGCTATGGGAAACACAAGACAAAAGCTTTTTCAGTTCAAAAGAAGACCCAGAGTTATGTGTTTTAAGAGTAGTTCCGCAATCCATTAAGCTaatg                                                  >  1:5366142/1‑100 (MQ=255)                                               tgattgGCTATGGGAAACACAAGACAAAAGCTTTTTCAGTTCAAAAGAAGACCCAGAGTTATGTGTTTTAAGAGTAGTTCCGCAATCCATTAAGCTaatg                                                  >  1:5414840/1‑100 (MQ=255)                                                 cttgGCTATGGGAAACACAAGACAAAAGCTTTTTCAGTTCAAAAGAAGACCCAGAGTTATGTGTTTTAAGAGTAGTTCCGCAATCCATTAAGCTaatgaa                                                >  1:2646601/2‑100 (MQ=255)                                                  ctgGCTATGGGAAACACAAGACAAAAGCTTTTTCAGTTCAAAAGAAGACCCAGAGTTATGTGTTTTAAGAGTAGTTCCGCAATCCATTAAGCTaatgaat                                               >  1:166419/2‑100 (MQ=255)                                                   tgGCTATGGGAAACACAAGACAAAAGCTTTTTCAGTTCGAAAGAAGACCCAGAGTTATGTGTTTTAAGAGTAGTTCCGCAATCCATTAAGCTaatgaatg                                              >  1:1206307/1‑100 (MQ=255)                                                   tgGCTATGGGAAACACAAGACAAAAGCTTTTTCAGTTCAAAAGAAGACCCAGAGTTATGTGTTTTAAGAGTAGTTCCGCAATCCATTAAGCTaatgaatg                                              >  1:1790428/1‑100 (MQ=255)                                                   tgGCTATGGGAAACACAAGACAAAAGCTTTTTCAGTTCAAAAGAAGACCCAGAGTTATGTGTTTTAAGAGTAGTTCCGCAATCCATTAAGCTaatgaatg                                              >  1:1734871/1‑100 (MQ=255)                                                   tgGCTATGGGAAACACAAGACAAAAGCTTTTTCAGTTCAAAAGAAGACCCAGAGTTATGTGTTTTAAGAGTAGTTCCGCAATCCATTAAGCTaatgaatg                                              >  1:5605149/1‑100 (MQ=255)                                                    ggCTATGGGAAACACAAGACAAAAGCTTTTTCAGTTCAAAAGAAGACCCAGAGTTATGTGTTTTAAGAGTAGTTCCGCAATCCATTAAGCTAatgaatga                                             >  1:2097211/1‑100 (MQ=255)                                                    ggCTATGGGAAACACAAGACAAAAGCTTTTTCAGTTCAAAAGAAGACCCAGAGTTATGTGTTTTAAGAGTAGTTCCGCAATCCATTAAGCTAatgaatga                                             >  1:133361/1‑100 (MQ=255)                                                     gCTATGGGAAACACAAGACAAAAGCTTTTTCAGTTCAAAAGAAGACCCAGAGTTATGTGTTTTAAGAGTAGTTCCGCAATCCATTAAGCTAATGAATGAt                                            >  1:4772144/1‑100 (MQ=255)                                                       caTGGGAAACACAAGACAAAAGCTTTTTCAGTTCAAAAGAAGACCCAGAGTTATGTGTTTTAAGAGTAGTTCCGCAATCCATTAAGCTAATGAATGATaa                                          >  1:4532984/2‑100 (MQ=255)                                                        aTGGGAAACACAAGACAAAAGCTTTTTCAGTTCAAAAGAAGACCCAGAGTTATGTGTTTTAAGAGTAGTTCCGCAATCCATTAAGCTAATGAATGATaaa                                         >  1:4327748/1‑100 (MQ=255)                                                          gggAAACACAAGACAAAAGCTTTTTCAGTTCAAAAGAAGACCCAGAGTTATGTGTTTTAAGAGTAGTTCCGCAATCCATTAAGCTAATGAATGATAAATc                                       >  1:5126515/1‑100 (MQ=255)                                                               acacAAGACAAAAGCTTTTTCAGTTCAAAAGAAGACCCAGCGTTATGTGTTTTAAGAGTAGTTCCGCAATCCATTAAGCTAATGAATGATAAATCATTAg                                  >  1:5437526/1‑100 (MQ=255)                                                                                               agaCCCAGAGTTATGTGTTTTAAGAGAAGTTCCGCAATCCATTAAGCTAATGAATGATAAATCATTAGATACACCTATCAAAATCGTTTTATAACACaaa  >  1:5208063/1‑100 (MQ=255)                                                                                                  |                                                                                                 TTGAGATGGAAGCGACAATCGAAATCGTTACAGAACAAGAAGTGATTGATTGGCTATGGGAAACACAAGACAAAAGCTTTTTCAGTTCAAAAGAAGATCCAGAGTTATGTGTTTTAAGAGTAGTTCCGCAATCCATTAAGCTAATGAATGATAAATCATTAGATACACCTATCAAAATCGATTTATAACACAAA  >  NC_007793/2503738‑2503931

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