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,2980AT98.1% 113.9 / ‑5.9 36T15T (ACA→ACTfemA_3Aminoacyltransferase FemA
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base T (11/24);  minor base G (0/1);  total (11/25)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-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.

TGTATTTGTAATCAGTTTTTGATAAAATAAAAAAATAGCATTATTAATAAATTCAAAAATTATAGTTCAATTCAGAGAGGGATGGGATAAAATGAACTTTGTAACGTTGACTTCGGATGAGTTCAATGCGTTTACAACAAAGCATTTTTCACATTACACACAATCAGCTATTCATTACAATCATAGAGTTGATTTAAAAGGCGATGTGCATCTTGTAGGGGTTAAAGATGACAATGGTCAAGTGATTGCAGGATGCTTATTGACA  >  USA300TCH1516_ALE/2532163‑2532427
                                                                                                                                       |                                                                                                                                 
tgtATTTGTAATCAGTTTTTGATAAAATAAAAAAATAGCATTATTAATAAATTCAAAAATTATAGTTCAATTCAGAGAGGGATGGGATAAAATGAACTTTGTAACGTTGACTTCGGATGAGTTCAATGCGTTTACTACaaa                                                                                                                              <  2:10902/141‑1 (MQ=255)
  tATTTGTAATCAGTTTTTGATAAAATAAAAAAATAGCATTATTAATAAATTCAAAAATTATAGTTCAATTCAGAGAGGGATGGGATAAAATGAACTTTGTAACGTTGACTTCGGATGAGTTCAATGCGTTTACTACAAAGc                                                                                                                            <  2:80344/141‑1 (MQ=255)
   aTTTGTAATCAGTTTTTGATAAAATAAAAAAATACCATTATTCATAAATTCAAAAATTATAGTTCAATTCAGAGAGGGATGGCATAAAATGAACTTTGTAACGTTGACTTCGGATGAGTTCAATGCGTTTACTACAAAGCa                                                                                                                           <  1:73295/141‑1 (MQ=255)
     ttGTAATCAGTTTTTGATAAAATAAAAAAATAGCATTATTAATAAATTCAAAAATTATAGTTCAATTCAGAGAGGGATGGGATAAAATGAACTTTGTAACGTTGACTTCGGATGAGTTCAATGCGTTTACTACAAAGCAtt                                                                                                                         <  1:312299/141‑1 (MQ=255)
      tGTAATCAGTTTTTCAGAAAATAAAAAAATAGCATTATTAATAAATTCAAAAATTATAGTTCAATTCAGAGAGGGATGGGATAAAATGAACTTTGTAACGTTGACTTCGGACGAGTTCAATGCGTTTACTACAAAGCAttt                                                                                                                        <  2:535/141‑1 (MQ=255)
            cAGTTTTTGATAAAATAAAAAAATAGCATTATTAATAAATTCAAAAATTATAGTTCAATTCAGAGAGGGATGGGATAAAATGAACTTTGTAACGTTGACTTCGGATGAGTTCAATGCGTTTACTACAAAACATTTTTcaca                                                                                                                  >  1:19234/1‑141 (MQ=255)
            cAGTTATTGATAAAATAAAAAAATAGCATTATTAATGAATTCAAGAATTATAGGTCAATTAAGAGAGGGATAGGATAAAAGGGACGTGGTAAGGTTAGCTTCGGATGAGTTCAATGCGTTTACTACAAAGCAtttttacca                                                                                                                  >  2:90966/1‑137 (MQ=255)
                ttttGATAAAATAAAAAAATAGCATTATTAATAAATTCAAAAATTATAGTTCAATTCAGAGAGGGATGGGATAAAATGAACTTTGTAACGTTGACTTCGGATGAGTTCAATGCGTTTACTACAAAGCATTTTTCACATTac                                                                                                              >  2:226846/1‑141 (MQ=255)
                  ttGATAAAATAAAAAAATAGCATTATTAATAAATTCAAAAATTATAGTTCAATTCAGAGAGGGATGGGATAAAATGAACTTTGTAACGTTGACTTCGGATGAGTTCAATGCGTTTACTACAAAGCATTTTTCACATTacac                                                                                                            <  2:354684/141‑1 (MQ=255)
                   taaaaaaataaaaAAATAGCATTATTAATAAATTCAAAAATTATAGTTCAATTCAGAGAGGGATGGGATAAAATGAACTTTGTAACGTTGACTTCGGATGAGTTCAATGCGTTTACTACAAAGCATTTTTCACATTacaca                                                                                                           <  2:119562/137‑1 (MQ=255)
                      taaaataaaaAAATAGCATTATTAATAAATTCAAAGATGATAGGTCAATTGAGAGAGGGATGGGATAATATGAACTTTGTAACGTTGAGTGCGTATGAGTTCAATGTGGTTAGTACAAAGCATTGTTCACATTACACACaa                                                                                                        <  2:101570/141‑1 (MQ=255)
                         aataaaaAAATAGCATTATTAATAAATTCAAAAATTATAGTTCAATTCAGAGAGGGATGGGATAAAATGAACTTTGTAACGTTGACTTCGGATGAGTTCAATGCGTTTACTACAACGCCTTTTTCACATTACACACAATCa                                                                                                     <  2:392325/141‑1 (MQ=255)
                         aataaaaAAATAGCATTATTAATAAATTCAAAAATTATACTTCAATTCAGAGATGGATTGGATAAAATGAACTTTGTAACGTTGACTTCGGATGAGTTCAATGCGTTTACTACAAAGCATTTTTCACATTACACACAATCa                                                                                                     <  2:161532/141‑1 (MQ=255)
                          ataaaaAAATAGCATTATTAATAACTTCAAAAATTATAGTTCAATTCAGAGAGGGATGGGATAAAATGAACTTTGTAACGTTGACTTCGGATGAGTTCAATGCGTTTACTACAAAGCATTTTTCATATTACACACAATCAg                                                                                                    >  2:159988/1‑141 (MQ=255)
                            aaaaaaaTAGCATTATTAATAAATTCAAAAATTATAGTTCAATTCAGAGAGGGATGGGATAAAATGAACTTTGTAACGTTGACTTCGGATGAGTTCAATGCGTTTACTACAAAGCATTTTTCACATTACACACAATCAGCt                                                                                                  <  2:252781/141‑1 (MQ=255)
                             aaaaaaTAGCATTATTAATAAATTCAAAAATTATAGTTCAATTCAGAGAGGGATGGGATAAAATGAACTTTGTAACGTTGACTTCGGATGAGTTCAATGCGTTTACTACAAAGCATTTTTCACATTACACACAATCAGCTa                                                                                                 >  1:59528/1‑141 (MQ=255)
                                 aaTAGCATTATTAATAAATTCAAAAATTATAGTTCAATTCAGAGAGGGATGGGATAAAATGAACTTTGTAACGTTGACTTCGGATGAGTTCAATGCGTTTACTACAAAGCATTTTTCACATTACACACAATCAGCTATTCa                                                                                             <  2:138240/141‑1 (MQ=255)
                                       attattAATAAATTCAAAAATTATAGTTCAATTCAGAGAGGGATGGGATAAAATGAACTTTGTAACGTTGACTTCGGATGAGTTCAATGCGTTTACTACAAAGCATTTTTCACATTACACACAATCAGCTATTCATTACaa                                                                                       >  1:366175/1‑141 (MQ=255)
                                       attattAATAAATTCAAAAATTATAGTTCAATTCAGAGAGGGATGGGATAAAATGAACTTTGTAACGTTGACTTCGGATGAGTTCAATGCGTTTACTACAAAGCATTTTTCACATTACACACAATCAGCTAGTCATTACaa                                                                                       >  1:235159/1‑141 (MQ=255)
                                            taataaATTCAAAAATTATAGTTCAATTCAGAGAGGGATGGGATAAAATGAACTTTGTAACGTTGACTTCGGATGAGTTCAATGCGTTTACTACAAAGCATTTTTCACATTACACACAATCAGCTATTCATTACAATCATa                                                                                  <  1:99213/141‑1 (MQ=255)
                                                              taGTTCAATTCAGAGAGGGATGGGATAAAATGAACTTTGTAACGTTGACTTCGGATGAGTTCAATGCGTTTACTACAAAGCATTTTTCACATTACACACAATCAGCTATTCATTACAATCATAGAGTTGATTTAAAAGGCg                                                                <  1:348972/141‑1 (MQ=255)
                                                               aGTTCAATTCAGAGAGGGATGGGATAAAATGAACTTTGTAACGTTGACTTCGGATGAGTTCAATGCGTTTACTACAAAGCATTTTTCACATTACACACAATCAGCTATTCATTACAATCATAGAGTTGATTTAAAAGGCGa                                                               <  2:222962/141‑1 (MQ=255)
                                                                 ttCAATTCAGAGAGGGATGGGATAAAATGAACTTTGTAACGTTGACTTCGGATGAGTTCAATGCGTTTACTACAAAGGATTTTTCACAGTACACACAATCAGCTATTCATTACAATCATAGAGTTGATTTAAAAGGCGAtg                                                             <  1:210635/141‑1 (MQ=255)
                                                                  tCAATTCAGAGAGGGATGGGATAAAATGAACTTTGTAACGTTGACTTCGGATGAGTTCAATGCGTTTACTACAAAGCATTTTTCACATTACACACAATCAGCTATTCATTACAATCATAGAGTTGATTTAAAAGGCGAtgt                                                            >  1:376459/1‑141 (MQ=255)
                                                                       tCAGAGAGGGATGGGATAAAATGAACTTTGTAACGTCGACTTCGGATGAGTTCAATGCGTTTACTACAAAGCATTTTTCACATTACACACAATCAGCTATTCATTACAATCATAGAGTTGATTTAAAAGGCGATGTGCATc                                                       <  1:128889/141‑1 (MQ=255)
                                                                        cAGAGAGGGATGGGATAAAATGAACTTTGTAACGTTGACTTCGGATGAGTTCAATGCGTTTACTACAAAGCATTTTTCACATTACACACAATCAGCTATTCATTACAATCATAGAGTTGATTTAAAAGGCGATGTGCATCt                                                      <  1:136909/141‑1 (MQ=255)
                                                                                gatgggatAAAATGAACTTTGCAACGTTGACCTCGGATGAGTTCAATGCGTTTACTACAAAGCATTTTTCACATTACACACAATCAGCTATTCATTACAATCATAGAGTTGATTTAAAAGGCGATGTGCATCTTGTAgggg                                              <  1:275509/141‑1 (MQ=255)
                                                                                              aaCTTTGTAACGTTGACTTCGGATGAGTTCAATGCGTTTACTACAAACCATTTTTCACATTACACAAAATCAGCTATTCATTACAATCATAGAGTTGATTTAAAAGGCGATGTGCATCTTGTAGGGGTTAAAGATGACAAt                                <  2:101230/141‑1 (MQ=255)
                                                                                                      aaCGTTGACTTCGGATGAGTTCAATGCGTTTACTACAAAGCATTTTTCACATTACACACAATCAGCTATTCATTACAATCATAGAGTTGATTTAAAAGGCGATGTGCATCTTGTAGGGGTTAAAGATGACAATGGTCAAGt                        <  2:349965/141‑1 (MQ=255)
                                                                                                           tGACTTCGGATGAGTTCAATGCGTTTACGAGAAAGCATTTTTCACAGTACACACAATCAGCTATTCATTACAATCATAGAGTTGATTTAAAAGGCGATGTGCATCTTGTAGGGGTTAAAGATGACAATGGTCAAGTGATTg                   <  2:257854/141‑1 (MQ=255)
                                                                                                               ttCGGATGAGTTCAATGCGTTTACTACAAAGCATTTTTCACATTACACACAATCAGCTATTCATTACAATCATAGAGTTGATTTAAAAGGCGATGTGCATCTTGTAGGGGTTAAAGATGACAATGGTCAAGTGATTGCAgg               <  2:32080/141‑1 (MQ=255)
                                                                                                                 cGGATGAGTTCAATGCGTTTACTACAAAGCATTTTTCACATTACACACAATCAGCTATTCATTACAATCATAGGGTTGATTTAAAAGGCGATGTGCATCTTGTAGGGGTTAAAGATAACAATGGTCAAGTGATTGCAGGAt             >  1:101935/1‑141 (MQ=255)
                                                                                                                 cGGATGAGTTCAATGCGTTTACTACAAAGCATTTTTCACATTACACACAATCAGCTATTCATTACAATCATAGAGTTGATTTAAAAGGCGATGTGCATCTTGTAGGGGTTAAAGATGACAATGGTCAAGTGATTGCAGGAt             >  1:72146/1‑141 (MQ=255)
                                                                                                                   gATGAGTTCAATGCGTTTACTACAAAGCATTTTTCACATTACACACAATCAGCTATTCATTACAATCATAGAGTTGATTTAAAAGGCGATGTGCATCTTGTAGGGGTTAAAGATGACAATGGTCAAGTGATTGCAGGGTGc           >  2:59537/1‑141 (MQ=255)
                                                                                                                      gAGTTCAATGCGTTTACTACAAAGCATTTTTCACATTACATACAATCAGCTATTCATTACAATCATAGAGTTGATTTAAAAGGCGATGTGCATCTTGTAGGGGTTAAAGATGACAATGGTCAAGTGATTGCAGGATGCTTa        <  2:201651/141‑1 (MQ=255)
                                                                                                                            aaTGCGTTTACTACAAAGCATTTTTCACATTACACACAATCAGCTATTCATTACAATAATGGATGAGGGGTAAAAGGCGATGTGCATCTTGTAGGGGTTAAAGATGACAATGGTCAAGTGATTGCAGGATGCTTATTGACa  <  1:389069/141‑1 (MQ=255)
                                                                                                                                       |                                                                                                                                 
TGTATTTGTAATCAGTTTTTGATAAAATAAAAAAATAGCATTATTAATAAATTCAAAAATTATAGTTCAATTCAGAGAGGGATGGGATAAAATGAACTTTGTAACGTTGACTTCGGATGAGTTCAATGCGTTTACAACAAAGCATTTTTCACATTACACACAATCAGCTATTCATTACAATCATAGAGTTGATTTAAAAGGCGATGTGCATCTTGTAGGGGTTAAAGATGACAATGGTCAAGTGATTGCAGGATGCTTATTGACA  >  USA300TCH1516_ALE/2532163‑2532427

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: