Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_007793	61,025	G→A	L62F (CTT→TTT)	SAUSA300_RS00265 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_007793	61,025	0	G	A	97.7%	128.6 / ‑3.5	43	coding (184/360 nt)	SAUSA300_RS00265	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (33/9);  minor base . (1/0);  total (34/9)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.02e-01

TTAAAAGAATGCTAAATATTAATAATTGATGATCCAAAGGTTTATCAAAGACTTCTTTAGAAAAAATAAATGCCAATCTACCTATAAAAAGACTAACCCAAGTCATAATATGATTAATGAATGCTTTTTTGCCATGTTCTGATGTTGGTTCATCAGCAAGTATCCATACAAATAAAGGTAAAA  >  NC_007793/60935‑61117                                                                                           |                                                                                             ttAAAAGAATGCTAAATATTAATAATTGATGATCCAAAGGTTTATCAAAGACTTCTTTAGAAAAAATAAATGCCAATCTACCTATAAAAAAACTAACCCa                                                                                     <  1:2279112/100‑1 (MQ=255)   aaaaGAATGCTAAATATTAATAATTGACGATCCAAAGGTTTATCAAAGACTTCTTTAGAAAAAATAAATGCCAATCTACCTATAAAAAAACTAACCCAAg                                                                                   <  1:5208666/100‑1 (MQ=255)      ttaaTGCTAAATATTAATAATTGATGATCCAAAGGTTTATCAAAGACTTCTTTAGAAAAAATAAATGCCAATCTACCTATAAAAAAACTAACCCAAgatc                                                                                >  1:1689514/3‑97 (MQ=255)      aGAATGCTAAATATTAATAATTGATGATCCAAAGGTTTATCAAAGACTTCTTTAGAAAAAATAAATGCCAATCTACCTATAAAAAAACTAACCCAAgatc                                                                                >  1:2655926/1‑97 (MQ=255)           gCTAAATATTAATAATTGATGATCCAAAGGTTTATCAAAGACTTCTTTAGAAAAAATAAATGCCAATCTACCTATAAAAAAACTAACCCAAgatcggaga                                                                           >  1:2499567/1‑92 (MQ=255)              aaaTATTAATAATTGATGATCCAAAGGTTTATCAAAGACTTCTTTAGAAAAAATAAATGCCAATCTACCTATAAAAAAACTAACCCAAGTCATAATATGa                                                                        >  1:1930282/1‑100 (MQ=255)                                   gtaaGGTTTATCAAAGACTTCTTTAGAAAAAATAAATGCCAATCTACCTATAAAAA‑ACTAACCCAAGTCATAATATGATTAATGAATGCTTTTTTGCCAt                                                  >  1:5758950/3‑100 (MQ=255)                                   ataaGGTTTATCAAAGACTTCTTTAGAAAAAATAAATGCCAATCTACCTATAAAAAAACTAACCCAAGTCATAATATGATTAATGAATGCTTTTTTGCCa                                                   >  1:3783944/3‑100 (MQ=255)                                      aGGTTTATCAAAGACTTCTTTAGAAAAAATAAATGCCAATCTACCTATAAAAAAACTAACCCAAGTCATAATATGATTAATGAATGCTTTTTTGCCATGt                                                <  1:1629589/100‑1 (MQ=255)                                         tttATCAAAGACTTCTTTAGAAAAAATAAATGCCAATCTACCTATAAAAAAACTAACCCAAGTCATAATATGATTAATGAATGC‑TTTTTGCCATGTTCta                                            <  1:3457417/100‑2 (MQ=255)                                          ttATCAAAGACTTCTTTAGAAAAAATAAATGCCAATCTACCTATAAAAAAGCTAACCCAAGTCATAATATGATTAATGAATGCTTTTTTGCCATGTTCTg                                            <  1:581762/100‑1 (MQ=255)                                                                  ataaatGCCAATCTACCTATAAAAAAACTAACCCAAGTCATAATATGATTAATGAATGCTTTTTTGCCATGTTCTGATGTTGGTTCATCAGCACGTATcc                    <  1:1696549/100‑1 (MQ=255)                                                                       ttccgaTCTGCCTATAAAAAAACTAACCCAAGTCATAATATGATTAATGAATGCTTTTTTGCCATGTTCTGATGTTGGTTCATCAGCAAGTATCCATACa               <  1:5752180/95‑1 (MQ=255)                                                                            aTCTACCTATAAAAAAACTAACCCAAGTCATAATATGATTAATGAATGCTTTTTTGCCATGTTCTGATGTTGGTTCATCAGCAAGTATCCATACAAATaa          <  1:54338/100‑1 (MQ=255)                                                                              acgccTATAAAAAAACTAACCCAAGTCATAATATGATTAATGAATGC‑TTTTTGCCATGTTCTGATGTTGGTTCATCAGCAAGTATCCATACAAATAAAgg       >  1:3671246/4‑100 (MQ=255)                                                                               caCCTATAAAAAAACTAACCCAAGTCATAATATGATTAATGAATGCTTTTTTGCCATGTTCTGATGTTGGTTCATCAGCAAGTATCCATACAAATAAAgg       >  1:2055937/2‑100 (MQ=255)                                                                                acctgcaaaaaaACTAACCCAAGTCATAATATGATTAATGAATGC‑TTTTTGCCATGTTCTGATGTTGGTTCATCAGCAAGTATCCATACAAATAAAGGTa     >  1:71757/7‑100 (MQ=255)                                                                                aCCTATAAAAAAACTAACCCAAGTCATAATATGATTAATGAATGCTTTTTTGCCATGTTCTGATGTTGGTTCATCAGCAAGTATCCTTACAAATAAAGGt      >  1:1146414/1‑100 (MQ=255)                                                                                aCCTATAAAAAAACTAACCCAAGTCATAATATGATTAATGAATGCTTTTTTGCCATGTTCTGATGTTGGTTCATCAGCAAGTATCCATACAAATAAAGGt      >  1:5124311/1‑100 (MQ=255)                                                                                aCCTATAAAAAAACTAACCCAAGTCATAATATGATTAATGAATGCTTTTTTGCCATGTTCTGATGTTGGTTCATCAGCAAGTATCCATACAAATAAAGGt      >  1:834150/1‑100 (MQ=255)                                                                                aCCTATAAAAAAACTAACCCAAGTCATAATATGATTAATGAATGCTTTTTTGCCATGTTCTGATGTTGGTTCATCAGCAAGTATCCATACAAATAAAGGt      >  1:91196/1‑100 (MQ=255)                                                                                aCCTATAAAAAAACTAACCCAAGTCATAATATGATTAATGAATGCTTTTTTGCCATGTTCTGATGTTGGTTCATCAGCAAGTATCCATACAAATAAAGGt      >  1:5794039/1‑100 (MQ=255)                                                                                aCCTATAAAAAAACTAACCCAAGTCATAATATGATTAATGAATGCTTTTTTGCCATGTTCTGATGTTGGTTCATCAGCAAGTATCCATACAAATAAAGGt      >  1:5760708/1‑100 (MQ=255)                                                                                aCCTATAAAAAAACTAACCCAAGTCATAATATGATTAATGAATGCTTTTTTGCCATGTTCTGATGTTGGTTCATCAGCAAGTATCCATACAAATAAAGGt      >  1:1620508/1‑100 (MQ=255)                                                                                aCCTATAAAAAAACTAACCCAAGTCATAATATGATTAATGAATGCTTTTTTGCCATGTTCTGATGTTGGTTCATCAGCAAGTATCCATACAAATAAAGGt      >  1:1626412/1‑100 (MQ=255)                                                                                aCCTATAAAAAAACTAACCCAAGTCATAATATGATTAATGAATGCTTTTTTGCCATGTTCTGATGTTGGTTCATCAGCAAGTATCCATACAAATAAAGGt      >  1:1895271/1‑100 (MQ=255)                                                                                aCCTATAAAAAAACTAACCCAAGTCATAATATGATTAATGAATGCTTTTTTGCCATGTTCTGATGTTGGTTCATCAGCAAGTATCCATACAAATAAAGGt      >  1:539979/1‑100 (MQ=255)                                                                                aCCTATAAAAAAACTAACCCAAGTCATAATATGATTAATGAATGCTTTTTTGCCATGTTCTGATGTTGGTTCATCAGCAAGTATCCATACAAATAAAGGt      >  1:2151377/1‑100 (MQ=255)                                                                                aCCTATAAAAAAACTAACCCAAGTCATAATATGATTAATGAATGCTTTTTTGCCATGTTCTGATGTTGGTTCATCAGCAAGTATCCATACAAATAAAGGt      >  1:5160030/1‑100 (MQ=255)                                                                                aCCTATAAAAAAACTAACCCAAGTCATAATATGATTAATGAATGCTTTTTTGCCATGTTCTGATGTTGGTTCATCAGCAAGTATCCATACAAATAAAGGt      >  1:2684012/1‑100 (MQ=255)                                                                                aCCTATAAAAAAACTAACCCAAGTCATAATATGATTAATGAATGCTTTTTTGCCATGTTCTGATGTTGGTTCATCAGCAAGTATCCATACAAATAAAGGt      >  1:4488776/1‑100 (MQ=255)                                                                                aCCTATAAAAAAACTAACCCAAGTCATAATATGATTAATGAATGCTTTTTTGCCATGTTCTGATGTTGGTTCATCAGCAAGTATCCATACAAATAAAGGt      >  1:4382327/1‑100 (MQ=255)                                                                                aCCTATAAAAAAACTAACCCAAGTCATAATATGATTAATGAATGCTTTTTTGCCATGTTCTGATGTTGGTTCATCAGCAAGTATCCATACAAATAAAGGt      >  1:3984971/1‑100 (MQ=255)                                                                                aCCTATAAAAAAACTAACCCAAGTCATAATATGATTAATGAATGCTTTTTTGCCATGTTCTGATGTTGGTTCATCAGCAAGTATCCATACAAATAAAGGt      >  1:3969430/1‑100 (MQ=255)                                                                                aCCTATAAAAAAACTAACCCAAGTCATAATATGATTAATGAATGCTTTTTTGCCATGTTCTGATGTTGGTTCATCAGCAAGTATCCATACAAATAAAGGt      >  1:893690/1‑100 (MQ=255)                                                                                aCCTATAAAAAAACTAACCCAAGTCATAATATGATTAATGAATGCTTTTTTGCCATGTTCTGATGTTGGTTCATCAGCAAGTATCCATACAAATAAAGGt      >  1:3884819/1‑100 (MQ=255)                                                                                aCCTATAAAAAAACTAACCCAAGTCATAATATGATTAATGAATGCTTTTTTGCCATGTTCTGATGTTGGTTCATCAGCAAGTATCCATACAAATAAAGGt      >  1:2321972/1‑100 (MQ=255)                                                                                aCCTATAAAAAAACTAACCCAAGTCATAATATGATTAATGAATGCTTTTTTGCCATGTTCTGATGTTGGTTCATCAGCAAGTATCCATACAAATAAAGGt      >  1:2417671/1‑100 (MQ=255)                                                                                aCCTATAAAAAAACTAACCCAAGTCATAATATGATTAATGAATGCTTTTTTGCCATGTTCTGATGTTGGTTCATCAGCAAGTATCCATACAAATAAAGGt      >  1:2630229/1‑100 (MQ=255)                                                                                aCCTATAAAAAAACTAACCCAAGTCATAATATGATTAATGAATGCTTTTTTGCCATGTTCTGATGTTGGTTCATCAGCAAGTATCCATACAAATAAAGGt      >  1:2858929/1‑100 (MQ=255)                                                                                aCCTATAAAAAAACTAACCCAAGTCATAATATGATTAATGAATGCCTTTTTGCCATGTTCTGATGTTGGTTCATCAGCAAGTATCCATACAAATAAAGGt      >  1:930599/1‑100 (MQ=255)                                                                                aCCTATAAAAAAACTAACCCAAGTCATAATATGATTAATGAATGC‑TTTTTGCCATGTTCTGATGTTGGTTCATCAGCAAGTATCCATACAAATAAAGGTa     >  1:396572/1‑100 (MQ=255)                                                                                 ccTATAAAAAAACTAACCCAAGTCATAATATGATTAATGAATGCTTTTTTGCCATGTTCTGATGTTGGTTCATCAGCAAGTATCCATACAAATAAAGGtc     <  1:1335930/100‑2 (MQ=255)                                                                                    ataAAAAAACTAACCCAAGTCATAATATGATTAATGAATGCTTTTTTGCCATGTTCTGATGTTGGTTCATCAGCAAGTATCCATACAAATAAAGGTaaag  >  1:1307845/1‑99 (MQ=255)                                                                                           |                                                                                             TTAAAAGAATGCTAAATATTAATAATTGATGATCCAAAGGTTTATCAAAGACTTCTTTAGAAAAAATAAATGCCAATCTACCTATAAAAAGACTAACCCAAGTCATAATATGATTAATGAATGCTTTTTTGCCATGTTCTGATGTTGGTTCATCAGCAAGTATCCATACAAATAAAGGTAAAA  >  NC_007793/60935‑61117

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