Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 2,532,298 A→T 100% T15T (ACA→ACT femA_3 → Aminoacyltransferase FemA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE2,532,2980AT97.9% 127.7 / ‑4.9 45T15T (ACA→ACTfemA_3Aminoacyltransferase FemA
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base T (15/29);  minor base G (0/1);  total (15/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.43e-01
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

TTGTAATCAGTTTTTGATAAAATAAAAAAATAGCATTATTAATAAATTCAAAAATTATAGTTCAATTCAGAGAGGGATGGGATAAAATGAACTTTGTAACGTTGACTTCGGATGAGTTCAATGCGTTTACAACAAAGCATTTTTCACATTACACACAATCAGCTATTCATTACAATCATAGAGTTGATTTAAAAGGCGATGTGCATCTTGTAGGGGTTAAAGATGACAATGGTCAAGTGATTGCAGGATGCTTATTGACAGAAGC  >  USA300TCH1516_ALE/2532168‑2532432
                                                                                                                                  |                                                                                                                                      
ttGTAATCATTTTTTGATAAAATAAAAAAATAGCATTATTAATAAATTCAAAAATTATAGTTCAATTCAGAGAGGGATGGGTTAAAATTAACTTTGTAACGTTGATTTCGGATGAGTTCAATGCGTTTACTACAAAGCAtt                                                                                                                              <  1:112955/141‑1 (MQ=255)
   taataaGTTTTTAATAAAATAAAAAAATAGCATTATTAATAAATTCAAAAATTATAGTTCAATTCAGAGAGGGATGGGATAAAATGAACTTTGTAACGTTGACTTCGGATGAGTTCAATGCGTTTACTACAAAGCAttttt                                                                                                                           <  2:3337/141‑1 (MQ=255)
     aTCAGATTTTGATAAAATAAAAAAATAGCATTATTAATAAATTCAAAAATTATAGTTCAATTCAGAGAGGGATGGGATAAAATGAACTTTGTAACGTTGACTTCGGATGAGTTCAATGCGTTTACTACACAACActttccc                                                                                                                         >  2:100464/1‑138 (MQ=255)
      tCAGTTTTTGATAAAATAAAAAAATAGCATTATTAATAAATTCAAAAATTATAGTTCAATTCAGAGAGGGATGGGATAAAATGAACTTTGTAACGTTTACTTCGGATGAGTTCAAAGCGTTTACTACAAAGCATTTTTcac                                                                                                                        >  2:140126/1‑141 (MQ=255)
        aGTTTTTGATAAAATAAAAAAATAGCATTATTAATAAATTCAAAAATTATAGTTAAATTCAGAGAGGGATGGGATAAAATGAACTTTGTAACGTTGACTTAGGATGAGTTCAATACGATCACTACAAAGCATTTTTAACAt                                                                                                                      >  2:310537/1‑141 (MQ=255)
         gTTTTTGATAAAATAAAAAAATAGCATTATTAATAAATTCAAAAATTATAGTTCAATTCAGAGAGGGATGGGATAAAATGAACTTTGTAACGTTGACTTCGGATGAGTTCAATGCGTTTACTACAAAGCATTTTTCACAtt                                                                                                                     >  1:292805/1‑141 (MQ=255)
           ttttGATAAAATAAAAAAATAGCATTATTAATAAATTCAAAAATTATAGTTCAATTCAGAGAGGGATGGGATAAAATGAACTTTGTAACGTTGACTTCGGATGAGTTCAATGCGTTTACTACAAAGCATTTTTCACATTac                                                                                                                   <  1:44560/141‑1 (MQ=255)
           ttatGATAAAATAAAAAAATAGCATTATTGATAAATTCAAAAATTATAGTTCAATTCAGAGAGGGATGGGATAAAATGAACTTTGTAACGTTGACTTCGGATGAGTTCAATGCGTTTACTACAAAGCATTTTTCACATTac                                                                                                                   <  2:332462/138‑1 (MQ=255)
           tcctGACAAAATAAAAAAATAGCATTATTAATAAATTCAAAAATTATAGTTCAATTCAGAGAGGGATGGGATAAAATGAACTTTGTAACGTTGACTTCGGATGAGTTCAATGCGTTTACTACAAAGCATTTTTCACATTat                                                                                                                   <  2:319428/138‑2 (MQ=255)
               gATAAAATAAAAAAATAGCATTATTAATAAATTCAAAAATTATAGTTCAATTCAGAGAGGGATGGGATAAAATGAACTTTGTAACGTTGACTTCGGATGAGTTCAATGCGTTTACTACAAAGCATTTTTCACATTacacac                                                                                                               >  2:333542/1‑141 (MQ=255)
                 aaaaataaaaAAATAGCATTATTAATAAATTCAAAAATTATAGTTCAATTCAGAGAGGGATGGGATAAAATGAACTTTGTAACGTTGACTTCGGATGAGTTCAATGCGTTTACTACAAAGCATTTTTCACATTACACACaa                                                                                                             <  2:246725/140‑1 (MQ=255)
                  aaaataaaaAAATAGCATTATTAATAAATTCAAAAATTATAGTTCAATTCAGAGAGGGATGGGATAAAATGAACTTTGTAACGTTGACTTCGGATGAGTTCAATGCGTTTACTACAAAGCATTTTTCACATTACACACAAt                                                                                                            <  1:208991/141‑1 (MQ=255)
                   aaataaaaAAATAGCATTATTAATAAATTCAAAAATTATAGTTCAATTCAGAGAGGGATGGGATAAAATGAACTTTGTAACGTTGACTTCCGATGAGTTCAATGCATTTACTACAAACCATTTTTCACATTACACACAATc                                                                                                           >  2:186758/1‑141 (MQ=255)
                       aaaaaaaaAGCAATATTAATAAATTCAAAAATTATAGTTCAATTCAGAGAGGGATGGGATAAAATGAACTTTGTAACGTTGACTTCGGATGAGTTCAATGCGTTTACTACAAAGCATTTTTCACATTACACACAATCAGCt                                                                                                       <  1:286794/141‑1 (MQ=255)
                          aaaaTAGCATTATTAATAAATTCAAAAATTATAGTTCAATTCAGAGAGGGATGGGATAAAATGAACTTTGTAACGTTGACTTCGGATGAGTTCAATGCGTTTACTACAAAGCATTTTTCACATTACACACAATCAGCTAtt                                                                                                    <  2:230402/141‑1 (MQ=255)
                          aaaaTAGCATTATTAATAAATTCAAAAATTATAGTTCAATTAAGAGAGGGATGAGATAAAATGAACTTTGTAACGTTGACTTCGGATGAGTTCAATGCGTTTACTACAAAGCATTTTTCACATTCCACACAATCAGCTAtt                                                                                                    >  2:59242/1‑141 (MQ=255)
                             aTAGCATTATTAATAAATTCAAAAATTATAGTTCAATTCAGAGAGGGATGGGATAAAATGAACTTTGTAACGTTGACTTCGGATGAGTTCAATGCGTTTACTACAAAGCATTTTTCACATTACACACAATCAGCTATTCAt                                                                                                 <  1:357104/141‑1 (MQ=255)
                             aTAGCATTATTAATAAATTCAAAAATTATAGGTCACTTCAGAGAGAGATGGGATAAAATGAACTTTGTAACGTTGACTTCGGATGAGTTCAATGCGTTTACTACAAAGCATTTTTCACATTACACACAATCAGCTATTCAt                                                                                                 <  1:54400/141‑1 (MQ=255)
                               aGCATTATTAATAAATTCAAAAATTATAGTTCAATTCAGAGAGGGATGGGATAAAATGAACTTTGTAACGTTGACTTCGGATGAGTTCAATGCGTTTACTACAAAGCATTTTTCACATTACACACAATCAGCTATTCATTa                                                                                               >  1:97891/1‑141 (MQ=255)
                                tcATTATTAATAAATTCAAAAATTATAGTTCAATTCAGAGAGGGATGGGATAAAATGAACTTTGTAACGTTGACTTCGGATGAGTTCAATGCGTTTACTACAAAGCATTTTTCACATTACACACAATCAGCTATTCATTAc                                                                                              <  1:137466/140‑1 (MQ=255)
                                            aaTTCAAAAATTATAGTTCAATTCAGAGAGGGATGGGATAAAATGAACTTTGTAACGTTGACTTCGGATGAGTTCAATGCGTTTACTACAAAGCATTTTTCACATTACACACAATCAGCTATTCATTACAATCATAGAGtt                                                                                  >  1:223920/1‑141 (MQ=255)
                                            aaTTCAAAAATTATAGTTCAATTCAGAGAGGGATGGGATAAAATGAACTTTGTAACGTTGACTTCGGATGAGTTAAATTCTTTTAATTCACAAAAATTTTCACATTACACACAATCAGCTATTCTTTACAATCATAGAGtt                                                                                  <  2:180032/141‑1 (MQ=255)
                                                 aaaaaTTATAGTGCAATTCAGAGAGGGATGGGATAAAATGAACTTTGGAACGTTGACGTCGGATGAGTTCAATGCGTTGACTACAAAGCATTTTTCACATTACACACAATCAGCTATTCATTACAATCATAGAGTTGAttt                                                                             <  1:254297/141‑1 (MQ=255)
                                                         taGTTCAATTCAGAGAGGGATGGGATAAAATGAACTTTGTAACGTTGACTTCGGATGAGTTCAATGCGTTTACTACAAAGCATTTTTCACATTACACACAATCAGCTATTCATTACAATCATAGAGTTGATTTAAAAGGCg                                                                     <  2:186831/141‑1 (MQ=255)
                                                            ttCAATTGAGAGAGGGATGGGATAAAATGAACTTTGTAACGTTGACTTCGGATGAGTTCAATGCGTTTACTACAAAGCATTTTTCACATTACACACAATCAGCTATTCATTACAATCATAGAGTTGATTTAAAAGGCGAtg                                                                  <  2:288415/141‑1 (MQ=255)
                                                             tCAATTCAGAGAGGGATGGGATAAAATGAACTTTGTAACGTTGACTTCGGATGAGTTCAATGCGTTTACTACAAAGCATTTTTCACATTACACACAATCAGCTATTCATTACAATCATAGAGTTGATTTAAAAGGCGAtgt                                                                 <  1:362267/141‑1 (MQ=255)
                                                              cAATTCAGAGAGGGATGGGATAAAATGAACTTTGTAACGTTGACTTCGGATGAGTTCAATGCGTTTACTACAAAGCATTTTTCACATTACACACAATCAGCTATTCATTACAATCATAGAGTTGATTTAAAAGGCGAtgtg                                                                >  2:269165/1‑141 (MQ=255)
                                                                  tCAGAGAGGGATGGGATAAAATGAACTTTGTAACGTTGACTTCGGATGAGTTCAATGCGTTTACTACAAAGCATTTTTCACATTACACACAATCAGCTATTCATTACAATCATAGAGTTGATTTAAAAGGCGATGTGCATc                                                            <  1:176173/141‑1 (MQ=255)
                                                                   tagagagGGATGGGATAAAATGAACTTTGTAACGTTGACTTCGGATGAGTTCAATGCGTTTACTACAAAGCATTTTTCACATTACACACAATCAGCTATTCATTACAATCATAGAGTTGATTTAAAAGGCGATGTGCATCt                                                           <  2:10970/140‑1 (MQ=255)
                                                                       gagGGATGGGATAAAATGAACTTTGTAACGTTGACTTCGGATGAGTTCAATGCGTTTACTACAAAGCATTTTTCACATTACACACAAGCAGGGATTCATTACAATCAGAGAGTTGATTTAAAAGGCGATGTGCATCTTGTa                                                       <  2:67303/141‑1 (MQ=255)
                                                                        agGGATGGGATAAAATGAACTTTGTAACGTTGACTTCGGATGAGTTCAATGCGTTTACTACAAAGCATTTTTCACATTACACACAATCAGCTATTCATTACAATCATAGAGTTGATTTAAAAGGCGATGTGCATCTTGTAg                                                      <  1:245026/141‑1 (MQ=255)
                                                                         ggtatgggagAAGATGAATTGTGTAACGTTGAGTGCGGATGAGGTCAATGCGTTGACGTTAAAGTGTTTTTCACATTACACACAATCAGCTATTCATTACAATCATAGAGTTGATTTAAAAGGCGATGTGCATCTTGTAgg                                                     <  2:298552/138‑1 (MQ=255)
                                                                         gggatgggatAAAATGAACTTTGTAACGTTGACTTCGGATGAGTTCAATGCGTTTACTACAAAGCATTTTTCACATTACACACAATCAGCTATTCATTACAATCATAGAGTTGATTTAAAAGGCGATGTGCATCTTGTAgg                                                     <  1:190853/141‑1 (MQ=255)
                                                                          ggatgggatAAAATGAACTTTTTAACGTTGACTTCGGATGAGTTCAATGCGTTTACTACAAAGCATTTTTCACATTACACACAATCAGCTATTCATTACAATCATAGAGTTGATTTAAAAGGCGATGTGCa                                                              <  1:63584/131‑1 (MQ=255)
                                                                          ggatgggatAAAATGAACTTTTTAACGTTGACTTCGGATGAGTTCAATGCGTTTACTACAAAGCATTTTTCACATTACACACAATCAGCTATTCATTACAATCATAGAGTTGATTTAAAAGGCGATGTGCa                                                              >  2:63584/1‑131 (MQ=255)
                                                                             gtggacAAAATGAACTTTGTAACGTTGACTTCGGATGAGTTCAATGCGTTTACTACAAAGCATTTTTCACATTACACACAATCAGCTATTCATTACAATCATAGAGTTGATTTAAAAGGCGATGTGCATCTTGTAGGGGtt                                                 <  1:366981/139‑1 (MQ=255)
                                                                               ggatAAAATGAACTTTGTAACGTTGACTTCGGATGAGTTCAATGCGTTTACTACAAAGCATTTTTCACATTACACACAATCAGCTATTCATTACAATCATAGACTTGATTTCAACGTCGATGTGCATCTTGTAGGGGttca                                               >  2:192334/1‑139 (MQ=255)
                                                                                    atattaaCTTTGTAACGTTGACTTCGAATGAGTTCAATGCGTTTACTACAAAGCATTTTTCACATTACACACAATCAGCTATTCATTACAATCATAGAGTTGATTTAAAAGGCGATGTGCATCTTGTAGGGGTTAAAGATg                                          <  1:153850/136‑1 (MQ=255)
                                                                                               gTAACGTTGACTTCGGATGAGTTCAATGCGTTTACTACAAAGCATTTTTCACATTACACACAATCAGCTATTCATTACAATCATAGAGTTGATTTAAAAGGCGATGTGCATCTTGTAGGGGTTAAAGATGACAATGGTCaa                               <  2:335271/141‑1 (MQ=255)
                                                                                                      ttatttCGGATGAGTTCAATGCGTTTACTACAAAGCATTTTTCACATTACACACAATCAGCTATTCATTACAATCATAGAGTTGATTTAAAAGGCGATGTGCATCTTGTAGGGGTTAAAGATGACAATGGTCAAGTGATTg                        <  2:84080/137‑1 (MQ=255)
                                                                                                       caattCGGATGAGTTCACTGCGTCTCATAAAATGAGTTTTTCACATTACACACAATCAGCAATTCATTTCAATCATAGAGTCGATTTAAAAGGCGATGCGCATCTCGTAGGGGTCAAAGATGACAATCGTCAAGTGATTGc                       <  2:192962/138‑1 (MQ=255)
                                                                                                           tCGGATGAGTTCAATGCGTTTACTACAAAGCATTTTTCACATTACACACAATCAGCTATTCATTACAATCATAGAGTTGATTTAAAAGGCGATGTGCATCTTGTAGGGGTTAAATATGACAATGGTGAAGTGATTGCATGa                   >  2:377163/1‑141 (MQ=255)
                                                                                                            cGGATGAGTTCAATGCGTTTACTACAAAGCATTTTTCACATTACACACAATCAGCTATTCATTACAATCATAGAGTTGATTTAAAAGGCGATGTGCATCTTGTAGGGATTAAAGATGACAATGGTCAAGTGATTGCAggct                  >  1:302975/1‑139 (MQ=255)
                                                                                                                      cAATGCGTTTACTACAAAGCATTTTTCACATTACACACAATCAGCTATTCATTACAATCATAGAGTTGATTTAAACGGCGCTGTGCATCTTGTAGGGGTTAAAGATGACAATGGTCAAGAGATAGCAGGATGCTTATTGAc        >  1:192023/1‑141 (MQ=255)
                                                                                                                            gTTTACTACAAAGCATTTTTCACATTACACACAATCAGCTATTCATTACAATCATAGAGTTGATTTAAAAGGCGATGTGCATCTTGTAGGGGTTAAAGATGACAATGGTCAAGTGATTGCAGGATGCTTATTGACAGAAgc  <  1:324472/141‑1 (MQ=255)
                                                                                                                                  |                                                                                                                                      
TTGTAATCAGTTTTTGATAAAATAAAAAAATAGCATTATTAATAAATTCAAAAATTATAGTTCAATTCAGAGAGGGATGGGATAAAATGAACTTTGTAACGTTGACTTCGGATGAGTTCAATGCGTTTACAACAAAGCATTTTTCACATTACACACAATCAGCTATTCATTACAATCATAGAGTTGATTTAAAAGGCGATGTGCATCTTGTAGGGGTTAAAGATGACAATGGTCAAGTGATTGCAGGATGCTTATTGACAGAAGC  >  USA300TCH1516_ALE/2532168‑2532432

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