Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_007793	1,579,434	A→G	P131P (CCT→CCC)	SAUSA300_RS07730 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_007793	1,579,434	0	A	G	95.8%	54.2 / ‑5.7	24	P131P (CCT→CCC)	SAUSA300_RS07730	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base G (8/15);  total (8/16)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

TTCCTTTTCATTTTCTTTTTGTATGCGTCAATGAGCTGGTCGATAGAATATAAGTTGAAAGCTATGTCTATCGCTATTACAATTGCCAATTGGTCAGGATAAAATTCTTTGAATATTATCTGTGGTGTGCTAACAACTGCGTCTTGAGCAAATTCTTTATCTTTAAAATTAAACATGTTGTGAA  >  NC_007793/1579339‑1579522                                                                                                |                                                                                         ttCCTTTTCATTTTCTTTTTGTATGGGTCAATGAGTTGGTCGATAGTATATAAGTTGTAAGCTATGTCTATCGCTATCACAATTGCCAATTGGTCgggta                                                                                      <  1:2115406/100‑3 (MQ=255) ttCCTTTTCATTTTCTTTTTGTATGCGTCAATGAGTTGGTCGATAGTATATAAGTTGTAAGCTATGTCTATCGCTATCACAATTGCCAATTGGTCGGGAt                                                                                      <  1:1915916/100‑1 (MQ=255)  tCCTTTTCATTTTCTTTTTGTATGCGTCAATGAGTTGGTCGATAGTATATAAGTTGTAAGCTATGTCTATCGCTATCACAATTGCCAATTGGTCGGGATa                                                                                     <  1:3549972/100‑1 (MQ=255)  tCCTTTTCATTTTCTTTTTGTATGCGTCAATGAGTTGGTCGATAGTATATAAGTTGTAAGCTATGTCTATCGCTATCACAATTGCCAATTGGTCGGGATa                                                                                     <  1:2646164/100‑1 (MQ=255)          aTTTTCTTTTTGTATGCGTCAATGAGTTGGTCGATAGTATATAAGTTGTAAGCTATGTCTATCGCTATCACAATTGCCAATTGGTCGGGATAAAATTCtg                                                                             <  1:3797489/100‑2 (MQ=255)                         gCGTCAATGAGTTGGTCGATAGTATATAAGTTGTAAGCTATGTCTATCGCTATCACAATTGCCAATTGGTCGGGATAAAATTCTTTGAATATTATCtgtg                                                              <  1:589753/100‑1 (MQ=255)                             cAATGAGTTGGTCGATAGTATATAAGTTGTAAGCTATGTCTATCGCTATCACGATTGCCAATTGGTCAGGATAAAATTCTTTGAATATTATCTGTGgtgt                                                          <  1:2942162/100‑1 (MQ=255)                             cAATGAGTTGGTCGATAGTATATAAGTTGTAAGCTATGTCTATCGCTATCACAATTGCCAATTGGTCGGGATAAAATTCTTTGAATATTATCTGTGgtgt                                                          >  1:1253634/1‑100 (MQ=255)                             cAATGAGTTGGTCGATAGTATATAAGTTGTAAGCTATGTCTATCGCTATCACAATTGCCAATTGGTCGGGATAAAATTCTTTGAATATTATCTGTGgtgt                                                          <  1:5429688/100‑1 (MQ=255)                              aaTGAGTTGGTCGATAGTATATAAGTTGTAAGCTATGTCTATCGCTATCACAATTGCCAATTGGTCGGGATAAAATTCTTTGAATATTATCTGTGgtgta                                                         <  1:2617567/100‑2 (MQ=255)                                 gAGTTGGTCGATAGTATATAAGTTGTAAGCTATGTCTATCGCTATCACAATTGCCAATTGGTCGGGATAAAATTCTTTGAATATTATCTGTGGTGTACTa                                                      <  1:3701064/100‑1 (MQ=255)                                           aTAGTATATAAGTTGTAAGCTATGTCTATCGCTATCACAATTGCCAATTGGTCGGGATAAAATTCTTTGAATATTATCTGTGGTGTACTAACAACTGCgc                                            <  1:190237/100‑2 (MQ=255)                                            tAGTATATAAGTTGTAAGCTATGTCTATCGCTATCACAATTGCCAATTGGTCGGGATAAAATTCTTTGAATATTATCTGTGGTGTACTAACAACTGCgcc                                           <  1:1149751/100‑3 (MQ=255)                                              gtatatAAGTTGTAAGCTATGTCTATCGCTATCACAATTGCCAATTGGTCGGGATAAAATTCTTTGAATATTATCTGTGGTGTACTAACAACTGCGCCtg                                         <  1:3847484/98‑2 (MQ=255)                                               tatatAAGTTGTAAGCTATGTCTATCGCTATCACAATTGCCAATTGGTCGGGATAAAATTCTTTGAATATTATCTGTGGTGTACTAACAACTGCGCCTTg                                        <  1:1733830/99‑1 (MQ=255)                                                atatAAGTTGTAAGCTATGTCTATCGCTATCACAATTGCCAATTGGTCGGGATAAAATTCTTTGAATATTATCTGTGGTGTACTAACAACTGCGCCTTGa                                       <  1:2293079/100‑1 (MQ=255)                                                  ataAGTTGTAAGCTATGTCTATCGCTATCACAATTGCCAATTGGTCGGGATAAAATTCTTTGAATATTATCTGTGGTGTACTAACAACTGCGCCttagat                                     <  1:2678909/100‑5 (MQ=255)                                                        tgtaaGCTATGTCTATCGCTATCACAATTGCCAATTGGTCGGGATAAAATTCTTTGAATATTATCTGTGGTGTACTAACAACTGCGCCTTGAGCAAATtc                               >  1:1861665/4‑100 (MQ=255)                                                                aTGTCTATCGCTATCACGATTGCCAATTGGTCGGGATAAAATTCTTTGAATATTATCTGTGGTGTACTAACAACTGCGCCTTGAGCAAATtctttatctt                       <  1:5723018/100‑1 (MQ=255)                                                                                  aTTGCCAATTGGTCGGGATAAAATTCTTTGAATATTATCTGTGGTGTACTAACAACTGCGCCTTGAGCAAATTCTTTATCTTTAAAATTAAACAttttgt     >  1:345257/1‑100 (MQ=255)                                                                                    tGCCAATTGGTCGGGATAAAATTCTTTGAATATTATCTGTGGTGTACTAACAACTGCGCCTTGAGCAAATTCTTTATCTTTAAAATTAAACATTTTGTGa   >  1:2331530/1‑100 (MQ=255)                                                                                    tGCCAATTGGTCGGGATAAAATTCTTTGAATATTATCTGTGGTGTACTAACAACTGCGCCTTGAGCAAATTCTTTATCTTTAAAATTAAACATTTTGTGa   >  1:5145958/1‑100 (MQ=255)                                                                                     tccAATTGGTCGGGATAAAATTCTTTGAATATTATCTGTGGTGTACTAACAACTGCGCCTTGAGCAAATTCTTTATCTTTAAAATTAAACATTTTGTGaa  >  1:1428089/2‑100 (MQ=255)                                                                                     gCCAATTGGTCGGGATAAAATTCTTTGAATATTATCTGTGGTGTACTAACAACTGCGCCTTGAGCAAATTCTTTATCTTTAAAATTAAACATTTTGTGaa  >  1:3859205/1‑100 (MQ=255)                                                                                     gCCAATTGGTCGGGATAAAATTCTTTGAATATTATCTGTGGTGTACTAACAACTGCGCCTTGAGCAAATTCTTTATCTTTAAAATTAAACATTTTGTGaa  >  1:5329101/1‑100 (MQ=255)                                                                                                |                                                                                         TTCCTTTTCATTTTCTTTTTGTATGCGTCAATGAGCTGGTCGATAGAATATAAGTTGAAAGCTATGTCTATCGCTATTACAATTGCCAATTGGTCAGGATAAAATTCTTTGAATATTATCTGTGGTGTGCTAACAACTGCGTCTTGAGCAAATTCTTTATCTTTAAAATTAAACATGTTGTGAA  >  NC_007793/1579339‑1579522

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