Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 18,890 C→G 100% N182K (AAC→AAG)  gdpP → Cyclic‑di‑AMP phosphodiesterase GdpP

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE18,8900CG97.7% 142.9 / ‑3.7 43N182K (AAC→AAG) gdpPCyclic‑di‑AMP phosphodiesterase GdpP
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base G (29/13);  total (30/13)
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.32e-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.

TCAGTATCATTTCCAAGTACGTTATTCTGAGAATGATCACTGCCTCTATTTCTTTGATATAACTGAACAAGTACAAACAAATGAACTATATGAAAATTCTAAACCAATCATTGCGACATTATTTTTAGATAACTACGATGAGATTACGCAAAATATGAATGATACGCAGCGTTCGGAAATCAACTCAATGGTAACGCGTGTCATTAGTCGATGGGCAACTGAGTATAATATATTTTTCAAAAGATACAGTTCCGATCAATTCGTA  >  USA300TCH1516_ALE/18758‑19022
                                                                                                                                    |                                                                                                                                    
tCAGTATCATTTCCAAGTACGTTATTCTGAGAATGATCACTGCCTCTATTTCTTTGATATAACTGAACAAGTACAAACAAATGAACTATATGAAAATTCTAAACCAATCATTGCGACATTATTTTTAGATAAGTACGATga                                                                                                                              <  1:96305/141‑1 (MQ=255)
tCAGTATCATTTCCAAGTACGTTATTCTGAGAATGATCACTGCCTCTATTTCTTTGATATAACTGAACAAGTACAAACAAATGAACTATATGAAAATTCTAAACCAATCATTGCGACATTATTTTTAGATAAGTACGATga                                                                                                                              <  1:795454/141‑1 (MQ=255)
  aGTATCATTTCCAAGTACGTTATTCTGAGAATGATCACTGCCTCTATTTCTTTGATATAACTGAACAAGTACAAACAAATGAACTATATGAAAATTCTAAACCAATCATTGCGCCATTATTTTTAGATAAGTACGATgaga                                                                                                                            >  1:379901/1‑141 (MQ=255)
  aGTATCATTTCCAAGTACGTTATTCTGAGAATGATCACTGCCTCTATTTCTTTGATATAACTGAACAAGTACAAACAAATGAACTATATGAAAATTCTAAACCAATCATTGCGACATTATTTTTAGATAAGTACGATgaga                                                                                                                            >  1:448713/1‑141 (MQ=255)
         tttCCAAGTACGTTATTCTGAGGCTAAGCACGTCCTCTATTTCTTTGATATAACTGAACAAGTACAAACAAATGAACTATATGAAAATTCTAAACCAATCATTGCGACATTATTTTTAGATAAGTACGATGAGATTACGCa                                                                                                                     >  1:208416/1‑141 (MQ=255)
              aaGTACGTTATTCTGAGAATGATCACTGCCTCTATTTCTTTGATATAACTGAACAAGTACAAACAAATGAACTATATGAAAATTCTAAACCAATCATTGCGACATTATTTTTAGATAAGTACGATGAGATTACGCAAAata                                                                                                                >  1:410408/1‑141 (MQ=255)
                            gagaATGATCACTGCCTCTATTTCTTTGATATAACTGAACAAGTACAAACAAATGAACTATAAGAAAATTCTAAACCAATCATTGCGACATTATTTTTAGATAAGTACGGAGCGATTACGCAAAAAATGAGTGAAACGCaa                                                                                                  >  1:154249/1‑140 (MQ=255)
                                aTGATCACTGCCTCTATTTCTTTGATATAACTGAACAAGTACAAACAAATGAACTATATGAAAATTCTAAACCAATCATTGCGACATTATTTTTAGATAAGTACGATGAGATTACGCAAAATATGAATGATATGCAACGtt                                                                                              >  1:439779/1‑141 (MQ=255)
                                    tCACTGCGTCTATTTCTTTGATATAACCGAACAAGTACAAACAAATGAACTATATGAAAATTCTAAACCAATCATTGCGACATTATTTTTAGATAAGTACCATGAGAGGAGGCAAAATATGAATGATACGCAGCGTTCGGa                                                                                          <  1:123138/141‑1 (MQ=255)
                                    tCACTGCCTCTATTTCTTTGATATAACTGAACAAGTACAAACAAATGAACTATATGAAAATTCTAAACCAATCATTGCGACATTATTTTTAGATAAGTACGATGAGATTACGCAAAATATGAATGATACGCAGCGTTCGGa                                                                                          <  1:730849/141‑1 (MQ=255)
                                    tCACTGCCTCTATTTCTTTGATATAACTGAACAAGTACAAACAAATGAACTATATGAAAATTCTAAACCAATCATTGCGACATTATTTTTAGATAAGTACGATGAGATTACGCAAAATATGAATGATACGCAGCGTTCGGa                                                                                          >  1:24353/1‑141 (MQ=255)
                                     cACTGCCTCTATTTCTTTGATATAACTGAACAAGTACAAACAAATGAACTATATGAAAATTCTAAACCAATCATTGCGACATTATTTTTAGATAAGTACGATGAGATTACGCAAAATATGAATGATACGCAGCGTTCGGaa                                                                                         >  1:277148/1‑141 (MQ=255)
                                             ctATTTCTTTGATATAACTGAACAAGTACAAACAAATGAACTATATGAAAATTCTAAACCAATCATTGCGACATTATTTTTAGATAAGTACGATGAGCTTACGCAAAATATGAATGATACGCAGCGTTCGGAAATCAACTc                                                                                 >  1:456744/1‑141 (MQ=255)
                                             ctATTTCTTTGATATAACTGAACAAGTACAAACAAATGAACTATATGAAAATTCTAAACCAATCATTGCGACATTATTTTTAGATAAGTACGATGAGATTACGCAAAATATGAATGATACGCAGCGTTCGGAAATCAACTc                                                                                 >  1:503037/1‑141 (MQ=255)
                                                  tCTTTGATATAACTGAACAAGTACAAACAAATGAACTATATGAAAATTCTAAACCAATCATTGCGACATTATTTTTAGATAAGTACGATGAGATTACGCAAAATATGAATGATACGCAGCGTTCGGAAATCAACTCAATgg                                                                            >  1:128737/1‑141 (MQ=255)
                                                       gATATAACTGAACAAGTACAAACAAATGAACTATATGAAAATTCTAAACCAATCATTGCGACATTATTTTTAGATAAGTACGATGAGATTACGCAAAATATGAATGATACGCAGCGTTCGGAAATCAACTCAATGGTAAcg                                                                       <  1:472653/141‑1 (MQ=255)
                                                          ctaACTGAACAAGGACAAACAAATGAACTATATGAAATTTCTAAACCAAGCATTGCGACCTTATTTTTGGATAAGTACGATGAGATTACGCAAAATATGAATGATACGCAGCGTTCGGAAATCAACTCAATGGTAACGCgt                                                                    <  1:115423/140‑1 (MQ=255)
                                                                    aaGTACAAACAAATGAACTATATGAAAATTCTAAACCAATCATTGCGACATTATTTTTAGATAAGTACGATGAGATTACGCAAAATATGAATGATACGCAGCGTTCGGAAATCAACTCAATGGTAACGCGTGTCATTAGtg                                                          >  1:307103/1‑140 (MQ=255)
                                                                    aaGTACAAACAAATGAACTATATGAAAATTCTAAACCAATCATTGCGACATTATTTTTAGATAAGTACGATGAGATTACGAAAAATATGAATGATACGCAGCGTTCGGAAATCAACTCAATGGTAACGCGTGTCATTAGTc                                                          >  1:418456/1‑141 (MQ=255)
                                                                       tACAAACAAATGAACTATATGAAAATTCTAAACCAATCATTGCGACATTATTTTTAGATAAGTACGATGAGATTACGCAAAATATGAATGATACGCAGCGTTCGGAAATCAACTCAATGGTAACGCGTGTCATTAGTCGAt                                                       >  1:346192/1‑141 (MQ=255)
                                                                        acaaacaaATGAACTATATGAAAATTCTAAACCAATCATTGCGACATTATTTTTAGATAAGTACGATGAGATTACGCCAAATATGAATGATACGCAGCGTTCGGAAATAAACTCAATGGTAACGCGTGTCATTAGTCGATg                                                      >  1:339782/1‑141 (MQ=255)
                                                                        acaaacaaATGAACTATATGAAAATTCTAAACCAATCATTGCGACATTATTTTTAGATAAGTACGATGAGATTACGCAAAATATGAATGATACGCAGCGTTCGGAAATCAACTCAATGGTAACGCGTGTCATTAGTCGATg                                                      >  1:436244/1‑141 (MQ=255)
                                                                         caaacaaaTGAACTATATGAAAATTCTAAACCAATCATTGCGACATTATTTTTAGATAAGTACGATGAGATTACGCAAAATATGAATGATACGCAGCGTTCGGAAATCAACTCAATGGTAACGCGTGTCATTAGTCGATgg                                                     <  1:308221/141‑1 (MQ=255)
                                                                         aaaaaaaaGGAACTATATGAAAATTCTAAACCAAGCATTGCGACATTAGTTTTAGATAAGTACGATGAGATTACGCAAAATATGAATGATACGCAGCGTTCGCAAATCAACTCAATGGTAACGCTTGTCATTAGTCGATgg                                                     <  1:486371/140‑1 (MQ=255)
                                                                          aaacaaaTGAATTATATGAACATTCTAAACCAATCATTGCGACATTATTTTTAGAGAAGTACGATGAGATTACGCAAAACATGAATGATACGCAGCGTTCGGAAATCAACTCAATGTTAACGCGTGTCATTAGTCGATggg                                                    <  1:558397/141‑1 (MQ=255)
                                                                              aaaTGAACTATATGAAAATTCTAAACCAATCATTGCGACATTATTTTTAGATAAGTACGATGAGATTACGCAAAATATGAATGATACGCAGCGTTCGGAAATCAACTCAATGGTAACGCGTGTCATTAGTCGATGGGCAAc                                                >  1:361229/1‑141 (MQ=255)
                                                                                aTGAACTATATGAAAATTCTAAACCAATCATTGCGACATTATTTTTAGATGAGTACGATGAGATTACGCAAAATATGAATGATACGCAGCGTTCGGAAATCAACTCAATGGTAACGCGTGTCATTAGTCGATGGGCGACTg                                              >  1:504801/1‑141 (MQ=255)
                                                                                aTGAACTATATGAAAATTCTAAACCAATCATTGCGACATTATTTTTAGATAAGTACGATGAGGTTACGCAAAATATGAATGATACGCAGCGTTCGGAAATCAGCTCAATGGTAACGCGTGTCATTAGTCGATGGGCGACTg                                              >  1:502686/1‑141 (MQ=255)
                                                                                aTGAACTATATGAAAATTCTAAACCAATCATTGCGACATTATTTTTAGATAAGTACGATGAGATTACGCAAAACATGAGTGATACGAAGCGCTCGGAAATCAACTCAATGGTAACGCGTGTCATTAGTCGATGGGCAAcgg                                              >  1:115431/1‑139 (MQ=255)
                                                                                  gAACTATATGAAAATTCTAAACCAATCATTGCGACATTATTTTTAGATAAGTACGATGAGATTACGCAAAATATGAATGATACGCAGCGTTCGGAAATCAACTCAATGGTAACGCGTGTCATTAGTCGATGGGCAACTGAg                                            <  1:82847/141‑1 (MQ=255)
                                                                                    aCTATATGAAAATTCTAAACCAATCATTGCGACATTATTTTTAGATAAGTACGATGAGATTACGCCAAATATGAATGATAGGCAGCGTGCGGAAATCAACTCAATGGTAACGCGGGTCATTAGTCGATGGGCAACTGAGta                                          >  1:163115/1‑141 (MQ=255)
                                                                                    aCTATATGAAAATTCTAAACCAATCATTGCGACATTATTTTTAGATAAGTACGATGAGATTACGCAAAGTATGAATGATACGCGGGTAGCGATAATCAACTCAATGGTAACGAGTGTGATTAGTCGATGGGCAACTGAGta                                          >  1:163633/1‑141 (MQ=255)
                                                                                       atatGAAAATTCTAAACCAATCATTGCGACATTATTTTTAGATAAGTACGATGAGATTACGCAAAATATGAATGATACGCAGCGTTCGGAAATCAACTCAATGGTAACGAGGGTCATTAATCGATGGGCAACTGAGTataa                                       >  1:400640/1‑141 (MQ=255)
                                                                                                  cTAAACCAATCATTGCGACATTATTTTTAGATAAGTACGATGAGATTACGCAAAATATGAATGATACGCAGCGTTCGGAAATCAACTCAATGGTAACGCGTGTCATTAGTCGATGGGCAACTGAGTATAATATATTTTTCa                            <  1:362063/141‑1 (MQ=255)
                                                                                                  cTAAACCAATCATTGCGACATTATTTTTAGATAAGTACGATGAGATTACGCAAAATATGAATGATACGCAGCGTTCGGAAATCAACTCAATGGTAACGCGTGTCATTAGTCGATGGGCAACTGAGTATAATATATTTTTCa                            >  1:815063/1‑141 (MQ=255)
                                                                                                           tCATTGCGACATTATTTTTAGATAACTACGATGAGATTACGCAAAATATGAATGATACGCAGCGTTCGGAAATCAACTCAATGGTAACGCGTGTCATTAGTCGATGGGCAACTGAGTATAATATATTTTTCAAAAGATACa                   >  1:188020/1‑141 (MQ=255)
                                                                                                             aTTGCGACATTATTTTTAGATAAGTACGATGAGATTACGCAAAATATGAATGATACGCAGCGTTCGGAAATCAACTCAATGGTAACGCGTGTCATTAGTCGATGGGCAACTGAGTATAATATATTTTTCAAAAGATACAGt                 >  1:745487/1‑141 (MQ=255)
                                                                                                               tGCGACATTATTTTTAGATAAGTACGATGAGATTACGCAAAATATGAATGATACGCAGCGTTCGGAAATCAACTCAATGGTAACGCGTGTCATTAGTCGATGGGCAACTGAGTATAATATATTTTTCAAAAGATACAGGTc               >  1:59756/1‑141 (MQ=255)
                                                                                                                gCGACATTATTTTTAGATAAGTACGATGAGATTACGCAAAATATGAATGATACGCAGCGTTCGGAAATCAACTCAATGGTAACGCGTGTCATTAGTCGATGGGCAACTGAGTATAATATATTTTTCAAAAGATACAGTTcc              <  1:388561/141‑1 (MQ=255)
                                                                                                                     attattTTTAGATAAGTACGATGAGATTACGCAAAATATGAATGATACGCAGCGTTCGGAAATCAACTCAGTGGTAACGCGTGTCATTAGTCGATGGGCAACTGAGTATAATATATTTTTCAAAAGATACGGTGCCGATCa         >  1:794344/1‑141 (MQ=255)
                                                                                                                          ttttAGATAAGTACGATGAGATTACGCAAAATATGAATGATACGCAGCGTTCGGAAATCAACTCAATGGTAACGCGTGTCATTAGTCGATGGGCAACTGAGTATAATATATTTTTCAAAAGATACAGTTCCGATCAATTCg    <  1:740307/141‑1 (MQ=255)
                                                                                                                            ttAGATAAGTACGATGAGATTACGCAAAATATGAATGATACGCAGCGTTCGGAAATCAGCTCAATGGTAACGCGTGTCATTAGTCGATGGGCAACTGAGTATAATATATTTTTCAAAAGATACAGTTCCGATCAATTCGTa  >  1:217774/1‑141 (MQ=255)
                                                                                                                            ttAGATAAGTACGATGAGATTACGCAAAATATGAATGATACGCAGCGTTCGGAAATCAACTCAATGGTAACGCGTGTCATTAGTCGATGGGCAACTGAGTATAATATATTTTTCAAAAGATACAGTTCCGATCAATTCGTa  >  1:459953/1‑141 (MQ=255)
                                                                                                                                    |                                                                                                                                    
TCAGTATCATTTCCAAGTACGTTATTCTGAGAATGATCACTGCCTCTATTTCTTTGATATAACTGAACAAGTACAAACAAATGAACTATATGAAAATTCTAAACCAATCATTGCGACATTATTTTTAGATAACTACGATGAGATTACGCAAAATATGAATGATACGCAGCGTTCGGAAATCAACTCAATGGTAACGCGTGTCATTAGTCGATGGGCAACTGAGTATAATATATTTTTCAAAAGATACAGTTCCGATCAATTCGTA  >  USA300TCH1516_ALE/18758‑19022

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