Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 111,516 C→A 28.8% Q196K (CAA→AAA)  USA300HOU_RS00520 → putative lipoprotein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE111,5160CA28.8% 89.7 / 44.1 59Q196K (CAA→AAA) USA300HOU_RS00520putative lipoprotein
Reads supporting (aligned to +/- strand):  ref base C (26/16);  new base A (6/11);  total (32/27)
Fisher's exact test for biased strand distribution p-value = 8.59e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.92e-01

AGGAAGTTAAAGAAGAAAACATAAAAAAAGAAATCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAGGTACAGGGAATTTAAAAGGTTCATCAGTTGGATATAAAAAAATTGAA  >  USA300TCH1516_ALE/111379‑111656
                                                                                                                                         |                                                                                                                                            
aGGAAGTTAAAGAAAAAAACATAAAAAAAGAAATCGAAAACTTTAAGTTTTTTGCGCCATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAat                                                                                                                                           <  1:3375/141‑1 (MQ=255)
 ggAAGTTAAAGAAGAAAACATAAAAAAAGAAATCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAata                                                                                                                                          <  1:171365/141‑1 (MQ=255)
  gAAGTTAAAGAAGAAAACATAAAAAAAGAAATCGACCTTTAGAAGTTTTTTGGGCCATGTAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAatat                                                                                                                                         >  1:204853/1‑141 (MQ=25)
  gAAGTTAAAGAAGAAAACATAAAAAAAGAAATCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATCTATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAatat                                                                                                                                         <  1:604447/141‑1 (MQ=255)
   aaGTTAAAGAAGAAAACATAATAAAAGAAATCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATc                                                                                                                                        >  1:449074/1‑141 (MQ=255)
   aaGTTAAAGAAGAAAACATAAAAAAAGAAATCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATc                                                                                                                                        >  1:129032/1‑141 (MQ=255)
   aaGTTAAAGAAGAAAACATAAAAAAAGAAATCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCAGATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATc                                                                                                                                        >  1:160518/1‑141 (MQ=255)
         aagaagAAAACATAAAAAAAGAAATCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTaa                                                                                                                                  <  1:269873/141‑1 (MQ=255)
           aaagaaAACATAAAAAAAGAAATCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACt                                                                                                                                <  1:754568/140‑1 (MQ=255)
             agaaAACATAAAAAAAGAAATCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTaa                                                                                                                              <  1:449060/141‑1 (MQ=255)
              gaaAACATAAAAAAAGAAATCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAat                                                                                                                             >  1:279688/1‑141 (MQ=255)
              gaaAACATAAAAAAAGAAATCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAat                                                                                                                             >  1:445129/1‑141 (MQ=255)
              gaaAACATAAAAAAAGAAATCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAat                                                                                                                             >  1:257460/1‑141 (MQ=255)
               aaaaCATAAAAAAAGAAATCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAatg                                                                                                                            >  1:780417/1‑141 (MQ=255)
                aaaCATAAAAAAAGAAATCGAAAACTTGAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAatga                                                                                                                           <  1:174552/141‑1 (MQ=255)
                  aCATAAAAAAAGAAATCGAAAACTTTAAGTTTTTTGCGCAATATCGCAATTTGAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAatgatg                                                                                                                         <  1:496234/141‑1 (MQ=255)
                     tAAAAAAAGAAATCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGAtt                                                                                                                      >  1:755535/1‑141 (MQ=255)
                     tAAAAAAAGAAATCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTGATTCGGCTAAATATCAATTAACTAATGATGAtt                                                                                                                      >  1:159419/1‑141 (MQ=255)
                       aaaaaaGAAATCGAAAACTTTAAGGTTTGTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATATAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATtat                                                                                                                    <  1:414182/141‑1 (MQ=255)
                        aaaaaGAAATCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTATCTAATGATGATtata                                                                                                                   <  1:648208/141‑1 (MQ=255)
                        aaaaaGAAATCGAAAACTTGAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATtata                                                                                                                   <  1:73528/141‑1 (MQ=255)
                          aaaGAAATCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAAt                                                                                                                 <  1:507345/141‑1 (MQ=255)
                          aaaGAAATCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAAt                                                                                                                 <  1:141027/141‑1 (MQ=255)
                           aaGAAATCGAAAACTTGAAGTATTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATg                                                                                                                >  1:445701/1‑141 (MQ=255)
                             gAAATCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTa                                                                                                              >  1:352568/1‑141 (MQ=255)
                              aaaTCGAAAACTTTAAGTTTTTTACGCAATATAGCAATTTTAAAGATTTAATGAATTAGAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTaa                                                                                                             <  1:549779/141‑1 (MQ=255)
                                aTCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTACATATCAATTAACTAATAATGATTATAATGTagaa                                                                                                           >  1:364472/1‑141 (MQ=18)
                                   gAAAACTTTAAGTNTTTTGCGCAAGATCGCNATTTTAAAGATTTAATGAAGTAGAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACaa                                                                                                        <  1:351081/141‑1 (MQ=255)
                                        cTTTAAGTTTTTTGCGCAATATAGCAATTTTAGAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGACTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACg                                                                                                   <  1:580742/141‑1 (MQ=255)
                                         tttAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAGATATCAATTAACTAATGATGATTATAATGTAAAACAATTACgc                                                                                                  <  1:813596/141‑2 (MQ=255)
                                              gTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAg                                                                                             >  1:489668/1‑141 (MQ=255)
                                                 ttttGCGCAATATAGCAATTTTAATTTTTTAATTAATTCCCTCGCCAGATGAATAGCAGATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGata                                                                                          <  1:629337/141‑1 (MQ=14)
                                                 ttttGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGata                                                                                          <  1:514287/141‑1 (MQ=255)
                                                 ttttGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGata                                                                                          <  1:514511/141‑1 (MQ=255)
                                                    tGCGCAATATAGCAATTTTAAAGATTAAATGAATTATAAAGAATGAGATAAATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGa                                                                                       <  1:419259/141‑1 (MQ=14)
                                                       gcAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATTATGATTATAATGTAAACAAATTAAGCAAAAGatataatat                                                                                    >  1:33399/1‑136 (MQ=17)
                                                                cAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGAGATGGTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATAccgacgag                                                                           <  1:633716/141‑7 (MQ=11)
                                                                 aaTTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATAcaaacaaaa                                                                          >  1:800136/1‑140 (MQ=255)
                                                                  aTTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATa                                                                         <  1:782155/141‑1 (MQ=255)
                                                                    tttAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATaaa                                                                       >  1:185283/1‑141 (MQ=255)
                                                                               aaTGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAGGCt                                                            >  1:283433/1‑141 (MQ=255)
                                                                               aaTGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAACt                                                            >  1:259248/1‑141 (MQ=255)
                                                                                aTGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCtg                                                           >  1:32163/1‑141 (MQ=255)
                                                                                aTGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCtg                                                           >  1:728558/1‑141 (MQ=255)
                                                                                 tGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCtgt                                                          <  1:185902/141‑1 (MQ=255)
                                                                                 tGAAGTATAAAGATGGCGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCtgt                                                          <  1:33014/141‑1 (MQ=255)
                                                                                             aTGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAGGTa                                              <  1:600545/141‑1 (MQ=255)
                                                                                              tGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAGGTAc                                             >  1:792677/1‑141 (MQ=255)
                                                                                                   atatatCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAGGTACAGGGa                                        >  1:140966/1‑141 (MQ=255)
                                                                                                     atatCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAGGTACAGGGAAt                                      >  1:481854/1‑141 (MQ=255)
                                                                                                         cATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCCAAGCTGTTGTTGAAAGGTACAGGGAATTTaa                                  >  1:765766/1‑141 (MQ=255)
                                                                                                                  cAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATAATAATTATAATGTAAAACTATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAATGTACAGGGAATTTAAAAGGTtcat                         >  1:381457/1‑141 (MQ=255)
                                                                                                                     agGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAGGTACAGGGAATTTAAAAGGTTCATCAg                      >  1:721654/1‑141 (MQ=255)
                                                                                                                      ggTGCCGAGTTATTCAGCTCAATATCAAATAACTAATGATGATTATAATGTAAAACACTTACGTAAAAGATATGCCATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAGGAACAGGGAAAATAAACGGCTcagcaga                     >  1:384786/1‑140 (MQ=18)
                                                                                                                         gCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAGGTACAGGGAATTTAAAAGGTTCATCAGTTgg                  >  1:40647/1‑141 (MQ=255)
                                                                                                                           cGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAGGTACAGGGAATTTAAAAGGTTCATCAGTTGGat                >  1:671457/1‑141 (MQ=255)
                                                                                                                            gAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAGGTACAGGGAATTTAAAAGGTTCATCAGTTGGata               <  1:304531/141‑1 (MQ=255)
                                                                                                                             aGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAGGTACAGGGAATTTAAAAGGTTCATCAGTTGGatat              >  1:55827/1‑141 (MQ=255)
                                                                                                                             aGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAGGTACAGGGAATTTAAAAGGTTCATCAGTTGGatat              >  1:536166/1‑141 (MQ=255)
                                                                                                                                         cAATAGTAATTACTTAAAGACGATTATAAGGTAAAACAATTACGTAAAAGATATCACATACCAAAAAAAAAAGATCCGAAGCGGTTGGTGTAAGGTACAGGGAATTTAAAAGGTTCATCAGTTGGATATAAAAAAATTGaa  <  1:309252/141‑1 (MQ=14)
                                                                                                                                         |                                                                                                                                            
AGGAAGTTAAAGAAGAAAACATAAAAAAAGAAATCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAGGTACAGGGAATTTAAAAGGTTCATCAGTTGGATATAAAAAAATTGAA  >  USA300TCH1516_ALE/111379‑111656

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: