Predicted mutation
evidence seq id position mutation annotation gene description
RA USA300TCH1516_ALE 2,532,298 A→T 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.6% 125.2 / NA 41T15T (ACA→ACT) femA_3Aminoacyltransferase FemA
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base T (15/25);  minor base C (0/1);  total (15/26)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCAGTTTTTGATAAAATAAAAAAATAGCATTATTAATAAATTCAAAAATTATAGTTCAATTCAGAGAGGGATGGGATAAAATGAACTTTGTAACGTTGACTTCGGATGAGTTCAATGCGTTTACAACAAAGCATTTTTCACATTACACACAATCAGCTATTCATTACAATCATAGAGTTGATTTAAAAGGCGATGTGCATCTTGTAGGGGTTAAAGATGACAATGGTCAAGTGATTGCAGGATGCTTATTGACAGAA  >  USA300TCH1516_ALE/2532174‑2532430
                                                                                                                            |                                                                                                                                    
taaGTTTTTGATAAAATAAAAAAATAGCATTATTAATAAATTCAAAAATTATAGTTCAATTCAGAGAGGGATGGGATAAAATGAACTTTGTAACGTTGACTTCGGATGAGTTCAATGCGTTTACTACAAAGCATTTTTcac                                                                                                                      <  2:14651/139‑1 (MQ=255)
tCAGTTTTTGATAAAATAAAAAAATAGCATTATTCATAAATTCAAAAATTATAGTTCAATTCAGAGAGGGATGGGATAAAATGAACTTTGTAACGTTGACTTCGGATGAGTTCAATGCGTTTACTACAAAGCATTTTTcac                                                                                                                      <  1:241536/141‑1 (MQ=255)
  aGTTTTTGATAAAATAAAAAAATAGCATTATTAATAAATTCAAAAATTATAGTTCAATTCAGAGAGGGATGGGATAAAATGAACTTTGTAACGTTGACTTCGGATGAGTTCAATGCGTTTACTACAAAGCATTTTTCACAt                                                                                                                    <  2:281998/141‑1 (MQ=255)
  aGTTTTTGATAAAATAAAAAAATAGCATTATTAATAAATTCAAAAATTATAGTTCAATTCAGAGAGGGATGGGATAAAATGAACTTTGTAACGTTGACTTCGGATGAGTTCAATGCGGTTACTACAAAGCATTTTTCACAt                                                                                                                    <  1:231887/141‑1 (MQ=255)
   gTTTTTGATAAAATAAAAAAATAGCATTTTTAATAAATTCAAAAATTATAGTTCAATTCAGAGAGGGATGGGATAAAATGCACTGTGTAACATGTAGTTCGGATGAGTTCAATGACTTTACTACAACGCATTTTTTACAtt                                                                                                                   <  2:160374/141‑1 (MQ=255)
      tttGATAAAATAAAAAAATAGCATTATTAATAAATTCAAAAATTATAGTTCAATTCAGAGAGGGATGGGATAAAATGAACTTTGTAACGTTGACTTCGGATGAGTTCAATGCGTTTACTACAAAGCCTTTGTAGCATTaca                                                                                                                >  2:378108/1‑141 (MQ=255)
       ttGAGAAAATAAAAAAATAGCATTATTAATAAACTCAAAAATTATAGTTCAATTCAGAGAGGGATGGGATAAAATGAACTTTGTAACGTTGACTTCGGATGAGTTCAATGCGTTTACTACAAAGCATTTTTCACATTacac                                                                                                               <  1:75372/141‑1 (MQ=255)
               agaaaaaaaTAGCAGTATTAATAAATTCAAAAATTATAGTTCAATTCAGAGAATGATGTTATAAAATGAACTTTGTAACGTTGACTTCGGATGAGTTCAATGCGTTTACTACAAAGCATTTTTCACATTACACACAATCAg                                                                                                       <  2:110574/139‑1 (MQ=255)
                aaaaaaaaGAGCATTATTAATAAATTCAAAAATTATAGTTCAAGTCAGAGAGGGATGGGATAAAATGAACTTTGTAACGTTGACTTCGGATGAGTTCAATGCGTTTACTACAAAGCATTTTTCACATTACACACAATCAGc                                                                                                      <  2:143689/140‑1 (MQ=255)
                 aaaaaaaTAGCATTATTAATAAATTCAAAAATTATAGTTCAATTCAGAGAGGGATGGGATAAAATGAACTTTGTAACGTTGACTTCGGATGAGTTCAATGCGTTTACTACAAAGCATTTTTCACATTACACACAATCAGCt                                                                                                     <  2:218313/141‑1 (MQ=255)
                   aaaaaTAGCATTATTAATAAATTCAAAAATTATAGTTCAATTCAGAGAGGGATGGGATAAAATGAACTTTGTAACGTTGACTTCGGATGAGTTCAATGCGTTTACTACAAAGCATTTTTCACATTACACACAATCAGCTAt                                                                                                   >  1:43501/1‑141 (MQ=255)
                                      aaTTCAAAAATTATAGTTCAATTCAGAGAGGGATGGGATAAAATGAACTTTGTAACGTTGACTTCGGATGAGTTCAATGCGTTTACTACAAAGCATTTTTCACATTACACACAATCAGCTATTCATTACAATCATAGAGtt                                                                                >  1:329985/1‑141 (MQ=255)
                                            aaaaTTATAGTTCAATTCAGAGAGGGATGGGATAAAATGAACTTTGTAACGTTGACTTCGGATGAGTTCAATGCGTTTACTACAAAGCATTTTTCACATTACACACAATCAGCTATTCATTACAATCATAGAGTTGATTTa                                                                          >  1:438720/1‑141 (MQ=255)
                                            aaaaTTATAGTTCAATTCAGAGAGGGATGGGATAAAATGAACTTTGTAACGTTGACTTCGGATGAGTTCAATGCGTTTACTACAAAGCATTTTTCACATTACACACAATCAGCTATTCATTACAATCATAGAGTTGATTTa                                                                          >  1:400182/1‑141 (MQ=255)
                                                ttATCGTTCAATTCAGAGAGGTATGGGATAAAATGAACTTTGTAACGTTGACTTCGGATGAGTTCAATGCGTTTACTACAAAGCATTTTTCACATTACACACAATCAGCTATTCATTACAATCATAGAGTTGATTTAAAAg                                                                      <  1:354319/141‑1 (MQ=255)
                                                  ataGTTCAATTCAGAGAGGGATGGGATAAAATGAACTTTGTAACGTTGACTTCGGATGAGTTCAATGCGTTTACTACAAAGCATTTTTCACATTACACACAATCAGCTATTCATTACAATCATAGAGTTGATTTAAAAGGc                                                                    >  1:70127/1‑141 (MQ=255)
                                                     gTTCAATTCAGAGAGGGATGGGATAAAATGAACTTTGTAACGTTGACTTCGGATGAGTTCAATGCGTTTACTACAAAGCATTTTTCACATTACACACAATCAGCTATTCATTACAATCATAGAGTTGATTTAAAAGGCGAt                                                                 >  1:402798/1‑141 (MQ=255)
                                                      ttCAATTCATATAGGTATGGGATAAAATGAACTTTGTAACGTTGACTTCGGATGAGTTCAATGCGTTTACTACAAAGCATTTTTCACATTACACACAATCAGCTATTCATTACAATCATAGAGTTGATTTAAAAGGCGAtg                                                                <  1:396523/141‑1 (MQ=255)
                                                             aaaatagGGATGGGATAAAATGAACTTTGTAACGTTGACTTCGGATGAGTTCAATGCGTTTACTACAAAGCATTTTTCACATTACACACAATCAGCTATTCATTACAATCATAGAGTTGATTTAAAAGGCGATGTGCATCt                                                         <  2:375980/136‑1 (MQ=255)
                                                              agagagGGATGGGATAAAATGAACTTTGTAACGTTGACTTCGGATGAGTTCAATGCGTTTACTACAAAGCATTTTTCACATTACACACAATCAGCTATTCATTACAATCATAGAGTTGATTTAAAAGGCGATGTGCATCtt                                                        <  1:315232/141‑1 (MQ=255)
                                                               gagagGGATGGGATAAAATGAACTTTGTAACGTTGACTTCGGATGAGTTCAATGCGTTTACTACAAAGCATTTTTCACATTACACACAATCAGCTATTCATTACAATCATAGCGTTGATTTAAAAGGCGATGTGCAACTTg                                                       >  1:99239/1‑141 (MQ=255)
                                                               gagagGGATGGGATAAAATGAACTTTGTAACGTTGACTTCGGATGAGTTCAATGCGTTTACTACAAAGCATTTTTCACATTACACACAATCAGCTATTCATCACAATCATAGAGTTTATTTAAAAGGCGATGTGCATCTTg                                                       >  1:166936/1‑141 (MQ=255)
                                                                  agGGATGGGATAAAATGAACTTTGTAACGTTGACTTCGGATGAGTTCAATGCGTTTACTACAAAGCATTTTTCACATTACACACAATCAGCTATTCATTACAATCATAGAGTTGATTTAAAAGGCGATGTGCATCTTGTAg                                                    <  1:421769/141‑1 (MQ=255)
                                                                    gtatgggatAAAATGAACTTTGTAACGTTAACTTCGGATGAGTTCAATGCGTTTACTACAAAGCATTTTTCACATTACACACAATCAGCTATTCATTACAATCATAGAGTTGATTTAAAAGGCGATGTGCATCTTGTAggg                                                  <  1:386244/139‑1 (MQ=255)
                                                                    ggatgggatAAAATGAACTTTGTAACGTTGACTTCGGATGAGTTCAATGCGTTTACTACAAAGCATTTTTCACATTACACACAATCAGCTATTCATTACAATCATAGAGTTGATTTAAAAGGCGATGTGCATCTTGTAggg                                                  <  2:383166/141‑1 (MQ=255)
                                                                    ggataggataAAATGAACTTTGTAACGTTGACTTCGGATGAGTTCAATGCGTTTACTACAAAGCATTTTTCACATTACACACAATCAGCTATTCATTACAATCATAGAGTTGATTTAAAAGGCGATGTGCATCTTGTAggg                                                  <  1:2819/141‑1 (MQ=255)
                                                                       tgggatAAAATGAACTTTGTAACGTTGACTTCGGATGAGTTCAATGCGTTTACTACAAAGCATTTTTCACATTACACACAATCAGCTATTCATTACAATCATAGAGTTGATTTAAAAGGCGATGTGCATCTTGTAGGGGtt                                               >  1:27290/1‑141 (MQ=255)
                                                                       tgggatAAAATGAACTTTGTAACGTTGACTTCGGATGAGTTCAATGCGTTTACTACAAAGCATTTTTCACATTACACACAATCAGCTATTCATTACAATCATAGAGTTGATTTAAAAGGCGATGTGCATCTTGTAGGGGtt                                               <  1:111284/141‑1 (MQ=255)
                                                                                 tGAACTTTGTAACGTTGACTTCGGATGAGTTCAATGCGTTTACTACAAAGCATTTTTCACATTACACACAATCAGCTATTCATTACAATCATAGAGTTGATTTAAAAGGCGATGTGCATCTTGTAGGGGTTAAAGATGACa                                     <  2:75283/141‑1 (MQ=255)
                                                                                 tGAACTTTGTAACGTTGACTTCGGATGAGTTCAATGCGTTTACTACAAAGCAGTTTTCACATTACACACAATCAGCTATTCATTACAATCATAGAGTTGATTTAAAAGGCGATGTGCATCTTGTAGGGGTTAAAGATGACa                                     <  1:278624/141‑1 (MQ=255)
                                                                                  gAACTTTGTAACGTTGACTTCGGATGAGTTCAATGCGTTTACTACAAAGCATTTTTCACATTACACACAATCAGCTATTCATTACAATCATAGAGTTGATTTAAAAGGCGATGTGCATCTTGTAAGGGTTAAAGATGACaa                                    >  1:354771/1‑141 (MQ=255)
                                                                                     cTTTGTAACGTTGACTTCGGATGAGTTCAATGCGTTTACTACAAAGCATGTTTCACATTACACACAATCAGCTATTCATTACAATCATAGAGTTGATTTAAAAGGCGATGTGCATCTTGTAGGGGTTAAAGATGACAATgg                                 <  2:91847/141‑1 (MQ=255)
                                                                                          tAACGTTGACTTCGGATGAGTTCAATGCGTTTACTACAAAGCATTTTTCACATTACACACAATCAGCTATTCATTACAATCATAGAGTTGATTTAAAAGGCGATGTGCATCTTGTAGGGGTTAAAGATGACAATGGTCAAg                            <  1:438914/141‑1 (MQ=255)
                                                                                                   cTTCGGAGGAGTTCAATGCGTTTACTACAAAGCATTTCTCACATTACACCCAATCAGCTATTCATTACAATCATAGAGGTTATTTAAAAGGCGATGTGCATCTTGTAGGGGTTAAAGATGACAATGGTCAAGTGATTGCAg                   <  1:216292/141‑1 (MQ=255)
                                                                                                    ttCGGATGAGTTCAATGCGTCTACTACAAAGCATTTTTCACATTACACACAATCAGCTATTCATTACAATCATAGAGTTGATTTAAAAGGCGATGTGCATCTTGTAGGGGTTAAAGATGACAATGGTCAAGTGATTGCAgg                  <  1:165132/141‑1 (MQ=255)
                                                                                                    tagcgaggAGTTCAATGCGTTTACCACAAATCACTTTGCACACTACACACAATCAGCTATTCATTACAATCATAGAGTTGATTTAAAAGGCGATGTGCATCTTGTAGGGGTTAAAGATGACAATGGTCAAGTGATTGCAgg                  <  1:196299/134‑1 (MQ=255)
                                                                                                     tCGGATGAGTTCAATGCGTTTACTACAAAGCATTTTTCACATTACACACAATCAGCTATTCATTACAATCATAGAGTTGATTTAAAAGGCGCTGTGCATCTTGTAGGGGTTAAAGATGACAATGGTCAAGTGATTGCAGGa                 >  2:22696/1‑141 (MQ=255)
                                                                                                       ggATGAGTTCAATGCGTTTACTACAAAGCATTTTTCACATTACACACAATCAGCTATTCATTACAATCATAGAGTTGATTTAAAAGGCGATGTGCATCTTGTAGGGGTTAAAGATGACAATGGTCAAGTAATTTCAGGATg               >  2:320212/1‑141 (MQ=255)
                                                                                                             gTTCAATGCGTTTACTACAAAGCATTTTTCACATTACACACAATCAGCTATTCATTACAATCATAGAGTTGATTTAAAAGGCGATGTGCATCTTGTAGGGGTTAAAGATGACAATGGTCAAGCAATTGCAGGATGCTTAtt         >  2:203975/1‑141 (MQ=255)
                                                                                                                 aaTGCGTTTACTACAAAGCATTTTTCACATTACACACAATCAGCTATTCATTACAATCATAGAGTTGATTTAAAAGGCGATGTGCATCTTGTAGGGGTTAAAGATGACAATGGTCAAGTGATTGCAGGATGCTTATTGACa     >  1:313108/1‑141 (MQ=255)
                                                                                                                    tcGTTTACTACAAAGCATTTTTCACATTACACACAATCAGCTATTCATTACAATCATAGAGTTGATTTAAAAGGCGATGTGCATCTTGTAGGGGTTAAAGATGACAATGGTCAAGTGATTGCAGGATGCTTATTGACAGaa  <  2:85739/140‑1 (MQ=255)
                                                                                                                            |                                                                                                                                    
TCAGTTTTTGATAAAATAAAAAAATAGCATTATTAATAAATTCAAAAATTATAGTTCAATTCAGAGAGGGATGGGATAAAATGAACTTTGTAACGTTGACTTCGGATGAGTTCAATGCGTTTACAACAAAGCATTTTTCACATTACACACAATCAGCTATTCATTACAATCATAGAGTTGATTTAAAAGGCGATGTGCATCTTGTAGGGGTTAAAGATGACAATGGTCAAGTGATTGCAGGATGCTTATTGACAGAA  >  USA300TCH1516_ALE/2532174‑2532430

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