Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NZ_CP009273	1,423,602	T→C	R99R (CGT→CGC)	stfR →	prophage tail fiber N‑terminal domain‑containing protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NZ_CP009273	1,423,602	0	T	C	100.0%	129.4 / NA	41	R99R (CGT→CGC)	stfR	prophage tail fiber N‑terminal domain‑containing protein
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (17/24);  total (17/24)

CGGGACCATTACCGTGTATGAAGATTCTCAACCCGGTACGCTGAATGATTTTCTCGGTGCCATGACGGAGGATGATGCCCGTCCGGAGGCACTGCGCCGTTTTGAACTGATGGTGGAAGAGGTGGCGCGTAACGCGTCCGCGGTGGCACAGAACACGGCAG  >  NZ_CP009273/1423521‑1423681                                                                                  |                                                                                cGGGACCATTACCGTGTATGAAGATTCTCAACCCGGTACGCTGAATGATTTTCTCGGTGCCATGACGGAGGATGATGCCCGCCCGGAGGc                                                                         <  2:131853/90‑1 (MQ=255)  gggACCATTACCGTGTATGAAGATTCTCAACCCGGTACGCTGAATGATTTTCTCGGTGCCATGACGGAGGATGATGCCCGCCCGGAGGCa                                                                        >  1:1224621/1‑90 (MQ=255)  gggACCATTACCGTGTATGAAGATTCTCAACCCGGTACGCTGAATGATTTTCTCGGTGCCATGACGGAGGATGATGCCCGCCCGGAGGCa                                                                        >  1:139504/1‑90 (MQ=255)    gACCATTACCGTGTATGAAGATTCTCAACCCGGTACGCTGAATGATTTTCTCGGTGCCATGACGGAGGATGATGCCCGCCCGGAGGCACt                                                                      >  1:456511/1‑90 (MQ=255)        cctgccGGGTATGAAGATTCGCACCCCGGTCCGCTGGATGCCTTTCTCGGTGCCATGACGGAGGATGATGCCCGCCCGGAGGCACTgcgc                                                                  <  1:1435676/86‑1 (MQ=255)        aTTACCGTGTATGAAGATTCTCAACCCGGTACGCTGAATGATTTTCTCGGTGCCATGACGGAGGATGATGCCCGCCCGGAGGCACTgcgc                                                                  <  1:845789/90‑1 (MQ=255)        aTTACCGTGTATGAAGATTCTCAACCCGGTACGCTGAATGATTTTCTCGGTGCCATGACGGAGGATGATGCCCGCCCGGAGGCACTgcgc                                                                  <  1:98910/90‑1 (MQ=255)        aTTACCGTGTATGAAGATTCTCAACCCGGTACGCTGAATGATTTTCTCGGTGCCATGACGGAGGATGATGCCCGCCCGGAGGCACTgcgc                                                                  <  1:184541/90‑1 (MQ=255)        aTTACCGTGTATGAAGATTCTCAACCCGGTACGCTGAATGATTTTCTCGGTGCCATGACGGAGGATGATGCCCGCCCGGAGGCACTgcgc                                                                  <  2:412317/90‑1 (MQ=255)        aTTACCGTGTATGAAGATTCTCAACCCGGTACGCTGAATGATTTTCTCGGTGCCATGACGGAGGATGATGCCCGCCCGGAGGCACTgcgc                                                                  <  1:1159526/90‑1 (MQ=255)              gtgtATGAAGATTCTCAACCCGGTACGCTGAATGATTTTCTCGGTGCCATGACGGAGGATGATGCCCGCCCGGAGGCACTgcgc                                                                  >  1:1429672/1‑84 (MQ=255)              gtgtATGAAGATTCTCAACCCGGTACGCTGAATGATTTTCTCGGTGCCATGACGGAGGATGATGCCCGCCCGGAGGCACTgcgc                                                                  <  2:1429672/84‑1 (MQ=255)              gtgtATGAAGATTCTCAACCCGGTACGCTGAATGATTTTCTCGGTGCCATGACGGAGGATGATGCCCGCCCGGAGGCACTGCGCCGtttt                                                            >  2:13848/1‑90 (MQ=255)              gtgtATGAAGATTCTCAACCCGGTACGCTGAATGATTTTCTCGGTGCCATGACGGAGGATGATGCCCGCCCGGAGGCACTGCGCCGtttt                                                            >  1:183405/1‑90 (MQ=255)              gtgtATGAAGATTCTCAACCCGGTACGCTGAATGATTTTCTCGGTGCCATGACGGAGGATGATGCCCGCCCGGAGGCACTGCGCCGtttt                                                            >  1:1276132/1‑90 (MQ=255)                gtATGAAGATTCTCAACCCGGTACGCTGAATGATTTTCTCGGTGCCATGACGGAGGATGATGCCCGCCCGGAGGCACTGCGCCGTTTTGa                                                          >  2:140056/1‑90 (MQ=255)                gtATGAAGATTCTCAACCCGGTACGCTGAATGATTTTCTCGGTGCCATGACGGAGGATGATGCCCGCCCGGAGGCACTGCGCCGTTTTGa                                                          >  2:459138/1‑90 (MQ=255)                gtATGAAGATTCTCAACCCGGTACGCTGAATGATTTTCTCGGTGCCATGACGGAGGATGATGCCCGCCCGGAGGCACTGCGCCGTTTTGa                                                          >  1:725351/1‑90 (MQ=255)                  aTGAAGATTCTCAACCCGGTACGCTGAATGATTTTCTCGGTGCCATGACGGAGGATGATGCCCGCCCGGAGGCACTGCGCCGTTTTGAAc                                                        <  2:534187/90‑1 (MQ=255)                  aTGAAGATTCTCAACCCGGTACGCTGAATGATTTTCTCGGTGCCATGACGGAGGATGATGCCCGCCCGGAGGCACTGCGCCGTTTTGAAc                                                        <  1:482537/90‑1 (MQ=255)                   tGAAGATTCTCAACCCGGTACGCTGAATGATTTTCTCGGTGCCATGACGGAGGATGATGCCCGCCCGGAGGCACTGCGCCGTTTTGAACt                                                       <  1:666983/90‑1 (MQ=255)                        aTTCTCAACCCGGTACGCTGAATGATTTTCTCGGTGCCATGACGGAGGATGATGCCCGCCCGGAGGCACTGCGCCGTTTTGAACTGAtgg                                                  >  1:892647/1‑90 (MQ=255)                          tctcAACCCGGTACGCTGAATGATTTTCTCGGTGCCATGACGGAGGATGATGCCCGCCCGGAGGCACTGCGCCGTTTTGAACTGAtggtg                                                >  1:1536587/1‑90 (MQ=255)                          tctcAACCCGGTACGCTGAATGATTTTCTCGGTGCCATGACGGAGGATGATGCCCGCCCGGAGGCACTGCGCCGTTTTGAACTGAt                                                    >  2:1120272/1‑86 (MQ=255)                          tctcAACCCGGTACGCTGAATGATTTTCTCGGTGCCATGACGGAGGATGATGCCCGCCCGGAGGCACTGCGCCGTTTTGAACTGAt                                                    <  1:1120272/86‑1 (MQ=255)                            tcAACCCGGTACGCTGAATGATTTTCTCGGTGCCATGACGGAGGATGATGCCCGCCCGGAGGCACTGCGCCGTTTTGAACTGATGGTGGa                                              >  2:1273170/1‑90 (MQ=255)                             cAACCCGGTACGCTGAATGATTTTCTCGGTGCCATGACGGAGGATGATGCCCGCCCGGAGGCACTGCGCCGTTTTGAACTGATGGTGGaa                                             <  1:140056/90‑1 (MQ=255)                                   ggTACGCTGAATGATTTTCTCGGTGCCATGACGGAGGATGATGCCCGCCCGGAGGCACTGCGCCGTTTTGAACTGATGGTGGAAGAGGTg                                       <  2:892647/90‑1 (MQ=255)                                   ggTACGCTGAATGATTTTCTCGGTGCCATGACGGAGGATGATGCCCGCCCGGAGGCACTGCGCCGTTTTGAACTGATGGTGGAAGAGGTg                                       <  2:725351/90‑1 (MQ=255)                                    gTACGCTGAATGATTTTCTCGGTGCCATGACGGAGGATGATGCCCGCCCGGAGGCACTgcgc                                                                  >  1:550044/1‑62 (MQ=255)                                    gTACGCTGAATGATTTTCTCGGTGCCATGACGGAGGATGATGCCCGCCCGGAGGCACTgcgc                                                                  <  2:550044/62‑1 (MQ=255)                                         cTGAATGATTTTCTCGGTGCCATGACGGAGGATGATGCCCGCCCGGAGGCACTGCGCCGTTTTGAACTGATGGTGGAAGAGGTGGCGCGt                                 <  2:693599/90‑1 (MQ=255)                                            aaTGATTTTCTCGGTGCCATGACGGAGGATGATGCCCGCCCGGAGGCACTGCGCCGTTTTGAACTGATGGTGGAAGAGGTGGCGCGTAAc                              <  1:922589/90‑1 (MQ=255)                                                    tctcGGTGCCATGACGGAGGATGATGCCCGCCCGGAGGCACTGCGCCGTTTTGAACTGATGGTGGAAGAGGTGGCGCGTAACGCGTCcgc                      <  1:1363767/90‑1 (MQ=255)                                                       cGGTGCCATGACGGAGGATGATGCCCGCCCGGAGGCACTGCGCCGTTTTGAACTGATGGTGGAAGAGGTGGCGCGTAACGCGTCCGCGGt                   <  2:1545389/90‑1 (MQ=255)                                                       cGGTGCCATGACGGAGGATGATGCCCGCCCGGAGGCACTGCGCCGTTTTGAACTGATGGTGGAAGAGGTGGCGCGTAACGCGTCCGCGGt                   <  1:1521485/90‑1 (MQ=255)                                                       cGGTGCCATGACGGAGGATGATGCCCGCCCGGAGGCACTGCGCCGTTTTGAACTGATGGTGGAAGAGGTGGCGCGTAACGCGTCCGCGGt                   <  1:478295/90‑1 (MQ=255)                                                            ccATGACGGAGGATGATGCCCGCCCGGAGGCACTGCGCCGTTTTGAACTGATGGTGGAAGAGGTGGCGCGTAACGCGTCCGCGGTGGcac              <  2:709231/90‑1 (MQ=255)                                                            acATGACGGAGGATGATGCCCGCCCGGAGGCACTGCGCCGTTTTGAACTGATGGTGGAAGAGGTGGCGCGTAACGCGTCCGCGGTGGcac              <  2:1345866/89‑1 (MQ=255)                                                                        atgatgCCCGCCCGGAGGCACTGCGCCGTTTTGAACTGATGGTGGAAGAGGTGGCGCGTAACGCGTCCGCGGTGGCACAGAACACGGCAg  >  1:1414440/1‑90 (MQ=255)                                                                        atgatgCCCGCCCGGAGGCACTGCGCCGTTTTGAACTGATGGTGGAAGAGGTGGCGCGTAACGCGTCCGCGGTGGCACAGAACACGGCAg  >  1:18303/1‑90 (MQ=255)                                                                                  |                                                                                CGGGACCATTACCGTGTATGAAGATTCTCAACCCGGTACGCTGAATGATTTTCTCGGTGCCATGACGGAGGATGATGCCCGTCCGGAGGCACTGCGCCGTTTTGAACTGATGGTGGAAGAGGTGGCGCGTAACGCGTCCGCGGTGGCACAGAACACGGCAG  >  NZ_CP009273/1423521‑1423681

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