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% 51.3 / NA 17F65S (TTC→TCC) yqfAhemolysin III family protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (7/10);  total (7/10)

TCCAGCCCCACCAGCAAAAACGGCGTGTAGGTTCCGGCAATCAACAGGTAAATAGCGCAATGGTCAAATTTCTTCAGCCACATTTTTGCCCGTTGATGAGGAATGGCGTGATAGAGCGTCGAAGCGAGGAACAGCAGGATCATACTGCCGCCATAGAGGCTGTAGCTGGTTATCGCTGTGGCGCTGGCATTAAGATCCACCGCCTGAACCAGTAGCAACACCAGCCCAACGATACCAAACACCAACCCAATGCCGT  >  NZ_CP009273/3036185‑3036440
                                                                                                                                 |                                                                                                                              
tCCAGCCCCACCAGCAAAAACGGCGTGTAGGTTCCGGCAATCAACAGGTAAATAGCGCAACGGTCAAATTTGTTCAGCCACATTTTTGCCCGTTGATGATGAATGGCGTGATAGAGCGTCGAAACGGGGGACAGCAGGa                                                                                                                       >  1:413595/1‑139 (MQ=255)
   aGCCCCACCAGCAAAAACGGCGTGTAGGATCCGGCAATCAACAGGTAAATAGCGCAATGGTCAAATTTCTTCAGCCACATTTTTGCCCGTTGATGAGGAATGGCGTGATAGAGCGTCGAAGCGAGGGACAGCAGGATCa                                                                                                                    >  1:348661/1‑139 (MQ=255)
        caccaGCAAAAACGGCGTGTAGGTTCCGGCAATCAACAGGTAAATAGCGCAATGGTCAAATTTCTTCAGCCACATTTTTGCCCGTTGATGAGGAATGGCGTGATAGAGCGTCGAAGCGAGGGACAGCAGGATCATACTg                                                                                                               <  2:259020/139‑1 (MQ=255)
                                              ggTAAATAGCGCAATGGTCAAATTTCTTCAGCCACATTTTTGCCCGTTGATGAGGAATGGCGTGATAGAGCGTCGAAGCGAGGGACAGCAGGATCATACTGCCGCCATAGAGGCTGTAGCTGGTTATCGCTGTGGCGCt                                                                         <  2:314036/139‑1 (MQ=255)
                                                         cAATGGTCAAATTTCTTCAGCCACATTTTTGCCCGTTGATGAGGAATGGCGTGATAGAGCGTCGAAGCGAGGGACAGCAGGATCATACTGCCGCCATAGAGGCTGTAGCTGGTTATCGCTGTGGCGCTGGCAtt                                                                   >  2:320936/1‑134 (MQ=255)
                                                            tGGTCAAATTTCTTCAGCCACATTTTTGCCCGTTGATGAGGAATGGCGTGATAGAGCGTCGAAGCGAGGGACAGCAGGATCATACTGCCGCCATAGAGGCTGTAGCTGGTTATCGCTGTGGCGCTGGCATTAAGATCCa                                                           <  2:171810/139‑1 (MQ=255)
                                                                          cAGCCACATTTTTGCCCGTTGATGAGGAATGGCGTGATAGAGCGTCGAAGCGAGGGACAGCAGGATCATACTGCCGCCATAGAGGCTGTAGCTGGTTATCGCTGTGGCGCTGGCATTAAGATCCACCGCCTGAAcc                                                <  2:413595/136‑1 (MQ=255)
                                                                          cAGCCACATTTTTGCCCGTTGATGAGGAATGGCGTGATAGAGCGTCGAAGCGAGGGACAGCAGGATCATACTGCCGCCATAGAGGCTGTAGCTGGTTATCGCTGTGGCGCTGGCATTAAGATCCACCGCCTGAACCa                                               <  2:36525/137‑1 (MQ=255)
                                                                                              gatgaGGAATGGCGTGATAGAGCGTCGAAGCGAGGGACAGCAGGATCATACTGCCGCCATAGAGGCTGTAGCTGGTTATCGCTGTGGCGCTGGCATTAAGATCCACCGCCTGAACCAGTAGCAACACCAGCCCAACGa                          <  1:238436/138‑1 (MQ=255)
                                                                                                  aGGAATGGCGTGATAGAGCTTCGAAGCGAGGGACAGCAGGATCATACTGCCGCCATAGAGGCTGTAGCTGGTTATCGCTGTGGCGCTGGCATTAAGATCCACCGCCTGAACCAGTAGCAACACCAGCCCAACGATACCa                     <  2:348661/139‑1 (MQ=255)
                                                                                                  aGGAATGGCGTGATAGAGCGTCGAAGCGAGGGACAGCAGGATCATACTGCCGCCATAGAGGCTGTAGCTGGTTATCGCTGTGGCGCTGGCATTAAGATCCACCGCCTGAACCAGTAGCAACACCAGCCCAACGATACCa                     <  1:282639/139‑1 (MQ=255)
                                                                                                         gCGTGATAGAGCGTCGAAGCGAGGGACAGCAGGATCATACTGCCGCCATAGAGGCTGTAGCTGGTTATCGCTGTGGCGCTGGCATTAAGATCCACCGCCTGAACCAGTAGCAACACCAGCCCAACGATACCAAACAcc               >  1:64254/1‑138 (MQ=255)
                                                                                                                    cGTCGAAGCGAGGGACAGCAGGATCATACTGCCGCCATAGAGGCTGTAGCTGGTTATCGCTGTGGCGCTGGCATTAAGATc                                                             >  1:69787/1‑81 (MQ=255)
                                                                                                                    cGTCGAAGCGAGGGACAGCAGGATCATACTGCCGCCATAGAGGCTGTAGCTGGTTATCGCTGTGGCGCTGGCATTAAGATc                                                             <  2:69787/81‑1 (MQ=255)
                                                                                                                    cGTCGAAGCGAGGGACAGCAGGATCATACTGCCGCCATAGAGGCTGTAGCTGGTTATCGCTGTGGCGCTGGCATTAAGATCCACCGCCTGAAcc                                                <  2:147305/94‑1 (MQ=255)
                                                                                                                    cGTCGAAGCGAGGGACAGCAGGATCATACTGCCGCCATAGAGGCTGTAGCTGGTTATCGCTGTGGCGCTGGCATTAAGATCCACCGCCTGAAcc                                                >  1:147305/1‑94 (MQ=255)
                                                                                                                      tCGAAGCGAGGGACAGCAGGATCATACTGCCGCCATAGAGGCTGTAGCTGGTTATCGCTGTGGCGCTGGCATTAAGATCCACCGCCTGAACCAGTAGCAACACCAGCCCAACGATACCAAACACCAACCCAATGCCGt  >  2:203068/1‑138 (MQ=255)
                                                                                                                                 |                                                                                                                              
TCCAGCCCCACCAGCAAAAACGGCGTGTAGGTTCCGGCAATCAACAGGTAAATAGCGCAATGGTCAAATTTCTTCAGCCACATTTTTGCCCGTTGATGAGGAATGGCGTGATAGAGCGTCGAAGCGAGGAACAGCAGGATCATACTGCCGCCATAGAGGCTGTAGCTGGTTATCGCTGTGGCGCTGGCATTAAGATCCACCGCCTGAACCAGTAGCAACACCAGCCCAACGATACCAAACACCAACCCAATGCCGT  >  NZ_CP009273/3036185‑3036440

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