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,2980AT96.1% 64.2 / ‑5.3 23T15T (ACA→ACTfemA_3Aminoacyltransferase FemA
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base T (4/18);  minor base C (0/1);  total (4/19)
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.39e-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.

GTGTATTTGTAATCAGTTTTTGATAAAATAAAAAAATAGCATTATTAATAAATTCAAAAATTATAGTTCAATTCAGAGAGGGATGGGATAAAATGAACTTTGTAACGTTGACTTCGGATGAGTTCAATGCGTTTACAACAAAGCATTTTTCACATTACACACAATCAGCTATTCATTACAATCATAGAGTTGATTTAAAAGGCGATGTGCATCTTGTAGGGGTTAAAGATGACAATGGTCAAGTGATTGCAGGATGCTTATTGACAGAAGCA  >  USA300TCH1516_ALE/2532162‑2532433
                                                                                                                                        |                                                                                                                                       
gtgtATTTGTAATCAGTTTTTGATAAAATAAAAAAATAGCATTATTAATAAATTCAAAAATTATAGTTCAATTCAGAGAGGGATGGGATAAAATGAACTTTGTAACGTTGACTTCGGATGAGTTCAATGCGTTTActacaa                                                                                                                                     >  1:88738/1‑141 (MQ=255)
          aaTCAGTTTTTGATAAAATAAAAAAATAGCATTATTAATAAATTCAAAAATTATAGTTCAATTCAGAGAGGGATGGGATAAAATGAACTTTGTAACGTTGACTTCGGATGAGTTCAATGCGTTTACTACAAAGCATTTTTc                                                                                                                           <  1:323798/141‑1 (MQ=255)
              agattttGATAACATAAAAAAATCGCATTATTACTAAATTCAAAAATAATATATACAAACCGAGAGGGATGGGATAAAAAGAACTAAGTAACGTTGACTTCGGATGGGTTCAATGCGTTTACTACAAAGCATTTTTCAAAt                                                                                                                       >  2:87276/4‑141 (MQ=255)
                         aaataaaaAAATAGCATTATTAATAAATTCAAAAATTATAGTTCAATTCAGAGAGGGATGGGATAAAATGAACTTTGTAACGTTGACTTCGGATGAGTTCAATGCGTTTACTACAAAGCATTTTTCACATTACACACAATc                                                                                                            <  2:172140/141‑1 (MQ=255)
                            taaaaAAATAGCATTATTAATAAATTCAAAAATTATAGTTCAATTCAGAGAGGGATGGGATAAAATGAACTTTGTAACGTTGACTTCGGATGAGTTCAATGCGTTTACTACAAAGCATTTTTCACATTACACACAATCAGc                                                                                                         <  1:67189/141‑1 (MQ=255)
                            taaaaAAATAGCATTATTAATAAATTCAAAAATTATAGTTCAATTCAGAGAGGGATGGGATAAAATGAACTTTGTAACGTTGACTTCGGATGAGTTCAATGCGTTTACTACAAAGCATTTTTCACATTACACACAATCAGc                                                                                                         <  2:282099/141‑1 (MQ=255)
                                 aaaTAGCATTATTAATAAATTCAAAAATTATAGTTCAAGTCAGATAGGGATGGGATAAAATGAACTTTGTAACGTTGACTTCGGATGAGTTCAATGCGTTTACTACAAAGCATTTTTCACATTACATACAATCTGCTATTc                                                                                                    <  2:289679/141‑1 (MQ=255)
                                  aaTAGCATTATTAATAAATTCAAAAATTATAGTTCAATTCAGAGAGGGATGGGATAAAATGAACTTTGTAACGTTGACTTCGGATGATTTCAATGCGTTTACTACAAAGCATTTTTCACATTACACACAATCAGCTATTCa                                                                                                   <  1:324692/141‑1 (MQ=255)
                                    tAGCATTATTAATAAATTCAAAAATTATAGTTCAATTCAGAGAGGGATGGGATAAAATGAACTTTGTAACGTTGACTTCGGATGAGTTCAATGCGTTTACTACAAAGCATTTTTCACATTACACACAATCAGCTATTCAtt                                                                                                 <  2:319572/141‑1 (MQ=255)
                                                 aaaTTCAAAAATTATAGTTCAATTCAGAGAGGGATGGGATAAAATGAACTTTGTAACGTTGACTTCGGATGAGTTCAATGCGTTTACTACAAAGCATTTTTCACATTACACACAATCAGCTATTCATTACAATCATAGAGt                                                                                    <  1:176870/141‑1 (MQ=255)
                                                  aactCAAAAAGTATAGTTTAATTTAGAGAGGGATGGGATAAAATGAACTTTGTAACGTTGACTTCGGATGAGTTCAATGCGTTTACTACAAAGTATTTTTCACATTACACACAATCAGCTATTCATTACAATCATAGAGtt                                                                                   <  2:223229/138‑1 (MQ=255)
                                                  aaTTCAATAATTATAGTCAAATTCAGAGAGGGATGGGATAAAATGAACTTTGTAACGTTGACTTCGGATGAGTTCAATGCGTTTACTACAAAGCATTTTTCACATTACACACAATCAGCTATTCATTACAATCATAGAGtt                                                                                   <  2:58397/141‑1 (MQ=255)
                                                                  ttCAATTCAGAGAGGGATGGGATAAAATGAACTTTGTAACGTTGACTTCGGATGAGTTCAATGCGTTTACTACAAAGCATTTTTCACATTACACACAATCAGCTATTCATTACAATCATAGAGTTGATTTAAAAGGCGAtg                                                                   <  2:175659/141‑1 (MQ=255)
                                                                           gagagGGATGGGATAAAATGAACTTTGTAACGTTGACTTCGGATGAGTTCAATGCGTTTACTACAAAGCATTTTTCACATTACACACAATCAGCTATTCATTACAATCATAGAGTTGATTTAAAAGGCGATGTGCATCTTg                                                          <  2:101875/141‑1 (MQ=255)
                                                                               gtgatgggatAAAAATAACTTTGTAACGTTGACTTCGGATGAGTTCAATGCGTTTACTACAAAGCATTTTTCACATTACACACAATCAGCTATTCATTACAATCATAGAGTTGATTTAAAAGGCGATGTGCATCTTGTAgg                                                      <  1:231270/139‑1 (MQ=255)
                                                                                   agggctACAATGAACTTTGTAACGTGGACTTCGGAATAGTTCAATGCGTTTACTACAAAGCAATTGTCACATTACACACAATCAGCTATTCATTACAATCATAGAGTTGATTTAAAACGCGATGTGCATCTTGTAGGGGtt                                                  <  2:88709/140‑1 (MQ=255)
                                                                                     ggatAAAATGAACTTTGTAACGTTGACTTCGGATGAGTTCAATGCGTTTACTACAAAGCATTTTTCACATTACACACAATCAGCTATTCATTACAATCATAGAGTTGATTTAAAAGGCGGTGTGCATCTTGTAGGGGttga                                                >  1:265552/1‑139 (MQ=255)
                                                                                               aaCTTTGTAACGTTGACTTCGGATGAGTTCAATGCGTTTACTACAAAGCATTTTTCACATTACACACAATCAGCTATTCATTACAATCATAGAGTTGCTTTAAAAGGCGCTGTGCATCTTGTAGGGGTTAAATATGACAAt                                      >  2:102048/1‑141 (MQ=255)
                                                                                                 cTTTGTAACGGTGACGTCGGATGAGTTCAATGTGTTTACTACAAGGCTTTTTTCACATTACACACAATCAGCTATTCATTACAATCATAGAGTTGATTTAAAAGGCGATGTGCATCTTGTAGGGGTTAAAGATGACAATgg                                    <  1:87276/141‑1 (MQ=255)
                                                                                                              aCTTCGGATGAGTTCAATGCGTTTACTACAAAGCATTTTTCACATTACACACAATCAGCTATTCATTACAATCATAGAGTTGATTTAAAAGGCGATGTGCATCTTGTAGGGGTTAAAGATGACAATGGTCAAGTGATTGCa                       >  1:276139/1‑141 (MQ=255)
                                                                                                                  cGGATGAGTTCAATGCGTCTACTACAAAGCATTTTTCACATTACACACAATCAGCTATTCATTACAATCATAGAGTTGATTTAAAAGGCGATGTGCATCTTGTAGGGGTTAAAGATGACAATGGTCAAGTGATTGCAGGAt                   <  1:90748/141‑1 (MQ=255)
                                                                                                                                  gTTTACCACAAAGCATTTTTCACATTACACACAATCAGATGTTCATTAAGAACTGAGACTTGATTTAAAAGGCGATGTGCATCTTGTAGGGGTTAAAGATGACAATTGTCAAGTGATTGCAGGATGCTTATTGACAGAAgc   <  1:9789/141‑1 (MQ=255)
                                                                                                                                   tttACTACAAAGCTTTTTTCACATTACACACAATCAGCTATTCATTACAATCATAGAGTTGATTTAAAAGGCGATGGTCATCTTGTAGGGGTTAAAGATGACAATGGTCAAGTGATTGCAGGATGCTTATTGACAGAAgca  <  2:102008/141‑1 (MQ=255)
                                                                                                                                   tttACTACAAAGCATTTTACACATTACACACAATCAACTATTCATTACAATCATAGAGTTGATTTAAAAGGCGATGTGCATCTTGTAGGGGTTAAAGATGACAATGGTCAAGTGATTGCAGGATGCTTATTGACAGAAgca  <  1:322145/141‑1 (MQ=255)
                                                                                                                                        |                                                                                                                                       
GTGTATTTGTAATCAGTTTTTGATAAAATAAAAAAATAGCATTATTAATAAATTCAAAAATTATAGTTCAATTCAGAGAGGGATGGGATAAAATGAACTTTGTAACGTTGACTTCGGATGAGTTCAATGCGTTTACAACAAAGCATTTTTCACATTACACACAATCAGCTATTCATTACAATCATAGAGTTGATTTAAAAGGCGATGTGCATCTTGTAGGGGTTAAAGATGACAATGGTCAAGTGATTGCAGGATGCTTATTGACAGAAGCA  >  USA300TCH1516_ALE/2532162‑2532433

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: