Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NZ_CP009273	3,981,822	A→G	I202T (ATC→ACC)	hemX ←	uroporphyrinogen‑III C‑methyltransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NZ_CP009273	3,981,822	0	A	G	100.0%	171.4 / NA	45	I202T (ATC→ACC)	hemX	uroporphyrinogen‑III C‑methyltransferase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (14/31);  total (14/31)

GTCTGAATCCATCGGCGAACCATCGCTGTCATTATCGGCCAGACGCAGGTTATCTACCTGATTTGAAAGCTGATTAAGCTTAAGGATGATGCCGTCATAATCCACCTGCGATACTGCAGAAAGGCTGGCGATATCATCGGTAATTGCCCGACGAACGGTAATCAGACTCGGGT  >  NZ_CP009273/3981734‑3981906                                                                                         |                                                                                     gTCTGAATCCATCGGCGAACCATCGCTGTCATTATCGGCCAGACGCAGGTTATCTACCTGATTTGAAAGCTGATTAAGCTTAAGGatggtg                                                                                    >  2:766664/1‑91 (MQ=255) gTCTGAATCCATCGGCGAACCATCGCTGTCATTATCGGCCAGACGCAGGTTATCTACCTGATTTGAAAGCTGATTAAGCTTAAGGatggtg                                                                                    >  2:13881/1‑91 (MQ=255) gTCTGAATCCATCGGCGAACCATCGCTGTCATTATCGGCCAGACGCAGGTTATCTACCTGATTTGAAAGCTGATTAAGCTTAAGGatggtg                                                                                    >  2:1062570/1‑91 (MQ=255) gTCTGAATCCATCGGCGAACCATCGCTGTCATTATCGGCCAGACGCAGGTTATCTACCTGATTTGAAAGCTGATTAAGCTTAAGGatggtg                                                                                    >  1:682700/1‑91 (MQ=255)  tCTGAATCCATCGGCGAACCATCGCTGTCATTATCGGCCAGACGCAGGTTATCTACCTGATTTGAAAGCTGATTAAGCTTAAGGATGGTGc                                                                                   <  1:878266/91‑1 (MQ=255)  tCTGAATCCATCGGCGAACCATCGCTGTCATTATCGGCCAGACGCAGGTTATCTACCTGATTTGAAAGCTGATTAAGCTTAAGGATGGTGc                                                                                   <  1:297661/91‑1 (MQ=255)  tCTGAATCCATCGGCGAACCATCGCTGTCATTATCGGCCAGACGCAGGTTATCTACCTGATTTGAAAGCTGATTAAGCTTAAGGATGGTGc                                                                                   <  1:535637/91‑1 (MQ=255)   cTGAATCCATCGGCGAACCATCGCTGTTATTATCGGCCAGACGCAGGTTATCTACCTGATTTGAAAGCTGATTAAGCTTAAGGATGGTGcc                                                                                  <  1:347814/91‑1 (MQ=255)   cTGAATCCATCGGCGAACCATCGCTGTCATTATCGGCCAGACGCAGGTTATCTACCTGATTTGAAAGCTGATTAAGCTTAAGGATGGTGcc                                                                                  <  1:154241/91‑1 (MQ=255)     gAATCCATCGGCGAACCATCGCTGTCATTATCGGCCAGACGCAGGTTATCTACCTGATTTGAAAGCTGATTAAGCTTAAGGATGGTGCCGt                                                                                >  2:443322/1‑91 (MQ=255)       aTCCATCGGCGAACCATCGCTGTCATTATCGGCCAGACGCAGGTTATCTACCTGATTTGAAAGCTGATTAAGCTTAAGGATGGTGCCGTCa                                                                              >  2:1121351/1‑91 (MQ=255)          cATCGGCGAACCATCGCTGTCATTATCGGCCAGACGCAGGTTATCTACCTGATTTGAAAGCTGATTAAGCTTAAGGATGGTGCCGTCATaa                                                                           <  1:867112/91‑1 (MQ=255)             cGGCGAACCATCGCTGTCATTATCGGCCAGACGCAGGTTATCTACCTGATTTGAAAGCTGATTAAGCTTAAGGATGGTGCCGTCATAATcc                                                                        <  2:785477/91‑1 (MQ=255)              ggCGAACCATCGCTGTCATTATCGGCCAGACGCAGGTTATCTACCTGATTTGAAAGCTGATTAAGCTTAAGGATGGTGCCGTCATAATCCa                                                                       <  1:1062570/91‑1 (MQ=255)              ggCGAACCATCGCTGTCATTATCGGCCAGACGCAGGTTATCTACCTGATTTGAAAGCTGATTAAGCTTAAGGATGGTGCCGTCATAATCCa                                                                       <  1:13881/91‑1 (MQ=255)                 gAACCATCGCTGTCATTATCGGCCAGACGCAGGTTATCTACCTGATTTGAAAGCTGATTAAGCTTAAGGATGGTGCCGTCATAATCCACCt                                                                    <  1:931305/91‑1 (MQ=255)                   aCCATCGCTGTCATTATCGGCCAGACGCAGGTTATCTACCTGATTTGAAAGCTGATTAAGCTTAAGGATGGTGCCGTCATAATCCACCTGc                                                                  <  1:994767/91‑1 (MQ=255)                      aTCGCTGTCATTATCGGCCAGACGCAGGTTATCTACCTGATTTGAAAGCTGATTAAGCTTAAGGATGGTGCCGTCATAATCCACCTGCGAt                                                               <  1:888404/91‑1 (MQ=255)                      aTCGCTGTCATTATCGGCCAGACGCAGGTTATCTACCTGATTTGAAAGCTGATTAAGCTTAAGGATGGTGCCGTCATAATCCACCTGCGAt                                                               <  1:829460/91‑1 (MQ=255)                      aTCGCTGTCATTATCGGCCAGACGCAGGTTATCTACCTGATTTGAAAGCTGATTAAGCTTAAGGATGGTGCCGTCATAATCCACCTGCGAt                                                               <  1:754242/91‑1 (MQ=255)                        cGCTGTCATTATCGGCCAGACGCAGGTTATCTACCTGATTTGAAAGCTGATTAAGCTTAAGGATGGTGCCGTCATAATCCACCTGCGATAc                                                             <  1:536009/91‑1 (MQ=255)                         gCTGTCATTATCGGCCAGACGCAGGTTATCTACCTGATTTGAAAGCTGATTAAGCTTAAGGATGGTGCCGTCATAATCCACCTGCGATACt                                                            <  2:403066/91‑1 (MQ=255)                               aTTATCGGCCAGACGCAGGTTATCTACCTGATTTGAAAGCTGATTAAGCTTAAGGATGGTGCCGTCATAATCCACCTGCGATACTGCAGaa                                                      <  1:256512/91‑1 (MQ=255)                               aTTATCGGCCAGACGCAGGTTATCTACCTGATTTGAAAGCTGATTAAGCTTAAGGATGGTGCCGTCATAATCCACCTGCGATACTGCAGaa                                                      <  1:808896/91‑1 (MQ=255)                                ttATCGGCCAGACGCAGGTTATCTACCTGATTTGAAAGCTGATTAAGCTTAAGGATGGTGCCGTCATAATCCACCTGCGATACTGCAGaaa                                                     <  2:341960/91‑1 (MQ=255)                                    cGGCCAGACGCAGGTTATCTACCTGATTTGAAAGCTGATTAAGCTTAAGGATGGTGCCGTCATAATCCACCTGCGATACTGCAGAAAGGCt                                                 >  2:529108/1‑91 (MQ=255)                                      gCCAGACGCAGGTTATCTACCTGATTTGAAAGCTGATTAAGCTTAAGGATGGTGCCGTCATAATc                                                                         >  1:13211/1‑65 (MQ=255)                                      gCCAGACGCAGGTTATCTACCTGATTTGAAAGCTGATTAAGCTTAAGGATGGTGCCGTCATAATc                                                                         <  2:13211/65‑1 (MQ=255)                                         aGACGCAGGTTATCTACCTGATTTGAAAGCTGATTAAGCTTAAGGATGGTGCCGTCATAATCCACCTGCGATACTGCAGAAAGGCTGGCGa                                            <  1:207107/91‑1 (MQ=255)                                          gACGCAGGTTATCTACCTGATTTGAAAGCTGATTAAGCTTAAGGATGGTGcc                                                                                  >  2:119614/1‑52 (MQ=255)                                          gACGCAGGTTATCTACCTGATTTGAAAGCTGATTAAGCTTAAGGATGGTGcc                                                                                  <  1:119614/52‑1 (MQ=255)                                                ggTTATCTACCTGATTTGAAAGCTGATTAAGCTTAAGGATGGTGCCGTCATAATCCACCTGCGATACTGCAGAAAGGCTGGCGATatcatc                                     <  1:961468/91‑1 (MQ=255)                                                ggTTATCTACCTGATTTGAAAGCTGATTAAGCTTAAGGATGGTGCCGTCATAATCCACCTGCGATACTGCAGAAAGGCTGGCGATatcatc                                     >  1:339470/1‑91 (MQ=255)                                                    aTCTACCTGATTTGAAAGCTGATTAAGCTTAAGGATGGTGCCGTCATAATCCACCTGCGATACTGCAGAAAGGCTGGCGATATCATCGGTa                                 >  1:678471/1‑91 (MQ=255)                                                     tCTACCTGATTTGAAAGCTGATTAAGCTTAAGGATGGTGcc                                                                                  >  1:174070/1‑41 (MQ=255)                                                     tCTACCTGATTTGAAAGCTGATTAAGCTTAAGGATGGTGcc                                                                                  <  2:174070/41‑1 (MQ=255)                                                      cTACCTGATTTGAAAGCTGATTAAGCTTAAGGATGGTGCCGTCATAATCCACCTGCGATACTGCAGAAAGGCTGGCGATATCATCGGTAAt                               <  1:835971/91‑1 (MQ=255)                                                       tACCTGATTTGAAAGCTGATTAAGCTTAAGGATGGTGCCGTCATAATCCACCTGCGATACTGCAGAAAGGCTGGCGATATCATCGGTAAtt                              <  1:674890/91‑1 (MQ=255)                                                          cTGATTTGAAAGCTGATTAAGCTTAAGGATGGTGCCGTCATAATCCACCTGCGATACTGCAGAAAGGCTGGCGATATCATCGGTAATTGcc                           >  1:59523/1‑91 (MQ=255)                                                          cTGATTTGAAAGCTGATTAAGCTTAAGGATGGTGCCGTCATAATCCACCTGCGATACTGCAGAAAGGCTGGCGATATCATCGGTAATTGcc                           >  2:752402/1‑91 (MQ=255)                                                           tGATTTGAAAGCTGATTAAGCTTAAGGATGGTGCCGTCATAATCCACCTGCGATACTGCAGAAAGGCTGGCGATATCATCGGTAATTGccc                          >  2:405716/1‑91 (MQ=255)                                                              tttGAAAGCTGATTAAGCTTAAGGATGGTGCCGTCATAATCCACCTGCGATACTGCAGAAAGGCTGGCGATATCATCGGTAATTGCCcgac                       >  2:921451/1‑91 (MQ=255)                                                                  aaaGCTGATTAAGCTTAAGGATGGTGCCGTCATAATCCACCTGCGATACTGCAGAAAGGCTGGCGATATCATCGGTAATTGCCCGACGAAc                   <  2:591802/91‑1 (MQ=255)                                                                   aaGCTGATTAAGCTTAAGGATGGTGCCGTCATAATCCACCTGCGATACTGCAGAAAGGCTGGCGATATCATCGGTAATTGCCCGACGAACg                  <  2:678471/91‑1 (MQ=255)                                                                      cTGATTAAGCTTAAGGATGGTGCCGTCATAATCCACCTGCGATACTGCAGAAAGGCTGGCGATATCATCGGTAATTGCCCGACGAACGGTa               >  1:739966/1‑91 (MQ=255)                                                                        gATTAAGCTTAAGGATGGTGCCGTCATAATCCACCTGCGATACTGCAGAAAGGCTGGCGATATCATCGGTAATTGCCCGACGAACGGTAAt             >  2:871048/1‑91 (MQ=255)                                                                              gCTTAAGGATGGTGCCGTCATAATCCACCTGCGATACTGCAGAAAGGCTGGCGATATCATCGGTAATTGCCCGACGAACGGTAATCAGACt       <  1:443322/91‑1 (MQ=255)                                                                              gCTTAAGGATGGTGCCGTCATAATCCACCTGCGATACTGCAGAAAGGCTGGCGATATCATCGGTAATTGCCCGACGAACGGTAATCAGACt       <  1:103939/91‑1 (MQ=255)                                                                                   aGGATGGTGCCGTCATAATCCACCTGCGATACTGCAGAAAGGCTGGCGATATCATCGGTAATTGCCCGACGAACGGTAATCAGACTCGGGt  <  2:809935/91‑1 (MQ=255)                                                                                   aGGATGGGGCCGTCATAATCCACCTGCGATACTGCAGAAAGGCTGGCGATATCATCGGTAATTGCCCGACGAACGGTAATCAGACTCGGGt  <  2:360226/91‑1 (MQ=255)                                                                                         |                                                                                     GTCTGAATCCATCGGCGAACCATCGCTGTCATTATCGGCCAGACGCAGGTTATCTACCTGATTTGAAAGCTGATTAAGCTTAAGGATGATGCCGTCATAATCCACCTGCGATACTGCAGAAAGGCTGGCGATATCATCGGTAATTGCCCGACGAACGGTAATCAGACTCGGGT  >  NZ_CP009273/3981734‑3981906

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