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_CP0092733,036,3140AG100.0% 115.3 / NA 31F65S (TTC→TCC) yqfAhemolysin III family protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (12/19);  total (12/19)

CAGGTAAATAGCGCAATGGTCAAATTTCTTCAGCCACATTTTTGCCCGTTGATGAGGAATGGCGTGATAGAGCGTCGAAGCGAGGAACAGCAGGATCATACTGCCGCCATAGAGGCTGTAGCTGGTTATCGCTGTGGCGCTGGCATTAAGATCCACCGCCTGA  >  NZ_CP009273/3036229‑3036391
                                                                                     |                                                                             
cAGGTAAATAGCGCAATGGTCAAATTTCTTCAGCCACATTTTTGCCCGTTGATGAGGAATGGCGTGATAGAGCGTCGAAGCGAGGGAcagc                                                                          <  2:386621/91‑1 (MQ=255)
cAGGTAAATAGCGCAATGGTCAAATTTCTTCAGCCACATTTTTGCCCGTTGATGAGGAATGGCGTGATAGAGCGTCGAAGCGAGGGAcagc                                                                          <  2:845180/91‑1 (MQ=255)
  ggTAAATAGCGCAATGGTCATATTTCTTCAGCCACATTTTTGCCCGTTGATGAGGAATGGCGTGATAGAGCGTCGAAGCGAGGGAcagcag                                                                        <  2:343252/91‑1 (MQ=255)
  ggTAAATAGCGCAATGGTCAAATTTCTTCAGCCACATTTTTGCCCGTTGATGAGGAATGGCGTGATAGAGCGTCGAAGCGAGGGAcagcag                                                                        <  1:833249/91‑1 (MQ=255)
          gcgcAATGGTCAAATTTCTTCAGCCACATTTTTGCCCGTTGATGAGGAATGGCGTGATAGAGCGTCGAAGCGAGGGACAGCAGGATCATAc                                                                >  1:579442/1‑91 (MQ=255)
            gcAATGGTCAAATTTCTTCAGCCACATTTTTGCCCGTTGATGAGGAATGGCGTGATAGAGCGTCGAAGCGAGGGACAGCAGGATCATACTg                                                              >  1:68983/1‑91 (MQ=255)
            gcAATGGTCAAATTTCTTCAGCCACATTTTTGCCCGTTGATGAGGAATGGCGTGATAGAGCGTCGAAGCGAGGGACAGCAGGATCATACTg                                                              >  2:302538/1‑91 (MQ=255)
            gcAATGGTCAAATTTCTTCAGCCACATTTTTGCCCGTTGATGAGGAATGGCGTGATAGAGCGTCGAAGCGAGGGACAGCAGGATCATACTg                                                              >  1:492711/1‑91 (MQ=255)
              aaTGGTCAAATTTCTTCAGCCACATTTTTGCCCGTTGATGAGGAATGGCGTGATAGAGCGTCGAAGCGAGGGACAGCAGGATCATACTgcc                                                            <  2:579442/91‑1 (MQ=255)
               aTGGTCAAATTTCTTCAGCCACATTTTTGCCCGTTGATGAGGAATGGCGTGATAGAGCGTCGAAGCGAGGGACAGCAGGATCATACTgccg                                                           >  1:50619/1‑91 (MQ=255)
                                    caTTTTTGCCCGTTGATGAGGAATGGCGTGATAGAGCGTCGAACCGAGGGACAGCAGGATCATACTGCCGCCATAGAGGCTGTAGCTGGtt                                      <  1:675507/91‑1 (MQ=255)
                                                ttGATGAGGAATGGCGTGATAGAGCGTCGAAGCGAGGGACAGCAGGATCATACTGCCGCCATAGAGGCTGTAGCTGGTTATCGCTGTGgcg                          >  1:412949/1‑91 (MQ=255)
                                                  gatgaGGAATGGCGTGATAGAGCGTCGAAGCGAGGGACAGCAGGATCATACTGCCGCCATAGAGGCTGTAGCTGGTTATCGCTGTGGCGCt                        <  2:212071/91‑1 (MQ=255)
                                                   gggaGGAATGGAAGGATAGAAGGTCGAGGCGAGGGACAGAAGGATCATCCTGCCGCCATAGAGGCTGTAGCTGGGTATTGCTGTGGCGCTg                       <  2:837411/89‑1 (MQ=255)
                                                   atgaGGAATGGCGTGATAGAGCGTCGAAGCGAGGGACAGCAGGATCATACTGCCGCCATAGAGGCTGTAGCTGGTTATCGCTGTGGCGCTg                       <  1:445625/91‑1 (MQ=255)
                                                       ggAATGGCGTGATAGAGCGTCGAAGCGAGGGACAGCAGGATCATACTGCCGCCATAGAGGCTGTAGCTGGTTATCGCTGTGGCGCTGGCAt                   <  1:393040/91‑1 (MQ=255)
                                                          aTGGCGTGATAGAGCGTCGAAGCGAGGGACAGCAGGATCATACTGCCGCCATAGAGGCTGTAGCTGGTTATCGCTGTGGCGCTGGCATTaa                >  2:740823/1‑91 (MQ=255)
                                                          aTGGCGTGATAGAGCGTCGAAGCGAGGGACAGCAGGATCATACTGCCGCCATAGAGGCTGTAGCTGGTTATCGCTGTGGCGCTGGCATTaa                >  2:818460/1‑91 (MQ=255)
                                                           tGGCGTGATAGAGCGTCGAAGCGAGGGACGGCAGGATCAGACTGCCGCCATAGAGGCTGTAGGTGGTTATCGCTGTGGCGCTGGCATTAAg               <  2:752220/91‑1 (MQ=255)
                                                            ggCGTGATAGAGCGTCGAAGCGAGGGACAGCAGGATCATACTGCCGCCATAGAGGCTGTAGCTGGTTATCGCTGTGGCGCTGGCATTAAGa              >  2:435689/1‑91 (MQ=255)
                                                            ggCGTGATAGAGCGTCGAAGCGAGGGACAGCAGGATCATACTGCCGCCATAGAGGCTGTAGCTGGTTATCGCTGTGGCGCTGGCATTAAGa              <  1:351673/91‑1 (MQ=255)
                                                                 gATAGAGCGTCGAAGCGAGGGACAGCAGGATCATACTGCCGCCATAGAGGCTGTAGCTGGtt                                      >  2:201932/1‑62 (MQ=255)
                                                                 gATAGAGCGTCGAAGCGAGGGACAGCAGGATCATACTGCCGCCATAGAGGCTGTAGCTGGtt                                      <  1:201932/62‑1 (MQ=255)
                                                                 gATAGAGCGTCGAAGCGAGGGACAGCAGGATCATACTGCCGCCATAGAGGCTGTAGCTGGTTATCGCTGTGGCGCTGGCATTAAGATCCAc         <  2:308104/91‑1 (MQ=255)
                                                                 gATAGAGCGTCGAAGCGAGGGACAGCAGGATCATACTGCCGCCATAGAGGCTGTAGCTGGTTATCGCTGTGGCGCTGGCATTAAGATCCAc         <  2:412949/91‑1 (MQ=255)
                                                                 gATAGAGCGTCGAAGCGAGGGACAGCAGGATCATACTGCCGCCATAGAGGCTGTAGCTGGTTATCGCTGTGGCGCTGGCATTAAGATCCAc         <  2:108347/91‑1 (MQ=255)
                                                                 gATAGAGCGTCGAAGCGAGGGACAGCAGGATCATACTGCCGCCATAGAGGCTGTAGCTGGTTATCGCTGTGGCGCTGGCATTAAGATCCAc         <  2:105791/91‑1 (MQ=255)
                                                                 gATAGAGCGTCGAAGCGAGGGACAGCAGGATCATACTGCCGCCATAGAGGCTGTAGCTGGTTATCGCTGTGGCGCTGGCATTAAGATCCAc         <  2:808073/91‑1 (MQ=255)
                                                                       gCGTCGAAGCGAGGGACAGCAGGATCATACTGCCGCCATAGAGGCTGTAGCTGGTTATCGCTGTGGCGCt                        >  2:295513/1‑70 (MQ=255)
                                                                       gCGTCGAAGCGAGGGACAGCAGGATCATACTGCCGCCATAGAGGCTGTAGCTGGTTATCGCTGTGGCGCt                        <  1:295513/70‑1 (MQ=255)
                                                                       gCGTCGAAGCGAGGGACAGCAGGATCATACTGCCGCCATAGAGGCTGTAGCTGGTTATCGCTGTGGCGCTGGCATTAAGATCCACCGCCTg   >  1:258258/1‑91 (MQ=255)
                                                                        cGTCGAAGCGAGGGACAGCAGGATCATACTGCCGCCATAGAGGCTGTAGCTGGTTATCGCTGTGGCGCTGGCATTAAGATCCACCGCCTGa  <  2:211708/91‑1 (MQ=255)
                                                                                     |                                                                             
CAGGTAAATAGCGCAATGGTCAAATTTCTTCAGCCACATTTTTGCCCGTTGATGAGGAATGGCGTGATAGAGCGTCGAAGCGAGGAACAGCAGGATCATACTGCCGCCATAGAGGCTGTAGCTGGTTATCGCTGTGGCGCTGGCATTAAGATCCACCGCCTGA  >  NZ_CP009273/3036229‑3036391

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