Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NZ_CP009273	3,036,314	A→G	F65S (TTC→TCC)	yqfA ←	hemolysin III family protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NZ_CP009273	3,036,314	0	A	G	100.0%	175.2 / NA	50	F65S (TTC→TCC)	yqfA	hemolysin III family protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (24/26);  total (24/26)

GAGAATCCAGCCCCACCAGCAAAAACGGCGTGTAGGTTCCGGCAATCAACAGGTAAATAGCGCAATGGTCAAATTTCTTCAGCCACATTTTTGCCCGTTGATGAGGAATGGCGTGATAGAGCGTCGAAGCGAGGAACAGCAGGATCATACTGCCGCCATAGAGGCTGTAGCTGGTTATCGCTGTGGCGCTGGCATTAAGATCCACCGCCTGAACCAGTAGCAACACCAGCCCAACGATACCAAACACCAACCCAATGCCGTGACTGACG  >  NZ_CP009273/3036180‑3036448                                                                                                                                       |                                                                                                                                       gagaATCCAGCCCCACCAGCAAAAACGGCGTGTAGGTTCCGGCAATCAACAGGTAAATAGCGCAATGGTCAAATTTCTTCAGCCACATTTTTGCCCGTTGATGAGGAATGGCGTGATAGAGCGTCGAAGCGAGGGAcag                                                                                                                                    <  2:712684/139‑1 (MQ=255)   gaATCCAGCCCCACCAGCAAAAACGGCGTGTAGGTTCCGGCAATCAACAGGTAAATAGCGCAATGGTCAAATTTCTTCAGCCACATTTTTGCCCGTTGATGAGGAATGGCGTGATAGAGCGTCGAAGCGAGGGAcagca                                                                                                                                  <  1:407526/139‑1 (MQ=255)    aaTCCAGCCCCACCAGCAAAAACGGCGTGTAGGTTCCGGCAATCAACAGGTAAATAGCGCAATGGTCAAATTTCTTCAGCCACATTTTTGCCCGTTGATGAGGAATGGCGTGATAGAGCGTCGAAGCGAGGGAcagcag                                                                                                                                 <  1:436433/139‑1 (MQ=255)      tCCAGCCCCACCAGCAAAAACGGCGTGTAGGTTCCGGCAATCAACAGGTAAATAGCGCAATGGTCAAATTTCTTCAGCCACATTTTTGCCCGTTGATGAGGAATGGCGTGATAGAGCGTCGAAGCGAGGGACAGCAGGa                                                                                                                               >  2:118770/1‑139 (MQ=255)         aGCCCCACCAGCAAAAACGGCGTGTAGGTTCCGGCAATCAACAGGTAAATAGCGCAATGGTCAAATTTCTTCAGCCACATTTTTGCCCGTTGATGGGGAATGGCGTGATAGAGCGTCGAAGCGAGGGACAGCAGGATCa                                                                                                                            <  1:454737/139‑1 (MQ=255)          gCCCCACCAGCAAAAACGGCGTGTAGGTTCCGGCAATCAAAAGGTAAATAGCGCAATGGACAAATTTCTTCAGCCACATTTTTGACCGTTGATGCGGAATGGCGTGATAGAGCGTCGAAGCGAGGGACAGAAGGATCa                                                                                                                            >  2:61566/1‑138 (MQ=255)                        aaCGGCGTGTAGGTTCCGGCAATCAACAGGTAAATAGCGCAATGGTCAAATTTCTTCAGCCACATTTTTGCCCGTTGATGAGGAATGGCGTGATAGAGCGTCGAAGCGAGGGACAGCAGGATCATACTGCCGCCATa                                                                                                               >  1:74997/1‑137 (MQ=255)                                 tAGGTTCCGGCAATCAACAGGTAAATAGCGCAATGGTCAAATTTCTTCAGCCACATTTTTGCCCGTTGATGAGGAATGGCGTGATAGAGCGTCGAAGCGAGGGACAGCAGGATCATACTGCCGCCATAGAGGCTGTAg                                                                                                     <  1:399408/138‑1 (MQ=255)                                     ttCCGGCAATCAACAGGTAAATAGCGCAATGGTCAAATTTCTTCAGCCACATTTTTGCCCGTTGATGAGGAATGGCGTGATAGAGCGTCGAAGCGAGGGACAGCAGGATCATACTGCCGCCATAGAGGCTGTAGCTGGt                                                                                                <  1:235064/139‑1 (MQ=255)                                      tCCGGCAATCAACAGGTAAATAGCGCAATGGTCAAATTTCTTCAGCCACATTTTTGCCCGTTGATGAGGAATGGCGTGATAGAGCGTCGAAGCGAGGGAcag                                                                                                                                    <  2:608312/102‑1 (MQ=255)                                      tCCGGCAATCAACAGGTAAATAGCGCAATGGTCAAATTTCTTCAGCCACATTTTTGCCCGTTGATGAGGAATGGCGTGATAGAGCGTCGAAGCGAGGGAcag                                                                                                                                    >  1:608312/1‑102 (MQ=255)                                       ccGGCAATCAACAGGTAAATAGCGCAATGGTCAAATTTCTTCAGCCACATTTTTGCCCGTTGATGAGGAATGGCGTGATAGAGCGTCGAAGCGAGGGACAGCAGGATCGTACTGCCGCCATAGAGGCTGTAGCTGGTTa                                                                                              >  2:830050/1‑139 (MQ=255)                                          gCAATCAACAGGTAAATAGCGCAATGGTCAAATTTCTTCAGCCACATTTTTGCCCGTTGATGAGGAATGGCGTGATAGAGCGTCGAAGCGAGGGACAGCAGGATCATACTGCCGCCATAGAGGCTGTAGCTGGTTATCg                                                                                           >  1:390321/1‑139 (MQ=255)                                                  cAGGTAAATAGCGCAATGGTCAAATTTCTTCAGCCACATTTTTGCCCGTTGATGAGGAATGGCGTGATAGAGCGTCGAAGCGAGGGACAGCAGGATCATACTGCCGCCATAGAGGCTGTAGCTGGTTATCGCTGTGgcg                                                                                   <  1:785617/139‑1 (MQ=255)                                                  cAGGTAAATAGCGCAATGGTCAAATTTCTTCAGCCACATTTTTGCCCGTTGATGAGGAATGGCGTGATAGAGCGTCGAAGCGAGGGACAGCAGGATCATACTGCCGCCATAGAGGCTGTAGCTGGTTATCGCTGTGgcg                                                                                   <  2:74997/139‑1 (MQ=255)                                                    ggaaaaTAGCGCAATGGTCAAATTTCTTCAGCCACATTTTTGCCCGTTGATGAGGAATGGCGTGATAGAGCGTCGAAGCGAGGGACAGCAGGATCATACTGCCGCCATAGAGGCTGTAGCTGGTTATCGCTGTGGCGCt                                                                                 <  2:390321/136‑1 (MQ=255)                                                           aGCGCAATGGTCAAATTTCTTCAGCCACATTTTTGCCCGTTGATGAGGAATGGCGTGATAGAGCGTCGAAGCGAGGGACAGCAGGATCATACTGCCGCCATAGAGGCTGTAGCTGGTTATCGCTGTGGCGCTGGCAtt                                                                           >  2:117600/1‑138 (MQ=255)                                                               cAATGGTCAAATTTCTTCAGCCACATTTTTGCCCGTTGATGAGGAATGGCGTGATAGAGCGTCGAAGCGAGGGACAGCAGGATCATACTGCCGCCATAGAGGCTGTAGCTGGTTATCGCTGTGGCGCTGGCAtt                                                                           >  2:155786/1‑134 (MQ=255)                                                                  tGGTCAAATTTCTTCAGCCACATTTTTGCCCGTTGATGAGGAATGGCGTGATAGAGCGTCGAAGCGAGGGACAGCAGGATCATACTGCCGCCATAGAGGCTGTAGCTGGTTATCGCTGTGGCGCTGGCATTAAGATCCa                                                                   >  1:631373/1‑139 (MQ=255)                                                                    gTCAAATTTCTTCAGCCACATTTTTGCCCGTTGATGAGGAATGGCGTGATAGAGCGTCGAAGCGAGGGACAGCAGGATCATACTGCCGCCATAGAGGCTGTAGCTGGTTATCGCTGTGGCGCt                                                                                 <  1:455883/123‑1 (MQ=255)                                                                    gTCAAATTTCTTCAGCCACATTTTTGCCCGTTGATGAGGAATGGCGTGATAGAGCGTCGAAGCGAGGGACAGCAGGATCATACTGCCGCCATAGAGGCTGTAGCTGGTTATCGCTGTGGCGCt                                                                                 >  2:455883/1‑123 (MQ=255)                                                                    gTCAAATTTCTTCAGCCACATTTTTGCCCGTTGATGAGGAATGGCGTGATAGAGCGTCGAAGCGAGGGACAGCAGGATCATACTGCCGCCATAGAGGCTGTAGCTGGTTATCGCTGTGGCGCTGGCATTAAGATCCAcc                                                                 >  2:573119/1‑139 (MQ=255)                                                                          tttCTTCAGCCACATTTTTGCCCGTTGATGAGGAATGGCGTGATAGAGCGTCGAAGCGAGGGACAGCAGGATCATACTGCCGCCATAGAGGCTGTAGCTGGTTATCGCTGTGGCGCTGGCATTAAGATCCACCGCCTg                                                            >  1:37217/1‑138 (MQ=255)                                                                              ttcAGCCACATTTTTGCCCGTTGATGAGGAATGGCGTGATAGAGCGTCGAAGCGAGGGACAGCAGGATCATACTGCCGCCATAGAGGCTGTAGCTGGTTATCGCTGTGgcg                                                                                   >  2:463726/1‑111 (MQ=255)                                                                              ttcAGCCACATTTTTGCCCGTTGATGAGGAATGGCGTGATAGAGCGTCGAAGCGAGGGACAGCAGGATCATACTGCCGCCATAGAGGCTGTAGCTGGTTATCGCTGTGgcg                                                                                   <  1:463726/111‑1 (MQ=255)                                                                                cAGCCACATTTTTGTCCGTTGATGAGGAATGGCGTGATAGAGCGTCGAAGCGAGGGACAGCAGGATCATACTGCCGCCATAGAGGCTGTAGCTGGTTATCGCTGTGGCGCTGGCATTAAGATCCACCGCCTGAAcc                                                        <  1:721257/136‑1 (MQ=255)                                                                                   ccACATTTTTGCCCGTTGATGAGGAATGGCGTGATAGAGCGTCGAAGCGAGGGACAGCAGGATCATACTGCCGCCATAGAGGCTGTAGCTGGTTATCGCTGTGGCGCTGGCATTAAGATCCACCGCCTGAACCAGTAGc                                                  <  1:624205/139‑1 (MQ=255)                                                                                        tttttGCCCGTTGATGAGGAATGGCGTGATAGAGCGTCGAAGCGAGGGACAGCAGGATCATACTGCCGCCATAGAGGCTGTAGCTGGTTATCGCTGTGGCGCTGGCATTAAGATCCACCGCCTGAACCAGTAGCAacac                                             <  1:118770/139‑1 (MQ=255)                                                                                            tGCCCGTTGATGAGGAATGGCGTGATAGAGCGTCGAAGCGAGGGACAGCAGGATCATACTGCCGCCATAGAGGCTGTAGCTGGTTATCGCTGTGGCGCTGGCATTAAGATCCACCGCCTGAACCAGTAGCAACACCAg                                          <  1:44545/138‑1 (MQ=255)                                                                                                        aGGAATGGCGTGATAGAGCGTCGAAGCGAGGGACAGCAGGATCATACTGCCGCCATAGAGGCTGTAGCTGGTTATCGCTGTGGCGCTGGCATTAAGATCCACCGCCTGAACCAGTAGCAACACCAGCCCAACGATACCa                             <  1:155786/139‑1 (MQ=255)                                                                                                             tGGCGTGATAGAGCGTCGAAGCGAGGGACAGCAGGATCATACTGCCGCCATAGAGGCTGTAGCTGGTTATCGCTGTGGCGCTGGCATTAAGATCCACCGCCTGAACCAGTAGCAACACCAGCCCAACGATACCa                             <  2:77649/134‑1 (MQ=255)                                                                                                             tGGCGTGATAGAGCGTCGAAGCGAGGGACAGCAGGATCATACTGCCGCCATAGAGGCTGTAGCTGGTTATCGCTGTGGCGCTGGCATTAAGATCCACCGCCTGAACCAGTAGCAACACCAGCCCAACGATACCa                             >  1:77649/1‑134 (MQ=255)                                                                                                             tGGCGTGATAGAGCGTCGAAGCGAGGGACAGCAGGATCATACTGCCGCCATAGAGGCTGTAGCTGGTTATCGCTGTGGCGCTGGCATTAAGATCCACCGCCTGAACCAGTAGCAACACCAGCCCAACGATACCAAacac                        <  1:583312/139‑1 (MQ=255)                                                                                                               gCGTGATAGAGCGTCGAAGCGAGGGACAGCAGGATCATACTGCCGCCATAGAGGCTGTAGCTGGTTATCGCTGTGGCGCTGGCATTAAGATCCACCGCCTGAACCAGTAGCAACACCAGCCCAACGATACCAAACAcc                       >  1:800811/1‑138 (MQ=255)                                                                                                                cGTGATAGAGCGTCGAAGCGAGGGACAGCAGGATCATACTGCCGCCATAGAGGCTGTAGCTGGTTATCGCTGTGGCGCTGGCATTAAGATCCACCGCCTGAACCAGTAGCAACACCAGCCCAACGATACCAAacgc                        <  1:42268/136‑3 (MQ=255)                                                                                                                cGTGATAGAGCGTCGAAGCGAGGGACAGCAGGATCATACTGCCGCCATAGAGGCTGTAGCTGGTTATCGCTGTGGCGCTGGCATTAAGATCCACCGCCTGAACCAGTAGCAACACCAGCCCAACGATACCAAacac                        >  2:42268/1‑136 (MQ=255)                                                                                                                  tGATAGAGCGTCGAAGCGAGGGACAGCAGGATCATACTgc                                                                                                                      >  1:607080/1‑40 (MQ=255)                                                                                                                  tGATAGAGCGTCGAAGCGAGGGACAGCAGGATCATACTgc                                                                                                                      <  2:607080/40‑1 (MQ=255)                                                                                                                    aTAGAGCGTCGAAGCGAGGGACAGCAGGATCATACTGCCGCCATAGAGGCTGTAGCTGGTTATCGCTGTGgcg                                                                                   >  2:89538/1‑73 (MQ=255)                                                                                                                    aTAGAGCGTCGAAGCGAGGGACAGCAGGATCATACTGCCGCCATAGAGGCTGTAGCTGGTTATCGCTGTGgcg                                                                                   <  1:89538/73‑1 (MQ=255)                                                                                                                    aTAGAGCGTCGAAGCGAGGGACAGCAGGATCATACTGCCGCCATAGAGGCTGTAGCTGGTTATCGCTGTGGCGCTGGCATTAAGATCCACCGCCTGAACCAGTAGCAACACCAGCCCAACGATACCAAACACCAAccc                  >  2:760709/1‑138 (MQ=255)                                                                                                                       gagCGTCGAAGCGAGGGACAGCAGGATCATACTgcc                                                                                                                     >  2:358677/1‑36 (MQ=255)                                                                                                                       gagCGTCGAAGCGAGGGACAGCAGGATCATACTgcc                                                                                                                     <  1:358677/36‑1 (MQ=255)                                                                                                                          cGTCGAAGCGAGGGACAGCAGGATCATACTGCCGCCATAGAGGCTGTAGCTGGt                                                                                                <  2:107621/54‑1 (MQ=255)                                                                                                                          cGTCGAAGCGAGGGACAGCAGGATCATACTGCCGCCATAGAGGCTGTAGCTGGt                                                                                                >  1:107621/1‑54 (MQ=255)                                                                                                                          cGTCGAAGCGAGGGACAGCAGGATCATACTGCCGCCATAGAGGCTGTAGCTGGTTATCGCTGTGGCGCTGGCATTAAGATCCACCGCCTGAACCAGTAGCAACACCAGCCCAACGATACCAAACACCAACCCAATGCCg           >  2:314699/1‑139 (MQ=255)                                                                                                                          cGTCGAAGCGAGGGACAGCAGGATCATACTGCCGCCATAGAGGCTGTAGCTGGTTATCGCTGTGGCGCTGGCATTAAGATCCACCGCCTGAACCAGTAGCAACACCAGCCCAACGATACCAAACACCAACCCAATGCCg           >  2:31399/1‑139 (MQ=255)                                                                                                                            tCGAAGCGAGGGACAGCAGGATCATACTGCCGCCATAGAGGCTGTAGCTGGTTATCGCTGTGGCGCTGGCATTAAGACCCACCGCCTGAAcc                                                        >  1:851597/1‑92 (MQ=255)                                                                                                                            tCGAAGCGAGGGACAGCAGGATCATACTGCCGCCATAGAGGCTGTAGCTGGTTATCGCTGTGGCGCTGGCATTAAGACCCACCGCCTGAAcc                                                        <  2:851597/92‑1 (MQ=255)                                                                                                                                   gAGGGACAGCAGGATCGTACTGCCGCCATAGAGGCTGTAGCTGGTTATCGCTGTGGCGCTGGCATTAAGATCCACCGCCTGAACCAGTAGCAACACCAGCCCAACGATACCAAACACCAACCCAATGCCGTGACTGACg  <  1:830050/139‑1 (MQ=255)                                                                                                                                       |                                                                                                                                       GAGAATCCAGCCCCACCAGCAAAAACGGCGTGTAGGTTCCGGCAATCAACAGGTAAATAGCGCAATGGTCAAATTTCTTCAGCCACATTTTTGCCCGTTGATGAGGAATGGCGTGATAGAGCGTCGAAGCGAGGAACAGCAGGATCATACTGCCGCCATAGAGGCTGTAGCTGGTTATCGCTGTGGCGCTGGCATTAAGATCCACCGCCTGAACCAGTAGCAACACCAGCCCAACGATACCAAACACCAACCCAATGCCGTGACTGACG  >  NZ_CP009273/3036180‑3036448

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