Predicted mutation
evidence seq id position mutation annotation gene description
RA USA300TCH1516_ALE 849,203 A→G Y110C (TAT→TGT)  USA300HOU_RS04145 → Nucleotide‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE849,2030AG90.0% 73.5 / 0.4 30Y110C (TAT→TGT) USA300HOU_RS04145Nucleotide‑binding protein
Reads supporting (aligned to +/- strand):  ref base A (0/3);  new base G (19/8);  total (19/11)
Fisher's exact test for biased strand distribution p-value = 4.06e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.66e-01

GCAATTGCAATTGATTTAAGAGGTAAGGAACTATTTAATTCATTAGTTGCAGTAGTGGATAAAGTCAAAAGTGAAAGTGACGTCATCATTGATGTTATGTTTTTAGAAGCAAGTACTGAAAAATTAATTTCAAGATATAAGGAAACGCGTCGTGCACATCCTTTGATGGAACAAGGTAAAAGATCGTTAATCAATGCAATTAATGATGAGCGAGAGCATTTGTCTCAAATTAGAAGTATAGCTAATTTTGTTAT  >  USA300TCH1516_ALE/849067‑849320
                                                                                                                                        |                                                                                                                     
gcaattgcaattGATTTAAGAGGAAAGGAACTATTTAATTCATTAGTTGCAGTAGTGGATAAAGTCAAAAGTGAAAGTGACGTCATCATTGATGTTATGTTTTTAGAAGCAAGTACTGAAAAATTAATTTCAAGATGTAAg                                                                                                                   >  2:393217/1‑141 (MQ=255)
           tgATTTAAGAGGTAAGGAACTATTTAATTCATTAGTTGCAGTAGTGGATAAAGTCAAAAGTGAAAGTGACGTCATCATTGATGTTATGTTTTTAGAAGCAAGTACTGAAAAATTAATTTCAAGATGTAAGGAAACGcgtcg                                                                                                        >  2:273490/1‑141 (MQ=255)
              tttAAGAGGTAAGGAACTATTTAATTCATTAGTTGCAGTAGTGGATAAAGTCAAAAGTGAAAGTGACGTCATCATTGATGTTATGTTTTTAGAAGCAAGTACTGAAAAATTAATTTCAAGATGTAAGGAAACGCGTCGTGc                                                                                                     <  1:409065/141‑1 (MQ=255)
                tAAGAGGTAAGGAACTATTTAATTCATTAGTTGCAGTAGTGGATAAAGTCAAAAGTGAAAGTGACGTCATCATTGATGTTATGTTTTTAGAAGCAAGTACTGAAAAATTAATTTCAAGATGTAAGGAAACGCGTCGTGcac                                                                                                   <  2:137984/141‑1 (MQ=255)
                  agagGTAAGGAACTATTTAATTCATTAGTTGCAGTAGTGGATAAAGTCAAAAGTGAAAGTGACGTCATCATTGATGTTATGTTTTTAGAAGCAAGTACTGAAAAATTAATTTCAAGATGTAAGGAAACGCGTCGTGCACAt                                                                                                 >  1:208362/1‑141 (MQ=255)
                  agagGTAAGGAACTATTTAATTCATTAGTTGCAGTAGTGGATAAAGTCAAAAGTGAAAGTGACGTCATCATTGATGTTATGTTTTTAGAAGCAAGTACTGAAAAATTAATTTCAAGATGTAAGGAAACGCGTCGTGCACAt                                                                                                 >  1:323382/1‑141 (MQ=255)
                     ggTAAGGAACTATTTAATTCATTAGTTGCAGTAGTGGATAAAGTCAAAAGTGAAAGTGACGTCATCATTGATGTTATGTTTTTAGAAGCAAGTACTGAAAAATTAATTTCAAGATGTAAGGAAACGCGTCGTGCACATCCt                                                                                              <  1:321747/141‑1 (MQ=255)
                               tATTTAATTCATTAGTTGCAGTAGTGGATAAAGTTAAAAGTGAAAGTGACGTCATCATTGATGTTATGTTTTTAGAAGCAAGTACTGAAAAATTAATTTCAAGATATAAAGAACCGCGTCGTGCACATCCTTTGATGGAAc                                                                                    <  1:373056/141‑1 (MQ=255)
                               tATTTAATTCATTAGTTGCAGTAGTGGATAAAGTTAAAAGTGAAAGTGACGTCATCATTGATGTTATGTTTTTAGAAGCAAGTACTGAAAAATTAATTTCAAGATATAAAGAAACGCGTCGTGCACATCCTTTGATGGAAc                                                                                    <  2:32831/141‑1 (MQ=255)
                               tATTTAATTCATTAGTTGCAGTAGTGGATAAAGTTAAAAGTGAAAGTGACGTCATCATTGATGTTATGTTTTTAGAAGCAAGTACTGAAAAATTAATTTCAAGATATAAAGAAACGCGTCGTGCACATCCTTTGATGGAAc                                                                                    <  1:373059/141‑1 (MQ=255)
                                 tttAATTCATTAGTTGCAGTAGTGGATAAAGTCAAAAGTGAAAGTGACGTCATCATTGATGTTATGTTTTTAGAAGCAAGTACTGAAAAATTAATTTCAAGATGTAAGGAAACGCGTCGTGCACATCCTTTGATGGAACaa                                                                                  >  2:193619/1‑141 (MQ=255)
                                    aaTTCATTAGTTGCAGTAGTGGATAAAGTCAAAAGTGAAAGTGACGTCATCATTGATGTTATGTTTTTAGAAGCAAGTACTGAAAAATTAATTTCAAGATGTAAGGAAACGCGTCGTGCACATCCTTTGATGGAACAAGGt                                                                               <  2:297434/141‑1 (MQ=255)
                                     aTTCATTAGTTGCAGTAGTGGATAAAGTCAAAAGTGAAAGTGACGTCATCATTGATGTTATGTTTTTAGAAGCAAGTACTGAAAAATTAATTTCAAGATGTAAGGAAACGCGTCGTGCACATCCTTTGATGGAACAAGGTa                                                                              >  1:5254/1‑141 (MQ=255)
                                     aTTCATTAGTTGCAGTAGTGGATAAAGTCAAAAGTGAAAGTGACGTCATCATTGATGTTATGTTTTTAGAAGCAAGTACTGAAAAATTAATTTCAAGATGTAAGGAAACGCGTCGTGCACATCCTTTGATGGAACAAGGTa                                                                              >  1:98777/1‑141 (MQ=255)
                                     aTTCATTAGTTGCAGTAGTGGATAAAGTCAAAAGTGAAAGTGACGTCATCATTGATGTTATGTTTTTAGAAGCAAGTACTGAAAAATTAATTTCAAGATGTAAGGAAACGCGTCGTGCACATCCTTTGATGGAACAAGGTa                                                                              >  2:103356/1‑141 (MQ=255)
                                     aTTCATTAGTTGCAGTAGTGGATAAAGTCAAAAGTGAAAGTGACGTCATCATTGATGTTATGTTTTTAGAAGCAAGTACTGAAAAATTAATTTCAAGATGTAAGGAAACGCGTCGTGCACATCCTTTGATGGAACAAGGTa                                                                              >  2:272650/1‑141 (MQ=255)
                                               tGCAGTAGTGGATAAAGTCAAAAGTGAAAGTGACGTCATCATTGATGTTATGTTTTTAGAAGCAAGTACTGAAAAATTAATTTCAAGATGTAAGGAAACGCGTCGTGCACATCCTTTGATGGAACAAGGTAAAAGATCGtt                                                                    >  2:157764/1‑141 (MQ=255)
                                                 cAGTAGTGGATAAAGTCAAAAGTGAAAGTGACGTCATCATTGATGTTATGTTTTTAGAAGCAAGTACTGAAAAATTAATTTCAAGATGTAAGGAAACGCGTCGTGCACATCCTTTGATGGAACAAGGTAAAAGATCGTTaa                                                                  <  1:193560/141‑1 (MQ=255)
                                                               gTCAAAAGTGAAAGTGACGTCATCATTGATGTTATGTTTTTAGAAGCAAGTACTGAAAAATTAATTTCAAGATGTAAGGAAACGCGTCGTGCACATCCTTTGATGGAACAAGGTAAAAGATCGTTAATCAATGCAATTAat                                                    <  2:393529/141‑1 (MQ=255)
                                                                  aaaaGTGAAAGTGACGTCATCATTGATGTTATGTTTTTAGAAGCAAGTACTGAAAAATTAATTTCAAGATGTAAGGAAACGCGTCGTGCACATCCTTTGATGGAACAAGGTAAAAGATCGTTAATCAATGCAATTAatgat                                                 >  2:25031/1‑141 (MQ=255)
                                                                           agtgaCGTCATCATTGATGTTATGTTTTTAGAAGCAAGTACTGAAAAATTAATTTCAAGATGTAAGGAAACGCGTCGTGCACATCCTTTGATGGAACAAGGTAAAAGATCGTTAATCAATGCAATTAATGATGAGCGagag                                        >  1:43514/1‑141 (MQ=255)
                                                                           agtgaCGTCATCATTGATGTTATGTTTTTAGAAGCAAGTACTGAAAAATTAATTTCAAGATGTAAGGAAACGCGTCGTGCACATCCTTTGATGGAACAAGGTAAAAGATCGTTAATCAATGCAATTAATGATGAGCGagag                                        >  1:150812/1‑141 (MQ=255)
                                                                                        ttGATGTTATGTTTTTAGAAGCAAGTACTGAAAAATTAATTTCAAGATGTAAGGAAACGCGTCGTGCACATCCTTTGATGGAACAAGGTAAAAGATCGTTAATCAATGCAATTAATGATGAGCGAGAGCATTTGTCTCaaa                           <  2:97594/141‑1 (MQ=255)
                                                                                                           aGCAAGTACTGAAAAATTGATTTCAAGATGTAAGGAAACGCGTCGTGCACATCCTTTGATGGAACAAGGTAAAAGATCGTTAATCAATGCAATTAATGATGAGCGAGAGCATTTGTATCAAATTAGAAGTATAGCTAAttt        <  1:256762/141‑1 (MQ=255)
                                                                                                            gCAAGTACTGAAAAATTAATTTCAAGATGTAAGGAAACGCGTCGTGCACATCCTTTGATGGAACAAGGTAAAAGATCGTTAATCAATGCAATTAATGATGAGCGAGAGCATTTGTCTCAAATTAGAAGTATAGCTAAtttt       >  2:159761/1‑141 (MQ=255)
                                                                                                            gCAAGTACTGAAAAATTAATTTCAAGATGTAAGGAAACGCGTCGTGCACATCCTTTGATGGAACAAGGTAAAAGATCGTTAATCAATGCAATTAATGATGAGCGAGAGCATTTGTCTCAAATTAGAAGTATAGCTAAtttt       >  2:126248/1‑141 (MQ=255)
                                                                                                            gCAAGTACTGAAAAATTAATTTCAAGATGTAAGGAAACGCGTCGTGCACATCCTTTGATGGAACAAGGTAAAAGATCGTTAATCAATGCAATTAATGATGAGCGAGAGCATTTGTCTCAAATTAGAAGTATAGCTAAtatt       >  2:52460/1‑141 (MQ=255)
                                                                                                             cAAGTACTGAAAAATTAATTTCAAGATGTAAGGAAACGCGTCGTGCACATCCTTTGATGGAACAAGGTAAAAGATCGTTAATCAATGCAATTAATGATGAGCGAGAGCATTTGTCTCAAATTAGAAGTATAGCTAATTTTg      >  1:129711/1‑141 (MQ=255)
                                                                                                             cAAGTACTGAAAAATTAATTTCAAGATGTAAGGAAACGCGTCGTGCACATCCTTTGATGGAACAAGGTAAAAGATCGTTAATCAATGCAATTAATGATGAACGAGAGCATTTGTCTCAAATTAGAAGTATAGCTAATTTTg      >  1:184605/1‑141 (MQ=255)
                                                                                                                 tACTGAAAAATTAATTTCAAGATGTAAGGAAACGCGTCGTGCACATCCTTTGATGGAACAAGGTAAAAGATCGTTAATCAATGCAATTAATGATGAGCGAGAGCATTTGTCTCAAATTAGAAGTATAGCTAATTTTGTtat  >  2:124801/1‑141 (MQ=255)
                                                                                                                                        |                                                                                                                     
GCAATTGCAATTGATTTAAGAGGTAAGGAACTATTTAATTCATTAGTTGCAGTAGTGGATAAAGTCAAAAGTGAAAGTGACGTCATCATTGATGTTATGTTTTTAGAAGCAAGTACTGAAAAATTAATTTCAAGATATAAGGAAACGCGTCGTGCACATCCTTTGATGGAACAAGGTAAAAGATCGTTAATCAATGCAATTAATGATGAGCGAGAGCATTTGTCTCAAATTAGAAGTATAGCTAATTTTGTTAT  >  USA300TCH1516_ALE/849067‑849320

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: ‑