Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 2,165,416 A→T 7.8% intergenic (‑384/‑94) tsaE ← / → ilvD tRNA threonylcarbamoyladenosine biosynthesis protein TsaE/Dihydroxy‑acid dehydratase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE2,165,4160AT 7.8% 306.9 / 13.6 90intergenic (‑384/‑94)tsaE/ilvDtRNA threonylcarbamoyladenosine biosynthesis protein TsaE/Dihydroxy‑acid dehydratase
Reads supporting (aligned to +/- strand):  ref base A (39/44);  new base T (5/2);  total (44/46)
Fisher's exact test for biased strand distribution p-value = 2.61e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.83e-01

Only 100 of 103 total aligned reads displayed.

TATACTCTTTAGGACTCGAACGTTAGTAAATATTTACTAAACGCTTTAAGTCCTATTTCTGTTTGAATGGGACTTGTAAACGTCCCAATAATATTGGGACGTTTTTTTATGTTTTATCTTTCAATTACTTATTTTTATTACTATAAAACATGATTAATCATTAAAATTTACGGGGGAATTT  >  USA300TCH1516_ALE/2165326‑2165506
                                                                                          |                                                                                          
tataCTCTTTAGGACTCGAACGTTAGTAAATATTTACTAAACGCTTTAAGTCCTATTTCTGTTTGAATGGGACTTGTAAACGTCCCaataa                                                                                            <  1:121021/91‑1 (MQ=255)
 ataCTCTTTAGGACTCGAACGTTAGTAAATATTTACTAAACGCTTTAAGTCCTATTTCTGTTTGAATGGGACTTGTAAACGTCCCaataat                                                                                           >  1:3499867/1‑91 (MQ=255)
 ataCTCTTTAGGACTCGAACGTTAGTAAATATTTACTAAACGCTTTAAGTCCTATTTCTGTTTGAATGGGACTTGTAAACGTCCCaataat                                                                                           >  1:3500217/1‑91 (MQ=255)
   aCTCTTTAGGACTCGAACGTTAGTAAATATTTACTAAACGCTTTAAGTCCTATTTCTGTTTGAATGGGACTTGTAAACGTCCCAAtattat                                                                                         <  1:1876955/91‑1 (MQ=255)
   aCTCTTTAGGACTCGAACGTTAGTAAATATTTACTAAACGCTTTAAGTCCTATTTCTGTTTGAATGGGACTTGTAAACGTCCCAATAatat                                                                                         <  1:1516435/91‑1 (MQ=255)
   aCTCTTTAGGACTCGAACGTTAGTAAATATTTACTAAACGCTTTAAGTCCTATTTCTGTTTGAATGGGACTTGTAAACGTCCCAATAatat                                                                                         >  1:4103406/1‑91 (MQ=255)
    ctctTTAGGACTCGAACGTTAGTAAATATTTACTAAACGCTTTAAGTCCTATTTCTGTTTGAATGGGACTTGTAAACGTCCCAATattatt                                                                                        <  1:1829822/91‑1 (MQ=255)
     tctTTAGGACTCGAACGTTAGTAAATATTTACTAAACGCTTTAAGTCCTATTTCTGTTTGAATGGGACTTGTAAACGTCCCAATATTATTg                                                                                       <  1:708802/91‑1 (MQ=255)
      ctTTAGGACTCGAACGTTAGTAAATATTTACTAAACGCTTTAAGTCCTATTTCTGTTTGAATGGGACTTGTAAACGTCCCAATAATATTgg                                                                                      <  1:1836543/91‑1 (MQ=255)
       tttAGGACTCGAACGTTAGTAAATATTTACTAAACGCTTTAAGTCCTATTTCTGTTTGAATGGGACTTGTAAACGTCCCAATAATATTggg                                                                                     >  1:3118303/1‑91 (MQ=255)
       tttAGGACTCGAACGTTAGTAAATATTTACTAAACGCTTTAAGTCCTATTTCTGTTTGAATGGGACTTGTAAACGTCCCAATAATATTggg                                                                                     <  1:1386118/91‑1 (MQ=255)
         tAGGACTCGAACGTTAGTAAATATTTACTAAACGCTTTAAGTCCTATTTCTGTTTGAATGGGACTTGTAAACGTCCCAATAATATTGGGAc                                                                                   <  1:1131578/91‑1 (MQ=255)
          aGGACTCGAACGTTAGTAAATATTTACTAAACGCTTTAAGTCCTATTTCTGTTTGAATGGGACTTGTAAACGTCCCAATATTATTGGGACg                                                                                  <  1:179490/91‑1 (MQ=255)
          aGGACTCGAACGTTAGTAAATATTTACTAAACGCTTTAAGTCCTATTTCTGTTTGAATGGGACTTGTAAACGTCCCAATAATATTGGGACg                                                                                  >  1:701594/1‑91 (MQ=255)
          aGGACTCGAACGTTAGTAAATATTTACTAAACGCTTTAAGTCCTATTTCTGTTTGAATGGGACTTGTAAACGTCCCAATAATATTGGGACg                                                                                  >  1:2605903/1‑91 (MQ=255)
          aGGACTCGAACGTTAGTAAATATTTACTAAACGCTTTAAGTCCTATTTCTGTTTGAATGGGACTTGTAAACGTCCCAATAATATTGGGACg                                                                                  >  1:2606042/1‑91 (MQ=255)
           ggACTCGAACGTTAGTAAATATTTACTAAACGCTTTAAGTCCTATTTCTGTTTGAATGGGACTTGTAAACGTCCCAATAATATTGGGACGt                                                                                 >  1:276333/1‑91 (MQ=255)
               tCGAACGTTAGTAAATATTTACTAAACGCTTTAAGTCCTATTTCTGTTTGAATGGGACTTGTAAACGTCCCAATAATATTGGGACGttttt                                                                             >  1:2750368/1‑91 (MQ=255)
                cGAACGTTAGTAAATATTTACTAAACGCTTTAAGTCCTATTTCTGTTTGAATGGGACTTGTAAACGTCCCAATAATATTGGGACGtttttt                                                                            >  1:367343/1‑91 (MQ=255)
                cGAACGTTAGTAAATATTTACTAAACGCTTTAAGTCCTATTTCTGTTTGAATGGGACTTGTAAACGTCCCAATAATATTGGGACGtttttt                                                                            >  1:770864/1‑91 (MQ=255)
                 gAACGTTAGTAAATATTTACTAAACGCTTTAAGTCCTATTTCTGTTTGAATGGGACTTGTAAACGTCCCAATAATATTGGGACGttttttt                                                                           <  1:2114696/91‑1 (MQ=255)
                  aaCGTTAGTAAATATTTACTAAACGCTTTAAGTCCTATTTCTGTTTGAATGGGACTTGTAAACGTCCCAATAATATTGGGACGTTTTTTTa                                                                          <  1:2319451/91‑1 (MQ=255)
                    gtttAGTAAATATTTACTAAACGCTTTAAGTCCTATTTCTGTTTGAATGGGACTTGTAAACGTCCCAATAATATTGGGACGTTTTTTTATg                                                                        >  1:1933586/3‑91 (MQ=255)
                     gTTAGTATATATTTACTAAACGCTTTAAGTCCTATTTCTGTTTGAATGGGACTTGTAAACGTCCCAATAATATTGGGACGTGTTTTTATGt                                                                       <  1:2086545/91‑1 (MQ=255)
                       tAGTAAATATTTACTAAACGCTTTAAGTCCTATTTCTGTTTGAATGGGACTTGTAAACGTCCCAATAATATTGGGACGTTTTTTTATGttt                                                                     >  1:3254370/1‑91 (MQ=255)
                       tAGTAAATATTTACTAAACGCTTTAAGTCCTATTTCTGTTTGAATGGGACTTGTAAACGTCCCAATAATATTGGGACGTTTTTTTATGttt                                                                     >  1:3831541/1‑91 (MQ=255)
                         gTAAATATTTACTAAACGCTTTAAGTCCTATTTCTGTTTGAATGGGACTTGTAAACGTCCCAATAATATTGGGACGTTTTTTTATGTTTTa                                                                   <  1:1690697/91‑1 (MQ=255)
                         gTAAATATTTACTAAACGCTTTAAGTCCTATTTCTGTTTGAATGGGACTTGTAAACGTCCCAATAATATTGGGACGTTTTTTTATGTTTTa                                                                   <  1:2347234/91‑1 (MQ=255)
                          tAAATATTTACTAAACGCTTTAAGTCCTATTTCTGTTTGAATGGGACTTGTAAACGTCCCAATAATATTGGGACGTTTTTTTATGTTTTAt                                                                  >  1:2841840/1‑91 (MQ=255)
                          tAAATATTTACTAAACGCTTTAAGTCCTATTTCTGTTTGAATGGGACTTGTAAACGTCCCAATAATATTGGGACGTTTTTTTATGTTTTAt                                                                  >  1:2841580/1‑91 (MQ=255)
                          tAAATATTTACTAAACGCTTTAAGTCCTATTTCTGTTTGAATGGGACTTGTAAACGTCCCAATAATATTGGGACGTTTTTTTATGTTTTAt                                                                  <  1:3675846/91‑1 (MQ=255)
                           aaaTATTTACTAAACGCTTTAAGTCCTATTTCTGTTTGAATGGGACTTGTAAACGTCCCAATAATATTGGGACGTTTTTTTATGTTTTATc                                                                 >  1:125985/1‑91 (MQ=255)
                           aaaTATTTACTAAACGCTTTAAGTCCTATTTCTGTTTGAATGGGACTTGTAAACGTCCCAATAATATTGGGACGTTTTTTTATGTTTTATc                                                                 <  1:457845/91‑1 (MQ=255)
                           aaaTATTTACTAAACGCTTTAAGTCCTATTTCTGTTTGAATGGGACTTGTAAACGTCCCAATAATATTGGGACGTTTTTTTATGTTTTATc                                                                 >  1:2841919/1‑91 (MQ=255)
                           aaaTATTTACTAAACGCTTTAAGTCCTATTTCTGTTTGAATGGGACTTGTAAACGTCCCAATAATATTGGGACGTTTTTTTATGTTTTATc                                                                 <  1:3618959/91‑1 (MQ=255)
                           aaaTATTCACTAAACGCTTTAAGTCCTATTTCTGTTTGAATGGGACTTGTAAACGTCCCAATAATATTGGGACGTTTTTTTATGTTTTATc                                                                 <  1:3620138/91‑1 (MQ=255)
                                  tACTAAACGCTTTAAGTCCTATTTCTGTTTGAATGGGACTTGTAAACGTCCCAATAATATTGGGACGTTTTTTTATGTTTTATCTTTCAAt                                                          <  1:212307/91‑1 (MQ=255)
                                  tACTAAACGCTTTAAGTCCTATTTCTGTTTGAATGGGACTTGTAAACGTCCCAATAATATTGGGACGTTTTTTTATGTTTTATCTTTCAAt                                                          >  1:3492971/1‑91 (MQ=255)
                                  tACTAAACGCTTTAAGTCCTATTTCTGTTTGAATGGGACTTGTAAACGTCCCAATAATATTGGGACGTTTTTTTATGTTTTATCTTTCAAt                                                          <  1:3961252/91‑1 (MQ=255)
                                  tACTAAACGCTTTAAGTCCTATTTCTGTTTGAATGGGACTTGTAAACGTCCCAATAATATTGGGACGTTTTTTTATGCTTTATCTTTCAAt                                                          >  1:1431288/1‑91 (MQ=255)
                                   aCTAAACGCTTTAAGTCCTATTTCTGTTTGAATGGGACTTGTAAACGTCCCAATAATATTGGGACGTTTTTTTATGTTTTATCTTTCAAtt                                                         >  1:1023990/1‑91 (MQ=255)
                                   aCTAAACGCTTTAAGTCCTATTTCTGTTTGAATGGGACTTGTAAACGTCCCAATAATATTGGGACGTTTTTTTATGTTTTATCTTTCAAtt                                                         <  1:15264/91‑1 (MQ=255)
                                    cTAAACGCTTTAAGTCCTATTTCTGTTTGAATGGGACTTGTAAACGTCCCAATAATATTGGGACGTTTTTTTATGTTTTATCTTTCAATTa                                                        <  1:1825458/91‑1 (MQ=255)
                                      aaaCGCTTTAAGTCCTATTTCTGTTTGAATGGGACTTGTAAACGTCCCAATAATATTGGGACGTTTTTTTATGTTTTATCTTTCAATTACt                                                      <  1:67566/91‑1 (MQ=255)
                                        aCGCTTTAAGTCCTATTTCTGTTTGAATGGGACTTGTAAACGTCCCAATAATATTGGGACGTTTTTTTATGTTTTATCTTTCAATTACTTa                                                    <  1:251199/91‑1 (MQ=255)
                                         cGCTTTAAGTCCTATTTCTGTTTGAATGGGACTTGTAAACGTCCCAATAATATTGGGACGTTTTTTTATGTTTTATCTTTCAATTACTTAt                                                   >  1:1815901/1‑91 (MQ=255)
                                          gCTTTAAGTCCTATTTCTGTTTGAATGGGACTTGTAAACGTCCCAATAATATTGGGACGTTTTTTTATGTTTTATCTTTCAATTACTTAtt                                                  <  1:3390807/91‑1 (MQ=255)
                                           cTTTAAGTCCTATTTCTGTTTGAATGGGACTTGTAAACGTCCCAATAATATTGGGACGTTTTTTTATGTTTTATCTTTCAATTACTTAttt                                                 >  1:4033758/1‑91 (MQ=255)
                                           cTTTAAGTCCTATTTCTGTTTGAATGGGACTTGTAAACGTCCCAATAATATTGGGACGTTTTTTTATGTTTTATCTTTCAATTACTTAttt                                                 >  1:96486/1‑91 (MQ=255)
                                           cTTTAAGTCCTATTTCTGTTTGAATGGGACTTGTAAACGTCCCAATAATATTGGGACGTTTTTTTATGTTTTATCTTTCAATTACTTAttt                                                 <  1:1548909/91‑1 (MQ=255)
                                            tttAAGTCCTATTTCTGTTTGAATGGGACTTGTAAACGTCCCAATAATATTGGGACGTTTTTTTATGTTTTATCTTTCAATTACTTAtttt                                                >  1:762489/1‑91 (MQ=255)
                                              tAAGTCCTATTTCTGTTTGAATGGGACTTGTAAACGTCCCAATAATATTGGGACGTTTTTTTATGTTTTATCTTTCAATTACTTATTTtta                                              <  1:329621/91‑1 (MQ=255)
                                               aaGTCCTATTTCTGTTTGAATGGGACTTGTAAACGTCCCAATAATATTGGGACGTTTTTTTATGTTTTATCTTTCAATTACTTATTTttat                                             >  1:1210719/1‑91 (MQ=255)
                                                aGTCCTATTTCTGTTTGAATGGGACTTGTAAACGTCCCAATAATATTGGGACGTTTTTTTATGTTTTATCTTTCAATTACTTATTTttatt                                            >  1:1732549/1‑91 (MQ=255)
                                                   ccTATTTCTGTTTGAATGGGACTTGTAAACGTCCCAATAATATTGGGACGTTTTTTTATGTTTTATCTTTCAATTACTTATTTTTATTACt                                         <  1:504988/91‑1 (MQ=255)
                                                   ccTATTTCTGTTTGAATGGGACTTGTAAACGTCCCAATAATATTGGGACGTTTTTTTATGTTTTATCTTTCAATTACTTATTTTTATTACt                                         <  1:3941980/91‑1 (MQ=255)
                                                   ccTATTTCTGTTTGAATGGGACTTGTAAACGTCCCAATAATATTGGGACGTTTTTTTATGTTTTATCTTTCAATTACTTATTTTTATTACt                                         <  1:503717/91‑1 (MQ=255)
                                                       tttCTGTTTGAATGGGACTTGTAAACGTCCCAATAATATTGGGACGTTTTTTTATGTTTTATCTTTCAATTACTTATTTTTATTACTATaa                                     <  1:4055483/91‑1 (MQ=255)
                                                       tttCTGTTTGAATGGGACTTGTAAACGTCCCAATAATATTGGGACGTTTTTTTATGTTTTATCTTTCAATTACTTATTTTTATTACTATaa                                     >  1:866537/1‑91 (MQ=255)
                                                       tttCTGTTTGAATGGGACTTGTAAACGTCCCAATAATATTGGGACGTTTTTTTATGTTTTATCTTTCAATTACTTATTTTTATTACTATaa                                     >  1:2509017/1‑91 (MQ=255)
                                                             tttGAATGGGACTTGTAAACGTCCCAATAATATTGGGACGTTTTTTTATGTTTTATCTTTCAATTACTTATTTTTATTACTATAAAACATg                               >  1:2544557/1‑91 (MQ=255)
                                                             tttGAATGGGACTTGTAAACGTCCCAATAATATTGGGACGTTTTTTTATGTTTTATCTTTCAATTACTTATTTTTATTACTATAAAACATg                               >  1:3002657/1‑91 (MQ=255)
                                                               tGAATGGGACTTGGAAACGTCCCAATAATATTGGGACGTTTTTGTATGTTTTATCTTTCAATTACTTATTTTTATTACTATAAAACATGAt                             <  1:4145568/91‑1 (MQ=255)
                                                                  aTGGGACTTGTAAACGTCCCAATAATATTGGGACGTTTTTTTATGTTTTATCTTTCAATTACTTATTTTTATTACTATAAAACATGATTaa                          >  1:3356439/1‑91 (MQ=255)
                                                                     ggACTTGTAAACGTCCCAATAATATTGGGACGTTTTTTTATGTTTTATCTTTCAATTACTTATTTTTATTACTATAAAACATGATTAATCa                       <  1:3915764/91‑1 (MQ=255)
                                                                     ggACTTGTAAACGTCCCAATAATATTGGGACGTTTTTTTATGTTTTATCTTTCAATTACTTATTTTTATTACTATAAAACATGATTAATCa                       <  1:1789905/91‑1 (MQ=255)
                                                                     ggACTTGTAAACGTCCCAATAATATTGGGACGTTTTTTTATGTTTTATCTTTCAATTACTTATTTTTATTACTATAAAACATGATTAATCa                       <  1:1789463/91‑1 (MQ=255)
                                                                      gACTTGTAAACGTCCCAATAATATTGGGACGTTTTTTTATGTTTTATCTTTCAATTACTTATTTTTATTACTATAAAACATGATTAATCAt                      >  1:4036296/1‑91 (MQ=255)
                                                                      gACTTGTAAACGTCCCAATAATATTGGGACGTTTTTTTATGTTTTATCTTTCAATTACTTATTTTTATTACTATAAAACATGATTAATCAt                      <  1:911312/91‑1 (MQ=255)
                                                                       aCTTGTAAACGTCCCAATAATATTGGGACGTTTTTTTATGTTTTATCTTTCAATTACTTATTTTTATTACTATAAAACATGATTAATCAtt                     <  1:173108/91‑1 (MQ=255)
                                                                        cTTGTAAACGTCCCAATAATATTGGGACGTTTTTTTATGTTTTATCTTTCAATTACTTATTTTTATTACTATAAAACATGATTAATCATTa                    <  1:3035646/91‑1 (MQ=255)
                                                                        cTTGTAAACGTCCCAATAATATTGGGACGTTTTTTTATGTTTTATCTTTCAATTACTTATTTTTATTACTATAAAACATGATTAATCATTa                    <  1:3035446/91‑1 (MQ=255)
                                                                        cTTGTAAACGTCCCAATAATATTGGGACGTTTTTTTATGTTTTATCTTTCAATTACTTATTTTTATTACTATAAAACATGATTAATCATTa                    <  1:3034744/91‑1 (MQ=255)
                                                                         ttGTAAACGTCCCAATAATATTGGGACGTTTTTTTATGTTTTATCTTTCAATTACTTATTTTTATTACTATAAAACATGATTAATCATTaa                   >  1:3476104/1‑91 (MQ=255)
                                                                         ttGTAAACGTCCCAATAATATTGGGACGTTTTTTTATGTTTTATCTTTCAATTACTTATTTTTATTACTATAAAACATGATTAATCATTaa                   >  1:3453155/1‑91 (MQ=255)
                                                                         ttGTAAACGTCCCAATAATATTGGGACGTTTTTTTATGTTTGATCTTTCAATTACTTATTTTTATTACTATAAAACATGATTAATCATTaa                   >  1:1079849/1‑91 (MQ=255)
                                                                          aaaaaaCGTCCCAATATTATTGGGACGTTTTTTTATGTTTTATCTTTCAATTACTTATTTTTATTACTATAAAACATGATTAATCATTaaa                  >  1:605598/4‑91 (MQ=255)
                                                                           gTAAACGTCCCAATAATATTGGGACGTTTTTTTATGTTTTATCTTTCAATTACTTATTTTTATTACTATAAAACATGATTAATCATTaaaa                 <  1:1375698/91‑1 (MQ=255)
                                                                           gTAAACGTCCCAATAATATTGGGACGTTTTTTTATGTTTTATCTTTCAATTACTTATTTTTATTACTATAAAACATGATTAATCATTaaaa                 <  1:2782349/91‑1 (MQ=255)
                                                                             aaaCGTCCCAATATTATTGGGACGTTTTTTTATGTTTTATCTTTCAATTACTTATTTTTATTACTATAAAACATGATTAATCATTAAAAtt               >  1:1091998/1‑91 (MQ=255)
                                                                              aaCGTCCCAATAATATTGGGACGTTTTTTTATGTTTTATCTTTCAATTACTTATTTTTATTACTATAAAACATGATTAATCATTAAAAttt              <  1:2411871/91‑1 (MQ=255)
                                                                               aCGTCCCAATAATATTGGGACGTTTTTTTATGTTTTATCTTTCAATTACTTATTTTTATTACTATAAAACATGATTAATCATTAAAATTTa             >  1:3076584/1‑91 (MQ=255)
                                                                               aCGTCCCAATAATATTGGGACGTTTTTTTATGTTTTATCTTTCAATTACTTATTTTTATTACTATAAAACATGATTAATCATTAAAATTTa             <  1:2958537/91‑1 (MQ=255)
                                                                                cGTCCCAATAATATTGGGACGTTTTTTTATGTTTTATCTTTCAATTACTTATTTTTATTACTATAAAACATGATTAATCATTAAAATTTAc            <  1:2196537/91‑1 (MQ=255)
                                                                                  tCCTAATAATATTGGGACGTTTTTTTATGTTTTATCTTTCAATTACTTATTTTTATTACTATAAAACATGATTAATCATTAAAATTTACgg          <  1:3214175/91‑1 (MQ=255)
                                                                                  tCCCAATATTATTGGGACGTTTTTTTATGTTTTATCTTTCAATTACTTATTTTTATTACTATAAAACATGATTAATCATTAAAATTTACgg          >  1:1615352/1‑91 (MQ=255)
                                                                                  tCCCAATATTATTGGGACGTTTTTTTATGTTTTATCTTTCAATTACTTATTTTTATTACTATAAAACATGATTAATCATTAAAATTTACgg          >  1:802979/1‑91 (MQ=255)
                                                                                  tCCCAATAATATTGGGACGTTTTTTTATGTTTTATCTTTCAATTACTTATTTTTATTACTATAAAACATGATTAATCATTAAAATTTACgg          >  1:2847445/1‑91 (MQ=255)
                                                                                  tCCCAATAATATTGGGACGTTTTTTTATGTTTTATCTTTCAATTACTTATTTTTATTACTATAAAACATGATTAATCATTAAAATTTACgg          >  1:1292707/1‑91 (MQ=255)
                                                                                  tCCCAATAATATTGGGACGTTTTTTTATGTTTTATCTTTCAATTACTTATTTTTATTACTATAAAACATGATTAATCATTAAAATTTACgg          <  1:1308903/91‑1 (MQ=255)
                                                                                  tCCCAATAATATTGGGACGTTTTTTTATGTTTTATCTTTCAATTACTTATTTTTATTACTATAAAACATGATTAATCATTAAAATTTACgg          >  1:3677620/1‑91 (MQ=255)
                                                                                    ccAATATTATTGGGACGTTTTTTTATGTTTTATCTTTCAATTACTTATTTTTATTACTATAAAACATGATTAATCATTAAAATTTACgggg        >  1:1396698/1‑91 (MQ=255)
                                                                                    ccAATAATATTGGGACGTTTTTTTATGTTTTATCTTTCAATTACTTATTTTTATTACTATAAAACATGATTAATCATTAAAATTTACgggg        >  1:1845327/1‑91 (MQ=255)
                                                                                    ccAATAATATTGGGACGTTTTTTTATGTTTTATCTTTCAATTACTTATTTTTATTACTATAAAACATGATTAATCATTAAAATTTACgggg        <  1:939852/91‑1 (MQ=255)
                                                                                       ataataTTGGGACGTTTTTTTATGTTTTATCTTTCAATTACTTATTTTTATTACTATAAAACATGATTAATCATTAAAATTTACGGGGGaa     <  1:3685132/91‑1 (MQ=255)
                                                                                       ataataTTGGGACGTTTTTTTATGTTTTATCTTTCAATTACTTATTTTTATTACTATAAAACATGATTAATCATTAAAATTTACGGGGGaa     <  1:3253444/91‑1 (MQ=255)
                                                                                        taataTTGGGACGTTTTTTTATGTTTTATCTTTCAATTACTTATTTTTATTACTATAAAACATGATTAATCATTAAAATTTACGGGGGAAt    >  1:1844431/1‑91 (MQ=255)
                                                                                          atatTGGGACGTTTTTTTATGTTTTATCTTTCAATTACTTATTTTTATTACTATAAAACATGATTAATCATTAAAATTTCCGGGGGAAttt  >  1:3758024/1‑91 (MQ=255)
                                                                                          atatTGGGACGTTTTTTTATGTTTTATCTTTCAATTACTTATTTTTATTACTATAAAACATGATTAATCATTAAAATTTACGGGGGAAttt  <  1:3471555/91‑1 (MQ=255)
                                                                                          atatTGGGACGTTTTTTTATGTTTTATCTTTCAATTACTTATTTTTATTACTATAAAACATGATTAATCATTAAAATTTACGGGGGAAttt  >  1:1002849/1‑91 (MQ=255)
                                                                                          |                                                                                          
TATACTCTTTAGGACTCGAACGTTAGTAAATATTTACTAAACGCTTTAAGTCCTATTTCTGTTTGAATGGGACTTGTAAACGTCCCAATAATATTGGGACGTTTTTTTATGTTTTATCTTTCAATTACTTATTTTTATTACTATAAAACATGATTAATCATTAAAATTTACGGGGGAATTT  >  USA300TCH1516_ALE/2165326‑2165506

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