Predicted mutation
evidence seq id position mutation annotation gene description
RA USA300TCH1516_ALE 2,086,905 (A)8→7 intergenic (‑56/‑629) scn_3 ← / → chp Staphylococcal complement inhibitor/Chemotaxis inhibitory protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE2,086,8980A.97.2% 155.2 / ‑1.6 36intergenic (‑49/‑636)scn_3/chpStaphylococcal complement inhibitor/Chemotaxis inhibitory protein
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base . (21/14);  total (21/15)
Fisher's exact test for biased strand distribution p-value = 4.17e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.75e-01

CTAGCTTGTGCCTCATTTTTATCTAGATTAGTTACTAGTGGTGATGCTAAAACGATTGCTAAAGTTCCCGCAAGTATAGATTTTCTAATTTTCATTTTTTATTCTCCTTAGTAATTAGTACTAATTTATTTTAATGACTGACTAAAAAAAAGTACATTAATTATTCATTAACTATATAGATCTTCTTATTAAATGAAAATTAATAAATTATCAAATTGACATACCTTATTGTATTTATAGAAAAACAAAAAAAGGTAAGCACCGAAATGCTTACCTACTTCCCATAA  >  USA300TCH1516_ALE/2086755‑2087041
                                                                                                                                               |                                                                                                                                               
cTAGCTTGTGCATCATGTTTATCTAGATTAGTTACTAGTGGTGATGCTAAAACGATTGCTAAAGTTCCCGCAAGTATAGATTTTCTAATTTTCATTTTTTATTCTCCTTAGTAATTAGTACTAATTTATTTTAATGACTGACTaaaaa                                                                                                                                             <  4:47894/148‑1 (MQ=255)
    cTTGTGCCTCATTTTTATCTAGATTAGTTACTAGTGGTGATGCTAAAACGATTGCTAAAGTTCCCGCAAGTATAGATTTTCTAATTTTCATTTTTTATTCTCCTTAGTAATTAGTACTAATTTATTTTAATGACTGACTaa                                                                                                                                                <  1:64904/141‑1 (MQ=255)
       gtgCCTCATTTTTATCTAGATTAGTTACTAGTGGTGATGCTAAAACGATTGCTAAAGTTCCCGCAAGTATAGATTTTCTAATTTTCATTTTTTATTCTCCTTAGTAATTAGTACTAATTTATTTTAATGACTGACTAAAAAAAGTACa                                                                                                                                     >  4:1579/1‑148 (MQ=255)
                  ttATATAGATTAGTTAATAGAGGAGATGCTAAAATGATTGTTAAAGTTCCCGAAAGTATAGATTTTCTAATTTTCATTTTTTATTCAAATTAGTAATTAGTACTAATTTATTTTAAAAAAAAAAAAAAAAAAGTACATTAATTATTCa                                                                                                                          <  4:258644/148‑1 (MQ=255)
                          aTTAGTTACTAGTGGTGATGCTAAAACGATTGCTAAAGTTCCCGCAAGTATAGATTTTCTAATTTTCATTTTTTCTTCTCCTTAGTAATTAGTACTAATTTCTTTTTATGACTGACTTAAAAAAGTACATTAATTATCCATCAACtat                                                                                                                  >  4:164354/1‑148 (MQ=255)
                                 aCTAGTGGTGATGCTAAAACGATTGCTAAAGTTCCCGCAAGTATAGATTTTCTAATTTTCATTTTTTATTCTCCTTAGTAATTAGTACTAATTTATTTTAATGACTGACTAAAAAAAGTACATTAATTATTCATTAACTATATAGAtc                                                                                                           >  4:30092/1‑148 (MQ=255)
                                 aCTAGTGGTGATGCTAAAACGATTGCTAAAGTTCCCGCAAGTATAGATTTTCTAATTTTCATTTTTTATTCTCCTTAGTAATTAGTACTAATTTATTTTAATGACTGACTAAAAAAAGTACATTAATTATACATTAACtat                                                                                                                  >  1:34835/1‑141 (MQ=255)
                                    aGTGGTGATGCTAAAACGATTGCTAAAGTTCCCGCAAGTATAGATTTTCTAATTTTCATTTTTTATTCTCCTTAGTAATTAGTACTAATTTATTTTAATGACTGACTAAAAAAAGTACATTAATTATTCATTAACTATATAGAtcttc                                                                                                        <  3:179561/148‑1 (MQ=255)
                                    aGTGGTGATGCTAAAACGATTGCTAAAGTTCCCGCAAGTATAGATTTTCTAATTTTCATTTTTTATTCTCCTTAGTAATTAGTACTAATTTATTTTAATGACTGACTAAAAAAAGTACATTAATTATTCATTAACTATATAAAtcttc                                                                                                        >  4:133286/1‑148 (MQ=255)
                                    aGTGGNGATGCTAAAACGATTGCTAAAGTTCCCGCAAGTATAGAATTTCTAATTTTCATTTTTTATTCTCCTTAGTAATTAGTACTAATTTATTTTAATGACTGACTNAAAAAAGTACATTAATTATTCATTAACTATATAGAtcttc                                                                                                        <  3:72797/148‑1 (MQ=255)
                                     gtggtgATGCTAAAACGATTGCTAAAGTTCCCGCAAGTATAGATTTTCTAATTTTCATTTTTTATTCTCCTTAGTAATTAGTACTAATTTATTTTAATGACTGACTAAAAAAAGTACATTAATTATTCATTAACTATATAGAtcttct                                                                                                       >  3:20604/1‑148 (MQ=255)
                                          taTGCTAAAACTATTGCTAAAGTGCCCTCAATTATAGATTTTCTAATTTTCATTTTTTAGTCGCCTTAGGAATTAGTACTAATTTATTTTAATGACTGACTAAAAAAAGTACATTAATTATTCATTAACTATATAGATCTTCttatta                                                                                                  <  3:331272/147‑1 (MQ=255)
                                          gATGCTAAAACGATTGCTAAAGTTCCCGCAAGTATAGATTTTCTAATTTTCATTTTTTATTCTCCTTAGTAATTAGTACTAATTTATTTTAATGACTGTATATAAAAAGTACATTAATTATTCATTAACTATATAGATCTTCttatta                                                                                                  <  3:303871/148‑1 (MQ=255)
                                          gATGCTAAAACGATTGCTAAAGTTCCCGCAAGTATAGATTTTCTAATTTTCATTTTTTATTCTCCTTAGTAATTAGTACTAATTTATTTTAATGACTGACTAAAAAAAGTACATTAATTATTCATTAACTATATAGATCTTCttatta                                                                                                  <  3:166801/148‑1 (MQ=255)
                                                    cGATTGCTAAAGTTCCCGCAAGTATAGATTTTCTAATTTTCATTTTTTATTCTCCTTAGTAATTAGTACTAATTTATTTTAATGACTGACTAAAAAAAGTACATTAATTATTCATTAACTATATAGATCTTCTTATTAAATGAAAAtt                                                                                        <  4:163047/148‑1 (MQ=255)
                                                        tGCTAAAGTTCCCGCAAGTATAGATTTTCTAATTTTCATTTTTTATTCTCCTTAGTAATTAGTACTAATTTATTTTAATGACTGACTAAAAAAAGTACATTAATTATTCATTAACTATATAGATCTTCTTATTAAATGAAAATtaata                                                                                    <  4:196097/148‑1 (MQ=255)
                                                         gCTAAAGTTCCCGCAAGTATAGATTTTCTAATTTTCATTTTTTATTCTCCTTAGTAATTAGTACTAATTTATTTTAATGACTGACTAAAAAAAGTACATTAATTATTCATTAACTATATAGATCTTCTTATTAAATGAAAATtaataa                                                                                   >  4:395835/1‑148 (MQ=255)
                                                           tAAAGTTCCCGCAAGTATAGATTTTCTAATTTTCATTTTTTATTCTCCTTAGTAATTAGTACTAATTTATTTTAATGACTGACTAAAAAAAGTACATTAATTATTCATTAACTACAGAGCTCCTTGTATGAAATGTAAATTAATAAAt                                                                                 >  4:187512/1‑148 (MQ=255)
                                                                    cGCAAGTATAGATTTTCTAATTTTCATTTTTTATTCTCCTTAGTAATTAGTACTAATTTATTTTAATGACTGACTAAAAAAAGTACATTAATTATTCATTAACTATATAGATCTTCTTATTAAATGAAAATTAATAAATTATCAAAtt                                                                        <  3:90470/148‑1 (MQ=255)
                                                                    cGCAAGTATAGATTTTCTAATTTTCATTTTTTATTCTCCTTAGTAATTAGTACTAATTTATTTTAATGACTGACTAAAAAAAGTACATTAATTATTCATTAACTATATAGATCTTCTTATTAAATGAAAATTAATAAATTATCAAAtt                                                                        <  3:312170/148‑1 (MQ=255)
                                                                    cGCAAGTATAGATTTTCTAATTTTCATTTTTTATTCTCCTTAGTAATTAGTACTAATTTATTTTAATGACTGACTAAAAAAAGTACATTAATTATTCATTAACTATATAGATCTTCTTATTAAATGAAAATTAATAAATTATCAAAtt                                                                        <  3:1579/148‑1 (MQ=255)
                                                                                 tttCTAATTTTCATTTTTTATTCTCCTTAGTAATTAGTACTAATTTATTTTAATGACTGACTAAAAAAAGTACATTAATTATTCATTAACTATATAGATCTTCTTATTAAATGAAAATTAATAAATTATCAAATTGACATACCTTAtt                                                           <  3:219589/148‑1 (MQ=255)
                                                                                         tttCATTTTTTTTTTTTTTTTGTAATTAGAACAAATTTATTTTAATTACTGACTAAAACAAGTAAATAAATTATTCATTAACTATATAGATTTTCTTATTAAATTAAAATTAATAAATTATAAAATTTACATATTTTATTGTAttttt                                                   >  4:84315/1‑146 (MQ=255)
                                                                                                          cTTAGTAATTAGTACTAATTTATTTTAATGACTGACTAAAAAAAGTACATTAATTATTCATTAACTATATAGATCTTCTTATTAAATGAAAATTAATAAATTATCAAATTGACATACCTTATTGTATTTATAGAAAAACAAAAAAAgg                                  >  4:127286/1‑148 (MQ=255)
                                                                                                                aaTTAGTACTAATTTATTTTAATGACTGACTAAAAAAAGTACATTAATTATTCATTAACTATATAGATCTTCTTATTAAATGAAAATTAATAAATTATCAAATTGACATACCTTATTGTATTTATAGAAAAACAAAAAAAGGTAAGCa                            >  3:361092/1‑148 (MQ=255)
                                                                                                                aaTTAGTACTAATTTATTTTAATGACTGACTAAAAAAAGTACATTAATTATTCATTAACTATATAGATCTTCTTATTAAATGAAAATTAATAAATTATCAAATTGACATACCTTATTGTATTTATAGAAAAACAAAAAAAGGTAAGCa                            >  4:79416/1‑148 (MQ=255)
                                                                                                                       aCTAATTTATTTTAATGACTGACTAAAAAAAGTACATTAATTATTCATTAACTATATCGATCTTCTTATTAAATGAAAATTAATAAATTATCAAATTGACATACCTTATTGTATATATAGAAAAACAAAAAAAGGTAAGAACCGAAAt                     >  4:135256/1‑148 (MQ=255)
                                                                                                                       aCTAATTTATTTTAATGACTGACTAAAAAAAGTACATTAATTATTCATTAACTATATAGATCTTCTTATTAAATGAAAATTAATAAATTATCAAATTGACATACCTTATTGTATTTATAGAAAAACAAAAAAAGGTAAGCACCGAAAt                     >  4:392484/1‑148 (MQ=255)
                                                                                                                        cTAATTTATTTTAATGACTGACTAAAAAAAGTACATTAATTATTCATTAACTATATAGATCTTCTTATTAAATGAAAATTAATAAATTATCAAATTGACATACCTTATTGTATTTATAGAAAAACAAAAAAAGGTAAGCACCGAAATg                    >  4:250521/1‑148 (MQ=255)
                                                                                                                           atttatttTAATGACTGACTAAAAAAAGTACATTAATTATTCATTAACTATATAGATCTTCTTATTAAATGAAAATTAATAAATTATCAAATTGACATACCTTATTGTATTTATAGAAAAACAAAAAAAGGTAAGCACCGAAATGCtt                 >  3:185256/1‑148 (MQ=255)
                                                                                                                           atttatttTAATGACTGACTAAAAAAAGTACATTAATTATTCATTAACTATATAGATCTTCTTATTAAATGAAAATTAATAAATTATCAAATTGACATACCTTATTGTATTTATAGAAAAACAAAAAAAGGTAAGCACCGAAATGCtt                 <  3:337822/148‑1 (MQ=255)
                                                                                                                              tatttTAATGACTGACTAAAAAAAGTACATTAATTATTCATTAACTATATAGATCTTCTTATTAAATGAAAATTAATAAATTATCAAATTGACATACCTTATTGTATTTATAGAAAAACAAAAAAAGGTAAGCACCGAAAt                     >  1:11642/1‑141 (MQ=255)
                                                                                                                                  ttAATGACTGACTAAAAAAAGTACATTAATTATTCATTAACTATATAGATCTTCTTATTAAATGAAAATTAATAAATTATCAAATTGACATACCTTATTGTATTTATAGAAAAACAAAAAAAGGTAAGCACCGAAATGCTTACCTACt          >  4:137581/1‑148 (MQ=255)
                                                                                                                                    aaTGACTGACTAAAAAAAGTACATTAATTATTCATTAACTATATAGATCTTCTTATTAAATGAAAATTAATAAATTATCAAATTGACATACCTTATTGTATTTATAGAAAAACAAAAAAAGGTAAGCACCGAAATGCTTAc               <  2:8102/141‑1 (MQ=255)
                                                                                                                                      tgactgacTAAAAAAAGTACATTAATTATTCATTAACTATATAGATCTTCTTATTAAATGAAAATTAATAAATTATCAAATTGACATACCTTATTGTATTTATAGAAAAACAAAAAAAGGTAAGCACCGAAATGCTTACCTACTTccc      <  3:274375/148‑1 (MQ=255)
                                                                                                                                      tgactgacTAAAAAAAGTACATTAATTATTCATTAACTATATAGATCTTCTTATTAAATGAAAATTAATAAATTATCAAATTGACATACCTTATTGTATTTATAGAAAAACAAAAAAAGGTAAGCACCGAAATGCTTACCTACTTccc      >  4:98366/1‑148 (MQ=255)
                                                                                                                                      tgactgacTAAAAAAAGTACATTAATTATTCATTAACTATATAGATCCTCTTATTAAATGAAAATTAATAAATTATCAAATTGACATACCTTATTGTATTTATAGAAAAACAAAAAAAGGTAAGCACCGAAATGCTTACCTACTTccc      >  4:358354/1‑148 (MQ=255)
                                                                                                                                         ctgactAAAAAAAGTACATTAATTATTCATTAACTATATAGATCTTCTTATTAAATGAAAATTAATAAATTATCAAATTGACATACCTTATTGTATTTATAGAAAAACAAAAAAAGGTAAGCACCGAAATGATTACCTACt          >  2:47821/1‑141 (MQ=255)
                                                                                                                                          tgactAAAAAAAGTACATTAATTATTCATTAACTATATAGATCTTCTTATTAAATGAAAATTAATAAATTATCAAATTGACATACCTTATTGTATTTATAGAAAAACAAAAAAAGGTAAGCACCGAAATGCTTACCTACTTCCCATaa  <  3:30092/148‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                               
CTAGCTTGTGCCTCATTTTTATCTAGATTAGTTACTAGTGGTGATGCTAAAACGATTGCTAAAGTTCCCGCAAGTATAGATTTTCTAATTTTCATTTTTTATTCTCCTTAGTAATTAGTACTAATTTATTTTAATGACTGACTAAAAAAAAGTACATTAATTATTCATTAACTATATAGATCTTCTTATTAAATGAAAATTAATAAATTATCAAATTGACATACCTTATTGTATTTATAGAAAAACAAAAAAAGGTAAGCACCGAAATGCTTACCTACTTCCCATAA  >  USA300TCH1516_ALE/2086755‑2087041

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