Predicted mutation
evidence seq id position mutation annotation gene description
RA USA300TCH1516_ALE 665,392 T→G 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,3920TG79.0% 145.7 / 27.4 57L79* (TTA→TGA) paiA_1Spermidine/spermine N(1)‑acetyltransferase
Reads supporting (aligned to +/- strand):  ref base T (10/2);  new base G (26/19);  total (36/21)
Fisher's exact test for biased strand distribution p-value = 1.77e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

CAACAGTTAATCACTTGTTAGAAGCTTCATATAATGAACATCATTTAAAGAAAAGACTTCAAGAACAATTATTCTTAGTCGTTGAAGAAGGTAATGACATCGTTGGCTTTGCTAACTTTATTTACGGTGAAGAATTATATTTATCAGCTCATTATGTTAAACCAGAATCGCAACATACAGGTTATGGTACAGCATTGTTAAATGAAGGATTATCACGTTTTGAAGATAAATTTGAAGGTGTTTACTTAGAAGTAGATAATAAAAATGAA  >  USA300TCH1516_ALE/665257‑665525
                                                                                                                                       |                                                                                                                                     
caacaGTTAATCACTTGTTAGAAGCTTCATATAATGAACATCATTTAAAGAAAAGACTTCAAGAACAATTATTCTTAGTCGTTGAAGAAGGTAATGACATCGTTGGCTTTGCTAACTTTATTTACGGTGAAGAATGATAtt                                                                                                                                  >  1:430390/1‑141 (MQ=255)
caacaGTTAATCACTTGTTAGAAGCTTCATATAATGAACATCATTTAAAGAAAAGACTTCAAGAACAATTATTCTTAGTCGTTGAAGAAGGTAATGACATCGTTGGCTTTGCTAACTTTATTTACGGTGAAGAATGATAtt                                                                                                                                  >  1:87660/1‑141 (MQ=255)
  acaGTTAATCACTTGTTAGAAGCTTCATATAATGAACATCATTTAAAGAAAAGACTTCAAGAACAATTATTCTTAGTCGTTGAAGAAGGTAATGACATCGTTGGCTTTGCTAACTTTATTTACGGTGAAGAATTATATTTa                                                                                                                                >  1:916365/1‑141 (MQ=255)
  acaGTTAATCACTTGTTAGAAGCTTCATATAATGAACATCATTTAAAGAAAAGACTTCAAGAACAATTATTCTTAGTCGTTGAAGAAGGTAATGACATCGTTGGCTTTGCTAACTTTAGTTACGGTGAAGAATGATATTTa                                                                                                                                >  2:351902/1‑141 (MQ=255)
               tGTTAGAAGCTTCATATAATGAACATCATTTAAAGAAAAGACTTCAAGAACAATTATTCTTAGTCGTTGAAGAAGGTAATGACATCGTTGGCTTTGCTAACTTTATTTACGGTGAAGAATGATATTTATCAGCTCATTATg                                                                                                                   <  1:178763/141‑1 (MQ=255)
               tGTTAGAAGCTTCATATAATGAACATCATTTAAAGAAAAGACTTCAAGAACAATTATTCTTAGTCGTTGAAGAAGGTAATGACATCGTTGGCTTTGCTAACTTTATTTACGGTGAAGAATGATATTTATCAGCTCATTATg                                                                                                                   >  1:865764/1‑141 (MQ=255)
                gTTAGAAGCTTCATATAATGAACATCATTTAAAGAAAAGACTTCAAGAACAATTATTCTTAGTCGTTGAAGAAGGTAATGACATCGTTGGCTTTGCTAACTTTATTTACGGTGAAGAATGATATTTATCAGCTCATTATGt                                                                                                                  <  1:537178/141‑1 (MQ=255)
                 ttAGAAGCTTCATATAATGAACATCATTTAAAGAAAAGACTTCAAGAACAATTATTCTTAGTCGTTGAAGAAGGTAATGACATCGTTGGCTTTGCTAACTTTATTTACGGTGAAGAATGATATTTATCAGCTCATTATGtt                                                                                                                 >  1:955085/1‑141 (MQ=255)
                      aGCTTCATATAATGAACATCATTTAAAGAAAAGACTTCAAGAACAATTATTCTTAGTCGTTGAAGAAGGTAATGACATCGTTGGCTTTGCTAACTTTATTTACGGTGAAGAATGATATTTATCAGCTCATTATGTTAAAcc                                                                                                            <  1:233892/141‑1 (MQ=255)
                       gCTTCATATAATGAACATCATTTAAAGAAAAGACTTCAAGAACAATTATTCTTAGTCGTTGAAGAAGGTAATGACATCGTTGGCTTTGCTAACTTTATTTACGGTGAAGAATGATATTTATCAGCTCATTATGTTAAACCa                                                                                                           <  2:223119/141‑1 (MQ=255)
                          tCATATAATGAACATCATTTAAAGAAAAGACTTCAAGAACAATTATTCTTAGTCGTTGAAGAAGGTAATGACATCGTTGGCTTTGCTAACTTTATTTACGGTGAAGAATGATATTTATCAGCTCATTATGTTAAACCAGaa                                                                                                        >  1:402520/1‑141 (MQ=255)
                            atatAATGAACATCATTTAAAGAAAAGACTTCAAGAACAATTATTCTTAGTCGTTGAAGAAGGTAATGACATCGTTGGCTTTGCTAACTTTATTTACGGTGAAGAATGATATTTATCAGCTCATTATGTTAAACCAGAATc                                                                                                      <  1:585845/141‑1 (MQ=255)
                            atatAATGAACATCATTTAAAGAAAAGACTTCAAGAACAATTATTCTTAGTCGTTGAAGAAGGTAATGACATCGTTGGCTTTGCTAACTTTATTTACGGTGAAGAATGATATTTATCAGCTCATTATGTTAAACCAGAATc                                                                                                      >  1:523967/1‑141 (MQ=255)
                             tataATGAACATCATTTAAAGAAAAGACTTCAAGAACAATTATTCTTAGTCGTTGAAGAAGGTAATGACATCGTTGGCTTTGCTAACTTTATTTACGGTGAAGAATGATATTTATCAGCTCATTATGTTAAACCAGAATCg                                                                                                     >  1:185392/1‑141 (MQ=255)
                              ataATGAACATCATTTAAAGAAAAGACTTCAAGAACAATTATTCTTAGTCGTTGAAGAAGGTAATGACATCGTTGGCTTTGCTAACTTTATTTACGGTGAAGAATTATATTTATCAGCTCATTATGTTAAACCAGAATCGc                                                                                                    >  1:948488/1‑141 (MQ=255)
                              ataATGAACATCATTTAAAGAAAAGACTTCAAGAACAATTATTCTTAGTCGTTGAAGAAGGTAATGACATCGTTGGCTTTGCTAACTTTATTTACGGTGAAGAATGATATTTATCAGCTCATTATGTTAAACCAGAATCGc                                                                                                    >  2:157181/1‑141 (MQ=255)
                                  tGAACATCATTTAAAGAAAAGACTTCAAGAACAATTATTCTTAGTCGTTGAAGAAGGTAATGACATCGTTGGCTTTGCTAACTTTATTTACGGTGAAGAATTATATTTATCAGCTCATTATGTTAAACCAGAATCGCAACa                                                                                                >  1:969305/1‑141 (MQ=255)
                                             tAAAGAAAAGACTTCAAGAACAATTATTCTTAGTCGTTGAAGAAGGTAATGACATCGTTGGCTTTGCTAACTTTATTTACGGTGAAGAATGATATTTATCAGCTCATTATGTTAAACCAGAATCGCAACATACAGGTTATg                                                                                     >  2:138108/1‑141 (MQ=255)
                                              aaagaaaagaCTTCAAGAACAATTATTCTTAGTCGTTGAAGAAGGTAATGACATCGTTGGCTTTGCTAACTTTATTTACGGTGAAGAATGATATTTATCAGCTCATTATGTTAAACCAGAATCGCAACATACAGGTTATgg                                                                                    <  1:413042/141‑1 (MQ=255)
                                              aaagaaaagaCTTCAAGAACAATTATTCTTAGTCGTTGAAGAAGGTAATGACATCGTTGGCTTTGCTAACTTTATTTACGGTGAAGAATGATATTTATCAGCTCATTATGTTAAACCAGAATCGCAACATACAGGTTATgg                                                                                    <  2:250747/141‑1 (MQ=255)
                                                     agaCTTCAAGAACAATTATTCTTAGTCGTTGAAGAAGGTAATGACATCGTTGGCTTTGCTAACTTTATTTACGGTGAAGAATGATATTTATCAGCTCATTATGTTAAACCAGAATCGCAACATACAGGTTATGGTACAGCa                                                                             >  2:384848/1‑141 (MQ=255)
                                                     agaCTTCAAGAACAATTATTCTTAGTCGTTGAAGAAGGTAATGACATCGTTGGCTTTGCTAACTTTATTTACGGTGAAGAATGATATTTATCAGCTCATTATGTTAAACCAGAATCGCAACATACAGGTTATGGTACAGCa                                                                             >  2:384991/1‑141 (MQ=255)
                                                     agaCTTCAAGAACAATTATTCTTAGTCGTTGAAGAAGGTAATGACATCGTTGGCTTTGCTAACTTTATTTACGGTGAAGAATGATATTTATCAGCTCATTATGTTAAACCAGAATCGCAACATACAGGTTATGGTACAGCa                                                                             <  1:663484/141‑1 (MQ=255)
                                                            aaGAACAATTATTCTTAGTCGTTGAAGAAGGTAATGACATCGTTGGCTTTGCTAACTTTATTTACGGTGAAGAATTATATTTATCAGCTCATTATGTTAAACCAGAATCGCAACATACAGGTTATGGTACAGCATTGTTaa                                                                      >  2:29507/1‑141 (MQ=255)
                                                             aGAACAATTATTCTTAGTCGTTGAAGAAGGTAATGACATCGTTGGCTTTGCTAACTTTATTTACGGTGAAGAATGATATTTATCAGCTCATTATGTTAAACCAGAATCGCAACATACAGGTTATGGTACAGCATTGTTaaa                                                                     >  2:403810/1‑141 (MQ=255)
                                                                aCAATTATTCTTAGTCGTTGAAGAAGGTAATGACATCGTTGGCTTTGCTAACTTTATTTACGGTGAAGAATGATATTTATCAGCTCATTATGTTAAACCAGAATCGCAACATACAGGTTATGGTACAGCATTGTTAAATGa                                                                  <  2:312233/141‑1 (MQ=255)
                                                                aCAATTATTCTTAGTCGTTGAAGAAGGTAATGACATCGTTGGCTTTGCTAACTTTATTTACGGTGAAGAATGATATTTATCAGCTCATTATGTTAAACCAGAATCGCAACATACAGGTTATGGTACAGCATTGTTAAATGa                                                                  <  2:312221/141‑1 (MQ=255)
                                                                  aaTTATTCTTAGTCGTTGAAGAAGGTAATGACATCGTTGGCTTTGCTAACTTTATTTACGGTGAAGAATTATATTTATCAGCTCATTATGTTAAACCAGAATCGCAACATACAGGTTATGGTACAGCATTGTTAAATGAAg                                                                >  1:74999/1‑141 (MQ=255)
                                                                  aaTTATTCTTAGTCGTTGAAGAAGGTAATGACATCGTTGGCTTTGCTAACTTTATTTACGGTGAAGAATGATATTTATCAGCTCATTATGTTAAACCAGAATCGCAACATACAGGTTATGGTACAGCATTGTTAAATGAAg                                                                >  2:252133/1‑141 (MQ=255)
                                                                  aaTTATTCTTAGTCGTTGAAGAAGGTAATGACATCGTTGGCTTTGCTAACTTTATTTACGGTGAAGAATGATATTTATCAGCTCATTATGTTAAACCAGAATCGCAACATACAGGTTATGGTACAGCATTGTTAAATGAAg                                                                >  1:449026/1‑141 (MQ=255)
                                                                  aaTTATTCTTAGTCGTTGAAGAAGGTAATGACATCGTTGGCTTTGCTAACTTTATTTACGGTGAAGAATGATATTTATCAGCTCATTATGTTAAACCAGAATCGCAACATACAGGTTATGGTACAGCATTGTTAAATGAAg                                                                >  1:86325/1‑141 (MQ=255)
                                                                   attattCTTAGTCGTTGAAGAAGGTAATGACATCGTTGGCTTTGCTAACTTTATTTACGGTGAAGAATGATATTTATCAGCTCATTATGTTAAACCAGAATCGCAACATACAGGTTATGGTACAGCATTGTTAAATGAAgg                                                               >  2:408046/1‑141 (MQ=255)
                                                                       ttCTTAGTCGTTGAAGAAGGTAATGACATCGTTGGCTTTGCTAACTTTATTTACGGTGAAGTATTATATATATCGTCTTGGTATGTTAGATCAGATTCGCAACATACAGGTTATGGTACAGCATTGTTAAATGAAGGATTa                                                           >  2:357844/1‑141 (MQ=255)
                                                                        tCTTAGTCGTTGAAGAAGGTAATGACATCGTTGGCTTTGCTAACTTTATTTACGGTGAAGAATGATATTTATCAGCTCATTATGTTAAACCAGAATCGCAACATACAGGTTATGGTACAGCATTGTTAAATGAAGGATTAt                                                          >  1:408828/1‑141 (MQ=255)
                                                                         cTTAGTCGTTGAAGAAGGTAATGACATCGTTGGCTTTGCTAACTTTATTTACGGTGAAGAATTATATTTATCAGCTCATTATGTTAAACCAGAATCGCAACATACAGGTTATGGTACAGCATTGTTAAATGAAGGATTATc                                                         >  1:167649/1‑141 (MQ=255)
                                                                          ttAGTCGTTGAAGAAGGTAATGACATCGTTGGCTTTGCTAACTTTATTTACGGTGAAGAATTATATTTATCAGCTCATTATGTTAAACCAGAATCGCAACATACAGGTTATGGTACAGCATTGTTAAATGAAGGATTATCa                                                        <  1:718589/141‑1 (MQ=255)
                                                                          ttAGTCGTTGAAGAAGGTAATGACATCGTTGGCTTTGCTAACTTTATTTACGGTGAAGAATGATATTTATCAGCTCATTATGTTAAACCAGAATCGCAACATACAGGTTATGGTACAGCATTGTTAAATGAAGGATTATCa                                                        <  1:481696/141‑1 (MQ=255)
                                                                                   gaagaaGGTAATGACATCGTTGGCTTTGCTAACTTTATTTACGGTGAAGAATGATATTTATCAGCTCATTATGTTAAACCAGAATCGCAACATACAGGTTATGGTACAGCATTGTTAAATGAAGGATTATCACGTTTTGaa                                               >  2:414733/1‑141 (MQ=255)
                                                                                    aagaagGTAATGACATCGTTGGCTTTGCTAACTTTATTTACGGTGAAGAATGATATTTATCAGCTCATTATGTTAAACCAGAATCGCAACATACAGGTTATGGTACAGCATTGTTAAATGAAGGATTATCACGTTTTGAAg                                              >  1:444248/1‑141 (MQ=255)
                                                                                     agaagGTAATGACATCGTTGGCTTTGCTAACTTTATTTACGGTGAAGAATGATATTTATCAGCTCATTATGTTAAACCAGAATCGCAACATACAGGTTATGGTACAGCATTGTTAAATGAAGGATTATCACGTTTTGAAGa                                             <  1:405715/141‑1 (MQ=255)
                                                                                      gaagGTAATGACATCGTTGGCTTTGCTAACTTTATTTACGGTGAAGAATTATATTTATCAGCTCATTATGTTAAACCAGAATCGCAACATACAGGTTATGGTACAGCATTGTTAAATGAAGGATTATCACGTTTTGAAGAt                                            >  1:164827/1‑141 (MQ=255)
                                                                                      gaagGTAATGACATCGTTGGCTTTGCTAACTTTATTTACGGTGAAGAATGATATTTATCAGCTCATTATGTTAAACCAGAATCGCAACATACAGGTTATGGTACAGCATTGTTAAATGAAGGATTATCACGTTTTGAAGAt                                            <  1:748582/141‑1 (MQ=255)
                                                                                         ggTAATGACATCGTTGGCTTTGCTAACTTTATTTACGGTGAAGAATGATATTTATCAGCTCATTATGTTAAACCAGAATCGCAACATACAGGTTATGGTACAGCATTGTTAAATGAAGGATTATCACGTTTTGAAGATaaa                                         <  1:494206/141‑1 (MQ=255)
                                                                                          gTAATGACATCGTTGGCTTTTCTAACTTTATTTACGGTGAAGAATGATATTTATCAGCTCATTATGTTAAACCAGAATCGCAACATACAGGTTATGGTACAGCATTGTTAAATGAAGGATTATCACGTTTTGAAGATAAAt                                        >  1:184535/1‑141 (MQ=255)
                                                                                             aTGACATCGTTGGCTTTGCTAACTTTATTTACGGTGAAGAATTATATTTATCAGCTCATTATGTTAAACCAGAATCGCAACATACAGGTTATGGTACAGCATTGTTAAATGAAGGATTATCACGTTTTGAAGATAAATTTg                                     >  1:207917/1‑141 (MQ=255)
                                                                                              tGACATCGTTGGCTTTGCTAACTTTATTTACGGTGAAGAATTATATTTATCAGCTCATTATGTTAAACCAGAATCGCAACATACAGGTTATGGTACAGCATTGTTAAATGAAGGATTATCACGTTTTGAAGATAAATTTGa                                    >  1:486313/1‑141 (MQ=255)
                                                                                              tGACATCGTTGGCTTTGCTAACTTTATTTACGGTGAAGAATGATATTTATCAGCTCATTATGTTAAACCAGAATCGCAACATACAGGTTATGGTACAGCATTGTTAAATGAAGGATTATCACGTTTTGAAGATAAATTTGa                                    <  1:13546/141‑1 (MQ=255)
                                                                                                   tCGTTGGCTTTGCTAACTTTATTTACGGTGAAGAATGATATTTAGCAGCTCATTATGTTAAACCAGAATCGCAACATACAGGTTATGGTACAGCATTGTTAAATGAAGGCTTATCACGTTTTGAAGATAAATTTGGAGgtg                               >  1:774681/1‑141 (MQ=255)
                                                                                                     gTTGGCTTTGCTAACTTTATTTACGGTGAAGAATGATATTTATCAGCTCATTATGTTAAACCAGAATCGCAACATACAGGTTATGGTACAGCATTGTTAAATGAAGGATTATCACGTTTTGAAGATAAATTTGAAGGTGtt                             <  1:486683/141‑1 (MQ=255)
                                                                                                       tGGCTTTGCTAACTTTATTTACGGTGAAGAATTATATTTATCAGCTCATTATGTTAAACCAGAATCGCAACATACAGGTTATGGTACAGCATTGTTAAATGAAGGATTATCACGTTTTGAAGATAAATTTGAAGGTGTTTa                           <  1:113047/141‑1 (MQ=255)
                                                                                                       tGGCTTTGCTAACTTTATTTACGGTGAAGAATGATATTTATCAGCTCATTATGTTAAACCAGAATCGCAACATACAGGTTATGGTACAGCATTGTTAAATGAAGGATTATCACGTTTTGAAGATAAATTTGAAGGTGTTTa                           >  1:790714/1‑141 (MQ=255)
                                                                                                              gCTAACTTTATTTACGGTGAAGAATGATATTTATCAGCTCATTATGTTAAACCAGAATCGCAACATACAGGTTATGGTACAGCATTGTTAAATGAAGGATTTTCACGTTTTGAAGATAAATTTGAAGGTGTTTACTTAGaa                    >  2:358917/1‑141 (MQ=255)
                                                                                                                   cTTTATTTACGGTGAAGAATGATATTTATCAGCTCATTATGTTAAACCAGAATCGCAACATACAGGTTATGGTACAGCATTGTTAAATGAAGGATTATCACGTTTTGAAGATAAATTTGAAGGTGTTTACTTAGAAGTAGa               <  1:290814/141‑1 (MQ=255)
                                                                                                                         ttaCGGTGAAGAATGANATTTATCAGCTCNTTATGTTAAACCAGAATCGCAACATACAGGTTATGGTACAGCATTGTTAAATGAAGGATTATCACGTTTTGAAGATAAATTTGAAGGTGTTTACTTAGAAGTAGAtaataa         >  2:256176/1‑141 (MQ=255)
                                                                                                                           aCGGTGAAGAATGATATTTATCAGCTCATTATGTTAAACCAGAATCGCAACATACAGGTTATGGTACAGCATTGTTAAATGAAGGATTATCACGTTTTGAAGATAAATTTGAAGGTGTTTACTTAGAAGTAGATAATaaaa       <  1:340163/141‑1 (MQ=255)
                                                                                                                               tGAAGAATGATATTTATCAGCTCATTATGTTAAACCAGAATCGCAACATACAGGTTATGGTACAGCATTGTTAAATGAAGGATTATCACGTTTTGAAGATAAATTTGAAGGTGTTTACTTAGAAGTAGATAATAAAAATga   <  1:784999/141‑1 (MQ=255)
                                                                                                                                gaagaaTGATATTTATCAGCTCATTATGTTAAACCAGAATCGCAACATACAGGTTATGGTACAGCATTGTTAAATGAAGGATTATCACGTTTTGAAGATAAATTTGAAGGTGTTTACTTAGAAGTAGATAATAAAAATgaa  >  1:249922/1‑141 (MQ=255)
                                                                                                                                       |                                                                                                                                     
CAACAGTTAATCACTTGTTAGAAGCTTCATATAATGAACATCATTTAAAGAAAAGACTTCAAGAACAATTATTCTTAGTCGTTGAAGAAGGTAATGACATCGTTGGCTTTGCTAACTTTATTTACGGTGAAGAATTATATTTATCAGCTCATTATGTTAAACCAGAATCGCAACATACAGGTTATGGTACAGCATTGTTAAATGAAGGATTATCACGTTTTGAAGATAAATTTGAAGGTGTTTACTTAGAAGTAGATAATAAAAATGAA  >  USA300TCH1516_ALE/665257‑665525

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