Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 2,454,626 A→T 44.3% L48L (CTA→CTT lyrA → Lysostaphin resistance protein A

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE2,454,6260AT44.3% 40.1 / 161.6 106L48L (CTA→CTTlyrALysostaphin resistance protein A
Reads supporting (aligned to +/- strand):  ref base A (35/24);  new base T (22/25);  total (57/49)
Fisher's exact test for biased strand distribution p-value = 2.41e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.10e-01

Only 100 of 107 total aligned reads displayed.

TTGTCATTACAATGGCGTTATCCATTATGCTCAGAGATTTCCAGTCTATAATTGGTGTCAAACACTTTATATTTGAAGTTACAGATCTAGCACCATTAATTGCTGCAATCATTTGTATACTCGTTTTCAAATATAAAAAGGTCCAACTTGCAGGTTTAAAATTCTCAATCAGCCTGAAA  >  USA300TCH1516_ALE/2454538‑2454716
                                                                                        |                                                                                           
ttGTCATTACAATGGCGTTATCCATTATGCTCAGAGATTTCCAGTCTATAATTGGTGTCAAACACTTTATATTTGAAGTTACAGATCTAGc                                                                                           <  1:3420234/91‑1 (MQ=255)
ttGTCATTACAATGGCGTTATCCATTATGCTCAGAGATTTCCAGTCTATAATTGGTGTCAAACACTTTATATTTGAAGTTACAGATCTAGc                                                                                           <  1:1731718/91‑1 (MQ=255)
  gTCATTACAATGGCGTTATCCATTATGCTCAGAGATTTCCAGTCTATAATTGGTGTCAAACACTTTATATTTGAAGTTACAGATCTAGCAc                                                                                         <  1:2394954/91‑1 (MQ=255)
  gTCATTACAATGGCGTTATCCATTATGCTCAGAGATTTCCAGTCTATAATTGGTGTCAAACACTTTATATTTGAAGTTACAGATCTAGCAc                                                                                         <  1:2394849/91‑1 (MQ=255)
  gTCATTACAATGGCGTTATCCATTATGCTCAGAGATTTCCAGTCTATAATTGGTGTCAAACACTTTATATTTGAAGTTACAGATCTAGCAc                                                                                         >  1:3710730/1‑91 (MQ=255)
   tCATTACAATGGCGTTATCCATTATGCTCAGAGATTTCCAGTCTATAATTGGTGTCAAACACTTTATATTTGAAGTTACAGATCTTGCAcc                                                                                        >  1:1204596/1‑91 (MQ=255)
    cATTACAATGGCGTTATCCATTATGCTCAGAGATTTCCAGTCTATAATTGGTGTCAAACACTTTATATTTGAAGTTACAGATCTTGCACCa                                                                                      <  1:785920/91‑1 (MQ=255)
    cATTACAATGGCGTTATCCATTATGCTCAGAGATTTCCAGTCTATAATTGGTGTCAAACACTTTATATTTGAAGTTACAGATCTTGCACCa                                                                                      >  1:397449/1‑91 (MQ=255)
    cATTACAATGGCGTTATCCATTATGCTCAGAGATTTCCAGTCTATAATTGGTGTCAAACACTTTATATTTGAAGTTACAGATCTTGCACCa                                                                                      <  1:3124731/91‑1 (MQ=255)
    cATTACAATGGCGTTATCCATTATGCTCAGAGATTTCCAGTCTATAATTGGTGTCAAACACTTTATATTTGAAGTTACAGATCTTGCACCa                                                                                      <  1:1273309/91‑1 (MQ=255)
    cATTACAATGGCGTTATCCATTATGCTCAGAGATTTCCAGTCTACAATTGGTGTCAAACACTTTATATTTGAAGTTACAGATCTTGCACCa                                                                                      <  1:3124869/91‑1 (MQ=255)
     aTTACAATGGCGTTATCCATTATGCTCAGAGATTTCCAGTCTATAATTGGTGTCAAACACTTTATATTTGAAGTTACAGATCTTGCACCAt                                                                                     <  1:3671624/91‑1 (MQ=255)
     aTTACAATGGCGTTATCCATTATGCTCAGAGATTTCCAGTCTATAATTGGTGTCAAACACTTTATATTTGAAGTTACAGATCTAGCACCAt                                                                                     <  1:167100/91‑1 (MQ=255)
     aTTACAATGGCGTTATCCATTATGCTCAGAGATTTCCAGTCTATAATTGGTGTCAAACACTTTATATTTGAAGTTACAGATCTAGCACCAt                                                                                     <  1:712256/91‑1 (MQ=255)
      ttACAATGGCGTTATCCATTATGCTCAGAGATTTCCAGTCTATAATTGGTGTCAAACACTTTATATTTGAAGTTACAGATCTAGCACCact                                                                                    >  1:523788/1‑89 (MQ=255)
      ttACAATGGCGTTATCCATTATGCTCAGAGATTTCCAGTCTATAATTGGTGTCAAACACTTTATATTTGAAGTTACAGATCTAGCACCAtt                                                                                    >  1:523890/1‑91 (MQ=255)
       tACAATGGCGTTATCCATTATGCTCAGAGATTTCCAGTCTATAATTGGTGTCAAACACTTTATATTTGAAGTTACAGATCTTGCACCATTa                                                                                   >  1:659458/1‑91 (MQ=255)
       tACAATGGCGTTATCCATTATGCTCAGAGATTTCCAGTCTATAATTGGTGTCAAACACTTTATATTTGAAGTTACAGATCTTGCACCATTa                                                                                   >  1:660260/1‑91 (MQ=255)
       tACAATGGCGTTATCCATTATGCTCAGAGATTTCCAGTCTATAATTGGTGTCAAACACTTTATATTTGAAGTTACAGATCTTGCACCATTa                                                                                   >  1:660304/1‑91 (MQ=255)
         cAATGGCGTTATCCATTATGCTCAGAGATTTCCAGTCTATAATTGGTGTCAAACACTTTATATTTGAAGTTACAGATCTAGCACCATTAAt                                                                                 <  1:326768/91‑1 (MQ=255)
             ggCGTTATCCATTATGCTCAGAGATTTCCAGTCTATAATTGGTGTCAAACACTTTATATTTGAAGTTACAGATCTTGCACCATTAATtgct                                                                             <  1:2072778/91‑1 (MQ=255)
             ggCGTTATCCATTATGCTCAGAGATTTCCAGGCTATAATTGGTGTCAAACACTTTATATTTGAAGTTACAGATCTTGCAGCAATTAAtggt                                                                              <  1:2033005/91‑4 (MQ=255)
             ggCGTTAACCATTATGCTCAGAGATTTCCAGTCTATAATTGGTGTCAAACACTTTATATTTGAAGTTACAGATCTTGCACCATTAATtgct                                                                             <  1:2072682/91‑1 (MQ=255)
              gCGTTATCCATTATGCTCAGAGATTTCCAGTCTATAATTGGTGTCAAACACTTTATATTTGAAGTTACAGATCTTGCACCATTAATtgctg                                                                            <  1:688650/91‑1 (MQ=255)
              gCGTTATCCATTATGCTCAGAGATTTCCAGTCTATAATTGGTGTCAAACACTTTATATTTGAAGTTACAGATCTAGCACCATTAATtgctg                                                                            <  1:108660/91‑1 (MQ=255)
                  tATCCATTATGCTCAGAGATTTCCAGTCTATAATTGGTGTCAAACACTTTATATTTGAAGTTACAGATCTTGCACCATTAATTGCTGCAAt                                                                        >  1:2423178/1‑91 (MQ=255)
                  tATCCATTATGCTCAGAGATTTCCAGTCTATAATTGGTGTCAAACACTTTATATTTGAAGTTACAGATCTTGCACCATTAATTGCTGCAAt                                                                        >  1:12014/1‑91 (MQ=255)
                   aTCCATTATGCTCAGAGATTTCCAGTCTATAATTGGTGTCAAACACTTTATATTTGAAGTTACAGATCTTGCACCATTAATTGCTGCAATc                                                                       <  1:2865209/91‑1 (MQ=255)
                   aTCCATTATGCTCAGAGATTTCCAGTCTATAATTGGTGTCAAACACTTTATATTTGAAGTTACAGATCTTGCACCATTAATTGCTGCAATc                                                                       >  1:1054176/1‑91 (MQ=255)
                   aTCCATTATGCTCAGAGATTTCCAGTCTATAATTGGTGTCAAACACTTTATATTTGAAGTTACAGATCTAGCACCATTAATTGCTGCAATc                                                                       >  1:1506351/1‑91 (MQ=255)
                    tCCATTATGCTCAGAGATTTCCAGTCTATAATTGGTGTCAAACACTTTATATTTGAAGTTACAGATCTAGCACCATTAATTGCTGCAATCa                                                                      >  1:1879155/1‑91 (MQ=255)
                    tCCATTATGCTCAGAGATTTCCAGTCTATAATTGGTGTCAAACACTTTATATTTGAAGTTACAGATCTAGCACCATTAATTGCTGCAATCa                                                                      >  1:2645200/1‑91 (MQ=255)
                       aTTATGCTCAGAGATTTCCAGTCTATAATTGGTGTCAAACACTTTATATTTGAAGTTACAGATCTAGCACCATTAATTGCTGCAATCAttt                                                                   >  1:2637477/1‑91 (MQ=255)
                        ttATGCTCAGAGATTTCCAGTCTATAATTGGTGTCAAACACTTTATATTTGAAGTTACAGATCTAGCACCATTAATTGCTGCAATCATTTg                                                                  >  1:1945480/1‑91 (MQ=255)
                          aTGCTCAGAGATTTCCAGTCTATAATTGGTGTCAAACACTTTATATTTGAAGTTACAGATCTTGCACCATTAATTGCTGCAATCATTTGta                                                                <  1:2424672/91‑1 (MQ=255)
                            gCTCAGAGATTTCCAGTCTATAATTGGTGTCAAACACTTTATATTTGAAGTTACAGATCTAGCACCATTAATTGCTGCAATCATTTGtata                                                              <  1:3756726/91‑1 (MQ=255)
                            gCTCAGAGATTTCCAGTCTATAATTGGTGTCAAACACTTTATATTTGAAGTTACAGATCTAGCACCATTAATTGCTGCAATCATTTGtata                                                              >  1:1338067/1‑91 (MQ=255)
                            gCTCAGAGATTTCCAGTCTATAATTGGTGTCAAACACTTTATATTTGAAGTTACAGATCTAGCACCATTAATTGCTGCAATCATTTGtata                                                              >  1:1337493/1‑91 (MQ=255)
                             cTCAGAGATTTCCAGTCTATAATTGGTGTCAAACACTTTATATTTGAAGTTACAGATCTAGCACCATTAATTGCTGCAATCATTTGTATAc                                                             <  1:1949904/91‑1 (MQ=255)
                              tCAGAGATTTCCAGTCTATAATTGGTGTCAAACACTTTATATTTGAAGTTACAGATCTTGCACCATTAATTGCTGCAATCATTTGTATACt                                                            <  1:2554760/91‑1 (MQ=255)
                               cAGAGATTTCCAGTCTATAATTGGTGTCAAACACTTTATATTTGAAGTTACAGATCTTGCACCATTAATTGCTGCAATCATTTGTATACTc                                                           <  1:622127/91‑1 (MQ=255)
                               cAGAGATTTCCAGTCTATAATTGGTGTCAAACACTTTATATTTGAAGTTACAGATCTTGCACCATTAATTGCTGCAATCATTTGTATACTc                                                           >  1:809100/1‑91 (MQ=255)
                               cAGAGATTTCCAGTCTATAATTGGTGTCAAACACTTTATATTTGAAGTTACAGATCTAGCACCATTAATTGCTGCAATCATTTGTATACTc                                                           >  1:3643741/1‑91 (MQ=255)
                               cAGAGATTTCCAGTCTATAATTGGTGTCAAACACTTTATATTTGAAGTTACAGATCTAGCACCATTAATTGCTGCAATCATTTGTATACTc                                                           <  1:3694221/91‑1 (MQ=255)
                               cAGAGATTTCCAGTCTATAATTGGTGTCAAACACTTTATATTTGAAGCTACAGATCTAGCACCATTAATTGCTGCAACCATTTGTATACTc                                                           >  1:4145778/1‑91 (MQ=255)
                                agagATTTCCAGTCTATAATTGGTGTCAAACACTTTATATTTGAAGTTACAGATCTTGCACCATTAATTGCTGCAATCATTTGTATACTCg                                                          <  1:522561/91‑1 (MQ=255)
                                agagATTTCCAGTCTATAATTGGTGTCAAACACTTTATATTTGAAGTTACAGATCTTGCACCATTAATTGCTGCAATCATTTGTATACTCg                                                          <  1:675563/91‑1 (MQ=255)
                                agagATTTCCAGTCTATAATTGGTGTCAAACACTTTATATTTGAAGTTACAGATCTAGCACCATTAATTGCTGCAATCATTTGTATACTCg                                                          >  1:4142053/1‑91 (MQ=255)
                                agagATTTCCAGTCTATAATTGGTGTCAAACACTTTATATTTGAAGTTACAGATCTAGCACCATTAATTGCTGCAATCATTTGTATACTCg                                                          >  1:4142003/1‑91 (MQ=255)
                                     tttCCAGTCTATAATTGGTGTCAAACACTTTATATTTGAAGTTACAGATCTTGCACCATTAATTGCTGCAATCATTTGTATACTCGTTTTc                                                     >  1:2771890/1‑91 (MQ=255)
                                      ttCCAGTCTATAATTGGTGTCAAACACTTTATATTTGAAGTTACAGATCTAGCACCATTAATTGCTTCAATCATTTGTATACTCGTTTTCa                                                    >  1:1772260/1‑91 (MQ=255)
                                       tCCAGTCTATAATTGGTGTCAAACACTTTATATTTGAAGTTACAGATCTTGCACCATTAATTGCTGCAATCATTTGTATACTCGTTTTCaa                                                   <  1:1289419/91‑1 (MQ=255)
                                       tCCAGTCTATAATTGGTGTCAAACACTTTATATTTGAAGTTACAGATCTAGCACCATTAATTGCTGCAATCATTTGTATACTCGTTTTCaa                                                   >  1:3987362/1‑91 (MQ=255)
                                       tCCAGTCTATAATTGGTGTCAAACACTTTATATTTGAAGTTACAGATCTAGCACCATTAATTGCTGCAATCATTTGTATACTCGTTTTCaa                                                   >  1:3987475/1‑91 (MQ=255)
                                        ccAGTCTATAATTGGTGTCAAACACTTTATATTTGAAGTTACAGATCTAGCACCATTAATTGCTGCAATCATTTGTATACTCGTTTTCaaa                                                  >  1:2803116/1‑91 (MQ=255)
                                        ccAGTCTATAATTGGTGTCAAACACTTTATATTTGAAGTTACAGATCTAGCACCATTAATTGCTGCAATCATTTGTATACTCGTTTTCaaa                                                  >  1:1546771/1‑91 (MQ=255)
                                         cAGTCTATAATTGGTGTCAAACACTTTATATTTGAAGTTACAGATCTTGCACCATTAATTGCTGCAATCATTTGTATACTCGTTTTCAAat                                                 >  1:128133/1‑91 (MQ=255)
                                         cAGTCTATAATTGGTGTCAAACACTTTATATTTGAAGTTACAGATCTTGCACCATTAATTGCTGCAATCATTTGTATACTCGTTTTCAAat                                                 >  1:80638/1‑91 (MQ=255)
                                         cAGTCTATAATTGGTGTCAAACACTTTATATTTGAAGTTACAGATCTATCACCATTACTTGCTGCAATCTTTTGTATACTCGTTTTCAAat                                                 >  1:1255523/1‑91 (MQ=255)
                                           taaaataATTGGTGTCAAACACTTTATATTTGGAGTTACAGATCTAGCACCATTAATTGCTGCAATCATTTGTATGCTCGTTTTCAAatat                                               <  1:3543072/87‑1 (MQ=255)
                                           gTCTATAATTGGTGTCAAACACTTTATATTTGAAGTTACAGATCTAGCACCATTAATTGCTGCAATCATTTGTATACTCGTTTTCAAatat                                               <  1:2583292/91‑1 (MQ=255)
                                           gTCTATAATTGGTGTCAAACACTTTATATTTGAAGTTACAGATCTAGCACCATTAATTGCTGCAATCATTTGTATACTCGTTTTCAAatat                                               >  1:1765681/1‑91 (MQ=255)
                                           gTCTATAATTGGTGTCAAACACTTTATATTTGAAGTTACAGATCTAGCACCATTAATTGCTGCAATCATTTGTATACTCGTTTTCAAatat                                               <  1:1634162/91‑1 (MQ=255)
                                             cTATAATTGGTGTCAAACACTTTATATTTGAAGTTACAGATCTTGCACCATTAATTGCTGCAATCATTTGTATACTCGTTTTCAAATATaa                                             <  1:2010814/91‑1 (MQ=255)
                                               ataATTGGTGTCAAACACTTTATATTTGAAGTTACAGATCTAGCACCATTAATTGCTGCAATCATTTGTATACTCGTTTTCAAATATaaaa                                           <  1:476700/91‑1 (MQ=255)
                                                  aTTGGTGTCAAACACTTTATATTTGAAGTTACAGATCTAGCACCATTAATTGCTGCAATCATTTGTATACTCGTTTTCAAATATAAAAAgg                                        >  1:1850765/1‑91 (MQ=255)
                                                   ttGGTGTCAAACACTTTATATTTGAAGTTACAGATCTTGCACCATTAATTGCTGCAATCATTTGTATACTCGTTTTCAAATATAAAAAGGt                                       <  1:1132722/91‑1 (MQ=255)
                                                   ttGGTGTCAAACACTTTATATTTGAAGTTACAGATCTAGCACCATTAATTGCTGCAATCATTTGTATACTCGTTTTCAAATATAAAAAGGt                                       <  1:3547713/91‑1 (MQ=255)
                                                    tGGTGTCAAACACTTTATATTTGAAGTTACAGATCTAGCACCATTAATTGCTGCAATCATTTGTATACTCGTTTTCAAATATAAAAAGGTc                                      >  1:2764274/1‑91 (MQ=255)
                                                     ggTGTCAAACACTTTATATTTGAAGTTACAGATCTTGCACCATTAATTGCTGCAATCATTTGTATACTCGTTTTCAAATATAAAAAGGTcc                                     >  1:2215338/1‑91 (MQ=255)
                                                     ggTGTCAAACACTTTATATTTGAAGTTACAGATCTTGCACCATTAATTGCTGCAATCATTTGTATACTCGTTTTCAAATATAAAAAGGTcc                                     <  1:3773746/91‑1 (MQ=255)
                                                       tgtCAAACACTTTATATTTGAAGTTACAGATCTAGCACCATTAATTGCTGCAATCATTTGTATACTCGTTTTCAAATATAAAAAGGTCCaa                                   >  1:6457/1‑91 (MQ=255)
                                                       tgtCAAACACTTTATATTTGAAGTTACAGATCTAGCACCATTAATTGCTGCAATCATTTGTATACTCGTTTTCAAATATAAAAAGGTCCaa                                   >  1:6495/1‑91 (MQ=255)
                                                        gtCAAACACTTTATATTTGAAGTTACAGATCTTGCACCATTAATTGCTGCAATCATTTGTATACTCGTTTTCAAATATAAAAAGGTCCAAc                                  >  1:2066004/1‑91 (MQ=255)
                                                        gtCAAACACTTTATATTTGAAGTTACAGATCTTGCACCATTAATTGCTGCAATCATTTGTATACTCGTTTTCAAATATAAAAAGGTCCAAc                                  >  1:2748635/1‑91 (MQ=255)
                                                        gtCAAACACTTTATATTTGAAGTTACAGATCTAGCACCATTAATTGCTGCAATCATTTGTATACTCGTTTTCAAATATAAAAAGGTCCAAc                                  >  1:2142703/1‑91 (MQ=255)
                                                        gtCAAACACTTTATATTTGAAGTTACAGATCTAGCACCATTAATTGCTGCAATCATTTGTATACTCGTTTTCAAATATAAAAAGGTCCAAc                                  >  1:598479/1‑91 (MQ=255)
                                                         tCAAACACTTTATATTTGAAGTTACAGATCTTGCACCATTAATTGCTGCAATCATTTGTATACTCGTTTTCAAATATAAAAAGGTCCAACt                                 <  1:3886540/91‑1 (MQ=255)
                                                          cAAACACTTTATATTTGAAGTTACAGATCTTGCACCATTAATTGCTGCAATCATTTGTATACTCGTTTTCAAATATAAAAAGGTCCAACtt                                <  1:2121989/91‑1 (MQ=255)
                                                          cAAACACTTTATATTTGAAGTTACAGATCTTGCACCATTAATTGCTGCAATCATTTGTATACTCGTTTTCAAATATAAAAAGGTCCAACtt                                >  1:1353270/1‑91 (MQ=255)
                                                           aaaCACTTTATATTTGAAGTTACAGATCTAGCACCATTAATTGCTGCAATCATTTGTATACTCGTTTTCAAATATAAAAAGGTCCAACTTg                               <  1:4057857/91‑1 (MQ=255)
                                                           aaaCACTTTATATTTGAAGTTACAGATCTAGCACCATTAATTGCTGCAATCATTTGTATACTCGTTTTCAAATATAAAAAGGTCCAACTTg                               <  1:1232808/91‑1 (MQ=255)
                                                            aaCACTTTATATTTGAAGTTACAGATCTAGCACCATTAATTGCTGCAATCATTTGTATACTCGTTTTCAAATATAAAAAGGTCCAACTTGc                              >  1:4125942/1‑91 (MQ=255)
                                                            aaCACTTTATATTTGAAGTTACAGATCTAGCACCATTAATTGCTGCAATCATTTGTATACTCGTTTTCAAATATAAAAAGGTCCAACTTGc                              <  1:3351729/91‑1 (MQ=255)
                                                             acacTTTATATTTGAAGTTACAGATCTAGCACCATTAATTGCTGCAATCATTTGTATACTCGTTTTCAAATATAAAAAGGTCCAACTTGCa                             <  1:1595116/91‑1 (MQ=255)
                                                               acTTTATATTTGAAGTTACAGATCTTGCACCATTAATTGCTGCAATCATTTGTATACTCGTTTTCAAATATAAAAAGGTCCAACTTGCAgg                           >  1:1879916/1‑91 (MQ=255)
                                                               acTTTATATTTGAAGTTACAGATCTAGCACCATTAATTGCTGCAATCATTTGTATACTCGTTTTCAAATATAAAAAGGTCCAACTTGCAgg                           >  1:1006799/1‑91 (MQ=255)
                                                                  ttATATTTGAAGTTACAGATCTAGCACCATTAATTGCTGCAATCATTTGTATACTCGTTTTCAAATATAAAAAGGTCCAACTTGCAGGttt                        <  1:1581633/91‑1 (MQ=255)
                                                                   tataTTTGAAGTTACAGATCTTGCACCATTAATTGCTGCAATCATTTGTATACTCGTTTTCAAATATAAAAAGGTCCAACTTGCAGGTTTa                       <  1:2748156/91‑1 (MQ=255)
                                                                    atatTTGAAGTTACAGATCTAGCACCATTAATTGCTGCAATCATTTGTATACTCGTTTTCAAATATAAAAAGGTCCAACTTGCAGGTTTaa                      >  1:4116546/1‑91 (MQ=255)
                                                                    atatTTGAAGTTACAGATCTAGCACCATTAATTGCTGCAATCATTTGTATACTCGTTTTCAAATATAAAAAGGTCCAACTTGAAGGTTTaa                      >  1:2300303/1‑91 (MQ=255)
                                                                        ttGAAGTTACAGATCTTGCACCATTAATTGCTGCAATCATTTGTATACTCGTTTTCAAATATAAAAAGGTCCAACTTGCAGGTTTAAAAtt                  >  1:2979251/1‑91 (MQ=255)
                                                                        ttGAAGTTACAGATCTTGCACCATTAATTGCTGCAATCATTTGTATACTCGTTTTCAAATATAAAAAGGTCCAACTTGCAGGTTTAAAAtt                  >  1:2979185/1‑91 (MQ=255)
                                                                         tGAAGTTACAGATCTTGCACCATTAATTGCTGCAATCATTTGTATACTCGTTTTCAAATATAAAAAGGTCCAACTTGCAGGTTTAAAATtc                 <  1:2814788/91‑1 (MQ=255)
                                                                           aaGTTACAGATCTAGCACCATTAATTGCTGCAATCATTTGTATACTCGTTTTCAAATATAAAAAGGTCCAACTTGCAGGTTTAAAATtctc               <  1:3317628/91‑1 (MQ=255)
                                                                             gTTACAGATCTAGCACCATTAATTGCTGCAATCATTTGTATACTCGTTTTCAAATATAAAAAGGTCCAACTTGCAGGTTTAAAATTCTCaa             <  1:2387853/91‑1 (MQ=255)
                                                                              ttACAGATCTAGCACCATTAATTGCTGCAATCATTTGTATACTCGTTTTCAAATATAAAAAGGTCCAACTTGCAGGTTTAAAATTCTCAAt            <  1:679237/91‑1 (MQ=255)
                                                                              ttACAGATCTAGCACCATTAATTGCTGCAATCATTTGTATACTCGTTTTCAAATATAAAAAGGTCCAACTTGCAGGTTTAAAATTCTCAAt            >  1:2545555/1‑91 (MQ=255)
                                                                               tACAGATCTTGCACCATTAATTGCTGCAATCATTTGTATACTCGTTTTCAAATATAAACCGGTCCAACTTGCAGGTTTAAAATTCTCAATc           >  1:1235588/1‑91 (MQ=255)
                                                                                        aGCACCATTAATTGCTGCAATCATTTGTATACTCGTTTTCAAATATAAAAAGGTCCAACTTGCAGGTTTAAAATTCTCAATCAGCCTGaaa  <  1:3057100/91‑1 (MQ=255)
                                                                                        |                                                                                           
TTGTCATTACAATGGCGTTATCCATTATGCTCAGAGATTTCCAGTCTATAATTGGTGTCAAACACTTTATATTTGAAGTTACAGATCTAGCACCATTAATTGCTGCAATCATTTGTATACTCGTTTTCAAATATAAAAAGGTCCAACTTGCAGGTTTAAAATTCTCAATCAGCCTGAAA  >  USA300TCH1516_ALE/2454538‑2454716

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