Predicted mutation
evidence seq id position mutation annotation gene description
RA USA300TCH1516_ALE 1,232,148 G→C 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,1480GC98.4% 189.0 / NA 63G18A (GGC→GCC) prkCSerine/threonine‑protein kinase PrkC
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base C (31/31);  total (31/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGTGATCAATTAATGCAATTGGCATTAGATAACCATTCGAAAGATAACGTTACTTTCATACTCGCGGCTATTGAAGGTGATAAAGTATGATAGGTAAAATAATAAATGAACGATATAAAATTGTAGATAAGCTTGGCGGCGGTGGCATGAGTACCGTTTATCTTGCTGAAGATACGATACTTAACATTAAAGTTGCAATTAAGGCGATTTTTATACCACCTAGAGAAAAAGAAGAAACATTAAAACGTTTTGAACGAGAAGTACATAACTCATCAC  >  USA300TCH1516_ALE/1232010‑1232285
                                                                                                                                          |                                                                                                                                         
ggTGATCAATTAATGCAATTGGCATTAGATAACCATTCGAAAGATAACGTTACTTTCATACTCGCGGCTATTGAAGGTGATAAAGTATGATAGGTGAAATAATAAATGAACGATATAACATCGTAGGTAAGCTTGgcgccg                                                                                                                                         >  1:769059/1‑141 (MQ=255)
ggTGATCAATTAATGCAATTGGCATTAGATAACCATTCGAAAGATAACGTTACTTTCATACTCGCGGCTATTGAAGGTGATAAAGTATGATAGGTAAAATAATAAATGAACGATATAAAATTGTGGATAAGCTTGgcgccg                                                                                                                                         >  1:538266/1‑141 (MQ=255)
ggTGATCAATTAATGCAATTGGCATTAGATAACCATTCGAAAGATAACGTTACTTTCATACTCGCGGCTATTGAAGGTGATAAAGTATGATAGGTAAAATAATAAATGAACGATATAAAATTGTAGATAAGCTTGgcgccg                                                                                                                                         >  1:115395/1‑141 (MQ=255)
ggTGATCAATTAATGCAATTGGCATTAGATAACCATTCGAAAGATAACGTTACTTTCATACTCGCGGCTATTGAAGGTGATAAAGTATGATAGGTAAAATAATAAATGAACGATATAAAATTGTAGATAAGCTTGgcgccg                                                                                                                                         >  1:214636/1‑141 (MQ=255)
   gATCAATTAATGCAATTGGCATTAGATAACCATTCGAAAGATAACGTTACTTTCATACTCGCGGCTATTGAAGGTGATAAAGTATGATAGGTAAAATAATAAATGAACGATATAAAATTGTAGATAAGCTTGGCGCCGGTg                                                                                                                                      >  2:188520/1‑141 (MQ=255)
   gATCAATTAATGCAATTGGCATTAGATAACCATTCGAAAGATAACGTTACTTTCATACTCGCGGCTATTGAAGGTGATAAAGCACGAAACGTAAAATAATAAATGAACGATATAAAATTGTAGATAAGCTTGGCGCAGGTg                                                                                                                                      >  1:428793/1‑141 (MQ=255)
      cAATTAATGCAATTGGCATTAGATAACCATTCGAAAGATAACGTTACTATCATACTCGCGGCAATTGAAGGTGATAAAGTATGATAGGTAAACTAATAAATGAACGATATAAAAGTGTAGATAAGCTTGGCGCGGGTGGCa                                                                                                                                   >  2:277581/1‑141 (MQ=255)
         ttAATGCAATTGGCATTAGATAACCATTCGAAAGATAACGTTACTTTCATACTCGCGGCTATTGAAGGTGATAAAGTATGATAGGTAAAATAATAAATGAACGATATAAAGTTGTCGATACGCGTGGCGCCGGTGTCATgg                                                                                                                                >  2:488550/1‑140 (MQ=255)
                    ggaaTTAGATAACCATTCGAAAGATAACGTTACTTTCATACTCGCGGCTATTGAAGGTGATAAAGTATGATAGGTAAAATAATAAATGAACGATATAAAATTGTAGATAAGGTGTGCTCCGGGGGCATGAGTACCGTTTAt                                                                                                                     <  1:236519/138‑1 (MQ=255)
                    ggCTTTAGATAACCATTCGAAAGATAACGTTACTTTCATACTCGCGGCTATTGAAGGTGATAAAGTATGATAGGTAAAATAATAAATGAACGATATAAAATTGTAGATAAGCTTGGCGCCGGTGGCATGAGTACCGTTTAt                                                                                                                     <  1:289234/141‑1 (MQ=255)
                    ggCATTAGATAACCATTCGAAAGATAACGTTACTTTCATACTCGCGGCTATTGAAGGTGATAAAGTATGATAGGTAAAATAATAAATGAACGATATAAAATTGTAGATAAGCTTGGCGCCGGTGGCATGAGTACCGTTTAt                                                                                                                     <  1:832403/141‑1 (MQ=255)
                    ggCATTAGATAACCATTCGAAAGATAACGTTACTTTCATACTCGCGGCTATTGAAGGTGATAAAGTATGATAGGTAAAATAATAAATGAACGATATAAAATTGTAGATAAGCTTGGCGCCGGTGGCATGAGTACCGTTTAt                                                                                                                     <  1:155910/141‑1 (MQ=255)
                    ggCATTAGATAACCATTCGAAAGATAACGTTACTTTCATACTCGCGGCTATTGAAGGTGATAAAGTATGATAGGTAAAATAATAAATGAACGATATAAAATTGTAGATAAGCTTGGCGCCGGTGGCATGAGTACCGTTTAt                                                                                                                     <  1:488601/141‑1 (MQ=255)
                      gaTTGGATTACCATTGGAAATGGAACGGTACTTCCGTACTCGCGGCTATTGAAGGTGATAAAGTATGATAGGTAAAATAATAAATGAACGATTTAAACTTATACTTTAGGTTGGGGCCGGTGGCGTGAGTACCGCTTATCt                                                                                                                   >  2:205201/2‑141 (MQ=255)
                      cATTAGATAACCATTCGAAAGATAACGTTACTTTCATACTCGCGGCTATTGAAGGTGATAAAGTATGATAGGTAAAATAATAAATGAACGATATAAAATTGTAGATAAGCTTGGCGCCGGTGGCATGGGTACCGTTTATCt                                                                                                                   >  1:893264/1‑141 (MQ=255)
                              aaCCATTCGAAAGATAACGTTACTTTCATACTCGCGGCTATTGAAGGTGATAAAGTATGATAGGTAAAATAATAAATGAACGATATAAAATTGTAGATAAGCTTGGCGCCGGTGGCATGAGTACCGTTTATCTTGCTGAAg                                                                                                           <  2:183472/141‑1 (MQ=255)
                                ccATTCGAAAGATAACGTTACTTTCATACTCGCGGCTATTGAAGGTGATAAAGTATGATAGGTAAAATAATAAATGAAAGCTCAAAAATGGTAGATAAGCTTGGCGCCGGTGGCATGAGTACCGTTTATCTTGCTGAAgat                                                                                                         >  1:718111/1‑141 (MQ=255)
                                  aTTCGAAAGATAACGTTACTTTCATACTCGCGGCTATTGAAGGTGATAAAGTATGATAGGTAAAATAATAAATGAACGATATAAAATTGTAGATAAGCTTGGCGCCGGTGGCATGAGTACCGTTTATCTTGCTGAAgatac                                                                                                       >  2:686408/1‑141 (MQ=255)
                                  aTTCGAAAGATAACGTTACTTTCATACTCGCGGCTATTGAAGGTGATAAAGTATGATAGGTAAAATAATAAATGAACGATATAAAATTGTAGATAAGCTTGGCGCCGGTGGCATGAGTACCGTTTATCTTCGTGAAgaaac                                                                                                       >  1:339136/1‑141 (MQ=255)
                                   ttCGAAAGATAACGTTACTTTCATACTCGCGGCTATTGAAGGTGATAAAGTATGATAGGTAAAATAATAAATGAACGATATAAAATTGTAGATAAGCTTGGCGCCGGTGGCATGAGTACCGTTTATCTTGCTGAAgatacg                                                                                                      <  2:349763/141‑1 (MQ=255)
                                       aaaGATAACGTTACTTTCATACTCGCGGCTATTGAAGGTGATAAAGTATGATAGGTAAAATAATAAATGAACGATATAAAATAGTAGATAAGCTTGGCGCCGGTGGCATGAGTACCGTTTATCTTGCTGAAgatacgatac                                                                                                  >  2:761829/1‑141 (MQ=255)
                                        aaGATAACGTTACTTTCATACTCGCGGCTATTGAAGGTGATAAAGTATGATAGGTAAAATAATAAATGAACGATATAAAATTGTAGATAAGCTTGGCGCCGGTGGCATGAGTACCGTTTATCTTGCTGAAGATACGATACt                                                                                                 <  2:854792/141‑1 (MQ=255)
                                          gATAACGTTACTTTCATACTCGCGGCTATTGAAGGTGATAAAGTATGATAGGTAAAATAATAAATGAACGATATAAAATTGTAGATAAGCTTGGCGCCGGTGGCATGAGTGCCGTTTATCTTGCTGAAGATACGATACTTa                                                                                               >  2:889567/1‑141 (MQ=255)
                                            tAACGTTACTTTCATACTCGCGGCTATTGAAGGTGATAAAGTATGATAGGTAAAATAATAAATGAACGATATAAAATTGTAGATAAGCTTGGCGCCGGTGGCATGAGTACCGTTTATCTTGCTGAAGATACGATACTTAAc                                                                                             <  1:897120/141‑1 (MQ=255)
                                            aaaCGTTACTTTCATACTCGCGGCTATTGAAGGTGATAAAGTATGATAGGTAAAATAATAAATGAACGATATAAAATTGTAGATAAGCTTGGCGCCGGTGGCATGAGTACCGTTTATCTTGCTGAAGATACGATACTTAAc                                                                                             <  1:897148/140‑1 (MQ=255)
                                               cGTTACTTTCATACTCGCGGCTATTGAAGGTGATAAAGTATGATAGGTAAAATAATAAATGAACGATATAAAATTGTAGATAAGCTTGGCGCCGGTGGCATGAGTACCGTTTATCTTGCTGAAGATACGATACTTAACAtt                                                                                          >  1:128131/1‑141 (MQ=255)
                                                    cTTTCATACTCGCGGCTATTGAAGGTGATAAAGTATGATAGGTAAAATAATAAATGAACGATATAAAATTGTAGATAAGCTTGGCGCCGGTGGCATGAGTACCGTTTATCTTGCTGAAGATACGATACTTAACATTAAAGt                                                                                     <  2:719471/141‑1 (MQ=255)
                                                            cTCGCGGCTATTGAAGGTGATAAAGTATGATAGGTAAAATAATAAATGAACGATATAAAATTGTAGATAAGCTTGGCGCCGGTGGCATGAGTACCGTTTATCTTGCTGAAGATACGATACTTAACATTAAAGTTGCAATTa                                                                             <  1:138091/141‑1 (MQ=255)
                                                              cgcgGCTATTGAAGGTGATAAAGTATGATAGGTAAAATAATAAATGAACGATATAAAATTGTAGATAAGCTTGGCGCCGGTGGCATGAGTACCGTTTATCTTGCTGAAGATACGATACTTAACATTAAAGGTGCAATTAAg                                                                           >  1:273139/1‑141 (MQ=255)
                                                                  gCTATAGAAGGTGATAAAGTATGATAGGTAAAATAATAAATGAACGATATAAAATTGTAGATAAGCTTGGCGCCGGTGGCATGAGTACCGTTTATCTTGCTGAAGATACGATACTTAACATTAAAGTTGCAATTAAGGCGa                                                                       <  2:718111/141‑1 (MQ=255)
                                                                       tGAAGGTGATAAAGTATGATAGGTAAAATAATAAATGAACGATATAAAATTGTAGATAAGCTTGGCGCCGGTGGCATGAGTACCGTTTATCTTGCTGAAGATACGATACTTAACATTAAAGTTGCAATTAAGGCGAttttt                                                                  <  1:417218/141‑1 (MQ=255)
                                                                        gAAGGTGATAAAGTATGATAGGTAAAATAATAAATGAACGATATAAAATTGTAGATAAGCTTGGCGCCGGTGGCATGAGTACCGTTTATCTTGCTGAAGATACGATACTTAACATTAAAGTTGCAAGTAAGGCGATTTTta                                                                 <  1:538818/141‑1 (MQ=255)
                                                                           ggTGATAAAGTATGATAGGTAAAATAATAAATGAACGATATAAAATTGTAGATAAGCTAGGCGCCGGAGGCATGAGTACCGATTATCTTGCTGAAGATACGAAACTTAACATAAAAGTTGCAATTAAGGCGAGTTTTATac                                                              >  2:255745/1‑141 (MQ=255)
                                                                            gTGATAAAGTATGATAGGTAAAATAATAAATGAACGATATAAAATTGTAGATAAGCTTGGCGCCGGTGGCATGAGTAACGATTATCTTGCTGCAGATACGATACTTAACATTAAAGGTGCAATTAAGGCGTTTCTTATaca                                                             >  2:638289/1‑140 (MQ=255)
                                                                              gATAAAGTATGATAGGTAAAATAATAAATGAACGATATAAAATTGTAGATAAGCTTGGCGCCGGTGGCATGAGTACCGTTTATCTTGCTGAAGATACGATACTTAACATTAAAGTTGCAATTAAGGCGATTTTTATaccac                                                           <  2:616313/141‑1 (MQ=255)
                                                                              gATAAAGTATGATAGGTAAAATAATAAATGAACGATATAAAATTGTAGATAAGCTTGGCGCCGGTGGCATGAGTACCGTTTATCTTGCTGAAGATACGATACTTAACATTAAAGTTGCAATTAAGGCGATTTTTATaccac                                                           <  2:416473/141‑1 (MQ=255)
                                                                                 aatgTATGATAGGTAAAATAATCACTGAACGACATAAATTTGTAGATTAGCTTGGCGCAGTTGGCAGGAGTACCGTTTATCTGGCTGAAGATACGATACTTAACATTAAAGTTGCAATTAAGGCGATTTTTATACCACCTa                                                        <  2:85805/138‑1 (MQ=255)
                                                                                       tGATAGGTAAAATAATAAATGAACGATATAAAATTGTAGATAAGCTTGGCGCCGGTGGCATGAGTACCGTTTATCTTGCTGAAGATACGATACTTAACATTAAAGTTGCAATTAAGGCGATTTTTATACCACCTAGAGaaa                                                  <  1:853540/141‑1 (MQ=255)
                                                                                       tGATAGGTAAAATAATAAATGAACGATATAAAATTGTAGATAAGCTTGGCGCCGGTGGCATGAGTACCGTTTATCTTGCTGAAGATACAATACTTAACATTAAAGTTGCAATTAAGGCGATTTTTATACCACCTAGAGaaa                                                  <  1:650222/141‑1 (MQ=255)
                                                                                       tGATAGGTAAAATAAAAAAAGAACGATATAAAATTGTAGATAAGCTTGGCGCCGGTGGCATGAGTACCGTTTATCTTGCTGAAGATACGATACTTAACATTAAAGTTGCAATTAAGGCGATTTTTATACCACCTAGAGaaa                                                  <  1:608546/141‑1 (MQ=255)
                                                                                        gATAGGTAAAATAATAAATGAACGATATAAAATTGTAGATAAGCTTGGCGCCGGTGGCATGAGTACCGTTTATCTTGCTGAGGATACGATACTTAACATTAAAGTTGCAATTAAGGCGATTTTTATACCACCTAGAGaaaa                                                 >  2:921436/1‑141 (MQ=255)
                                                                                        gATAGGTAAAATAATAAATGAACGATATAAAATTGTAGATAAGCTTGGCGCCGGTGGCATGAGTACCGTTTATCTTGCTGAAGATACGATACTTAACATTAAAGTTGCAATTAAGGCGATTTTTATACCACCTAGAGaaaa                                                 >  1:935116/1‑141 (MQ=255)
                                                                                          tAGGTAAAATAATAAATGAACGATATAAAATTGTAGATAAGCTTGGCGCCGGTGGCATGAGTACCGTTTATCTTGCTGAAGATACGATACTTAACATTAAAGTTGCAATTAAGGCGATTTTTATACCACCTAGAGACAaag                                               >  1:487339/1‑141 (MQ=255)
                                                                                          tAGGTAAAATAATAAATGAACGATATAAAATTGTAGATAAGCTTGGCGCCGGTGGCATGAGTACCGTTTATCTTGCTGAAGATACGATACTTAACATTAAAGTTGCAATTAAGGCGATTTTTATACCACCTAGAGAAAaag                                               >  1:910216/1‑141 (MQ=255)
                                                                                            ggTAAAATAATAAATGAACGATATAAAATTGTAGATAAGCTTGGCGCCGGTGGCATGAGTACCGTTTATCTTGCTGAAGATACGATACTTAACATTAAAGTTGCAATTAAGGCGATTTTTATACCACCTAGAGAAAaagaa                                             >  1:359392/1‑141 (MQ=255)
                                                                                            ggTAAAATAATAAATGAACGATATAAAATTGTAGATAAGCTTGGCGCCGGTGGCATGAGTACCGTTTATCTTGCTGAAGATACGATACTTAACATTAAAGTTGCAATTAAGGCGATTTTTATACCACCTAGAGAAAaagaa                                             >  1:239017/1‑141 (MQ=255)
                                                                                             gTAAAATAATAAATGAACGATATAAAATTGTAGATAAGCTTGGCGCCGGTGGCATGAGTACCGTTTATCTTGCTGAAGATACGATACTTAACATTAAAGTTGCAATTAAGGCGATTTTTATACCACCTAGAGAAAaagaag                                            >  1:154393/1‑141 (MQ=255)
                                                                                              tAAAATAATAAATGAACGATATAAAATTGTAGCTAAGCTTGGCGCCGGTGGCATGAGTACCGTTTATCTTGCTGAAGATACGATACTCAACATTAAAGTTGCAATTAAGGCGATTTTTATACCACCTAGAGAAAAagaaga                                           <  1:663208/141‑1 (MQ=255)
                                                                                                   aaataaATGAACGATATAAAATTGTAGATAAGCTTGGCGCCGGTGGCATGAGTACCGTTTATCTTGCTGAAGATACGATACTTAACATTAAAGTTGCAATTAAGGCGATTTTTATACCACCTAGAGAAAAAGAAGAAACAt                                      <  1:834110/140‑1 (MQ=255)
                                                                                                           gAACGATATAAAATTGTAGATAAGCTTGGCGCCGGTGGCATGAGTACCGTTTATCTTGCTGAAGATACGATACTTAACATTAAAGTTGCAATTAAGGCGATTTTTATACCACCTAGAGAAAAAGAAGAAACATTAAAACGt                              <  2:397675/141‑1 (MQ=255)
                                                                                                           gAACGATATAAAATTGTAGATAAGCTTGGCGCCGGTGGCATGAGTACCGTTTATCTTGCTGAAGATACGATACTTAACATTAAAGTTGCAATTAAGGCGATTTTTATACCACCTAGAGAAAAAGAAGAAACATTAAAACGt                              <  1:31415/141‑1 (MQ=255)
                                                                                                           gAACGATATAAAATTGTAGATAAGCTTGGCGCCGGTGGCATGAGTACCGTTTATCTTGCTGAAGATACGATACTTAACATTAAAGTTGCAATTAAGGCGAGTTTTATACCACCTAGAGAAAAAGAAGAAACATTAAAACGt                              <  2:699782/141‑1 (MQ=255)
                                                                                                           gAACGATATAAAAGTGTAGATAAGCTTGGCTGCGGTGGCATGAGTACCGTTTATCATGCTGAAGATACGATACTTAAAAATAAAGTTGCAATTTCGGCGAGTTTTATACTACCTGGAGAAATAGAAGAAACATGAAAACGt                              <  2:683571/141‑1 (MQ=255)
                                                                                                             aCGATATAAAATTGTAGATAAGCTTGGCGCCGGTGGCATGAGTACCGTTTATCTTGCTGAAGATACGATACTTAACATTAAAGTTGCAATTAAGGCGATTTTTATACCACCTAGAGAAAAAGAAGAAACATTAAAActtct                            >  1:347569/1‑137 (MQ=255)
                                                                                                             aCGATATAAAATTGTAGATAAGCTTGGCGCCGGTGGCATGAGTACCGTTTATCTTGCTGAAGATACGATACTTAACATTAAAGTTGCAATTAAGGCGATTTTTATACCACCTAGAGAAAAAGAAGAAACATTAAAACGttt                            >  1:793102/1‑141 (MQ=255)
                                                                                                               gATATAAAATTGTAGATAAGCTTGGCGCCGGTGGCATGAGTACCGTTTATCTTGCTGAAGATACGATACTTAACATTAAAGTTGCAATTAAGGCGATTTTTATACCACCTAGAGAAAAAGAAGAAACATTAAAACGTTTTg                          >  2:500497/1‑141 (MQ=255)
                                                                                                                atatAAAATTGTAGATAAGCTTGGCGCCGGTGGCATGAGTACCGTTTATCTTGCTGAAGATACGATACTTAACATTAAAGTTGCAATTAAGGCGATTTTTATACCACCTAGAGAAAAAGAAGAAACATTAAAACGTTTTGa                         <  1:627924/141‑1 (MQ=255)
                                                                                                                   tgggcttGTAGATAGGCTTTGCGCCTGTGGCATGAGTACCGTTTATCTTGCTGAAGATACGATACTTAACATTAAAGTTGCAATTAAGGCGATTTTTATACCACCTAGAGAAAAAGAAGAAACATTAAAACGTTTTGAACg                      <  2:988551/136‑1 (MQ=255)
                                                                                                                   taAAATTGTAGATAAGCTTGGCGCCGGTGGCATGAGTACCGTTTATCTTGCTGAAGATACGATACTTAACATTAAAGTTGCAATTAAGGCGATTTTTATACCACCTAGAGAAAAAGAAGAAACATTAAAACGTTTTGACCg                      >  2:67581/1‑141 (MQ=255)
                                                                                                                    aaaaTTGTAGATAAGCTTGGCGCTGGTGGCTGGAATACCGTTTATCTTGCTGAAGATACGATACTTAACATTAAAGTTGCAATTAAGGCGATTTTTATACCACCTAGAGAAAAAGAAGAAACATTAAAACGTTTTGAACga                     >  2:908021/1‑141 (MQ=255)
                                                                                                                         tGTAGATAAGCGTGGCGCCGGTGGCATGAGTACCGTTTATCTTGCTGAAGATACGATACTTAACATTAAAGTTGCAATTAAGGCGATTTTTATACCACCTAGAGAAAAAGAAGAAACATTAAAACGTTTTGAACGAGAAGt                <  2:215041/141‑1 (MQ=255)
                                                                                                                             gATAAGCTTGGCGCCGGTGGCATGAGTACCGTTTATCTTGCTGAAGATACGATACTTAACGTTAAAGTTGCAATTAAGGCGATTTTTATACCACCTAGAGAAAAAGAAGAAACATTATAACGTTTTGAACGAGAAGTACAt            <  2:707936/141‑1 (MQ=255)
                                                                                                                             gATAAGCTTGGCGCCGGTGGCATGAGTACCGTGTATCTTGCTGAAGATACGATACTTAACATTAAAGTTGCAATTAAGGCGATTTTTATACCACCTAGAGAAAAAGAAGAAACATTAAAACGTTTTGAACGAGAAGTACAt            <  2:673066/141‑1 (MQ=255)
                                                                                                                              aTAAGCTTGGCGCCGGTGGCATGAGTACCGTTTATCTTGCTGAAGATACGATACTTAACATTAAAGTTGCAATTAAGGCGATTTTTATACCACCTAGAGAAAAAGAAGAAACATTAAAACGCTTTGAACGAGAAGTACATa           >  2:815118/1‑141 (MQ=255)
                                                                                                                                aaGCTTGGCGCCGGTGGCATGAGTACCGTTTATCTTGCTGAAGATACGATACTTAACATTAAAGTTGCAATTAAGGCGATTTTTATACCACCTAGAGAAAAAGAAGAAACATTAAAACGTTTTGAACGAGAAGTACATAAc         >  1:262180/1‑141 (MQ=255)
                                                                                                                                aaGCTTGGCGCCGGTGGCATGAGTACCGTTTATCTTGCTGAAGATACGATACTTAACATTAAAGTTGCAATTAAGGCGATTTTTATACCACCTAGAGAAAAAGAAGAAACATTAAAACGTTTTGAACGAGAAGTACATAAc         >  1:262321/1‑141 (MQ=255)
                                                                                                                                aaGCTTGGCGCCGGTGGCATGAGTACCGTTTATCTTGCTGAAGATACGATACTTAACATTAAAGTTGCAATTAAGGAGATTTTTATACCACCTAGAGAAAAAGCAGAAGCATGAAAACGTTTTGCACGAGAGGAACATAAc         >  1:316787/1‑141 (MQ=255)
                                                                                                                                      ggcgccGGTGGCATGAGTACCGTTTATCTTGCTGAAGATACGATACTTAACATTAAAGTTGCAATTAAGGCGATTTTTATACCACCTAGAGAAAAAGAAGAAACATTAAAACGTTTGGAACGAGAAGTGCATAACgcatca   >  2:341260/1‑141 (MQ=255)
                                                                                                                                       gcgccgGTGGCATGAGTACCGTTTATCTTGCTGAAGATACGATACTTAACATTAAAGTTGCAATTAAGGCGATTTGTATACCACCTAGAGAAAAAGAAGAAAAATTAACACGTCGTGAAGGAGAAACAGATAACTCATcac  >  2:470622/1‑141 (MQ=255)
                                                                                                                                          |                                                                                                                                         
GGTGATCAATTAATGCAATTGGCATTAGATAACCATTCGAAAGATAACGTTACTTTCATACTCGCGGCTATTGAAGGTGATAAAGTATGATAGGTAAAATAATAAATGAACGATATAAAATTGTAGATAAGCTTGGCGGCGGTGGCATGAGTACCGTTTATCTTGCTGAAGATACGATACTTAACATTAAAGTTGCAATTAAGGCGATTTTTATACCACCTAGAGAAAAAGAAGAAACATTAAAACGTTTTGAACGAGAAGTACATAACTCATCAC  >  USA300TCH1516_ALE/1232010‑1232285

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: