Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NZ_CP009273	1,079,189	T→C	S120P (TCC→CCC)	efeB →	iron uptake transporter deferrochelatase/peroxidase subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NZ_CP009273	1,079,189	0	T	C	100.0%	123.1 / NA	33	S120P (TCC→CCC)	efeB	iron uptake transporter deferrochelatase/peroxidase subunit
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (18/15);  total (18/15)

GCTTGTTGACTCAGCGTTTTGCTTTTCTGACTCAGGGCGGAGCAGCACCAGAAACGCCAAATCCGCGCCTGCCACCACTCGATTCCGGCATTCTTGGCGGCTACATTGCGCCCGATAATCTCACCATCACGTTATCGGTGGGTCACTCATTGTTTGATGAGCGCTTTGGCC  >  NZ_CP009273/1079106‑1079276                                                                                    |                                                                                        gCTTGTTGACTCAGCGTTTTGCTTTTCTGACTCAGGGCGGAGCAGCACCAGAAACGCCAAATCCGCGCCTGCCACCACTCGATCCCGGCAt                                                                                  >  1:285767/1‑91 (MQ=255) gCTTGTTGACTCAGCGTTTTGCTTTTCTGACTCAGGGCGGAGCAGCACCAGAAACGCCAAATCCGCGCCTGCCACCACTCGATCCCGGCAt                                                                                  >  1:673667/1‑91 (MQ=255)   ttgttgACTCAGCGTTTTGCTTTTCTGACTCAGGGCGGAGCAGCACCAGAAACGCCAAATCCGCGCCTGCCACCACTCGATCCCGGCATTc                                                                                >  1:393694/1‑91 (MQ=255)        gACTCAGCGTTTTGCTTTTCTGACTCAGGGCGGAGCAGCACCAGAAACGCCAAATCCGCGCCTGCCACCACTCGATCCCGGCATTCTTggc                                                                           <  1:585264/91‑1 (MQ=255)        gACTCAGCGTTTTGCTTTTCTGACTCAGGGCGGAGCAGCACCAGAAACGCCAAATCCGCGCCTGCCACCACTCGATCCCGGCATTCTTggc                                                                           <  2:621039/91‑1 (MQ=255)                gTTTTGCTTTTCTGACTCAGGGCGGAGCAGCACCAGAAACGCCAAATCCGCGCCTGCCACCACTCGATCCCGGCATTCTTGGCGGCTACAt                                                                   <  1:83462/91‑1 (MQ=255)                   ttGCTTTTCTGACTCAGGGCGGAGCAGCACCAGAAACGCCAAATCCGCGCCTGCCACCACTCGATCCCGGCATTCTTGGCGGCTACATTgc                                                                >  1:30930/1‑91 (MQ=255)                     tcTTTTCTGACTCAGGGCGGAGCAGCACCAGAAACGCCAAATCCGCGCCTGCCACCACTCGATCCCGGCATTCTTGGCGGCTACATTgcgc                                                              <  1:803985/90‑1 (MQ=255)                     gCTTTTCTGACTCAGGGCGGAGCAGCACCAGAAACGCCAAATCCGCGCCTGCCACCACTCGATCCCGGCATTCTTGGCGGCTACATTgcgc                                                              <  2:340206/91‑1 (MQ=255)                             gACTCAGGGCGGAGCAGCACCAGAAACGCCAAATCCGCGCCTGCCACCACTCGATCCCGGCATTCTTGGCGGCTACATTGCGCCCGATAAt                                                      >  1:87226/1‑91 (MQ=255)                              aCTCAGGGCGGAGCAGCACCAGAAACGCCAAATCCGCGCCTGCCACCACTCGATCCCGGCATTCTTGGCGGCTACATTGCGCCCGATAAtc                                                     <  1:735904/91‑1 (MQ=255)                                tCAGGGCGGAGCAGCACCAGAAACGCCAAATCCGCGCCTGCCACCACTCGATCCCGGCATTCTTGGCGGCTACATTGCGCCCGATAAtctc                                                   <  2:46928/91‑1 (MQ=255)                                tCAGGGCGGAGCAGCACCAGAAACGCCAAATCCGCGCCTGCCACCACTCGATCCCGGCATTCTTGGCGGCTACATTGCGCCCGATAAtctc                                                   <  1:697498/91‑1 (MQ=255)                                 cAGGGCGGAGCAGCACCAGAAACGCCAAATCCGCGCCTGCCACCACTCGATCCCGGCATTCTTGGCGGCTACATTGCGCCCGATAATCTCa                                                  >  2:188770/1‑91 (MQ=255)                                  aGGGCGGAGCAGCACCAGAAACGCCAAATCCGCGCCTGCCACCACTCGATCCCGGCATTCTTGGCGGCTACATTGCGCCCGATAATCTCAc                                                 >  1:202018/1‑91 (MQ=255)                                   gggCGGAGCAGCACCAGAAACGCCAAATCCGCGCCTGCCACCAGTCGATCCCGGCATTCTTGGCGGCTACATTGCGCCCGATAATCTCAcc                                                >  1:817327/1‑91 (MQ=255)                                   gggCGGAGCAGCACCAGAAACGCCAAATCCGCGCCTGCCACCACTCGATCCCGGCATTCTTGGCGGCTACATTGCGCCCGATAATCTCAcc                                                >  2:742401/1‑91 (MQ=255)                                        gAGCAGCACCAGAAACGCCAAATCCGCGCCTGCCACCACTCGATCCCGGCATTCTTGGCGGCTACATTGCGCCCGATAAtc                                                     <  1:427413/81‑1 (MQ=255)                                        gAGCAGCACCAGAAACGCCAAATCCGCGCCTGCCACCACTCGATCCCGGCATTCTTGGCGGCTACATTGCGCCCGATAAtc                                                     >  2:427413/1‑81 (MQ=255)                                             gcaCCAGAAACGCCAAATCCGCGCCTGCCACCACTCGATCCCGGCATTCTTGGCGGCTACATTGCGCCCGATAATCTCACCATCACGTTAt                                      >  2:541856/1‑91 (MQ=255)                                             gcaCCAGAAACGCCAAATCCGCGCCTGCCACCACTCGATCCCGGCATTCTTGGCGGCTACATTGCGCCCGATAATCTCACCATCACGTTAt                                      >  1:521840/1‑91 (MQ=255)                                             gcaCCAGAAACGCCAAATCCGCGCCTGCCACCACTCGATCCCGGCATTCTTGGCGGCTACATTGCGCCCGATAATCTCACCATCACGTTAt                                      >  2:352916/1‑91 (MQ=255)                                                ccAGAAACGCCAAATCCGCGCCTGCCACCACTCGATCCCGGCATTCTTGGCGGCTACATTGCGCCCGATAATCTCACCATCACGTTATCgg                                   <  2:173866/91‑1 (MQ=255)                                                 cAGAAACGCCAAATCCGCGCCTGCCACCACTCGATCCCGGCATTCTTGGCGGCTACATTGCGCCCGATAATCTCACCATCACGTTATCGGt                                  <  1:357361/91‑1 (MQ=255)                                                  aGAAACGCCAAATCCGCGCCTGCCACCACTCGATCCCGGCATTCTTGGCGGCTACATTGCGCCCGATAATCTCACCATCACGTTATCGGTg                                 >  1:517722/1‑91 (MQ=255)                                                     aaCGCCAAATCCGCGCCTGCCACCACTCGATCCCGGCATTCTTGGCGGCTACATTGCGCCCGATAATCTCACCATCACGTTATCGGTGGGt                              <  1:352916/91‑1 (MQ=255)                                                        gCCAAATCCGCGCCTGCCACCACTCGATCCCGGCATTCTTGGCGGCTACATTGCGCCCGATAATCTCACCATCACGTTATCGGTGGGTCAc                           >  1:466834/1‑91 (MQ=255)                                                             aTCCGCGCCTGCCACCACTCGATCCCGGCATTCTTGGCGGCTACATTGCGCCCGATAATCTCACCATCACGTTATCGGTGGGTCACTCAtt                      <  2:87226/91‑1 (MQ=255)                                                                 gcgcCTGCCACCACTCGATCCCGGCATTCTTGGCGGCTACATTGCGCCCGATAATCTCACCATCACGTTATCGGTGGGTCACTCATTGttt                  >  1:242791/1‑91 (MQ=255)                                                                       gCCACCACTCGATCCCGGCATTCTTGGCGGCTACATTGCGCCCGATAATCTCACCATCACGTTATCGGTGGGTCACTCATTGTTTGATGAg            >  1:593084/1‑91 (MQ=255)                                                                        ccaccaCTCGATCCCGGCATTCTTGGCGGCTACATTGCGCCCGATAATCTCACCATCACGTTATCGGTGGGTCACTCATTGTTTGATGAgc           <  2:48769/91‑1 (MQ=255)                                                                        ccaccaCTCGATCCCGGCATTCTTGGCGGCTACATTGCGCCCGATAATCTCACCATCACGTTATCGGTGGGTCACTCATTGTTTGATGAgc           <  1:660928/91‑1 (MQ=255)                                                                                 gATCCCGGCATTCTTGGCGGCTACATTGCGCCCGATAATCTCACCATCACGTTATCGGTGGGTCACTCATTGTTTGATGAGCGCTTTGGcc  >  1:421457/1‑91 (MQ=255)                                                                                    |                                                                                        GCTTGTTGACTCAGCGTTTTGCTTTTCTGACTCAGGGCGGAGCAGCACCAGAAACGCCAAATCCGCGCCTGCCACCACTCGATTCCGGCATTCTTGGCGGCTACATTGCGCCCGATAATCTCACCATCACGTTATCGGTGGGTCACTCATTGTTTGATGAGCGCTTTGGCC  >  NZ_CP009273/1079106‑1079276

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