Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	AM260479	2,594,128	Δ1 bp	coding (326/957 nt)	cR ←	Recombinational DNA repair protein (RecF pathway)

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	AM260479	2,594,128	0	C	.	100.0%	62.6 / NA	15	coding (326/957 nt)	cR	Recombinational DNA repair protein (RecF pathway)
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base . (9/6);  total (9/6)

GCTGGGCGGACTTGGGCCCGACGCCGGGCAGCACCCGCAGCGCTTCGACCAGTGCTTGCAGCGAGGTCGGGGACGATGCCTTCACAGTGGCGTTCAGAACGCGCAGCTTGAAGCCCGGGGGCAGCGGCAGGCCCGAGGTCATCGAGCCCATCTTTTCCTGCGAGGTGGCTTCGGCCTTGCGCACGGCGTCATTGAAGGCGGCGGCGACCAGGTCTTCCAGCAGGTC  >  AM260479/2594027‑2594252                                                                                                      |                                                                                                                             gcTGGGCGGACTTGGGCCCGACGCCGGGCAGCACCCGCAGCGCTTCGACCAGTGCCTGCAGCGAGGTCGGGGACGATGCCTTCACAGTGGCGTTCAGAACG‑GCAGCTTGAAGCCCGGGGGCAGCGGCAGGccc                                                                                              >  1:822157/1‑133 (MQ=255)         gACTTGGGCCCGACGCCGGGCAGCACCCGCAGCGCTTCGACCAGTGCCTGCAGCGAGGTCGGGGACGATGCCTTCACAGTGGCGTTCAGAACG‑GCAGCTTGAAGCCCGGGGGCAGCGGCAGGCCCGAGGTCATc                                                                                     >  2:1012569/1‑134 (MQ=255)         gACTTGGGCCCGACGCCGGGCAGCACCCGCAGCGCTTCGACCAGTGCCTGCAGCGAGGTCGGGGACGATGCCTTCACAGTGGCGTTCAGAACG‑GCAGCTTGAAGCCCGGGGGCAGCGGCAGGCCCGAGGTCATc                                                                                     >  2:935735/1‑134 (MQ=255)                   cGACGCCGGGCAGCACCCGCAGCGCTTCGACCAGTGCCTGCAGCGAGGTCGGGGACGATGCCTTCACAGTGGCGTTCAGAACG‑GCAGCTTGAAGCCCGGGGGCAGCGGCAGGCCCGAGGTCATCGAGCCCATCt                                                                           >  1:175958/1‑134 (MQ=255)                   cGACGCCGGGCAGCACCCGCAGCGCTTCGACCAGTGCCTGCAGCGAGGTCGGGGACGATGCCTTCACAGTGGCGTTCAGAACG‑GCAGCTTGAAGCCCGGGGGCAGCGGCAGGCCCGAGGTCATCGAGCCCATCt                                                                           >  1:224132/1‑134 (MQ=255)                    gACGCCGGGCAGCACCCGCAGCGCTTCGACCAGTGCCTGCAGCGAGGTCGGGGACGATGCCTTCACAGTGGCGTTCAGAACG‑GCAGCTTGAAGCCCGGGGGCAGCGGCAGGCCCGAGGTCATCGAGCCCATCtt                                                                          <  2:175958/134‑1 (MQ=255)                                              cGACCAGTGCCTGCAGCGAGGTCGGGGACGATGCCTTCACAGTGGCGTTCAGAACG‑GCAGCTTGAAGCCCGGGGGCAGCGGCAGGCCCGAGGTCATCGAGCCCATCTTTTCCTGCGAGGTGGCTTCGGCCTTg                                                 <  1:1012569/133‑1 (MQ=255)                                              cGACCAGTGCCTGCAGCGAGGTCGGGGACGATGCCTTCACAGTGGCGTTCAGAACG‑GCAGCTTGAAGCCCGGGGGCAGCGGCAGGCCCGAGGTCATCGAGCCCATCTTTTCCTGCGAGGTGGCTTCGGCCTTg                                                 <  2:224132/133‑1 (MQ=255)                                                                             tGCCTTCACAGTGGCGTTCAGAACG‑GCAGCTTGAAGCCCGGGGGCAGCGGCAGGCCCGAGGTCATCGAGCCCATCTTTTCCTGCGAGGTGGCTTCGGCCTTGCGCACGGCGTCATTGAAGGCGGCGGCGACCAg                 <  2:517106/134‑1 (MQ=255)                                                                                 ttCACAGTGGCGTTCAGAACG‑GCAGCTTGAAGCCCGGGGGCAGCGGCAGGCCCGAGGTCATCGAGCCCATCTTTTCCTGCGAGGTGGCTTCGGCCTTGCGCACGGCGTCATTGAAGGCGGCGGCGACCAGGTCt             >  1:842596/1‑134 (MQ=255)                                                                                    acaGTGGCGTTCAGAACG‑GCAGCTTGAAGCCCGGGGGCAGCGGCAGGCCCGAGGTCATCGAGCCCATCTTTTCCTGCGAGGTGGCTTCGGCCTTGTGCACGGCGTCATTGAAGGCGGCGGCGACCAGGTCTTc           <  2:842596/133‑1 (MQ=255)                                                                                           cGTTCAGAACG‑GCAGCTTGAAGCCCGGGGGCAGCGGCAGGCCCGAGGTCATCGAGCCCATCTTTTCCTGCGAGGTGGCTTCGGCCTTg                                                 >  1:537092/1‑88 (MQ=255)                                                                                           cGTTCAGAACG‑GCAGCTTGAAGCCCGGGGGCAGCGGCAGGCCCGAGGTCATCGAGCCCATCTTTTCCTGCGAGGTGGCTTCGGCCTTg                                                 <  2:537092/88‑1 (MQ=255)                                                                                           cGTTCAGAACG‑GCAGCTTGAAGCCCGGGGGCAGCGGCAGGCCCGAGGTCATCGAGCCCATCTTTTCCTGCGAGGTGGCTTCGGCCTTGCGCACGGCGTCATTGAAGGCGGCGGCGACCAGGTCTTCCAGCAGGt   >  1:946889/1‑134 (MQ=255)                                                                                            gTTCAGAACG‑GCAGCTTGAAGCCCGGGGGCAGCGGCAGGCCCGAGGTCATCGAGCCCATCTTTTCCTGCGAGGTGGCTTCGGCCTTGCGCACGGCGTCATTGAAGGCGGCGGCGACCAGGTCTTCCAGCAGGTc  >  1:437780/1‑134 (MQ=255)                                                                                                      |                                                                                                                             GCTGGGCGGACTTGGGCCCGACGCCGGGCAGCACCCGCAGCGCTTCGACCAGTGCTTGCAGCGAGGTCGGGGACGATGCCTTCACAGTGGCGTTCAGAACGCGCAGCTTGAAGCCCGGGGGCAGCGGCAGGCCCGAGGTCATCGAGCCCATCTTTTCCTGCGAGGTGGCTTCGGCCTTGCGCACGGCGTCATTGAAGGCGGCGGCGACCAGGTCTTCCAGCAGGTC  >  AM260479/2594027‑2594252

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