Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_007793	1,923,107	A→T	I33F (ATT→TTT)	SAUSA300_RS09510 →	calcium‑binding protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_007793	1,923,107	0	A	T	84.6%	15.3 / 0.7	13	I33F (ATT→TTT)	SAUSA300_RS09510	calcium‑binding protein
Reads supporting (aligned to +/- strand):  ref base A (2/0);  new base T (11/0);  total (13/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 = 3.39e-02

TGAATACAGAGAACAACAAGAATCAAAACCAATCTGTTAAAAATTCTGAAAGACGCGGCATGTTAAAAGGATGCGGCGGTTGCCTTATTTCTTTTATTTTATTAATAATCTTATTATCAGCCTGTTCAATGATGTTTAGTAATAATGACAATTCCACTAATAATCAATCATCAAAAACGCAATTAACTCAAA  >  NC_007793/1923012‑1923203                                                                                                |                                                                                                 tGAATACAGAGAACAACAAGAATCAAAACCAATCTGTTAAAAATTCTGAAAGACGCGGCATGTTAAAAGGATGCGGCGGTTTTCtttttttttttttttt                                                                                              >  1:222211/1‑100 (MQ=255) tGAATACAGAGAACAACAAGAATCAAAACCAATCTGTTAAAAATTCTGAAAGACGCGGCATGTTAAAAGGATGCGGCGGTTTGCtttttttttttttttt                                                                                              >  1:3252745/1‑100 (MQ=255) tGAATACAGAGAACAACAAGAATCAAAACCAATCTGTTAAAAATTCTGAAAGACGCGGCATGTTAAAAGGATGCGGCGGTTTCCtttttttttttttttt                                                                                              >  1:3706758/1‑100 (MQ=255) tGAATACAGAGAACAACAAGAATCAAAACCAATCTGTTAAAAATTCTGAAAGACGCGGCATGTTAAAAGGATGCGGCGGTTGCCCTCTtttttttttttt                                                                                              >  1:1996741/1‑100 (MQ=255) agAATACAGAGAACAACAAGGATCAAAACCAATCTGTTAAAAATTCTGAAAGACGCGGCATGTTAAAAGGATGCGGCGGTTTTCtttttttttttttttt                                                                                              >  1:2744919/2‑100 (MQ=255)  gAATACAGAGAACAACAAGAATCAAAACCAATCTGTTAAAAATTCTGAAAGACGCGGCATGTTAAAAGGATGCGGCGGTTTCCTTTTTTTttttatttta                                                                                             >  1:2360219/1‑100 (MQ=255)   aaTACAGAGAACAACAAGAATCAAAACCAATCTGTTAAAAATTCTGAAAGACGCGGCATGTTAAAAGGATGCGGCGGTTTCCCTtttttttttttttttt                                                                                            >  1:2862953/1‑98 (MQ=255)             acaacaAGAATCAAAACCAATCTGTTAAAAATTCTGAAAGACGCGGCATGTTAAAAGGATGCGGCGGTTGCCTTTTTTCTTTTTTTTTTTTATTAATAtt                                                                                  >  1:5323751/1‑100 (MQ=255)                acaaGAATCAAAACCAATCTGTTAAAAATTCTGAAAGACGCGGCATGTTAAAAGGATGCGGCGGTTGCCCTTTTTCTTTTTTTTTTTTAATAAACttttt                                                                               >  1:4955220/1‑100 (MQ=255)                        cAAAACCAATCTGTTAAAAATTCTGAAAGACGCGGCATGTTAAAAGGATGTGGCGGTTGTTTTTTTTTTTTTATTTTTTTAATAAACTTAATATCAAccc                                                                       >  1:3447400/1‑99 (MQ=255)                                    gTTAAAAATTCTGAAAGACGCGGCATGTTAAAAGGATGCGGCGGTTTCCTTTTTTTTTTTTTTTTTTTTATAAACTTTTTATTAACCTGTTTAATGATGt                                                           >  1:1305519/1‑100 (MQ=255)                                          aaTTCTGAAAGACGCGGCATGTTAAAAGGATGCGGCGGTTGCCTTTTTTTTTTTTTTTTTTTAATAAACTTTTTATTAACCTGTTTCATGATGTTTTGta                                                     >  1:2050925/1‑100 (MQ=255)                                             ccTGAAAGACGCGGCATGTTAAAAGGATGCGGCGGTTGCCCTTTTTTTTTTTTTTTTTTAATAAACTTATTTTCAACCCGTTCAAAGATGTTTATtaaaa                                                  >  1:5492092/2‑98 (MQ=255)                                                   aGACGCGGCATGTTAAAAGGATGCGGCGGTTGCCTTTTTTTTTTTTTTTTTTTAATAAATTTTTTTTTAACCTGTTCAATTATGTTTTGTaaaaattaaa                                            >  1:4514933/1‑96 (MQ=255)                                                    atggtggCATGTTAAAAGGATGCGGCGGTTTCCCTATTTTTTTTTTTTTATTTATAAACTTATTTTCAACCCGTTTAAAGATGTTTTTTAAAAATGACaa                                           >  1:548009/6‑100 (MQ=255)                                                      aacgGCATGTTAAAAGGATGCGGCGGTTGTCCTTTTTTTTTTTTTTTTTTAATAATCTTTTTTTTAACCCGTTTTATGATGTTTTATAATAAAAACAAtt                                         >  1:5593093/3‑100 (MQ=255)                                                        cgGCATGTTAAAAGGATGCGGCGGTTTCCTTTTTTTTTTTTTTTTTTTAATAATTTTTTTTTTAACCCGTTTAATGATGTTTTGTAATAATGACAATTcc                                       >  1:2533490/1‑100 (MQ=255)                                                        cgGCATGTTAAAAGGATGCGGCGGTTGCCTTTTTTTTTTTTTTTTTTTTATAAACTTATTATTAACCCGTTCTATGTTGTTTTGTAATAAATACAAtttc                                       >  1:948153/1‑98 (MQ=255)                                                         tgCATGTTAAAAGGATGCGGCGGTTGCCTTATTTTTTTTTTTTTTTTAATAAACTTATTTTTAACCTGTTCAATGATGTTTTGTAATAATGACAATTCCa                                      >  1:5125824/2‑100 (MQ=255)                                                         ggCATGTTAAAAGGATGCGGCGGTTGCCCTTTTTTTTTTTTTTTTTTTATAATCTTTTTATCACCCCGTTCTATTATGTTTTTTAATAATGACaaattca                                      >  1:2753908/1‑95 (MQ=255)                                                                                ctGCCTTATTTCTTTTATTTTATTAATAATCTTATTATCAGCCTGTTCAATGATGTTTAGTAATAATGACAATTCCACTAATAATCAATCATCAAAAACg               >  1:916465/2‑100 (MQ=255)                                                                                   ccTTATTTCTTTTATTTTATTAATAATCTTATTATCAGCCTGTTCAATGATGTTTAGTAATAATGACAATTCCACTAATAATCAATCATCAAAAACGCaa            >  1:3223076/1‑100 (MQ=255)                                                                                             tttattttatTAACAATCTTATTATCAGCCTGTTCAATGATGTTTAGTAATAATGACAATTCCACTAATAATCAATCATCAAAAACGCAATTAACTCaaa  >  1:13162/1‑100 (MQ=255)                                                                                                |                                                                                                 TGAATACAGAGAACAACAAGAATCAAAACCAATCTGTTAAAAATTCTGAAAGACGCGGCATGTTAAAAGGATGCGGCGGTTGCCTTATTTCTTTTATTTTATTAATAATCTTATTATCAGCCTGTTCAATGATGTTTAGTAATAATGACAATTCCACTAATAATCAATCATCAAAAACGCAATTAACTCAAA  >  NC_007793/1923012‑1923203

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