Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 1,221,653 G→A 10.5% intergenic (+68/‑148) rpoZ → / → coaBC DNA‑directed RNA polymerase subunit omega/Coenzyme A biosynthesis bifunctional protein CoaBC

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE1,221,6530GA10.5% 654.3 / 46.4 224intergenic (+68/‑148)rpoZ/coaBCDNA‑directed RNA polymerase subunit omega/Coenzyme A biosynthesis bifunctional protein CoaBC
Reads supporting (aligned to +/- strand):  ref base G (94/106);  new base A (2/22);  total (96/128)
Fisher's exact test for biased strand distribution p-value = 1.54e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 232 total aligned reads displayed.

GTGATTATTATGGTAAAGAATAGTTTCGACCATTAAAAATATGTGAAGCTAGACGTACATCAAATACGTCAAGGCATATCCAAATCAAACGCACCTATGTAGGTGTGTTTTTTTAGTTTATCCAATGATTGATGTTATAATAATACTAAATTTGTATCTATAAAAAAGTAATGAGCATTTG  >  USA300TCH1516_ALE/1221563‑1221743
                                                                                          |                                                                                          
gTGATTATTATGGTAAAGAATAGTGTCGACCATTAAAAATATGTGAAGCTAGACGTACATCAAATACGTCAAGGCATATCCAAATCAAACg                                                                                            <  2:8593/91‑1 (MQ=255)
 tGATTATTATGGTAAAGAATAGTTTCGACCATTAAAAATATGTGAAGCTAGACGTACATCAAATACGTCAAGGCATATCCAAATCAAACGc                                                                                           <  1:1266701/91‑1 (MQ=255)
  gATTATTATGGTAAAGAATAGTTTCGACCATTAAAAATATGTGAAGCTAGACGTACATCAAATACGTCAAGGCATATCCAAATCAAACGCa                                                                                          <  1:1026483/91‑1 (MQ=255)
  gATTATTATGGTAAAGAATAGTTTCGACCATTAAAAATATGTGAAGCTAGACGTACATCAAATACGTCAAGGCATATCCAAATCAAACGCa                                                                                          >  2:2411808/1‑91 (MQ=255)
   attattATGGTAAAGAATAGTTTCGACCATTAAAAATATGTGAAGCTAGACGTACATCAAATACGTCAAGGCATATCCAAATCAAACGCAc                                                                                         <  2:2348244/91‑1 (MQ=255)
    ttattaTGGTAAAGAATAGTTTCGACCATTAAAAATATGTGAAGCTAGACGTACATCAAATACGTCAAGGCATATCCAAATCAAACGCAcc                                                                                        <  1:45332/91‑1 (MQ=255)
      attagGGTAAAGAATAGTTTCGACCATTAAAAATATGTGAAGCTAGACGTACATCAAATACGTCAAGGCATATCCAAATCAAACACACCTa                                                                                      <  2:454153/91‑1 (MQ=255)
       ttatGGTAAAGAATAGTTTCGACCATTAAAAATATGTGAAGCTAGACGTACATCAAATACGTCAAGGCATATCCAAATCAAACGCACCTAt                                                                                     <  1:1024094/91‑1 (MQ=255)
        tatGGTAAAGAATAGTTTCGACCATTAAAAATATGTGAAGCTAGACGTACATCAAATACGTCAAGGCATATCCAAATCAAACGCACCTATg                                                                                    >  2:3326930/1‑91 (MQ=255)
        tatGGTAAAGAATAGTTTCGACCATTAAAAATATGTGAAGCTAGACGTACATCAAATACGTCAAGGCATATCCAAATCAAACACACCTaca                                                                                    <  1:1366214/91‑3 (MQ=255)
        tatGGTAAAGAATAGTTTCGACCATTAAAAACATGTGAAGCTAGACGTACATCAAATACGTCAAGGCATATCCAAATCAAACGCACCTATg                                                                                    >  1:3509017/1‑91 (MQ=255)
         atGGTAAAGAATAGTTTCGACCATTAAAAATATGTGAAGCTAGACGTACATCAAATACGTCAAGGCATATCCAAATCAAACGCACCTATGt                                                                                   >  2:2291441/1‑91 (MQ=255)
         atGGTAAAGAATAGTTTCGACCATTAAAAATATGTGAAGCTAGACGTACATCAAATACGTCAAGGCATATCCAAATCAAACACACCTacat                                                                                   <  2:1154097/91‑4 (MQ=255)
          tGGTAAAGAATAGTTTCGACCATTAAAAATATGTGAAGCTAGACGTACATCAAATACGTCAAGGCATATCCAAATCAAACGCACCTATGTa                                                                                  <  1:669293/91‑1 (MQ=255)
          tGGTAAAGAATAGTTTCGACCATTAAAAATATGTGAAGCTAGACGTACATCAAATACGTCAAGGCATATCCAAATCAAACACACCTacata                                                                                  <  1:192278/91‑5 (MQ=255)
            gTAAAGAATAGTTTCGACCATTAAAAATATGTGAAGCTAGACGTACATCAAATACGTCAAGGCATATCCAAATCAAACGCACCTATGTAgg                                                                                <  2:2069661/91‑1 (MQ=255)
             tAAAGAATAGTTTCGACCATTAAAAATATGTGAAGCTAGACGTACATCAAATACGTCAAGGCATATCCAAATCAAACGCACCTATGTAGgt                                                                               <  1:260937/91‑1 (MQ=255)
             tAAAGAATAGTTTCGACCATTAAAAATATGTGAAGCTAGACGTACATCAAATACGTCAAGGCATATCCAAATCAAACACACCTACATAGgt                                                                               <  2:3154093/91‑1 (MQ=255)
               aaGAATAGTTTCGACCATTAAAAATATGTGAAGCTAGACGTACATCAAATACGTCAAGGCATATCCAAATCAAACGCACCTATGTAGgtgt                                                                             >  2:2986055/1‑91 (MQ=255)
               aaGAATAGTTTCGACCATTAAAAATATGTGAAGCTAGACGTACATCAAATACGTCAAGGCATATCCAAATCAAACACACCTACATAGgtgc                                                                             <  1:991931/91‑2 (MQ=255)
                aGAATAGTTTCGACCATTAAAAATATGTGAAGCTAGACGTACATCAAATACGTCAAGGCATATCCAAATCAAACGCACCTATGTAGgtgtg                                                                            <  2:1029066/91‑1 (MQ=255)
                 gAATAGTTTCGACCATTAAAAATATGTGAAGCTAGACGTACATCAAATACGTCAAGGCATATCCAAATCAAACGCACCTATGTAGgtgtgt                                                                           <  1:3598009/91‑1 (MQ=255)
                  aaTAGTTTCGACCATTAAAAATATGTGAAGCTAGACGTACATCAAATACGTCAAGGCATATCCAAATCAAACGCACCTATGTAGGTGTGtt                                                                          <  1:2041659/91‑1 (MQ=255)
                   aTAGTTTCGACCATTAAAAATATGTGAAGCTAGACGTACATCAAATACGTCAAGGCATATCCAAATCAAACGCACCTATGTAGGTGTGttt                                                                         >  2:3573719/1‑91 (MQ=255)
                    tAGTTTCGACCATTAAAAATATGTGAAGCTAGACGTACATCAAATACGTCAAGGCATATCCAAATCAAACGCACCTATGTAGGTGTGtttt                                                                        <  1:1600933/91‑1 (MQ=255)
                      gTTTCGACCATTAAAAATATGTGAAGCTAGACGTACATCAAATACGTCAAGGCATATCCAAATCAAACGCACCTATGTAGGTGTGtttttt                                                                      <  1:867888/91‑1 (MQ=255)
                        ttCGACCATTAAAAATATGTGAAGCTAGACGTACATCAAATACGTCAAGGCATATCCAAATCAAACGCACCTATGTAGGTGTGTTTTTTTa                                                                    >  2:2814700/1‑91 (MQ=255)
                         tCGACCATTAAAAATATGTGAAGCTAGACGTACATCAAATACGTCAAGGCATATCCAAATCAAACGCACCTATGTAGGTGTGTTTTTTTAg                                                                   <  1:1164342/91‑1 (MQ=255)
                          cGACCATTAAAAATATGTGAAGCTAGACGTACATCAAATACGTCAAGGCATATCCAAATCAAACGCACCTATGTAGGTGTGTTTTTTTAGt                                                                  >  2:1802695/1‑91 (MQ=255)
                           gACCATTAAAAATATGTGAAGCTAGACGTACATCAAATACGTCAAGGCATATCCAAATCAAACGCACCTATGTAGGTGTGTTTTTTTAGtt                                                                 >  2:1430759/1‑91 (MQ=255)
                              cATTAAAAATATGTGAAGCTAGACGTACATCAAATACGTCAAGGCATATCCAAATCAAACGCACCTATGTAGGTGTGTTTTTTTAGTTTAt                                                              <  2:2141694/91‑1 (MQ=255)
                              cATTAAAAATATGTGAAGCTAGACGTACATCAAATACGTCAAGGCATATCCAAATCAAACGCACCTATGTAGGTGTGTTTTTTTAGTTTAt                                                              <  2:2865242/91‑1 (MQ=255)
                                ttAAAAATATGTGAAGCTAGACGTACATCAAATACGTCAAGGCATATCCACATCAAACGCACCTATGCAGGCGTGTTTTTTTAGTTTATcc                                                            >  2:2185277/1‑91 (MQ=255)
                                ttAAAAATATGTGAAGCTAGACGTACATCAAATACGTCAAGGCATATCCAAATCAAACGCACCTATGTAGGTGTGTTTTTTTAGTTTATcc                                                            >  1:2140931/1‑91 (MQ=255)
                                ttAAAAATATGTGAAGCTAGACGTACATCAAATACGTCAAGGCATATCCAAATCAAACGCACCTATGTAGGTGTGTTTTTTTAGTTTATcc                                                            >  1:186685/1‑91 (MQ=255)
                                 tAAAAATATGTGAAGCTAGACGTACATCAAATACGTCAAGGCATATCCAAATCAAACGCACCTATGTAGGTGTGTTTTTTTAGTTTATCCa                                                           >  1:180912/1‑91 (MQ=255)
                                 tAAAAATATGTGAAGCTAGACGTACATCAAATACGTCAAGGCATATCCAAATCAAACGCACCTATGTAGGTGTGTTTTTTTAGTTTATCCa                                                           <  1:1157909/91‑1 (MQ=255)
                                 tAAAAATATGTGAAGCTAGACGTACATCAAATACGTCAAGGCATATCCAAATCAAACGCACCTATGTAGGTGTGTTTTTTTAGTTTATCCa                                                           <  2:1493270/91‑1 (MQ=255)
                                   aaaaTATGTGAAGCTAGACGTACATCAAATACGTCAAGGCATATCCAAATCAAACGCACCTATGTAGGTGTGTTTTTTTAGTTTATCCAAt                                                         >  1:595399/1‑91 (MQ=255)
                                    aaaTATGTGAAGCTAGACGTACATCAAATACGTCAAGGCATATCCAAATCAAACGCACCTATGTAGGTGTGTTTTTTTAGTTTATTCAAtg                                                        >  2:110784/1‑91 (MQ=255)
                                     aaTATGTGAAGCTAGACGTACATCAAATACGTCAAGGCATATCCAAATCAAACGCACCTATGTAGGTGTGTTTTTTTAGTTTATCCAAtga                                                       >  2:877278/1‑91 (MQ=255)
                                     aaTATGTGAAGCTAGACGTACATCAAATACGTCAAGGCATATCCAAATCAAACGCACCTATGTAGGTGTGTTTTTTTAGTTTATCCAAtga                                                       <  2:2209528/91‑1 (MQ=255)
                                     aaTATGTGAAGCTAGACGTACATCAAATACGTCAAGGCATATCCAAATCAAACGCACCTATGTAGGTGTGTTTTTTTAGTTTATCCAAtga                                                       <  1:3275354/91‑1 (MQ=255)
                                       tatGTGAAGCTAGACGTACATCAAATACGTCAAGGCATATCCAAATCAAACGCACCTATGTAGGTGTGCTTTTTTAGTTTATCCAAtgatt                                                     >  2:3184380/1‑91 (MQ=255)
                                         tgtgAAGCTAGACGTACATCAAATACGTCAAGGCATATCCAAATCAAACGCACCTATGTAGGTGTGTTTTTTTAGTTTATCCAAtgattga                                                   >  1:704105/1‑91 (MQ=255)
                                         tgtgAAGCTAGACGTACATCAAATACGTCAAGGCATATCCAAATCAAACGCACCTATGTAGGTGTGTTTTTTTAGTTTATCCAAtgattga                                                   >  1:1755112/1‑91 (MQ=255)
                                         tgtgAAGCTAGACGTACATCAAATACGTCAAGGCATACCCAAATCAAACGCCCCTATGTAGGTGTATTTTTTTAGTTTATCCAAtcattga                                                   >  1:3679492/1‑91 (MQ=255)
                                          gtgAAGCTAGACGTACATCAAATACGTCAAGGCATATCCAAATCAAACGCACCTATGTAGGTGTGTTTTTTTAGTTTATCCAAtgattgat                                                  >  1:1107416/1‑91 (MQ=255)
                                             aaGCTAGACGTACATCTAATACGTCAAGGCATATCCAAATCAAACGCACCTATGTAGGTGTGTTTTTTTAGTTTATCCAATGATTGATGtt                                               <  2:3035465/91‑1 (MQ=255)
                                             aaGCTAGACGTACATCAAATACGTCAAGGCATATCCAAATCAAACGCACCTATGTAGGTGTGTTTTTTTAGTTTATCCAATGATTGATGtt                                               <  1:2002768/91‑1 (MQ=255)
                                               gCTAGACGTACATCAAATACGTCAAGGCATATCCAAATCAAACGCACCTATGTAGGTGTGTTTTTTTAGTTTATCCAATGATTGATGTtat                                             <  1:2249165/91‑1 (MQ=255)
                                                cTAGACGTACATCAAATACGTCAAGGCATATCCAAATCAAACGCACCTATGTAGGTGTGTTTTTTTAGTTTATCCAATGATTGATGTtata                                            <  2:3307094/91‑1 (MQ=255)
                                                   gACGTACATCAAATACGTCAAGGCATATCCAAATCAAACGCACCTATGTAGGTGTGTTTTTTTAGTTTATCCAATGATTGATGTTataata                                         <  1:881628/91‑1 (MQ=255)
                                                     cGTACATCAAATACGTCAAGGCATATCCAAATCAAAGGCACCTATGTAGGTGTGTTTTTTTAGTTTATCCAATGATTGATGTTataataat                                       <  2:1985264/91‑1 (MQ=255)
                                                     cGTACATCAAATACGTCAAGGCATATCCAAATCAAACGCACCTATGTAGGTGTGTTTTTTTAGTTTATCCAATGATTGATGTTataataat                                       <  1:1754194/91‑1 (MQ=255)
                                                       tACATCAAATACGTCAAGGCATATCCAAATCAAACGCACCTATGTAGGTGTGTTTTTTTAGTTTATCCAATGATTGATGTTATAATAATAc                                     <  2:501087/91‑1 (MQ=255)
                                                        aCATCAAATACGTCAAGGCATATCCAAATCAAACGCACCTATGTAGGTGTGTTTTTTTAGTTTATCCAATGATTGATGTTATAATAATACt                                    <  2:1010398/91‑1 (MQ=255)
                                                         cATCAAATACGTCAAGGCATATCCAAATCAAACGCACCTATGTAGGTGTGTTTTTTTAGTTTATCCAATGATTGATGTTATAATAATACTa                                   >  1:3814206/1‑91 (MQ=255)
                                                         cATCAAATACGTCAAGGCATATCCAAATCAAACGCACCTATGTAGGTGTGTTTTTTTAGTTTATCCAATGATTGATGTTATAATAATACTa                                   <  1:3392889/91‑1 (MQ=255)
                                                           gcAAATACGTCAAGGCATATCCAAATCAAACGCACCTATGTAGGTGCGTTATTTTAGTTTATCCAATGCTTGATGTTATAATAATACTaaa                                 >  2:732501/2‑91 (MQ=255)
                                                             aaaTACGTCAAGGCATATCCAAATCAAACGCACCTATGTAGGTGTGTTTTTTTAGTTTATCCAATGATTGATGTTATAATAATACTAAAtt                               <  2:1532526/91‑1 (MQ=255)
                                                             aaaTACGTCAAGGCATATCCAAATCAAACGCACCTATGTAGGTGTGTTTTTTTAGTTTATCCAATGATTGATGTTATAATAATACTAAAtt                               <  1:1124071/91‑1 (MQ=255)
                                                              aaTACGTCAAGGCATATCCAAATCAAACGCACCTATGTAGGTGTGTTTTTTTAGTTTATCCAATGATTGATGTTATAATAATACTAAAttt                              <  2:1682141/91‑1 (MQ=255)
                                                              aaTACGTCAAGGCATATCCAAATCAAACGCACCTATGTAGGTGTGTTTTTTTAGTTTATCCAATGATTGATGTTATAATAATACTAAAttt                              <  2:3822389/91‑1 (MQ=255)
                                                               aTACGTCAAGGCATATCCAAATCAAACGCACCTATGTAGGTGTGTTTTTTTAGTTTATCCAATGATTGATGTTATAATAATACTAAATTTg                             <  2:1714134/91‑1 (MQ=255)
                                                                tACGTCAAGGCATATCCAAATCAAACGCACCTATGTAGGTGTGTTTTTTTAGTTTATCCAATGATTGATGTTATAATAATACTAAATTTGt                            >  2:2945960/1‑91 (MQ=255)
                                                                 aCGTCAAGGCATATCCAAATCAAACGCACCTATGTAGGTGTGTTTTTTTAGTTTATCCAATGATTGATGTTATAATAATACTAAATTTGTa                           >  1:1144051/1‑91 (MQ=255)
                                                                 aCGTCAAGGCATATCCAAATCAAACGCACCTATGTAGGTGTGTTTTTTTAGTTTATCCAATGATTGATGTTATAATAATACTAAATTTGTa                           >  1:2181849/1‑91 (MQ=255)
                                                                  cGTCAAGGCATATCCAAATCAAACGCACCTATGTAGGTGTGTTTTTTTAGTTTATCCAATGATTGATGTTATAATAATACTAAATTTGTAt                          <  1:3090554/91‑1 (MQ=255)
                                                                   gTCAAGGCATATCCAAATCAAACGCACCTATGTAGGTGTGTTTTTTTAGTTTATCCAATGATTGATGTTATAATAATACTAAATTTGTATc                         >  1:1416952/1‑91 (MQ=255)
                                                                   gTCAAGGCATATCCAAATCAAACGCACCTATGTAGGTGTGTTTTTTTAGTTTATCCAATGATTGATGTTATAATAATACTAAATTTGTATc                         >  2:2256365/1‑91 (MQ=255)
                                                                   gTCAAGGCATATCCAAATCAAACGCACCTATGTAGGTGTGTTTTTTTAGTTTATCCAATGATTGATGTTATAATAATACTAAATTTGTATc                         <  1:1511459/91‑1 (MQ=255)
                                                                   gTCAAGGCATATCCAAATCAAACGCACCTATGTAGGTGTGTTTTTTTAGTTTATCCAATGATTGATGTTATAATAATACTAAATTTGTATc                         <  2:2554701/91‑1 (MQ=255)
                                                                   gTCAAGGCATATCCAAATCAAACGCACCTATGTAGGTGTGTTTTTTTAGTTTATCCAATGATTGATGTTATAATAATACTAAATTTGTATc                         <  1:3022370/91‑1 (MQ=255)
                                                                     cAAGGCATATCCAAATCAAACGCACCTATGTAGGTGTGTTTTTTTAGTTTATCCAATGATTGATGTTATAATAATACTAAATTTGTATCta                       >  2:2645753/1‑91 (MQ=255)
                                                                     cAAGGCATATCCAAATCAAACGCACCTATGTAGGTGTGTTTTTTTAGTTTATCCAATGATTGATGTTATAATAATACTAAATTTGTATCta                       >  2:685065/1‑91 (MQ=255)
                                                                      aaGGCATATCCAAATCAAACGCACCTATGTAGGTGTGTTTTTTTAGTTTATCCAATGATTGATGTTATAATAATACTAAATTTGTATCtat                      <  1:1895431/91‑1 (MQ=255)
                                                                      aaGGCATATCCAAATCAAACGCACCTATGTAGGTGTGTTTTTTTAGTTTATCCAATGATTGATGTTATAATAATACTAAATTTGTATCtat                      >  2:728826/1‑91 (MQ=255)
                                                                       aGGCATATCCAAATCAGACGCACCTATGTAGGTGTGTTTTTTTAGTTTATCCAATGATTGATGTTATAATAATACTAAATTTGTATCtata                     <  2:3777012/91‑1 (MQ=255)
                                                                        ggCATATCCAAATCAAACGCACCTATGTAGGTGTGTTTTTTTAGTTTATCCAATGATTGATGTTATAATAATACTAAATTTGTATCTATaa                    >  2:2169686/1‑91 (MQ=255)
                                                                         gCATATCCAAATCAAACGCACCTATGTAGGTGTGTTTTTTTAGTTTATCCAATGATTGATGTTATAATAATACTAAATTTGTATCTATaaa                   <  1:3240845/91‑1 (MQ=255)
                                                                         gCATATCCAAATCAAACGCACCTATGTAGGTGTGTTTTTTTAGTTTATCCAATGATTGATGTTATAATAATACTAAATTTGTATCTATaaa                   <  2:1363345/91‑1 (MQ=255)
                                                                          cATATCCAAATCAAACGCACCTATGTAGGTGTGTTTTTTTAGTTTATCCAATGATTGATGTTATAATAATACTAAATTTGTATCTATaaaa                  <  1:404891/91‑1 (MQ=255)
                                                                          cATATCCAAATCAAACGCACCTATGTAGGTGTGTTTTTTTAGTTTATCCAATGATTGATGTTATAATAATACTAAATTTGTATCTATaaaa                  >  1:308048/1‑91 (MQ=255)
                                                                           atatCCAAATCAAACGCACCTATGTAGGTGTGTTTTTTTAGTTTATCCAATGATTGATGTTATAATAATACTAAATTTGTATCTATaaaaa                 <  1:1542101/91‑1 (MQ=255)
                                                                           atatCCAAATCAAACGCACCTATGTAGGTGTGTTTTTTTAGTTTATCCAATGATTGATGTTATAATAATACTAAATTTGTATCTATaaaaa                 <  1:2823199/91‑1 (MQ=255)
                                                                            tatCCAAATCAAACGCACCTATGTAGGTGTGTTTTTTTAGTTTATCCAATGATTGATGTTATAATAATACTAAATTTGTATCTATaataaa                >  1:3188624/1‑91 (MQ=255)
                                                                             atCCAAATCAAACGCACCTATGTAGGTGTGTTTTTTTAGTTTATCCAATGATTGATGTTATAATAATACTAAATTTGTATCTATAAAAAAg               >  2:1562397/1‑91 (MQ=255)
                                                                              tCCAAATCAAACGCACCTATGTAGGTGTGTTTTTTTAGTTTATCCAATGATTGATGTTATAATAATACTAAATTTGTATCTATAAAAAAGt              >  2:3959545/1‑91 (MQ=255)
                                                                                  aaaaaaaCACACCTACATAGGTGCGTTTTTTTAGTTTATCCAATGATTGATGTTATAATAATACTAAATTTGTATCTATAAAAAAGTAATg          >  1:2916120/5‑91 (MQ=37)
                                                                                  aaTCAAACGCACCTATGTAGGTTTGTTTTTTTAGTTTATCCAATGATTGATGTTATAATAATACTAAATTTGTATCTATAAAAAAGTAATg          <  2:2603898/91‑1 (MQ=255)
                                                                                  aaTCAAACGCACCTATGTAGGTGTGTTTTTTTAGTTTATCCAATGATTGATGTTATAATAATACTAAATTTGTATCTATAAAAAAGTAATg          <  2:46691/91‑1 (MQ=255)
                                                                                  aaTCAAACGCACCTATGTAGGTGTGTTTTTTTAGTTTATCCAATGATTGATGTTATAATAATACTAAATTTGTATCTATAAAAAAGTAATg          <  1:3512303/91‑1 (MQ=255)
                                                                                    tCAAACGCACCTATGTAGGTGTGTTTTTTTAGTTTATCCAATGATTGATGTTATAATAATACTAAATTTGTATCTATAAAAAAGTAATGAg        >  2:780416/1‑91 (MQ=255)
                                                                                      aaaCGCACCTATGTAGGTGTGTTTTTTTAGTTTATCCAATGATTGATGTTATAATAATACTAAATTTGTATCTATAAAAAAGTAATGAGCa      >  2:3065246/1‑91 (MQ=255)
                                                                                      aaaCGCACCTATGTAGGTGTGTTTTTTTAGTTTATCCAATGATTGATGTTATAATAATACTAAATTTGTATCTATAAAAAAGTAATGAGCa      <  2:1595411/91‑1 (MQ=255)
                                                                                       aaCGCACCTATGTAGGTGTGTTTTTTTAGTTTATCCAATGATTGATGTTATAATAATACTAAATATGTATCTATAAAAAAGTAATGAGCAt     >  1:386474/1‑91 (MQ=255)
                                                                                        aCGCACCTATGTAGGTGTGTTTTTTTAGTTTATCCAATGATTGATGTTATAATAATACTAAATTTGTATCTATAAAAAAGTAATGAGCAtt    >  1:125663/1‑91 (MQ=255)
                                                                                        aCGCACCTATGTAGGTGTGTTTTTTTAGTTTATCCAATGATTGATGTTATAATAATACTAAATTTGTATCTATAAAAAAGTAATGAGCAtt    >  2:1955499/1‑91 (MQ=255)
                                                                                         cGCACCTATGTAGGTGTGTTTTTTTAGTTTATCCAATGATTGATGTTATAATAATACTAAATTTGTATCTATAAAAAAGTAATGAGCAttt   >  1:264958/1‑91 (MQ=255)
                                                                                          |                                                                                          
GTGATTATTATGGTAAAGAATAGTTTCGACCATTAAAAATATGTGAAGCTAGACGTACATCAAATACGTCAAGGCATATCCAAATCAAACGCACCTATGTAGGTGTGTTTTTTTAGTTTATCCAATGATTGATGTTATAATAATACTAAATTTGTATCTATAAAAAAGTAATGAGCATTTG  >  USA300TCH1516_ALE/1221563‑1221743

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