Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NZ_CP009273	3,980,918	(T)5→4	coding (325/1197 nt)	hemY ←	protoheme IX biogenesis protein HemY

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NZ_CP009273	3,980,914	0	T	.	100.0%	130.8 / NA	31	coding (329/1197 nt)	hemY	protoheme IX biogenesis protein HemY
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base . (15/16);  total (15/16)

TGGTTGGCGCGTGCTTCATCACCACGTTGTTGCGCGGCTTCGGCAGCCAGTAGATAGTTCACCACCGGTTGTTCCGCGTGATCGGCATTTTTCGCCATCAGCTTTTCAACTTGCTGATAATCGCCTTCCGCCAGTTTCAGCAGCGCCTGTTCGGTCTGCTTACGTGCAC  >  NZ_CP009273/3980827‑3980995                                                                                        |                                                                                  tggttggCGCGTGCTTCATCACCACGTTGTTGCGCGGCTTCGGCAGCCAGTAGATAGTTCACCACCGGTTGTTCCGCGTGATCGGCAttt                                                                                 >  2:340992/1‑90 (MQ=255)     tggCGCGTGCTTCATCACCACGTTGTTGCGCGGCTTCGGCAGCCAGTAGATAGTTCACCACCGGTTGTTCCGCGTGATCGGCAttttcgc                                                                             >  1:8918/1‑87 (MQ=255)      ggCGCGTGCTTCATCACCACGTTGTTGCGCGGCTTCGGCAGCCAGTAGATAGTTCACCACCGGTTGTTCCGCGTGATCGGCAttttcgcc                                                                            >  1:317512/1‑86 (MQ=255)         gcgTGCTTCATCACCACGTTGTTGCGCGGCTTCGGCAGCCAGTAGATAGTTCACCACCGGTTGTTCCGCGTGATCGGCA‑TTTTCGCCATc                                                                        <  2:272433/90‑1 (MQ=255)          cgTGCTTCATCACCACGTTGTTGCGCGGCTTCGGCAGCCAGTAGATAGTTCACCACCGGTTGTTCCGCGTGATCGGCA‑TTTTCGCCATCa                                                                       <  1:105586/90‑1 (MQ=255)          cgTGCTTCATCACCACGTTGTTGCGCGGCTTCGGCAGCCAGTAGATAGTTCACCACCGGTTGTTCCGCGTGATCGGCA‑TTTTCGCCATCa                                                                       <  1:504779/90‑1 (MQ=255)          cgTGCTTCATCACCACGTTGTTGCGCGGCTTCGGCAGCCAGTAGATAGTTCACCACCGGTTGTTCCGCGTGATCGGCA‑TTTTCGCCATCa                                                                       <  1:482028/90‑1 (MQ=255)           gTGCTTCATCACCACGTTGTTGCGCGGCTTCGGCAGCCAGTAGATAGTTCACCACCGGTTGTTCCGCGTGATCGGCA‑TTTTCGCCATCAg                                                                      >  1:446726/1‑90 (MQ=255)                     accacGTTGTTGCGCGGCTTCGGCAGCCAGTAGATAGTTCACCACCGGTTGTTCCGCGTGATCGGCA‑TTTTCGCCATCAGCTTTTCAACt                                                            <  2:386543/90‑1 (MQ=255)                      ccacGTTGTTGCGCGGCTTCGGCAGCCAGTAGATAGTTCACCACCGGTTGTTCCGCGTGATCGGCA‑TTTTCGCCATCAGCTTTTCAACtt                                                           >  1:159121/1‑90 (MQ=255)                        acGTTGTTGCGCGGCTTCGGCAGCCAGTAGATAGTTCACCACCGGTTGTTCCGCGTGATCGGCA‑TTTTCGCCATCAGCTTTTCAACTTGc                                                         >  2:422194/1‑90 (MQ=255)                         cGTTGTTGCGCGGCTTCGGCAGCCAGTAGATAGTTCACCACCGGTTGTTCCGCGTGATCGGCA‑TTTTCGCCATCAGCTTTTCAACTTGCt                                                        <  1:271792/90‑1 (MQ=255)                                       ttCGGCAGCCAGTAGATAGTTCACCACCGGTTGTTCCGCGTGATCGGCA‑TTTTCGCCATCAGCTTTTCAACTTGCTGATAATCGCCTTcc                                          >  2:139682/1‑90 (MQ=255)                                         cGGCAGCCAGTAGATAGTTCACCACCGGTTGTTCCGCGTGATCGGCA‑TTTTCGCCATCAGCTTTTCAACTTGCTGATAATCGCCTTCCGc                                        <  2:446726/90‑1 (MQ=255)                                                   tAGATAGTTCACCACCGGTTGTTCCGCGTGATCGGCA‑TTTTCGCCATCAGCTTTTCAACTTGCTGATAATCGCCTTCCGCCAGTTTcagc                              <  2:205297/90‑1 (MQ=255)                                                    aGATAGTTCACCACCGGTTGTTCCGCGTGATCGGCA‑TTTTCGCCATCAGCTTTTCAACTTGCTGATAATCGCCTTCCGCCAGTTTcagca                             <  1:449371/90‑1 (MQ=255)                                                        aGTTCACCACCGGTTGTTCCGCGTGATCGGCA‑TTTTCGCCATCAGCTTTTCAACTTGCTGATAATCGCCTTCCGCCAGTTTCAGCAgcgc                         >  1:427707/1‑90 (MQ=255)                                                        aGTTCACCACCGGTTGTTCCGCGTGATCGGCA‑TTTTCGCCATCAGCTTTTCAACTTGCTGATAATCGCCTTCCGCCAGTTTCAGCAgcgc                         >  1:248338/1‑90 (MQ=255)                                                         gTTCACCACCGGTTGTTCCGCGTGATCGGCA‑TTTTCGCCATCAGCTTTTCAACTTGCTGATAATCGCCTTCCGCCAGTTTCAGCAGCGcc                        >  1:523212/1‑90 (MQ=255)                                                         gTTCACCACCGGTTGTTCCGCGTGATCGGCA‑TTTTCGCCATCAGCTTTTCAACTTGCTGATAATCGCCTTCCGCCAGTTTCAGCAGCGcc                        >  2:390517/1‑90 (MQ=255)                                                          ttCACCACCGGTTGTTCCGCGTGATCGGCA‑TTTTCGCCATCAGCTTTTCAACTTGCTGATAATCGCCTTCCGCCAGTTTCAGCAGCGCCt                       >  2:191704/1‑90 (MQ=255)                                                              ccaccGGTTGTTCCGCGTGATCGGCA‑TTTTCGCCATCAGCTTTTCAACTTGCTGATAATCGCCTTCCGCCAGTTTCAGCAGCGCCTGTTc                   <  1:409519/90‑1 (MQ=255)                                                                   ggTTGTTCCGCGTGATCGGCA‑TTTTCGCCATCAGCTTTTCAACTTGCTGATAATCGCCTTCCGCCAGTTTCAGCAGCGCCTGTTCGGTCt              >  1:36956/1‑90 (MQ=255)                                                                   ggTTGTTCCGCGTGATCGGCA‑TTTTCGCCATCAGCTTTTCAACTTGCTGATAATCGCCTTCCGCCAGTTTCAGCAGCGCCTGTTCGGTCt              >  2:191934/1‑90 (MQ=255)                                                                     ttgttCCGCGTGATCGGCA‑TTTTCGCCATCAGCTTTTCAACTTGCTGATAATCGCCTTCCGCCAGTTTCAGCAGCGCCTGTTCGGTCTGc            <  1:414501/90‑1 (MQ=255)                                                                       gttCCGCGTGATCGGCA‑TTTTCGCCATCAGCTTTTCAACTTGc                                                         >  2:271761/1‑43 (MQ=38)                                                                       gttCCGCGTGATCGGCA‑TTTTCGCCATCAGCTTTTCAACTTGc                                                         <  1:271761/43‑1 (MQ=38)                                                                       gttCCGCGTGATCGGCA‑TTTTCGCCATCAGCTTTTCAACTTGCTGATAATCGCCTTCCGc                                        <  2:280058/60‑1 (MQ=255)                                                                       gttCCGCGTGATCGGCA‑TTTTCGCCATCAGCTTTTCAACTTGCTGATAATCGCCTTCCGc                                        >  1:280058/1‑60 (MQ=255)                                                                        ttCCGCGTGATCGGCA‑TTTTCGCCATCAGCTTTTCAACTTGCTGATAATCGCCTTCCGCCAGTTTCAGCAGCGCCTGTTCGGTCTGCTTa         <  1:476025/90‑1 (MQ=255)                                                                            gcgTGATCGGCA‑TTTTCGCCATCAGCTTTTCAACTTGCTGATAATCGCCTTCCGCCAGTTTCAGCAGCGCCTGTTc                   <  2:152557/76‑1 (MQ=255)                                                                            gcgTGATCGGCA‑TTTTCGCCATCAGCTTTTCAACTTGCTGATAATCGCCTTCCGCCAGTTTCAGCAGCGCCTGTTc                   >  1:152557/1‑76 (MQ=255)                                                                            gcgTGATCGGCA‑TTTTCGCCATCAGCTTTTCAACTTGCTGATAATCGCCTTCCGCCAGTTTCAGCAGCGCCTGTTCGGTCTGCTTACGTg     >  2:422298/1‑90 (MQ=255)                                                                               tGATCGGCA‑TTTTCGCCATCAGCTTTTCAACTTGCTGATAATCGCCTTCCGCCAGTTTCAGCAGCGCCTGTTCGGTCTGCTTACGTGCAc  <  1:340992/90‑1 (MQ=255)                                                                                        |                                                                                  TGGTTGGCGCGTGCTTCATCACCACGTTGTTGCGCGGCTTCGGCAGCCAGTAGATAGTTCACCACCGGTTGTTCCGCGTGATCGGCATTTTTCGCCATCAGCTTTTCAACTTGCTGATAATCGCCTTCCGCCAGTTTCAGCAGCGCCTGTTCGGTCTGCTTACGTGCAC  >  NZ_CP009273/3980827‑3980995

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