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

AGAGGAAACAAAATTAACGCAGAATCTTCTTCTCAGCCAAATCCAGCGCAAAGTAGCTGAAAATCAGATCCGCACCCGCACGCTTAATCGAACCTAAGCTTTCGAGCACGACTTTCTCTTCATCTATAGCACCCGCCAGCGCGGCGAACTTAATCATCGCATACTCACCGCTCACCTGATACGCGCCAATCGGCAATTCAGTACGTTCACGCAGCTCACGCACGATGTCGAGGTACGCTCCAGCAGGTTTAACCATCAGG  >  NZ_CP009273/384196‑384455                                                                                                                                   |                                                                                                                                  agagGAAACAAAATTAACGCAGAATCTTCTTCTCAGCCAAATCCAGCGCAAAGTAGCTGAAAATCAGATCCGCACCCGCACGCTTAATCGAACCTAAGCTTTCGAGCACGACTTTCTCTTCATCTATAGCCCCCGCCAg                                                                                                                           >  2:424618/1‑139 (MQ=255)        aCAAAATTAACGCAGAATCTTCTTCTCAGCCAAATCCAGCGCAAAGTAGCTGAAAATCAGATCCGCACCCGCACGCTTAATCGAACCTAAGCTTTCGAGCACGACTTTCTCTTCATCTATAGCCCCCGCCAGCgcggcg                                                                                                                    <  2:62024/139‑1 (MQ=255)        aCAAAATTAACGCAGAATCTTCTTCTCAGCCAAATCCAGCGCAAAGTAGCTGAAAATCAGATCCGCACCCGCACGCTTAATCGAACCTAAGCTTTCGAGCACGACTTTCTCTTCATCTATAGCCCCCGCCAGCgcggcg                                                                                                                    <  2:530726/139‑1 (MQ=255)        aCAAAATTAACGCAGAATCTTCTTCTCAGCCAAATCCAGCGCAAAGTAGCTGAAAATCAGATCCGCACCCGCACGCTTAATCGAACCTAAGCTTTCGAGCACGACTTTCTCTTCATCTATAGCCCCCGCCAGCgcggcg                                                                                                                    <  1:678967/139‑1 (MQ=255)          aaaaTTAACGCAGAATCTTCTTCTCAGCCAAATCCAGCGCAAAGTAGCTGAAAATCAGATCCGCACCCGCACGCTTAATCGAACCTAAGCTTTCGAGCACGACTTTCTCTTCATCTATAGCCCCCGCCAGCGCGGCGaa                                                                                                                  <  1:32973/139‑1 (MQ=255)                 aCGCAGAATCTTCTTCTCAGCCAAATCCAGCGCAAAGTAGCTGAAAATCAGATCCGCACCCGCACGCTTAATCGAACCTAAGCTTTCGAGCACGACTTTCTCTTCATCTATAGCCCCCGCCAGCGCGGCGAACTTAat                                                                                                            <  2:209230/138‑1 (MQ=255)                   gCAGAATCTTCTTCTCAGCCAAATCCAGCGCAAAGTAGCTGAAAATCAGATCCGCACCCGCACGCTTAATCGAACCTAAGCTTTCGAGCACGACTTTCTCTTCATCTATAGCCCCCGCCAGCGCGGCGAACTTAatcat                                                                                                         <  1:424618/139‑1 (MQ=255)                      gAATCTTCTTCTCAGCCAAATCCAGCGCAAAGTAGCTGAAAATCAGATCCGCACCCGCACGCTTAATCGAACCTAAGCTTTCGAGCACGACTTTCTCTTCATCTATAGCCCCCGCCAGCGCGGCGAACTTAATCATCGc                                                                                                      <  1:708333/139‑1 (MQ=255)                            tcttctCAGCCAAATCCGGCGCAAAGTAGCTGAAAATCAGATCCGCACCCGCACGCTTAATCGAACCTAAGCTTTCGAGCACGACTTTCTCTTCATCTATAGCCCCCGCCAGCGCGGCGAACTTAATCATCGCATACt                                                                                                 >  1:452793/1‑138 (MQ=255)                             cttctCAGCCAAATCCAGCGCAAAGTAGCTGAAAATCAGATCCGCACCCGCACGCTTAATCGAACCTAAGCTTTCGAGCACGACTTTCTCTTCATCTATAGCCCCCGCCAGCGCGGCGAACTTAATCATCGCATACt                                                                                                 <  1:288714/137‑1 (MQ=255)                               tctcAGCCAAATCCAGCGCAAAGTAGCTGAAAATCAGATCCGCACCCGCACGCTTAATCGAACCTAAGCTTTCGAGCACGACTTTCTCTTCATCTATAGCCCCCGCCAGCGCGGCGAACTTAATCATCGCATACTCAcc                                                                                             >  1:377122/1‑139 (MQ=255)                                   aGCCAAATCCAGCGCAAAGTAGCTGAAAATCAGATCCGCACCCGCACGCTTAATCGAACCTAAGCTTTCGAGCACGACTTTCTCTTCATCTATAGCCCCCGCCAGCGCGGCGAACTTAATCATCg                                                                                                       >  2:40394/1‑125 (MQ=255)                                   aGCCAAATCCAGCGCAAAGTAGCTGAAAATCAGATCCGCACCCGCACGCTTAATCGAACCTAAGCTTTCGAGCACGACTTTCTCTTCATCTATAGCCCCCGCCAGCGCGGCGAACTTAATCATCg                                                                                                       <  1:40394/125‑1 (MQ=255)                                   aGCCAAATCCAGCGCAAAGTAGCTGAAAATCAGATCCGCACCCGCACGCTTAATCGAACCTAAGCTTTCGAGCACGACTTTCTCTTCATCTATAGCCCCCGCCAGCGCGGCGAACTTAATCATCGCATACTCACCg                                                                                            <  2:691792/136‑1 (MQ=255)                                   aGCCAAATCCAGCGCAAAGTAGCTGAAAATCAGATCCGCACCCGCACGCTTAATCGAACCTAAGCTTTCGAGCACGACTTTCTCTTCATCTATAGCCCCCGCCAGCGCGGCGAACTTAATCATCGCATACTCACCg                                                                                            <  1:336832/136‑1 (MQ=255)                                   aGCCAAATCCAGCGCAAAGTAGCTGAAAATCAGATCCCCCCCCGCACGCTTAACCGAACCTAAGCTTTCGAGCACGACTTTCTCTTCATCTATAGCCCCCGCCAGCGCGGCGAACTTAATCATCGCATACTCACCGCt                                                                                          <  2:508622/138‑1 (MQ=255)                                     ccAAATCCAGCGCAAAGTAGCTGAAAATCAGATCCGCACCCGCACGCTTAATCGAACCTAAGCTTTCGAGCACGACTTTCTCTTCATCTATAGCCCCCGCCAGCGCGGCGAACTTAATCATCGCATACTCACCGCTCAc                                                                                       >  2:270235/1‑139 (MQ=255)                                     ccAAATCCAGCGCAAAGTAGCTGAAAATCAGATCCGCACCCGCACGCTTAATCGAACCTAAGCTTTCGAGCACGACTTTCTCTTCATCTATAGCCCCCGCCAGCGCGGCGAACTTAATCATCGCATACTCACCGCTCAc                                                                                       <  2:148880/139‑1 (MQ=255)                                      cAAATCCAGCGCAAAGTAGCTGAAAATCAGATCCGCACCCGCACGCTTAATCGAACCTAAGCTTTCGAGCACGACTTTCTCTTCATCTATAGCCCCCGCCAGCGCGGCGAACTTAATCATCGCATACTCACCGCTCAcc                                                                                      <  1:393990/139‑1 (MQ=255)                                                  aaaGTAGCTGAAAATCAGATCCGCACCCGCACGCTTAATCGAACCTAAGCTTTCGAGCACGACTTTCTCTTCATCTATAGCCCCCGCCAGCGCGGCGAACTTAATCATCGCATACTCACCGCTCACCTGATACGCGCCa                                                                          <  2:494488/139‑1 (MQ=255)                                                        gCTGAAAATCAGATCCGCACCCGCACGCTTAATCGAACCTAAGCTTTCGAGCACGACTTTCTCTTCATCTATAGCCCCCGCCAGCGCGGCGAACTTAATCATCGCATACTCACCGCTCACCTGATACGCGCCAATCgg                                                                     <  2:67080/138‑1 (MQ=255)                                                        gCTGAAAATCAGATCCGCACCCGCACGCTTAATCGAACCTAAGCTTTCGAGCACGACTTTCTCTTCATCTATAGCCCCCGCCAGCGCGGCGAACTTAATCATCGCATACTCACCGCTCACCTGATACGCGCCAATCgg                                                                     <  2:377122/138‑1 (MQ=255)                                                        gCTGAAAATCAGATCCGCACCCGCACGCTTAATCGAACCTAAGCTTTCGAGCACGACTTTCTCCTCATCTATAGCCACCGCCAGCGCGGCGAACTTAATCATCGCATACTCACCGCTCACCTGATACGCGCTAATCgg                                                                     <  1:687724/138‑1 (MQ=255)                                                                         cacccgcacGCTTAATCGAACCTAAGCTTTCGAGCACGACTTTCTCTTCATCTATAGCCCCCGCCAGCgcgg                                                                                                                      <  1:757856/72‑1 (MQ=255)                                                                         cacccgcacGCTTAATCGAACCTAAGCTTTCGAGCACGACTTTCTCTTCATCTATAGCCCCCGCCAGCgcgg                                                                                                                      >  2:757856/1‑72 (MQ=255)                                                                          acccgcacGCTTAATCGAACCTAAGCTTTCGAGCACGACTTTCTCTTCATCTATAGCCCCCGCCAGCGCGGCGAACTTAATCATCGCATACTCACCGCTCACCTGATACGCGCCAATCGGCAATTCAGTACGTTCACGc                                                  <  1:601811/139‑1 (MQ=255)                                                                          acccgcacGCTTAATCGAACCTAAGCTTTCGAGCACGACTTTCTCTTCATCTATAGCCCCCGCCAGCGCGGCGAACTTAATCATCGCATACTCACCGCTCACCTGATACGCGCCAATCGGCAATTCAGTACGTTCACGc                                                  <  2:82638/139‑1 (MQ=255)                                                                              gcacGCTTAATCGAACCTAAGCTTTCGAGCACGACTTTCTCTTCATCTATAGCCCCCGCCAGCGCGGCGAACTTAATCATCGCATACTCACCGCTCACCTGATACGCGCCAATCGGCAATTCAGTACGTTCACGCAGCt                                              >  1:741996/1‑139 (MQ=255)                                                                               cacGCTTAATCGAACCTAAGCTTTCGAGCACGACTTTCTCTTCATCTATAGCCCCCGCCAGCGCGGCGAACTTAATCATCGCATACt                                                                                                 <  2:98552/87‑1 (MQ=255)                                                                               cacGCTTAATCGAACCTAAGCTTTCGAGCACGACTTTCTCTTCATCTATAGCCCCCGCCAGCGCGGCGAACTTAATCATCGCATACt                                                                                                 >  1:98552/1‑87 (MQ=255)                                                                                      aaTCGATCCTAAGCTTTCGAGAACGACTATCACTACATCTAGAGACCCCGCCAGCGCGGCGAACGTAAGCATCGCATACTCACCGCTCACGTGATACGCGCCAATCGGCAATTCAGTACGTTCACGCAGCTcacgcacg                                      <  2:452793/139‑1 (MQ=255)                                                                                      aaTCGAACCTAAGCTTTCGAGCACGACTTTCTCTTCATCTATAGCCCCCGCCAGCGCGGCGAACTTAATCATCGCATACTCACCGCTCACCTGATACGCGCCAATCGGCAATTCAGTACGTTCACGCAGCTcacgcac                                       >  1:385903/1‑138 (MQ=255)                                                                                      aaTCGAACCTAAGCTTTCGAGCACGACTTTCTCTTCATCTATAGCCCCCGCCAGCGCGGCGAACTTAATCATCGCATACTCACCGCTCACCTGATACGCGCCAATCGGCAATTCAGTACGTTCACGCAGCTcacgca                                        <  2:352009/137‑1 (MQ=255)                                                                                              cTAAGCTTTCGAGCACGACTTTCTCTTCATCTATAGCCCCCGCCAGCGCGGCGAACTTAATCATCGCATACTCACCGCTCACCTGATACGCGCCAATCGGCAATTCAGTACGTTCACGCAGCTCACGCACGATGTCGAg                              >  2:745042/1‑139 (MQ=255)                                                                                                 aGCTTTCGAGCACGACTTTCTCTTCATCTATAGCCCCCGCCAGCGCGGCGAACTTAATCATCGCATACTCACCGCTCACCTGATACGCGCCAATCGGCAATTCAGTACGTTCACGCAGCTCACGCACGATGTCGAgg                             <  2:613817/137‑1 (MQ=255)                                                                                                      tCGAGCACGACTTTCTCTTCATCTATAGCCCCCGCCAGCGCGGCGAACTTAATCATCGCATACTCACCGCTCACCTGATACGCGCCAATCGGCAATTCAGTACGTTCACGCAGCTCACGCACGATGTCGAGGTACGCTc                      >  2:479042/1‑139 (MQ=255)                                                                                                      tCGAGCACGACTTTCTCTTCATCTATAGCCCCCGCCAGCGCGGCGAACTTAATCATCGCATACTCACCGCTCACCTGATACGCGCCAATCGGCAATTCAGTACGTTCACGCAGCTCACGCACGATGTCGAGGTACGCTc                      >  1:681144/1‑139 (MQ=255)                                                                                                         aGCACGACTTTCTCTTCATCTATAGCCCCCGCCAGCGCGGCGAACTTAATCATCGCATACTCACCGCTCACCTGATACGCGCCAATCGGCAATTCAGTACGTTCACGCAGCTCACGCACGATGTCGAGGTACGCTCcag                   <  1:270235/139‑1 (MQ=255)                                                                                                          gCACGACTTTCTCTTCATCTATAGCCCCCGCCAGCGCGGCGAACTTAATCATCGCATACTCACCGCTCACCTGATACGCGCCAATCGGCAATTCAGTACGTTCACGCAGCTCACGCACGATGTCGAGGTACGCTCcagc                  <  1:745042/139‑1 (MQ=255)                                                                                                           cACGACTTTCTCTTCATCTATAGCCCCCGCCAGCGCGGCGAACTTAATCATCGCATACTCACCGCTCACCTGATACGCGCCAATCGGCAATTCAGTACGTTCACGCAGCTCACGCACGATGTCGAGGTACGCTCcagca                 <  1:6752/139‑1 (MQ=255)                                                                                                                   tctcTTCATCTATAGCCCCCGCCAGCGCGGCGAACTTAATCATCGCATACTCACCGCTCACCTGATACGCGCCAATCGGCAATTCAGTACGTTCACGCAGCTCACGCACGATGTCGAGGTACGCTCCAGCAGGTTTaa          <  2:376474/138‑1 (MQ=255)                                                                                                                          aTCTATAGCCCCCGCCAGCGCGGCGAACTTAATCATCGCATACTCACCGCTCACCTGATACGCGCCAATCGGCAATTCAGTACGTTCACGc                                                  <  2:478828/91‑1 (MQ=255)                                                                                                                          aTCTATAGCCCCCGCCAGCGCGGCGAACTTAATCATCGCATACTCACCGCTCACCTGATACGCGCCAATCGGCAATTCAGTACGTTCACGc                                                  >  1:478828/1‑91 (MQ=255)                                                                                                                          aTCTATAGCCCCCGCCAGCGCGGCGAACTTAATCATCGCATACTCACCGCTCACCTGATACGCGCCAATCGGCAATTCAGTACGTTCACGCAGCTCACGCACGATGTCGAGGTACGCTCCAGCAGGTTTAACCATCAgg  <  2:513525/139‑1 (MQ=255)                                                                                                                               taGCCCCCGCCAGCGCGGCGAACTTAATCATCGCATACTCACCGCTCACCTGATACGCGCCAATCGGCAATTCAGTACGTTCACGCAGCTCACGCACGa                                     <  2:50950/99‑1 (MQ=255)                                                                                                                               taGCCCCCGCCAGCGCGGCGAACTTAATCATCGCATACTCACCGCTCACCTGATACGCGCCAATCGGCAATTCAGTACGTTCACGCAGCTCACGCACGa                                     >  1:50950/1‑99 (MQ=255)                                                                                                                                   |                                                                                                                                  AGAGGAAACAAAATTAACGCAGAATCTTCTTCTCAGCCAAATCCAGCGCAAAGTAGCTGAAAATCAGATCCGCACCCGCACGCTTAATCGAACCTAAGCTTTCGAGCACGACTTTCTCTTCATCTATAGCACCCGCCAGCGCGGCGAACTTAATCATCGCATACTCACCGCTCACCTGATACGCGCCAATCGGCAATTCAGTACGTTCACGCAGCTCACGCACGATGTCGAGGTACGCTCCAGCAGGTTTAACCATCAGG  >  NZ_CP009273/384196‑384455

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