Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NZ_CP009273	1,503,965	T→C	D67D (GAT→GAC)	hicB →	type II toxin‑antitoxin system antitoxin HicB

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NZ_CP009273	1,503,965	0	T	C	100.0%	148.3 / NA	45	D67D (GAT→GAC)	hicB	type II toxin‑antitoxin system antitoxin HicB
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (25/20);  total (25/20)

GCAAAAGATGCTTTACTGACCGCATTTGATTTTTATTTTGAAGATAACGAGCTTATCCCTTTACCTTCGCCATTAAATAGTCACGATCACTTTATTGAAGTACCTTTGAGCGTCGCCTCTAAGGTATTGCTGTTAAATGCTTTTTTACAGTCAGAAATCACTCAGCAAGAG  >  NZ_CP009273/1503879‑1504049                                                                                       |                                                                                     gCAAAAGATGCTTTACTGACCGCATTTGATTTTTATTTTGAAGATAACGAGCTTATCCCTTTACCTTCGCCATTAAATAGTCACGACCAc                                                                                   <  1:1388284/90‑1 (MQ=255) gCAAAAGATGCTTTACTGACCGCATTTGATTTTTATTTTGAAGATAACGAGCTTATCCCTTTACCTTCGCCATTAAATAGTCACGACCAc                                                                                   <  1:359032/90‑1 (MQ=255)    aaaGATGCTTTACTGACCGCATTTGATTTTTATTTTGAAGATAACGAGCTTATCCCTTTACCTTCGCCATTAAATAGTCACGACCACttt                                                                                <  2:1179657/90‑1 (MQ=255)      aGATGCTTTACTGACCGCATTTGATTTTTATTTTGAAGATAACGAGCTTATCCCTTTACCTTCGCCATTAAATAGTCACGACCACTTTAt                                                                              <  1:576224/90‑1 (MQ=255)           cTTTACTGACCGCATTTGATTTTTATTTTGAAGATAACGAGCTTATCCCTTTACCTTCGCCATTAAATAGTCACGACCACTTTATTGAAg                                                                         <  1:1447564/90‑1 (MQ=255)                 tGACCGCATTTGATTTTTATTTTGAAGATAACGAGCTTATCCCTTTACCTTCGCCATTAAATAGTCACGACCACTTTATTGAAGTACCtt                                                                   >  1:966400/1‑90 (MQ=255)                  gACCGCATTTGATTTTTATTTTGAAGATAACGAGCTTATCCCTTTACCTTCGCCATTAAATAGTCACGACCACTTTATTGAAGTACCttt                                                                  >  2:131275/1‑90 (MQ=255)                  gACCGCATTTGATTTTTATTTTGAAGATAACGAGCTTATCCCTTTACCTTCGCCATTAAATAGTCACGACCACTTTATTGAAGTACCttt                                                                  >  1:765811/1‑90 (MQ=255)                   aCCGCATTTGATTTTTATTTTGAAGATAACGAGCTTATCCCTTTACCTTCGCCATTAAATAGTCACGACCACTTTATTGAAGTACCTTTg                                                                 >  2:814298/1‑90 (MQ=255)                    ccGCATTTGATTTTTATTTTGAAGATAACGAGCTTATCCCTTTACCTTCGCCATTAAATAGTCACGACCACTTTATTGAAGTACCTTTGa                                                                >  2:796798/1‑90 (MQ=255)                    ccGCATTTGATTTTTATTTTGAAGATAACGAGCTTATCCCTTTACCTTCGCCATTAAATAGTCACGACCACTTTATTGAAGTACCTTTGa                                                                >  2:776286/1‑90 (MQ=255)                       cATTTGATTTTTATTTTGAAGATAACGAGCTTATCCCTTTACCTTCGCCATTAAATAGTCACGACCACTTTATTGAAGTACCTTTGAGCg                                                             >  2:757974/1‑90 (MQ=255)                        aTTTGATTTTTATTTTGAAGATAACGAGCTTATCCCTTTACCTTCGCCATTAAATAGTCACGACCACTTTATTGAAGTACCTTTGAGCGt                                                            <  1:757974/90‑1 (MQ=255)                         tttGATTTTTATTTTGAAGATAACGAGCTTATCCCTTTACCTTTGCCATTAAATAGTCACGACCACTTTATTGAAGTACCTTTGAGCGTc                                                           <  1:51306/90‑1 (MQ=255)                         tttGATTTTTATTTTGAAGATAACGAGCTTATCCCTTTACCTTCGCCATTAAATAGTCACGACCACTTTATTGAAGTACCTTTGAGCGTc                                                           <  1:814298/90‑1 (MQ=255)                         tttGATTTTTATTTTGAAGATAACGAGCTTATCCCTTTACCTTCGCCATTAAATAGTCACGACCACTTTATTGAAGTACCTTTGAGCGTc                                                           <  1:1472774/90‑1 (MQ=255)                         tttGATTTTTATTTTGAAGATAACGAGCTTATCCCTTTACCTTCGCCATTAAATAGTCACGACCACTTTATTGAAGTACCTTTGAGCGTc                                                           <  2:331591/90‑1 (MQ=255)                         tttGATTTTTATTTTGAAGATAACGAGCTTATCCCTTTACCTTCGCCATTAAATAGTCACGACCACTTTATTGAAGTACCTTTGAGCGTc                                                           <  1:1383293/90‑1 (MQ=255)                         tttGATTTTTATTTTGAAGATAACGAGCTTATCCCTTTACCTTCGCCATTAAATAGTCACGACCACTTTATTGAAGTACCTTTGAGCGTc                                                           <  1:1272421/90‑1 (MQ=255)                          ttGATTTTTATTTTGAAGATAACGAGCTTATCCCTTTACCTTCGCCATTAAATAGTCACGACCACTTTATTGAAGTACCTTTGAGCGTCg                                                          <  2:1011762/90‑1 (MQ=255)                               ttttATTTTGAAGATAACGAGCTTATCCCTTTACCTTCGCCATTAAATAGTCACGACCACTTTATTGAAGTACCTTTGAGCGTCGCctct                                                     <  2:987570/90‑1 (MQ=255)                                tttATTTTGAAGATAACGAGCTTATCCCTTTACCTTCGCCATTAAATAGTCACGACCACTTTATTGAAGTACCTTTGAGCGTCGCCTCTa                                                    >  1:1565504/1‑90 (MQ=255)                                 ttATTTTGAAGATAACGAGCTTATCCCTTTACCTTCGCCATTAAATAGTCACGACCACTTTATTGAAGTACCTTTGAGCGTCGCCTCTaa                                                   <  1:335865/90‑1 (MQ=255)                                                     ttATCCCTTTACCTTCGCCATTAAATAGTCACGACCACTTTATTGAAGTACCTTTGAGCGTCGCCTCTAAGGTATTg                                            >  1:1000109/1‑77 (MQ=255)                                                     ttATCCCTTTACCTTCGCCATTAAATAGTCACGACCACTTTATTGAAGTACCTTTGAGCGTCGCCTCTAAGGTATTg                                            <  2:1000109/77‑1 (MQ=255)                                                      tATCCCTTTACCTTCGCCATTAAATAGTCACGACCACTTTATTGAAGTACCTTTGAGCGTc                                                           >  1:1040447/1‑61 (MQ=255)                                                      tATCCCTTTACCTTCGCCATTAAATAGTCACGACCACTTTATTGAAGTACCTTTGAGCGTc                                                           <  2:1040447/61‑1 (MQ=255)                                                      tATCCCTTTACCTTCGCCATTAAATAGTCACGACCACTTTATTGAAGTACCTTTGAGCGTCGCCTCTAAGGTATTGCTGTTAAATGCttt                              >  2:1475746/1‑90 (MQ=255)                                                             ttACCTTCGCCATTAAATAGTCACGACCACTTTATTGAAGTACCTTTGAGCGTCGCCTCTAAGGTATTGCTGTTAAATGCTTTTTTACAg                       >  2:598556/1‑90 (MQ=255)                                                                  ttCGCCATTAAATAGTCACGACCACTTTATTGAAGTACCTTTGAGCGTCGCCTCTAAGGCATTGCTGTTAAATGCTTTTTTACAGTCAGa                  >  1:1359638/1‑90 (MQ=255)                                                                   tCGCCATTAAATAGTCACGACCACTTTATTGAAGTACCTTTGAGCGTCGc                                                         >  1:426951/1‑50 (MQ=255)                                                                   tCGCCATTAAATAGTCACGACCACTTTATTGAAGTACCTTTGAGCGTCGc                                                         <  2:426951/50‑1 (MQ=255)                                                                   tCGCCATTAAATAGTCACGACCACTTTATTGAAGTACCTTTGAGCGTCGCCTCTAAGGTATTGCTGTTAAATGCTTTTTTACAGTCAGaa                 >  2:1090699/1‑90 (MQ=255)                                                                   tCGCCATTAAATAGTCACGACCACTTTATTGAAGTACCTTTGAGCGTCGCCTCTAAGGTATTGCTGTTAAATGCTTTTTTACAGTCAGaa                 >  1:843247/1‑90 (MQ=255)                                                                   tCGCCATTAAATAGTCACGACCACTTTATTGAAGTACCTTTGAGCGTCGCCTCTAAGGTATTGCTGTTAAATGCTTTTTTACAGTCAGaa                 >  1:642995/1‑90 (MQ=255)                                                                   tCGCCATTAAATAGTCACGACCACTTTATTGAAGTACCTTTGAGCGTCGCCTCTAAGGTATTGCTGTTAAATGCTTTTTTACAGTCAGaa                 >  2:469080/1‑90 (MQ=255)                                                                    cGCCATTAAATAGTCACGACCACTTTATTGAAGTACCTTTGAGCGTCGCCTCTAAGGTATTGCTGTTAAATGCTTTTTTACAGTCAGaaa                >  2:1478021/1‑90 (MQ=255)                                                                    cGCCATTAAATAGTCACGACCACTTTATTGAAGTACCTTTGAGCGTCGCCTCTAAGGTATTGCTGTTAAATGCTTTTTTACAGTCAGaaa                >  1:191940/1‑90 (MQ=255)                                                                              tAGTCACGACCACTTTATTGAAGTACCTTTGAGCGTCGCCTCTAAGGTATTGCTGTTAAATGCTTTTTTACAGTCAGAAATCACTCAGCa      <  2:1045400/90‑1 (MQ=255)                                                                              tAGTCACGACCACTTTATTGAAGTACCTTTGAGCGTCGCCTCTAAGGTATTGCTGTTAAACGCTTTTTTACAGTCAGAAATCACTCAGCa      <  2:167818/90‑1 (MQ=255)                                                                                gTCACGACCACTTTATTGAAGTACCTTTGAGCGTCGCCTCTAAGGTATTGCTGTTAAATGCTTTTTTACAGTCAGAAATCACTCAGCAag    >  1:29611/1‑90 (MQ=255)                                                                                gTCACGACCACTTTATTGAAGTACCTTTGAGCGTCGCCTCTAAGGTATTGCTGTTAAATGCTTTTTTACAGTCAGAAATCACTCAGCAag    >  2:1012079/1‑90 (MQ=255)                                                                                  cACGACCACTTTATTGAAGTACCTTTGAGCGTCGCCTCTAAGGTATTGCTGTTAAATGCTTTTTTACAGTCAGAAATCACTCAGCAagag  >  2:254056/1‑90 (MQ=255)                                                                                  cACGACCACTTTATTGAAGTACCTTTGAGCGTCGCCTCTAAGGTATTGCTGTTAAATGCTTTTTTACAGTCAGAAATCACTCAGCAagag  >  1:364546/1‑90 (MQ=255)                                                                                  cACGACCACTTTATTGAAGTACCTTTGAGCGTCGCCTCTAAGGTATTGCTGTTAAATGCTTTTTTACAGTCAGAAATCACTCAGCAagag  >  1:1305111/1‑90 (MQ=255)                                                                                       |                                                                                     GCAAAAGATGCTTTACTGACCGCATTTGATTTTTATTTTGAAGATAACGAGCTTATCCCTTTACCTTCGCCATTAAATAGTCACGATCACTTTATTGAAGTACCTTTGAGCGTCGCCTCTAAGGTATTGCTGTTAAATGCTTTTTTACAGTCAGAAATCACTCAGCAAGAG  >  NZ_CP009273/1503879‑1504049

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