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,2980AT100.0% 46.1 / NA 19T15T (ACA→ACTfemA_3Aminoacyltransferase FemA
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base T (6/13);  total (6/13)
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

TTTGTAATCAGTTTTTGATAAAATAAAAAAATAGCATTATTAATAAATTCAAAAATTATAGTTCAATTCAGAGAGGGATGGGATAAAATGAACTTTGTAACGTTGACTTCGGATGAGTTCAATGCGTTTACAACAAAGCATTTTTCACATTACACACAATCAGCTATTCATTACAATCATAGAGTTGATTTAAAAGGCGATGTGCATCTTGTAGGGGTTAAAGATGACAATGGTCAAGTGATTGCAGGATGCTTATTGACA  >  USA300TCH1516_ALE/2532167‑2532427
                                                                                                                                   |                                                                                                                                 
tttGTAATCAGTTTTTGATAAAATAAAAAAATAGCATTATTAATAAATTCAAAAATTATAGTTCAATTCAGAGAGGGATGGGATAAAATGAACTTTGTAACGTTGACTTCGGATGAGTTCAATGCGTTTActacaa                                                                                                                               <  2:121108/136‑1 (MQ=255)
                   aaaataaaaAAATAGCATTATTAATAAATTCAAAAATTATAGTTCAATTCAGAGAGGGATGGGATAAAATGAACTTTGTAACGTTGACTTCGTATGAGTTCAATTTTTTTACTATAAAGCATTTTTCACAATacac                                                                                                            <  1:58116/136‑1 (MQ=255)
                            aaaTAGCATTATTAATAAATTCAAAAATTATAGTTCAATTCAGAGAGGGATGGGATAAAATGAACTTTGTAACGTTGACTTCGGATGAGTTCAATGCGTTTACTACAAAGCATTTTTCACATTACACACAATCAGc                                                                                                   <  2:61080/136‑1 (MQ=255)
                                      attaATAAATTTAAGAATTATAGTTCAATTCAGAGAGGGATGGGATAAAATGAACTTTGTAACGTTGACTTCGGATGAGTTCAATGCGTTTACTACAAAGCATTTTTCACATTACACACAATCAGCTATTCATTAc                                                                                         <  2:129981/136‑1 (MQ=255)
                                                              tCAATTCAGAGAGGGATGGGATAAAATGAACTTTGTAACGTTGACTTCGGATGAGTTCAATGCGTTTACTACAAAGCATTTTTCACATTACACACAATCAGCTATTCATTACAATCATAGAGTTGATTTAAAAGGc                                                                 <  2:21052/136‑1 (MQ=255)
                                                                        gagGGATAGGATAAAATACACTTTGTAACGTTGACTTCAGATGAGTTCAATGCGTTTACTACAAAGAATTTTTCACATTACACACAATCATCTATTCATTACAATCATATAGTTTATTTAAAAGGCGATGTGCATc                                                       <  2:44596/136‑1 (MQ=255)
                                                                            gatgggatAAAATGAACTTTGTAACGTTGACTTCGGATGAGTTCAATGCGTTTACTACAAAGCATTTTTCACATTACACACAATCAGTTATTCATTATTATTTTAGTGTTTATTTAAAAGGCGATGTGCATCTTGt                                                   <  1:106779/136‑1 (MQ=255)
                                                                                               tGTAACGTTGACTTCGGATGAGTTCAATGCGTTTACTACAAAGCATTTTTCACATTACACACAATCAGCTATTCATTACAATCATAGAGTTGATTTAAAAGGCGATGTGCATCTTGTAGGGGTTAAAGATGACAAt                                >  1:25942/1‑136 (MQ=255)
                                                                                                gTAACGTTGACTTCGGATGAATTAAATACAAAAACTACAAAACATAATTCACATTACACACAATCAGCTATTCATTACAATCATAGAGTTGATTTAAAAAGCGATGTGCATCTTGTAGGGGTTAAAGATGACAATg                               >  1:6423/1‑136 (MQ=255)
                                                                                                    cGTTGACTTCGGATGAGTTCAATGCGTTTACTACAAAGCATTTTTCACATTACACACAATCAGCTATTCATTACAATCATAGAGTTGATTTAAAAGGCGATGTGCATCTTGTAGGGGTTAAAGATGACAATGGTCa                           >  2:75663/1‑136 (MQ=255)
                                                                                                      tttaCTTCGGATGAGTTCAATGCGTTTACTACAAAGCATTTTTCACATTACACACAATCAGCTATTCATTACAATCATAGAGTTGATTTAAAAGGCGTTGTGCATTTTTTTGGGGTTAAAGATGACAATGGTCAAg                         <  1:143090/133‑1 (MQ=255)
                                                                                                       tGACTTCGGATGAGTTCAATGCGTTTACTACAAAGCATTTTTCACATTACACACAATCAGCTATTCATTACAATCATAGAGTTGATTTAAAAGGCGATGTGCATCTTGTAGGGGTTAAAGATGACAATGGTCAAGt                        <  2:58567/136‑1 (MQ=255)
                                                                                                        gACTTCGGATGAGTTCAAGGTGTTTACTACAAAGCATTTTTCACATTACACACAATCAGCTATTCATTACAATCATAGAGTTGATTTAAAAGGCGATGTGCATCTTGTAGGGGTTAAAGATGACAATGGTCAAGTg                       <  2:156990/136‑1 (MQ=255)
                                                                                                          cTTCGGATGAGTTCAATGCGTTTACTACAAAGCATTTTTCACATTACACACAATCAGCTATTCATTACAATCATAGAGTTGATTTAAAATGCGATGTGTATCTTGTATGGGTTAAAGATGACAATGGTCAAGTGAt                     <  1:51429/136‑1 (MQ=255)
                                                                                                           ttCGGATGAGTTCAATGCGTTTACTACAAAGCATTTTTCACATTACACACAATCAGCTATTCATTACAATCATAGAGTTGATTTAAAAGGCGATGTGTATCTTGTTGTGTTTAAAGATGACAATGGTCAAGTGAtt                    <  1:95087/136‑1 (MQ=255)
                                                                                                           ttCGGATGAGTTCAATGCGTTTACTACAAAGCATTTTTCACATTACACACAATCAGCTATTCATTACAATCATAGAGTTGATTTAAAAGGCGATGTGTATCTTGTAGGGTTTAAAGATGACAATGGTCAAGTGAtt                    <  1:157932/136‑1 (MQ=255)
                                                                                                           ttCGGATGAGTTCAATGCGTTTACTACAAAGCATATTTCACATTACACACAATCAGCTATTCATTACAATCATAGAGTTGATTTAAAAGGCGATGTGCATCTTGTAGGGGTTAAAGATGACAATGGTCAAGTGAtt                    >  1:149905/1‑136 (MQ=255)
                                                                                                                     ttCAATGCGTTTACTACAAAGCATTTTTCACATTACACACAATCAGCTATTCACTACAATCATAGAGTTGATTTAAAAGGCGATGTGCATCTTGTAAGGGTTAAAGATGACAATGGTCAAGTGATTGCAGGACGCt          >  2:122997/1‑136 (MQ=255)
                                                                                                                        actGCGTTTACTACAAAGCATTTTTCACATTACACACAATCAGCTATTCATTACAATCATAGAGTTGATTTAAAAGGCGATGTGCATCTAGTAGGGGTTAAAGATGACAATGGTCAAGTGATTGCAGGATGCTTAt       >  1:7106/3‑136 (MQ=255)
                                                                                                                             gTTTATTACAAAGCATTTTTCAGATTACACACAATCAGCTATTCATTACAATCATAGAGTTGATTTAAAAGGCGATGTGCATCTTGTAGGGGTTAAAGATGACAATGGTCAAGTGATTGCAGGATGCTTATTGACa  <  2:141467/136‑1 (MQ=255)
                                                                                                                                   |                                                                                                                                 
TTTGTAATCAGTTTTTGATAAAATAAAAAAATAGCATTATTAATAAATTCAAAAATTATAGTTCAATTCAGAGAGGGATGGGATAAAATGAACTTTGTAACGTTGACTTCGGATGAGTTCAATGCGTTTACAACAAAGCATTTTTCACATTACACACAATCAGCTATTCATTACAATCATAGAGTTGATTTAAAAGGCGATGTGCATCTTGTAGGGGTTAAAGATGACAATGGTCAAGTGATTGCAGGATGCTTATTGACA  >  USA300TCH1516_ALE/2532167‑2532427

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: