Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 1,435,065 C→G 7.6% A141G (GCT→GGT)  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,0650CG 7.6% 714.7 / 9.1 216A141G (GCT→GGT) dapH2,3,4,5‑tetrahydropyridine‑2,6‑dicarboxylate N‑acetyltransferase
Reads supporting (aligned to +/- strand):  ref base C (66/131);  new base G (16/2);  total (83/133)
Fisher's exact test for biased strand distribution p-value = 5.37e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 224 total aligned reads displayed.

CTGTCGTTATGATGGGCGCAACAATTAATATTGGCGCAGTCGTTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAACCCCCTAGTGCTTCACCGGTTATAATCGAGGA  >  USA300TCH1516_ALE/1434975‑1435155
                                                                                          |                                                                                          
cTGTCGTTATGATGGGCGCAACAATTAATATTGGCGCAGTCGTTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGc                                                                                            >  1:504747/1‑91 (MQ=255)
cTGTCGTTATGATGGGCGCAACAATTAATATTGGCGCAGTCGTTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGc                                                                                            <  1:2851122/91‑1 (MQ=255)
cTGTCGTTATGATGGGCGCAACAATTAATATTGGCGCAGTCGTTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGc                                                                                            >  2:1858074/1‑91 (MQ=255)
 tGTCGTTATGATGGGCGCAACAATTAATATTGGCGCAGTCGTTGGCGAAGGTTCAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCt                                                                                           <  2:3305155/91‑1 (MQ=255)
     gTTATGATGGGCGCAACAATTAATATTGGCGCAGTCGTTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTAcca                                                                                       >  2:4242837/1‑89 (MQ=255)
     gTTATGATGGGCGCAACAATTAATATTGGCGCAGTCGTTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACaa                                                                                       >  1:4766576/1‑91 (MQ=255)
     gTTATGATGGGCGCAACAATTAATATTGGCGCAGTCGTTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACaa                                                                                       <  1:4920367/91‑1 (MQ=255)
      ttATGATGGGCGCAACAATTAATATTGGCGCAGTCGTTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAAc                                                                                      <  1:4056777/91‑1 (MQ=255)
      ttATGATGGGCGCAACAATTAATATTGGCGCAGTCGTTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAAc                                                                                      >  1:1791372/1‑91 (MQ=255)
       tATGATGGGCGCAACAATTAATATTGGCGCAGTCGTTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGCGCTACCACt                                                                                     >  2:192162/1‑91 (MQ=255)
         tgatgGGCGCAACAATTAATATTGGCGCAGTCGTTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTgg                                                                                   <  1:2226006/91‑1 (MQ=255)
           atgGGCGCAACAATTAATATTGGCGCAGTCGTTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTa                                                                                 >  2:5184762/1‑91 (MQ=255)
           atgGGCGCAACAATTAATATTGGCGCAGTCGTTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTa                                                                                 >  1:9092/1‑91 (MQ=255)
           atgGGCGCAACAATTAATATTGGCGCAGTCGTTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTCGTACCACtaaaa                                                                                 >  2:625953/1‑87 (MQ=255)
            tgGGCGCAACAATTAATATTGGCGCAGTCGTTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTaa                                                                                <  2:1205268/91‑1 (MQ=255)
             gggCGCAACAATTAATATTGGCGCAGTCGTTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTaaa                                                                               <  2:1393955/91‑1 (MQ=255)
              ggCGCAACAATTAATATTGGCGCAGTCGTTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTaaaa                                                                              >  1:4561027/1‑91 (MQ=255)
               gcgcAACAATTAATATTGGCGCAGTCGTTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTCGTACCACTGGaaaaaa                                                                             >  2:4908373/1‑91 (MQ=255)
                cgcAACAATTAATATTGGCGCAGTCGTTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGGTACCATTGGTaaaaaa                                                                            >  2:694356/1‑90 (MQ=255)
                cgcAACAATTAATATTGGCGCAGTCGTTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGCGCTACAACTGGaaaaaaa                                                                            >  2:3241124/1‑90 (MQ=255)
                  caacaaTTAATATTGGCGCAGTCGTTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACCACTGGTAAAAATGt                                                                          >  2:5334791/1‑91 (MQ=255)
                   aacaaTTAATATTGGCGCAGTCGTTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTa                                                                         >  1:509077/1‑91 (MQ=255)
                    acaaTTAATATTGGCGCAGTCGTTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTAc                                                                        <  2:1495000/91‑1 (MQ=255)
                    acaaTTAATATTGGCGCAGTCGTTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTAc                                                                        >  1:1011069/1‑91 (MQ=255)
                     caaTTAATATTGGCGCAGTCGTTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACa                                                                       <  1:5408415/91‑1 (MQ=255)
                      aaTTAATATTGGCGCAGTCGTTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACAt                                                                      <  2:85070/91‑1 (MQ=255)
                       tttAATGTTGGCGAGGTCGTTGGCGAAGGTTCAGTGATTGTTATGAATGCTGCTCTCTGTGGTCGTGCTACAACTGGTAAAAATGTACATg                                                                     <  2:4502964/90‑1 (MQ=255)
                       aTTAATATTGGCGCAGTCGTTGGCGAAGGTAGAATGATTGATATGAATGCTACTCTCTGTGGTCGTGCTACAACTGGTAAGAATGTACATg                                                                     >  1:5166574/1‑91 (MQ=255)
                       aTTAATATTGGCGCAGTCGTTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATg                                                                     <  1:4815007/91‑1 (MQ=255)
                        ttAATATTGGCGCAGTCGTTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGt                                                                    <  1:59895/91‑1 (MQ=255)
                        ttAATATTGGCGCAGTCGTTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGt                                                                    >  1:651855/1‑91 (MQ=255)
                        ttAATATTGGCGCAGTCGTTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGt                                                                    <  1:1945124/91‑1 (MQ=255)
                          aaTATTGGCGCAGTCGTTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAg                                                                  <  2:4422996/91‑1 (MQ=255)
                            tatTGGGGCAGTCGTTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAggg                                                                <  2:1094926/91‑1 (MQ=255)
                            tatTGGCGCAGTCGTTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAggg                                                                <  1:329986/91‑1 (MQ=255)
                            taatGGCGCAGGCGTTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAggg                                                                <  2:5435541/88‑1 (MQ=255)
                             atTGGCGCAGTCATTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAgggg                                                               <  2:3655262/91‑1 (MQ=255)
                              ttGGCGCAGTCGTTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGc                                                              <  1:3231607/91‑1 (MQ=255)
                              ttGGCGCAGTCGTTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGc                                                              <  1:1247829/91‑1 (MQ=255)
                              ttGGCGCAGTCGTTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGc                                                              <  1:3404190/91‑1 (MQ=255)
                                ggCGCAGTCGTTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTg                                                            <  2:3825879/91‑1 (MQ=255)
                                   gcAGTCGTTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGgcg                                                         >  1:1596516/1‑91 (MQ=255)
                                   gcAGTCGTTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGgcg                                                         <  1:788423/91‑1 (MQ=255)
                                   gcAGTCGTTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGgcg                                                         <  1:1418142/91‑1 (MQ=255)
                                   gcAGTCGTTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGgcg                                                         <  1:4234945/91‑1 (MQ=255)
                                    cAGTCGTTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGgcgc                                                        <  2:544371/91‑1 (MQ=255)
                                    cAGTCGTTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGgcgc                                                        <  1:3793260/91‑1 (MQ=255)
                                      gTCGTTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAg                                                      >  1:1552948/1‑91 (MQ=255)
                                      gTCGTTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAg                                                      >  1:3045893/1‑91 (MQ=255)
                                      gTCGTTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAg                                                      <  1:2402414/91‑1 (MQ=255)
                                       tCGTTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTCCTACCACTGGAAAAAATGTACATGTAGGGGCTGGCGCAGt                                                     >  2:4114137/1‑91 (MQ=255)
                                        cGTTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTa                                                    <  2:4992819/91‑1 (MQ=255)
                                          ttGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTAtt                                                  <  2:1682761/91‑1 (MQ=255)
                                          ttGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGCCGTACCACTGGAAAAAAAGTACATGTAGGGGCTGGCGCAGTAtt                                                  >  2:1374437/1‑91 (MQ=255)
                                             gCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGc                                               <  1:1219627/91‑1 (MQ=255)
                                                  ggTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACCACTGGAAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGgtg                                          >  2:3873240/1‑91 (MQ=255)
                                                  ggTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGgtg                                          <  2:4453191/91‑1 (MQ=255)
                                                  ggTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGgtg                                          >  1:4036969/1‑91 (MQ=255)
                                                   gTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGgtgt                                         <  1:1103292/91‑1 (MQ=255)
                                                   gTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGgtgt                                         >  1:3903409/1‑91 (MQ=255)
                                                   gTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGgtgt                                         <  1:5110946/91‑1 (MQ=255)
                                                      cAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGAt                                      >  1:3490392/1‑91 (MQ=255)
                                                       aaTGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGAtt                                     <  1:3708137/91‑1 (MQ=255)
                                                        aTGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTg                                    >  1:4558805/1‑91 (MQ=255)
                                                        aTGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTg                                    <  1:2099360/91‑1 (MQ=255)
                                                        aTGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTg                                    <  2:2939089/91‑1 (MQ=255)
                                                        aTGATTGATATGAATGCTACTCTCGGTGGTCGCCGTATCATTAAAAAAAAAATCCATGTAGGCGCTGCCGCAGTATTAGCAGGTGTGAttc                                    >  2:2626236/1‑90 (MQ=255)
                                                         tgattgatATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGa                                   <  2:3189981/91‑1 (MQ=255)
                                                           attgatATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAAc                                 <  2:1885573/91‑1 (MQ=255)
                                                            ttgatATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAAcc                                <  1:3591978/91‑1 (MQ=255)
                                                            ttgatATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAAcc                                >  1:1488167/1‑91 (MQ=255)
                                                            ttgatATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAAcc                                <  1:2932521/91‑1 (MQ=255)
                                                            ttgatATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGCGCCGGCGCAGTATTAGCAGGTTTGATTGAAcc                                >  1:850338/1‑91 (MQ=255)
                                                             tgatATGAATGCTACTCTCGGTGGTCGCCGTATCATTAGAAAAAAAATACATGTAGGGGCCGGCGCAGTATTAGCAGGTGTGATTGAAccc                               >  2:2883796/1‑91 (MQ=255)
                                                               atatGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAAccccc                             <  1:4867620/91‑1 (MQ=255)
                                                                 ttgtaTGCTACTCTCGGTGGTCGTGCTACAACTGGTAAGAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAACCCCCTa                           <  2:2385926/87‑1 (MQ=255)
                                                                   gAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAACCCCCTAGt                         <  1:3657587/91‑1 (MQ=255)
                                                                   gAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAACCCCCTAGt                         <  1:3360920/91‑1 (MQ=255)
                                                                    aaTGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAACCCCCTAGTg                        <  2:3991927/91‑1 (MQ=255)
                                                                     aTGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGTAGTATTAGCAGGTGTGATTGAACCCCCTAGTGc                       <  1:5607382/91‑1 (MQ=255)
                                                                      tGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAACCCCCTAGTGCt                      <  2:1503340/91‑1 (MQ=255)
                                                                       gCTACTCTCGGTGGTCGTCCTACCACTGGAAAAAAAGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAACCCCCTAGTGCtt                     >  2:4886727/1‑91 (MQ=255)
                                                                        cTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAACCCCCTAGTGCTTc                    <  2:2263523/91‑1 (MQ=255)
                                                                          aCTCTCGGTGGTAGTGCTAAAACTGGTAAAGATGTACATGTAGGGGCTGGGGCAGTATTAGCAGGTGTGATTGAACCCCCTAGTGCTTCAc                  <  2:1652141/91‑1 (MQ=255)
                                                                            tctcGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAACCCCCTAGTGCTTCACCg                <  1:1766980/91‑1 (MQ=255)
                                                                            tctcGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAACCCCCTAGTGCTTCACCg                >  1:4754862/1‑91 (MQ=255)
                                                                            tctcGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAACCCCCTAGTGCTTCACCg                >  1:1128737/1‑91 (MQ=255)
                                                                             ctcGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAACCCCCTAGTGCTTCACCgg               <  1:1475855/91‑1 (MQ=255)
                                                                             ctcGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAACCCCCTAGTGCTTCACCgg               >  1:2712010/1‑91 (MQ=255)
                                                                              tcGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAACCCCCTAGTGCTTCACCGGt              >  1:1683534/1‑91 (MQ=255)
                                                                              tcGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAACCCCCTAGTGCTTCACCGGt              <  2:2345233/91‑1 (MQ=255)
                                                                              tcGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAACCCCCTAGTGCTTCACCGGt              <  2:454506/91‑1 (MQ=255)
                                                                              tcGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCGGGGGCAGTATTAGCAGGTGTGATTGAACCCCCTAGTGCTTCACCGGt              <  2:1011069/91‑1 (MQ=255)
                                                                               tggtggtCGTGCTAAAACTGGTATAAATGTATATGTTGGGGCTGGCGCAGTATTAGCATGTGTGATTGACCCTCCTAGTGCTTCACAGGtt             <  2:5528669/90‑1 (MQ=255)
                                                                               cGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAACCCCCTAGTGCTTCACCGGtt             >  1:5268834/1‑91 (MQ=255)
                                                                               cGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAACCCCCTAGTGCTTCACCGGtt             <  2:757178/91‑1 (MQ=255)
                                                                                     tCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAACCCCCTAGTGCTTCACCGGTTATAATc       <  2:3563482/91‑1 (MQ=255)
                                                                                     tCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCACGTGTGATTGAACCCCCTAGTGCTTCACCGGTTATAATc       >  1:4680273/1‑91 (MQ=255)
                                                                                        tGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAACCCCCTAGTGCTTCACCGGTTATAATCGAg    <  2:1571574/91‑1 (MQ=255)
                                                                                          cTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAACCCCCTAGTGCTTCACCGGTTATAATCGAGga  >  2:2838449/1‑91 (MQ=255)
                                                                                          |                                                                                          
CTGTCGTTATGATGGGCGCAACAATTAATATTGGCGCAGTCGTTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAACCCCCTAGTGCTTCACCGGTTATAATCGAGGA  >  USA300TCH1516_ALE/1434975‑1435155

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: