Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 1,232,148 G→C 68.1% G18A (GGC→GCC)  prkC → Serine/threonine‑protein kinase PrkC

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE1,232,1480GC68.1% 122.1 / 87.2 97G18A (GGC→GCC) prkCSerine/threonine‑protein kinase PrkC
Reads supporting (aligned to +/- strand):  ref base G (16/15);  new base C (29/36);  total (46/51)
Fisher's exact test for biased strand distribution p-value = 6.62e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.04e-01

GGTGATCAATTAATGCAATTGGCATTAGATAACCATTCGAAAGATAACGTTACTTTCATACTCGCGGCTATTGAAGGTGATAAAGTATGATAGGTAAAATAATAAATGAACGATATAAAATTGTAGATAAGCTTGGCGGCGGTGGCATGAGTACCGTTTATCTTGCTGAAGATACGATACTTAACATTAAAGTTGCAATTAAGGCGATTTTTATACCACCTAGAGAAAAAGAAGAAACATTAAAACGTTTTGAACGAGAAGTACATAACTCAT  >  USA300TCH1516_ALE/1232010‑1232282
                                                                                                                                          |                                                                                                                                      
ggTGATCAATTAATGCAATTGGCATTAGATAACCATTCGAAAGATAACGTTACTTTCATACTCGCGGCTATTGAAGGTGATAAAGTATGATAGGTAAAATAATAAATGAACGATATAAAATTGTAGATAAGCTTGgcgccg                                                                                                                                      <  2:803949/141‑1 (MQ=255)
  ggtacAATTAAGTCAATTGGCATTAGATAACCATTCGAAAGATAACGTTACTTTCATTCTCGCGGCTATTGAAGGTGATAAAGGATGATAGGTAAAATAATAAATGAACGATATAAAATTGTAGATAAGCTTGGCGGCGGt                                                                                                                                    <  2:559243/137‑1 (MQ=255)
   gATCAATTAATGCAATTGGCATTAGATAACCATTCGAAAGATAACGTTACTTTCATACTCGCGGCTATTGAAGGTGATAAAGTATGATAGGTAAAATAATAAATGAACGATATAAAGTTGTAGATAAGCGTGGcggcggag                                                                                                                                   >  1:637166/1‑139 (MQ=255)
      aaactAATGCAACTGGCATTAGATAACCATTCGAAACGTAACCCTACGTTCATACTCGCGTCTATTGAGGGTGATAAAGTATGATAGTTAATATTATAAATGAACGTTATCAAATTGTAGATAAGCTTGGCGCCGGTGGCa                                                                                                                                <  2:192527/137‑1 (MQ=255)
       aaTTAATGCCATTGGCAGTAGATAACCATTCGAAAGATAACGTTACTTGCATACTCGCGGCTATTGAAGGTGATAAAGTATGATAGGTAAAATAATAAATGAACGATATAAAATTGTAGATAAGCTTGGCGCCGGTGGCAt                                                                                                                               <  2:271067/141‑1 (MQ=255)
        aTTAATGCAATTGGCATTAGATAACCATTCGAAAGATAACGTTACTTTCATACTCGCGGCTATTGAAGGTGATAAAGTATGATAGGTAAAATAATAAATGAACGATATAAAATTGTAGATAAGCTTGGCGGCGGTGGCATg                                                                                                                              <  1:872336/141‑1 (MQ=255)
         ttAATGCAATTGGCATTAGATAACCATTCGAAAGATAACGTTACTTTCATACTCGCGGCTATTGAAGGTGATAAAGTATGATAGGTAAAATAATAAATGAACGATATAAAATTGTAAATAAGCTTGGCGGCGGTGGCATGa                                                                                                                             >  1:895508/1‑141 (MQ=255)
          tAATGCAATTGGCATTAGATAACCATTCGAAAGATAACGTTACTTTCATACTCGCGGCTATTGAAGGTGATAAAGTATGATAGGTAAAATAATAAATGAACGATATAAAATTGTAGATAAGCTTGGCGCCGGTGGCATGAg                                                                                                                            >  1:686332/1‑141 (MQ=255)
           aaTGCAATTGGCATTAGATAACCATTCGAAAGATAACGTTACTTTCATACTCGCGGCTATTGAAGGTGATAAAGTATGATAGGTAAAATAATAAATGAACGATATAAAATTGTAGATAAGCTTGGCGCCGGTGGCATGact                                                                                                                           >  2:270151/1‑139 (MQ=255)
            aTGCAATTGGCATTAGATAACCATTCGAAAGATAACGTTACTTTCATACTCGCGGCTATTGAGGGTGATAAAGAATAATAGGTAAAATAATAAATGAACGATATAAAATTGTAGATAAGCGTGGCGGCGGATGAATGAGTa                                                                                                                          >  2:991510/1‑141 (MQ=255)
             tGCAATTGGCATTAGATAACCATTCGAAAGATAACGTTACTTTCATACTCGCGGCTATTGAAGGTGATAAAGTATGATAGGTAAAATAATAAATGAACGATATAAAATTGTAGATAAGCTTGGCGCCGGTGGCATGAGTAc                                                                                                                         >  1:959852/1‑141 (MQ=255)
              gCAATTGGCATTAGATAACCATTCGAAAGATAACGTTACTTTCATACTCGCGGCTATTGAAGGTGATAAAGTATGATAGGTAAAATAATAAATGAACGATATAAAATTGTAGATAAGCTTGGCGGCGGTGGGATGAGTAct                                                                                                                        >  2:644184/1‑140 (MQ=255)
              gCAATTGGCATTAGATAACCATTCGAAAGATAACGTTACTTTCATACGCGCGGCTATTGAAGGTGATAAAGTATGATAGGTAAAATAATAAATGAACGATATAAAATTGTAGATAAGCTTGGCGCCGGTGACATGAGTAcc                                                                                                                        <  2:810435/141‑1 (MQ=255)
               taaTTGGCAGTATATAACCACTCGAAAGATAACGTTACTTTCATACTCGCGGCTATTGAAGGTGATAAAGTATGATAGGTAAAATAATAAATTAACGATATAAAATTGTAGATAAGCTTGGCGGCGGTGGCATGAGTACCg                                                                                                                       <  1:510274/140‑1 (MQ=255)
                 aTTGGCATTAGATAACCATTCGAAAGATAACGTTACTTTCATACTCGCGGCTATTGAAGGTGATAAAGTATGATAGGTAAAATAATAAATGAACGATATAAAATTGTAGATAAGCTTGGCGCCGGTGGCATGAGTACCGtt                                                                                                                     >  2:109796/1‑141 (MQ=255)
                 aTTGGCATTAGATAACCATTCGAAAGATAACGTTACTTAAAAACTCGCGGCTATGGAAGGTGAAAAAGTATGATAGGTAAAATAATAAATAAACGATATAAAATTGTCGATAAGCAAGGCGCCGGGAGCATGAGTACCGtt                                                                                                                     >  2:1133544/1‑141 (MQ=255)
                   tGGAATTAGATAACCATTCGAAAGATAACGTTACTTTCATACTCGCGGCTATTGAAGGTGATAAAGTATGATAGGTAAAATAATAAATGAACGATATAAAATTGTAGATAAGCTTGGCGCCGGTGGCATGAGTACCGTTTa                                                                                                                   <  2:1040358/141‑1 (MQ=255)
                    ggCATTAGATAACCATTCGAAAGATAACGTTACTTTCATACTCGCGGCTATTGAAGGTGATAAAGTATGATAGGTAAAATAATAAATGAACGATATAAAATTGTAGATAAGCTTGGCGGCGGTGGCATGAGTACCGTTTAt                                                                                                                  <  1:907154/141‑1 (MQ=255)
                     gCATGAGATAACCATTCGAAAGATAACGTTACTTTCATACTCGCGGCAAATGAAGGTGATAAAGTATGATAGGTAAAATAATAAATGAACGATATAAAATTGTAGATAAGCTTGGCGCCGGGGGGCTGGGTACCGCTTATc                                                                                                                 >  2:481457/1‑141 (MQ=255)
                     acATTAGATAACCATTCGAAAGATAACGTTACTTTCATACTCGCGGCTATTGAAGGTGATAAAGTATGATAGGTAAAATAATAAATGAACGATATAAAATTGTATATAAGCTTGACGACGGGGGCATGGGTAACGCTTATc                                                                                                                 >  2:311943/2‑141 (MQ=255)
                         tAGATAACCATTCGAAAGATAACGTTACTTTCATACTCGCGGCTATTGAAGGTGATAAAGTATGATAGGTAAAATAATAAATGAACGATATAAAATTGTAGATAAGCTTGGCGCCGGTGGCATGAGTACCGTTTATCTTGc                                                                                                             <  1:802032/141‑1 (MQ=255)
                           gATAACCATTCGAAAGATGACGTTACTTTCATACTCGCGGCTATTGAAGGTGATAAAGTATGATAGGTAAAATAATAAATGAACGGTATAAAATTGTAGATAAGCTTGGCGCCGGTGGCATGAGTACCGTTTATCTTGCTg                                                                                                           <  1:775505/141‑1 (MQ=255)
                           gATAACCATTCGAAAGATAACGTTACTTTCATACTCGCGGCTATTGAAGGTGATAAAGTATGATAGGTAAAATAATAAATGAACGATATAAAATTGTAGATAAGCTTGGCGGCGGTGGCATGAGTACCGTTTATCTTGCTg                                                                                                           <  2:1046351/141‑1 (MQ=255)
                                ccATTCGAAAGATAACGTTACTTTCATACTCGCGGCTATTGAAGGTGATAAAGTATGATAGGTAAAATAATAAATGAACGATATAAAATTGTAGATAAGCTTGGCGCCGGTGGCATGAGTACCGTTTATCTTGCTGAAgat                                                                                                      <  1:1056493/141‑1 (MQ=255)
                                ccATTCGAAAGATAACGTTACTTTCATACTCGCGACTATTGAAGGTGATAAAGTATGATAGGTAAAATAATAAATGAACGATATAAAATTGTAGATAAGCTTGGCGCCGGTGGCATGAGTACCGTTTATCTTGCTGAAgat                                                                                                      <  1:994643/141‑1 (MQ=255)
                                  aTTCGAAGCCTAACGTTACTTCAATTCGCGCCGCTATGGAAGGTGATAAAGTAAGATAGGTAAAATAATAAATGAACGATATAAAAGTGTAGATACGCTTGGCGGCGGTGGCATGAGTACCGTTTATCTTGCTGAAgatac                                                                                                    <  1:108953/141‑1 (MQ=255)
                                  aTTCGAAAGATAACGTTACTTTCATACTCGCGGCTATTGAAGGTGATAAAGTATGATAGGTAAAATAATAAATGAACGATATAAAATTGTAGATAAGCTTGGCGCCGGTGGCATGAGTACCGTTTATCTTGCTGAAgatac                                                                                                    <  2:907962/141‑1 (MQ=255)
                                  aTTCGAAAGATAACGTTACTTTCATACTCGCGGCTATTGAAGGTGATAAAGTATGATAGGTAAAATAATAAATGAACGATATAAAATTGTAGATAAGCTTGGCGCCGGTGGCATGAGTACCGTTTATCTTGCTGAAgatac                                                                                                    <  1:1021579/141‑1 (MQ=255)
                                   ttCGAAAGATAACGTTACTTTCATACTCGCGGCTATTGAAGGTGATAAAGTATGATAGGTAAAATAATAAATGAACGATATAAAATTGTAGATAAGCTTGGCGCCGGTGGCATGAGTACCGTTTATCTTGCTGAAgatacg                                                                                                   <  2:428087/141‑1 (MQ=255)
                                   ttCGAAAGATAACGTGACTTTCATACTCGCGGCTATTGAAGGTGATAAAGTATGATAGGTAAAATAATAAATGAACGATATAAAATTGTAGATAAGCTTGGCGCCGGTGGCATGAGTACCGTTTATCTTGCTGAAgttacg                                                                                                   >  1:887783/1‑141 (MQ=255)
                                    tCGAAAGATAACGTTACTTTCATACTCGCGGCTATTGAAGGTGATAAAGTATGATAGGTAAAATAATAAATGAACGATATAAAATTGTAGATAAGCTTGGCGCCGGTGGCATGAGTACCGTTTATCTTGCTGAAgatacga                                                                                                  <  2:458311/141‑1 (MQ=255)
                                          gaaaatgTTACTTTCATACTCGCGGCCAGTGAAGGTGATAAAGTATGACAGGTAAAATAATAAATGAACGATATAAAATTGTAGATAAGCTTGGCGCCGGTGGCATGAGTACCGTTTATCTTGCTGCAGATACGGTACGTa                                                                                            >  2:911528/7‑141 (MQ=255)
                                          gATAACGTTACTTTCATACTCGCGGCTATTGAAGGTGATAAAGTATGATAGGTAAAATAATAAATGAACGATATAAAATTGTGGATGAGCTCGGCGCCGGTGGCATGAGTACCGTTTATCTTGCTGAAGATACGATACTTa                                                                                            >  1:461313/1‑141 (MQ=255)
                                            tAACGTTACTTTCATACTCGCGGCTATTGAAGGTGATAAAGTATGATAGGTAAAATAATAAATGAACGATATAAAATTGTAGATAAGCTTGGCGCCGGTGGCATGAGTACCGTTTATCTTGCTGAAGATACGATACTTAAc                                                                                          >  1:974002/1‑141 (MQ=255)
                                             aaCGTTACTTTCATACTCGCGGCTATTGAAGGTGATAAAGTATGATAGGTAAAATAATAAATGAACGATATAAAATTGTAGATAAGCTTGGCACCGGTGGCATGAGTACCGTTTATCTTGCTGAAGATACGATACTTAACa                                                                                         <  2:778689/141‑1 (MQ=255)
                                               cGTTACTTTCATACTCGCGGCTATTGAAGGTGAGAAAGTATGATAGGTAAAATAATAAATGAACGCTGTAAAATTGTAGGTAAGCTTGGCGCCGGTGGCATGAGTACCGTTTATCTTGCTGAAGATACGATGCTTAACAtt                                                                                       >  1:30660/1‑141 (MQ=255)
                                                gTTACTTTCATACTCGCGGCTATTGAAGGTGATAAAGTATGATAGGTAAAATAATAAATGAACGATATAAAATTGTAGATAAGCTTGGCGCCGGTGGCATGAGTACCGTTTATCTTGCTGAAGATACGATACTTAACATTa                                                                                      >  1:744671/1‑141 (MQ=255)
                                                  tACTTTCATACTCGCGGCTAATGAAGGTGATAAAGTATGATAGGTAAAATAATAAATGAACGATATAAAATTGTAGATAAGCTTGGCGCCGGGGGCATGGGTACGGTTTATCTTGCTGAAGATCACATACTTAACATTaaa                                                                                    >  1:69008/1‑141 (MQ=255)
                                                   aCTTTCATACTCGCGGCTATTGAAGGTGATAAAGTATGATAGGTAAAATAATAAATGAACGATATAAAATTGTAGATAAGCTTGGCGCCGGTGGCATGAGTACCGCTTATCTTGCTGAAGATACGATACTTAACAtagcag                                                                                   >  1:1022297/1‑136 (MQ=255)
                                                    cTTTCATACTCGCGGCTATTGAAGGTGATAAAGTATGATAGGTAAAATAATAAATGAACGATATAAAATTGTAGATAAGCTTGGCGGCGGTGGCATGAGTACCGTTTATCTTGCTGAAGATACGATACTTAACATTAAAGt                                                                                  <  1:1085524/141‑1 (MQ=255)
                                                      ttgctACTCGCGGCTATTCAGGGTCATGAAGTATGATAGGTAAAATAATAAATGAACGATATAAAATTGTAGATAAGCTTGGCGGCGGTGGCATGAGTACCGTTTATCTTGCTGAAGATACGATACTTAACATTAAAGTTg                                                                                <  2:835913/137‑1 (MQ=255)
                                                         aTACTCGCGGCTATTGAAGGTGATAAAGTATGATAGGTAAAATAATAAATGAACGATATAAAATTGTAGATAAGCTTGGCGCCGGTGGCATGAGTACCGTTTATCTTGCTGAAGATACGATACTTAACATTAAAGTTGCaa                                                                             >  1:181438/1‑141 (MQ=255)
                                                          tACTCGCGGCTATTGAAGGTGATAAAGTATGATAGGTAAAATAATAAATGAACGATATAAAATTGTAGATAAGCTTGGCGCCGGTGGCATGAGTACCGTTTATCTTGCTGAAGATGCGATACTTAACATTAAAGTTGCatt                                                                            >  2:708871/1‑139 (MQ=255)
                                                              tcggggtATTGAAGGTGATAAAGTATGATAGGTAAAATAATAAATGAACGATATAAAAGTGTAGATAAGCTTGGCGCCGGTGGCATGAGTACCGTTTATCTTGCTGAAGATACGATACTTAACATTAAAGTTGCAATTAAg                                                                        <  1:537384/135‑1 (MQ=255)
                                                              cgcgGCTATTGAAGGTGATAAAGTATGATAGGTAAAATAATAAATGAACGATATAAAATTGTAGATAAGCTTGGCGCCGGTGGCATGAGTACCGTTTATCTTGCTGAAGATACGATACTTAACATTAAAGTTGCAATTAAg                                                                        >  2:614690/1‑141 (MQ=255)
                                                               gcgGCTATTGAAGGTGATAAAGTATGATAGGTAAAATAATAAATGAACGATATAAAATTGTAGATAAGCTTGGCGCCGGTGGCATGAGTACCGTTTATCTTGCTGAAGATACGATACTTAACATTAAAGTTGCAATTAAgg                                                                       >  2:435571/1‑141 (MQ=255)
                                                                 ggCTATTGAAGGTGATAAAGTATGATAGGTAAAATAATAAATGAACGATATAAAATTGTAGATAAGCTTGGCGCCGGTGGCATGAGTACCGTTTATCTTGCTGAAGATACGATACTTAACATTAAAGTTCGAATGAAGGCg                                                                     >  2:835280/1‑141 (MQ=255)
                                                                    tATTGAATGTGATAAAGTAGGATAGGGAAAATAATAAATGAACGATATAAAATTGTAGATAAGCTTGGCGGCGGTGGCATGAGTACCGTTTATCTTGCTGAAGATACGATACTTAACATTAAAGTTGCAATTAAGGCGAtt                                                                  <  2:169451/141‑1 (MQ=255)
                                                                     aTTGAAGGTGATAAAGTATGATAGGTAAAATAATAAATGAACGATATCCCCTTGTAGATGAGCTTGGCGCCGGTGGCATGAGTACCGTTTATCTTGCTGAAGATACGATACTTAACATTAAAGTTGCAATTAAGGCGAttt                                                                 <  2:1149590/141‑1 (MQ=255)
                                                                      ttGAAGGTGATAAAGTATGATAGGTAAAATAATAAATGAACGATATAAAATTGTAGATAAGCTTGGCACCGGTGGCATGAGTACCGTTTATCTTGCTGAAGATACGATACTTAACATTAAAGTTGCAATTAAGGCGAtttt                                                                >  2:333707/1‑141 (MQ=255)
                                                                      ttGAAGGGTATAAAGTATGATAGGTAAAATAATAAATGAACGATATAAAATTGTAGATAAGCTTGGCGCCGGTGGCATGAGTCCCGTTTATCTTGCTGAAGATACGATACTTCACATTAAAGTTGCAATTAAGGCGAtttt                                                                <  2:937960/141‑1 (MQ=255)
                                                                      ttGAAGGGGATAAAGTATTATAGGTAAAATAATAAATGAACGATATAAAATTGTAGATAAGCTTGGCGGCGGTGGCATGAGTACCGTTTATCTTGCTGAAGATACGATACTTAACATTAAAGTTGCAATTAAGGCGAtttt                                                                <  2:180551/141‑1 (MQ=255)
                                                                       tGAAGGTGATACAGCAATATAGGTAAAATAATAATGGAACGTCTTCGGTTCATACGCGAGCTTGGCGCCGGTGGCATTAGTACCGTTTATCTTGCTGAAGATACGATACTTAACATTAAAGTTGCAATTAAGGCGAttttt                                                               <  1:270151/141‑1 (MQ=255)
                                                                       tGAAGGTGATAAAGTATGATAGGTAAAATAATAAATGAACGATATAAAATTGTAGATAAGCTTGGCGGCGGTGGCATGAGTACCGTTTATCTTGCTGAAGATACGATACTTAACATTAAAGTTGCAATTAAGGCGAttttt                                                               <  2:939499/141‑1 (MQ=255)
                                                                        gAAGGTGATAAAGTATGATAGGTAAAATAATAAATGAACGATATAAAATTGTAGATAAGCTTGGCGGCGGTGGCATGAGTACCGTTTATCTTGCTGAAGATACGATACTTAACATTAAAGTTGCAATTAAGGCGATTTTta                                                              >  1:709984/1‑141 (MQ=255)
                                                                        gAAGGTGATAAAGTATGATAGGTAAAATAATAAATGAACGATATAAAATTGTAGATAAGCTTGGCGGCGGTGGCATGAGTACCGTTTATCTTGCTGAAGATACGATACTTAACATTAAAGTGGCAATTAAGGCGAtttttt                                                              >  2:1115033/1‑140 (MQ=255)
                                                                         aaGGTGATAAAGTATGATAGGTAAAATAATAAATGAACGATATAAAATTGTAGATAAGCTTGGCGCCGGTGGCATGAGTACCGTTTATCTTGCTGAAGATACGATACTTAACATTAAAGTTGCAATTAAGGCGATTTTtat                                                             <  1:897055/141‑1 (MQ=255)
                                                                         aaGGTGATAAAGTATGATAGGTAAAATAATAAATGAACGATATAAAATTGTAGATAAGCTTGGCGCCGGTGGCATGAGTACCGTTTATCTTGCTGAAGATACGATACTTAACATTAAAGTTGCAATTAAGGCGATTTTtat                                                             <  1:179216/141‑1 (MQ=255)
                                                                          aGGTGATAAAGTATGATAGGTAAAATAATAAATGAACGATATAAAATTGTAGATAAGCTTGGCGGCGGTGGCATGAGTACCGTTTATCTTGCTGAAGATACGATACTTAACATTAAAGTTGCAATTAAGGCGATTTTtata                                                            >  2:306359/1‑141 (MQ=255)
                                                                          aGGTGATAAAGTATGATAGGTAAAATAATAAATGAACGATATAAAATTGTAGATAAGCTTGGCGCCGGTGGCATGAGTACCGTTTATCTTGCTGAAGATACGATACTTAACATTAAAGTTGCAATTAAGGCGATTTTtata                                                            <  1:166635/141‑1 (MQ=255)
                                                                                tAAAGTATGATAGGTAAAATAATAAATGAACGATACAAAATTGTAGATAAGCTTGGCGCCGGTGGCATGAGTACCGTTTATCTTGCTGAAGATACGATACTTAACATTAAAGTTGCAATTAAGGCGATTTTTATACCACCt                                                      <  2:651577/141‑1 (MQ=255)
                                                                                tAAAGTAAGATAGGTAAAATAATAAATGAACGATATAAAATTGTAGATAAGCTTGGCGGCGGTGGCATGAGTACCGTTTATCTTGCTGAAGATACGATACTTAACATTAAAGTTGCAATTAAGGCGATTTTTATACCACCt                                                      >  2:1188367/1‑141 (MQ=255)
                                                                                tAAAGGATGATAGGTAAAATAATAAATGAACGATATAAAATTGTAGATAAGCTTGGCGCCGGTGGCATGAGTACCGTTTATCTTGCTGAAGATACGATACTTAACATTAAAGTTGCAATTAAGGCGATTTTTATACCACCt                                                      <  2:119725/141‑1 (MQ=255)
                                                                                  aaGTATGATAGGTAAAATAATAAATGAACGATATAAAATTGTAGATAAGCTTGGCGCCGGTGGCATGAGTACCGTTTATCTTGCTGAAGATACGATACTTAACATTAAAGTTGCAATTAAGGCGATTTTTATACCACCTag                                                    <  1:1073131/141‑1 (MQ=255)
                                                                                     tATGATAGGTAAAATAATAAATGAACGATATAAAATTGTAGATAAGCTTGGCGGCGGTGGCATGAGTACCGTTTATATTGCTGAAGATACGATACTTAACATTAAAGTTGCAATTAAGGCGATTTTTATACCACCTAgaga                                                 >  2:51813/1‑141 (MQ=255)
                                                                                      atgtaaGGCAAAAGAATAAATCAACGACATAAAATTGTAGATACGCTTGGCGGCGGTGGCATGAGTACCGTTTATCTTGCTGAAGATACGATACTAAACATTAAAGTTGCAATTAAGGCGATTTTTATACCACCTAGAGaa                                                <  1:444806/136‑1 (MQ=255)
                                                                                          tAGGTAAAATAATAAATGAACGATATAAAATTGTAGATAAGCTTGGCGGCGGTGGCATGAGTACCGTTTATCTTGCTGAAGATACGATACTTAACATTAAAGTTGCAGTTAAGGCGATTTTTATACCACCTAGAGAAAaag                                            >  1:927759/1‑141 (MQ=255)
                                                                                          tAGGTAAAATAATAAATGAACGATATAAAATTGTAGATAAGCTTGGCGCCGGTGGCATGAGTACCGTTTATCTTGCTGAAGATACGATACTTAACATTAAAGTTGCAATTAAGGCGATTTTTATACCACCTAGAGAAAaag                                            >  2:917554/1‑141 (MQ=255)
                                                                                          tAGGTAAAATAATAAATGAACGATATAAAATTGTAGATAAGCTTGGCGCCGGTGGCATGAGTACCGTTTATCTTGCTGAAGATACGATACTTAACATTAAAGTTGCAATTAAGGCGATTTTTATACCACCTAGAGAAAaag                                            <  1:155488/141‑1 (MQ=255)
                                                                                            ggTAAAATAATAAATGAACGATATAAAATTGTAGATAAGCTTGGCGCCGGTGGCATGAGTACCGTTTATCTTGCTGAAGATACTATACTTAACATTAAAGTTGCAATTAAGGCGATTTTTATACCACCTAGAGAAAaagaa                                          >  2:483506/1‑141 (MQ=255)
                                                                                             gTAAAATAATAAATGAACGATATAAAATTGTAGATAAGCTTGGCGCCGGTGGCATGAGTACCGTTTATCTTGCTGAAGATACGATACTTAACATTAAAGTTGCAATTAAGGCGATTTTTATACCACCTAGAGAAAaagaag                                         <  2:481052/141‑1 (MQ=255)
                                                                                                  ataataAATGAACGATATAAAATTGTAGATAAGCTTGGCGGCGGTGGCATGAGTACCGTTTATCTTGCTGAAGATGCGATACTTAACATTAAAGTTGCAATTAAGGCGATTTTTCTACCCCCTAGGGACAAGGAAGAAGCa                                    >  2:1021127/1‑141 (MQ=255)
                                                                                                  ataataAATGAACGATATAAAATTGTAGATAAGCTTGGCGCCGGTGGCATGAGTACCGTTTATCTTGCTGAAGATACGATACTTAACATTAAAGTTGCAATTAAGGCGATTTTTATACCACCTAGAGAAAAAGAAGAAACa                                    <  2:959852/141‑1 (MQ=255)
                                                                                                     ataaATGAACGATATAAAATTGTAGATAAGCTTGGCGGCGGTGGCATGAGTACCGTTTATCTTGCTGAAGATACGATACTTAACATTAAAGTTGCAATTAAGGCGATTTTTATACCACCTAGAGAAAAAGAAGAAACATTa                                 <  1:1193255/141‑1 (MQ=255)
                                                                                                      taaATGAACGATATAAGATTGTAGATAAGCTTGGCGCCGGTGGCATGAGTACCGTTTATCTTGCTGAAGATACGATACTTAACATTAAAGTTGCAATTAAGGCGATTTTTATCCCACCTAGCGAAAAAGAAGAAACATTaa                                >  1:713592/1‑141 (MQ=255)
                                                                                                      taaATGAACGATATAAAATTGTAGATAAGCTTGGCGGCGGTGGCATGAGTACCGTTTATCTTGCTGAAGATACGATACTTAACATTAAAGTTGCAATTAAGGCGATTTTTATACCACCTAGAGAAAAAGAAGAAACATTaa                                >  1:587542/1‑141 (MQ=255)
                                                                                                      taaATGAACGATATAAAATTGTAGATAAGCTTGGCGCCGGTGGCATGAGTACCGTTTATCTTGCTGAAGATACGATACTTAACATTAAAGTTGCAATTAAGGCGATTTTTATACCACCTAGAGAAAAAGAAGAAACATTaa                                >  1:990024/1‑141 (MQ=255)
                                                                                                       aaaTGAACGATATAAGATTGTAGATAAGCTTGGCGGCGGTGGCATGAGTACCGTTTATCTTGCTGAAGATACGATGCTTAACATTAAAGTTGCAATTAAGGCGATTTTTATACCACCTAGAGAAAAAGAAGAAACGTTaaa                               >  2:642976/1‑141 (MQ=255)
                                                                                                        aacgAACGATATAAAATTGTAGATAAACTTAGAGCCGGTGGCATGAGTACCGTATATCTTGAAAAAGATACAATACTTAACATTAAAATTGCAATTAAGGCTATTTTTATACCACCTAGAGAAAAAGAAGAAACATTaaaa                              <  2:235734/138‑1 (MQ=255)
                                                                                                        aaTGAACGATATAAAATTGTAGATAAGCTTGGCGCCGGTGGCATGAGTACCGTTTATCTTGCTGAAGATACGATACTTAACATTAAAGTTGCAATTAAGGCGATTTTTATACCACCTAGAGAAAAAGAAGAAACATTaaaa                              <  1:109796/141‑1 (MQ=255)
                                                                                                           gAACGATATAAAATTGTAGATAAGCTTGGCGCCGGTGGCATGAGTACCGTTTATCTTGCTGAAGATACGATACTTAACATTAAAGTTGCAATTAAGGCGATTTTTATACCACCTAGAGAAAAAGAAGAAACATTAAAACGt                           <  1:930763/141‑1 (MQ=255)
                                                                                                             aCGATATAAAATTGTAGATAAGCTTGGCGCCGGTGGCATGAGTACCGTTTATCTTGCTGAAGATACGATACTTAACATTAAAGTTGCAATTAAGGCGATTTTTATACCAGATTGAGGAAATGAAGAAACATTAAAACGttt                         >  1:681154/1‑141 (MQ=255)
                                                                                                             aCGATATAAAATTGTAGATAAGCTTGGCGCCGGTGGCATGAGTACCGTTTATCTTGCTGAAGATACGATACTTAACATTAAAGTTGCAATTAAGGCGATTTTTATACCACCTAGAGAAAAAGAAGAAACATTAAAACGttt                         <  1:389196/141‑1 (MQ=255)
                                                                                                             aCGATATAAAATTGTAGATAAGCTTGGCGCCGGTGGCATGAGTACCGTTTATCTTGCTGAAGATACGATACTTAACATTAAAGTTGCAATTAAGGCGATTTTTATACCACCTAGAGAAAAAGAAGAAACATTAAAACGttt                         >  1:688946/1‑141 (MQ=255)
                                                                                                                 tataAAATTGTAGATAAGCTTGGCGCCGGTGGCATGAGTACCATTTATCTTGCTGAAGATACGATACTTAACATTAAAGTTGCAATTAAGGCGATTTTTATACCACCTAGAGAAAAAGAAGAAACATTAAAACGTTTTGaa                     >  2:530544/1‑141 (MQ=255)
                                                                                                                 tataAAATTGTAGATAAGCTCGGTTCCGGTGGCATGAGTACCGTTTATCTTGCTGAAGATACGATACTTAACATTAAAGTTGCAATTAAGGCGATTTTTATACCACCTAGAGAAAAAGAAGAAACATTAAAACGTTTTGaa                     >  1:1189233/1‑141 (MQ=255)
                                                                                                                  ataAAATTGTAGATAAGCTTGGCGCCGGTGGCATGAGTACCGTTTATCTTGCTGAAGATACGATACTTAACATTAAAGTTGCAATTAAGGCGATTTTTATACCACCTAGAGAAAAAGAAGAAACATTAAAACGTTTTGAAc                    <  2:439071/141‑1 (MQ=255)
                                                                                                                    aaaaTTGTAGATAAGCTTGGCGCCGGTGGCATGAGTACCGTTTATCTTGCTGAAGATACGATACTTAACATTAAAGTTGCAATTAAGGCGATTTTTATACCACCTAGAGAAAAAGAAGAAACATTAAAACGTTTTGAACga                  <  2:154009/141‑1 (MQ=255)
                                                                                                                    aaaaTGGTAGATAAGCTTGGCGGCGGTGGCATGGGTACCATTTATCTTGCTGAAGATACGATACTTAACATTAAAGTTGCAATTAAGGCGATTGTTATACCACCTAGAGAAAAAGAAGAAACATTAAAACGTTTTGAACga                  <  2:471553/141‑1 (MQ=255)
                                                                                                                     aaaTTGTAGATAAGCTTGGGCCCGGTGCCATGAGTACCGTTTATCTTGCTGAAGATACGATACTTAACATTAAAGTTGCAATTAAGGCGATTTTTATACCACCTAGAGAAAAAGAAGAAACATTAAAACGTTTTGAACgag                 <  2:1022297/141‑1 (MQ=255)
                                                                                                                     aaaTTGTAGATAAGCTTGGCGGCGGTGGCATGAGTACCGTTTATCTTGCTGAAGATACGATACTTAACATTAAAGTTGCAATTAAGGCGATTTTTATACCACCTAGAGAAAAAGAAGAAACATTAAAACGTTTTGAACgag                 >  1:402639/1‑141 (MQ=255)
                                                                                                                     aaaTTGTAGATAAGCTTGGCGGCGGTGGCATGAGTACCGTTTATCTTGCTGAAGATACGATACTTAACATTAAAGTTGCAATTAAGGCGATTTTTATACCACCTAGAGAAAAAGAAGAAACATTAAAACGTTTTGAACgag                 >  1:363870/1‑141 (MQ=255)
                                                                                                                     aaaTTGTAGATAAGCTTGGCGGCGGTGGCATGAGTACCGTTTATCTTGCTGAAGATACGATACTTAACATTAAAGTTGCAATTAAGGCGATTTTTATACCACCTAGAGAAAAAGAAGAAACATTAAAACGTTTTGAACgag                 >  1:672235/1‑141 (MQ=255)
                                                                                                                       aTTGTAGATAAGCTTGGCGCCGGTGGCATGAGTACCGTTTATCTTGCTGAAGATACGATACTTAACATTAAAGTTGCAATTAAGGCGATTTTTATACCACCTAGAGAAAAAGAAGAAACATTAAAACGTTTTGAACGAGaa               >  2:645705/1‑141 (MQ=255)
                                                                                                                           tAGATAAGCTTGGCGCCGGTGGCATGAGTACCGTTTATCTTGCTGAAGATACGATACTTAACATTAAAGTTGCAATTAAGGCGATTTTTATACCACCTAGAGAAAAAGAAGAAACATTAAAACGTTTTGAACGAGAAGTAc           <  1:1094528/141‑1 (MQ=255)
                                                                                                                             gATAAGCTTGGCGGCGGTGGCATGAGTACCGTTTATCTTGCTGAAGATACGATACTTAACATTAAAGTTGCAATTAAGGCGATTTTTATACCACCTAGAGAAAAAGAAGAAACATTAAAACGTTTTGAACGAGAAGTACAt         <  1:918386/141‑1 (MQ=255)
                                                                                                                                aaGCTGGGCGCCGGTGGCATGAGTACCGTTTATCTTGCTGAAGATACGATACTTAACATTAAAGTTGCAATTAAGGCGATTTTTATACCACCTAGAGAAAAAGAAGAAACATTAAAACGTTTTGAACGAGAAGTACATAAc      <  2:675527/141‑1 (MQ=255)
                                                                                                                                   cTTGGCGCCGGTGGCATGAGTACCGTTTATCTTGCTGAAGATACGATACTTAACATTAAAGTTGCAATTAAGGCGATTTTTATACCACCTAGAGAAAAAGAAGAAACATTAAAACGTTTTGAACGAGAAGTACATAACtca   <  1:1200118/141‑1 (MQ=255)
                                                                                                                                    ttGGCGGCGGTGGCATGAGTACCGTTTATCTTGCTGAAGATACGATACTTAACATTAAAGTTGCAATTAAGGCGATTTTTATACCACCTAGAGAAAAAGAAGAAACATTAAAACGTTTTGAACGAGAAGTACATAACtcat  >  1:1138492/1‑141 (MQ=255)
                                                                                                                                    taggcgccTGTGGCATGAGTACCGTTTATCTTGCTGAAGATACGATACGTAACATTAAAGTTGCAATTAAGGCGATTTTTATACCACCTAGAGAAAAAGAAGAAACATTAAAACGTTTTGAACGAGAAGTACATAACtcat  <  1:1133544/139‑1 (MQ=255)
                                                                                                                                          |                                                                                                                                      
GGTGATCAATTAATGCAATTGGCATTAGATAACCATTCGAAAGATAACGTTACTTTCATACTCGCGGCTATTGAAGGTGATAAAGTATGATAGGTAAAATAATAAATGAACGATATAAAATTGTAGATAAGCTTGGCGGCGGTGGCATGAGTACCGTTTATCTTGCTGAAGATACGATACTTAACATTAAAGTTGCAATTAAGGCGATTTTTATACCACCTAGAGAAAAAGAAGAAACATTAAAACGTTTTGAACGAGAAGTACATAACTCAT  >  USA300TCH1516_ALE/1232010‑1232282

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