Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 1,685,149 T→A 7.1% intergenic (‑161/+59) USA300HOU_RS08395 ← / ← rpsU hypothetical protein/30S ribosomal protein S21

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE1,685,1490TA 7.1% 596.1 / 24.6 182intergenic (‑161/+59)USA300HOU_RS08395/rpsUhypothetical protein/30S ribosomal protein S21
Reads supporting (aligned to +/- strand):  ref base T (88/81);  new base A (8/5);  total (96/86)
Fisher's exact test for biased strand distribution p-value = 5.75e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.87e-01

Only 100 of 188 total aligned reads displayed.

ATATCAATTTGATTATTTATAAACTATATAACTTTTTTTAATTAAAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACTTGTGATGTTTGTTTTATATAATTAATATTCGTGTTGAGGGAGTCAACAGAGGTATTAATTATTTGAATTTACGTTTACGTGCAGCTTCTG  >  USA300TCH1516_ALE/1685059‑1685239
                                                                                          |                                                                                          
atatCAATTTGATTATTTATAAACTATATAACTTTTTTTAATTAAAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACt                                                                                            <  1:2196352/91‑1 (MQ=255)
 tatCAATTTGATTATTTATAAACTATATAACTTTTTTTAATTAAAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACtt                                                                                           >  1:2200398/1‑91 (MQ=255)
   tCAATTTGATTATTTATAAACTATATAACTTTTTTTAATTAAAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACTtgt                                                                                         >  2:2377465/1‑91 (MQ=255)
   tCAATTTGACTATTTATAAACTATATAACTTTTTTTAATTAAAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACTtgt                                                                                         <  1:3711886/91‑1 (MQ=255)
   tCAATTTGACTATTTATAAACTATATAACTTTTTTTAATTAAAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACTtgt                                                                                         >  1:1955380/1‑91 (MQ=255)
    cAATTTGATTATTTATAAACTATATAACGTTTTTTAATTAAAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACTtgtg                                                                                        <  1:82816/91‑1 (MQ=255)
     aaTTTGATTATTTATAAACTATATAACTTTTTTTAATTAAAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACTTGTGa                                                                                       >  2:184032/1‑91 (MQ=255)
     aaTTTGATTATTTATAAACTATATAACTTTTTTTAATTAAAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACTTGTGa                                                                                       <  2:2686234/91‑1 (MQ=255)
      aTTTGATTATTTATAAACTATATAACTTTTTTTAATTAAAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACTTGTGAt                                                                                      >  2:2452031/1‑91 (MQ=255)
      aTTTGATTATTTATAAACTATATAACTTTTTTTAATTAAAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACTTGTGAt                                                                                      <  1:1973016/91‑1 (MQ=255)
       tttGATTATTTATAAACTATATAACTTTTTTTAATTAAAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACTTGTGAtg                                                                                     >  2:1078387/1‑91 (MQ=255)
        ttGATTATTTATAAACTATATAACTTTTTTTAATTAAAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACTTGTGAtgt                                                                                    >  1:2646473/1‑91 (MQ=255)
        ttGATTATTTATAAACTATATAACTTTTTTTAATTAAAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACTTGTGAtgt                                                                                    >  1:1103001/1‑91 (MQ=255)
        ttGATTATTTATAAACTATATAACTTTTTTTAATTAAAAATTGACAACAAAAAAACATCACAAGTTAGTGTCTGACACTAATATGTGAtgt                                                                                    <  2:1621450/91‑1 (MQ=21)
        ataattattTATAAACTATATAACTTTTTTTAATTAAAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACTTGTGAtgt                                                                                    <  1:257267/88‑1 (MQ=255)
         tGATTATTTATAAACTATATAACTTTTTTTAATTAAAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACTTGTGAtgtt                                                                                   <  1:2484165/91‑1 (MQ=255)
           attattTATAAACTATATAACTTTTTTTAATTAAAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACTTGTGAtgtttg                                                                                 >  2:2201339/1‑91 (MQ=255)
             tatttatAAACTATATAACTTTTTTTAATTAAAAATTGACAACAAACAAACATCACAAGTTAGTGTCTGACACTAATATGTGAtgtttttt                                                                               <  2:3031866/91‑4 (MQ=17)
             tatttatAAACTATATAACTTTTTTTAATTAAAAATTGACAACAAACAAACATCACAAGTTAGTGTCTGACACTAATATGTGAtgtttttt                                                                               <  2:2371864/91‑4 (MQ=17)
             tatttatAAACTATATAACTTTTTTTAATTAAAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACTTGTGAtgtttgtt                                                                               <  2:100752/91‑1 (MQ=255)
             tatttatAAACTATATAACTTTTTTTAATTAAAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACTTGTGAtgtttgtt                                                                               >  1:424290/1‑91 (MQ=255)
                ttatAAACTATATAACTTTTTTTAATTAAAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACTTGTGATGTTTGTTTta                                                                            <  1:3520068/91‑1 (MQ=255)
                  ataAACTATATAACTTTTTTTAATTAAAAATTGACAACAAAAAAACATCACAAGTTAGTGTCTGACACTAATATGTGATgtttttttgttg                                                                          <  1:2887705/91‑5 (MQ=17)
                   taAACTATATAACTTTTTTTAATTAAAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACTTGTGATGTTTGTTTtatat                                                                         <  2:2271397/91‑1 (MQ=255)
                    aaaCTATATAACTTTTTTTAATTAAAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACTTGTGATGTTTGTTTtatata                                                                        >  1:1028962/1‑91 (MQ=255)
                    aaaCTATATAACTTTTTTTAATTAAAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACTTGTGATGTTTGTTTtatata                                                                        >  2:3630612/1‑91 (MQ=255)
                      aCTATATAACTTTTTTTAATTAAAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACTTGTGATGTTTGTTTTATAtaat                                                                      >  2:3689631/1‑91 (MQ=255)
                      aCTATATAACTTTTTTTAATTAAAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACTTGTGATGTTTGTTTTATAtaat                                                                      >  2:3689793/1‑91 (MQ=255)
                       cTATATAACTTTTTTTAATTAAAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACTTGTGATGTTTGTTTTATAtaatt                                                                     <  2:2865002/91‑1 (MQ=255)
                       cTATATAACTTTTTTTAATTAAAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACTTGTGATGTTTGTTTTATAtaatt                                                                     >  1:959862/1‑91 (MQ=255)
                         atataACTTTTTTTAATTAAAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACTTGTGATGTTTGTTTTATAtaattaa                                                                   <  2:1058499/91‑1 (MQ=255)
                          tataACTTTTTTTAATTAAAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACTTGTGATGTTTGTTTTATAtaattaat                                                                  <  2:1758067/91‑1 (MQ=255)
                          tataACTTTTTTTAATTAAAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACTTGTGATGTTTGTTTTATAtaattaat                                                                  <  1:2367910/91‑1 (MQ=255)
                           ataACTTTTTTTAATTAAAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACTTGTGATGTTTGTTTTATATAATTAata                                                                 >  1:3189337/1‑91 (MQ=255)
                             aaCTTTTTTTAATTAAAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACTTGTGATGTTTGTTTTATATAATTAATAtt                                                               <  1:3759663/91‑1 (MQ=255)
                                tttttttAATTAAAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACTTGTGATGTTTGTTTTATATAATTAATATTCgt                                                            >  1:2492539/1‑91 (MQ=255)
                                tttttttAATTAAAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACTTGTGATGTTTGTTTTATATAATTAATATTCgt                                                            <  2:690710/91‑1 (MQ=255)
                                 ttttttAATTAAAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACTTGTGATGTTTGTTTTATATAATTAATATTCgtg                                                           >  1:1555053/1‑91 (MQ=255)
                                  tttttAATTAAAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACTTGTGATGTTTGTTTTATATAATTAATATTCgtgt                                                          <  2:731145/91‑1 (MQ=255)
                                     ttaattaaAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACTTGTGATGTTTGTTTTATATAATTAATATTCGTGTTGa                                                       <  2:672020/91‑1 (MQ=255)
                                     ttaattaaAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACTTGTGATGTTTGTTTTATATAATTAATATTCGTGTTGa                                                       >  1:1423821/1‑91 (MQ=255)
                                      taattaaAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACTTGTGATGTTTGTTTTATATAATTAATATTCGTGTTGAg                                                      >  2:2745943/1‑91 (MQ=255)
                                        attaaAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACTTGTGATGTTTGTTTTATATAATTAATATTCGTGTTGactc                                                    <  1:1518933/91‑4 (MQ=255)
                                        attaaAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACTTGTGATGTTTGTTTTATATAATTAATATTCGTGTTGAggg                                                    <  2:153674/91‑1 (MQ=255)
                                        attaaAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACTTGTGATGTTTGTTTTATATAATTAATATTCGTGTTGAggg                                                    >  1:549567/1‑91 (MQ=255)
                                         ttaaAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACTTGTGATGTTTGTTTTATATAATTAATATTCGTGTTGAGGGa                                                   >  1:3318758/1‑91 (MQ=255)
                                           aaaaaTTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACTTGTGATGTTTGTTTTATATAATTAATATTCGTGTTGAGGGAGt                                                 <  2:179959/91‑1 (MQ=255)
                                                ttGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACTTGTGATGTTTGTTTTATATAATTAATATTCGTGTTGAGGGAGTCAACa                                            <  2:200169/91‑1 (MQ=255)
                                                 tGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACTTGTGATGTTTGTTTTATATAATTAATATTCGTGTTGAGGGAGTCAACag                                           <  1:3646827/91‑1 (MQ=255)
                                                  gACAACAAAAAAACATCACATATTAGTGTCAGACACTAACTTGTGATGTTTGTTTTATATAATTAATATTCGTGTTGAGGGAGTCAACaga                                          >  1:3402684/1‑91 (MQ=255)
                                                  gACAACAAAAAAACATCACATATTAGTGTCAGACACTAACTTGTGATGTTTGTTTTATATAATTAATATTCGTGTTGAGGGAGTCAACaga                                          >  1:3401826/1‑91 (MQ=255)
                                                   acaacaAAAAAACATCACATATTAGTGTCAGACACTAACTTGTGATGTTTGTTTTATATAATTAATATTCGTGTTGAGGGAGTCAACagag                                         <  1:2412331/91‑1 (MQ=255)
                                                    caacaaAAAAACATCACATATTAGTGTCAGACACTAACTTGTGATGTTTGTTTTATATAATTAATATTCGTGTTGAGGGAGTCAACAGAgg                                        >  1:351905/1‑91 (MQ=255)
                                                     aacaaAAAAACATCACATATTAGTGTCAGACACTAACTTGTGATGTTTGTTTTATATAATTAATATTCGTGTTGAGGGAGTCGACAGAGGt                                       >  2:943697/1‑91 (MQ=255)
                                                      acaaAAAAACATCACATATTAGTGTCAGACACTAACTTGTGATGTTTGTTTTATATAATTAATATTCGTGTTGAGGGAGTCAACAGAGGTa                                      >  1:2512042/1‑91 (MQ=255)
                                                      acaaAAAAACATCACATATTAGTGTCAGACACTAACTTGTGATGTTTGTTTTATATAATTAATATTCGTGTTGAGGGAGTCAACAGAGGTa                                      >  2:2852698/1‑91 (MQ=255)
                                                       caaAAAAACATCACATATTAGTGTCAGACACTAACTTGTGATGTTTGTTTTATATAATTAATATTCGTGTTGAGGGAGTCAACAGAGGTat                                     <  2:211369/91‑1 (MQ=255)
                                                         aacaaaCATCACAAGTTAGTGTCTGACACTAATATGTGATGTTTGTTTTATATAATTAATATTCGTGTTGAGGGAGTCAACAGAGGTatta                                   >  2:3477440/4‑91 (MQ=18)
                                                         aaaaaaCATCACATATTAGTGTCAGACACTAACTTGTGATGTTTGTTTTATATAATTAATATTCGTGTTGAGGGAGTCAACAGAGGTatta                                   <  1:2243172/91‑1 (MQ=255)
                                                            aaaCATCACATATTAGTGTCAGACACTAACTTGTGATGTTTGTTTTATATAATTAATATTCGTGTTGAGGGAGTCAACAGAGGTattaatt                                >  2:112340/1‑91 (MQ=255)
                                                            aaaCATCACAACTTAGTGCCACACACTAACTTGTGACGATTGTTCTATATAAGTAATATTCGTGTTGAGGGAGGCAACAGAGGTattaatt                                >  2:3275759/1‑91 (MQ=255)
                                                               cATCACATATTAGTGTCAGACACTAACTTGTGATGTTTGTTTTATATAATTAATATTCGTGTTGAGGGAGTCAACAGAGGTATTAattatt                             <  1:243586/91‑1 (MQ=255)
                                                               cATCACATATTAGTGTCAGACACTAACTTGTGATGTTTGTTTTATATAATTAATATTCGTGTTGAGGGAGTCAACAGAGGTATTAattatt                             <  1:169656/91‑1 (MQ=255)
                                                               cATCACATATTAGTGTCAGACACTAACTTGTGATGTTTGTTTTATATAATTAATATTCGTGTTGAGGGAGTCAACAGAGGTATTAattatt                             <  1:3401009/91‑1 (MQ=255)
                                                               cATCACAAGTTAGTGTCTGACACTAATATGTGATGTTTGTTTTATATAATTAATATTCGTGTTGAGGGAGTCAACAGAGGTATTAattatt                             >  2:3007986/1‑91 (MQ=21)
                                                                 tCACATATTAGTGTCAGACACTAACTTGTGATGTTTGTTTTATATAATTAATATTCGTGTTGAGGGAGTCAACAGAGGTATTAATTATTTg                           >  1:1595153/1‑91 (MQ=255)
                                                                 tCACATATTAGTGTCAGACACTAACTTGTGATGTTTGTTTTATATAATTAATATTCGTGTTGAGGGAGTCAACAGAGGTATTAATTATTTg                           >  2:2854926/1‑91 (MQ=255)
                                                                 tCACATATTAGTGTCAGACACTAACTTGTGATGTTTGTTTTATATAATTAATATTCGTGTTGAGGGAGTCAACAGAGGTATTAATTATTTg                           >  1:2149400/1‑91 (MQ=255)
                                                                 tCACATATTAGTGTCAGACACTAACTTGTGATGTTTGTTTTATATAATTAATATTCGTGTTGAGGGAGTCAACAGAGGTATTAATTATTTg                           >  2:1042001/1‑91 (MQ=255)
                                                                  cacaTATTAGTGTCAGACACTAACTTGTGATGTTTGTTTTATATAATTAATATTCGTGTTGAGGGAGTCAACAGAGGTATTAATTATTTGa                          >  2:3798167/1‑91 (MQ=255)
                                                                   acaTATTAGTGTCAGACACTAACTTGTGATGTTTGTTTTATATAATTAATATTCGTGTTGAGGGAGTCAACAGAGGTATTAATTATTTGaa                         <  2:1888294/91‑1 (MQ=255)
                                                                   acaTATTAGTGTCAGACACTAACTTGTGATGTTTGTTTTATATAATTAATATTCGTGTTGAGGGAGTCAACAGAGGTATTAATTATTTGaa                         <  2:2627262/91‑1 (MQ=255)
                                                                    caagttAGTGTCTGACACTAATATGTGATGTTTGTTTTATATAATTAATATTCGTGTTGAGGGAGTCAACAGAGGTATTAATTATTTGAAt                        >  2:1195865/5‑91 (MQ=37)
                                                                    caagttAGTGTCTGACACTAATATGTGATGTTTGTTTTATATAATTAATATTCGTGTTGAGGGAGTCAACAGAGGTATTAATTATTTGAAt                        >  1:2538557/5‑91 (MQ=37)
                                                                    caTATTAGTGTCAGACACTAACTTGTGATGTTTGTTTTATATAATTAATATTCGTGTTGAGGGAGTCAACAGAGGTATTAATTATTTGAAt                        <  2:2797635/91‑1 (MQ=255)
                                                                     atatTAGTGTCAGACACTAACTTGTGATGTTTGTTTTATATAATTAATATTCGTGTTGAGGGAGTCAACAGAGGTATTAATTATTTGAAtt                       <  2:512812/91‑1 (MQ=255)
                                                                     atatTAGTGTCAGACACTAACTTGTGATGTTTGTTTTATATAATTAATATTCGTGTTGAGGGAGTCAACAGAGGTATTAATTATTTGAAtt                       <  1:772784/91‑1 (MQ=255)
                                                                       gttAGTGTCTGACACTAATATGTGATGTTTGTTTTATATAATTAATATTCGTGTTGAGGGAGTTAACAGAGGTATTATTTATTTGAAttta                     >  2:2475829/2‑91 (MQ=37)
                                                                        ttAGTGTCAGACACTAACTTGTGATGTTTGTTTTATATAATTAATATTCGTGTTGAGGGAGTCAACAGAGGTATTAATTATTTGAAtttac                    <  1:1022267/91‑1 (MQ=255)
                                                                         tAGTGTCAGACACTAACTTGTGATGTTTGTTTTATATAATTAATATTCGTGTTGAGGGAGTCAACAGAGGTATTAATTATTTGAAtttacg                   <  1:3832775/91‑1 (MQ=255)
                                                                         tAGTGTCAGACACTAACTTGTGATGTTTGTTTTATATAATTAATATTCGTGTTGAGGGAGTCAACAGAGGTATTAATTATTTGAAtttacg                   >  1:3542647/1‑91 (MQ=255)
                                                                          aGTGTCAGACACTAACTTGTGATGTTTGTTTTATATAATTAATATTCGTGTTGAGGGAGTCAACAGAGGTATTAATTATTTGAAtttacgt                  <  1:3417153/91‑1 (MQ=255)
                                                                            tgtCAGACACTAACTTGTGATGTTTGTTTTATATAATTAATATTCGTGTTGAGGGAGTCAACAGAGGTATTAATTATTTGAAtttatgttt                <  1:3782131/91‑1 (MQ=255)
                                                                             gtCAGACACTAACTTGTGATGTTTGTTTTATATAATTAATATTCGTGTTGAGGGAGTCAACAGAGGTATTAATTATTTGAAtttacgttta               >  1:2705389/1‑91 (MQ=255)
                                                                              tctgACACTAACTTGTGATGTTTGTTTTATATAATTAATATTCGTGTTGAGGGAGTCAACAGAGGTATTAATTATTTGAAtttacgtttac              >  1:1204641/4‑91 (MQ=255)
                                                                              tCAGACACTAACTTGTGATGTTTGTTTTATATAATTAATATTCGTGTTGAGGGAGTCAACCGAGGTATTAATTATTTGAAtttacgtttac              >  1:593275/1‑91 (MQ=255)
                                                                              tCAGACACTAACTTGTGATGTTTGTTTTATATAATTAATATTCGTGTTGAGGGAGTCAACAGAGGTATTAATTATTTGAAtttacgtttac              <  1:428333/91‑1 (MQ=255)
                                                                               cAGACACTAACTTGTGATGTTTGTTTTATATAATTAATATTCGTGTTGAGGGAGTCAACAGAGGTATTAATTATTTGAAtttacgtttacg             <  2:2938585/91‑1 (MQ=255)
                                                                                aGACACTAACTTGTGATGTTTGTTTTATATAATTAATATTCGTGTTGAGGGAGTCAACAGAGGTATTAATTATTTGAATttacgtttacgt            >  1:309008/1‑91 (MQ=255)
                                                                                aGACACTAACTTGTGATGTTTGTTTTATATAATTAATATTCGTGTTGAGGGAGTCAACAGAGGTATTAATTATTTGAATttacgtttacgt            >  1:1833643/1‑91 (MQ=255)
                                                                                aGACACTAACTTGTGATGTTTGTTTTATATAATTAATATTCGTGTTGAGGGAGTCAACAGAGGTATTAATTATTTGAATttacgtttacgt            >  2:1400111/1‑91 (MQ=255)
                                                                                aGACACTAACTTGTGATGTTTGTTTTATATAATTAATATTCGTGTTGAGGGAGTCAACAGAGGTATTAATTATTTGAATttacgtttacgt            >  1:820568/1‑91 (MQ=255)
                                                                                aGACACTAACTCGTGATGTTTGTGTTATATAATTAATATTCGTGTTGAGGGAGTCAACAGAGGTATTAATTATTTGAATttacgtttacgt            <  2:2963836/91‑1 (MQ=255)
                                                                                  acacTAGCTTGTGATGTTTGTTTTATATAATTAATATTCGTGTTGAGGGAGTCAACAGAGGTATTAATTATTTGAATTTACGTTTACGTGc          <  1:3816698/91‑1 (MQ=255)
                                                                                  acacTAACTTGTGATGTTTGTTTTATATAATTAATATTCGTGTTGAGGGAGTCAACAGAGGTATTAATTATTTGAATTTACGTTTACGTGc          >  1:1846293/1‑91 (MQ=255)
                                                                                   cacTAACTTGTGATGTTTGTTTTATATAATTAATATTCGTGTTGAGGGAGTCAACAGAGGTATTAATTATTTGAATTTACGTTTACGTGCa         >  2:450022/1‑91 (MQ=255)
                                                                                   cacTAACTTGTGATGTTTGTTTTATATAATTAATATTCGTGTTGAGGGAGTCAACAGAGGTATTAATTATTTGAATTTACGTTTACGTGCa         <  1:2326578/91‑1 (MQ=255)
                                                                                    acTAACTTGTGATGTTTGTTTTATATAATTAATATTCGTGTTGAGGGAGTCAACAGAGGTATTAATTATTTGAATTTACGTTTACGTGCAg        >  2:144918/1‑91 (MQ=255)
                                                                                     cTAACTTGTGATGTTTGTTTTATATAATTAATATTCGTGTTGAGGGAGTCAACAGAGGTATTAATTATTTGAATTTACGTTTACGTGCAGc       >  1:665580/1‑91 (MQ=255)
                                                                                       aaCTTGTGATGTTTGTTTTATATAATTAATATTCGTGTTGAGGGAGTCAACAGAGGTATTAATTATTTGAATTTACGTTTACGTGCAGCtt     >  1:1958137/1‑91 (MQ=255)
                                                                                          |                                                                                          
ATATCAATTTGATTATTTATAAACTATATAACTTTTTTTAATTAAAAATTGACAACAAAAAAACATCACATATTAGTGTCAGACACTAACTTGTGATGTTTGTTTTATATAATTAATATTCGTGTTGAGGGAGTCAACAGAGGTATTAATTATTTGAATTTACGTTTACGTGCAGCTTCTG  >  USA300TCH1516_ALE/1685059‑1685239

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