Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	AM260479	3,301,576	+G	intergenic (+45/‑146)	lspA → / → dfp	lipoprotein signal peptidase II/phosphopantothenoylcysteine synthetase/decarboxylase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	AM260479	3,301,574	1	.	G	100.0%	46.0 / NA	13	intergenic (+43/‑148)	lspA/dfp	lipoprotein signal peptidase II/phosphopantothenoylcysteine synthetase/decarboxylase
Reads supporting (aligned to +/- strand):  ref base . (0/0);  new base G (4/9);  total (4/9)

CCGGCCTTCAACGTGGCCGACTGCGCGATCACGGTCGGCGCGGTGCTGCTGATCGTGGACGAGCTGCGCAGGGTGCGGCGCCACTGAGCGCCCTGATTGTGTTCCGCTCTCCCCTCACGGGGAGAGGCCA‑GGTGAGGGGCGGTTCAGCAAGGCACCACGCCACCAGAAC  >  AM260479/3301445‑3301613                                                                                                                                   |                                        ccggccTTCAACGTGGCCGACTGCGCGATCACGGTCGGCGCGGTGCTGCTGATCGTGGACGAGCTGCGCAGGGTGCGGCGCCACTGAGCGCCCTGATTGTGTTCCGCTCTCCCCTCACGGGGAGAGGCCA‑gggt                                     <  1:979381/134‑3 (MQ=255)   ggccTTCAACGTGGCCGACTGCGCGATCACGGTCGGCGCGGTGCTGCTGATCGTGGACGAGCTGCGCAGGGTGCGGCGCCACTGAGCGCCCTGATTGTGTTCCGCTCTCCCCTCACGGGGAGAGGCCA‑gggtga                                   <  1:826422/134‑5 (MQ=255)                        cgcgATCACGGTCGGCGCGGTGCTGCTGATCGTGGACGAGCTGCGCAGGGTGCGGCGCCACTGAGCGCCCTGATTGTGTTCCGCTCTCCCCTCACGGGGAGAGGCCAGGGTGAGGGGCGGTTCAGCAAGGca                 <  2:112311/132‑1 (MQ=255)                        cgcgATCACGGTCGGCGCGGTGCTGCTGATCGTGGACGAGCTGCGCAGGGTGCGGCGCCACTGAGCGCCCTGATTGTGTTCCGCTCTCCCCTCACGGGGAGAGGCCAGGGTGAGGGGCGGTTCAGCAAGGca                 <  2:493544/132‑1 (MQ=255)                                  gTCGGCGCGGTGCTGCTGATCGTGGACGAGCTGCGCAGGGTGCGGCGCCACTGAGCGCCCTGATTGTGTTCCGCTCTCCCCTCACGGGGAGAGGCCAGGGTGAGGGGCGGTTCAGCAAGGCACCACGCCACCAg     <  1:720292/134‑1 (MQ=255)                                  gTCGGCGCGGTGCTGCTGATCGTGGACGAGCTGCGCAGGGTGCGGCGCCACTGAGCGCCCTGATTGTGTTCCGCTCTCCCCTCACGGGGAGAGGCCAGGGTGAGGGGCGGTTCAGCAAGGCACCACGCCACCAg     <  2:125881/134‑1 (MQ=255)                                     ggCGCGGTGCTGCTGATCGTGGACGAGCTGCGCAGGGTGCGGCGCCACTGAGCGCCCTGATTGTGTTCCGCTCTCCCCTCACGGGGAGAGGCCAGGGTGAGGGGCGGTTCAGCAAGGCACCACGCCACCAGAAc  <  2:407234/134‑1 (MQ=255)                                      gcgcGGTGCTGCTGATCGTGGACGAGCTGCGCAGGGTGCGGCGCCACTGAGCGCCCTGATTGTGTTCCGCTCTCCCCTCACGGGGAGAGGCCAGGGTGAGGGGc                               >  1:1067866/1‑104 (MQ=255)                                      gcgcGGTGCTGCTGATCGTGGACGAGCTGCGCAGGGTGCGGCGCCACTGAGCGCCCTGATTGTGTTCCGCTCTCCCCTCACGGGGAGAGGCCAGGGTGAGGGGc                               <  2:1067866/104‑1 (MQ=255)                                           gTGCTGCTGATCGTGGACGAGCTGCGCAGGGTGCGGCGCCACTGAGCGCCCTGATTGTGTTCCGCTCTCCCCTCACGGGGAGAGGCCAGGGTGAggg                                 >  1:452767/1‑97 (MQ=255)                                           gTGCTGCTGATCGTGGACGAGCTGCGCAGGGTGCGGCGCCACTGAGCGCCCTGATTGTGTTCCGCTCTCCCCTCACGGGGAGAGGCCAGGGTGAggg                                 <  2:452767/97‑1 (MQ=255)                                           gTGCTGCTGATCGTGGACGAGCTGCACAGGGTGCGGCGCCACTGAGCGCCCTGATTGTGTTCCGCTCTCCCCTCACGGGGAGAGGCCAGGGTGAggg                                 >  1:125783/1‑97 (MQ=255)                                           gTGCTGCTGATCGTGGACGAGCTGCACAGGGTGCGGCGCCACTGAGCGCCCTGATTGTGTTCCGCTCTCCCCTCACGGGGAGAGGCCAGGGTGAggg                                 <  2:125783/97‑1 (MQ=255)                                                                                 ccACTGAGCGCCCTGATTGTGTTCCGCTCTCCCCTCACGGGGAGAGGCCAGGGTGAGGGGCGGTTCAGCAAGGca                 >  1:1118487/1‑75 (MQ=255)                                                                                 ccACTGAGCGCCCTGATTGTGTTCCGCTCTCCCCTCACGGGGAGAGGCCAGGGTGAGGGGCGGTTCAGCAAGGca                 <  2:1118487/75‑1 (MQ=255)                                                                                                                                   |                                        CCGGCCTTCAACGTGGCCGACTGCGCGATCACGGTCGGCGCGGTGCTGCTGATCGTGGACGAGCTGCGCAGGGTGCGGCGCCACTGAGCGCCCTGATTGTGTTCCGCTCTCCCCTCACGGGGAGAGGCCA‑GGTGAGGGGCGGTTCAGCAAGGCACCACGCCACCAGAAC  >  AM260479/3301445‑3301613

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