Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_007793	787,567	T→C	S61P (TCA→CCA)	SAUSA300_RS03810 →	5'‑3'‑deoxyribonucleotidase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_007793	787,567	0	T	C	100.0%	45.8 / NA	19	S61P (TCA→CCA)	SAUSA300_RS03810	5'‑3'‑deoxyribonucleotidase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (17/2);  total (17/2)

GCGGATTTAGGTATTAAAATGGAAGCTTTGAATGGTCAAAAACTTAAACATGTTATTCCTGAACATGATGGATTAATTACAGAAGTATTGAGAGAATCAGGCTTCTTCAGACATCTTAAAGTGATGCCGTATGCACAAGAAGTTGTGAAAAAATTAACTGAACATTATGATGTATATATTGCTACAGCAGC  >  NC_007793/787471‑787661                                                                                                 |                                                                                               gCGGATTTAGGTATTAAAATGGCAGCTTTGAATGGTCAAAAACTTAAACATGTTATTCCTGAACATGATGGATTAATTACAGAAGTATTGAGAGAACCAg                                                                                             >  1:929132/1‑100 (MQ=255)    gATTTAGGTATTAAAATGGAAGCTTTGAATGGTCAAAAACTTAAACATGTTATTCCTGAACATGATGGATTAATTACAGAAGTATTGAGAGAACCAGGct                                                                                          <  1:4342825/100‑1 (MQ=255)      gctAGGTATTAAAATGGAAGCTTTGAATGGTCAAAAACTTAAACATGTTATTCCTGAACATGATGGATTAATTACAGAAGTATTGAGAGAACCAGGcttc                                                                                        >  1:3051565/3‑100 (MQ=255)       ttAGGTATTAAAATGGAAGCTTTGAATGGTCAAAAACTTAAACATGTTATTCCTGAACATGATGGATTAATTACAGAAGTATTGAGAGAACCAGGcttct                                                                                       >  1:3937954/1‑100 (MQ=255)       ttAGGTATTAAAATGGAAGCTTTGAATGGTCAAAAACTTAAACATGTTATTCCTGAACATGATGGATTAATTACAGAAGTATTGAGAGAACCAGGcttct                                                                                       >  1:5337480/1‑100 (MQ=255)        tAGGTATTAAAATGGAAGCTTTGAATGGTCAAAAACTTAAACATGTTATTCCTGAACATGATGGATTAATTACAGAAGTATTGAGGGAACCAGGcttctt                                                                                      >  1:4851563/1‑100 (MQ=255)         aGGTATTAAAATGGAAGCTTTGAATGGTCAAAAACTTAAACATGTTATTCCTGAACATGATGGATTAATTACAGAAGTATTGAGAGAACCAGGCttcttc                                                                                     >  1:3937483/1‑100 (MQ=255)         aGGTATTAAAATGGAAGCTTTGAATGGTCAAAAACTTAAACATGTTATTCCTGAACATGATGGATTAATTACAGAAGTATTGAGAGAACCAGGCttcttc                                                                                     >  1:101261/1‑100 (MQ=255)         aGGTATTAAAATGGAAGCTTTGAATGGTCAAAAACTTAAACATGTTATTCCTGAACATGATGGATTAATTACAGAAGTATTGAGAGAACCAGGCttcttc                                                                                     >  1:4027346/1‑100 (MQ=255)         aGGTATTAAAATGGAAGCTTTGAATGGTCAAAAACTTAAACATGTTATTCCTGAACATGATGGATTAATTACAGAAGTATTGAGAGAACCAGGCttcttc                                                                                     >  1:5388576/1‑100 (MQ=255)         aGGTATTAAAATGGAAGCTTTGAATGGTCAAAAACTTAAACATGTTATTCCTGAACATGATGGATTAATTACAGAAGCATTGAGAGCACCAGGCttcttc                                                                                     >  1:3423026/1‑100 (MQ=255)                     ggAAGCTTTGAATGGTCAAAAACTTAAACATGTTATTCCTGAACATGATGGATTAATTACAGAAGTATTGAGAGAACCAGGCTTCTTCAGACATCTTaaa                                                                         >  1:1818485/1‑100 (MQ=255)                                                             gAACATGATGGATTAATTACAGAAGTATTGAGAGAACCAGGCTTCTTCAGACATCTTAAAGTGATGCCGTATGCACAAGAAGTTGTGAAAAAATTAACTg                                 >  1:3653698/1‑100 (MQ=255)                                                             gAACATGATGGATTAATTACAGAAGTATTGAGAGAACCAGGCTTCTTCAGACATCTTAAAGTGATGCCGTATGCACAAGAAGTTGTGAAAAAATTAACTg                                 >  1:2769498/1‑100 (MQ=255)                                                             gAACATGATGGATTAATTACAGAAGTATTGAGAGAACCAGGCTTCTTCAGACATCTTAAAGTGATGCCGTATGCACAAGAAGTTGTGAAAAAATTAACTg                                 >  1:887786/1‑100 (MQ=255)                                                                cATGATGGATTAATTACAGAAGTATTGAGAGAACCAGGCTTCTTCAGACATCTTAAAGTGATGCCGTATGCACAAGAAGATGTGAAAAAATTAACTGAAc                              >  1:5155908/1‑100 (MQ=255)                                                                                    aGTATTGAGAGAACCAGGCTTCTTCAGACATCTTAAAGTGATGCCGTATGCACAAGAAGTTGTGAAAAAATTAACTGAACATTATGATGTATATATTGCt          >  1:4070708/1‑100 (MQ=255)                                                                                       aTTGAGAGAACCAGGCTTCTTCAGACATCTTAAAGTGATGCCGTATGCACAAGAAGTTGTGAAAAAATTAACTGAACATTATGATGTATATATTGCTAca       >  1:3202261/1‑100 (MQ=255)                                                                                            gagaACCAGGCTTCTTCAGACATCTTAAAGTGATGCCGTATGCACAAGAAGTTGTGAAAAAATTAACTGAACATTATGATGTATATATTGCTAcagcagg  <  1:135559/100‑2 (MQ=255)                                                                                                 |                                                                                               GCGGATTTAGGTATTAAAATGGAAGCTTTGAATGGTCAAAAACTTAAACATGTTATTCCTGAACATGATGGATTAATTACAGAAGTATTGAGAGAATCAGGCTTCTTCAGACATCTTAAAGTGATGCCGTATGCACAAGAAGTTGTGAAAAAATTAACTGAACATTATGATGTATATATTGCTACAGCAGC  >  NC_007793/787471‑787661

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