Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 1,273,930 A→G 11.7% T111T (ACA→ACG codY → GTP‑sensing transcriptional pleiotropic repressor CodY

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE1,273,9300AG11.7% 580.9 / 57.2 212T111T (ACA→ACGcodYGTP‑sensing transcriptional pleiotropic repressor CodY
Reads supporting (aligned to +/- strand):  ref base A (93/94);  new base G (20/5);  total (113/99)
Fisher's exact test for biased strand distribution p-value = 5.02e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.39e-01

Only 100 of 220 total aligned reads displayed.

AGTTAAACAAACAGAATCAAATATTGATATCGACAATGTATTAACAGTATTCCCACCTGAAAACAGAGAATTATTCATAGATAGTCGTACAACTATCTTCCCAATTTTAGGTGGAGGGGAAAGATTAGGTACATTAGTACTTGGTCGAGTACATGATGATTTTAATGAAAATGATTTGGTA  >  USA300TCH1516_ALE/1273840‑1274020
                                                                                          |                                                                                          
aGTTAAACAAACAGAATCAAATATTGATATCGACAATGTATTAACAGTATTCCCACCTGAAAACAGAGAATTATTCATAGATAGTCGTACa                                                                                            <  2:484786/91‑1 (MQ=255)
aGTTAAACAAACAGAATCAAATATTGATATCGACAATGTATTAACAGTATTCCCACCTGAAAACAGAGAATTATTCATAGATAGTCGTACa                                                                                            <  2:3152491/91‑1 (MQ=255)
    aaacaaacAGAATCAAATATTGATATCGACAATGTATTAACAGTATTCCCACCTGAAAACAGAGAATTATTCATAGATAGTCGTACAACTa                                                                                        >  2:2738640/1‑91 (MQ=255)
     aacaaacaGAATCAAATATTGATATCGACAATGTATTAACAGTATTCCCACCTGAAAACAGAGAATTATTCATAGATAGTCGTACAACGAt                                                                                       <  1:3672100/91‑1 (MQ=255)
      acaaacaGAATCAAATATTGATATCGACAATGTATCAACAGTATTCCCACCTGAAAACAGAGAATTATTCATAGATAGTTGTACGACTATc                                                                                      <  1:3110550/91‑1 (MQ=255)
         aacaGAATCAAATATTGATATCGACAATGTATTAACAGTATTCCCACCTGAAAACAGAGAATTATTCATAGATAGTCGTACAACTATCTTc                                                                                   <  1:3019574/91‑1 (MQ=255)
          acaGAATCAAATATTGATATCGACAATGTATTAACAGTATTCCCACCTGAAAACAGAGAATTATTCATAGATAGTCGTACAACTATCTTcc                                                                                  >  1:2867433/1‑91 (MQ=255)
           caGAATCAAATATTGATATCGACAATGTATTAACAGTATTCCCACCTGAAAACAGAGAATTATTCATAGATAGTCGTACAACTATCTTccc                                                                                 <  2:1200009/91‑1 (MQ=255)
             gAATCAAATATTGATATCGACAATGTATTAACAGTATTCCCACCTGAAAACAGAGAATTATTCATAGATAGTCGTACAACTATCTTCCCaa                                                                               >  1:614157/1‑91 (MQ=255)
              aaTCAAATATTGATATCGACAATGTATTAACAGTATTCCCACCTGAAAACAGAGAATTATTCATAGATAGTCGTACAACTATCTTCCCAAt                                                                              <  1:2292798/91‑1 (MQ=255)
               aTCAAATATTGATATCGACAATGTATTAACAGTATTCCCACCTGAAAACAGAGAATTATTCATAGATAGTCGTACAACTATCTTCCCAAtt                                                                             >  2:2130547/1‑91 (MQ=255)
               aTCAAATATTGATATCGACAATGTATTAACAGTATTCCCACCTGAAAACAGAGAATTATTCATAGATAGTCGTACAACTATCTTCCCAAtt                                                                             >  2:1843113/1‑91 (MQ=255)
                tCAAATATTGATATCGACAATGTATTAACAGTATTCCCACCTGAAAACAGAGAATTATTCATAGATAGTCGTACAACTATCTTCCCAAttt                                                                            >  2:1097880/1‑91 (MQ=255)
                 cAAATATTGATATCGACAATGTATTAACAGTATTCCCACCTGAAAACAGAGAATTATTCATAGATAGTCGTACAACTATCTTCCCAAtttt                                                                           <  1:1149847/91‑1 (MQ=255)
                 cAAATATTGATATCGACAATGTATTAACAGTATTCCCACCTGAAAACAGAGAATTATTCATAGATAGTCGTACAACTATCTTCCCAAtttt                                                                           <  1:3168068/91‑1 (MQ=255)
                 cAAATATTGATATCGACAATGTATTAACAGTATTCCCACCTGAAAACAGAGAATTATTCATAGATAGTCGTACAACTATCTTCCCAAtttt                                                                           <  2:56551/91‑1 (MQ=255)
                  aaaTATTGATATCGACAATGTATTAACAGTATTCCCACCTGAAAACAGAGAATTATTCATAGATAGTCGTACAACTATCTTCCCAATTTTa                                                                          <  2:1675871/91‑1 (MQ=255)
                   aaTATTGATATCGACAATGTATTAACAGTATTCCCACCTGAAAACAGAGAATTATTCATAGATAGTCGTACAACTATCTTCCCAATTTTAg                                                                         >  1:2051183/1‑91 (MQ=255)
                        tGATATCGACAATGTATTAACAGTATTCCCACCTGAAAACAGAGAATTATTCATAGATAGTCGTACAACTATCTTCCCAATTTTAGGTGGa                                                                    >  1:33301/1‑91 (MQ=255)
                         gATATCGACAATGTATTAACAGTATTCCCACCTGAAAACAGAGAATTATTCATAGATAGTCGTACAACTATCTTCCCAATTTTAGGTGGAg                                                                   <  1:2684248/91‑1 (MQ=255)
                          atatCGACAATGTATTAACAGTATTCCCACCTGAAAACAGAGAATTATTCATAGATAGTCGTACAACTATCTTCCCAATTTTAGGTGGAgg                                                                  <  2:3468693/91‑1 (MQ=255)
                          atatCGACAATGTATTAACAGTATTCCCACCTGAAAACAGAGAATTATTCATAGATAGTCGTACAACTATCTTCCCAATTTTAGGTGGAgg                                                                  <  1:3540882/91‑1 (MQ=255)
                          atatCGACAATGTATTAACAGTATTCCCACCTGAAAACAGAGAATTATTCATAGATAGTCGTACAACTATCTTCCCAATTTTAGGTGGAgg                                                                  <  1:389611/91‑1 (MQ=255)
                          atatCGACAATGTATTAACAGTATTCCCACCTGAAAACAGAGAATTATTCATAGATAGTCGTACAACTATCTTCCCAATTTTAGGTGGAgg                                                                  <  1:508691/91‑1 (MQ=255)
                           tatCGACAATGTATTAACAGTATTCCCACCTGAAAACAGAGAATTATTCATAGATAGTCGTACAACTATCTTCCCAATTTTAGGTGGAggg                                                                 >  1:1322485/1‑91 (MQ=255)
                            atCGACAATGTATTAACAGTATTCCCACCTGAAAACAGAGAATTATTCATAGATAGTCGTACAACTATCTTCCCAATTTTAGGTGGAgggg                                                                >  1:2269811/1‑91 (MQ=255)
                            atCGACAATGTATTAACAGTATTCCCACCTGAAAACAGAGAATTATTCATAGATAGTCGTACAACTATCTTCCCAATTTTAGGTGGAgggg                                                                <  2:2874430/91‑1 (MQ=255)
                            atCGACAATGTATTAACAGTATTCCCACCTGAAAACAGAGAATTATTCATAGATAGTCGTACAACTATCTTCCCAATTTTAGGTGGAgggg                                                                >  1:189504/1‑91 (MQ=255)
                            atCGACAATGTATTAACAGTATTCCCACCTGAAAACAGAGAATTATTCATAGATAGTCGTACAACTATCTTCCCAATTTTAGGTGGAgggg                                                                >  1:3830103/1‑91 (MQ=255)
                             tCGACAATGTATTAACAGTATTCCCACCTGAAAACAGAGAATTATTCATAGATAGTCGTACAACTATCTTCCCAATTTTAGGTGGAGGGGa                                                               <  1:3495495/91‑1 (MQ=255)
                                  aaTGTATTAACAGTATTCCCACCTGAAAACAGAGAATTATTCATAGATAGTCGTACAACTATCTTCCCAATTTTAGGTGGAGGGGAAAGAt                                                          <  2:1938377/91‑1 (MQ=255)
                                   aTGTATTAACAGTATTCCCACCTGAAAACAGAGAATTATTCATAGATAGTCGTACAAGTATCTTCCCAATTTTAGGTGGAGGGGAAAGAtt                                                         >  1:3448982/1‑91 (MQ=255)
                                   aTGTATTAACAGTATTCCCACCTGAAAACAGAGAATTATTCATAGATAGTCGTACAACTATCTTCCCAATTTTAGGTGGAGGGGAAAGAt                                                          <  1:2949919/90‑1 (MQ=255)
                                   aTGTATTAACAGTATTCCCACCTGAAAACAGAGAATTATTCATAGATAGTCGTACAACTATCTTCCCAATTTTAGGTGGAGGGGAAAGAt                                                          >  2:2949919/1‑90 (MQ=255)
                                     gTATTAACAGTATTCCCACCTGAAAACAGAGAATTATTCATAGATAGTCGTACAACTATCTTCCCAATTTTAGGTGGAGGGGAAAGATTAg                                                       >  2:1730088/1‑91 (MQ=255)
                                     gTATTAACAGTATTCCCACCTGAAAACAGAGAATTATTCATAGATAGTCGTACAACTATCTTCCCAATTTTAGGTGGAGGGGAAAGATTAg                                                       >  2:1235683/1‑91 (MQ=255)
                                         tAACAGTATTCCCACCTGAAAACAGAGAATTATTCATAGATAGTCGTACAACTATCTTCCCAATTTTAGGTGGAGGGGAAAGATTAGGTAc                                                   >  2:1322360/1‑91 (MQ=255)
                                          aaCAGTATTCCCACCTGAAAACAGAGAATTATTCATAGATAGTCGTACAACTATCTTCCCAATTTTAGGTGGAGGGGAAAGATTAGGTACa                                                  >  1:923952/1‑91 (MQ=255)
                                            cAGTATTCCCACCTGAAAACAGAGAATTATTCATAGATAGTCGTACAACTATCTTCCCAATTTTAGGTGGAGGGGAAAGATTAGGTACAtt                                                <  1:1088036/91‑1 (MQ=255)
                                              gTATTCCCACCTGAAAACAGAGAATTATTCATAGATAGTCGTACAACTATCTTCCCAATTTTAGGTGGAGGGGAAAGATTAGGTACATTAg                                              <  2:3546964/91‑1 (MQ=255)
                                              gTATTCCCACCTGAAAACAGAGAATTATTCATAGATAGTCGTACAACTATCTTCCCAATTTTAGGTGGAGGGGAAAGATTAGGTACATTAg                                              >  1:642815/1‑91 (MQ=255)
                                               tATTCCCACCTGAAAACAGAGAATTATTCATAGATAGTCGTACAACTATCTTCCCAATTTTAGGTGGAGGGGAAAGATTAGGTACATTAGt                                             >  1:1459087/1‑91 (MQ=255)
                                               tATTCCCACCTGAAAACAGAGAATTATTCATAGATAGTCGTACAACTATCTTCCCAATTTTAGGTGGAGGGGAAAGATTAGGTACATTAGt                                             >  2:3517709/1‑91 (MQ=255)
                                               tATTCCCACCTGAAAACAGAGAATTATTCATAGATAGTCGTACAACTATCTTCCCAATTTTAGGTGGAGGGGAAAGATTAGGTACATTAGt                                             >  2:3757455/1‑91 (MQ=255)
                                                aTTCCCACCTGAAAACAGAGAATTATTCATAGATAGTCGTACAACTATCTTCCCAATTTTAGGTGGAGGGGAAAGATTAGGTACATTAGTa                                            <  2:794766/91‑1 (MQ=255)
                                                 ttCCCACCTGAAAACAGAGAATTATTCATAGATAGTCGTACAACTATCTTCCCAATTTTAGGTGGAGGGGAAAGATTAGGTACATTAGTAc                                           <  1:1117085/91‑1 (MQ=255)
                                                 ttCCCACCTGAAAACAGAGAATTATTCATAGATAGTCGTACAACTATCTTCCCAATTTTAGGTGGAGGGGAAAGATTAGGTACATTAGTAc                                           >  1:779147/1‑91 (MQ=255)
                                                  tCCCACCTGAAAACAGAGAATTATTCATAGATAGTCGTACAACTATCTTCCCAATTTTAGGTGGAGGGGAAAGATTAGGTACATTAGTACt                                          <  2:2456245/91‑1 (MQ=255)
                                                  tCCCACCTGAAAACAGAGAATTATTCATAGATAGTCGTACAACTATCTTCCCAATTTTAGGTGGAGGGGAAAGATTAGGTACATTAGTACt                                          <  1:1698221/91‑1 (MQ=255)
                                                   cccACCTGAAAACAGAGAATTATTCATAGATAGTCGTACAACTATCTTCCCAATTTTAGGTGGAGGGGAAAGATTAGGTACATTAGTACtt                                         >  1:989805/1‑91 (MQ=255)
                                                   cccACCTGAAAACAGAGAATTATTCATAGATAGTCGTACAACTATCTTCCCAATTTTAGGTGGAGGGGAAAGATTAGGTACATTAGTACtt                                         >  1:591021/1‑91 (MQ=255)
                                                   cccACCTGAAAACAGAGAATTATTCATAGATAGTCGTACAACTATCTTCCCAATTTTAGGTGGAGGGGAAAGATTAGGTACATTAGTACtt                                         >  1:2206407/1‑91 (MQ=255)
                                                   cccACCTGAAAACAGAGAATTATTCATAGATAGTCGTACAACTATCTTCCCAATTTTAGGTGGAGGGGAAAGATTAGGTACATTAGTACtt                                         <  1:1286217/91‑1 (MQ=255)
                                                    ccACCTGAAAACAGAGAATTATTCATAGATAGTCGTACAACTATCTTCCCAATTTTAGGTGGAGGGGAAAGATTAGGTACATTAGTACTTg                                        >  2:2205257/1‑91 (MQ=255)
                                                    ccACCTGAAAACAGAGAATTATTCATAGATAGTCGTACAACTATCTTCCCAATTTTAGGTGGAGGGGAAAGATTAGGTACATTAGTACTTg                                        >  1:2335767/1‑91 (MQ=255)
                                                    ccACCTGAAAACAGAGAATTATTCATAGATAGTCGTACAACTATCTTCCCAATTTTAGGTGGAGGGGAAAGATTAGGTACATTAGTACTTg                                        >  2:1888028/1‑91 (MQ=255)
                                                     cACCTGAAAACAGAGAATTATTCATAGATAGTCGTACAACTATCTTCCCAATTTTAGGTGGAGGGGAAAGATTAGGTACATTAGTACTTgg                                       >  1:2638973/1‑91 (MQ=255)
                                                     cACCTGAAAACAGAGAATTATTCATAGATAGTCGTACAACTATCTTCCCAATTTTAGGTGGAGGGGAAAGATTAGGTACATTAGTACTTgg                                       <  1:278686/91‑1 (MQ=255)
                                                      aCCTGAAAACAGAGAATTATTCATAGATAGTCGTACAACTATCTTCCCAATTTTAGGTGGAGGGGAAAGATTAGGTACATTAGTACTTGGt                                      >  1:3179058/1‑91 (MQ=255)
                                                          gAAAACAGAGAATTATTCATAGATAGTCGTACAACTATCTTCCCAATTTTAGGTGGAGGGGAAAGATTAGGTACATTAGTACTTGGTCGAg                                  <  2:2832690/91‑1 (MQ=255)
                                                          gAAAACAGAGAATTATTCATAGATAGTCGTACAACTATCTTCCCAATTTTAGGTGGAGGGGAAAGATTAGGTACATTAGTACTTGGTCGAg                                  <  1:2713727/91‑1 (MQ=255)
                                                           aaaaCAGAGAATTATTCATAGATAGTCGTACAACTATCTTCCCAATTTTAGGTGGAGGGGAAAGATTAGGTACATTAGTACTTGGTCGAGt                                 <  1:1349490/91‑1 (MQ=255)
                                                             aaCAGAGAATTATTCATAGATAGTCGTACAACTATCTTCCCAATTTTAGGTGGAGGGGAAAGATTAGGTACATTAGTACTTGGTCGAGTAc                               >  1:1413366/1‑91 (MQ=255)
                                                              aCAGAGAATTATTCATAGATAGTCGTACAACTATCTTCCCAATTTTAGGTGGAGGGGAAAGATTAGGTACATTAGTACTTGGTCGAGTACa                              >  2:136923/1‑91 (MQ=255)
                                                                agagAATTATTCATAGATAGTCGTACAACTATCTTCCCAATTTTAGGTGGAGGGGAAAGATTAGGTACATTAGTACTTGGTCGAGTACatg                            >  1:2310591/1‑91 (MQ=255)
                                                                 gagaATTATTCATAGATAGTCGTACAACTATCTTCCCAATTTTAGGTGGAGGGGAAAGATTAGGTACATTAGTACTTGGTCGAGTACatga                           >  2:410191/1‑91 (MQ=255)
                                                                     attattCATAGATAGTCGTACAACTATCTTCCCAATTTTAGGTGGAGGGGAAAGATTAGGTACATTAGTACTTGGTCGAGTACATgatgat                       <  1:1097880/91‑1 (MQ=255)
                                                                        attCATAGATAGTCGTACAACTATCTTCCCAATTTTAGGTGGAGGGGAAAGATTAGGTACATTAGTACTTGGTCGAGTACATGATGAtttt                    <  2:1933978/91‑1 (MQ=255)
                                                                         ttCATAGATAGTCGTACAACTATCTTCCCAATTTTAGGTGGAGGGGAAAGATTAGGTACATTAGTACTTGGTCGAGTACATGATGATTTTa                   <  2:1558631/91‑1 (MQ=255)
                                                                         ttCATAGATAGTCGTACAACTATCTTCCCAATTTTAGGTGGAGGGGAAAGATTAGGTACATTAGTACTTGGTCGAGTACATGATGATTTTa                   >  2:3421095/1‑91 (MQ=255)
                                                                         ttCATAGATAGTCGTACAACTATCTTCCCAATTTTAGGTGGAGGGGAAAGATTAGGTACATTAGTACTTGGTCGAGTACATGATGATTTTa                   >  2:694067/1‑91 (MQ=255)
                                                                          tCATAGATAGTCGTACAACTATCTTCCCAATTTTAGGTGGAGGGGAAAGATTAGGTACATTAGTACTTGGTCGAGTACATGATGATTTTaa                  <  2:2615598/91‑1 (MQ=255)
                                                                            gaagatagTTGTACGACTATCTTCCCAATTTTAGGTGGAGGGGAAAGATTAGGTACATTAGTACTTGGTCGAGTACATGATGATTTTAATg                >  2:215148/3‑91 (MQ=255)
                                                                            atagatagTCGTACAACTATCTTCCCAATTTTAGGTGGAGGGGAAAGATTAGGTACATTAGTACTTGGTCGAGTACATGATGATTTTAATg                <  2:1686046/91‑1 (MQ=255)
                                                                            atagatagTCGTACAACTATCTTCCCAATTTTAGGTGGAGGGGAAAGATTAGGTACATTAGTACTTGGTCGAGTACATGATGATTTTAATg                <  1:1459570/91‑1 (MQ=255)
                                                                            atagatagTCGTACAACTATCTTCCCAATTTTAGGTGGAGGGGAAAGATTAGGTACATTAGTACTTGGTCGAGTACATGATGATTTTAATg                <  2:3775494/91‑1 (MQ=255)
                                                                                atagTTGTACGACTATCTTCCCAATTTTAGGTGGAGGGGAAAGATTAGGTACATTAGTACTTGGTCGAGTACATGATGATTTTAATGaaaa            >  1:3296674/1‑91 (MQ=255)
                                                                                atagTCGTACAACTATCTTCCCAATTTTAGGTGGAGGGGAAAGATTAGGTACATTAGTACTTGGTCGAGTACATGATGATTTTAATGaaaa            >  2:1506791/1‑91 (MQ=255)
                                                                                atagTCGTACAACTATCTTCCCAATTTTAGGTGGAGGGGAAAGATTAGGTACATTAGTACTTGGTCGAGTACATGATGATTTTAATGaaaa            <  2:779147/91‑1 (MQ=255)
                                                                                 tagTTGTACGACTATCTTCCCAATTTTAGGTGGAGGGGAAAGATTAGGTACATTAGTACTTGGTCGAGTACATGATGATTTTAATGAAAAt           >  1:1378131/1‑91 (MQ=255)
                                                                                 tagTTGTACGACTATCTTCCCAATTTTAGGTGGAGGGGAAAGATTAGGTACATTAGTACTTGGTCGAGTACATGATGATTTTAATGAAAAt           >  1:1057982/1‑91 (MQ=255)
                                                                                 tagTTGTACGACTATCTTCCCAATTTTAGGTGGAGGGGAAAGATTAGGTACATTAGTACTTGGTCGAGTACATGATGATTTTAATGAAAAt           >  2:2551631/1‑91 (MQ=255)
                                                                                   gttgTACGACTATCTTCCCAATTTTAGGTGGAGGGGAAAGATTAGGTACATTAGTACTTGGTCGAGTACATGATGATTTTAATGAAAATGa         >  2:3051601/4‑91 (MQ=255)
                                                                                   gttgTACGACTATCTTCCCAATTTTAGGTGGAGGGGAAAGATTAGGTACATTAGTACTTGGTCGAGTACATGATGATTTTAATGAAAATGa         >  1:2492526/4‑91 (MQ=255)
                                                                                   gTCGTACAACTATCTTCCCAATTTTAGGTGGAGGGGAAAGATTAGGTACATTAGTACTTGGTCGAGTACATGATGATTTTAATGAAAATGa         <  2:1334479/91‑1 (MQ=255)
                                                                                   gTCGTACAACTATCTTCCCAATTTTAGGTGGAGGGGAAAGATTAGGTACATTAGTACTTGGTCGAGTACATGATGATTTTAATGAAAATGa         <  2:3865836/91‑1 (MQ=255)
                                                                                   gTCGTACAACTATCTTCCCAATTTTAGGTGGAGGGGAAAGATTAGGTACATTAGTACTTGGTCGAGTACATGATGATTTTAATGAAAATGa         >  2:2330657/1‑91 (MQ=255)
                                                                                    ttgTACGACTATCTTCCCAATTTTAGGTGGAGGGGAAAGATTAGGTACATTAGTACTTGGTCGAGTACATGATGATTTTAATGAAAATGAt        >  2:2728818/3‑91 (MQ=255)
                                                                                     tgTACGACTATCTTCCCAATTTTAGGTGGAGGGGAAAGATTAGGTACATTAGTACTTGGTCGAGTACATGATGATTTTAATGAAAATGAtt       >  2:3322039/2‑91 (MQ=255)
                                                                                     cGTACAACTATCTTCCCAATTTTAGGTGGAGGGGAAAGATTAGGTACATTAGTACTTGGTCGAGTACATGATGATTTTAATGAAAATGAtt       >  1:1771052/1‑91 (MQ=255)
                                                                                      gTACAACTATCTTCCCAATTTTAGGTGGAGGGGAAAGATTAGGTACATTAGTACTTGGTCGAGTACATGATGATTTTAATGAAAATGAttt      <  1:2212075/91‑1 (MQ=255)
                                                                                       tACAACTATCTTCCCAATTTTAGGTGGAGGGGAAAGATTAGGTACATTAGTACTTGGTTGAGTACATGATGATTTTAATGAAAATGATTTg     <  2:3803296/91‑1 (MQ=255)
                                                                                       tACAACTATCTTCCCAATTTTAGGTGGAGGGGAAAGATTAGGTACATTAGTACTTGGTCGAGTACATGATGATTTTAATGAAAATGATTTg     >  2:1580044/1‑91 (MQ=255)
                                                                                       tACAACTATCTTCCCAATTTTAGGTGGAGGGGAAAGATTAGGTACATTAGTACTTGGTCGAGTACATGATGATTTTAATGAAAATGATTTg     >  1:1026509/1‑91 (MQ=255)
                                                                                        aCAACTATCTTCCCAATTTTAGGTGGAGGGGAAAGATTAGGTACATTAGTACTTGGTCGAGTACATGATGATTTTAATGAAAATGATTTgg    >  1:1180113/1‑91 (MQ=255)
                                                                                        aCAACTATCTTCCCAATTTTAGGTGGAGGGGAAAGATTAGGTACATTAGTACTTGGTCGAGTACATGATGATTTTAATGAAAATGATTTgg    <  2:1137143/91‑1 (MQ=255)
                                                                                         cAACTATCTTCCCAATTTTAGGTGGAGGGGAAAGATTAGGTACATTAGTACTTGGTCGAGTACATGATGATTTTAATGAAAATGATTTGGt   >  2:2059081/1‑91 (MQ=255)
                                                                                         cAACTATCTTCCCAATTTTAGGTGGAGGGGAAAGATTAGGTACATTAGTACTTGGTCGAGTACATGATGATTTTAATGAAAATGATTTGGt   <  1:3510522/91‑1 (MQ=255)
                                                                                          aaCTATCTTCCCAATTTTAGGTGGAGGGGAAAGATTAGGTACATTAGTACTTGGTCGAGTACATGATGATTTTAATGAAAATGATTTGGTa  >  2:1717294/1‑91 (MQ=255)
                                                                                          aaCTATCTTCCCAATTTTAGGTGGAGGGGAAAGATTAGGTACATTAGTACTTGGTCGAGTACATGATGATTTTAATGAAAATGATTTGGTa  <  2:2650514/91‑1 (MQ=255)
                                                                                          |                                                                                          
AGTTAAACAAACAGAATCAAATATTGATATCGACAATGTATTAACAGTATTCCCACCTGAAAACAGAGAATTATTCATAGATAGTCGTACAACTATCTTCCCAATTTTAGGTGGAGGGGAAAGATTAGGTACATTAGTACTTGGTCGAGTACATGATGATTTTAATGAAAATGATTTGGTA  >  USA300TCH1516_ALE/1273840‑1274020

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: