Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_007793	1,578,801	T→C	K18K (AAA→AAG)	SAUSA300_RS07715 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_007793	1,578,801	0	T	C	94.1%	39.9 / ‑5.6	17	K18K (AAA→AAG)	SAUSA300_RS07715	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base C (7/9);  minor base A (1/0);  total (8/9)
Fisher's exact test for biased strand distribution p-value = 4.71e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTATCTTTAATTTCTTTAATAAAATCATCTTCAATCTTAGCCTTTTCTTCAGGTGCTTTATTTTTATATTTTTTGATTATTTCAGTGTACTTTTTCGGGAATTTCATTTTAGGTATGTTAATCATCGTTTGCCTCCTTAATAAATGTAAATGATTCAATCTCATCTCT  >  NC_007793/1578730‑1578897                                                                        |                                                                                                 ttatcattaATTTCTTTAATGAAATCATCTTCAATCTTAGCTTTTTCTTCAGGTGTTTTATTTTTATATTTCTTGATTATTTCAGTGTACTTTTTCGGGa                                                                      >  1:1740296/1‑100 (MQ=255)                            cttccgaTCTTAGCTTTTTCTTCAGGTGTTTTATTTTTATATTTCTTGATTATTTCAGTGTACTTTTTCGGGAATTTCATTTTAGGTATGTTAatcgtca                                           <  1:602314/94‑2 (MQ=255)                              tccgaTCTTAGCTCTTTCTTCAGGTGTTTTATCTTTATATTTCTTGATTATTTCAGTGTACTTTTTCGGGAATTTCATTTTAGGTATGTTAATCGTCGtt                                         <  1:331796/96‑1 (MQ=255)                               tCAATCTTAGCTTTTTCTTCAGGTGTTTTATTTTTATATTTCTTGATTATTTCAGTGTACTTTTTCGGGAATTTCATTTTAGGTATGTTAATCGTCGttt                                        <  1:769779/100‑1 (MQ=255)                                   ctatAGCTTTTTCTTCAGGTGTTTTATTTTTATATTTCTTGATTATTTCAGTGTACTTTTTCGGGAATTTCATTTTAGGTATGTTAATCGTCGTTTGcct                                    <  1:5756763/97‑1 (MQ=255)                                   ccTTAGCTCTTTCTTCAGGTGTTTTATTTTTATATTGCTTGATTAATTCAGTGTACTTTTTCGGGAATGTCATTTTATGTATGTTAATCGTCGTTTGcct                                    <  1:3136575/99‑1 (MQ=255)                                     ctAGCTTTTTCTTCAGGTGTTTTATTTTTATATTTCTTGATTATTTCAGTGTACTTTTTCGGGAATTTCATTTTAGGTATGTTAATCGTCGTTTGcctcc                                  >  1:2253350/2‑100 (MQ=255)                                      tggcgttttCTTCAGGTGTTTTATTTTTATATTTCTTGATTATTTCAGTGTACTTTTTCGGGAATTTCATTTTAGGTATGTTAATCGTCGTTTGcctcct                                 >  1:5618552/6‑100 (MQ=255)                                      tAGCTTTTTCTTCAGGTGTTTTATTTTTATATTTCTTGATTATTTCAGTGTACTTTTTCGGGAATTTCATTTTAGGTATGTTAATCGTCGTTTGcctcct                                 >  1:6168786/1‑100 (MQ=255)                                      tAGCTTTTTCTTCAGGTGTTTTATTTTTATATTTATTGATTATTTCAGTGTACTTTTTCGGGAATTTCATTTTAGGTATGTTAATCGTCGTTTGcctcct                                 >  1:4335275/1‑100 (MQ=255)                                      acgctttttCTTCAGGTGTTTTATTTTTATATTTCTTGATTATTTCAGTGTACTTTTTCGGGAATTTCATTTTAGGTATGTTAATCGTCGTTTGcctcct                                 >  1:6395019/6‑100 (MQ=255)                                       ggctttttCTTCAGGTGTTTTATTTTTATATTTCTTGATTATTTCAGTGTACTTTTTCGGGAATTTCATTTTAGGTATGTTAATCGTCGTTTGCCTCCtt                                <  1:1892265/96‑1 (MQ=255)                                        gctttttCTTCAGGTGTTTTATTTTTATATTTCTTGATTATTTCAGTGTACTTTTTCGGGAATTTCATTTTAGGTATGTTAATCGTCGTTTGCCTCCTta                               <  1:5063034/97‑1 (MQ=255)                                                ctcAGGTGTTTTATTTTTATATTTCTTGATTATTTCAGTGTACTTTTTCGGGAATTTCATTTTAGGTATGTTAATCGTCGTTTGCCTCCTTAAtaaatgt                       >  1:5422518/2‑100 (MQ=255)                                                 tttggTGTTTTATTTTTATATTTCTTGATTATTTCAGTGTACTTTTTCGGGAATTTCATCTTAGGTATGTTAATCGTCGTTTGCCTCCTTAAtaaatgta                      >  1:6271426/4‑100 (MQ=255)                                                 tcAGGTGTTTTATTTTTATATTTCTTGATTATTTCAGTGTACTTTTTCGGGAATTTCATTTTAGGTATGTTAATCGTCGTTTGCCTCCTTAAtaaatgta                      <  1:444880/100‑1 (MQ=255)                                                           tATTTTTATATTTCTTGATTATTTCAGTGTACTTTTTCGGGAATTTCATTTTAGGTATGTTAATCGTCGTTTGCCTCCTTAATAAATGTAAATGATTCac            <  1:6180366/100‑2 (MQ=255)                                                                     tttcttGATTATTTCAGTGTACTTTTTCGGGAATTTCATTTTAGGTATGTTAATCGTCGTTTGCCTCCTTAATAAATGTAAATGATTCAatctcatctcg  <  1:3764962/100‑2 (MQ=255)                                                                        |                                                                                                 TTATCTTTAATTTCTTTAATAAAATCATCTTCAATCTTAGCCTTTTCTTCAGGTGCTTTATTTTTATATTTTTTGATTATTTCAGTGTACTTTTTCGGGAATTTCATTTTAGGTATGTTAATCATCGTTTGCCTCCTTAATAAATGTAAATGATTCAATCTCATCTCT  >  NC_007793/1578730‑1578897

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: ‑