Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NZ_CP009273	384,326	A→C	G286G (GGT→GGG)	hemB ←	porphobilinogen synthase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NZ_CP009273	384,326	0	A	C	100.0%	70.1 / NA	20	G286G (GGT→GGG)	hemB	porphobilinogen synthase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base C (8/12);  total (8/12)

GAGGAAACAAAATTAACGCAGAATCTTCTTCTCAGCCAAATCCAGCGCAAAGTAGCTGAAAATCAGATCCGCACCCGCACGCTTAATCGAACCTAAGCTTTCGAGCACGACTTTCTCTTCATCTATAGCACCCGCCAGCGCGGCGAACTTAATCATCGCATACTCACCGCTCACCTGATACGCGCCAATCGGCAATTCAGTACGTTCACGCAGCTCACGCACGATGTCGAGGTACGCTCCAGCAGGTTTAACCATCAGGCAGTC  >  NZ_CP009273/384197‑384460                                                                                                                                  |                                                                                                                                       gagGAAACAAAATTAACGCAGAATCTTCTTCTCAGCCAAATCCAGCGCAAAGTAGCTGAAAATCAGATCCGCACCCGCACGCTTAATCGAACCTAAGCTTTCGAGCACGACTTTCTCTTCATCTATAGCCCCCGCCAgc                                                                                                                               >  1:246778/1‑139 (MQ=255)  agGAAACAAAATTAACGCAGAATCTTCTTCTCAGCCAAATCCAGCGCAAAGTAGCTGAAAATCAGATCCGCACCCGCACGCTTAATCGAACCTAAGCTTTCGAGCACGACTTTCTCTTCATCTATAGCCCCCGCCAgcg                                                                                                                              >  1:626835/1‑139 (MQ=255)       aCAAAATTAACGCAGAATCTTCTTCTCAGCCAAATCCAGCGCAAAGTAGCTGAAAATCAGATCCGCACCCGCACGCTTAATCGAACCTAAGCTTTCGAGCACGACTTTCTCTTCATCTATAGCCCCCGCCAGCgcggcg                                                                                                                         <  2:26364/139‑1 (MQ=255)                  gCAGAATCTTCTTCTCAGCCAAATCCAGCGCAAAGTAGCTGAAAATCAGATCCGCACCCGCACGCTTAATCGAACCTAAGCTTTCGAGCACGACTTTCTCTTCATCTATAGCCCCCGCCAGCGCGGCGAACTTAatcat                                                                                                              <  1:489165/139‑1 (MQ=255)                                  aGCCAAATCCAGCGCAAAGTAGCTGAAAATCAGATCCGCACCCGCACGCTTAATCGAACCTAAGCTTTCGAGCACGACTTTCTCTTCATCTATAGCCCCCGCCAGCGCGGCGAACTTAATCATCGCATACTCACCg                                                                                                 <  2:311637/136‑1 (MQ=255)                                  aGCCAAATCCAGCGCAAAGTAGCTGAAAATCAGATCCGCACCCGCACGCTTAATCGAACCTAAGCTTTCGAGCACGACTTTCTCTTCATCTATAGCCCCCGCCAGCGCGGCGAACTTAATCATCGCATACTCACCg                                                                                                 <  2:192853/136‑1 (MQ=255)                                                                             gcacGCTTAATCGAACCTAAGCTTTCGAGCACGACTTTCTCTTCATCTATAGCCCCCGCCAGCGCGGCGAACTTAATCATCGCATACTCACCGCTCACCTGATACGCGCCAATCGGCAATTCAGTACGTTCACGCAGCt                                                   >  2:179124/1‑139 (MQ=255)                                                                                     aaTCGAACCTAAGCTTTCGAGCACGACTTTCTCTTCATCTATAGCCCCCGCCAGCGCGGCGAACTTAATCATCGCATACTCACCGCTCACCTGATACGCGCCAATCGGCAATTCAGTACGTTCACGCAGCTcacgca                                             <  2:457287/137‑1 (MQ=255)                                                                                        cGAACCTAAGCTTTCGACCACGACTTTCTCTTCATCTATAGCCCCCGCCAGCGCGGCGAACTTAATCATCGCATACTCACCGCTCACCTGATACGCGCCAATCGGCAATTCAGTACGTTCACGCAGCTCACGCACGATg                                        <  2:354030/139‑1 (MQ=255)                                                                                           aCCTAAGCTTTCGAGCACGACTTTCTCTTCATCTATAGCCCCCGCCAGCGCGGCGAACTTAATCATCGCATACTCACCg                                                                                                 >  1:281513/1‑79 (MQ=255)                                                                                           aCCTAAGCTTTCGAGCACGACTTTCTCTTCATCTATAGCCCCCGCCAGCGCGGCGAACTTAATCATCGCATACTCACCg                                                                                                 <  2:281513/79‑1 (MQ=255)                                                                                           aCCTAAGCTTTCGAGCACGACTTTCTCTTCATCTATAGCCCCCGCCAGCGCGGCGAACTTAATCATCGCATACTCACCGCTCACCTGATACGCGCCAATCgg                                                                          <  1:117628/102‑1 (MQ=255)                                                                                           aCCTAAGCTTTCGAGCACGACTTTCTCTTCATCTATAGCCCCCGCCAGCGCGGCGAACTTAATCATCGCATACTCACCGCTCACCTGATACGCGCCAATCgg                                                                          >  2:117628/1‑102 (MQ=255)                                                                                                aGCTTTCGAGCACGACTTTCTCTTCATCTATAGCCCCCGCCAGCGCGGCGAACTTAATCATCGCATACTCACCGCTCACCTGATACGCGCCAATCGGCAATTCAGTACGTTCACGCAGCTCACGCACGATGTCGAGGTa                                <  1:179124/139‑1 (MQ=255)                                                                                                        aGCACGACTTTCTCTTCATCTATAGCCCCCGCCAGCGCGGCGAACTTAATCATCGCATACTCACCGCTCACCTGATACGCGCCAATCGGCAATTCAGTACGTTCACGCAGCTCACGCACGATGTCGAGGTACGCTCcag                        <  1:193849/139‑1 (MQ=255)                                                                                                               ctctctcTTCATCTATAGCCCCCGCCAGCGCGGCGAACTTAATCATCGCATACTCACCGCTCACCTGATACGCGCCAATCGGCAATTCAGTACGTTCACGCAGCTCACGCACGATGTCGAGGTACGCTCCAGCAGGttt                 <  2:626835/136‑1 (MQ=255)                                                                                                                         aTCTATAGCCCCCGCCAGCGCGGCGAACTTAATCATCGCATACTCACCGCTCACCTGATACGCGCCAATCGGCAATTCAGTACGTTCACg                                                        >  1:636673/1‑90 (MQ=255)                                                                                                                         aTCTATAGCCCCCGCCAGCGCGGCGAACTTAATCATCGCATACTCACCGCTCACCTGATACGCGCCAATCGGCAATTCAGTACGTTCACg                                                        <  2:636673/90‑1 (MQ=255)                                                                                                                         aTCTATAGCCCCCGCCAGCGCGGCGAACTTAATCATCGCATACTCACCGCTCACCTGATACGCGCCAATCGGCAATTCAGTACGTTCACGCAGCTCACGCACGATGTCGAGGTACGCTCCAGCAGGTTTAACCATCAgg       >  1:369087/1‑139 (MQ=255)                                                                                                                              taGCCCCCGCCAGCGCGGCGAACTTAATCATCGCATACTCACCGCTCACCTGATACGCGCCAATCGGCAATTCAGTACGTTCACGCAGCTCACGCACGATGTCGAGGTACGCTCCAGCAGGTTTAACCATCAGGCAGTc  >  1:451631/1‑139 (MQ=255)                                                                                                                                  |                                                                                                                                       GAGGAAACAAAATTAACGCAGAATCTTCTTCTCAGCCAAATCCAGCGCAAAGTAGCTGAAAATCAGATCCGCACCCGCACGCTTAATCGAACCTAAGCTTTCGAGCACGACTTTCTCTTCATCTATAGCACCCGCCAGCGCGGCGAACTTAATCATCGCATACTCACCGCTCACCTGATACGCGCCAATCGGCAATTCAGTACGTTCACGCAGCTCACGCACGATGTCGAGGTACGCTCCAGCAGGTTTAACCATCAGGCAGTC  >  NZ_CP009273/384197‑384460

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