Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 1,252,446 A→G 13.4% K207K (AAA→AAG ffh → Signal recognition particle protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE1,252,4460AG13.4% 65.9 / 3.6 30K207K (AAA→AAGffhSignal recognition particle protein
Reads supporting (aligned to +/- strand):  ref base A (11/15);  new base G (2/2);  total (13/17)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.74e-01

GTGAAGGAGATCAAGTAAAGCCACAACAAATTGTAACTAATGCATTAAAACATGCTAAAGAAGAACATTTAGACTTTGTAATCATTGATACAGCAGGTCGATTACACATCGATGAAGCATTGATGAACGAATTAAAAGAAGTAAAAGACATTGCTAAACCAAACGAAATTATGTTAGTTGTCGATTCAATGACGGGTCAAGATGCTGTCAATGTTGCAGAATCTTTTGACGATCAACTTGATGTCACAGGTGTTACCTTAACTAAA  >  USA300TCH1516_ALE/1252310‑1252575
                                                                                                                                        |                                                                                                                                 
gTGAAGTAGAGAAAGTAAAGGCACAACAAATTGTAACTAATGCATTAAAACGTGCTAAAGAAGAACATTTAAACTTTATAATCATGTATACAGCAGGTTGATTACACATCGATGAAGCATTGATGAACGAATTAAaagaag                                                                                                                               <  2:187984/141‑1 (MQ=255)
    aggatctCAAGTAAACCCACAACAAATTGTAACTAATGCATTAAAACATGCTAAAGAAGAACATTTAGACTTTGTAATCATTGATACAGCATGTCGATTACACATCTATGAATCATTGATTAACGAATTAAAAGAAGTaaa                                                                                                                           <  1:115993/135‑1 (MQ=255)
    aGGAGATCAAGTAAAGCCACAACAAATTGTGACGATTCAAGCAGAAGATGCTAAAGAAGAACATTTAGACTTTGTAATCATTGATACAGCAGGTCGATTACACATCGATGAAGCATTGATGAACGAATTAAAAGAAGTaaa                                                                                                                           <  1:175602/141‑1 (MQ=255)
    aGGAGATCAAGTAAAGCCACAACAAATTGTAACTAATGCATTAAAACATGCTAAAGAAGAACATTTAGACTTTGTAATCATTGATACAGCAGGTCGATTACACATCGATGAAGCATTGATGAACGAATTAAAAGAAGTaaa                                                                                                                           <  1:312061/141‑1 (MQ=255)
      gagaTCAAGTAAAGCCACAACAAATTGTAACTAATGCATTAAAACATGCTAAAGAAGAACATTTAGACTTTGTAATCATTGATACAGCAGGTCGATTACACATCGATGAAGCATTGATGAACGAATTAAAAGAAGTAAAAg                                                                                                                         >  1:155228/1‑141 (MQ=255)
        gaTCAAGTAAAGCCACAACAAATTGTAACTAATGCATTAAAACATGCTAAAGAAGAACATTTAGACTTTGTAATCATTGATACAGCAGGTCGATTACACATCGATGAAGCATTGATGAACGAATTAAAAGAAGTAAAAGAc                                                                                                                       <  1:164216/141‑1 (MQ=255)
           cAAGTAAAGCCACAACAAATTGTAACTAATGCATTAAAACATGCTAAAGAAGAACATTTAGACTTTGTAATCATTGATACAGCAGGTCGATTACACATCGATGAAGCATTGATGAACGAATTAAAAGAAGTAAAAGACAtt                                                                                                                    >  2:38365/1‑141 (MQ=255)
             aGTAAAGCTACAGCAAAGTGTAACTAATGCATTAGAACATGCTAAAGAAGAACATTTAGACTTTGTAATCATTGATACAGCAGGTCGATTACACATCGATGAAGCATTGATGAACGAATTAAAAGAAGTAAAAGACATTGc                                                                                                                  <  2:155177/141‑1 (MQ=255)
             aGTAAAGCCACAACAAATTGTAACTAATGCATTAAAACATGCTAAAGAAGAACATTTAGACTTTGTAATCATTGATACAGCAGGTCGATTACACATCGATGAAGCATTGATGAACGAATTAAAAGAAGTAAAAGACATTGc                                                                                                                  >  2:13113/1‑141 (MQ=255)
              gTAAAGCCACAACAAATTGTAACTAATGCATTAAAACATGCTCCAGTCAAAGATTTAGACTTTGTAATCATTGATACAGCAGGTCGATTACACATCGATGAAGCATTGATGAACGAATTAAAAGAAGTAAAAGACATTGCt                                                                                                                 <  1:134117/141‑1 (MQ=255)
              gTAAAGCCACAACAAATTGTAACTAATGCATTAAAACATGCTAAAGAAGAACATTTAGACTTTGTAATCATTGATACAGCAGGTCGATTACACATCGATGAAGCATTGATGAACGAATTAAAGGAAGTAAAAGACATTGCt                                                                                                                 >  1:166048/1‑141 (MQ=255)
              gTAAAGCCACAACAAATTGTAACTAATGCATTAAAACATGCTAAAGAAGAACATTTAGACTTTGTAATCATTGATACAGCAGGTCGATTACACATCGATGAAGCATTGATGAACGAATTAAAGGAAGTAAAAGACATTGCt                                                                                                                 >  1:389158/1‑141 (MQ=255)
              gTAAAGACACAACAAATTGTTTCAATCTCATTAAAACATGCTAAAGAACAACATTTAGACTTTGTAATCATTGATACAGCAGGTCGATTACACATCGATGAAGCATTGATGAACGAATTAAAAGAAGTAAAAGACATTGCt                                                                                                                 >  2:403054/1‑141 (MQ=255)
                  aGCCACAACAAATTGTAACTAATGCATTAAAACATGCTAAAGAAGAACATTTAGACTTTGTAATCATTGATACAGCAGGTCGATTACACATCGATGAAGCATTGATGAACGAATTAAAAGAAGTAAAAGACATTGCTAAAc                                                                                                             >  1:331859/1‑141 (MQ=255)
                      accacaaTTTGAAACTAATGCAAGAAAACATGCGAAAGAAGAACATGTAGACTTTGTAATCATCGATACAGCAGGGCGATTACACATCGATGAAGCATTGATGAACGAATTAAAAGAAGTAAAAGACATTGCTAAACCaaa                                                                                                         <  1:127546/138‑1 (MQ=255)
                            aaTTGTAACTAATGCATTAAAACATGCTAAAGAAGAACATTTAGACTTTGTAATCATTGATACAGCAGGTCGATTACACATCGATGAAGCATTGATGAACGAATTAAAAGAAGTAAAAGACATTGCTAAACCAAACGAAAt                                                                                                   >  1:79108/1‑141 (MQ=255)
                               tGTAACTAATGCATTAAAACATGCTAAATAAAAACAGTTAGACTTTGTAATCATTGATACAGCAGGTCGATTACACATCGATGAAGCATTGATAAACGAATTAAAAGAAGTAAAAGACATTGCTAAACCAAACGAAATTAt                                                                                                >  2:90432/1‑141 (MQ=255)
                                        tGCATTAAAACATGCTAAAGAAGAACATTTAGACTTTGTAATCATTGATACAGCAGGTCGATTACACATCGATGAAGCATTGATGAACGAATTAAAAGAAGTAAAAGACATTGCTAAACCAAACGAAATTATGTTAGTTGt                                                                                       <  1:30678/141‑1 (MQ=255)
                                             tAAAACATGCTAAAGAAGAACATTTAGACTTTGTAATCATTGATACAGCAGGTCGATTACACATCGATGAAGCATTGATGAACGAATTAAAAGAAGTAAAAGACATTGCTAAACCAAACGAAATTATGTTAGTTGTCGAtt                                                                                  <  2:221356/141‑1 (MQ=255)
                                                       tAAAGAAGAACATTTAGACTTTGTAATCATTGATACAGCAGGTCGATTACACATCGATGAAGCATTGTTGAACGAATTAAAAGAAGTAAAAGACATTGCTAAACCAAACGAAGTTCTGTTAGATGTCGATTCATTGAAggg                                                                        >  2:9323/1‑141 (MQ=255)
                                                          agaagaACATTTAGACTTTGTAATCATTGATACAGCAGGTCGATTACACATCGATGAAGGATTGATGAACGAATTAAAAGAAGTAAAAGACATTGCTAAACCAAACGAAATTATGTTAGTTGTCGATTCAATGACGGGtaa                                                                     >  1:368538/1‑139 (MQ=255)
                                                                 cATTTAGACTTTGTAATGATTGATACAGAAGGTCGATTACACATAGAAGAAGCATTGATGAAAGAATTAAAAGAAGTAAAAGACATTGCTAAACCAAACGAAATTATGTTAGTTGTCGATTCAATGACGGGTCAAGATGCt                                                              <  2:197618/141‑1 (MQ=255)
                                                                                 ttaTTGATACAGTAGGTTGATTACACATAGATGAAGCATTGATGAACGAATTAAAAGTAGTAAAAGACATTGCTAAACCAAACGAATTTATGTTAGTTGTTGATTCAATGACGGGTCAAGATGCTGTCAATGTTGCAGAAt                                              <  2:300083/139‑1 (MQ=255)
                                                                                      gATACAGCAGGTCGATTACACATCGATGAAGCATTGATGAACGAATTAAAAGAAGTAAAAGACATTGCTAAACCAAACGAAATTATGTTAGTTGTCGATTCAATGACGGGTCAAGATGCTGTCAATGTTGCAGAATCtttt                                         >  1:399290/1‑141 (MQ=255)
                                                                                         aCAGCAGGTCGATTACACATCGATGAAGCATTGATGAACGAATTAAAAGAAGTAAAAGACATTGCTAAACCAAACGAAATTATGTTAGTTGTCGATTCAATGACGGGTCAAGATGCTGTCAATGTTGCAGAATCTTTTGAc                                      >  1:352500/1‑141 (MQ=255)
                                                                                           agcagGCCAATTACACATCGATGAAGCATTGATGAACGAATTAAAAGAAGTAAAAGACAATGCTAAACCAAACGAAATTTTGTTAGTTGTCGATTCAATGACGGGTCAAGATGCTGTCAATGTTGCAGAATCTTTTGACGa                                    <  1:236996/141‑1 (MQ=255)
                                                                                                     ttACACATCGATGAAGCATTGATGAACGAATTAAAGGAAGTAAAAGACATTGCTAAACCAAACGAAATTATGTTAGTTGTCGATTCAATGACGGGTCAAGATGCTGTCAATGTTGCAGAATCTTTTGACGATCAACTTGAt                          <  2:166048/141‑1 (MQ=255)
                                                                                                     atACACATCGATGAAGCATTGATGAACGAATTAAAGGAAGAAAAAGAAATTGCTAAACCAAGCGAAATTATGTTAGTTGTCGATTCAATGACGGGTCAAGACTCTGTCAATGGTGCAGAATCTTTTGACGATCAACTTGAt                          <  2:389158/140‑1 (MQ=255)
                                                                                                          caTGGAGGAAGGATTGATGATGGAATTAAAAGAAGTAAAAGACATTGCTAAACAACACGAAATTATGTTAGTTGTCGAGTCAATGACGGGTCAAGATGCTGTCAATGTTGCAGAATCTTTTGACGATCAACTTGATGTcac                     <  2:122894/141‑1 (MQ=255)
                                                                                                                     cATTGATGAACGAATTAAAAGAAGTAAAAGACATTGCTAAACCAAACGAAATTATGTTAGTTGTCGATTCAATGACGGGTCAAGATGCTGTCAATGTTGCAGAATCTTTTGACGATCAACTTGATGTCACAGGTGTTACCt          <  1:381168/141‑1 (MQ=255)
                                                                                                                             aaCGAATTAAAAGAAGTAAAAGACATTGCTAAACCAAACGAAATTATGTTAGTTGTCGATTCAATGACGGGTCAAGATGCTGTCAATGTTGCAGAATCTTTTGACGATCAACTTGATGTCACAGGTGTTACCTTAACTaaa  <  1:163398/141‑1 (MQ=255)
                                                                                                                                        |                                                                                                                                 
GTGAAGGAGATCAAGTAAAGCCACAACAAATTGTAACTAATGCATTAAAACATGCTAAAGAAGAACATTTAGACTTTGTAATCATTGATACAGCAGGTCGATTACACATCGATGAAGCATTGATGAACGAATTAAAAGAAGTAAAAGACATTGCTAAACCAAACGAAATTATGTTAGTTGTCGATTCAATGACGGGTCAAGATGCTGTCAATGTTGCAGAATCTTTTGACGATCAACTTGATGTCACAGGTGTTACCTTAACTAAA  >  USA300TCH1516_ALE/1252310‑1252575

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