Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 1,435,073 G→A 5.3% G144S (GGT→AGT) ‡ dapH → 2,3,4,5‑tetrahydropyridine‑2,6‑dicarboxylate N‑acetyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE1,435,0730GA 5.3% 829.9 / 6.7 269G144S (GGT→AGT) ‡dapH2,3,4,5‑tetrahydropyridine‑2,6‑dicarboxylate N‑acetyltransferase
Reads supporting (aligned to +/- strand):  ref base G (82/169);  new base A (16/2);  total (98/171)
Fisher's exact test for biased strand distribution p-value = 3.07e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.92e-01

Only 100 of 276 total aligned reads displayed.

TGATGGGCGCAACAATTAATATTGGCGCAGTCGTTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAACCCCCTAGTGCTTCACCGGTTATAATCGAGGATGATGTAT  >  USA300TCH1516_ALE/1434984‑1435163
                                                                                         |                                                                                          
tgatgGGCGCAAAAATTAATATTGGCGCAGTCGTTGGCGAGGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTgg                                                                                           <  2:6280235/91‑1 (MQ=255)
 gatgGGCGCAACAATTAATATTGGCGCAGTCGTTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGt                                                                                          <  1:1555099/91‑1 (MQ=255)
 gatgGGCGAAACAATTAATATTGGCGCAGTCGTTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGt                                                                                          <  2:2039089/91‑1 (MQ=255)
    gggCGCAACAATTAATATTGGCGCAGTCGTTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTaaa                                                                                       <  2:4992812/91‑1 (MQ=255)
    gggCGCAACAATTAATATTGGCGCAGTCGTTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTAGTaaa                                                                                       <  1:6643/91‑1 (MQ=255)
        gcAACAATTAATATTGGCGCAGTCGTTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGGTACAACTGGTAAAAATg                                                                                   >  2:2334034/1‑91 (MQ=255)
        gcAACAATTAATATTGGCGCAGTCGTTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATg                                                                                   <  1:5304791/91‑1 (MQ=255)
        gcAACAATTAATATTGGCGCAGTCGTTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATg                                                                                   <  2:3513246/91‑1 (MQ=255)
         caacaaTTAATATTGGCGCAGTCGTTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGt                                                                                  <  2:5255673/91‑1 (MQ=255)
          aacaaTTAATATTGGCGCAGTCGTTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGGTACAACTGGTAAAAATGTa                                                                                 >  2:6076422/1‑91 (MQ=255)
           acaaTTAATATTGGCGCAGTCGTTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTAc                                                                                >  1:3571138/1‑91 (MQ=255)
           acaaTTAATATTGGCGCAGTCGTTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTAc                                                                                <  2:2874795/91‑1 (MQ=255)
           acaaTTAATATTGGCGCAGTCGTTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTAc                                                                                >  1:3782672/1‑91 (MQ=255)
            caaTTAATATTGGCGCAGTCGTTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACa                                                                               >  1:279592/1‑91 (MQ=255)
              aTTAATATTGGCGCAGTCGTTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATg                                                                             <  1:2984370/91‑1 (MQ=255)
              aTTAATATTGGCGCAGTCGTTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATg                                                                             <  2:4067882/91‑1 (MQ=255)
              aTTAATATTGGCGCAGTCGTTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATg                                                                             >  1:4565857/1‑91 (MQ=255)
               ttAATATTGGCGCAGTCGTTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACCACTGGTAAAAATGTACATGt                                                                            >  2:6837964/1‑91 (MQ=255)
               ttAATATTGGCGCAGTCGTTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGt                                                                            <  2:3526064/91‑1 (MQ=255)
                tAATATTGGCGCAGTCGTTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGCGCTACCACTGGTAAAAATGTACATGTa                                                                           >  2:759938/1‑91 (MQ=255)
                 aaTATTGGCGCAGTCGTTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAg                                                                          >  1:1054134/1‑91 (MQ=255)
                 aaTATTGGCGCAGTCGTTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAg                                                                          >  1:6683438/1‑91 (MQ=255)
                   tatTGGCGCAGTCGTTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAggg                                                                        <  2:4631376/91‑1 (MQ=255)
                        gcgcAGTCGTTGGCGAAGGTACGATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTgg                                                                   <  2:1738697/91‑1 (MQ=255)
                        gcgcAGTCGTTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTgg                                                                   <  1:6169678/91‑1 (MQ=255)
                        gcgcAGTCGTTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTgg                                                                   <  2:3615906/91‑1 (MQ=255)
                        gcgcAGTCGTTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTgg                                                                   <  2:2731824/91‑1 (MQ=255)
                         cgcAGGCGTTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGgc                                                                  <  2:5613095/91‑1 (MQ=255)
                          gcAGTCGTTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGgcg                                                                 <  1:1363424/91‑1 (MQ=255)
                           cgggcGTTGGCGAAGGTATAATGATTGATATGAATGCTGCTCTCGGTGGTGGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGgcgc                                                                <  2:2753477/87‑1 (MQ=255)
                             gTCGTTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAg                                                              >  2:3114872/1‑91 (MQ=255)
                             gTCGTTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAg                                                              <  1:6860268/91‑1 (MQ=255)
                              tCGTTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGt                                                             <  1:3467714/91‑1 (MQ=255)
                               cGTTGGCGAAGGTACAATGATTGATATGAATGCTGCTCTCGGGGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTa                                                            <  2:4348401/91‑1 (MQ=255)
                               cGTTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTa                                                            <  1:3816650/91‑1 (MQ=255)
                                 ttGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTAtt                                                          <  1:6930784/91‑1 (MQ=255)
                                 ttGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTAtt                                                          >  1:4377976/1‑91 (MQ=255)
                                 ttGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGCGGTACCACTGGTAAAAAAGTACATGTAGGGGCTGGCGCAGTAtt                                                          >  2:490190/1‑91 (MQ=255)
                                 ttGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGCCGTATCATTAAAAAAAAAATACATGTAGGGGCGGGCGCAGTAtt                                                          >  2:6550136/1‑91 (MQ=255)
                                   ggCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGCCGTACCACTAGAAAAAATGTACATGTAGGGGCTGGCGCAGTATTAg                                                        >  2:1269231/1‑91 (MQ=255)
                                     cGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCa                                                      >  1:1220937/1‑91 (MQ=255)
                                        aGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGgt                                                   <  2:8469/91‑1 (MQ=255)
                                         ggTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGgtg                                                  >  1:2816247/1‑91 (MQ=255)
                                          ttACAATAATTAATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGCGCTGGCGCAGTATTAGCAGgtgt                                                 <  2:3238530/90‑1 (MQ=255)
                                          gTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGgtgt                                                 <  1:1975505/91‑1 (MQ=255)
                                          gTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGgtgt                                                 <  2:1445674/91‑1 (MQ=255)
                                            aCAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGa                                               <  1:6215147/91‑1 (MQ=255)
                                            aCAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGa                                               <  1:3134478/91‑1 (MQ=255)
                                             cAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGAt                                              <  2:5945009/91‑1 (MQ=255)
                                             cAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGAt                                              <  1:3177417/91‑1 (MQ=255)
                                                 gattgatATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGaa                                          <  1:1639821/91‑1 (MQ=255)
                                                 gattgatATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGaa                                          <  2:2092327/91‑1 (MQ=255)
                                                  attgatATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTGGCAGGTGTGATTGAAc                                         >  1:1931159/1‑91 (MQ=255)
                                                   ttggtATGAATGCTACTCTCGGTGGTCGTGCGACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAAcc                                        <  1:6504304/91‑1 (MQ=255)
                                                   ttgatATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAAcc                                        <  1:1000535/91‑1 (MQ=255)
                                                     gatATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAAcccc                                      <  2:5400292/91‑1 (MQ=255)
                                                     gatATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAAcccc                                      <  2:1825719/91‑1 (MQ=255)
                                                       tatGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAACCCCCt                                    >  2:3891519/1‑91 (MQ=255)
                                                       tatGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAACCCCCt                                    >  1:3673549/1‑91 (MQ=255)
                                                       tatGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAACCCCCt                                    >  1:5248279/1‑91 (MQ=255)
                                                         tGAATGCTACTCTCGGTGGTCGCCGTACCATTGGAAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAACCCCCTAg                                  >  2:265063/1‑91 (MQ=255)
                                                          gAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAACCCCCTAGt                                 <  1:2280713/91‑1 (MQ=255)
                                                           aaTGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAACCCCCTAGTg                                <  2:2423853/91‑1 (MQ=255)
                                                           aaTGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAACCCCCTAGTg                                <  1:5539532/91‑1 (MQ=255)
                                                           aaTGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAACCCCCTAGTg                                <  2:1149717/91‑1 (MQ=255)
                                                              gCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAACCCCCTAGTGCtt                             <  1:2178832/91‑1 (MQ=255)
                                                                tACTCTCGGTGGTCGTCGTACCATTAAAAAAAAAGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAACCCCCTAGTGCTTCa                           >  2:3367989/1‑91 (MQ=255)
                                                                 aCTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAACCCCCTAGTGCTTCAc                          <  1:5735736/91‑1 (MQ=255)
                                                                  ctctCGGTGGTCGTGCTACAATTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAACCCCCTAGTGCTTCAcc                         >  2:4691403/1‑91 (MQ=255)
                                                                  ctctCGGTGGTCGCGGTACCATTAGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAACCCCCTAGTGCTTCAcc                         >  2:185078/1‑91 (MQ=255)
                                                                   tctcGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAACCCCCTAGTGCTTCACCg                        >  2:1497362/1‑91 (MQ=255)
                                                                   tctcGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAACCCCCTAGTGCTTCACCg                        <  1:2635264/91‑1 (MQ=255)
                                                                   tctcGGTGGTCGTGCAACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAACCCCCTAGTGCTTCACCg                        <  2:1013955/91‑1 (MQ=255)
                                                                    ctcGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAACCCCCTAGTGCTTCACCgg                       >  1:5170150/1‑91 (MQ=255)
                                                                    ctaggtgggCGTGCTACAACAGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAACCCCCTAGTGCTTCACCgg                       <  2:679050/88‑1 (MQ=255)
                                                                     tcGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAACCCCCTAGTGCTTCACCGGt                      <  1:3891519/91‑1 (MQ=255)
                                                                       ggtggtCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAACCCCCTAGTGCTTCACCGGTta                    >  1:4473386/1‑91 (MQ=255)
                                                                         tggtCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAACCCCCTAGTGCTTCACCGGTtata                  >  1:541486/1‑91 (MQ=255)
                                                                         tggtCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAACCCCCTAGTGCTTCACCGGTtata                  >  1:4224640/1‑91 (MQ=255)
                                                                         tggtCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAACCCCCTAGTGCTTCACCGGTtata                  >  1:2736098/1‑91 (MQ=255)
                                                                           gtCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAACCCCCTAGTGCTTCACCGGTTATAAt                <  1:617983/91‑1 (MQ=255)
                                                                            tCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAACCCCCTAGTGCTTCACCGGTTATAATc               >  1:4977733/1‑91 (MQ=255)
                                                                            tCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAACCCCCTAGTGCTTCACCGGTTATAATc               <  1:4333998/91‑1 (MQ=255)
                                                                             cGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAACCCCCTAGTGCTTCACCGGTTATAATCg              >  1:4694634/1‑91 (MQ=255)
                                                                               tGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAACCCCCTAGTGCTTCACCGGTTATAATCGAg            <  1:2379019/91‑1 (MQ=255)
                                                                                gCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAACCCCCTAGTGCTTCACCGGTTATAATCGAgg           <  1:863587/91‑1 (MQ=255)
                                                                                gCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAACCCCCTAGTGCTTCACCGGTTATAATCGAgg           <  2:5775668/91‑1 (MQ=255)
                                                                                gCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAACCCCCTAGTGCTTCACCGGGTATAATCGAgg           <  2:6947665/91‑1 (MQ=255)
                                                                                gCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGAAGGTGTGATTGAACTCCCTAGTGCTTCACCGGTTATAATCGAgg           <  2:5150037/91‑1 (MQ=255)
                                                                                 cTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAACCCCCTAGTGCTTCACCGGTTATAATCGAGga          >  2:2508374/1‑91 (MQ=255)
                                                                                  tACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAACCCCCTAGTGCTTCACCGGTTATAATCGAGgat         >  1:72229/1‑91 (MQ=255)
                                                                                    cAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAACCCCCTAGTGCTTCACCGGTTATAATCGAGgatga       <  2:3257841/91‑1 (MQ=255)
                                                                                    cAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAACCCCCTAGTGCTTCACCGGTTATAATCGAGgatga       <  2:1628670/91‑1 (MQ=255)
                                                                                      aCTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAACCCCCTAGTGCTTCACCGGTTATAATCGAGGatgatg     <  2:5459162/91‑1 (MQ=255)
                                                                                      aCTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAACCCCCTAGTGCTTCACCGGTTATAATCGAGGatgatg     <  2:3038798/91‑1 (MQ=255)
                                                                                      aCTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAACCCCCTAGTGCTTCACCGGTTATAATCGAGGatgatg     <  1:3389075/91‑1 (MQ=255)
                                                                                        tGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAACCCCCTAGTGCTTCACCGGTTATATTCGAGGATGATGTa   >  1:5805165/1‑91 (MQ=255)
                                                                                        tGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAACCCCCTAGTGCTTCACCGGTTATAATCGAGGATGATGTa   <  2:1210308/91‑1 (MQ=255)
                                                                                         ggTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAACCCCCTAGTGCTTCACCGGTTATAATCGAGGATGATGTAt  <  2:4197375/91‑1 (MQ=255)
                                                                                         ggTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAACCCCCTAGTGCTTCACCGGTTATAATCGAGGATGATGTAt  >  1:3479157/1‑91 (MQ=255)
                                                                                         |                                                                                          
TGATGGGCGCAACAATTAATATTGGCGCAGTCGTTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAACCCCCTAGTGCTTCACCGGTTATAATCGAGGATGATGTAT  >  USA300TCH1516_ALE/1434984‑1435163

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