Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_007793	1,579,752	T→C	P25P (CCA→CCG)	SAUSA300_RS07730 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_007793	1,579,752	0	T	C	95.0%	40.6 / ‑4.0	20	P25P (CCA→CCG)	SAUSA300_RS07730	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base C (2/17);  total (2/18)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.77e-01

TTCCAGTTCTTAAACGTTTCCAATGTATTAAACCATTCAAAGAATTCAACTACATATGCTATTTTGCTATCTCGTAAGTTCAGCGTTGGTATTCTATCGTCGAACTCCTTTTGTATTTGTAATAACTCTTGTAACTGATCAATTGTTAATGTATTAGTCATTTTCCTGATCCTCCTCATATTT  >  NC_007793/1579666‑1579848                                                                                       |                                                                                                 ttcCAGTTCTTAAACGTTTCCAATGTATTAAACCATTCAAAGAATTCAACTACATATGCAATTTTGCTATCTCCTAAGTTTAGCGTCGGTATTCTAtcgt                                                                                     <  1:5363549/100‑1 (MQ=255)   ccAGTTCTTAAACGTTTCCAATGTATTAAACCATTCAAAGAATTCAACTACATATGCAATTTTGCTATCTCCTAAGTTTAGCGTCGGTATTCTAtcgtcg                                                                                   >  1:2129334/1‑100 (MQ=255)   ccAGTTCTTAAACGTTTCCAATGTATTAAACCATTCAAAGAATTCAACTACATATGCAATTTTGCTATCTCCTAAGTTTAGCGTCGGTATTCTAtcgtcg                                                                                   <  1:3851165/100‑1 (MQ=255)                 tttCCAATGTATTAAACCATTCAAAGAATTCAACTACATATGCAATTTTGCTATCTCCTAAGTTTAGCGTCGGTATTCTATCGTCGAACTCCTtttgtat                                                                     <  1:3896818/100‑1 (MQ=255)                   tCCAATGTATTAAACCATTCAAAGAATTCAACTACATATGCAATTTTGCTATCTCCTAAGTTTAGCGTCGGTATTCTATCGTCGAACTCCTtctgtattt                                                                   <  1:442589/100‑1 (MQ=255)                     cAATGTATTAAACCATTCAAAGAATTCAACTACATATGCAATTTTGCTATCTCCTAAGTTTAGCGTCGGTATTCTATCGTCGAACTCCTtttgtatttgt                                                                 <  1:3363503/100‑1 (MQ=255)                     cAATGTATTAAACCATTCAAAGAATTCAACTACATATGCAATTTTGCTATCTCCTAAGTTTAGCGTCGGTATTCTATCGTCGAACTCCTtttgtatttgt                                                                 <  1:3997210/100‑1 (MQ=255)                      aaTGTATTAAACCATTCAAAGAATTCAACTACATATGAAATTTTGCTATCTCCTAAGTTTAGCGTCGGTATTCTATCGTCGAACTCCTtttgtatttgta                                                                <  1:3076615/100‑1 (MQ=255)                                 ccATTCAAAGAATTCAACTACATATGCAATTTTGCTATCTCCTAAGTTTAGCGTCGGTATACTATCGTCGAACTCCTTTTGTATTTGTAATAACTCTTGt                                                     <  1:475983/100‑1 (MQ=255)                                  cATTCAAAGAATTCAACTACATATGCAATTTTGCTATCTCCTAAGTTTAGCGTCGGTATTCTATCGTCGAACTCCTTTTGTATTTGTAATAACTCTTGTa                                                    <  1:1816431/100‑1 (MQ=255)                                    ttCAAAGAATTCAACTACATATGCAATTTTGCTATCTCCTAAGTTTAGCGTCGGTATTCTATCGTCGAACTCCTTTTGTATTTGTAATAACTCTTGTAAc                                                  <  1:3362928/100‑1 (MQ=255)                                      cAAAGAATTCAACTACATATGCAATTTTGCTATCTCCTAAGTTTAGCGTCGGTATTCTATCGTCGAACTCCTTTTGTATTTGTAATAACTCTTGTAACTg                                                <  1:4751227/100‑1 (MQ=255)                                          gAATTCAACTACATATGCAATTTTGCTATCTCCTAAGTTTAGCGTCGGTATTCTATCGTCGAACTCCTTTTGTATTTGTAATAACTCTTGTAACTGATcg                                            <  1:1613239/100‑2 (MQ=255)                                                        atGCAATTTTGCTATCTCCTAAGTTTAGCGTCGGTATTCTATCGTCGAACTCCTTTTGTATTTGTAATAACTCTTGTAACTGATCGATTGTTAATGtgtc                              <  1:2170858/100‑4 (MQ=255)                                                                         cctAAGTTTAGCGTCGGTATTCTATCGTCGAACTCCTTTTGTATTTGTAATAACTCTTGTAACTGATCGATTGTTAATGTGTTATTCATTATcgttgtgt             <  1:485701/98‑8 (MQ=255)                                                                         cctAAGTTTAGCGTCGGTATTCTATCGTCGAACTCCTTTTGTATTTGTAATAACTCTTGTAACTGATCGATTGTTAATGTGTTATTCATTATcgttgtcc             >  1:1415375/3‑93 (MQ=255)                                                                              atttagtgTTGGTATTCTATCGTCGAACTCCTTTTGTATTTGTAATAACCCTTGTAACTGATCGATTGTTAATGTGTTATTCATTATCGTTGtcctcccg        <  1:5190440/93‑3 (MQ=255)                                                                                ttaGCGTCGGTATTCTATCGTCGAACTCCTTTTGTATTTGTAATAACTCTTGTAACTGATCGATTGTTAATGTGTTATTCATTATCGTTGtcctcccatg      <  1:194809/98‑5 (MQ=255)                                                                                    cGTCGGTATTCTATCGTCGAACTCCTTTTGTATTTGTAATAACTCTTGTAACTGATCGATTGTTAATGTGTTATTCATTATCGTTGtcctcccattgcag  <  1:1789366/100‑9 (MQ=255)                                                                                    cGTCGGTATTCTATCGTCGAACTCCTTTTGTATTTGTAATAACTCTTGTAACTGATCGATTGTTAATGTGTTATTCATTATCGTTGtcctcccattgcaa  <  1:3289268/100‑9 (MQ=255)                                                                                       |                                                                                                 TTCCAGTTCTTAAACGTTTCCAATGTATTAAACCATTCAAAGAATTCAACTACATATGCTATTTTGCTATCTCGTAAGTTCAGCGTTGGTATTCTATCGTCGAACTCCTTTTGTATTTGTAATAACTCTTGTAACTGATCAATTGTTAATGTATTAGTCATTTTCCTGATCCTCCTCATATTT  >  NC_007793/1579666‑1579848

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