Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 1,435,065 C→G 19.9% 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,0650CG19.9% 418.5 / 33.4 155A141G (GCT→GGT) dapH2,3,4,5‑tetrahydropyridine‑2,6‑dicarboxylate N‑acetyltransferase
Reads supporting (aligned to +/- strand):  ref base C (42/80);  new base G (30/2);  total (73/82)
Fisher's exact test for biased strand distribution p-value = 5.24e-10
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 160 total aligned reads displayed.

CTGTCGTTATGATGGGCGCAACAATTAATATTGGCGCAGTCGTTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAACCCCCTAGTGCTTCACCGGTTATAATCGAGGA  >  USA300TCH1516_ALE/1434975‑1435155
                                                                                          |                                                                                          
cTGTCGTTATGATGGGCGCAACAATTAATATTGGCGCAGTCGTTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGc                                                                                            <  1:1273525/91‑1 (MQ=255)
cTGTCGTTATGATGGGCGCAACAATTAATATTGGCGCAGTCGTTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGc                                                                                            <  1:1273684/91‑1 (MQ=255)
  gTCGTTATGATGGGCGCAACAATTAATATTGGCGCAGTCGTTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTa                                                                                          >  1:540275/1‑91 (MQ=255)
   tCGTTATGATGGGCGCAACAATTAATATTGGCGCAGTCGTTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGGTAc                                                                                         >  2:1184525/1‑91 (MQ=255)
     gTTATGATGGGCGCAACAATTAATATTGGCGCAGTCGTTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACaa                                                                                       <  2:2538999/91‑1 (MQ=255)
       tATGATGGGCGCAACAATTAATATTGGCGCAGTCGTTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGCGCTACCACt                                                                                     >  2:579864/1‑91 (MQ=255)
       tATGATGGGCGCAACAATTAATATTGGCGCAGTCGTTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGCCCTAccatt                                                                                     >  2:1850929/1‑87 (MQ=255)
           atgGGCGCAACAATTAATATTGGCGCAGTCGTTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTa                                                                                 <  1:2745091/91‑1 (MQ=255)
            tgGGCGCAACAATTAATATTGGCGCAGTCGTTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTaa                                                                                >  1:2354325/1‑91 (MQ=255)
            tgGGCGCAACAATTAATATTGGCGCAGTCGTTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTaa                                                                                >  1:3071949/1‑91 (MQ=255)
             gggCGCAACAATTAATATTGGCGCAGTCGTTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTaaa                                                                               <  1:624390/91‑1 (MQ=255)
             gggCGCAACAATTAATATTGGCGCAGTCGTTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTaaa                                                                               <  1:1058504/91‑1 (MQ=255)
              ggCGCAACAATTAATATTGGCGCGGTCGTTGGCGAGGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTaaaa                                                                              <  2:2137890/91‑1 (MQ=255)
               gcgcAACAATTAATATTGGCGCAGTCGTTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTaaaaa                                                                             >  1:966433/1‑91 (MQ=255)
               gcgcAACAATTAATATTGGCGCAGTCGTTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAAATGGTaaaaa                                                                             <  2:507894/91‑1 (MQ=255)
               gcgcAACAATTAATATTGGCGCAGTCGTTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTCCTACAATTGGaaaaaa                                                                             >  2:2333268/1‑91 (MQ=255)
                cgcAACAATTAATATTGGCGCAGTCGTTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGCGGTACCACTGGaaaaaaa                                                                            >  2:640502/1‑90 (MQ=255)
                 gcAACAATTAATATTGGCGCAGTCGTTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATg                                                                           <  1:171892/91‑1 (MQ=255)
                  caacaaTTAATATTGGCGCAGTCGTTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGt                                                                          <  1:2542643/91‑1 (MQ=255)
                    acaaTTAATATTGGCGCAGTCGTTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGCGGTACCATTAGAAAAAAAGTAc                                                                        >  2:1781026/1‑91 (MQ=255)
                     caaTTAATATTGGCGCAGTCGTTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACa                                                                       >  1:3409460/1‑91 (MQ=255)
                       aTTAATATTGGCGCGGTAGTTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATg                                                                     <  2:3823836/91‑1 (MQ=255)
                       aTTAATATTGGCGCAGTCGTTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATg                                                                     >  2:3128213/1‑91 (MQ=255)
                        ttAATATTGGCGCAGTCGTTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGt                                                                    <  2:244787/91‑1 (MQ=255)
                        ttAATATTGGCGCAGTCGTTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGCGGTACCACTAGAAAAAAAGTACATGt                                                                    >  2:2422796/1‑91 (MQ=255)
                        ttAATATTGGCGCAGTCGTTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGCCGTACCACTAAAAAAAAAATACATGt                                                                    >  2:2112308/1‑91 (MQ=255)
                          aaTATTGGCGCAGTCGTTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTCGTACCATTGGAAAAAATGTACATGTAg                                                                  >  2:1286132/1‑91 (MQ=255)
                           atatTGGCGCAGTCGTTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAgg                                                                 >  1:3329086/1‑91 (MQ=255)
                            tatTGGCGCAGTCGTTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAggg                                                                <  2:3015757/91‑1 (MQ=255)
                              ttGGCGCAGTCGTTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGc                                                              <  2:3481765/91‑1 (MQ=255)
                              ttGGCGCAGTCGTTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGc                                                              <  2:287050/91‑1 (MQ=255)
                                ggCGCAGTCGTTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTg                                                            <  1:3141622/91‑1 (MQ=255)
                                  cgcAGTCGTTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGAGCTACAACTGGTAAAAATGTACATGTAGGGGCTGgc                                                          <  2:2391037/91‑1 (MQ=255)
                                  cgcAGGCGTTGGCGAAGGTACAATGTTTGATATGAATGCGACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGgc                                                          <  1:3679075/91‑1 (MQ=255)
                                    taGTCGGTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGgcgc                                                        <  2:2778364/90‑1 (MQ=255)
                                    cAGTCGTTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGgcgc                                                        <  1:1067082/91‑1 (MQ=255)
                                      gTCGTTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAg                                                      <  2:3108153/91‑1 (MQ=255)
                                      gTCGTTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAg                                                      <  1:2402871/91‑1 (MQ=255)
                                      gTCGTTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAg                                                      <  1:1815709/91‑1 (MQ=255)
                                       tcattGGCGAAGGTACAATGATTGATGTGAATGCTGCTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTAGATGTAGGGGCTGGCGCAGt                                                     <  2:717659/88‑1 (MQ=255)
                                        cGTTGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTa                                                    >  1:855170/1‑91 (MQ=255)
                                         gTTGGCGAAGGTAAAATGATTGATATGAGTGCTAATCTAGGTGGTCGTGATACAACTGGTAACACTGTAAACGTAGGGGCTGGCGCAGTAt                                                   >  1:3144057/1‑91 (MQ=255)
                                          ttGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGTTGGCGCAGTAtt                                                  <  1:3036965/91‑1 (MQ=255)
                                          ttGGCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTCGTATCATTGAAAAAAATGTACATGTAGGGGCTGGCGCAGTAtt                                                  >  2:919286/1‑91 (MQ=255)
                                           tGGCGAGGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGGTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTa                                                 <  2:3057552/91‑1 (MQ=255)
                                             gCGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGCCGTATCATTAAAAAAAAAATACATGTAGGGGCTGGCGCAGTATTAGc                                               >  2:738258/1‑91 (MQ=255)
                                              cGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCa                                              >  1:2952451/1‑91 (MQ=255)
                                              cGAAGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGCCGTACCACTAGAAAAAAAGTACATGCAGGGGCTGGCGCAGTATTAGCa                                              >  2:3183299/1‑91 (MQ=255)
                                                aaGGTACAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAgg                                            >  1:3479949/1‑91 (MQ=255)
                                                     aCAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGa                                       <  2:1860981/91‑1 (MQ=255)
                                                     aCAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGa                                       <  1:2058092/91‑1 (MQ=255)
                                                      cAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGAt                                      >  1:2984180/1‑91 (MQ=255)
                                                      cAATGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGAt                                      <  1:2647807/91‑1 (MQ=255)
                                                      cAATGATTGATATGAATGCTACTCTCGGTGGTCGCCGTATCATTAAAAAAAAAATACCTGTAGGGGCTGGCGCAGTATTAGCAGGAGTGAt                                      >  2:1559672/1‑91 (MQ=255)
                                                        aTGATTGATATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTg                                    <  1:2456131/91‑1 (MQ=255)
                                                         tgattgatATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGa                                   >  1:1469080/1‑91 (MQ=255)
                                                         tgattgatATGAATGCTACTCTCGGTGGTCGCCCTACCACTGGTAAAAATGTACATGCAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGa                                   >  2:3050395/1‑91 (MQ=255)
                                                            ttgatATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAAcc                                <  1:2481432/91‑1 (MQ=255)
                                                            ttgatATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAAcc                                >  1:2826438/1‑91 (MQ=255)
                                                            ttgatATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAAcc                                >  1:2869934/1‑91 (MQ=255)
                                                            ttgatATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAAcc                                >  1:2405524/1‑91 (MQ=255)
                                                            ttgatATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAAcc                                >  1:492968/1‑91 (MQ=255)
                                                            ttgatATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAAcc                                >  1:3912088/1‑91 (MQ=255)
                                                            ttgatATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAAcc                                >  1:493049/1‑91 (MQ=255)
                                                            ttgatATGAATGCTACTCTCGGTGGTCGCCGTATCATTAAAAAAAAAATACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAAcc                                >  2:1470045/1‑91 (MQ=255)
                                                            ttgatATGAATGCTACTCTCGGTGGTCGCCGTATCACTAAAAAAAAAGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAAcc                                >  2:1466110/1‑91 (MQ=255)
                                                              gatATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAAcccc                              <  1:2483709/91‑1 (MQ=255)
                                                              gatATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAAcccc                              <  1:3553712/91‑1 (MQ=255)
                                                              gatATGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAAcccc                              <  1:33914/91‑1 (MQ=255)
                                                              gatATGAATGCTACTCTCGGTGGTCGCCGTATCATTAAAAAAAAAATACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAAcccc                              >  2:893146/1‑91 (MQ=255)
                                                              gatATGAATGCTACTCTCGGTGGTCGCCGTATCATTAAAAAAAAAATACATGTAGGGGCTGGCGCAGTATTAGCAGGGATGACTGAAcccc                              >  2:893034/1‑91 (MQ=255)
                                                              gatATGAATGCTACTCTCGGTGGTCGCCGTACCATTAAAAAAAAAATACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAAcccc                              >  2:426362/1‑91 (MQ=255)
                                                                tatGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAACCCCCt                            >  1:1785473/1‑91 (MQ=255)
                                                                tatGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAACCCCCt                            >  1:66129/1‑91 (MQ=255)
                                                                tatGAATGCTACTCTCGGTGGTCGCCGTATCATTAAAAAAAAAATACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAAccccca                            >  2:1934781/1‑90 (MQ=255)
                                                                tatGAATGCTACTCTCGGTGGTCGCCGTATCATTAAAAAAAAAATACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAACCCCCt                            >  2:2462157/1‑91 (MQ=255)
                                                                tatGAATGCTACTCTCGGTGGTCGCCGTATCATTAAAAAAAAAATACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAACCCCCt                            >  2:3606647/1‑91 (MQ=255)
                                                                 atGAATGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGACCCCCCTa                           >  2:3847968/1‑91 (MQ=255)
                                                                    aaTGCTACTGTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAACCCCCTAGTg                        <  1:3785607/91‑1 (MQ=255)
                                                                    aaTGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAACCCCCTAGTg                        <  1:174811/91‑1 (MQ=255)
                                                                    aaTGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAACCCCCTAGTg                        <  1:75272/91‑1 (MQ=255)
                                                                     aTGCTACTCTCGGTGGTCGTGGTACAATTAATAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAACCCCCTAGTGc                       >  2:1723173/1‑91 (MQ=255)
                                                                     aTGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAACCCCCTAGTGc                       <  1:2595357/91‑1 (MQ=255)
                                                                     aTGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACAGGTGGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAACCCCCTAGTGc                       <  2:511539/91‑1 (MQ=255)
                                                                      tGCTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAACCCCCTAGTGCt                      <  1:2294789/91‑1 (MQ=255)
                                                                        cTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAACCCCCTAGTGCTTc                    >  1:3265069/1‑91 (MQ=255)
                                                                        cTACTCTCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAACCCCCTAGTGCTTc                    >  1:2541919/1‑91 (MQ=255)
                                                                        cTACTCTCGGTGGTCGTCGTACCATTGAAAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAACCCCCTAGTGCTTc                    >  2:1990524/1‑91 (MQ=255)
                                                                           ctctCGGTGGTCGTGCTACAACTGGTAAAAGTGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAACCCCCTAGTGCTTCAcc                 <  2:575104/91‑1 (MQ=255)
                                                                           ctctCGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAACCCCCTAGTGCTTCAcc                 <  1:3140105/91‑1 (MQ=255)
                                                                            tctcGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAACCCCCTAGTGCTTCACCg                >  1:2927771/1‑91 (MQ=255)
                                                                             ctcGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAACCCCCTAGTGCTTCACCgg               >  1:2893317/1‑91 (MQ=255)
                                                                              tcGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAACCCCCTAGTGCTTCACCGGt              <  2:2070009/91‑1 (MQ=255)
                                                                              tcGGTGGTCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAACCCCCTAGTGCTTCACCGGt              <  1:3808245/91‑1 (MQ=255)
                                                                                  tggtCGTCGTACCATTGGTAAAAATATACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAACCCCCTAGTGCTTCACCGGTtata          >  2:3743472/1‑91 (MQ=255)
                                                                                   ggtCGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAACCCCCTAGTGCTTCACCGGTTATaa         <  1:2718319/91‑1 (MQ=255)
                                                                                     tCGTGGTAAAAATGGAAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAACCCCCTAGTGCTTCACCGGTTATAATc       <  2:1302906/91‑1 (MQ=255)
                                                                                     tCGTGCTACAACTGGTAAGAATGTACAGGGAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAACCCCCTAGTGCTTCACCGGTTATAATc       <  2:1302819/91‑1 (MQ=255)
                                                                                      cGTGCTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAACCCCCTAGTGCTTCACCGGTTATAATCg      <  2:2730650/91‑1 (MQ=255)
                                                                                          cTACAACTGGTAAAAATGTACATGTAGGGGCTGGCGCAGTATTAGCAGGTGTGATTGAACCCCCTAGTGCTTCACCGGTTATAATCGAGga  <  1:2763778/91‑1 (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: