Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NZ_CP009273	3,349,491	T→C	S470P (TCC→CCC)	gltB →	glutamate synthase large subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NZ_CP009273	3,349,491	0	T	C	100.0%	90.2 / NA	25	S470P (TCC→CCC)	gltB	glutamate synthase large subunit
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (15/10);  total (15/10)

ATGAAGAAGTGGGTAGCCGCGAACTGGACGACGACACGCTTGCCAGCTACCAGAAACA‑GTTTAACTACAGCGCGGAAGAGCTGGACTCCGTAATTCGCGTACTGGGCGAAAACGGTCAGGAAGCGGTCGGTTCGATGGGCGATGATACCCCATTCGCCG  >  NZ_CP009273/3349405‑3349563                                                                                        |                                                                         aTGAAGAAGTGGGTAGCCGCGAACTGGACGACGACACGCTTGCCAGCTACCAGAAACA‑GTTTAACTACAGCGCGGAAGAGCTGGACCCCGt                                                                      >  2:845670/1‑91 (MQ=255) aTGAAGAAGTGGGTAGCCGCGAACTGGACGACGACACGCTTGCCAGCTACCAGAAACA‑GTTTAACTACAGCGCGGAAGAGCTGGACCCCGt                                                                      >  2:463530/1‑91 (MQ=255)           gggTAGCCGCGAACTGGACGACGACACGCTTGCCAGCTACCAGAAACA‑GTTTAACTACAGCGCGGAAGAGCTGGACCCCGTa                                                                     <  1:27902/82‑1 (MQ=255)           gggTAGCCGCGAACTGGACGACGACACGCTTGCCAGCTACCAGAAACA‑GTTTAACTACAGCGCGGAAGAGCTGGACCCCGTa                                                                     >  2:27902/1‑82 (MQ=255)           gggTAGCCGCGAACTGGACGACGACACGCTTGCCAGCTACCAGAAACA‑GTTTAACTACAGCGCGGAAGAGCTGGACCCCGTAATTCGCGTa                                                            >  2:5917/1‑91 (MQ=255)            ggTAGCCGCGAACTGGACGACGACACGCTTGCCAGCTACCAGAAACA‑GTTTAACTACAGCGCGGAAGAGCTGGACCCCGTAATTCGCGTAc                                                           <  1:723563/91‑1 (MQ=255)             gTAGCCGCGAGCTGGACGACGACACGCTTGCCAGCTACCAGAAACA‑GTTTAACTACAGCGCGGAAGAGCTGGACCCCGTAATTCGCGTACt                                                          >  2:18117/1‑91 (MQ=255)                gCCGCGAACTGGACGACGACACGCTTGCCAGCTACCAGAAACAGGTTTAACTACAGCGCGGAAGAGCTGGACCCCGTAATTCGCGTACTgg                                                        >  2:125882/1‑91 (MQ=255)                       aCTGGACGACGACACGCTTGCCAGCTACCAGAAACA‑GTTTAACTACAGCGCGGAAGAGCTGGACCCCGTAATTCGCGTACTGGGCGAAAAc                                                <  2:467155/91‑1 (MQ=255)                       aCTGGACGACGACACGCTTGCCAGCTACCAGAAACA‑GTTTAACTACAGCGCGGAAGAGCTGGACCCCGTAATTCGCGTACTGGGCGAAAAc                                                <  2:177924/91‑1 (MQ=255)                       aCTGGACGACGACACGCTTGCCAGCTACCAGAAACA‑GTTTAACTACAGCGCGGAAGAGCTGGACCCCGTAATTCGCGTACTGGGCGAAAAc                                                <  2:211731/91‑1 (MQ=255)                             cgacgacACGCTTGCCAGCTACCAGAAACA‑GTTTAACTACAGCGCGGAAGAGCTGGACCCCGTAATTCGCGTACTGGGCGAAAACGGTCAg                                          >  1:560834/1‑91 (MQ=255)                                       cTTGCCAGCTACCAGAAACA‑GTTTAACTACAGCGCGGAAGAGCTGGACCCCGTAATTCGCGTACTGGGCGAAAACGGTCAGGAAGcggtcg                                >  2:350012/1‑91 (MQ=255)                                         tGCCAGCTACCAGAAACA‑GTTTAACTACAGCGCGGAAGAGCTGGACCCCGTAATTCGCGTACTGGGCGAAAACGGTCAGGAAGcggtcggt                              >  1:613418/1‑91 (MQ=255)                                                  ccAGAAACA‑GTTTAACTACAGCGCGGAAGAGCTGGACCCCGTAATTCGCGTACTGGGCGAAAACGGTCAGGAAGCGGTCGGTTCGATGGGc                     <  1:463530/91‑1 (MQ=255)                                                  ccAGAAACA‑GTTTAACTACAGCGCGGAAGAGCTGGACCCCGTAATTCGCGTACTGGGCGAAAACGGTCAGGAAGCGGTCGGTTCGATGGGc                     <  2:792077/91‑1 (MQ=255)                                                  ccAGAAACA‑GTTTAACTACAGCGCGGAAGAGCTGGACCCCGTAATTCGCGTACTGGGCGAAAACGGTCAGGAAGCGGTCGGTTCGATGGGc                     <  1:218441/91‑1 (MQ=255)                                                            gTTTAACTACAGCGCGGAAGAGCTGGACCCCGTAATTCGCGTACTGGGCGAAAACGGTCAGGAAGCGGTCGGTTCGATGGGCGATGATAcc            >  1:418739/1‑91 (MQ=255)                                                                      aGCGCGGAAGAGCTGGACCCCGTAATTCGCGTACTGGGCGAAAACGGTCAGGaa                                       >  1:729196/1‑54 (MQ=255)                                                                      aGCGCGGAAGAGCTGGACCCCGTAATTCGCGTACTGGGCGAAAACGGTCAGGaa                                       <  2:729196/54‑1 (MQ=255)                                                                      aGCGCGGAAGAGCTGGACCCCGTAATTCGCGTACTGGGCGAAAACGGTCAGGAAGCGGTCGGTTCGATGGGCgatga                <  1:815699/77‑1 (MQ=255)                                                                      aGCGCGGAAGAGCTGGACCCCGTAATTCGCGTACTGGGCGAAAACGGTCAGGAAGCGGTCGGTTCGATGGGCgatga                >  2:815699/1‑77 (MQ=255)                                                                      aGCGCGGAAGAGCTGGACCCCGTAATTCGCGTACTGGGCGAAAACGGTCAGGAAGCGGTCGGTTCGATGGGCGATGATACCCCATTCGCCg  >  2:205378/1‑91 (MQ=255)                                                                      aGCGCGGAAGAGCTGGACCCCGTAATTCGCGTACTGGGCGAAAACGGTCAGGAAGCGGTCGGTTCGATGGGCGATGATACCCCATTCGCCg  >  2:219630/1‑91 (MQ=255)                                                                      aGCGCGGAAGAGCTGGACCCCGTAATTCGCGTACTGGGCGAAAACGGTCAGGAAGCGGTCGGTTCGATGGGCGATGATACCCCATTCGCCg  >  1:792648/1‑91 (MQ=255)                                                                                        |                                                                         ATGAAGAAGTGGGTAGCCGCGAACTGGACGACGACACGCTTGCCAGCTACCAGAAACA‑GTTTAACTACAGCGCGGAAGAGCTGGACTCCGTAATTCGCGTACTGGGCGAAAACGGTCAGGAAGCGGTCGGTTCGATGGGCGATGATACCCCATTCGCCG  >  NZ_CP009273/3349405‑3349563

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑