Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 1,343,619 A→C 14.5% intergenic (+32/‑98) USA300HOU_RS06640 → / → USA300TCH1516_01248 hypothetical protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE1,343,6190AC14.5% 921.1 / 153.8 344intergenic (+32/‑98)USA300HOU_RS06640/USA300TCH1516_01248hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (133/160);  new base C (13/37);  total (147/197)
Fisher's exact test for biased strand distribution p-value = 1.28e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.95e-01

Only 100 of 358 total aligned reads displayed.

TATACATAAATACTAGTAAAGGGCTATTTCCTGAATCTGAACTCGTAGCATGGGGATAAAGAAAGTTTACTTTAAAAATTTTAACCTTCAAATGTGAAGGTTATTTTTTATTGCCCAAAACGTGCTGATGACATTTTAAAAGCAAGTATGGAATATCAGTCGACAGACTATAAACGGAGGT  >  USA300TCH1516_ALE/1343529‑1343709
                                                                                          |                                                                                          
tataCATAAATACTAGTAAAGGGCTATTTCCTGAATCTGAACTCGTAGCATGGGGATAAAGAAAGTTTACTTTAAAAATTTTAACCTTCaa                                                                                            <  2:3157142/91‑1 (MQ=255)
   aCATAAATACTAGTAAAGGGCTATTTCCTGAATCTGAACTCGTAGCATGGGGATAAAGAAAGTTTACTTTAAAAATTTTAACCTTCAAAtg                                                                                         <  1:3217208/91‑1 (MQ=255)
   aCATAAATACTAGTAAAGGGCTATTTCCTGAATCTGAACTCGTAGCATGGGGATAAAGAAAGTTTACTTTAAAAATTTTAACCTTCAAAtg                                                                                         >  1:2964432/1‑91 (MQ=255)
    cATAAATACTAGTAAAGGGCTATTTCCTGAATCTGAACTCGTAGCATGGGGATAAAGAAAGTTTACTTTAAAAATTTTAACCTTCAAAtgt                                                                                        <  1:370618/91‑1 (MQ=255)
     aTAAATACTAGTAAAGGGCTATTTCCTGAATCTGAACTCGTAGCATGGGGATAAAGAAAGTTTACTTTAAAAATTTTAACCTTCAAAtgtg                                                                                       <  1:1462641/91‑1 (MQ=255)
      tAAATACTAGTAAAGGGCTATTTCCTGAATCTGAACTCGTAGCATGGGGATAAAGAAAGTTTACTTTAAAAATTTTAACCTTCACATTTGa                                                                                      <  1:1395848/91‑1 (MQ=255)
      tAAATACTAGTAAAGGGCTATTTCCTGAATCTGAACTCGTAGCATGGGGATAAAGAAAGTTTACTTTAAAAATTTTAACCTTCACATTTGa                                                                                      <  1:682807/91‑1 (MQ=255)
      tAAATACTAGTAAAGGGCTATTTCCTGAATCTGAACTCGTAGCATGGGGATAAAGAAAGTTTACTTTAAAAATTTTAACCTTCAAATGTGa                                                                                      <  1:3831424/91‑1 (MQ=255)
       aaaTACTAGTAAAGGGCTATTTCCTGAATCTGAACTCGTAGCATGGGGATAAAGAAAGTTTACTTTAAAAATTTTAACCTTCAAATGTGaa                                                                                     >  2:2290921/1‑91 (MQ=255)
       aaaTACTAGTAAAGGGCTATTTCCTGAATCTGAACTCGTAGCATGGGGATAAAGAAAGTTTACTTTAAAAATTTTAACCTTCAAATGTGaa                                                                                     <  2:3372217/91‑1 (MQ=255)
        aaTACTAGTAAAGGGCTATTTCCTGAATCTGAACTCGTAGCATGGGGATAAAGAAAGTTTACTTTAAAAATTTTAACCTTCAAATGTGAAg                                                                                    <  1:2423302/91‑1 (MQ=255)
         aTACTAGTAAAGGGCTATTTCCTGAATCTGAACTCGTAGCATGGGGATAAAGAAAGTTTACTTTAAAAATTTTAACCTTCACATTTGAAgg                                                                                   <  2:1051584/91‑1 (MQ=255)
          tACTAGTAAAGGGCTATTTCCTGAATCTGAACTCGTAGCATGGGGATAAAGAAAGTTTACTTTAAAAATTTTAACCTTCACATTTGAAGGt                                                                                  <  2:2140504/91‑1 (MQ=255)
          tACTAGTAAAGGGCTATTTCCTGAATCTGAACTCGTAGCATGGGGATAAAGAAAGTTTACTTTAAAAATTTTAACCTTCACATTTGAAGGt                                                                                  <  1:2361812/91‑1 (MQ=255)
          tACTAGTAAAGGGCTATTTCCTGAATCTGAACTCGTAGCATGGGGATAAAGAAAGTTTACTTTAAAAATTTTAACCTTCAAATGTGAAGGt                                                                                  <  1:1637284/91‑1 (MQ=255)
           aCTAGTAATGGGCTATTGCCTGAATCTTAACTCGTAGCATGGGGATAGAGAAAGTTGACTTTAAAAATTTTAACCTTCACATTTGATGGtt                                                                                 <  2:3454506/91‑1 (MQ=255)
           aCTAGTAAAGGGCTATTTCCTGAATCTGAACTCGTAGCATGGGGATAAAGAAAGTTTACTTTAAAAATTTTAACCTTCACATTTGAAGGtt                                                                                 <  1:608520/91‑1 (MQ=255)
           aCTAGTAAAGGGCTATTTCCTGAATCTGAACTCGTAGCATGGGGATAAAGAAAGTTTACTTTAAAAATTTTAACCTTCAAATGTGAAGGtt                                                                                 >  2:3800572/1‑91 (MQ=255)
            cTAGTAAAGGGCTATTTCCTGAATCTGAACTCGTAGCATGGGGATAAAGAAAGTTTACTTTAAAAATTTTAACCTTCACATTTGAAGGTTa                                                                                <  1:1951325/91‑1 (MQ=255)
            cTAGTAAAGGGCTATTTCCTGAATCTGAACTCGTAGCATGGGGATAAAGAAAGTTTACTTTAAAAATTTTAACCTTCACATTTGAAGGTTa                                                                                <  2:1256188/91‑1 (MQ=255)
            cTAGTAAAGGGCTATTTCCTGAATCTGAACTCGTAGCATGGGGATAAAGAAAGTTTACTTTAAAAATTTTAACCTTCAAATGTGAAGGTTa                                                                                >  1:2082653/1‑91 (MQ=255)
            cTAGTAAAGGGCTATTTCCTGAATCTGAACTCGTAGCATGGGGATAAAGAAAGTTTACTTTAAAAATTTTAACCTTCAAATGTGAAGGTTa                                                                                >  1:3305767/1‑91 (MQ=255)
             tAGTAAAGGGCTATTTCCTGAATCTGAACTCGTAGCATGGGGATAAAGAAAGTTTACTTTAAAAATTTTAACCTTCACATTTGAAGGTTaa                                                                               <  1:506538/91‑2 (MQ=255)
              aGTAAAGGGCTATTTCCTGAATCTGAACTCGTAGCATGGGGATAAAGAAAGTTTACTTTAAAAATTTTAACCTTCACATTTGAAGGTTaaa                                                                              <  1:1224414/91‑3 (MQ=39)
              aGTAAAGGGCTATTTCCTGAATCTGAACTCGTAGCATGGGGATAAAGAAAGTTTACTTTAAAAATTTTAACCTTCAAATGTGAAGGTTAtt                                                                              >  1:3153369/1‑91 (MQ=255)
              aGTAAAGGGCTATTTCCTGAATCTGAACTCGTAGCATGGGGATAAAGAAAGTTTACTTTAAAAATTTTAACCTTCAAATGTGAAGGTTAtt                                                                              >  1:405200/1‑91 (MQ=255)
              aGTAAAGGGCTATTTCCTGAATCTGAACTCGTAGCATGGGGATAAAGAAAGTTTACTTTAAAAATTTTAACCTTCAAATGTGAAGGTTAtt                                                                              >  2:1425714/1‑91 (MQ=255)
              aGTAAAGGGCTATTTCCTGAATCTGAACTCGTAGCATGGGGATAAAGAAAGTTTACTTTAAAAATTTTAACCTTCAAATGTGAAGGTTAtt                                                                              <  1:2499162/91‑1 (MQ=255)
                tAAAGGGCTATTTCCTGAATCTGAACTCGTAGCATGGGGATAAAGAAAGTTTACTTTAAAAATTTTAACCTTCAAATGTGAAGGTTAtttt                                                                            >  1:379772/1‑91 (MQ=255)
                 aaaGGGCTATTTCCTGAATCTGAACTCGTAGCATGGGGATAAAGAAAGTTTACTTTAAAAATTTTAACCTTCAAATGTGAAGGTTAttttt                                                                           <  1:105410/91‑1 (MQ=255)
                   aGGGCTATTTCCTGAATCTGAACTCGTAGCATGGGGATAAAGAAAGTTTACTTTAAAAATTTTAACCTTCAAATGTGAAGGTTATTTTTTa                                                                         <  1:968395/91‑1 (MQ=255)
                      gCTATTTCCTGAATCTGAACTCGTAGCATGGGGATAAAGAAAGTTTACTTTAAAAATTTTAACCTTCAAATGTGAAGGTTATTTTTTATTg                                                                      >  2:1285239/1‑91 (MQ=255)
                      gCTATTTCCTGAATCTGAACTCGTAGCATGGGGATAAAGAAAGTTTACTTTAAAAATTTTAACCTTCAAATGTGAAGGTTATTTTTTATTg                                                                      >  1:3205266/1‑91 (MQ=255)
                       cTATTTCCTGAATCTGAACTCGTAGCATGGGGATAAAGAAAGTTTACTTTAAAAATTTTAACCTTCAAATGTGAAGGTTATTTTTTATTGc                                                                     <  2:1453460/91‑1 (MQ=255)
                       cTATTTCCTGAATCTGAACTCGTAGCATGGGGATAAAGAAAGTTTACTTTAAAAATTTTAACCTTCAAATGTGAAGGTTATTTTTTATTGc                                                                     <  1:3870261/91‑1 (MQ=255)
                       cTATTTCCTGAATCTGAACTCGTAGCATGGGGATAAAGAAAGTTTACTTTAAAAATTTTAACCTTCAAATGTGAAGGTTATTTTTTATTGc                                                                     >  2:1211928/1‑91 (MQ=255)
                        tATTTCCTGAATCTGAACTCGTAGCATGGGGATAAAGAAAGTTTACTTTAAAAATTTTAACCTTCAAATGTGAAGGTTATTTTTTATTGcc                                                                    <  1:207355/91‑1 (MQ=255)
                        tATTTCCTGAATCTGAACTCGTAGCATGGGGATAAAGAAAGTTTACTTTAAAAATTTTAACCTTCAAATGTGAAGGTTATTTTTTATTGcc                                                                    <  2:1558036/91‑1 (MQ=255)
                        tATTTCCTGAATCTGAACTCGTAGCATGGGGATAAAGAAAGTTTACTTTAAAAATTTTAACCTTCAAATGTGAAGGTTATTTTTTATTGcc                                                                    >  2:3959639/1‑91 (MQ=255)
                         atgtCCTGAATCTGAACTCGTAGCAGGGGGATAAAGAAAGTTTACTTTAAAAATTTTAGCCTTCAAATGTGAAGGTTATTTTTTATTGccc                                                                   <  2:279280/88‑1 (MQ=255)
                         aTTTCCTGAATCTGAACTCGTAGCATGGGGATAAAGAAAGTTTACTTTAAAAATTTTAACCTTCAAATGTGAAGGTTATTTTTTATTGccc                                                                   >  1:3470496/1‑91 (MQ=255)
                         aTTTCCTGAATCTGAACTCGTAGCATGGGGATAAAGAAAGTTTACTTTAAAAATTTTAACCTTCAAATGTGAAGGTTATTTTTTATTGccc                                                                   >  2:3699583/1‑91 (MQ=255)
                         aTTTCCTGAATCAGAACTCGTAGCATGGGGATAAAGAAAGTTTACTTTAAAAATTTTAACCTTCAAATGTGAAGGTTATTTTTTATTGccc                                                                   <  2:783859/91‑1 (MQ=255)
                          tttCCTGAATCTGAACTCGTAGCATGGGGATAAAGAAAGTTTACTTTAAAAATTTTAACCTTCAAATGTGAAGGTTATTTTTTATTGCCCa                                                                  <  1:547479/91‑1 (MQ=255)
                            tCCTGAATCTGAACTCGTAGCATGGGGATAAAGAAAGTTTACTTTAAAAATTTTAACCTTCAAATGTGAAGGTTATTTTTTATTGCCCaaa                                                                <  1:1900722/91‑1 (MQ=255)
                             ccTGAATCTGAACTCGTAGCATGGGGATAAAGAAAGTTTACTTTAAAAATTTTAACCTTCAAATGTGAAGGTTATTTTTTATTGCCCaaaa                                                               >  2:1691584/1‑91 (MQ=255)
                               tGAATCTGAACTCGTAGCATGGGGATAAAGAAAGTTTACTTTAAAAATTTTAACCTTCAAATGTGAAGGTTATTTTTTATTGCCCAAAACg                                                             <  1:2209977/91‑1 (MQ=255)
                                    cTGATCTCGTAGCATGGGGATAAAGAAAGTTTACTTTAAAAATTTTAACCTTCAAATGTGAAGGTTATTTTTTATTGCCCAAAACGTGCtg                                                        <  1:3295353/91‑1 (MQ=255)
                                    cTGAACTCGTAGCATGGGGATAAAGAAAGTTTACTTTAAAAATTTTAACCTTCAAATGTGAAGGTTATTTTTTATTGCCCAAAACGTGCtg                                                        <  1:2168883/91‑1 (MQ=255)
                                    cTGAACTCGTAGCATGGGGATAAAGAAAGTTTACTTTAAAAACTTTAACCTTCAAATGTGAAGGTTATTTTTTATTGCCCAAAACGTGCtg                                                        <  2:1082213/91‑1 (MQ=255)
                                     tGAACTCGTAGCATGGGGATAAAGAAAGTTTACTTTAAAAATTTTAACCTTCAAATGTGAAGGTTATTTTTTATTGCCCAAAACGTGCtga                                                       >  1:1611628/1‑91 (MQ=255)
                                       aaCTCGTAGCATGGGGATAAAGAAAGTTTACTTTAAAAATTTTAACCTTCAAATGTGAAGGTTATTTTTTATTGCCCAAAACGTGCtgatg                                                     <  1:3562603/91‑1 (MQ=255)
                                       aaCTCGTAGCATGGGGATAAAGAAAGTTTACTTTAAAAATTTTAACCTTCAAATGTGAAGGTTATTTTTTATTGCCCAAAACGTGCtgatg                                                     <  1:2042885/91‑1 (MQ=255)
                                        acgcGTAGCATGGGGATAAAGAAAGTTTACTTTAAAAATTTTAACCTTCAAATGTGAAGGTTATTTTTTATTGCCCAAAACGTGCtgatga                                                    <  2:3002406/88‑1 (MQ=255)
                                        aCTCGTAGCATGGGGATAAAGAAAGTTTACTTTAAAAATTTTAACCTTCAAATGTGAAGGTTATTTTTTATTGCCCAAAACGTGCtgatga                                                    <  1:2891645/91‑1 (MQ=255)
                                        aCTCGTAGCATGGGGATAAAGAAAGTTTACTTTAAAAATTTTAACCTTCAAATGTGAAGGTTATTTTTTATTGCCCAAAACGTGCtgatga                                                    >  1:2734842/1‑91 (MQ=255)
                                        aCTCGTAGCATGGGGATAAAGAAAGTTTACTTTAAAAATTTTAACCTTCAAATGTGAAGGTTATTTTTTATTGCCCAAAACGTGCtgatga                                                    >  1:2307667/1‑91 (MQ=255)
                                          tCGTAGCATGGGGATAAAGAAAGTTTACTTTAAAAATTTTAACCTTCAAATGTGAAGGTTATTTTTTATTGCCCAAAACGTGCTGATGACa                                                  <  1:2253462/91‑1 (MQ=255)
                                           cGTAGCATGGGGATAAAGAAAGTTTACTTTAAAAATTTTAACCTTCAAATGTGAAGGTTATTTTTTATTGCCCAAAACGTGCTGATGACAt                                                 >  2:2457473/1‑91 (MQ=255)
                                           cGTAGCATGGGGATAAAGAAAGTTTACTTTAAAAATTTTAACCTTCAAATGTGAAGGTTATTTTTTATTGCCCAAAACGTGCTGATGACAt                                                 >  2:503897/1‑91 (MQ=255)
                                            gTAGCATGGGGATAAAGAAAGTTTACTTTAAAAATTTTAACCTTCAAATGTGAAGGTTATTTTTTATTGCCCAAAACGTGCTGATGACAtt                                                >  2:2871166/1‑91 (MQ=255)
                                            gTAGCATGGGGATAAAGAAAGTTTACTTTAAAAATTTTAACCTTCAAATGTGAAGGTTATTTTTTATTGCCCAAAACGTGCTGATGACAtt                                                >  1:1540599/1‑91 (MQ=255)
                                             tAGCATGGGGATAAAGAAAGTTTACTTTAAAAATTTTAACCTTCAAATGTGAAGGTTATTTTTTATTGCCCAAAACGTGCTGATGACAttt                                               <  1:3942685/91‑1 (MQ=255)
                                             tAGCATGGGGATAAAGAAAGTTTACTTTAAAAATTTTAACCTTCAAATGTGAAGGTTATTTTTTATTGCCCAAAACGTGCTGATGACAttt                                               <  1:812231/91‑1 (MQ=255)
                                                  tGGGGATAAAGAAAGTTTACTTTAAAAATTTTAACCTTCAAATGTGAAGGTTATTTTTTATTGCCCAAAACGTGCTGATGACATTTTaaaa                                          <  1:2651763/91‑1 (MQ=255)
                                                  tGGGGATAAAGAAAGTTTACTTTAAAAATTTTAACCTTCAAATGTGAAGGTTATTTTTTATTGCCCAAAACGTGCTGATGACATTTTaaaa                                          >  2:2220865/1‑91 (MQ=255)
                                                   gggtATAAAGAAAGTTTACTTTAAAAATTTTAACCTTCAAATGTGAAGGTTATTTTTTATTGCCCAAAACGTGCTGATGACATTTTAAAAg                                         <  1:3066490/91‑1 (MQ=255)
                                                   ggggATAAAGAAAGTTTACTTTAAAAATTTTAACCTTCAAATGTGAAGGTTATTTTTTATTGCCCAAAACGTGCTGATGACATTTTAAAAg                                         <  1:2843472/91‑1 (MQ=255)
                                                    gggATAAAGAAAGTTTACTTTAAAAATTTTAACCTTCAAATGTGAAGGTTATTTTTTATTGCCCAAAACGTGCTGATGACATTTTAAAAGc                                        <  1:2573345/91‑1 (MQ=255)
                                                        tAAAGAAAGTTTACTTTAAAAATTTTAACCTTCAAATGTGAAGGTTATTTTTTATTGCCCAAAACGTGCTGATGACATTTTAAAAGCAAGt                                    <  1:741619/91‑1 (MQ=255)
                                                        tAAAGAAAGTTTACTTTAAAAATTTTAACCTTCAAATGTGAAGGTTATTTTTTATTGCCCAAAACGTGCTGATGACATTTTAAAAGCAAGt                                    >  1:2024400/1‑91 (MQ=255)
                                                        tAAAGAAAGTTTACTTTAAAAATTTTAACCTTCAAATGTGAAGGTTATTTTTTATTGCCCAAAACGTGCTGATGACATTTTAAAAGCAAGt                                    <  2:1942994/91‑1 (MQ=255)
                                                         aaagaaagTTTACTTTAAAAATTTTAACCTTCAAATGTGAAGGTTATTTTTTATTGCCCAAAACGTGCTGATGACATTTGAAAAGCAAGTa                                   <  2:3229560/91‑1 (MQ=255)
                                                          aagaaagTTTACTTTAAAAATTTTAACCTTCAAATGTGAAGGTTATTTTTTATTGCCCAAAACGTGCTGATGACATTTTAAAAGCAAGTAt                                  <  1:2775397/91‑1 (MQ=255)
                                                           agaaagTTTACTTTAAAAATTTTAACCTTCAAATGTGAAGGTTATTTTTTATTGCCCACAACGTGCTGATGACATTTTAAAAGCAAGTATg                                 >  2:2516493/1‑91 (MQ=255)
                                                            gaaagTTTACTTTAAAAATTTTAACCTTCAAATGTGAAGGTTATTTTTTATTGCCCAAAACGTGCTGATGACATTTTAAAAGCAAGTATgg                                >  1:2721266/1‑91 (MQ=255)
                                                            gaaagTTTACTTTAAAAATTTTAACCTTCAAATGTGAAGGTTATTTTTTATTGCCCAAAACGTGCTGATGACATTTTAAAAGCAAGTATgg                                >  1:383751/1‑91 (MQ=255)
                                                               agTTTACTTTAAAAATTTTAACCTTCAAATGTGAAGGTTATTTTTTATTGCCCAAAACGTGCTGATGACATTTTAAAAGCAAGTATGGAat                             <  1:1117274/91‑1 (MQ=255)
                                                               agTTTACTTTAAAAATTTTAACCTTCAAATGTGAAGGTTATTTTTTATTGCCCAAAACGTGCTGATGACATTTTAAAAGCAAGTATGGAat                             <  2:159328/91‑1 (MQ=255)
                                                                   tACTTTAAAAATTTTAACCTTCAAATGTGAAGGTTATTTTTTATTGCCCAAAACGTGCTGATGACATTTTAAAAGCAAGTATGGAATATCa                         <  2:2352975/91‑1 (MQ=255)
                                                                    aCTTTAAAAATTTTAACCTTCAAATGTGAAGGTTATTTTTTATTGCCCAAAACGTGCTGATGACATTTTAAAAGCAAGTATGGAATATCAg                        >  1:1761914/1‑91 (MQ=255)
                                                                     cTTTAAAAATTTTAACCTTCAAATGTGAAGGTTATTTTTTATTGCCCAAAACGTGCTGATGACATTTTAAAAGCAAGTATGGAATATCAGt                       >  2:1640042/1‑91 (MQ=255)
                                                                       ttAAAAATTTTAACCTTCAAATGTGAAGGTTATTTTTTATTGCCCAAAACGTGCTGATGACATTTTAAAAGCAAGTATGGAATGTCAGTCg                     <  1:116776/91‑1 (MQ=255)
                                                                          aaaaTTTTAACCTTCAAATGTGAAGGTTATTTTTTATTGCCCAAAACGTGCTGTTGACATTTTAAAAGCAAGTATGGAATATCAGTCGACa                  >  1:3658741/1‑91 (MQ=255)
                                                                          aaaaTTTTAACCTTCAAATGTGAAGGTTATTTTTTATTGCCCAAAACGTGCTGATGACATTTTAAAAGCAAGTATGGAATATCAGTCGACa                  >  2:3587743/1‑91 (MQ=255)
                                                                             aTTTTAACCTTCAAATGTGAAGGTTATTTTTTATTGCCCAAAACGTGCTGATGACATTTTAAAAGCAAGTATGGAATATCAGTCGACAGAc               >  2:3344058/1‑91 (MQ=255)
                                                                               tttAACCTTCAAATGTGAAGGTTCTTTTTTATTGCCCAAAACGTGCTGATGATATTTCAAAAGCAAGTATGGAATATCAGTCCACAGAGta             >  2:1777891/1‑91 (MQ=255)
                                                                               tttAACCTTCAAATGTGAAGGTTATTTTTTATTGCCCAAAACGTGCTGATGACATTTTAAAAGCAAGTATGGAATATCAGTCGACAGACta             <  2:2610770/91‑1 (MQ=255)
                                                                               tttAACCTTCAAATGTGAAGGTTATTTTTTATTGCCCAAAACGTGCTGATGACATTTTAAAAGCAAGTATGGAATATCAGTCGACAGACta             <  2:2101747/91‑1 (MQ=255)
                                                                                atAACCTTCACATTTGAAGGTTATTTTTTATTGCCCAAAACGTGCTGATGACATTTTAAAAGCAAGTATGTAATATCAGTCGACAGACtat            >  2:698788/2‑91 (MQ=255)
                                                                                atAACCTTCACATTTGAAGGTTATTTTTTATTGCCCAAAACGTGCTGATGACATTTTAAAAGCAAGTATGGAATATCAGTCGACAGACtat            >  1:130404/2‑91 (MQ=255)
                                                                                  aaCCTTCACATTTGAAGGTTATTTTTTATTGCCCAAAACGTGCTGATGACATTTTAAAAGCAAGTATGGAATATCAGTCGACAGACTATaa          >  1:2525016/1‑91 (MQ=255)
                                                                                  aaCCTTCAAATGTGAAGGTTATTTTTTATTGCCCAAAACGTGCTGATGACATTTTAAAAGCAAGTATGGAATATCAGTCGACAGACTATaa          >  1:2382489/1‑91 (MQ=255)
                                                                                  aaCCTTCAAATGTGAAGGTTATTTTTTATTGCCCAAAACGTGCTGATGACATTTTAAAAGCAAGTATGGAATATCAGTCGACAGACTATaa          <  2:3895034/91‑1 (MQ=255)
                                                                                    ccTTCAAATGTGAAGGTTATTTTTTATTGCCCAAAACGTGCTGATGACATTTTAAAAGCAAGTATGGAATATCAGTCGACAGACTATAAAc        >  1:1341712/1‑91 (MQ=255)
                                                                                     cTTCACATTTGAAGGTTATTTTTTATTGCCCAAAACGTGCTGATGACATTTTAAAAGCAAGTATGGAATATCAGTCGACAGACTATAAACg       >  1:636699/1‑91 (MQ=255)
                                                                                     cTTCACATTTGAAGGATATTTTTTATTGCCCAAAACGTGCTGATGACATTTTAAAAGCAAGTATGGAATATCAGTCGACAGACTATAAACg       >  1:723909/1‑91 (MQ=255)
                                                                                         aaaTGTGAAGGTTATTTTTTATTGCCCAAAACGTGCTGATGACATTTTAAAAGCAAGTATGGAATATCAGTCGACAGACTATAAACGGAgg   <  2:60940/91‑1 (MQ=255)
                                                                                         aaaTGTGAAGGTTATTTTTTATTGCCCAAAACGTGCTGATGACATTTTAAAAGCAAGTATGGAATATCAGTCGACAGACTATAAACGGAgg   <  1:2303800/91‑1 (MQ=255)
                                                                                         aaaTGTGAAGGTTATTTTTTATTGCCCAAAACGTGCTGATGACATTTTAAAAGCAAGTATGGAATATCAGTCGACAGACTATAAACGGAgg   <  2:2014629/91‑1 (MQ=255)
                                                                                          |                                                                                          
TATACATAAATACTAGTAAAGGGCTATTTCCTGAATCTGAACTCGTAGCATGGGGATAAAGAAAGTTTACTTTAAAAATTTTAACCTTCAAATGTGAAGGTTATTTTTTATTGCCCAAAACGTGCTGATGACATTTTAAAAGCAAGTATGGAATATCAGTCGACAGACTATAAACGGAGGT  >  USA300TCH1516_ALE/1343529‑1343709

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