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%	93.0 / NA	26	G286G (GGT→GGG)	hemB	porphobilinogen synthase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base C (12/14);  total (12/14)

TCTGAGAGGAAACAAAATTAACGCAGAATCTTCTTCTCAGCCAAATCCAGCGCAAAGTAGCTGAAAATCAGATCCGCACCCGCACGCTTAATCGAACCTAAGCTTTCGAGCACGACTTTCTCTTCATCTATAGCACCCGCCAGCGCGGCGAACTTAATCATCGCATACTCACCGCTCACCTGATACGCGCCAATCGGCAATTCAGTACGTTCACGCAGCTCACGCACGATGTCGAGGTA  >  NZ_CP009273/384192‑384430                                                                                                                                       |                                                                                                         tCTGAGAGGAAACAAAATTAACGCAGAATCTTCTTCTCAGCCAAATCCAGCGCAAAGTAGCTGAAAATCAGATCCGCACCCGCACGCTTAATCGAACCTAAGCTTTCGAGCACGACTTTCTCTTCATCTATAGCCCCCg                                                                                                      >  1:141871/1‑139 (MQ=255)     agagGAAACAAAATTAACGCAGAATCTTCTTCTCAGCCAAATCCAGCGCAAAGTAGCTGAAAATCAGATCCGCACCCGCACGCTTAATCGAACCTAAGCTTTCGAGCACGACTTTCTCTTCATCTATAGCCCCCGCCAg                                                                                                  >  1:298428/1‑139 (MQ=255)           aaCAAAATTAACGCAGAATCTTCTTCTCAGCCAAATCCAGCGCAAAGTAGCTGAAAATCAGATCCGCACCCGCACGCTTAATCGAACCTAAGCTTTCGAGCACGACTTTCTCTTCATCTATAGCCCCCGCCAGCgcggc                                                                                            >  2:207504/1‑139 (MQ=255)            aCAAAATTAACGCAGAATCTTCTTCTCAGCCAAATCCAGCGCAAAGTAGCTGAAAATCAGATCCGCACCCGCACGCTTAATCGAACCTAAGCTTTCGAGCACGACTTTCTCTTCATCTATAGCCCCCGCCAGCgcggc                                                                                            >  1:12256/1‑138 (MQ=255)                    aaCGCAGAATCTTCTTCTCAGCCAAATCCAGCGCAAAGTAGCTGAAAATCAGATCCGCACCCGCACGCTTAATCGAACCTAAGCTTTCGAGCACGACTTTCTCTTCATCTATAGCCCCCGCCAGCGCGGCGAACTTaa                                                                                    <  2:298428/138‑1 (MQ=255)                     aCGCAGAATCTTCTTCTCAGCCAAATCCAGCGCAAAGTAGCTGAAAATCAGATCCGCACCCGCACGCTTAATCGAACCTAAGCTTTCGAGCACGACTTTCTCTTCATCTATAGCCCCCGCCAGCGCGGCGAACTTAatc                                                                                  <  2:486160/139‑1 (MQ=255)                           aaTCTTCTTCTCAGCCAAATCCAGCGCAAAGTAGCTGAAAATCAGATCCGCACCCGCACGCTTAATCGAACCTAAGCTTTCGAGCACGACTTTCTCTTCATCTATAGCCCCCGCCAGCGCGGCGAACTTAATCATCGCa                                                                            >  2:492111/1‑139 (MQ=255)                           aaTCTTCTTCTCAGCCAAATCCAGCGCAAAGTAGCTGAAAATCAGATCCGCACCCGCACGCTTAATCGAACCTAAGCTTTCGAGCACGACTTTCTCTTCATCTATAGCCCCCGCCAGCGCGGCGAACTTAATCATCGCa                                                                            <  2:141871/139‑1 (MQ=255)                           aaTCTTCTTCTCAGCCAAATCCAGCGCAAAGTAGCTGAAAATCAGATCCGCACCCGCACGCTTAATCGAACCTAAGCTTTCGAGCACGACTTTCTCTTCATCTATAGCCCCCGCCAGCGCGGCGAACTTAATCATCGCa                                                                            >  1:429027/1‑139 (MQ=255)                                 cttctCAGCCAAATCCAGCGCAAAGTAGCTGAAAATCAGATCCGCACCCGCACGCTTAATCGAACCTAAGCTTTCGAGCACGACTTTCTCTTCATCTATAGCCCCCGCCAGCGCGGCGAACTTAATCATCGCATAc                                                                         <  2:543394/136‑1 (MQ=255)                                 cttctCAGCCAAATCCAGCGCAAAGTAGCTGAAAATCAGATCCGCACCCGCACGCTTAATCGAACCTAAGCTTTCGAGCACGACTTTCTCTTCATCTATAGCCCCCGCCAGCGCGGCGAACTTAATCATCGCATACt                                                                        <  2:12256/137‑1 (MQ=255)                                       aGCCAAATCCAGCGCAAAGTAGCTGAAAATCAGATCCGCACCCGCACGCTTAATCGAACCTAAGCTTTCGAGCACGACTTTCTCTTCATCTATAGCCCCCGCCAGCGCGGCGAACTTAATCATCGCa                                                                            <  2:183730/127‑1 (MQ=255)                                       aGCCAAATCCAGCGCAAAGTAGCTGAAAATCAGATCCGCACCCGCACGCTTAATCGAACCTAAGCTTTCGAGCACGACTTTCTCTTCATCTATAGCCCCCGCCAGCGCGGCGAACTTAATCATCGCa                                                                            >  1:183730/1‑127 (MQ=255)                                       aGCCAAATCCAGCGCAAAGTAGCTGAAAATCAGATCCGCACCCGCACGCTTAATCGAACCTAAGCTTTCGAGCACGACTTTCTCTTCATCTATAGCCCCCGCCAGCGCGGCGAACTTAATCATCGCATACTCACCg                                                                   <  2:499062/136‑1 (MQ=255)                                                             cTGAAAATCAGATCCGCACCCGCACGCTTAATCGAACCTAAGCTTTCGAGCACGACTTTCTCTTCATCTATAGCCCCCGCCAGCGCGGCGAACTTAATCATCGCATACTCACCg                                                                   >  1:69027/1‑114 (MQ=255)                                                             cTGAAAATCAGATCCGCACCCGCACGCTTAATCGAACCTAAGCTTTCGAGCACGACTTTCTCTTCATCTATAGCCCCCGCCAGCGCGGCGAACTTAATCATCGCATACTCACCg                                                                   <  2:69027/114‑1 (MQ=255)                                                             cTGAAAATCAGATCCGCACCCGCACGCTTAATCGAACCTAAGCTTTCGAGCACGACTTTCTCTTCATCTATAGCCCCCGCCAGCGCGGCGAACTTAATCATCGCATACTCACCGCTCACCTGATACGCGCCAATCGGc                                           >  1:312175/1‑138 (MQ=255)                                                                   aTCAGATCCGCACCCGCACGCTTAATCGAACCTAAGCTTTCGAGCACGACTTTCTCTTCATCTATAGCCCCCGCCAGCGCGGCGAACTTAATCATCGCATACTCACCGCTCACCTGATACGCGCCAATCGGCAATTCAg                                    <  1:410611/139‑1 (MQ=255)                                                                              acccgcacGCTTAATCGAACCTAAGCTTTCGAGCACGACTTTCTCTTCATCTATAGCCCCCGCCAGCGCGGCGAACTTAATCATCGCATACTCACCGCTCACCTGATACGCGCCAATCGGCAATTCAGTACGTTCACGc                         <  1:207504/139‑1 (MQ=255)                                                                                            tCGAACCTAAGCTTTCGAGCACGACTTTCTCTTCATCTATAGCCCCCGCCAGCGCGGCGAACTTAATCATCGCATACTCACCGCTCACCTGATACGCGCCAATCGGCAATTCAGTACGTTCACGCAGCTCACGCACGAt           >  1:62510/1‑139 (MQ=255)                                                                                                aCCTAAGCTTTCGAGCACGACTTTCTCTTCATCTATAGCCCCCGCCAGCGCGGCGAACTTAatcat                                                                                <  1:392451/66‑1 (MQ=255)                                                                                                aCCTAAGCTTTCGAGCACGACTTTCTCTTCATCTATAGCCCCCGCCAGCGCGGCGAACTTAatcat                                                                                >  2:392451/1‑66 (MQ=255)                                                                                                     aGCTTTCGAGCACGACTTTCTCTTCATCTATAGCCCCCGCCAGCGCGGCGAACTTAATCATCGCATACTCACCGCTCACCTGATACGCGCCAATCGGCAATTCAGTTCGTTCACGCAGCTCACGCACGATGTCGAGGTa  <  1:492111/139‑1 (MQ=255)                                                                                                     aGCTTGCGAGCTGGCCTTTCTCTTCATCTATAGCCCCCGCCAGCGCGGCGAACTTAATCATCGCATACTCACCGCTCACCTGATACGCGCCAATCGGCAATTCAGTACGTTCACGCAGCCCACGCACGATGTCGAGGTa  <  2:151529/139‑1 (MQ=255)                                                                                                          tCGAGCACGACTTTCTCTTCATCTATAGCCCCCGCCAGCgcggcg                                                                                           <  2:171270/45‑1 (MQ=255)                                                                                                          tCGAGCACGACTTTCTCTTCATCTATAGCCCCCGCCAGCgcggcg                                                                                           >  1:171270/1‑45 (MQ=255)                                                                                                                                       |                                                                                                         TCTGAGAGGAAACAAAATTAACGCAGAATCTTCTTCTCAGCCAAATCCAGCGCAAAGTAGCTGAAAATCAGATCCGCACCCGCACGCTTAATCGAACCTAAGCTTTCGAGCACGACTTTCTCTTCATCTATAGCACCCGCCAGCGCGGCGAACTTAATCATCGCATACTCACCGCTCACCTGATACGCGCCAATCGGCAATTCAGTACGTTCACGCAGCTCACGCACGATGTCGAGGTA  >  NZ_CP009273/384192‑384430

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