Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 493,809 G→T 11.4% V82L (GTA→TTA)  bltD_1 → Spermine/spermidine acetyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE493,8090GT11.4% 86.6 / 3.1 35V82L (GTA→TTA) bltD_1Spermine/spermidine acetyltransferase
Reads supporting (aligned to +/- strand):  ref base G (13/18);  new base T (2/2);  total (15/20)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.24e-01

CGCCAGAGGAAAATATTACAGCTGCAATGTCTGATAATGAAAGGTATCCCATCGTTGTAATGGATGGCAGGCAATGTGTGGCCTTTTTTACATTACATCGTGGAAAAGGGGTCGCACCATTTAGCGATAACCAAGATGCAGTATTTTTCAGGTCATTTAGTGTTGATCAACGTTATCGTAATAGAGGAATAGGTAAAGTGGTAATGGAAAAATTGGCGTCATTTATCACTTCAACATTTCAGGATATTAATGAGATTGTGTTAACGGTTAATACTGACAAT  >  USA300TCH1516_ALE/493669‑493949
                                                                                                                                            |                                                                                                                                            
ggCCATATGAAGATGTGACAGGTGCAATGTATGATAATGAAAGGTAGCCCATCGTTGTAATGGATGGCAGGCAATGTGTGGCCTTTGTTACATTACATCGGTGAAAAGGGGTCGCACCATTTAGCGATAACCAAGATGCAg                                                                                                                                              <  2:196963/140‑1 (MQ=255)
cGCCAGAGGAACATATTACAGCTGCAATGTCTGATAATGAATGGTATCCCATCGTTGTAATGGATGGCAGGCAATGTGTGGCCTTTTTTACATTACATCGTGGAAAAGGGGTCGCACCATTTAGCGATAACCAAGATGCAg                                                                                                                                              <  2:309830/141‑1 (MQ=255)
 gCCAGAGGAAAATATTACAGCTGCAATGTATGATAATGAAAGGTATCCCATCGTTGTAATGGATGGCAGGCAATGTGTGGCCTTTTTTACATTACATCGTGGAAAAGGGGTCGCACCATTTAGCGATAACCAAGATGCAGt                                                                                                                                             <  2:173619/141‑1 (MQ=255)
     gagGAAAATATTACAGCTGCAATGTCTGATGATGAAAGGTATCCCATCGTTGTAATGGATGGCAGGCAATGTGTGGCCTTTTTTACATTACATCGTGGAAAAGGGGTCGCACCATTTAGCGATAACCAAGCTGCAGTAttt                                                                                                                                         >  1:207113/1‑141 (MQ=255)
      agGAAAATATTACAGCTGCAATGTCTGATAATGAAAGGTATCCCATCGTTGTAATGGATGGCAGGCAATGTGTGGCCTTTTTTACATTACATCGTGGAAAAGGGGTCGCACCATTTAGCGATAACCAAGATGCAGTAtttt                                                                                                                                        <  2:336274/141‑1 (MQ=255)
           aaTATTACAGCTGCAATGTCTGATAATGAAAGGTATCCCATCGTTGTAATGGATGGCAGGCAATGTGTGGCCTTTTTTACATTACATCGTGGAAAAGGGGTCGCACCATTTAGCGATAACCAAGATGCAGTATTTTTCAgg                                                                                                                                   <  1:93390/141‑1 (MQ=255)
                    gCTGCAATGTCTGATAATGAAAGGTATCCCATCGTTGTAATGGATGGCAGGCAATGTGTGGCCTTTTTTACATTACATCGTGGAAAAGAGGGGTCGCCATTTAGCGATAACCAAGATGTAGTATTTTTAGGGTCATTTAga                                                                                                                          >  2:293641/1‑140 (MQ=255)
                          aTGTCTGATAATGAAAGGTATCCCATCGTTGTAATGGATGGCAGGCAATGTGTGGCCTTTTTTACATTACATCGTGGAAAAGGGGTCGCACCATTTAGCGATAACCAAGATGCAGTATTTTCCCGGTCATTTAGTGTTGAt                                                                                                                    >  2:164679/1‑141 (MQ=255)
                                      gAAAGGTATCCCATCGTGGTAATGGATGGCAGGCAATGTGTGGCCTTTTTTACATTACATCGTGGAAAAGGGGTCGCACCATTTAGCGATAACCAAGATGCAGTATTTTTCAGGTCATTTAGTGTTGATCAACGTTATCGt                                                                                                        <  2:67453/141‑1 (MQ=255)
                                        aaGGTATCCCATCGTTGTAATGGATGGCAGGCAATGTGTGGCCTTTTTTACATTACATCGTGGAAAAGGGGTCGCACCCTTTAGCGATAACCAAAATGCCGTATTTTTCCGGTCATTTAGCGTTGATCAACGTTATCGTaa                                                                                                      >  1:38046/1‑141 (MQ=255)
                                        aaGGTATCCCATCGTTGTAATGGATGGCAGGCAATGTGTGGCCTTTTTTACATTACATCGTGGAAAAGGGGTCGCACCATTTAGCGATAACCAAGATGCAGTATTTTTCAGGTCATTTAGTGTTGATCAACGTTATCGTaa                                                                                                      <  1:34885/141‑1 (MQ=255)
                                          ggTATCCCATCGTTGTAATGGATGGCAGGCAATGTGTGGCCTTTTTTTTTTTTTTTCGTGGAAAAGAGGTCGCACCATTTAGCGATATCGAAGCTGCAGTATTTTTCAGGTCATTTAGTGTTGATCAACGTTATCGtctta                                                                                                    >  2:26359/1‑137 (MQ=255)
                                             tgccgttCGTTGTAATGGATGGTAGGCAATGTGTGTCCTGTTTTAGATTACATCGTGGAAAAGGGGTCGCACCATTTAGCGATAACCAAGATGCAGTATTTTTCAGTTTCTTAGGTATTTATCTACGTTATCGTAATagag                                                                                                 <  2:135557/135‑1 (MQ=255)
                                                 cATCGTTGTAATGGATGGCAGGCAATGTGTGGCCTTTTTTACATTACATCGTGGAAAAGGGGTCGCACCATTTAGCGATAACCAAGATGCAGTATTTTTCAGGTCATTTAGTGTTGATCAACGTTATCGTAATAGAGGAAt                                                                                             <  1:217732/141‑1 (MQ=255)
                                                    cGTTGTAATGGATGGCAGGCAATGTGTGGCCTTTTTTACATTACATCGTGGAAAAGGGGTCGCACCATTTAGCGATAACCAAGATGCAGTATTTTTCAGGTCATTTAGTGTTGATCAACGCTACCGTACTAGAGGCATAgg                                                                                          >  1:58510/1‑141 (MQ=255)
                                                        gTAATGGATGGCAGGCAATGTGTGGCCTTTTTTACATTACATCGTGGAAAAGGGGTCGCACCATTTAGCGATAACCAAGATGCAGTATTTTTCAGGTCATTTAGTGTTGATCAACGTTATCGTAATAGAGGAATAGGTaaa                                                                                      <  1:419103/141‑1 (MQ=255)
                                                                tggCAGGCAATGTGTGGCCTTTTTTACATTACATCGTGGAAAAGGGGTCGCAACATTTAGCGATAACCAAGATGCAGTTTTTTTCAGGTCATTTAGTGTTGATCAACGTTATCGTAATAGAGGAATAGGTAAAGTGGTAAt                                                                              <  1:336041/141‑1 (MQ=255)
                                                                  gcaggcaATGTGTGGCCTTTTTTACATTACATCGTGGAAAAGGGGTCGCACCATTTAGCGATAACCAAGATGCAGTATTTTTCAGGTCATTTAGTGTTGATCAACGTTATCGTAATAGAGGAATAGGTAAAGTGGTAATgg                                                                            >  1:393883/1‑141 (MQ=255)
                                                                   caggcaATGTGTGGCCTTTTTTACATTACATCGTGGAAAAGGGGTCGGATCATTTAGCGATAACCAAGATGCAGTATTTTTCAGGTCATTTAGTGTTGATCAACGTTATCGTAATAGAGGAATAGGTAAAGTGGTAATGGa                                                                           <  2:234871/141‑1 (MQ=255)
                                                                       caATGTGTGGCGTTTTTTACATTACATCGGGGAAAAGGGGTCGTACCATTTAGCGATAACCAAGATGCAGTATTTTTCAGGTCATTTAGTGTTGATCAACGTTATCGTAATAGAGGAATAGGTAAAGTGGTAATGGaaaaa                                                                       <  1:393267/141‑1 (MQ=255)
                                                                              tgGCCTTTTTTACATTACATCGTGGAAAAGGGGTCGCACCATTTAGCGATAACCAAGATGCATTATTTTTCAGGTCATTTAGTGTTGATCAACGTTATCGTAATAGAGGAATAGGTAAAGTGGTAATGGAAAAATTGGCGt                                                                >  1:123209/1‑141 (MQ=255)
                                                                              tgGCCTTTTTTACATTACATCGTGGAAAAGGGGTCGCACCATTTAGCGATAACCAAGATGCATTATTTTTCAGGTCATTTAGTGTTGATCAACGTTATCGTAATAGAGGAATAGGTAAAGTGGTAATGGAAAAATTGGCGt                                                                >  1:414269/1‑141 (MQ=255)
                                                                                    tttttACATTACATCGTGGAAAAGGGGTCGCACCATTTAGCGATAACCAAGATGCAGTATTTTTCAGGTCATTTAGTGTTGATCAACGTTATCGTAATAGAGGAATAGGTACAGTGGTAATGGAAAAATTGGCGTCATTTa                                                          >  2:360042/1‑141 (MQ=255)
                                                                                      tttACATTACATCGTGGAAAAGGGGTCGCACCATTTAGCGATAACCAAGATGCAGTATTTTTCAGGTCATTTAGTGTTGATCAACGTTATCGTAATAGAGGAATAGGTAAAGTGGTAATGGAAAAATTGGCGTCATTTATc                                                        <  1:86230/141‑1 (MQ=255)
                                                                                          aattacgtAGTGGAAAATGTGTCGCAACATTTAGTGATAACGAAGAAGTAGTATTTTTCAGGTTATTTAGTGTTGATCAAAGTTATCGTAATAGTGGAATTGGTAAAGTGGTAATGGAAAAATTGTCATCATTTAACACtt                                                    <  2:266879/140‑1 (MQ=255)
                                                                                           attacatCGTGGAAAAGGGGTCGCACCATTTAGCGATAACCAAGATGCAGTATTTTTCAGGTCATTTAGTGTTGATCAACGTTATCGTAATAGAGGAATAGGTAAAGTGGTAATGGAAAAATTGGCGTCATTTATCACTTc                                                   <  2:294774/141‑1 (MQ=255)
                                                                                            ttacatCGTGGAAAAGGGGTCGCACCATTTAGCGATAACCAAGATGCAGTATTTTTCAGGTCATTTAGTGTTGATCAACGTTATCGTAATAGGGGCATAGGTACAGTGGTAATGGAACAATTGGCGTCATTTATCACTTCa                                                  >  1:148514/1‑141 (MQ=255)
                                                                                                 ggggggAAGAGGGGGCGGAACATGTAGCGAGAACCAAGGGGCAGTATTGTTCAGGTGATTTAGGGCTGATCAAGGTTAGCGTAATAGAGGAAGAGGTAAAGTGGTAATGGAAAAATTGGCGTCATTTATCACTTCAACAtt                                             <  2:127391/137‑1 (MQ=255)
                                                                                                       aaaaGGGGTCGCACCATTTAGCGATAACCAAGATGCAGTATTTTTAAGGTCATTTAGTGTTGATCAACGTTATCGTAATAGAGGAATCGGTAAAGTGGTACTGGAAAAATTGGAGTCATACATCACCTCAACATTTCAAGa                                       >  2:267621/1‑141 (MQ=255)
                                                                                                       aaaaGGGGTCGCACCATGTAGCGATAACCAAGATGCAGTATTTTTCAGGTCATTTAGTGTTGATCAACGTTATCGTAATAGAGGAATAGGTAAAGTGGTAATGGAAAAATTGGCGTCATTTATCACTTCAACATTTCAGGa                                       <  1:39591/141‑1 (MQ=255)
                                                                                                                     cATTTAGCGATAACCAAGATGCAGTATTTTTCAGGTCATTTAGTGTTGATCAACGTTATCGTAATAGAGGAATAGGTAAAGTGGTAATGGAAAAATTGGCGTCATTTATCACTTCAACATTTCAGGATATTAATGAGATtg                         <  1:360042/141‑1 (MQ=255)
                                                                                                                          aGCGATAACCAAGATGCAGTATTTGTCAGGTCATTTAGTGTTGATCAACGTTATCGTAATAGAGGAATAGGTAAAGTGGTAATGGAAAAATTGGCGTCATTTATCACTTCAACATTTCAGGATATTAATGAGATTGTGTTa                    >  1:63301/1‑141 (MQ=255)
                                                                                                                           gCGATAACCAAGATGCAGTATTTTTCAGGTCATTTAGTGTTGATCAACGTTATCGTAATAGAGGAATAGGTAAAGTGGTAATGGAAAAATTGGCGTCATTTATCACTTCAACATTTCAGGATATTAATGAGATTGTGTTaa                   >  1:355974/1‑141 (MQ=255)
                                                                                                                            cGATAACCAAGATGCAGTATTTTTCAGGTCATTTAGTGTTGATCAACGTTATCGTAATAGAGGAATAGGTAAAGTGGTAATGGAAAAATTGGCGTCATTTATCACTTCAACATTTCAGGATATTAATGAGATTGTGTTAAc                  <  1:192264/141‑1 (MQ=255)
                                                                                                                                 aCCAAGATGCATTATTTTTCAGGTCATTTAGTGTTGATCAACGTTATCGTAATAGAGGAATAGGTAAAGTGGTAATGGAAAAATTGGCGTCATTTATCACTTCAACATTTCAGGATATTAATGAGATTGTGTTAACGGTTa             <  2:414269/141‑1 (MQ=255)
                                                                                                                                 aCCAAGATGCATTATTTTTCAGGTCATTTAGTGTTGATCAACGTTATCGTAATAGAGGAATAGGTAAAGTGGTAATGGAAAAATTGGCGTCATTTATCACTTCAACATTTCAGGATATTAATGAGATTGTGTTAACGGTTa             <  2:123209/141‑1 (MQ=255)
                                                                                                                                      gATGCAGTATTTTTCAGGTCATTTAGTGTTGATCAACGTTATCGTAATAGAGGAATAGGTAAAGTGGTAATGGAAAAATTGGCGTCATTTATCACTTCAACATTTCAGGATATTAATGAGATTGTGTTAACGGTTAATACt        >  1:426742/1‑141 (MQ=255)
                                                                                                                                            gTATTTTTCAGGTCATTTAGTGTTGATGAATGTTATCGTAATAGAGGAATAGGTAAAGTGGTAATGGAAAAATTGGCGTCATTTATCACTTCAACATTTCAGGATATTAATGAGATTGTGTTAACGGTTAATACTGACAAt  >  2:375053/1‑141 (MQ=255)
                                                                                                                                            |                                                                                                                                            
CGCCAGAGGAAAATATTACAGCTGCAATGTCTGATAATGAAAGGTATCCCATCGTTGTAATGGATGGCAGGCAATGTGTGGCCTTTTTTACATTACATCGTGGAAAAGGGGTCGCACCATTTAGCGATAACCAAGATGCAGTATTTTTCAGGTCATTTAGTGTTGATCAACGTTATCGTAATAGAGGAATAGGTAAAGTGGTAATGGAAAAATTGGCGTCATTTATCACTTCAACATTTCAGGATATTAATGAGATTGTGTTAACGGTTAATACTGACAAT  >  USA300TCH1516_ALE/493669‑493949

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: