Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 849,656 G→A 12.2% G261D (GGT→GAT)  USA300HOU_RS04145 → Nucleotide‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE849,6560GA12.2% 134.6 / 9.7 49G261D (GGT→GAT) USA300HOU_RS04145Nucleotide‑binding protein
Reads supporting (aligned to +/- strand):  ref base G (14/29);  new base A (3/3);  total (17/32)
Fisher's exact test for biased strand distribution p-value = 4.05e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.94e-01

ATAAAGACGTTTATAATTATGTTATGAAATGGAAAGAGACGGAGATTTTCTTTGAAAAATTAACTGATTTGTTAGATTTTATGATACCCGGGTATAAAAAAGAAGGGAAATCTCAATTAGTAATTGCCATCGGTTGTACGGGTGGACAACATCGATCTGTAGCATTAGCAGAACGACTAGGTAATTATCTAAATGAAGTATTTGAATATAATGTTTATGTGCATCATAGGGACGCACATATTGAAAGTGGCGAGAAAAAATGAGACAAATAA  >  USA300TCH1516_ALE/849524‑849795
                                                                                                                                    |                                                                                                                                           
aTAAAGACGTTTATAATTATGTTATGAAATGGAAAGAGACGGAGATTTTCTTTGAAAAATTAACTGATTTGTTAGATTTTATGATACCCGGGTATAAAAAAGAAGGGAAATCTCAATTAGTAATTGCCATCGATTGTACgg                                                                                                                                     >  2:3246/1‑141 (MQ=255)
aTAAAGACGTTTATAATTATGTTATGAAATGGAAAGAGACGGAGATTTTCTTTGAAAAATTAACTGATTTGTTAGATTTTATGATACCCGGGTATAAAAAAGAAGGGAAATCTCAATTAGTAATTGCCATCGATTGTACgg                                                                                                                                     <  1:220856/141‑1 (MQ=255)
 tAAAGACGTTTATAATTATGTTATGAAATGGAAAGAGACGGAGATTTTCTTTGAAAAATTAACTGATTTGTTAGATTTTATGATACCCGGGTATAAAAAAGAAGGGAAATCTCAATTAGTAATTGCCATCGGTTGTACggg                                                                                                                                    <  1:336646/141‑1 (MQ=255)
      aagTTTATAATTATGTTATGAAATGGAAAGAGACGGAGATTTTCTTTGAAAAATTAACTGATTTGTTAGATTTTATGATACCCGGGTATAAAAAAGAAGGGAAATCTCAATTAGTAATTGCCATCGGTTGTACGGGTGGac                                                                                                                               <  2:50177/139‑1 (MQ=255)
      aCGTTTATAATTATGTTATGAAATGGAAAGAGACGGAGATTTTCTTTGAAAAATTAACTGATTTGTTAGATTTTATGATACCCGGGTATAAAAAAGAAGGGAAATCTCAATTAGTAATTGCCATCGGTTGTACGGGTGGac                                                                                                                               <  1:595313/141‑1 (MQ=255)
        gTTTATAATTATGTTATGAAATGGAAAGAGACGGAGATTTTCTTTGAAAAATTAACTGATTTGTTAGATTTTATGATACCCGGGTATAAAAAAGAAGGGAAATCTCAATTAGTAATTGCCATCGGTTGTACGGGTGGacaa                                                                                                                             >  2:130584/1‑141 (MQ=255)
          ttATAATTATGTTATGAAATGGAAAGAGACGGAGATTTTCTTTGAAAAATTAACTGATTTGTTAGATTTTATGATACCCGGGTATAAAAAAGAAGGGAAATCTCAATTAGTAATTGCCATCGGTTGTACGGGTGGacaaca                                                                                                                           <  1:514970/141‑1 (MQ=255)
            ataATTATGTTATGAAATGGAAAGAGACGGAGATTTTCTTTGAAAAATTAACTGATTTGTTAGATTTTATGATACCCGGGTATAAAAAAGAAGGGAAATCTCAATTAGTAATTGCCATCGGTTGTACGGGTGGACAACATc                                                                                                                         <  1:16546/141‑1 (MQ=255)
             taATTATGTTATGAAATGGAAAGAGACGGAGATTTTCTTTGAAAAATTAACTGATTTGTTAGATTTTATGATACCCGGGTATAAAAAAGAAGGGAAATCTCAATTAGTAATTGCCATCGGTTGTACGGGTGGACAACATCg                                                                                                                        >  2:182959/1‑141 (MQ=255)
             taATTATGTTATGAAATGGAAAGAGACGGAGATTTTCTTTGAAAAATTAACTGATTTGTTAGATTTTATGATACCCGGGTATAAAAAAGAAGGGAAATCTCAATTAGTAATTGCCATCGATTGTACGGGTGGACAACATCg                                                                                                                        >  1:466146/1‑141 (MQ=255)
              aaTTATGTTATGAAATGGAAAGAGACGGAGATTTTCTTGGAAAAATTAACTGATTTGTTAGATTTTATGATACCCGGGTATAAAAAAGAAGGGAAATCTCAATTAGTAATTGCCATCGGTTGTACGGGTGGACAACATCGa                                                                                                                       >  1:417383/1‑141 (MQ=255)
                 tatgttatgAAATGGAAAGAGACGGAGATTTTCTTTGAAAAATTAACTGATTTGTTAGATTTTATGATACCCGGGTATAAAAAAGAAGGGAAATCTCAATTAGTAATTGCCATCGGTTGTACGGGTGGACAACATCGATCt                                                                                                                    <  1:532266/141‑1 (MQ=255)
                  atgttatgAAATGGAAAGAGACGGAGATTTTCTTTGAAAAATTAACTGATTTGTTAGATTTTATGATACCCGGGTATAAAAAAGAAGGGAAATCTCAATTAGTAATTGCCATCGGTTGTACGGGTGGACAACATCGATCTg                                                                                                                   >  1:215593/1‑141 (MQ=255)
                   tgttatgAAATGGAAAGAGACTGAGATTTTCTTTGAAAAATTAACTGATTTGTTAGATTTTATGATACCCGGGTATAAAAAAGAAGGGAAATCTCAATTAGTAATTGCCATCGGTTGTACGGGTGGACAACATCGATCTGt                                                                                                                  <  1:422408/141‑1 (MQ=255)
                            aTGGAAAGAGACGGAGATTTTCTTTGAAAAATTAACTGATTTGTTAGATTTTATGATACCCGGGTATAAAAAAGAAGGGAAATCTCAATTAGTAATTGCCATCGGTTGTACGGGTGGACAACATCGATCTGTAGCATTAGc                                                                                                         <  2:113644/141‑1 (MQ=255)
                            aTGGAAAGAGACGGAGATTTTCTTTGAAAAATTAACTGATTTGTTAGATTTTATGATACCCGGGTATAAAAAAGAAGGGAAATCTCAATTAGTAATTGCCATCGGTTGTACGGGTGGACAACATCGATCTGTAGCATTAGc                                                                                                         <  2:305400/141‑1 (MQ=255)
                                aaaGAGACGGAGATTTTCTTTGAAAAATTAACTGATTTGTTAGATTTTATGATACCCGGGTATAAAAAAGAAGGGAAATCTCAATTAGTAATTGCCATCGGTTGTACGGGTGGACAACATCGATCTGTAGCATTAGCAGaa                                                                                                     <  2:197930/141‑1 (MQ=255)
                                  agagACGGAGATTTTCTTTGAAAAATTAACTGATTTGTTAGATTTTATGATACCCGGGTATAAAAAAGAAGGGAAATCTCAATTAGTAATTGCCATCGGTTGTACGGGTGGACAACATCGATCTGTAGCATTAGCAGAACg                                                                                                   <  1:217440/141‑1 (MQ=255)
                                   gagaCGGAGATTTTCTTTGAAAAATTAACTGATTTGTTAGATTTTATGATACCCGGGTATAAAAAAGAAGGGAAATCTCAATTAGTAATTGCCATCGATTGTACGGGTGGACAACATCGATCTGTAGCATTAGCAGAACGa                                                                                                  >  1:13165/1‑141 (MQ=255)
                                   gagaCGGAGATTTTCTTTGAAAAATTAACTGATTTGTTAGATTTTATGATACCCGGGTATAAAAAAGAAGGGAAATCTCAACTAGTAATTGCCATCGGGTGTACGGGTGGACAACATCGATCTGTAGCATTAGTAGAACGa                                                                                                  >  1:371788/1‑141 (MQ=255)
                                     gaCGGAGATTTTCTTTGAAAAATTAACTGATTTGTTAGATTTTATGATACCCGGGTATAAAAAAGAAGGGAAATCTCAATTAGTAATTGCCATCGGTTGTACGGGTGGACAACATCGATCTGTAGCATTAGCAGAACGACt                                                                                                >  1:177867/1‑141 (MQ=255)
                                        ggAGATTTTCTTTGAAAAATTAACTGATTTGTTAGATTTTATGATACCCGGGTATAAAAAAGAAGGGAAATCTCAATTAGTAATTGCCATCGGTTGTACGGGTGGACAACATCGATCTGTAGCATTAGCAGAACGACTAgg                                                                                             >  1:187622/1‑141 (MQ=255)
                                           gaTTTTCTTTGAAAAATTAACTGATTTGTTAGATTTTATGATACCCGGGTATAAAAAAGAAGGGAAATCTCAATTAGTAATTGCCATCGGTTGTACGGGTGGACAACATCGATCTGTAGCATTAGCAGAACGACTAGGTaa                                                                                          <  1:541256/141‑1 (MQ=255)
                                            aTTTTCTTTGAAAAATTAACTGATTTGTTAGATTTTATGATACCCGGGTATAAAAAAGAAGGGAAATCTCAATTAGTAATTGCCATCGGTTGTACGGGTGGACAACATCGATCTGTAGCATTAGCAGAACGACTAGGTAAt                                                                                         <  1:331260/141‑1 (MQ=255)
                                              tttCTTTGAAAAATTAACTGATTTGTTAGATTTTATGATACCCGGGTATAAAAAAGAAGGGAAATCTCAATTAGTAATTGCCATCGGTTGTACGGGTGGACAACATCGATCTGTAGCATTAGCAGAACGACTAGGTAATTa                                                                                       <  1:549338/141‑1 (MQ=255)
                                              tttCTTTGAAAAATTAACTGATTTGTTAGATTTTATGATACCCGGGTATAAAAAAGAAGGGAAATCTCAATTAGTAATTGCCATCGGTTGTACGGGTGGACAACATCGATCTGTAGCATTAGCAGAACGACTAGGTAATTa                                                                                       <  1:452259/141‑1 (MQ=255)
                                               ttCTTTGAAAAATTAACTGATTTGTTAGATTTTATGATACCCGGGTATAAAAAAGAAGGGAAATCTCAATTAGTAATTGCCATCGGTTGTACGGGTGGACAACATCGATCTGTAGCATTAGCAGAACGACTAGGTAATTAt                                                                                      <  2:3233/141‑1 (MQ=255)
                                               cttcatGAAAAATTAACTGATTTGTTAGATTTTATGATACCCGGGTATAAAAAAGAAGGGAAATCTCAATTAGTAATTGCCATCGGTTGTACGGGTGGACAACATCGATCTGTAGCATTAGCAGAACGACTAGGTAATTAt                                                                                      <  2:57797/136‑1 (MQ=255)
                                                 cTTTGAAAAATTAACTGATTTGTTAGATTTTATGATACCCGGGTATAAAAAAGAAGGGAAATCTCAATTAGTAATTGCCATCGGTTGTACGGGTGGACAACATCGATCTGTAGCATTAGCAGAACGACTAGGTAATTATCt                                                                                    >  1:573631/1‑141 (MQ=255)
                                                         aaTTAACTGATTTGTTAGATTTTATGATACCCGGGTATAAAAAAGAAGGGAAATCTCAATTAGTAATTGCCATCGGTTGTACGGGTGGACAACATCGATCTGTAGCATTAGCAGAACGACTACGTAATTATCTAAATGAAg                                                                            >  1:483355/1‑141 (MQ=255)
                                                           ttAACTGATTTGTTAGATTTTATGATACCCGGGTATAAAAAAGAAGGGAAATCTCAATTAGTAATTGCCATCGGTTGTACGGGTGGACAACATCGATCTGTAGCATTAGCAGAACGACTAGGTAATTATCTAAATGAAGTa                                                                          <  1:390350/141‑1 (MQ=255)
                                                                 gATTTGTTAGATTTTATGATCCCCGGGTATAAAAAAGAAGGTAAATCTCAATTAGTAATTGCCATCGATTGTACGGGTGGACAACATCGATCTGTAGCATTAGCAGAACGACTAGGTAATTATCTAAATGAAGTATTTGaa                                                                    <  1:10610/141‑1 (MQ=255)
                                                                    ttGTTAGATTTTATGATACCCGGGTATAAAAAAGAAGGGAAATCTCAATTAGTAATTGCCATCGGTTGTACGGGTGGACAACATCGATCTGTAGCATTAGCAGAACGACTAGGTAATTATCTAAATGAAGTATTTGAatat                                                                 <  1:329631/141‑1 (MQ=255)
                                                                      gTTAGATTTTATGATACCCGGGTATAAAAAAGAAGGGAAATCTCAATTAGTAATTGCCATCGGTTGTACGGGTGGACAACATCGATCTGTAGCATTAGCAGAACAACTAGGTAATTATCTAAATGAAGTATTTGAATATaa                                                               <  1:184488/141‑1 (MQ=255)
                                                                      gTTAGATTTTATGATACCCGGGTATAAAAAAGAAGGGAAATCTCAATTAGTAATTGCCATCGATTGTACGGGTGGACAACATCGATCTGTAGCATTAGCAGAACGACTAGGTAATTATCTAAATGAAGTATTTGAATATaa                                                               <  1:191309/141‑1 (MQ=255)
                                                                              ttATGATACCCGGGTATAAAAAAGAAGGGAAATCTCAATTAGTAATTGCCATCGGTTGTACGGGTGGACAACATCGATCTGTAGCATTAGCAGAACGACTAGGTAATTATCTAAATGAAGTATTTGAATATAATGTTTAtg                                                       <  1:230190/141‑1 (MQ=255)
                                                                                aTGATACCCGGGTATAAAAAAGAAGGGAAATCTCAATTAGTAATTGCCATCGGTTGTACGGGTGGACAACATCGATCTGTAGCATTAGCAGAACGACTCGGTAATTATCTAAATGAAGTATTTGAATATAATGTTTAtgtg                                                     >  1:158583/1‑141 (MQ=255)
                                                                                aTGATACCCGGGTATAAAAAAGAAGGGAAATCTCAATTAGTAATTGCCATCGGTTGTACGGGTGGACAACATCGATCTGTAGCATTAGCAGAACGACTAGGTAATTATCTAAATGAAGTATTTGAATATAATGTTTAtgtg                                                     >  1:782590/1‑141 (MQ=255)
                                                                                  gATACCCGGGTATAAAAAAGAAGGGAAGTCTCAATTAGTAATTGCCATCGGTTGTACGGGTGGACAACATCGATCTGTAGCATTAGCAGAACGACTAGGTAATTATCTAAATGAAGTATTTGAATATAATGTTTATGTGca                                                   <  1:420672/141‑1 (MQ=255)
                                                                                       ccGGGTATAAAAAAGAAGGGAAATCTCAATTAGTAATTGCCATCGGTTGTACGGGTGGACAACATCGATCTGTAGCATTAGCAGAACGACTAGGTAATTATCTAAATGAAGTATTTGAATATAATGTTTATGTGCATCATa                                              >  1:682413/1‑141 (MQ=255)
                                                                                               aaaaaaGAAGGGAAATCTCAATTAGTAATTGCCATCGGTTGTACGGGTGGACAACATCGATCTGTAGCATTAGCAGAACGACTAGGTAATTATCTAAATGAAGTATTTGAATATAATGTTTATGTGCATCATAGGGACGca                                      <  1:118630/141‑1 (MQ=255)
                                                                                                   aagaagGGAAATCTCAATTAGTAATTGCCATCGGTTGTACGGGTGGACAACATCGATCTGTAGCATTAGCAGAACGACTAGGTAATTATCTAAATGAAGTATTTGAATATAATGTTTATGTGCATCATAGGGACGCACata                                  <  1:55638/141‑1 (MQ=255)
                                                                                                    agaagGGAAATCTCAATTAGTAATTGCCATCGGTTGTACGGGTGGACAACATCGATCTGTAGCATTAGCAGAACGACTAGGTAATTATCTAAATGAAGTATTTGAATATAATGTTTATGTGCATCATAGGGTCGCACatat                                 <  1:692007/141‑1 (MQ=255)
                                                                                                     gaagGGAAATCTCAATTAGTAATTGCCATCGGTTGTACGGGTGGACAACATCGATCTGTAGCATTAGCAGAACGACTAGGTAATTATCTAAATGAAGTATTTGAATATAATGTTTATGTGCATCATAGGGACGCACATAtt                                <  1:3458/141‑1 (MQ=255)
                                                                                                       agGGAAATCTCAATTAGTAATTGCCATCGGTTGTACGGGTGGACAACATCGATCTGTAGCATTAGCAGAACGACTAGGTAATTATCTAAATGAAGTATTTGAATATAATGTTTATGTGCATCATAGGGACGCACATATTGa                              <  1:767736/141‑1 (MQ=255)
                                                                                                            aaTCTCAATTAGTAATTGCCATCGGTTGTACGGGTGGACAACATCGATCTGTAGCATTAGCAGAACGACTAGGTAATTATCTAAATGAAGTATTTGAATATAATGTTTATGTGCATCATAGGGACGCACATATTGAAAGTg                         >  1:361664/1‑141 (MQ=255)
                                                                                                               ctcAATTAGTAATTGCCATCGGTTGTACGGGTGGACAACATCGATCTGTAGCATTAGCAGAACGACTAGGTAATTATCTAAATGAAGTATTTGAATATAATGTTTATGTGCATCATAGGGACGCACATATTGAAAGTGGCg                      <  2:279506/141‑1 (MQ=255)
                                                                                                                             gCCATCGGTTGTACGGGTGGACAACATCGATCTGTAGCATTAGCAGAACGACTAGGTAATTATCTAAATGAAGTATTTGAATATAATGTTTATGTGCATCATAGGGACGCACATATTGAAAGTGGCGAGAAAAAATGAGAc        <  1:113430/141‑1 (MQ=255)
                                                                                                                                aTCGGTTGTACGGGTGGACAACATCGATCTGTAGCATTAGCAGAACGACTAGGTAATTATCTAAATGAAGTATTTGAATATAATGTTTATGTGCATCATAGGGACGCACATATTGAAAGTGGCGAGAAAAAATGAGACaaa     >  1:446662/1‑141 (MQ=255)
                                                                                                                                   ggTTGTACGGGTGGACAACATCGATCTGTAGCATTAGCAGAACGACTAGGTAATTATCTAAATGAAGTATTTGAATATAATGTTTATGTGCATCATAGGGACGCACATATTGAAAGTGGCGAGAAAAAATGAGACAAATaa  <  2:242772/141‑1 (MQ=255)
                                                                                                                                    |                                                                                                                                           
ATAAAGACGTTTATAATTATGTTATGAAATGGAAAGAGACGGAGATTTTCTTTGAAAAATTAACTGATTTGTTAGATTTTATGATACCCGGGTATAAAAAAGAAGGGAAATCTCAATTAGTAATTGCCATCGGTTGTACGGGTGGACAACATCGATCTGTAGCATTAGCAGAACGACTAGGTAATTATCTAAATGAAGTATTTGAATATAATGTTTATGTGCATCATAGGGACGCACATATTGAAAGTGGCGAGAAAAAATGAGACAAATAA  >  USA300TCH1516_ALE/849524‑849795

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