Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 2,172,081 G→C 11.6% G44R (GGT→CGT)  leuB → 3‑isopropylmalate dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE2,172,0810GC11.6% 561.3 / 48.4 172G44R (GGT→CGT) leuB3‑isopropylmalate dehydrogenase
Reads supporting (aligned to +/- strand):  ref base G (78/74);  new base C (12/8);  total (90/82)
Fisher's exact test for biased strand distribution p-value = 4.87e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.62e-01

Only 100 of 196 total aligned reads displayed.

GGTCCAGAAATTTTGAACGGATCTCTATCATTGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGAGCACCACGAATTTGGTGGTGCCTCTATTGATACATTCGGCGAGCCTTTAACTGAGAAAACCTTAAATGCGTGTAAAAGAGCAGATGCTATTTTACTGGGTGCAA  >  USA300TCH1516_ALE/2171991‑2172171
                                                                                           |                                                                                           
ggTCCAGAAATTTTGAACGGATCTCTATCATTGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGAGCACCACGAATTtg                                                                                             <  2:964900/91‑1 (MQ=255)
 gTCCAGAAATTTTGAACGGATCTCTATCATTGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGAGCACCACGAATTtgg                                                                                            >  1:880680/1‑91 (MQ=255)
 gTCCAGAAATTTTGAACGGATCTCTATCATTGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGAGCACCACGAATTtgg                                                                                            <  1:2439174/91‑1 (MQ=255)
 gTCCAGAAATTTTGAACGGATCTCTATCATTGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGAGCACCACGAATTtgg                                                                                            <  1:3361522/91‑1 (MQ=255)
  tCCAGAAATTTTGAACGGATCTCTATCATTGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGAGCACCACGAATTtggt                                                                                           >  2:2853442/1‑91 (MQ=255)
  tCCAGAAATTTTGAACGGATCGCGATCATTGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGAGCACCACGAATTtggt                                                                                           <  2:306170/91‑1 (MQ=255)
   ccAGAAATTTTGAACGGATCTCTATCATTGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGAGCACCACGAATTtggtg                                                                                          >  2:2949171/1‑91 (MQ=255)
   ccAGAAATTTTGAACGGATCTCTATCATTGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGAGCACCACGAATTtggtg                                                                                          >  1:3021734/1‑91 (MQ=255)
    cAGAAATTTTGAACGGATCTCTATCATTGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGAGCACCACGAATTtggtgg                                                                                         <  1:1022387/91‑1 (MQ=255)
    cAGAAATTTTGAACGGATCTCTATCATTGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGAGCACCACGAATTtggtgg                                                                                         <  1:1506857/91‑1 (MQ=255)
     aGAAATTTTGAACGGATCTCTATCATTGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGAGCACCACGAATTTggtggt                                                                                        >  1:1708944/1‑91 (MQ=255)
     aGAAATTTTGAACGGATCTCTATCATTGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGAGCACCACGAATTTggtggt                                                                                        <  1:3523713/91‑1 (MQ=255)
      gAAATTTTGAACGGATCTCTATCATTGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGAGGCACCACCAAATTcgtggt                                                                                        <  2:3825896/91‑1 (MQ=38)
      gAAATTTTGAACGGATCTCTATCATTGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGAGCACCACGAATTTGgtggtg                                                                                       <  2:905593/91‑1 (MQ=255)
      gAAATTTTGAACGGATCTCTATCATTGCCTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGAGGCACCACCAAATTcgtggt                                                                                        <  2:510795/91‑1 (MQ=37)
        aaTTTTGAACGGATCTCTATCATTGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGAGGCACCACCAAATTCGTGGTGc                                                                                      <  1:175118/91‑1 (MQ=37)
        aaTTTTGAACGGATCTCTATCATTGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGAGCACCACGAATTTGGTGGTGcc                                                                                     <  2:2253535/91‑1 (MQ=255)
        aaTTTTGAACGGATCTCTATCATTGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGAGCACCACGAATTTGGTGGTGcc                                                                                     >  1:2939986/1‑91 (MQ=255)
         aTTTTGAACGGATCTCTATCATTGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGAGCACCACGAATTTGGTGGTGCct                                                                                    <  2:1474128/91‑1 (MQ=255)
             tGAACGGATCTCTATCATTGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGAGCACCACGAATTTGGTGGTGCCTCTAt                                                                               >  2:3180745/1‑91 (MQ=255)
              gAACGGATCTCTATCATTGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGAGGCACCACCAAATTCGTGGTGCTCTAtt                                                                              <  2:3280407/91‑1 (MQ=25)
              gAACGGATCTCTATCATTGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGAGCACCACGAATTTGGTGGTGCCTCTAtt                                                                              >  1:2581337/1‑91 (MQ=255)
               aaCGGATCTCTATCATTGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGAGCACCACGAATTTGGTGGTGCCTCTATTg                                                                             >  2:2197157/1‑91 (MQ=255)
                aCGGATCTCTATCATTGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGAGCACCACGAATTTGGTGGTGCCTCTATTGa                                                                            <  1:728662/91‑1 (MQ=255)
                aCGGATCTCTATCATTGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGAGCACCACGAATTTGGTGGTGCCTCTATTGa                                                                            >  2:2717917/1‑91 (MQ=255)
                aCGGATCTCTATCATTGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGAGCACCACGAATTTGGTGGTGCCTCTATTGa                                                                            >  1:1444419/1‑91 (MQ=255)
                aCGGATCTCTATCACTGCTTGAAGCTATAAGTAATAAATATAAGTTTAATTATCGAATAGAGCACCACCAAATTCGTGGTGCTCTAtttga                                                                           <  2:3677822/91‑4 (MQ=16)
                   gATCTCTATCATTGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGAGCACCACGAATTTGGTGGTGCCTCTATTGATAc                                                                         >  1:125484/1‑91 (MQ=255)
                    aTCTCTATCATTGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGAGCACCACGAATTTGGTGGTGCCTCTATTGATACa                                                                        <  1:3920109/91‑1 (MQ=255)
                     tctcTATCATTGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGAGCACCACGAATTTGGTGGTGCCTCTATTGATACAt                                                                       >  1:3860170/1‑91 (MQ=255)
                     tctcTATCATTGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGAGCACCACGAATTTGGTGGTGCCTCTATTGATACAt                                                                       <  1:996770/91‑1 (MQ=255)
                     tctcTATCATTGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGAGCACCACGAATTTGGTGGTGCCTCTATTGATACAt                                                                       >  2:409652/1‑91 (MQ=255)
                     tctcTATCATTGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGAGCACCACGAATTTGGTGGTGCCTCTATTGATACAt                                                                       <  1:3598111/91‑1 (MQ=255)
                         tATCATTGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGAGCACCACGAATTTGGTGGTGCCTCTATTGATACATTCgg                                                                   >  1:2523004/1‑91 (MQ=255)
                         tATCATTGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGAGCACCACGAATTTGGTGGTGCCTCTATTGATACATTCgg                                                                   >  2:2928689/1‑91 (MQ=255)
                          aTCATTGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGAGCACCACGAATTTGGTGGTGCCTCTATTGATACATTCGGc                                                                  >  2:389062/1‑91 (MQ=255)
                          aTCATTGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGAGCACCACGAATTTGGTGGTGCCTCTATTGATACATTCGGc                                                                  >  2:2319687/1‑91 (MQ=255)
                             aTTGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGAGCACCACGAATTTGGTGGTGCCTCTATTGATACATTCGGCGAg                                                               >  2:3387012/1‑91 (MQ=255)
                              ttGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGAGCACCACGAATTTGGTGGTGCCTCTATTGATACATTCGGCGAGc                                                              >  2:2309772/1‑91 (MQ=255)
                                gCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGAGCACCACGAATTTGGTGGTGCCTCTATTGATACATTCGGCGAGCCt                                                            <  1:2577048/91‑1 (MQ=255)
                                gCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGAGCACCACGAATTTGGTGGTGCCTCTATTGATACATTCGGCGAGCCt                                                            <  1:804225/91‑1 (MQ=255)
                                  ttGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGAGCACCACGAATTTGGTGGTGCCTCTATTGATACATTCGGCGAGCCttt                                                          >  1:463257/1‑91 (MQ=255)
                                   tGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGAGCACCACGAATTTGGTGGTGCCTCTATTGATACATTCGGCGAGCCTTTa                                                         >  1:2305761/1‑91 (MQ=255)
                                     aaaTTATAAGTAATAAATATAACTTTAATTATCAAATAGAGCACCACGAATTTGGTGGTGCCTCTATTGATACATTCGGCGAGCCTTTAAc                                                       >  2:1749700/1‑91 (MQ=255)
                                      aaTTATAAGTAATAAATATAACTTTAATTATCAAATAGAGCACCACGAATTTGGTGGTGCCTCTATTGATACATTCGGCGAGCCTTTAACt                                                      <  2:862877/91‑1 (MQ=255)
                                      aaTTATAAGTAATAAATATAACTTTAATTATCAAATAGAGCACCACGAATTTGGTGGTGCCTCTATTGATACATTCGGCGAGCCTTTAACt                                                      <  1:938891/91‑1 (MQ=255)
                                       aTTATAAGTAATAAATATAACTTTAATTATCAAATAGAGCACCACGAATTTGGTGGTGCCTCTATTGATACATTCGGCGAGCCTTTAACTg                                                     <  1:366608/91‑1 (MQ=255)
                                        ttATAAGTAATAAATATAACTTTAATTATCAAATAGAGCACCACGAATTTGGTGGTGCCTCTATTGATACATTCGGCGAGCCTTTAACTga                                                    <  1:1657558/91‑1 (MQ=255)
                                          ataAGTAATAAATATAACTTTAATTATCAAATAGAGCACCACGAATTTGGTGGTGCCTCTATTGATACATTCGGCGAGCCTTTAACTgaga                                                  <  1:3626866/91‑1 (MQ=255)
                                            aaGTAATAAATATAACTTTAATTATCAAATAGAGCACCACGAATTTGGTGGTGCCTCTATTGATACATTCGGCGAGCCTTTAACTGAGaaa                                                >  1:2160224/1‑91 (MQ=255)
                                             aGTAATAAATATAACTTTAATTATCAAATAGAGCACCACGAATTTGGTGGTGCCTCTATTGATACATTCGGCGAGCCTTTAACTGAGcaaa                                               >  2:3848706/1‑91 (MQ=255)
                                             aGTAATAAATATAACTTTAATTATCAAATAGAGCACCACGAATTTGGTGGTGCCTCTATTGATACATTCGGCGAGCCTTTAACTGAGaaaa                                               >  1:2492590/1‑91 (MQ=255)
                                                aataaataTAACTTTAATTATCAAATAGAGCACCACGAATTTGGTGGTGCCTCTATTGATACATTCGGCGAGCCTTTAACTGAGAAAACCt                                            <  1:2216033/91‑1 (MQ=255)
                                                   aaataTAACTTTAATTATCAAATAGAGCACCACGAATTTGGTGGTGCCTCTATTGATACATTCGGCGAGCCTTTAACTGAGAAAACCTTaa                                         <  1:1060265/91‑1 (MQ=255)
                                                      tataACTTTAATTATCAAATAGAGCACCACGAATTTGGTGGTGCCTCTATTGATACATTCGGCGAGCCTTTAACTGAGAAAACCTTAAATg                                      <  2:1203800/91‑1 (MQ=255)
                                                      tataACTTTAATTATCAAATAGAGCACCACGAATTTGGTGGTGCCTCTATTGATACATTCGGCGAGCCTTTAACTGAGAAAACCTTAAATg                                      <  2:1359328/91‑1 (MQ=255)
                                                       ataACTTTAATGATCAAATAGAGCACCACGAATTTGGTGGTGCCTCTATTGATACATTCGGCGAGCCTTTAACTGAGAAAACCTTAAATGc                                     <  1:3577538/91‑1 (MQ=255)
                                                        taACTTTAATTATCAAATAGAGCACCACGAATTTGGTGGTGCCTCTATTGATACATTCGGCGAGCCTTTAACTGAGAAAACCTTAAATGCg                                    <  2:3039409/91‑1 (MQ=255)
                                                        taACTTTAATTATCAAATAGAGCACCACGAATTTGGTGGTGCCTCTATTGATACATTCGGCGAGCCTTTAACTGAGAAAACCTTAAATGCg                                    <  1:1993048/91‑1 (MQ=255)
                                                         aaCTTTAATTATCAAATAGAGCACCACGAATTTGGTGGTGCCTCTATTGATACATTCGGCGAGCCTTTAACTGAGAAAACCTTAAATGCgt                                   <  1:3761062/91‑1 (MQ=255)
                                                          aCTTTAATTATCAAATAGAGCACCACGAATTTGGTGGTGCCTCTATTGATACATTCGGAGAGCCTCTAACTGAGAAAACCTTAAATGCgtg                                  >  2:1623694/1‑91 (MQ=255)
                                                           cTTTAATTATTAAATAGAGCACCACGAATTTGGTGGTGCCTCTATTGATACATTCGGCGAGCCTTTAACTGAGAAAACCTTAAATGCgtgt                                 >  2:1055697/1‑91 (MQ=255)
                                                           cTTTAATTATCAAATAGAGCACCACGAATTTGGTGGTGCCTCTATTGATACATTCGGCGAGCCTTTAACTGAGAAAACCTTAAATGCgtgt                                 >  1:1276026/1‑91 (MQ=255)
                                                            tttAATTATCAAATAGAGCACCACGAATTTGGTGGTGCCTCTATTGATACATTCGGCGAGCCTTTAACTGAGAAAACCTTAAATGCGTGTa                                >  1:137995/1‑91 (MQ=255)
                                                               aaTTATCAAATAGAGCACCACGAATTTGGTGGTGCCTCTATTGATACATTCGGCGAGCCTTTAACTGAGAAAACCTTAAATGCGTGTaaaa                             >  1:1145181/1‑91 (MQ=255)
                                                               aaTTATCAAATAGAGCACCACGAATTTGGTGGTGCCTCTATTGATACATTCGGCGAGCCTTTAACTGAGAAAACCTTAAATGCGTGTaaaa                             <  1:910506/91‑1 (MQ=255)
                                                                aTTATCAAATAGAGCACCACGAATTTGGTGGTGCCTCTATTGATACATTCGGCGAGCCTTTAACTGAGAAAACCTTAAATGCGTGGAAAag                            <  2:2434247/91‑1 (MQ=255)
                                                                 ttATCAAATAGAGCACCACGAATTTGGTGGTGCCTCTATTGATACATTCGGCGAGCCTTTAACTGAGAAAACCTTAAATGCGTGTAAAaga                           <  1:394434/91‑1 (MQ=255)
                                                                  tATCAAATAGAGCACCACGAATTTGTTGGTGCCGCTATTGTTACATTCGGCGAGCCTTTAACTGAAAAAACCTTAAATGCGTGTAAAagag                          <  1:3507782/91‑1 (MQ=255)
                                                                  tATCAAATAGAGCACCACGAATTTGGTGGTGCCTCTATTGATACCTTCGGCGAGCCTTTAACTGAGAAAACCTTAAATGCGTGTAAAcgag                          >  2:1957109/1‑91 (MQ=255)
                                                                     tcaaTAGAGGCACCACCAAATTCGTGGTGCTCTATTTGATACATTCGGCGAGCCTTTAACTGAGAAAACCTTAAATGCGTGTAAAAGAGCa                        >  1:711013/3‑91 (MQ=17)
                                                                     cAAATAGAGCACCACGAATTTGGTGGTGCCTCTATTGATACATTCGGCGAGCCTTTAACTGAGAAAACCTTAAATGCGTGTAAAAGAGCAg                       >  2:1114447/1‑91 (MQ=255)
                                                                     cAAATAGAGCACCACGAATTTGGTGGTGCCTCTATTGATACATTCGGCGAGCCTTTAACTGAGAAAACCTTAAATGCGTGTAAAAGAGCAg                       >  1:49484/1‑91 (MQ=255)
                                                                     cAAATAGAGCACCACGAATTTGGTGGTGCCTCTATTGATACATTCGGCGAGCCTTTAACTGAGAAAACCTTAAATGCGTGTAAAAGACCAg                       >  1:3194800/1‑91 (MQ=255)
                                                                      caaTAGAGGCACCACCAAATTCGAGGTGCTCTATTTGATACATCCGGCGAGCCTTTAACTGAGAAACCCTTAAATGCCTGTAAAAGAGCAg                       >  2:2145209/2‑91 (MQ=9)
                                                                        aTAGAGGCACCACCAAATTCGTGGTGCTCTATTTGATACATTCGGCGAGCCTTTAACTGAGAAAACCTTAAATGCGTGTAAAAGAGCAGAt                     >  2:3107851/1‑91 (MQ=18)
                                                                        aTAGAGCACCACGAATTTGGTGGTGCCTCTATTGATACATTCGGCGAGCCTTTAACTGAGAAAACCTTAAATGCGTGTAAAATAGCatctg                    <  1:571683/91‑5 (MQ=255)
                                                                         tAGAGCACCACGAATTTGGTGGTGCCTCTATTGATACATTCGGCGAGCCTTTAACTGAGAAAACCTTAAATGCGTGTAAAAGAGCAGATGc                   <  2:3771425/91‑1 (MQ=255)
                                                                          agagCACCACGAATTTGGTGGTGCCTCTATTGATACATTCGGCGAGCCTTTAACTGAGAACCCCTTAAATGCGTGTAAAATAGCATCTGCt                  <  2:2651595/91‑1 (MQ=255)
                                                                           gagCACCACGAATTTGGTGGTGCCTCTATTGATACATTCGGCGAGCATTTAACTGAGAAAGCCTTAAATGCGTGTAAAAGAGCAGATGCTa                 <  2:3873192/91‑1 (MQ=255)
                                                                            agCACCACGAATTTGGTGGTGCCTCTATTGATACATTCGGCGAGCCTTTAACTGAGAAAACCTTAAATGCGTGTAAAAGAGCAGATGCTAt                <  1:2338974/91‑1 (MQ=255)
                                                                             ggCACCACCAAATTCGTGGTGCTCTATTGATACATTCGGCGAGCCTTTAACTGAGAAAACCTTAAATGCGTGTAAAAGAGCAGATGCTAtt               >  1:1842209/2‑91 (MQ=37)
                                                                              gCACCACGAATTTGGTGGTGCCTCTATTGCTACATTCGGCGAGCCTTTAACTGAGAAAACCTTAAATGCGTGTAAAAGAGCAGATGCTAtt               >  2:3446694/1‑91 (MQ=255)
                                                                              gCACCACGAATTTGGTGGTGCCTCTATTGATACATTCGGCGAGCCTTTAACTGAGAAAACCTTAAATGCGTGTAAAAGAGCAGATGCTAtt               <  2:626283/91‑1 (MQ=255)
                                                                              gCACCACCAAATTCGTGGTGCTCTATTGATACATTCGGCGAGCCTTTAACTGAGAAAACCTTAAATGCGTGTAAAAGAGCAGATGCTAttt              >  1:2808540/1‑91 (MQ=38)
                                                                               caccacGAATTTGGTGGTGCCTCTATTGATACATTCGGCGAGCCTTTAACTGAGAAAACCTTAAATGCGTGTAAAAGAGCAGATGCTAttt              <  2:2373912/91‑1 (MQ=255)
                                                                                accacGAATTTGGTGGTGCCTCTATTGATACATTCGGCGAGCCTTTAACTGAGACAACCTTAAATGCGTGTAAAAGAGCAGATGCTAtttt             >  2:3118295/1‑91 (MQ=255)
                                                                                accacGAATTTGGTGGTGCCTCTATTGATACATTCGGCGAGCCTTTAACTGAGAAAACCTTAAATGCGTGTAAAAGAGCAGATGCTAtttt             >  2:2972482/1‑91 (MQ=255)
                                                                                 ccacGAATTTGGTGGTGCCTCTATTGATACATTCGGCGAGCCTTTAACTGAGAAAACCTTAAATGCGTGTAAAAGAGCAGATGCTATTTTa            >  1:1323796/1‑91 (MQ=255)
                                                                                 ccacGAATTTGGTGGTGCCTCTATTGATACATTCGGCGAGCCTTTAACTGAGAAAACCTTAAATGCGTGTAAAAGAGCAGATGCTATTTTa            >  1:2289436/1‑91 (MQ=255)
                                                                                 ccacGAATTTGGTGGTGCCTCTATTGATACATTCGGCGAGCCTTTAACTGAGAAAACCTTAAATGCGTGTAAAAGAGCAGATGCTATTTTa            >  2:1097410/1‑91 (MQ=255)
                                                                                  cacGAATTTGGTGGTGCCTCTATTGATACATTCGGCGAGCCTTTAACTGAGAAAACCTTAAATGCGTGTAAAAGAGCAGATGCTATTTTAc           >  1:3548333/1‑91 (MQ=255)
                                                                                    cGAATTTGGTGGTGCCTCTATTGATACATTCGGCGAGCCTTTAACTGAGAAAACCTTAAATGCGTGTAAAAGAGCAGATGCTATTTTACTg         >  2:688512/1‑91 (MQ=255)
                                                                                       aTTTGGTGGTGCCTCTATTGATACATTCGGCGAGCCTTTAACTGAGAAAACCTTAAATGCGTGTAAAAGAGCAGATGCTATTTTACTGGGt      <  1:2868693/91‑1 (MQ=255)
                                                                                        tttGGTGGTGCCTCTATTGATACATTCGGCGAGCCTTTAACTGAGAAAACCTTAAATGCGTGTAAAAGAGCAGATGCTATTTTACTGGGTg     <  2:1590935/91‑1 (MQ=255)
                                                                                        tttGGTGGTGCCTCTATTGATACATTCGGCGAGCCTTTAACTGAGAAAACCTTAAATGCGTGTAAAAGAGCAGATGCTATTTTACTGGGTg     <  2:1218244/91‑1 (MQ=255)
                                                                                         ttGGTGGAGCCTCTATTGATACATTCGGCGAGCCTTTAACTGAGAAAACCTTAAATGCGTGTAAAAGAGCAGATGCTATTTTACTGGGTGc    <  1:1955536/91‑1 (MQ=255)
                                                                                           ggtggtGCCTCTATTGATACATTCGGCGAGCCTTTAACTGAGAAAACCTTAAATGCGTGTAAAAGAGCAGATGCTATTTTACTGGGTGCaa  >  2:2020686/1‑91 (MQ=255)
                                                                                           ggtggtGCCTCTATTGATACATTCGGCGAGCCTTTAACTGAGAAAACCTTAAATGCGTGTAAAAGAGCAGATGCTATTTTACTGGGTGCaa  >  2:1715869/1‑91 (MQ=255)
                                                                                           ggtggtGCCTCTATTGATACATTCGGCGAGCCTTTAACTGAGAAAACCTTAAATGCGTGTAAAAGAGCAGATGCTATTTTACTGGGTGCaa  <  1:987840/91‑1 (MQ=255)
                                                                                           |                                                                                           
GGTCCAGAAATTTTGAACGGATCTCTATCATTGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGAGCACCACGAATTTGGTGGTGCCTCTATTGATACATTCGGCGAGCCTTTAACTGAGAAAACCTTAAATGCGTGTAAAAGAGCAGATGCTATTTTACTGGGTGCAA  >  USA300TCH1516_ALE/2171991‑2172171

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