Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 1,366,121 A→T 6.2% intergenic (+420/‑34) rpmG2_1 → / → rpsN2 50S ribosomal protein L33 2/Alternate 30S ribosomal protein S14

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE1,366,1210AT 6.2% 1021.7 / 40.5 304intergenic (+420/‑34)rpmG2_1/rpsN250S ribosomal protein L33 2/Alternate 30S ribosomal protein S14
Reads supporting (aligned to +/- strand):  ref base A (123/161);  new base T (3/16);  total (126/178)
Fisher's exact test for biased strand distribution p-value = 2.79e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.55e-01

Only 100 of 328 total aligned reads displayed.

TAAATATGATGGTTTAGAAAATGATTTAACAAAGAAATGAAACTTTACTGTTGAATTATGTGAGGATTGTGTTATTATATAAATCGTAATAATTACGATTTGATAAAAAGTGAGGTAACTATATATGGCTAAGAAATCTAAAATAGCAAAAGAGAGAAAAAGAGAAGAGTTAGTAAATAAA  >  USA300TCH1516_ALE/1366031‑1366211
                                                                                          |                                                                                          
 aaaTATGATGGTTTAGAAAATGATTTAACAAAGAAATGAAACTTTACTGTTGAATTATGTGAGGATTGTGTTATTATATAAATCGtaataa                                                                                           <  1:2215716/91‑1 (MQ=255)
 aaaTATGATGGTTTAGAAAATGATTTAACAAAGAAATGAAACTTTACTGTTGAATTATGTGAGGATTGTGTTATTATATAAATCGtaataa                                                                                           <  2:2085027/91‑1 (MQ=255)
   atatGATGGTTTAGAAAATGATTTAACAAAGAAATGAAACTTTACTGTTGAATTATGTGAGGATTGTGTTATTATATAAATCGTAattatt                                                                                         <  2:3853756/91‑1 (MQ=255)
   atatGATGGTTTAGAAAATGATTTAACAAAGAAATGAAACTTTACTGTTGAATTATGTGAGGATTGTGTTATTATATAAATCGTAATAAtt                                                                                         <  2:2208419/91‑1 (MQ=255)
   atatGATGGTTTAGAAAATGATTTAACAAAGAAATGAAACTTTACTGTTGAATTATGTGAGGATTGTGTTATTATATAAATCGTAATAAtt                                                                                         <  1:1638586/91‑1 (MQ=255)
    tatGATGGCTTAGAAAATGATTTAACAAAGAAATGAACCTTCACTGTTGAATTATGTGAGGATTGTGTTATCATATAAATCGTAATAATTa                                                                                        >  2:595980/1‑91 (MQ=255)
     atgatgGTTTAGAAAATGATTTAACAAAGAAATGAAACTTTACTGTTGAATTATGTGAGGATTGTGTTATTATATAAATCGTAATAATTAc                                                                                       <  1:2821730/91‑1 (MQ=255)
       gatgGTTTAGAAAATGATTTAACAAAGAAATGAAACTTTACTGTTGAATTATGTGAGGATTGTGTTATTATATAAATCGTAATAATTACGa                                                                                     >  1:1003972/1‑91 (MQ=255)
       gatgGTTTAGAAAATGATTTAACAAAGAAATGAAACTTTACTGTTGAATTATGTGAGGATTGTGTTATTATATAAATCGTAATAATTACGa                                                                                     >  1:61687/1‑91 (MQ=255)
        atgGTTTAGAAAATGATTTAACAAAGAAATGAAACTTTACTGTTGAATTATGTGAGGATTGTGTTATTATATAAATCGTAATAATTACGAt                                                                                    >  1:1370201/1‑91 (MQ=255)
         tgGTTTAGAAAATGATTTAACAAAGAAATGAAACTTTACTGTTGAATTATGTGAGGATTGTGTTATTATATAAATCGTAATAATTACGAtt                                                                                   >  1:1388003/1‑91 (MQ=255)
         tgGTTTAGAAAATGATTTAACAAAGAAATGAAACTTTACTGTTGAATTATGTGAGGATTGTGTTATTATATAAATCGTAATAATTACGAtt                                                                                   <  1:1888757/91‑1 (MQ=255)
         tgGTTTAGAAAATGATTTAACAAAGAAATGAAACTTTACTGTTGAATTATGTGAGGATTGTGTTATTATATAAATCGTAATAATTACGAtt                                                                                   <  1:2988216/91‑1 (MQ=255)
         tgGTTTAGAAAAGGATTTAACAAAGAGATGAAACTTTACTGTTGAATGATGTGAGGATTGTGTTATTAGATAAATCGTAATAATTACGAtt                                                                                   <  2:290397/91‑1 (MQ=255)
          ggTTTAGAAAATGATTTAACGAAGAAATGAAACTTTACTGTTGAATTATGTGAGGATTGTGTTATTATATAAATCGTAATTATTACGAttt                                                                                  <  2:3168242/91‑1 (MQ=255)
          ggTTTAGAAAATGATTTAACAAAGAAATGAAACTTTACTGTTGAATTATGTGAGGATTGTGTTATTATATAAATCGTAATAATTACGAttt                                                                                  <  1:2953122/91‑1 (MQ=255)
           gTTTAGAAAATGATTTGACAAAGAAATGAAACTTTACTGTTGAATTATGTGAGGATTGTGTTATTATATAAATCGTAATAATTACGATTTg                                                                                 <  2:3782138/91‑1 (MQ=255)
           gTTTAGAAAATGATTTAACAAAGAAATGAAACTTTACTGTTGAATTATGTGAGGGTTGTGTTATTATATAAATCGTAATAATTACGAttta                                                                                 >  2:3920942/1‑90 (MQ=255)
           gTTTAGAAAATGATTTAACAAAGAAATGAAACTTTACTGTTGAATTATGTGAGGATTGTGTTATTATATAAATCGTAATAATTACGATTTg                                                                                 >  2:482091/1‑91 (MQ=255)
            tttAGAAAATGATTTAACAAAGAAATGAAACTTTACTGTTGAATTATGTGAGGATTGTGTTATTATATAAATCGTAATTATTACGATTTGa                                                                                <  1:3035819/91‑1 (MQ=255)
             ttAGAAAATGATTTAACAAAGAAATGAAACTTTACTGTTGAATTATGTGAGGATTGTGTTATTATATAAATCGTAATAATTACGATTTGAt                                                                               <  1:880739/91‑1 (MQ=255)
             ttAGAAAATGATTTAACAAAGAAATGAAACTTTACTGTTGAATTATGTGAGGATTGTGTTATTATATAAATCGTAATAATTACGATTTGAt                                                                               <  1:2794470/91‑1 (MQ=255)
              tAGAAAATGATTTAACAAAGAAATGAAACTTTACTGTTGAATTATGTGAGGATTGTGTTATTATATAAATCGTAATAATTACGATTTGATa                                                                              <  2:1476330/91‑1 (MQ=255)
                 aaaaTGATTTAACAAAGAAATGAAACTTTACTGTTGAATTATGTGAGGATTGTGTTATTATATAAATCGTAATAATTACGATTTGATaaaa                                                                           <  2:3743403/91‑1 (MQ=255)
                  aaaTGATTTAACAAAGAAATGAAACTTTACTGTTGAATTATGTGAGGATTGTGTTATTATATAAATCGTAATAATTACGATTTGATaaaaa                                                                          <  1:1695507/91‑1 (MQ=255)
                  aaaTGATTTAACAAAGAAATGAAACTTTACTGTTGAATTATGTGAGGATTGTGTTATTATATAAATCGTAATAATTACGATTTGATaaaaa                                                                          >  1:868039/1‑91 (MQ=255)
                   aaTGATTTAACAAAGAAATGAAACTTTACTGTTGAATTATGTGAGGATTGTGTTATTATATAAATCGTAATAATTACGATTTGATAAAAAg                                                                         <  2:198547/91‑1 (MQ=255)
                     tGATTTAACAAAGAAATGAAACTTTACTGTTGAATTATGTGAGGATTGTGTTATTATATAAATCGTAATAATTACGATTTGATAAAAAGTg                                                                       >  1:1924263/1‑91 (MQ=255)
                     tGATTTAACAAAGAAATGAAACTTTACTGTTGAATTATGTGAGGATTGTGTTATTATATAAATCGTAATAATTACGATTTGATAAAAAGTg                                                                       <  2:3491282/91‑1 (MQ=255)
                     tGATTTAACAAAGAAATGAAACTTTACTGTTGAATTATGTGAGGATTGTGTTATTATATAAATCGTAATAATTACGATTTGATAAAAAGTg                                                                       >  2:665660/1‑91 (MQ=255)
                     tGATTTAACAAAGAAATGAAACTTTACTGTTGAATTATGTGAGGATTGTGTTATTATATAAATCGTAATAATTACGATTTGATAAAAAGTg                                                                       >  2:1569377/1‑91 (MQ=255)
                         ttAACAAAGAAATGAAACTTTACTGTTGAATTATGTGAGGATTGTGTTATTATATAAATCGTAATAATTACGATTTGATAAAAAGTGAGGt                                                                   >  1:3888415/1‑91 (MQ=255)
                         ttAACAAAGAAATGAAACTTTACTGTTGAATTATGTGAGGATTGTGTTATTATATAAATCGTAATAATTACGATTTGATAAAAAGTGAGGt                                                                   >  1:2405477/1‑91 (MQ=255)
                          tAACAAAGAAATGAAACTTTACTGTTGAATTATGTGAGGATTGTGTTATTATATAAATCGTAATAATTACGATTTGATAAAAAGTGAGGTa                                                                  >  2:1971991/1‑91 (MQ=255)
                          tAACAAAGAAATGAAACTTTACTGTTGAATTATGTGAGGATTGTGTTATTATATAAATCGTAATAATTACGATTTGATAAAAAGTGAGGTa                                                                  >  1:1250969/1‑91 (MQ=255)
                            aCAAAGAAATGAAACTTTACTGTTGAATTATGTGAGGATTGTGTTATTATATAAATCGTAATAATTACGATTTGATAAAAAGTGAGGTAAc                                                                <  1:2040733/91‑1 (MQ=255)
                            aCAAAGAAATGAAACTTTACTGTTGAATTATGTGAGGATTGTGTTATTATATAAATCGTAATAATTACGATTTGATAAAAAGTGAGGTAAc                                                                <  1:795723/91‑1 (MQ=255)
                             cAAAGAAATGAAACTTTACTGTTGAATTATGTGAGGATTGTGTTATTATATAAATCGTAATAATTACGATTTGATAAAAAGTGAGGTAACt                                                               >  2:3496275/1‑91 (MQ=255)
                              aaaGAAATGAAACTTTACTGTTGAATTATGTGAGGATTGTGTTATTATATAAATCGTAATAATTACGATTTGATAAAAAGTGAGGTAACta                                                              <  2:1389567/91‑1 (MQ=255)
                                aGAAATGAAACTTTACTGTTGAATTATGTGAGGATTGTGTTATTATATAAATCGTAATAATTACGATTTGATAAAAAGTGAGGTAACtata                                                            <  1:957455/91‑1 (MQ=255)
                                aGAAATGAAACTTTACTGTTGAATTATGTGAGGATTGTGTTATTATATAAATCGTAATAATTACGATTTGATAAAAAGTGAAGTAACtata                                                            >  1:1769128/1‑91 (MQ=255)
                                 gAAATGAAACTTTACTGTTGAATTATGTGAGGATTGTGTTATTATATAAATCGTAATAATTACGATTTGATAAAAAGTGAGGTAACtatat                                                           >  1:3180347/1‑91 (MQ=255)
                                     tGAAACTTTACTGTTGAATTATGTGAGGATTGTGTTATTATATAAATCGTAATAATTCCGATTTGATAAAAAGTGAGGTAACTATATATgg                                                       >  1:3504393/1‑91 (MQ=255)
                                     tGAAACTTTACTGTTGAATTATGTGAGGATTGTGTTATTATATAAATCGTAATAATTACGATTTGATAAAAAGTGAGGTAACTATATATgg                                                       <  2:3816698/91‑1 (MQ=255)
                                     tGAAACTTTACTGTTGAATTATGTGAGGATTGTGTTATTATATAAATCGTAATAATTACGATTTGATAAAAAGTGAGGTAACTATATATgg                                                       <  1:1991267/91‑1 (MQ=255)
                                     tGAAACTTTACTGTTGAATTATGTGAGGATTGTGTTATTATATAAATCGTAATAATTACGATTTGATAAAAAGTGAGGTAACTATATATgg                                                       <  1:2679843/91‑1 (MQ=255)
                                     tGAAACTTTACTGTTGAATTATGTGAGGATTGTGTTATTATATAAATCGTAATAATTACGATTTGATAAAAAGTGAGGTAACTATATATgg                                                       <  2:2385941/91‑1 (MQ=255)
                                        aaCTTTACTGTTGAATTATGTGAGGATTGTGTTATTATATAAATCGTAATAATTACGATTTGATAAAAAGTGAGGTAACTATATATGGCTa                                                    >  2:2819279/1‑91 (MQ=255)
                                        aaCTTTACTGTTGAATTATGAGAGGATCGTGTTATTATATAAATCGTAATAATTACGATTTGATAAAAAGTGAGGTAACTATATATGGCTa                                                    >  2:1910851/1‑91 (MQ=255)
                                          cTTTACTGTTGAATTACGTGAGGATTGTGTTATTATATAAATCGTAATAATTACGATTTGATAAAAAGTGAGGTAACTATATATGGCTaat                                                  <  2:3034332/91‑2 (MQ=255)
                                           tttACTGTTGAATTATGTGAGGATTGTGTTATTATATAAATCGTAATAATTACGATTTGATAAAAAGTGAGGTAACTATATATGGCTAAGa                                                 >  2:1925270/1‑91 (MQ=255)
                                           tttACTGTTGAATTATGCGAGGATTGTGTTATTATATAAATCGTAATAATTACGATTTGATAAAAAGTGAGGTAACTATATATGGCTAAGa                                                 <  2:2587077/91‑1 (MQ=255)
                                             tACTGTTGAATTATGTGAGGATTGTGTTATTATATAAATCGTAATAATTACGATTTGATAAAAAGTGAGGTAACTATATATGGCTAAGaaa                                               <  1:1303127/91‑1 (MQ=255)
                                              aCTGTTGAATTATGTGAGGATTGTGTTATTATATAAATCGTAATAATTACGATTTGATAAAAAGTGAGGTAACTATATATGGCTAAGAAAt                                              <  1:3344432/91‑1 (MQ=255)
                                                tGTTGAATTATGTGAGGATTGTGTTATTATATAAATCGTAATAATTACGATTTGATAAAAAGTGAGGTAACTATATATGGCTAAGAAATCt                                            <  1:3744692/91‑1 (MQ=255)
                                                tGTTGAATTATGTGAGGATTGTGTTATTATATAAATCGTAATAATTACGATTTGATAAAAAGTGAGGTAACTATATATGGCTAAGAAATCt                                            >  1:2751998/1‑91 (MQ=255)
                                                 gTTGAATTATGTGAGGATTGTGTTATTATATAAATCGTAATAATTACGATTTGATAAAAAGTGAGGTAACTATATATGGCTAAGAAATCTa                                           <  2:2371601/91‑1 (MQ=255)
                                                 gTTGAATTATGAGAGGATCGGGTTATTATCTAAATCGTAATAATTACGATTTGCTAAAAAGTGAGGGATCTATATATGGCTAAGACATCtc                                           >  2:1071356/1‑90 (MQ=255)
                                                  ttGAATTATGTGAGGATTGTGTTATTATATAAATCGTAATAATTACGATTTGATAAAAAGTGAGGTAACTATATATGGCTAAGAAATCTaa                                          <  1:680914/91‑1 (MQ=255)
                                                   tGAATTATGTGAGGATTGTGTTATTATATAAATCGTAATAATTACGATTTGATAAAAAGTGAGGTAACTATATATGGCTAAGAAATCTaaa                                         >  2:973212/1‑91 (MQ=255)
                                                    gAATTATGTGAGGATTGTGTTATTATATAAATCGTAATAATTACGATTTGATAAAAAGTGAGGTAACTATATATGGCTAAGAAATCTaaaa                                        <  1:3447585/91‑1 (MQ=255)
                                                    gAATTATGTGAGGATTGTGTTATTATATAAATCGTAATAATTACGATTTGATAAAAAGTGAGGTAACTATATATGGCTAAGAAATCTaaaa                                        <  2:2968021/91‑1 (MQ=255)
                                                       ttATGTGAGGATTGTGTTATTATATAAATCGTAATAATTACGATTTGATAAAAAGTGAGGTAACTATATATGGCTAAGAAATCTAAAATAg                                     <  1:1453763/91‑1 (MQ=255)
                                                           gtgAGGATTGTGTTATTATATAAATCGTAATAATTACGATTTGATAAAAAGTGAGGTAACTATATATGGCTAAGAAATCTAAAATAGCaaa                                 >  2:1716976/1‑91 (MQ=255)
                                                           gtgAGGATTGTGTTATTATATAAATCGTAATAATTACGATTTGATAAAAAGTGAGGTAACTATATATGGCTAAGAAATCTAAAATAGCaaa                                 <  2:1699722/91‑1 (MQ=255)
                                                             gAGGATTGTGTTATTATATAAATCGTAATAATTACGATTTGATAAAAAGTGAGGTAACTATATATGGCTAAGAAATCTAAAATAGCAAAag                               <  2:528913/91‑1 (MQ=255)
                                                             gAGGATTGTGTTATTATATAAATCGTAATAATTACGATTTGATAAAAAGTGAGGTAACTATATATGGCTAAGAAATCTAAAATAGCAAAag                               <  1:2118990/91‑1 (MQ=255)
                                                             gAGGATTGTGTTACTATATAAATCGTAATAATTACGATTTGATGAAAAGTGAGGTAACTATATATGGCTAAGAAATCTAAAATAGCAAAag                               <  1:3962090/91‑1 (MQ=255)
                                                               ggATTGTGTTATTATATAAATCGTAATAATTACGATTTGATAAAAAGTGAGGTAACTATATATGGCTAAGAAATCTAAAATAGCAAAagag                             >  2:224224/1‑91 (MQ=255)
                                                               ggATTGTGTTATTATATAAATCGTAATAATTACGATTTGATAAAAAGTGAGGTAACTATATATGGCTAAGAAATCTAAAATAGCAAAagag                             <  2:3219206/91‑1 (MQ=255)
                                                               ggATTGTGTTATTATATAAATCGTAATAATTACGATTTGATAAAAAGTGAGGTAACTATATATGGCTAAGAAATCTAAAATAGCAAAagag                             >  1:1427599/1‑91 (MQ=255)
                                                               ggATTGTGTTATTATATAAATCGTAATAATTACGATTTGATAAAAAGTGAGGTAACTATATATGGCTAAGAAATCTAAAATAGCAAAagag                             >  2:3958598/1‑91 (MQ=255)
                                                                gATTGTGTTATTATATAAATCGTAATAATTACGATTTGATAAAAAGTGAGGTAACTATATATGGCTAAGAAATCTAAAATAGCAAAagaga                            <  2:892185/91‑1 (MQ=255)
                                                                  ttGTGTTATTATATAAATCGTAATAATTACGATTTGATAAAAAGTGAGGTAACTATATATGGCTAAGAAATCTAAAATAGCAAAAgagaga                          <  2:2673183/91‑1 (MQ=255)
                                                                   tgtgTTATTATATAAATCGTAATAATTACGATTTGATAAAAAGTGAGGTAACTATATATGGCTAAGAAATCTAAAATAGCAAAAGAGAGaa                         <  1:117738/91‑1 (MQ=255)
                                                                    gtgtTATTATATAAATCGTAATAATTACGATTTGATAAAAAGTGAGGTAAGTATATATGGCTAAGAAATATAAAATAGCAAAAGAGAGaaa                        >  2:2473466/1‑91 (MQ=255)
                                                                    gtgtTATTATATAAATCGTAATAATTACGATTTGATAAAAAGTGAGGTAACTATATATGGCTAAGAAATCTAAAATAGCAAAAGAGAGaaa                        <  1:2591452/91‑1 (MQ=255)
                                                                     tgtTATTATATAAATCGTAATAATTACGATTTGATAAAAAGTGAGGTAACTATATATGGCTAAGAAATCTAAAATAGCAAAAGAGAGaaaa                       >  2:2996392/1‑91 (MQ=255)
                                                                      gtTATTATATAAATCGTAATAATTACGATTTGATAAAAAGTGAGGTAACTATATATGGCTAAGAAATCTAAAATAGCAAAAGAGAGaaaaa                      <  1:51913/91‑1 (MQ=255)
                                                                      gtTATTATATAAATCGTAATAATTACGATTTGATAAAAAGTGAGGTAACTATATATGGCTAAGAAATCTAAAATAGCAAAAGAGAGaaaaa                      <  2:403062/91‑1 (MQ=255)
                                                                         attatATAAATCGTAATAATTACGATTTGATAAAAAGTGAGGTAACTATATATGGCTAAGAAATCTAAAATAGCAAAAGAGAGAAAaagag                   <  1:1109074/91‑1 (MQ=255)
                                                                         attatATAAATCGTAATAATTACGATTTGATAAAAAGTGAGGTAACTATATATGGCTAAGAAATCTAAAATAGCAAAAGAGAGAAAaagag                   >  1:3660534/1‑91 (MQ=255)
                                                                         attatATAAATCGTAATAACTACGATTTGATAAAAAGTGAGGTAACTATATATGGCAAAGAAATCTAAAATAGCAAAAGAGAGAAAaagag                   >  2:1023205/1‑91 (MQ=255)
                                                                          ttatATAAATCGTAATAATTACGATTTGATAAAAAGTGAGGTAACTATATATGGCTAAGAAATCTAAAATAGCAAAAGAGAGAAAaagaga                  <  2:3305747/91‑1 (MQ=255)
                                                                          ttatATAAATCGTAATAATTACGATTTGATAAAAAGTGAGGTAACTATATATGGCTAAGAAATCTAAAATAGCAAAAGAGAGAAAaagaga                  <  2:1647019/91‑1 (MQ=255)
                                                                          ttatATAAATCGTAATAATTACGATTTGATAAAAAGTGAGGTAACTATATATGGCTAAGAAATCTAAAATAGCAAAAGAGAGAAAaagaga                  <  2:814619/91‑1 (MQ=255)
                                                                          ttatATAAATCGTAATAATTAAGATTGGATAAAAAGTGAGGTAACTATATATGGCTAAGAAAGCTAAAATAGCAAAAGAGAGAAAaagaga                  <  2:3646418/91‑1 (MQ=255)
                                                                            atataAATCGTAATAATTACGATTTGATAAAAAGTGAGGTAACTATATATGGCTAAGAAATCTAAAATAGCAAAAGAGAGAAAaagagaag                >  1:2364265/1‑91 (MQ=255)
                                                                            atataAATCGTAATAATTACGATTTGATAAAAAGTGAGGTAACTATATATGGCTAAGAAATCTAAAATAGCAAAAGAGAGAAAaagagaag                >  2:1140081/1‑91 (MQ=255)
                                                                               taAATCGTAATAATTACGATTTGATAAAAAGTGAGGTAACTATATATGGCTAAGAAATCTAAAATAGCAAAAGAGAGAAAAAGAGAAGAGt             <  1:152844/91‑1 (MQ=255)
                                                                                  aTCGTAATAATTACGATTTGATAAAAAGTGAGGTAACTATATATGGCTAAGAAATCTAAAATAGCAAAAGAGAGAAAAAGAGAAGAGTTAg          <  2:3481703/91‑1 (MQ=255)
                                                                                  aTCGTAATAATTACGATTTGATAAAAAGTGAGGTAACTATATATGGCTAAGAAATCTAAAATAGCAAAAGAGAGAAAAAGAGAAGAGTTAg          >  2:2292606/1‑91 (MQ=255)
                                                                                   tCGTAATAATTACGATTTGATAAAAAGTGAGGTAACTATATATGGCTAAGAAATCTAAAATAGCAAAAGAGAGAAAAAGAGAAGAGTTAGt         >  1:2549334/1‑91 (MQ=255)
                                                                                   tCGTAATAATTACGATTTGATAAAAAGTGAGGTAACTATATATGGCTAAGAAATCTAAAATAGCAAAAGAGAGAAAAAGAGAAGAGTTAGt         <  2:3118146/91‑1 (MQ=255)
                                                                                    cGTAATAATTACGATTTGATAAAAAGTGAGGTAACTATATATGGCTAAGAAATCTAAAATAGCAAAAGAGAGAAAAAGAGAAGAGTTAGta        >  1:3308322/1‑91 (MQ=255)
                                                                                      taattaTTACGATTTGATAAAAAGTGAGGTAACTATATATGGCTAAGAAATCTAAAATAGCAAAAGAGAGAAAAAGAGAAGAGTTAGtaaa      >  2:1812732/1‑91 (MQ=255)
                                                                                      taataaTTACGATTTGATAAAAAGTGAGGTAACTATATATGGCTAAGAAATCTAAAATAGCAAAAGAGAGAAAAAGAGAAGAGTTAGtaaa      <  1:2281736/91‑1 (MQ=255)
                                                                                       aataatTACGATTTGATAAAAAGTGAGGTAACTATATATGGCTAAGAAATCTAAAATAGCAAAAGAGAGAAAAAGAGAAGAGTTAGtaaat     >  2:1441682/1‑91 (MQ=255)
                                                                                         taatTACGATTTGATAAAAAGTGAGGTAACTATATATGGCTAAGAAATCTAAAATAGCAAAAGAGAGAAAAAGAGAAGAGTTAGtaaataa   >  2:2739867/1‑91 (MQ=255)
                                                                                         taaaTACGATTTGATAAAAGTGAGGTAACTATATATGGCTAAGAAATCTAAAATAGCAAAAGAGAGAAAAAGAGAAGAGTTAGcaaataaa  >  2:3501714/1‑91 (MQ=255)
                                                                                          |                                                                                          
TAAATATGATGGTTTAGAAAATGATTTAACAAAGAAATGAAACTTTACTGTTGAATTATGTGAGGATTGTGTTATTATATAAATCGTAATAATTACGATTTGATAAAAAGTGAGGTAACTATATATGGCTAAGAAATCTAAAATAGCAAAAGAGAGAAAAAGAGAAGAGTTAGTAAATAAA  >  USA300TCH1516_ALE/1366031‑1366211

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