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,3140AG97.3% 131.2 / NA 37F65S (TTC→TCC) yqfAhemolysin III family protein
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base G (11/25);  total (12/25)
Fisher's exact test for biased strand distribution p-value = 3.24e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AACAGGTAAATAGCGCAATGGTCAAATTTCTTCAGCCACATTTTTGCCCGTTGATGAGGAATGGCGTGATAGAGCGTCGAAGCGAGGAACAGCAGGATCATACTGCCGCCATAGAGGCTGTAGCTGGTTATCGCTGTGGCGCTGGCATTAAGATCCACCGCCTGAACCA  >  NZ_CP009273/3036227‑3036395
                                                                                       |                                                                                 
aaCAGGTAAATAGCGCAATGGTCAAATTTCTTCAGCCACATTTTTGCCCGTTGATGAGGAATGGCGTGATAGAGCGTCGAAGCGAGGGAca                                                                                <  1:598762/91‑1 (MQ=255)
aaCAGGTAAATAGCGCAATGGTCAAATTTCTTCAGCCACATTTTTGCCCGTTGATGAGGAATGGCGTGATAGAGCGTCGAAGCGAGGGAca                                                                                <  1:583517/91‑1 (MQ=255)
  cAGGTAAATAGCGCAATGGTCAAATTTCTTCAGCCACATTTTTGCCCGTTGATGAGGAATGGCGTGATAGAGCGTCGAAGCGAGGGAcagc                                                                              <  1:643620/91‑1 (MQ=255)
  cAGGTAAATAGCGCAATGGTCAAATTTCTTCAGCCACATTTTTGCCCGTTGATGAGGAATGGCGTGATAGAGCGTCGAAGCGAGGGAcagc                                                                              <  1:286184/91‑1 (MQ=255)
  cAGGTAAATAGCGCAATGGTCAAATTTCTTCAGCCACATTTTTGCCCGTTGATGAGGAATGGCGTGATAGAGCGTCGAAGCGAGGGAcagc                                                                              <  2:1022657/91‑1 (MQ=255)
    ggTAAATAGCGCAATGGTCAAATTTCTTCAGCCACATTTTTGCCCGTTGATGAGGAATGGCGTGATAGAGCGTCGAAGCGAGGGAcagcag                                                                            >  2:507190/1‑91 (MQ=255)
    ggTAAATAGCGCAATGGTCAAATTTCTTCAGCCACATTTTTGCCCGTTGATGAGGAATGGCGTGATAGAGCGTCGAAGCGAGGGAcagcag                                                                            <  2:724709/91‑1 (MQ=255)
      tAAATAGCGCAATGGTCAAATTTCTTCAGCCACATTTTTGCCCGTTGATGAGGAATGGCGTGATAGAGCGTCGAAGCGAGGGACAGCAGGa                                                                          <  1:460652/91‑1 (MQ=255)
            gcgcAATGGTCAAATTTCTTCAGCCACATTTTTGCCCGTTGATGAGGAATGGCGTGATAGAGCGTCGAAGCGAGGGACAGCAGGATCATAc                                                                    >  1:155748/1‑91 (MQ=255)
             cgcAATGGTCAAATTTCTTCAGCCACATTTTTGCCCGTTGATGAGGAATGGCGTGATAGAGCGTCGAAGCGAGGGACAGCAGGATCATACt                                                                   <  1:984916/91‑1 (MQ=255)
              gcAATGGTCAAATTTCTTCAGCCACATTTTTGCCCGTTGATGAGGAATGGCGTGATAGAGCGTCGAAGCGAGGGACAGCAGGATCATACTg                                                                  >  2:1095389/1‑91 (MQ=255)
              gcAATGGTCAAATTTCTTCAGCCACATTTTTGCCCGTTGATGAGGAATGGCGTGATAGAGCGTCGAAGCGAGGAACAGCAGGATCATACTg                                                                  >  1:1024367/1‑91 (MQ=255)
               cAATGGTCAAATTTCTTCAGCCACATTTTTGCCCGTTGATGAGGAATGGCGTGATAGAGCGTCGAAGCGAGGGACAGCAGGATCATACTgc                                                                 <  2:330592/91‑1 (MQ=255)
                 aTGGTCAAATTTCTTCAGCCACATTTTTGCCCGTTGATGAGGAATGGCGTGATAGAGCGTCGAAGCGAGGGACAGCAGGATCATACTgccg                                                               >  1:413966/1‑91 (MQ=255)
                 aTGGTCAAATTTCTTCAGCCACATTTTTGCCCGTTGATGAGGAATGGCGTGATAGAGCGTCGAAGCGAGGGACAGCAGGATCATACTgccg                                                               <  2:441476/91‑1 (MQ=255)
                         aTTTCTTCAGCCACATTTTTGCCCGTTGATGAGGAATGGCGTGATAGAGCGTCGAAGCGAGGGACAGCAGGATCATACTGCCGCCATagag                                                       <  1:53305/91‑1 (MQ=255)
                                          tttGCCCGTTGATGAGGAATGGCGTGATAGAGCGTCGAAGCGAGGGACAGCAGGATCATACTGCCGCCATAGAGGCTGTAGCTGGTTATCg                                      >  1:885566/1‑91 (MQ=255)
                                             gCCCGTTGATGAGGAATGGCGTGATAGAGCGTCGAAGCGAGGGACAGCAGGATCATACTGCCGCCATAGAGGCTGTAGCTGGTTATCGCtg                                   <  2:868877/91‑1 (MQ=255)
                                              cccGTTGATGAGGAATGGCGTGATAGAGCGTCGAAGCGAGGGACAGCAGGATCATACTGCCGCCATAGAGGCTGTAGCTGGTTATCGCtgt                                  <  2:48154/91‑1 (MQ=255)
                                                   tgatgaGGAATGGCTTGATAGAGCGTCGAAGGGAGGGACAGCAGGATCATACTGCCGCCATAGAGGCTGTAGCTGGTTATCGCTGTGgcgc                             <  2:393295/91‑1 (MQ=255)
                                                   tgatgaGGAATGGCGTGATAGAGCGTCGAAGCGAGGGACAGCAGGATCATACTGCCGCCATAGAGGCTGTAGCTGGTTATCGCTGTGgcgc                             <  2:68620/91‑1 (MQ=255)
                                                   tgatgaGGAATGGCGTGATAGAGCGTCGAAGCGAGGGACAGCAGGATCATACTGCCGCCATAGAGGCTGTAGCTGGTTATCGCTGTGgcgc                             <  1:30439/91‑1 (MQ=255)
                                                   tgatgaGGAATGGCGTGATAGAGCGTCGAAGCGAGGGACAGCAGGATCATACTGCCGCCATAGAGGCTGTAGCTGGTTATCGCTGTGgcgc                             <  2:918146/91‑1 (MQ=255)
                                                           aaTGGCGTGATAGAGCGTCGAAGCGAGGGACAGCAGGATCATACTGCCGCCATAGAGGCTGTAGCTGGTTATCGCTGTGGCGCTGGCATTa                     <  1:955647/91‑1 (MQ=255)
                                                            aTGGCGTGATAGAGCGTCGAAGCGAGGGACAGCAGGATCATACTGCCGCCATAGAGGCTGTAGCTGGTTATCGCTGTGGCGCTGGCATTaa                    >  1:72249/1‑91 (MQ=255)
                                                            aTGGCGTGATAGAGCGTCGAAGCGAGGGACAGCAGGATCATACTGCCGCCATAGAGGCTGTAGCTGGTTATCGCTGTGGCGCTGGCATTaa                    >  2:239318/1‑91 (MQ=255)
                                                              ggCGTGATAGAGCGTCGAAGCGAGGGACAGCAGGATCATACTGCCGCCATAGAGGCTGTAGCTGGTTATCGCTGTGGCGCTGGCATTAAGa                  >  1:177431/1‑91 (MQ=255)
                                                                 gTGATAGAGCGTCGAAGCGAGGGACAGCAGGATCATACTGCCGCCATAGAGGCTGTAGCTGGTTATCGCTGTGGCGCTGGCATTAAGATcc               <  1:272101/91‑1 (MQ=255)
                                                                   gATAGAGCGTCGAAGCGAGGGACAGCAGGATCATACTGCCGCCATAGAGGCTGTAGCTGGTTATCGCTGTGGCGCTGGCATTAAGATCCAc             >  1:600899/1‑91 (MQ=255)
                                                                    aTAGAGCGTCGAAGCGAGGGACAGCAGGATCATACTGCCGCCATAGAGGCTGTAGCTGGTTATCGCTGTGGCGCTGGCATTAAGATCCAcc            >  1:934617/1‑91 (MQ=255)
                                                                       gagCGTCGAAGCGAGGGACAGCAGGATCATACTGCCGCCATAGAGGCTGTAGCTGGTTATCGCTGTGGCGCTGGCATTAAGATCCACCGcc         >  2:1054023/1‑91 (MQ=255)
                                                                       gagCGTCGAAGCGAGGGACAGCAGGATCATACTGCCGCCATAGAGGCTGTAGCTGGTTATCGCTGTGGCGCTGGCATTAAGATCCACCGcc         <  1:1069412/91‑1 (MQ=255)
                                                                        agCGTCGAAGCGAGGGACAGCAGGATCATACTGCCGCCATAGAGGCTGTAGCTGGTTATCGCTGTGGCGCTGGCATTAAGATCCACCGCCt        <  2:155748/91‑1 (MQ=255)
                                                                        agCGTCGAAGCGAGGGACAGCAGGATCATACTGCCGCCATAGAGGCTGTAGCTGGTTATCGCTGTGGCGCTGGCATTAAGATCCACCGCCt        <  2:661131/91‑1 (MQ=255)
                                                                         gCGTCGAAGCGAGGGACAGCAGGATCATACTGCCGCCATAGAGGCTGTAGCTGGTTATCGCTGTGGCGCTGGCATTAAGATCCACCGCCTg       <  2:606172/91‑1 (MQ=255)
                                                                           gTCGAAGCGAGGGACAGCAGGATCATACTGCCGCCATAGAGGCTGTAGCTGGTTATCGCTGTGGCGCTGGCATTAAGATCCACCGCCTGaa     <  1:529337/91‑1 (MQ=255)
                                                                              gAAGCGAGGGACAGCAGGATCATACTGCCGCCATAGAGGCTGTAGCTGGTTATCGCTGTGGCGCTGGCATTAAGATCCACCGCCTGAACCa  <  1:311028/91‑1 (MQ=255)
                                                                                       |                                                                                 
AACAGGTAAATAGCGCAATGGTCAAATTTCTTCAGCCACATTTTTGCCCGTTGATGAGGAATGGCGTGATAGAGCGTCGAAGCGAGGAACAGCAGGATCATACTGCCGCCATAGAGGCTGTAGCTGGTTATCGCTGTGGCGCTGGCATTAAGATCCACCGCCTGAACCA  >  NZ_CP009273/3036227‑3036395

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