Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NZ_CP009273	2,882,413	A→C	V349G (GTT→GGT)	cysI ←	assimilatory sulfite reductase (NADPH) hemoprotein subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NZ_CP009273	2,882,413	0	A	C	100.0%	80.0 / NA	24	V349G (GTT→GGT)	cysI	assimilatory sulfite reductase (NADPH) hemoprotein subunit
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base C (12/12);  total (12/12)

CCGGTTTTCAGCGGACGCGCCGGATAATCAAGGATGCGACCATTTTCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTAACCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGACGGATCGGTTCAAATTTGATCCCCGCGCGACGCTCCACTTCCGC  >  NZ_CP009273/2882326‑2882498                                                                                        |                                                                                      ccGGTTTTCAGCGGACGCGCCGGATAATCAAGGATGCGACCATTTTCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTAccc                                                                                     <  1:1048836/90‑1 (MQ=255) ccGGTTTTCAGCGGACGCGCCGGATAATCAAGGATGCGACCATTTTCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTAccc                                                                                     >  1:323154/1‑90 (MQ=255) ccGGTTTTCAGCGGACGCGCCGGATAATCAAGGATGCGACCATTTTCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTAccc                                                                                     >  2:61534/1‑90 (MQ=255)              gACGCGCCGGATAATCAAGGATGCGACCATTTTCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTACCCCAGCCAATACGAt                                                                        >  2:1043145/1‑90 (MQ=255)              gACGCGCCGGATAATCAAGGATGCGACCATTTTCGATAAACAGCGTCAGGGGCCAGTTATCATCAATGCCCTTACCCCAGCCAATACGAt                                                                        >  2:52262/1‑90 (MQ=255)                  cgcCGGATAATCAAGGATGCGACCATTTTCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTACCCCAGCCAATACGATCGcc                                                                    <  1:487213/90‑1 (MQ=255)                               aGGATGCGACCATTTTCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTGCCCCAGCCAATACGATCGCCTCGTCCGGTGAAc                                                       <  1:1118646/90‑1 (MQ=255)                               aGGATGCGACCATTTTCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTACCCCAGCCAATACGATCGCCTCGTCCGGTGAAc                                                       <  1:34913/90‑1 (MQ=255)                                 gATGCGACCATTTTAGATAAACAGCGTCAGGTGCCAGTTATCATCAATTCCCTTACCCCAGCCAATACGATCGCCTCGTCCGGTGAACTc                                                     >  2:943979/1‑90 (MQ=255)                                   tGCGACCATTTTCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTACCCCAGCCAATACGATCGCCTCGTCCGGTGAACTCat                                                   >  1:510579/1‑90 (MQ=255)                                    gCGACCATTTTCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTACCCCAGCCAATACGATCGCCTCGTCCGGTGAACTCata                                                  >  1:1126934/1‑90 (MQ=255)                                     cGACCATTTTCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTACCCCAGCCAATACGATCGCCTCGTCCGGTGAACTCatat                                                 >  2:125664/1‑90 (MQ=255)                                          aTTTTCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTACCCCAGCCAATACGATCGCCTCGTCCGGTGAAc                                                       <  2:916484/79‑1 (MQ=255)                                          aTTTTCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTACCCCAGCCAATACGATCGCCTCGTCCGGTGAAc                                                       >  1:916484/1‑79 (MQ=255)                                                 aTAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTACCCCAGCCAATACGATCGCCTCGTCCGGTGAACTCata                                                  >  1:316126/1‑77 (MQ=255)                                                 aTAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTACCCCAGCCAATACGATCGCCTCGTCCGGTGAACTCata                                                  <  2:316126/77‑1 (MQ=255)                                                     aCAGCGTCAGGTGCCAGTTATCATCAATGCCCTTACCCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGACGGATCGGTTCaa                                 <  1:478372/90‑1 (MQ=255)                                                                    aGTTATCATCAATGCCCTTACCCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGACGGATCGGTTCAAATTTGATCCCcgcgc                  <  2:998199/90‑1 (MQ=255)                                                                    aGTTATCATCAATGCCCTTACCCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGACGGATCGGTTCAAATTTGATCCCcgcgc                  <  2:510579/90‑1 (MQ=255)                                                                       tATCATCAATGCCCTTCCCCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGACGGATCGGTTCAAATTTGATCCCCGCGCGAc               <  1:515843/90‑1 (MQ=255)                                                                       tATCATCAATGCCCTTACCCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGACGGATCGGTTCAAATTTGATCCCCGCGCGAc               <  2:1126934/90‑1 (MQ=255)                                                                       tATCATCAATGCCCTTACCCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGACGGATCGGTTCAAATTTGATCCCCGCGCGAc               <  1:641337/90‑1 (MQ=255)                                                                            tcaATGCCCTTCCCCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGACGGATCGGTTCAAATTTGATCCCCGCGCGACGCTcc          <  1:747636/90‑1 (MQ=255)                                                                             caATGCCCTTACCCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGACGGATCGGTTCAAATTTGATCCCCGCGCGACGCTCCa         >  1:617042/1‑90 (MQ=255)                                                                               aTGCCCTTACCCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGACGGATCGGTTCAAATTTGATCCCCGCGCGACGCTCCACt       >  1:925115/1‑90 (MQ=255)                                                                                    cTTACCCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGACGGATCGGTTCAAATTTGATCCCCGCGCGACGCTCCACTTCCGc  >  2:961623/1‑90 (MQ=255)                                                                                    cTTACCCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGACGGATCGGTTCAAATTTGATCCCCGCGCGACGCTCCACTTCCGc  >  1:234440/1‑90 (MQ=255)                                                                                        |                                                                                      CCGGTTTTCAGCGGACGCGCCGGATAATCAAGGATGCGACCATTTTCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTAACCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGACGGATCGGTTCAAATTTGATCCCCGCGCGACGCTCCACTTCCGC  >  NZ_CP009273/2882326‑2882498

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