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.7% 122.3 / ‑4.8 41T15T (ACA→ACTfemA_3Aminoacyltransferase FemA
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base T (13/27);  minor base C (0/1);  total (13/28)
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.42e-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.

GTATTTGTAATCAGTTTTTGATAAAATAAAAAAATAGCATTATTAATAAATTCAAAAATTATAGTTCAATTCAGAGAGGGATGGGATAAAATGAACTTTGTAACGTTGACTTCGGATGAGTTCAATGCGTTTACAACAAAGCATTTTTCACATTACACACAATCAGCTATTCATTACAATCATAGAGTTGATTTAAAAGGCGATGTGCATCTTGTAGGGGTTAAAGATGACAATGGTCAAGTGATTGCAGGATGCTTATTGACAGAAGCAC  >  USA300TCH1516_ALE/2532164‑2532434
                                                                                                                                      |                                                                                                                                        
gtATTTGTAATCAGTTTTTGATAAAATAAAAAAATAGCATTATTAATAAATTCAAAAATTATAGTTCAATTCAGAGAGGGATGGGATAAAATGAACTTTGTAACGTTGACTTCGGATGAGTTCAATGCGTTTACTACaaa                                                                                                                                     <  2:392048/140‑1 (MQ=255)
 tATTTGTAATCAGTTTTTGATAAAATAAAAAAATAGCATTATTAATAAATTCAAAAATTATAGTTCAATTCAGAGAGGGATGGGATAAAATGAACTTTGTAACGTTGACTTCGGATGAGTTCAATGCGTTTACTACAAAg                                                                                                                                    <  1:391231/140‑1 (MQ=255)
  agtggTAATCAGGTGTGGAAAAAATAAAAAAATAGCATTATTAATAAATTCAAAAATGATAGTTCAATTCAGAGAGGGATGGGATAAAAGGAACTTTGTAACGTTGACTTCGGATGAGTTCAATGCGTTTACTACAAAGc                                                                                                                                   <  1:501237/136‑1 (MQ=255)
  aTTTGTAATCAGTTTTTGATAAAATAAAAAAATATCATTATTAATAAATTCAAAAATTATAGTTCAATTCAGAGAGGGATGGGATAAAATGAACTTTGTAACGTTGACTTCGGATGAGTTGAATGCGTTTACTACAAAGc                                                                                                                                   <  2:155642/140‑1 (MQ=255)
  aTTTGTAATCAGTTTTTGATAAAATAAAAAAATAGCATTATTAATAAATTCAAAAATTATAGTTCAATTCAGAGAGGGATGGGATAAAATGAACTTTGTAACGTTGACTTCGGATGAGTTCAATGCGTTTACTACAAAGc                                                                                                                                   <  2:519196/140‑1 (MQ=255)
  aTTTGTAATCAGTTTTTGATAAAATAAAAAAATAGCATTATTAATAAATTCAAAAATTATAGTTCAATTCAGAGAGGGATGGGATAAAATGAACTTTGTAACGTTGACTTCGGATGAGTTCAATGCGTTTACTACAAAGc                                                                                                                                   <  2:467883/140‑1 (MQ=255)
   tttGTAATCAGTTTTTGATAAAATAAAAAAATAGCATTATTAATAAATTCAAAAATTATAGTTCAATTCAGAGAGGGATGGGATAAAATGAACTTTGTAACGTTGACTTCGGATGAGTTCAATGCGTTTACTACAAAGCa                                                                                                                                  <  2:170065/140‑1 (MQ=255)
     tGTAATTAGTTTTTGATAAAATAAAAAAATAGCATTATTAATTAATCCAAAAATGATAGTTCAATTCAGAGAGGGATGGGATAAAATGAACTGTGTACCGTTGACTTCGGATGAGTTCAATGCGTTTACTACAAAGCAtt                                                                                                                                <  2:12947/140‑1 (MQ=255)
      gTAATCAGTTTTTGATAAAATAAAAAAATAGCATTATTAATAAATTCAAAAATTATAGTTCAATTCAGAGATGGATGGGATAAAATGAACTTTGTAACGTTGACTTCGGATGAGTTCAATGCGTTTACTACAAAGCAttt                                                                                                                               <  1:37372/140‑1 (MQ=255)
      gTAATCAGTTTTTGATAAAATAAAAAAATAGCATTATTAATAAATTCAAAAATTATAGTTCAATTCAGAGAGGGATGGGATAAAATGAACTTTGTAACGTTGACTTCGGATGAGTTCAATGCGTTTACTACAAAGCAttt                                                                                                                               <  1:396033/140‑1 (MQ=255)
         aTCAGTTTTTGATAAAATAAAAAAATAGCATTATTAATAAATTCAAAAATTATAGTTCAATTCAGAGAGGGATGGGATAAAATGAACTTTGTAACGTTGACTCCGGATGAGTTCAATGCGTTTACTACAAAGCATTTTTc                                                                                                                            >  2:238511/1‑140 (MQ=255)
              tttttGATAAAATAAAAAAATAGCCTTATTAATAACTTCAAAACTTATTGTTCAATTCAGAGAGGGATGGGATAAAATGAACTTTGTAACGTTGACTTCGGATGAGTTCAATGCGTTTACTACAAAGCATTTTTCATAtt                                                                                                                       <  2:2166/140‑1 (MQ=255)
                    aaaaaataaaaAAATACGATTATTAATAAATTCAAAAATTATAGTTAAATCCAGAGAGGTATGTGAAAAAATGAACTTTGCAACGTTGACTTCGGATGAGTTCAATGCGCTGACCACAACGCACTTTTCACATTacacac                                                                                                                 <  2:142458/138‑1 (MQ=255)
                      aaaataaaaAAATAGCATTATTAATAAATTCAAAAATTATAGTTCAATTCAGAGAGGGATGGGATAAAATGAACTTTGTAACGTTGACTTCGGATGAGTTCAATGCGTTTACTACAAAGCATTTTTCACATTACACACaa                                                                                                               <  1:270677/140‑1 (MQ=255)
                            aaaaaaTAGCATTATTAATAAATTCAAAAATTATAGTTCAATTCAGAGAGGGATGGGATAAAATGAACTTTGTAACGTTGACTTCGGATGAGTTCAATGCGTTTACTACAAAGCATTTTTCACATTACACACAATCAGCt                                                                                                         <  1:428637/140‑1 (MQ=255)
                             aaaaaTAGCATTATTAATAAATTCAAAAATTATAGTTCAATTCAGAGAGGGATGGGATAAAATGAACTTTGTAACGTTGACTTCGGATGAGTTCAATGCGTTTACTACAAAGCATTTTTCACATTACACACAAACAGCTa                                                                                                        >  1:94981/1‑140 (MQ=255)
                                        tattaATAAATTCAAAAATTATAGTTCAATTCAGAGAGGGATGGGATAAAATGAACTTTTTAACGTTGACTTCGGATGAGTTCAATGCGTTTTCTACAAAGCATTTTTCACATTACACACAATCAGCTATTCATTACAAt                                                                                             <  1:171268/140‑1 (MQ=255)
                                           taataaGTTCAAAAATTATAGTTCAATTCAGGGAGGGATGGGATAAAATGAACTTTGTAACGTTTACTTCGGATGAGTTCAATGCGTTTACTACAAAGCATTGTTCACATTACACACAATCAGCTATTCATTACAATCAt                                                                                          <  2:247978/140‑1 (MQ=255)
                                           taataaATTCAAAAATTATAGTTCAATTCAGAGAGGGATGGGATAAAATGAACTTTGTAACGTTGACTTCGGATGAGTTCAATGCGTTTACTACAAAGCATTTTTCACATTACACACAATCAGCTATTCATTACAATCAt                                                                                          <  2:304590/140‑1 (MQ=255)
                                                     aaaaaTTATAGTTCAATTCAGAGAGGGATGGGATAAAATGAACTTTGTAACGTTGACTTCGGATGAGTTCAATGCGTTTACTACAAAGCATTTTTCACATTACACACAATCAGCTATTCATTACAATCATAGAGTTGAtt                                                                                <  2:256746/140‑1 (MQ=255)
                                                         aTTATAGTTCAATTCAAAGAGGGATGGGATAAAATGAACTTTGTAACGTTGACTTCGGATGAGTTCAATGCGTTTACTACAAAGCATTTTTCTCGTCACGCACAATCAGCTATTCATTACAATCATAGAGTTGATTTaaa                                                                            <  1:378425/140‑1 (MQ=255)
                                                              aGTTCAATTCAGAGAGGGATGGGATAAAATGAACTTTGTAACGTTGACTTCGGATTAGTTCAATGCGTTTACTACAAAGCATTTTTCACATTACACACAATCAGCTATTCATTACAATCATAGAGTTGATTTAAAAGGCg                                                                       <  1:281220/140‑1 (MQ=255)
                                                              aGTTCAATTCAGAGAGGGATGGGATAAAATGAACTTTGTAACGTTGACTTCGGATGAGTTCAATGCGTTTACTACAAAGCATTTTTCACATTACACACAATCAGCTATTCATTACAATCATAGAGTTGAATTAAAAAGCg                                                                       >  1:156490/1‑140 (MQ=255)
                                                              aGTTCAATTCAGAGAGGGATGGGATAAAATGAACTTTGTAACGTTGACTTCGGATGAGTTAAATGCGTTTACTACAAAGCATTTTTTAAATTACACAGAATCAGCTATTCATTACAATCATAGAGTTGATTTaaaaagag                                                                       >  1:509639/1‑136 (MQ=255)
                                                                 tCAATTCAGAGAGGGATGGGATAAAATGAACTTTGTAACGTTGACTTCGGATGAGTTCAATGCGTTTACTACAAAGCATTTTTCACATTACACACAATCAGTTATTCATTACAATCATAGATTTGATTTAAAAGGCGAtg                                                                    <  2:254506/140‑1 (MQ=255)
                                                                              ggatgggatAAAATGAACTTTGTAACGTTGACTTCGGATGAGTTCAATGCGTTTACTACAAAGCATTTTTCACATTACACACAATCAGCTATTCATTCCAATCATAGAATAAATTTAAAATGCGATGTGCATCtgttaga                                                       >  1:258715/1‑134 (MQ=255)
                                                                                     atAAAATGAACTTTGTAACGTTGACTTCGGATGAGTTCAATGCGTTTACTACAAAGCATTTTTCACATTACACACAATCAGCTATTCATTACAATCATAGAGTTGATTTAAAAGGCGATGTGCATCTTGTAGGGGTTcaa                                                >  2:217295/1‑140 (MQ=255)
                                                                                     atAAAATGAACTTTGTAACGTTGACTTCGGATGAGTTCAATGCGTTTACTACAAAGCAGTTTTCACATTACACACAATCAGCTATTCATTACAATCATAGAGTTGATTTAAAAAGCCATGTGCATCTTGTAGGGGTTaaa                                                >  2:275151/1‑140 (MQ=255)
                                                                                            gAACTTTGTAACGTTGACTTCGGATGAGTTCAATGCGTTTACTACAAAGCATTTTTCACATTACACACAATCAGCTATTCATTACAATCATAGAGTTGATTTAAAAGGCGATGTGCATCTTGTAGGGGTTAAAGATGaaa                                         <  1:545921/140‑3 (MQ=255)
                                                                                             aaCTTTGTAACGTTGACTTCGGATGAGTTCAATGCGTTTACTACAAAGCATTTTTCACATTACACACAATCAGCTATTCCTTACAATCATAGAGTTGATTTAAAAGGCGATGTGCATCTTGTAGGGGTTAAAGATGACaa                                        <  2:13174/140‑1 (MQ=255)
                                                                                             aaCTTTGTAACGTTGACTTCGGATGAGTTCAATGCGTTTACTACAAAGCATTTTTCACATTACAAACAATCAGCTATTCATTACAATCATAGAGTTGATTTAAAAGGCGATGTGCATCTTGTAGGGGTTAAAGATGACaa                                        >  2:381818/1‑140 (MQ=255)
                                                                                                tttGTAACGTTGACTTCGGGTGAGTTCAATGCGTTTACTACAAAGCATTTTTCACATTACACACAATCAGATATTCATTACAATCATAGAGTTGATTTAAAAGGCGATGTGCATCTTGNANGGGTTGAAGATGACAAggg                                     >  1:282124/1‑140 (MQ=255)
                                                                                                tttGTAACGTTGACTTCGGATGAGTTCAATGCGTTTACTACAAAGCATTTTTCACATTACACACAATCAGCTATTCATTACAATCATAGAGTTGATTTAAAAGGCGATGTGCATCTTGTAGGGGTTAAAGATGACAGTgg                                     >  2:296622/1‑140 (MQ=255)
                                                                                                tttGTAACGTTGACTTCGGATGAGTTCAATGCGTTTACTACAAAGCATTTTTCACATTACACACAATCAGCTATTCATTACAATCATAGAGTTGATTTAAAAGGCGATGTGCATCTTGTAGGGGTTAAAGATGACAATgg                                     <  2:73409/140‑1 (MQ=255)
                                                                                                  tGTAACGTTGACTTCGGATGAGTTCAATGCGTTTACTACAAAGCATTTTTCACATTACACACAATCAGCTATTCATTACAATCATAGAGTTTATTTAAAAGGCGATGTGCATCTTGTAGGGGTTAAAGATGACAATGGTc                                   <  1:240095/140‑1 (MQ=255)
                                                                                                   gTAATGTTGACTTCGGATGAGTTCAATGCGTTTACTATAAAGCATTTTTCACATTACACACAATCAGCTATTCATTACAATCATAGAGTTGATTTAAAAGGCGATGTGCATCTTGTAGGGGTTAAAGATGACAATGGTCa                                  <  1:56217/140‑1 (MQ=255)
                                                                                                            aCTTCGGATGAGTTCAATGCGGTTACTAAAAAGCATTTTTCACATTACACACAATCAGCTATTCATTACAATCATAGAGTTGATTTAAAAGGCGATGTGCATCTTGTAGGGGTTAAAGATGACAATGGTCAAGTGATTGc                         <  2:417169/140‑1 (MQ=255)
                                                                                                                 ggATGAGTTCAATGCGTTTACTACAAAGCATTTTTCACATTACACACAATCAGCTATTCATTCCAATCATAGAGTTGATTTAAAAGGCGATGTGCATCTTGTAGGGGTTAAAGATGAAAATGGTCAAGTGATTGCAGGAt                    >  2:274906/1‑140 (MQ=255)
                                                                                                                 ggATGAGTTCAATGCGTTTACTACAAAGCATTTTTCACATTAAACACAATCAGCTATTCATTACAATCATAGAGTTGATTTAAAAGGCGATGTGCATCTTGTAGGGGTTAAAGATGACAATGGTCAAGTGATTGCAGGAt                    >  2:275013/1‑140 (MQ=255)
                                                                                                                      aGTTCAATGCGTTTACTACAAAGCATTTTTCACATTACACACAATCAGCTATTCATTACAATCATAGAGTTGATTTAAAAGGCGATGTGCATCTTGTAGGGGTTAAAGATGACAATGGTCAAGTGATTGCAGGATGATTa               >  2:320399/1‑140 (MQ=255)
                                                                                                                                  ttACTACAAAGCATTTTTCACATTACACACAATCAGCTATTCATTACAATCATAGAGTTGATTTAAAAGGCGATGTGCATCTTGTAGGGGTTAAAGATGACAATGGTCAAGTGATTGCAGGATGCTTATTGACAGAAgca   <  1:445530/140‑1 (MQ=255)
                                                                                                                                   tactacAAAGCATTTTTCACATTACACACAATCAGCTATTCATTACAATCATAGAGTTGATTTAAAAGGCGATGTGCATCTTGTAGGGGTTAAAGATGACAATGGTCAAGTGATTGCAGGATGCTTATTGACAGAAgcac  >  2:97141/1‑140 (MQ=255)
                                                                                                                                      |                                                                                                                                        
GTATTTGTAATCAGTTTTTGATAAAATAAAAAAATAGCATTATTAATAAATTCAAAAATTATAGTTCAATTCAGAGAGGGATGGGATAAAATGAACTTTGTAACGTTGACTTCGGATGAGTTCAATGCGTTTACAACAAAGCATTTTTCACATTACACACAATCAGCTATTCATTACAATCATAGAGTTGATTTAAAAGGCGATGTGCATCTTGTAGGGGTTAAAGATGACAATGGTCAAGTGATTGCAGGATGCTTATTGACAGAAGCAC  >  USA300TCH1516_ALE/2532164‑2532434

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: