Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	1	873,132	T→C	C565R (TGC→CGC)  P10P (CCT→CCC)	unknown → unknown →	CRISPR‑associated helicase Cas3 CRISPR‑associated protein, Csd2/Csh2 family

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	1	873,132	0	T	C	100.0%	24.8 / NA	11	C565R (TGC→CGC)  P10P (CCT→CCC)	unknown unknown	CRISPR‑associated helicase Cas3 CRISPR‑associated protein, Csd2/Csh2 family
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (4/7);  total (4/7)

GCGTCCATGTGTTTTGGCGGTGGTCAATACCCGCAAACATGCCCGTAGGCTTTTTGCCGCGCTGCCTTCAAACGGCATCAAACTTCATCTCTCTGCCAATATGTGTGCTACGCACCGTAGCGAAGTGATTGCACTGGTTCGCCGTTATTTGGAACTGTATCGCGAAGGTCGTCTGAAAAAGCCTTTATGGGTGGTCAGTACCCAACTGATTGAAGCAGGAGTCGATTTGGATTTTCCCTGCGTTTAAAGAAGTTGATAACCCAAATCCTGA  >  1/872894‑873164                                                                                                                                                                                                                                               |                                 gCGTCCATGTGTTTTGGCGGTGGTCAATACCCGCAAACATGCCCGTAGGCTTTTTGCCGCGCTGCCTTCAAACGGCATCAAACTTCATCTCTCTGCCAATATGTGTGCTACGCACCGTAGCGAAGTGATTGCACTGGTTCGCCGTTATTTGGAACTGTATCGCGAAGGTCGTCTGAAAAAGCCTTTATGGGTGGTCAGTACCCAACTGATTGAAGCAGGAGTCGATTTGGATTTTCCCCGCGTTTatcggg                      >  1:521913/1‑246 (MQ=255) gCGTCCATGTGTTTTGGCGGTGGTCAATACCCGCAAACATGCCCGTAGGCTTTTTGCCGCGCTGCCTTCAAACGGCATCAAACTTCATCTCTCTGCCAATATGTGTGCTACGCACCGTAGCGAAGTGATTGCACTGGTTCGCCGTTATTTGGAACTGTATCGCGAAGGTCGTCTGAAAAAGCCTTTATGGGTGGTCAGTACCCAACTGATTGAAGCAGGAGTCGATTTGGATTTTCCCCGCGTTTatcggg                      >  1:772801/1‑246 (MQ=255)             tttGGCGGTGGTCAATACCCGCAAACATGCCCGTAGGCTTTTTGCCGCGCTGCCTTCAAACGGCATCAAACTTCATCTCTCTGCCAATATGTGTGCTACGCACCGTAGCGAAGTGATTGCACTGGTTCGCCGTTATTTGGAACTGTATCGCGAAGGTCGTCTGAAAAAGCCTTTATGGGTGGTCAGTACCCAACTGATTGAAGCAGGAGTCGATTTGGATTTTCCCCGCGTTTatcgggcgatggcagggc          <  1:746052/251‑18 (MQ=255)                  cGGTGGTCAATACCCGCAAACATGCCCGTAGGCTTTTTGCCGCGCTGCCTTCAAACGGCATCAAACTTCATCTCTCTGCCAATATGTGTGCTACGCACCGTAGCGAAGTGATTGCACTGGTTCGCCGTTATTTGGAACTGTATCGCGAAGGTCGTCTGAAAAAGCCTTTATGGGTGGTCAGTACCCAACTGATTGAAGCAGGAGTCGATTTGGATTTTCCCCGCGTTTatcgggcgatggcagggcttgac     <  1:100906/251‑23 (MQ=255)                  cGGTGGTCAATACCCGCAAACATGCCCGTAGGCTTTTTGCCGCGCTGCCTTCAAACGGCATCAAACTTCATCTCTCTGCCAATATGTGTGCTACGCACCGTAGCGAAGTGATTGCACTGGTTCGCCGTTATTTGGAACTGTATCGCGAAGGTCGTCTGAAAAAGCCTTTATGGGTGGTCAGTACCCAACTGATTGAAGCAGGAGTCGATTTGGATTTTCCCCGCGTTTatcgggcgatggcagggcttgac     <  1:697399/251‑23 (MQ=255)                    gtggtCAATACCCGCAAACATGCCCGTAGGCTTTTTGCCGCGCTGCCTTCAAACGGCATCAAACTTCATCTCTCTGCCAATATGTGTGCTACGCACCGTAGCGAAGTGATTGCACTGGTTCGCCGTTATTTGGAACTGTATCGCGAAGGTCGTCTGAAAAAGCCTTTATGGGTGGTCAGTACCCAACTGATTGAAGCAGGAGTCGATTTGGATTTTCCCCGCGTTTatcgggcgatggcagggcttgacag   >  1:359762/1‑227 (MQ=255)                    gtggtCAATACCCGCAAACATGCCCGTAGGCTTTTTGCCGCGCTGCCTTCAAACGGCATCAAACTTCATCTCTCTGCCAATATGTGTGCTACGCACCGTAGCGAAGTGATTGCACTGGTTCGCCGTTATTTGGAACTGTATCGCGAAGGTCGTCTGAAAAAGCCTTTATGGGTGGTCAGTACCCAACTGATTGAAGCAGGAGTCGATTTGGATTTTCCCCGCGTTTatcgggcgatggcagggcttgaca    <  1:512109/250‑24 (MQ=255)                     tggtCAATACCCGCAAACATGCCCGTAGGCTTTTTGCCGCGCTGCCTTCAAACGGCATCAAACTTCATCTCTCTGCCAATATGTGTGCTACGCACCGTAGCGAAGTGATTGCACTGGTTCGCCGTTATTTGGAACTGTATCGCGAAGGTCGTCTGAAAAAGCCTTTATGGGTGGTCAGTACCCAACTGATTGAAGCAGGAGTCGATTTGGATTTTCCCCGCGTTTatcgggcgatggcagggcttgacagc  <  1:238173/251‑26 (MQ=255)                     tggtCAATACCCGCAAACATGCCCGTAGGCTTTTTGCCGCGCTGCCTTCAAACGGCATCAAACTTCATCTCTCTGCCAATATGTGTGCTACGCACCGTAGCGAAGTGATTGCACTGGTTCGCCGTTATTTGGAACTGTATCGCGAAGGTCGTCTGAAAAAGCCTTTATGGGTGGTCAGTACCCAACTGATTGAAGCAGGAGTCGATTTGGATTTTCCCCGCGTTTatcgggcgatggcagggcttgacag   <  1:873571/250‑25 (MQ=255)                                                                                               gCCAATATGTGTGCTACGCACCGTAGCGAAGTGATTGCACTGGTTCGCCGTTATTTGGAACTGTATCGCGAAGGTCGTCTGAAAAAGCCTTTATGGGTGGTCAGTACCCAACTGATTGAAGCAGGAGTCGATTTGGATTTTCCCCGCGTTTatcgggcgatggc               <  1:627508/164‑13 (MQ=255)                                                                                               gCCAATATGTGTGCTACGCACCGTAGCGAAGTGATTGCACTGGTTCGCCGTTATTTGGAACTGTATCGCGAAGGTCGTCTGAAAAAGCCTTTATGGGTGGTCAGTACCCAACTGATTGAAGCAGGAGTCGATTTGGATTTTCCCCGCGTTTatcg                        >  1:866076/1‑152 (MQ=255)                                                                                                                                                                                                                                               |                                 GCGTCCATGTGTTTTGGCGGTGGTCAATACCCGCAAACATGCCCGTAGGCTTTTTGCCGCGCTGCCTTCAAACGGCATCAAACTTCATCTCTCTGCCAATATGTGTGCTACGCACCGTAGCGAAGTGATTGCACTGGTTCGCCGTTATTTGGAACTGTATCGCGAAGGTCGTCTGAAAAAGCCTTTATGGGTGGTCAGTACCCAACTGATTGAAGCAGGAGTCGATTTGGATTTTCCCTGCGTTTAAAGAAGTTGATAACCCAAATCCTGA  >  1/872894‑873164

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