Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 2,541,351 A→G 9.5% L307L (TTA→TTG)  sbi → Immunoglobulin‑binding protein sbi

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE2,541,3510AG 9.5% 104.0 / 4.7 43L307L (TTA→TTG) sbiImmunoglobulin‑binding protein sbi
Reads supporting (aligned to +/- strand):  ref base A (21/17);  new base G (2/2);  total (24/19)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.86e-01

AGTGGCGCCAAAAGTTGAGGCTCCTCAAATTCAATCACCACAAATTGAAAAACCTAAAGTAGAATCACCAAAAGTTGAAGTCCCTCAAATTCAATCACCAAAAGTTGAGGTTCCTCAATCTAAATTATTAGGTTACTACCAATCATTAAAAGATTCATTTAACTATGGTTACAAGTATTTAACAGATACTTATAAAAGCTATAAAGAAAAATATGATACAGCAAAGTACTACTATAATACGTACTATAAATACAAAGGTGCGATTGA  >  USA300TCH1516_ALE/2541222‑2541488
                                                                                                                                 |                                                                                                                                         
agTGGCGCCAAAAGTTGAGGCTCCTCAAATTCTATCACCACAAATTGAAAAACCTAAAGTAGAATCACCAAAATTTGAAGTCCCTCAAATTCAATCACCAAAAGTTGAGGTTCCTCAATCTAAATTATTAGGTTACTACCa                                                                                                                                >  2:371826/1‑141 (MQ=255)
 gTGGCGCCAAAAGTTGCGGCTCCTCAAATTCAATCACCACAAATTGAAAAACCTAAAGTAGAATCACCAAAAGATGAAGTCCCTCAAATTCAATCACCAAAAATTGAAGTTCCTCAATCTAAATTATTAGCTTACTCCCaa                                                                                                                               >  2:81606/1‑141 (MQ=255)
 gTGGCGCCAAAAGTTGAGGCTCCTCAAATTCAATCACCACAAATTGAAAAACCTAAAGTAAAATCACCAAAAGTTGAAGTCCCTCAAATTCAATCACCAAAAGTTGAGGTTCCTCACTCTAAATAATTAGGTTACTACCaa                                                                                                                               >  2:37264/1‑141 (MQ=255)
   ggCGCCAAAAGTTGAGGCTCCTCAAATTCAATCACCACAAATTGAAAAACCTAAAGTAGAATCACCAAAAGTTGAAGTCCCTCAAATTCAATCACCAAAAGTTGCGGTTCCTCAATCTAAATTATTAGGTTACTACCAATc                                                                                                                             >  1:187675/1‑141 (MQ=255)
      gcCAAAAGTTGAGGCTCCTCAAATTCAATCACCACAAATTGAAAAACCTAAAGTAGAATCACCAAAAGTTGAAGTCCCTCAAATTCAATCACCAAAAGTTGAGGTTCCTCAATCTAAATTATTAGGTTACTACCAATCAtt                                                                                                                          <  2:24125/141‑1 (MQ=255)
      gcCAAAAGTTGAGGCTCCTCAAATTCAATAACCACAAATTGAAAAACCTAAAGTAGAATCACCAAAAGTTGAAGTCCCTCAAATTCAATCACCAAAAGTTGAGGTTCTCCATCCCAATTTACTATGTTACTACCAATCAtg                                                                                                                          >  2:385695/1‑140 (MQ=255)
             gTTGAGGCTCCTCAAATTCAATCACCACAAATTGAAAAACCTAAAGTAGAATCACCAACAGTTGAAGTCCCTCAAATTCAATCACCAAAAGTTGAGGTTCCTCAATCTAAATTATTAGGTTACTACCAATCATTAAAAGAt                                                                                                                   <  1:322000/141‑1 (MQ=255)
               tGAGGCTCCTCAAATTCAATCACCACAAATTGAAAAACCTAAAGTAGAATCACCAAAAGTTGAAGTCCCTCAAATTCAATCACCAAAAGTTGAGGTTCCTCAATCTAAATTATTCGGTTACTACCAATCATTAACAGATTc                                                                                                                 >  1:20035/1‑141 (MQ=255)
                 aGGCTCCTCAAATTCAATCACCACAAATTGAAAAACCTAAAGTAGAATCACCAAAAGTTGAAGTCCCTCAAATTCAATCACCAAAAGTTGAGGTTCCTCAATCTAAATTATTAGGTTACTACCAATCATTAAAAGATTCAt                                                                                                               <  2:153491/141‑1 (MQ=255)
                 aGGCTCCTCAAATTCAATCACCACAAATTGAAAAACCTAAAGTAGAATCACCAAAAGTTGAAGTCCCTCAAATTCAATCACCAAAAGTTGAGGTTCCTCAATCTAAATTATTAGGTTACTACCAATCATTAAAAGATTCAt                                                                                                               >  1:83305/1‑141 (MQ=255)
                 aGGCTCCTCAAATTCAATCACCACAAATTGAAAAACCTAAAGTAGAATCACCAAAAGTTGAAGTCCCTCAAATTCAATCACCAAAAGTTGAGGTTCCTCAATCTAAATTATTAGGTTACTACCAATCATTAAAAGATTCAt                                                                                                               >  1:5878/1‑141 (MQ=255)
                          aaaTTCAATCACCACAAATTGAAAAACCTAAAGTAGAATCACTAAAAGTTGAAGTCCTTAAAATGCAATCACCAAAAGTTGAGGTTCCTCAACCATAATTACTAGGGTACTACCAATCATTAAAAGATTCATTTAACTATg                                                                                                      <  2:278697/141‑1 (MQ=255)
                            aTTCAATCACCACAAATTGAAAAACCTAAAGTAGAATCACCAAAAGTTGAAGTCCCTCAAATTCAATCACCAAAAGTTGAGGTTCCTCAATCTAAATTATTAGGTTACTACCAATCATTAAAAGATTCATTTAACTATGGt                                                                                                    >  1:357191/1‑141 (MQ=255)
                            aTTCAATCACCACAAATTGAAAAACCTAAAGTAGAATCACCAAAAGTTGAAGTCCCTCAAATTCAATCACCAAAAGTTGAGGTTCCTCAATCTAAATTATTAGGTTACTACCAATAATTAAAAGATTCATTTAACTATGGt                                                                                                    >  1:89384/1‑141 (MQ=255)
                              tCAATCACCACAAATTGAAAAACCTAAAGTAGAATCACCAAAAGTTGAAGTCCCTCAAATTCAATCACCAAAAGTTGAGGTTCCTCAATCTAAATTATTAGGTTACTACCAATCATTAAAAGATTCATTTAACTATGGTTa                                                                                                  >  1:5799/1‑141 (MQ=255)
                                 aTCACCACAAATTGAAAAACCTAAAGTAGAATCACCAAAAGTTGAAGTCCCTCAAATTCAATCACCAAAAGTTGAGGTTCCTCAATCTAAATTATTAGGTTACTACCAATCATTAAAAGATTCATTTAACTATGGTTACaa                                                                                               >  1:211328/1‑141 (MQ=255)
                                  tCACCACAAATTGAAAAACCTAAAGTAGAATCACCAAAAGTTGAAGTCCCTCAAATTCAATCACCAAAAGTTGAGGTTCCTCAATCTAAATTATTAGGTTACTACCAATCATTAAAAGATTCATTTAACTATGGTTACAAg                                                                                              <  1:104847/141‑1 (MQ=255)
                                       acaAATTGAAAAACCTAAAGTAGAATCACCAAAAGTTGAAGTCCCTCAAATTCAATCAACAAAAGTTGAGGTTCCTCAATCTAAATTATTAGGTTACTACCAATCATTAAAAGATTCATTTAACTATGGTTACAAGTAttt                                                                                         <  2:83305/141‑1 (MQ=255)
                                          aaTTGAAAAACCTAAAGTAGAATCACCAAAAGTTGAAGTCCCTCAAATTCAATCACCAAAAGTTGAGGTTCCTCAATCTAAATTATTAGGTTACTACCAATCATTAAAAGATTCATTTAACTATGGTTACAAGTATTTAAc                                                                                      <  1:129566/141‑1 (MQ=255)
                                            ttGAAAAACCTAAAGTAGAATCACCAAAAGTTGAAGTCCCTCAAATTCAATCACCAAAAGTTGAGGTTCCTCAATCTAAATTATTAGGTTACTACCAATCATTAAAAGATTCATTTAACTATGGTTACAAGTATTTAACAg                                                                                    >  1:175584/1‑141 (MQ=255)
                                                aaaaCCTAAAGTAGAATCACCAAAAGTTGAAGTCCCTCAAATTCAATCACCAAAAGTTGAGGTTCCTCAATCTAAATTATTAGGTTACTACCAATCATTAAAAGATTCATTTAACTATGGTTACAAGTATTTAACAGATAc                                                                                >  1:315050/1‑141 (MQ=255)
                                                   aCCTAAAGTAGAATCACCAAAAGTTGAAGTCCCTCAAATTCAATCACCAAAAGTTGAGGTTCCTCAATCTAAATTATTAGGTTACTACCAATCATTAAAAGATTCATTTAACTATGGTTACAAGTATTTAACAAATACTta                                                                             >  2:356562/1‑141 (MQ=255)
                                                      tAAAGTAGAATCACCAAAAGTTGAAGTCCCTCAAATTCAATCACCAAAAGTTGAGGTTCCTCAATCTAAATTATTAGGTTACTACCAATCATTAAAAGATTCATTTAACTATGGTTACAAGTATTTAACAGATACTTATaa                                                                          <  2:19045/141‑1 (MQ=255)
                                                       aaaGTAGAATCACCAAAAGTTGAAGTCCCTCAAATTCAATCACCAAAAGTTGAGGTTCCTCAATCTAAATTATTAGGTTACTACCAATCCTTAAACAGCGCATTTAACTATGGTTACAAGTATTTAACAGATACTTATaaa                                                                         >  2:409007/1‑141 (MQ=255)
                                                              aaaaagcAAAAGTTGAAGCTGCTGAAATACACTCACCAAAAGTTGAGGTTCCTCAATCTATATTATTAGGTTACTACCAAGCATTAAAAGATTCATTTAACTATGGTTACAAGTATTTAACAGATACTTATAAAAGCtata                                                                  <  2:293571/135‑1 (MQ=255)
                                                               aaaaacAAAAGTTGAAGGTCATCAAAATAAATCACCAAAAGTTGAGGTTCATCAATCTAAATTATTAGGTTAATAAAAATCATTAAAAGATTCATTAAACTATGGTTACAAGTATTTAACAGATACATAAGAAAGCTATaa                                                                 <  2:8188/136‑1 (MQ=255)
                                                                  aCCAAAAGTTGAAGTCCCTCAAATTCAATCACCAAAAGTTGAGGTTCCTCAATCTAAATTATTAGGTTACTACCAATCATTAAAAGATTCATTTAACTATGGTTACAAGTATTTAACAGATACTTATAAAAGCTATAAAGa                                                              <  1:401201/141‑1 (MQ=255)
                                                                      aaaGTTGAAGTCCCTCAAATTCAATCACCAAAAGTTGAGGTTCCTCAATCTAAATTATGAGGTTACTACCAATAATGAAAAGATTCATTTAACTATGGTTACAAGTATTTAACAGATACTTATAAAAGCTATAAAGaaaaa                                                          <  2:233014/141‑1 (MQ=255)
                                                                         gTTGAAGTCCCTCAAATTCAATCACCAAAAGTTGAGGTTCCTCAATCTAAATTATTAGGTTACTACCAATCATTAAAAGATTCATTTAACTATGGTTACAAGTATTTAACAGATACTTATAAAAACTATAAAGAAAAatat                                                       >  2:357919/1‑141 (MQ=255)
                                                                             aaGTCCCTCAAATTCAATCACCAAAAGTTGAGGTTCCTCAATCTAAATTATTAGGTTACTACCAATCATTAAAAGATTCATTTAACTATGGTTACAAGTATTTAACAGATACTTATAAAAGCTATAAAGAAAAATATGATa                                                   >  2:206645/1‑141 (MQ=255)
                                                                                   cTCAAATTCAATCACCAAAAGTTGAGGTTCCTCAATCTAAATTATTAGGTTACTACCAATCATTAAAAGATTCATTTAACTATGGTTACAAGTATTTAACAGATACTTATAAAAGCTATAAAGAAAAATATGATACAGCaa                                             <  1:17322/141‑1 (MQ=255)
                                                                                        aTTCAATCACCAAAAGTTGAGGTTACTCAATCTAAATTATTAGGTTACTACCAATCATTAAAAGATTCATTTATCTATGTATACAAGTATTTAACAGATACTTATAAAAGCTATAAAGAAAAATATGATACAGCAAAGtac                                        <  2:7637/141‑1 (MQ=255)
                                                                                         ttCAATCACCAAAAGTTGAGGTTCCTCAATCTAAATTATTAGGTTACTACCAATCATTAAAAGATTCATTTAACTATGGTTACAAGTATTTAACAGATACTTATAAAAGCTATAAAGAAAAATATGATACAGGAAAAtaat                                       >  2:356157/1‑139 (MQ=255)
                                                                                           cAATCACCAAAAGTTGAGGTTCCTCAATCTAAATTATTGGGTTACTACCAATCATTAAAAGATTCATTTAACTATGGTTACAAGTATTTAACAGATACt                                                                               >  2:13733/1‑99 (MQ=255)
                                                                                           cAATCACCAAAAGTTGAGGTTCCTCAATCTAAATTATTGGGTTACTACCAATCATTAAAAGATTCATTTAACTATGGTTACAAGTATTTAACAGATACt                                                                               >  2:12562/1‑99 (MQ=255)
                                                                                           cAATCACCAAAAGTTGAGGTTCCTCAATCTAAATTATTGGGTTACTACCAATCATTAAAAGATTCATTTAACTATGGTTACAAGTATTTAACAGATACt                                                                               <  1:13733/99‑1 (MQ=255)
                                                                                           cAATCACCAAAAGTTGAGGTTCCTCAATCTAAATTATTGGGTTACTACCAATCATTAAAAGATTCATTTAACTATGGTTACAAGTATTTAACAGATACt                                                                               <  1:12562/99‑1 (MQ=255)
                                                                                             aTCACCAAAAGTTGAGGTTCCTCAATCTAAATTATTAGGTTACTACCAATCATTAAAAGATTCATTTAACTATGGTTACAAGTATTTAACAGATACTTATAAAAGCTATAAAGAAAAATATGATACAGCAAAGTActacta                                   >  1:225131/1‑141 (MQ=255)
                                                                                                acaaaaaTTTGAGATTCCTTAATCTAAATTATTAGGTTACTACCAATCATTAAAAGATTCATTTAACTATGGTTACAAGTATTTAACAGATACTTATAAAAGCTATAAAGAAAAATATGATACAGCAAAGTACTACTataa                                <  2:22458/138‑1 (MQ=255)
                                                                                                aCCAAAAGTTGAGGTTCCTCAATCTAAATTATTAGGTTACTACCAATCATTAAAAGATTCATTTAACTATGGTTACAAGTATTTAACAGATACTTATAAAAGCTATAAAGAAAAATATGATACAGCAAAGTACTACTataa                                <  1:110766/141‑1 (MQ=255)
                                                                                                       gTTGAGGTTCCTCAATCTAAATTATTAGGTTACTACCAATCATTAAAAGATTCATTTAACTATGGTTACAAGTATTTAACAGATACTTATAAAAGCTATAAAGAAAAATATGATACAGCAAAGTACTAATATAATACGTAc                         >  1:240605/1‑141 (MQ=255)
                                                                                                            ggTTCCTCAATCTAAATTATTAGGTTACTACCAATCATTAAAAGATTCATTTAACTATGGTTACAAGTATTTAACAGATACTTATAAAAGCTATAAAGAAAAATATGATACAGCAAAGTACTACTATAATACGTACTATaa                    >  1:26569/1‑141 (MQ=255)
                                                                                                                   aaaTCTAAATTATTAGGTTACTACCAATCATTAAAAGATTCATTTAACTATGGTTACAAGTATTTAACAGATACTTATAAAAGCTATAAAGAAAAATATGATACAGCAAAGTACTACTATAATACGTACTATAAATACaaa             <  2:219438/140‑1 (MQ=255)
                                                                                                                            ttattaGGTTACTACCAATCATTAAAAGATTCATTTAAATTTGGTAACGAGTATTTAACAGATACTTATAAAAGCTATAAAGAAAAATATGATACAGCAAAGTACTACTATAATACGTACTATAAATACAAAGGTGCGAtt    <  1:398413/141‑1 (MQ=255)
                                                                                                                              attaGGTTACTACCAATCATTAAAAGATTCATTTAACTATGGTTACAAGTATTTAACAGATACTTATAAAAGCTATAAAGAAAAATATGATACAGCAAAGTACTACTATAATACGTACTATAAATACAAAGGTGCGATTGa  <  1:206645/141‑1 (MQ=255)
                                                                                                                                 |                                                                                                                                         
AGTGGCGCCAAAAGTTGAGGCTCCTCAAATTCAATCACCACAAATTGAAAAACCTAAAGTAGAATCACCAAAAGTTGAAGTCCCTCAAATTCAATCACCAAAAGTTGAGGTTCCTCAATCTAAATTATTAGGTTACTACCAATCATTAAAAGATTCATTTAACTATGGTTACAAGTATTTAACAGATACTTATAAAAGCTATAAAGAAAAATATGATACAGCAAAGTACTACTATAATACGTACTATAAATACAAAGGTGCGATTGA  >  USA300TCH1516_ALE/2541222‑2541488

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