Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 665,392 T→G 100% L79* (TTA→TGA)  paiA_1 → Spermidine/spermine N(1)‑acetyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE665,3920TG100.0% 122.2 / NA 37L79* (TTA→TGA) paiA_1Spermidine/spermine N(1)‑acetyltransferase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base G (22/15);  total (22/15)
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

CACTTGTTAGAAGCTTCATATAATGAACATCATTTAAAGAAAAGACTTCAAGAACAATTATTCTTAGTCGTTGAAGAAGGTAATGACATCGTTGGCTTTGCTAACTTTATTTACGGTGAAGAATTATATTTATCAGCTCATTATGTTAAACCAGAATCGCAACATACAGGTTATGGTACAGCATTGTTAAATGAAGGATTATCACGTTTTGAAGATAAATTTGAAGGTGTTTACTTAGAAGTAGATAATAAAAATGAAGA  >  USA300TCH1516_ALE/665268‑665527
                                                                                                                            |                                                                                                                                       
cACTTGTTAGAAGCTTCATATAATGAACATCATTTAAAGAAAAGACTTCAAGAACAATTATTCTTAGTCGTTGAAGAAGGTAATGACATCGTTGGCTTTGCTAGCTTTATTTACGGTGAAGAATGATATTTATCAGCTCAt                                                                                                                         <  1:444181/141‑1 (MQ=255)
      ttAGAAGCTTCATATAATGAACATCATTTAAAGAAAAGACTTCAAGAACAATTATTCTTAGTCGTTTAAGAAGGTAATGACATCGTTGGCTTTGCTAACTTTATTTACGGTGAAGAATGATATTTATCAGCTCATGATGtt                                                                                                                   >  1:601745/1‑141 (MQ=255)
      ttAGAAGCTTCATATAATAAACATCGTTTAAAGAAAAGACTTCAAGAACAATTATTCTTAGTCGTTTAAGAAGGTAATGACATCGTTGGCTTTGCTAAATTTATTTACGGTGAAGAATGATATTTATCAGCTCATTATAtt                                                                                                                   >  1:370186/1‑141 (MQ=255)
       tAGAAGCTTCATATAATGAACATCATTTAAAGAAAAGACTTCAAGAACAATTATTCTTAGTCGTTGAAGAAGGTAATGACATCGTTGGCTTTGCTAACTTTATTTACGGTGAAGAATGATATTTATCAGCTCATTATGTTa                                                                                                                  <  1:5885/141‑1 (MQ=255)
        aGAAGCTTCATATAATGAACATCATTTAAAGAAAAGACTTCAAGAACAATTATTCTTAGTCGTTGAAGAAGGTAATGACATCGTTGGCTTTGCTAACTTTATTTACGGTGAAGAATGATATTTATCAGCTCATTATGTTaa                                                                                                                 <  1:527555/141‑1 (MQ=255)
          aaGCTTCATATAATGAACATCATTTAAAGAAAAGACTTCAAGAACAATTATTCTTAGGCGTTGAAGAAGGTAATGACATCGTTGGCTTTGCTAACTTTATTTACGGTGAAGAATGATGTTTATCAGCTCATTATGTTAAAc                                                                                                               >  1:275043/1‑141 (MQ=255)
                tatatAATGAACATCACTTAAAGAAAAGACTTCAAGAACAATTATTAGTACTCGTTAAAGAAGGTAATGACATCGTTGGCTTTGCCAACTTTATTTACGGTGACTAATGATATTTATCAGCTCATTATGTTAAACCAGAAt                                                                                                         <  1:303262/140‑1 (MQ=255)
                  tataATGAACATCATTTAAAGAAAAGACTTCAAGAACAATTATTCTTAGTCGTTGAAGAAGGTAATGACATCGTTGGCTTTGCTAACTTTATTTACGGTGGAGAATGATATTTATTCGCTCATTATGTTAAACCAGGATCg                                                                                                       >  1:28832/1‑141 (MQ=255)
                  tataATGAACATCATTTAAAGAAAAGACTTCAAGAACAATTATTCTTAGTCGTTGAAGAAGGTAATGACATCGTTGGCTTTGCTAACTTTATTTACGGTGAAGAATGATATTTATCAGCTCATTATGTTAAACCAGAATCg                                                                                                       >  1:684729/1‑141 (MQ=255)
                          aCATCATTTAAAGAAAAGACTTCAAGAACAATTATTCTTAGTCGTTGAAGAAGGTAATGACATCGTTGGCTTTGCTAACTTTATTTACGGTGAAGAATGATATTTATCAGCTCATTATGTTAAACCAGAATCGCAACATAc                                                                                               >  1:458168/1‑141 (MQ=255)
                             tcatTTAAAGAAAAGACTTCAAGAACAATTATTCTTAGTCGTTGAAGAAGGTAATGACATCGTTGGCTTTGCTAACTTTATTTACGGTGCTGCATGAGATTTATCAGCTCATTATGTTAAACCAGAATCGCAACATACAgg                                                                                            >  1:359713/1‑141 (MQ=255)
                             tcatTTAAAGAAAAGACTTCAAGAACAATTATTCTTAGTCGTTGAAGAAGGTAATGACATCGTTGGCTTTGCTAACTTTATTTACGGTGAAGAATGATATTTATCAGCTCATTATGTTAAACCAGAATCGCAACATACAgg                                                                                            >  1:157370/1‑141 (MQ=255)
                                tttAAAGAAAAGACTTCAAGAACAATTATTCTTAGTCGTTGAAGAAGGTAATGACATCGTTGGCTTTGCTAACTTTATTTACGGTGAAGAATGATATTTATCAGCTCATTATGTTAAACCAGAATCGCAACCCACAGGGTa                                                                                         >  1:748210/1‑141 (MQ=255)
                                 ttAAAGAAAAGACTTCAAGAACAATTATTCTTAGTCGTTGAAGAAGGTAATGACATCGTTGGCTTTGCTAACTTTATTTACGGTGAAGAATGATATTTATCAGCTCATTATGTTAAACCAGAATCGCAACATACAGGTTAt                                                                                        <  1:124158/141‑1 (MQ=255)
                                 ttAAAGAAAAGACTTCAAGAACAATTATTCTTAGTCGTTGAAGAAGGTAATGACATCGTTGGCTTTGCTAACTTTATTTACGGTGAAGAATGATATTTATCAGCTCATTATGTTAAACCAGAATCGCAACATACAGGTTAt                                                                                        <  1:29381/141‑1 (MQ=255)
                                          agaCTTCAAGAACAATTATTCTTAGTCGTTGAAGAAGGTAATGACATCGTTGGCTTTGCTAACTTTATTTACGGTGAAGAATGATATTTATCAGCTCATTATGTTAAACCAGAATCGCAACATACAGGTTATGGTACAGCa                                                                               >  1:574089/1‑141 (MQ=255)
                                            aCTTCAAGAACAATTATTCTTAGTCGTTGAAGAAGGTAATGACATCGTTGGCTTTGCTAACTTTATTTACGGTGAAGAATGATATTTATCAGCTCATTATGTTAAACCAGAATCGCAACATACAGGTTATGGTACAGCAtt                                                                             >  1:323057/1‑141 (MQ=255)
                                                  aGAACAATTATTCTTAGTCGTTGAAGAAGGTAATGACATCGTTGGCTTTGCTAACTTTATTTACGGTGAAGAATGATATTTATCAGCTCATTATGTTAAACCAGAATCGCAACATACAGGTTATGGTACAGCATTGTTaaa                                                                       >  1:663204/1‑141 (MQ=255)
                                                       aaTTATTCTTAGTCGTTGAAGAAGGTAATGACATCGTTGGCTTTGCTAACTTTATTTACGGTGAAGAATGATATTTATCAGCTCATTATGTTAAACCAGAATCGCAACATACAGGTTATGGTACAGCATTGTTAAATGAAg                                                                  >  1:148936/1‑141 (MQ=255)
                                                                tAGTCGTTGAAGAAGGTAATGACATCGTTGGCTTTGCTAACTTTATTTACGGTGAAGAATGATATTTATCAGCTCATTATGTTAAACCAGAATCGCAACATACAGGTTATGGTACAGCATTGTGAAATGAAGGATTATCAc                                                         >  1:616701/1‑141 (MQ=255)
                                                                      ttGAAGAAGGTAATGACATCGTTGGCTGTGCTAACTTTATTTACGGTGAAGAATGATATTTATCAGCTCATTATGTTAAACCAGAATCGCAACATACAGGTTATGGTACAGCATTGTTAAATGAAGGATTATCACGTTTTg                                                   <  1:695161/141‑1 (MQ=255)
                                                                            aagGTAATGACATCGTTGGCTTTGCTAACTTTATTTACGGTGAAGAATGATATTCATCAGCTCATTATGTTAAACCAGAATCGCAACATACAGGTTATGGTACAGCATTGTTAAATGAAGGATTATCACGTTTTGAAGATa                                             >  1:508045/1‑141 (MQ=255)
                                                                               gTAATGACATCGTTGGCTTTGCTAACTTTATTTACGGTGAAGAATGATATTTATCAGCTCATTATGTTAAACCAGAATCGCAACATACAGGGTATGGTACAGCATTGTTAAATGAAGGATTATCACGTTTTGAAGATAAAt                                          >  1:423630/1‑141 (MQ=255)
                                                                                  aTGACATCGTTGGCTTTGCTAACTTTATTTACGGTGAAGAATGATATTTATCAGCTCATTATGTTAAACCAGAATCGCAACATACAGGTTATGGTACAGCATTGTTAAATGAAGGATTATCACGTTTTGAAGATAAATTTg                                       <  1:460132/141‑1 (MQ=255)
                                                                                       aTCGTTGGCTTTGCTAACTTTATTTACGGTGAAGAATGATATTTATCAGCTCATTATGTTAAACCAGAATCGCAACATACAGGTTATGGTACAGCATTGTTAAATGAAGGATTATCACGTTGTGAAGATAAATTTGAAGgt                                  <  1:404556/141‑1 (MQ=255)
                                                                                               cTTTGCTAACTTTATTTACGGTGAAGAATGATATTTATCAGCTCATTATGTTAAACCAGAATCGCAACATACAGGTTATGGTACAGCATTGTTAAATGAAGGATTATCACGTTTTGAAGATAAATTTGAAGGTGTTTACtt                          <  1:265035/141‑1 (MQ=255)
                                                                                               cTTTGCTAACTTTATTTACGGTGAAGAATGATATTTATCAGCTCATTATGTTAAACCAGAATCGCAACATACAGGTTATGGTACAGCATTGTTAAATGAAGGATTATCACGTTTTGAAGATAAATTTGAAGGTGTTTACtt                          >  1:658025/1‑141 (MQ=255)
                                                                                                         tttatttaCGGTGAAGAATGATATTTATCAGCTCATTATGTTAAACCAGAATCGCAACATACAGGTTATGGTACAGCATTGTTAAATGAAGGATTATCACGTTTTGAAGATAAATTTGAAGGTGTTTACTTAGAAGTAGat                >  1:223540/1‑141 (MQ=255)
                                                                                                         tttatttaCGGTGAAGAATGATATTTATCAGCTCATTATGTTAAACCAGAATCGCAACATACAGGTTATGGTACAGCATTGTTAAATGAAGGATTATCACGTTTTGAAGATAAATTTGAAGGTGTTTACTTAGAAGTAGat                <  1:284807/141‑1 (MQ=255)
                                                                                                         tttatttaCGGTGAAGAATGATATTTATCAGCTCATTATGTTAAACCAGAATCGCAACATACAGGTTATGGTACAGCATTGTTAAATGAAGGATTATCACGTTTTGAAGATAAATTTGAAGGTGTTTACTTAGAAGTAGat                >  1:737694/1‑141 (MQ=255)
                                                                                                            atttaCGGTGAAGAATGATATTTATCAGCTCATTATGTTAAACCAGAATCGCAACATACAGGTTATGGTACAGCATTGTTAAATGAAGGATTATCACGTTTTGAAGATAAATTTGAAGGTGTTTACTTAGAAGTAGataat             <  1:265783/141‑1 (MQ=255)
                                                                                                               taCGGTGAAGAATGATATTAATCAGCTCATTATGTTAAACCAGAATCGCAACATACAGGTTATGGTACAGCATTGTTAAATGAAGGATTATCACGTTTTGAAGATAAATTTGAAGGTGTTTACTTAGAAGTAGATAATaaa          >  1:458062/1‑141 (MQ=255)
                                                                                                                 cGGTGAAGAATGATATTTATCAGCTCATTATGTTAAACCAGAATCGCAACATACAGGTTATGGTACAGCATTGTTAAATGAAGGATTATCACGTTTTGAAGATAAATTTGAAGGTGTTTACTTAGAAGTAGATAATaaaaa        >  1:510750/1‑141 (MQ=255)
                                                                                                                  ggTGAAGAATGATATTTATCAGCTCATTATGTTAAACCAGAATCGCAACATACAGGTTATGGTACAGCATTGTTAAATGAAGGATTATCACGTGTGGAAGATAAATTTGAAGGTGTTTACTTAGAAGTAGATAATAAAAAt       >  1:435790/1‑141 (MQ=255)
                                                                                                                     gaagaaTGATATTTATCAGCTCATTATGTTAAACCAGAATCGCAACATACAGGTTATGGTACAGCATTGTTAAATGAAGGATTATCACGTTTTGAAGATAAATTTGAAGGTGTTTACTTAGAAGTAGATAATAAAAATgaa    <  1:726858/141‑1 (MQ=255)
                                                                                                                      aagaaTGATATTTATCAGCTCATTATGTCAAACCAGAATCGCAACATACAGGTTATGGTACAGCATTGTTAAATGAAGGATTATCACGTTTTGAAGCTAAATTTGAAGGTGTTTACTTAGAAGTAGATAATAAAAATgaag   <  1:433048/141‑1 (MQ=255)
                                                                                                                       agaaTGATATTTATCAGCTCATTATGTTAAACCAGAATCGCAACATACAGGTTATGGTACAGCATTGTTAAATGAAGGATTATCACGTTTTGAAGATAAATTTGAAGGTGTTTACTTAGAAGTAGATAATAAAAATgaaga  <  1:649996/141‑1 (MQ=255)
                                                                                                                            |                                                                                                                                       
CACTTGTTAGAAGCTTCATATAATGAACATCATTTAAAGAAAAGACTTCAAGAACAATTATTCTTAGTCGTTGAAGAAGGTAATGACATCGTTGGCTTTGCTAACTTTATTTACGGTGAAGAATTATATTTATCAGCTCATTATGTTAAACCAGAATCGCAACATACAGGTTATGGTACAGCATTGTTAAATGAAGGATTATCACGTTTTGAAGATAAATTTGAAGGTGTTTACTTAGAAGTAGATAATAAAAATGAAGA  >  USA300TCH1516_ALE/665268‑665527

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