Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 2,424,608 A→G 8.1% intergenic (‑166/+62) USA300HOU_RS12490 ← / ← USA300HOU_RS12495 hypothetical protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE2,424,6080AG 8.1% 388.4 / 18.9 124intergenic (‑166/+62)USA300HOU_RS12490/USA300HOU_RS12495hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (47/67);  new base G (8/2);  total (55/69)
Fisher's exact test for biased strand distribution p-value = 2.25e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 128 total aligned reads displayed.

CCTGTCGATGATTGATTCACAATAGCACGTGATTTTTTTGTTTTTCAATATTATTCATTTATTCCATCAAAAACACCCTTTTTAATTTTTACAAAAATTAAAAAAAGTGCTCCTACACTGCTTGCATGTAGAAACACTTTTTCATTGTAATGTTATTCTTCTCGAGACATACCTTTTAGCA  >  USA300TCH1516_ALE/2424518‑2424698
                                                                                           |                                                                                          
ccTGTCGATGATTGATTCACAATAGCACGTGATTTTTTTGTCTTTCAATATTATTCATTTATTCCATCAAAAACACCCTTTTTAATTTTTa                                                                                            >  1:3599826/1‑91 (MQ=255)
 cTGTCGATGATTGATTCACAATAGCACGTGATTTTTTTGTTTTTCAATATTATTCATTTATTCCATCAAAAACACCCTTTTTAATTTTTAc                                                                                           <  1:3914682/91‑1 (MQ=255)
    tCGATGATTGATTCACAATAGCACGTGATTTTTTTTGTTTTTCAATATTATTCATTTATTCCATCAAAAACACCCTTTTTAATTTTTACaa                                                                                         <  1:4044036/91‑1 (MQ=255)
     cGATGATTGATTCACAATAGCACGTGATTTTTTTGTTTTTCAATATTATTCATTTATTCCATCAAAAACACCCTTTTTAATTTTTACaaaa                                                                                       <  1:2238558/91‑1 (MQ=255)
      gatgatTGATTCACAATAGCACGTGATTTTTTTGTTTTTCAATATTATTCATTTATTCCATCAAAAACACCCTTTTTAATTTTTACaaaaa                                                                                      >  1:590914/1‑91 (MQ=255)
      gatgatTGATTCACAATAGCACGTGATTTTTTTGTTTTTCAATATTATTCATTTATTCCATCAAAAACACCCTTTTTAATTTTTACaaaaa                                                                                      <  1:3088112/91‑1 (MQ=255)
      gatgatTGATTCACAATAGCACGTGATTTTTTTGTTTTTCAATATTATTCATTTATTCCATCAAAAACACCCTTTTTAATTTTTACaaaaa                                                                                      <  1:2690434/91‑1 (MQ=255)
      gatgatTGATTCACAATAGCACGTGATTTTTTTGTTTTTCAATATTATTCATTTATTCCATCAAAAACACCCTTTTTAATTTTTACaaaaa                                                                                      >  1:1549575/1‑91 (MQ=255)
      gatgatTGATTCACAATAGCACGTGATTTTTTTGTTTTTCAATATTATTCATTTATTCCATCAAAAACACCCTTTTTAATTTTTACaaaaa                                                                                      <  1:2138985/91‑1 (MQ=255)
        tgattgatTCACAATAGCACGTGATTTTTTTGTTTTTCAATATTATTCATTTATTCCATCAAAAACACCCTTTTTAATTTTTGTAAAAAtt                                                                                    <  1:1463510/91‑1 (MQ=255)
         gattgattCACAATAGCACGTGATTTTTTTGTTTTTCAATATTATTCATTTATTCCATCAAAAACACCCTTTTTAATTTTTACAAAAATTa                                                                                   <  1:4001426/91‑1 (MQ=255)
            tgattCACAATAGCACGTGATTTTTTTGTTTTTCAATATTATTCATTTATTCCATCAAAAACACCCTTTTTAATTTTTACAAAAATTaaaa                                                                                <  1:108757/91‑1 (MQ=255)
             gattCACAATAGCACGTGATTTTTTTGTTTTTCAATATTATTCATTTATTCCATCAAAAACACCCTTTTTAATTTTTACAAAAATTaaaaa                                                                               >  1:3464059/1‑91 (MQ=255)
                tCACAATAGCACGTGATTTTTTTGTTTTTCAATATTATTCATTTATTCCATCAAAAACACCCTTTTTAATTTTTACAAAAATTAAAAAAAg                                                                            >  1:1231150/1‑91 (MQ=255)
                tCACAATAGCACGTGATTTTTTTGTTTTTCAATATTATTCATTTATTCCATCAAAAACACCCTTTTTAATTTTTACAAAAATTAAAAAAAg                                                                            >  1:1192778/1‑91 (MQ=255)
                 cacaATAGCACGTGATTTTTTTGTTTTTCAATATTATTCATTTATTCCATCAAAAACACCCTTTTTAATTTTTGTAAAAATTaaaaagggt                                                                           <  1:535953/91‑5 (MQ=38)
                  acaATAGCACGTGATTTTTTTGTTTTTCAATATTATTCATTTATTCCATCAAAAACACCCTTTTTAATTTTTACAAAAATTAAAAAAAGTg                                                                          <  1:973526/91‑1 (MQ=255)
                    aaTAGCACGTGATTTTTTTGTTTTTCAATATTATTCATTTATTCCATCAAAAACACCCTTTTTAATTTTTACAAAAATTAAGAAAAGTGCt                                                                        >  1:3976707/1‑91 (MQ=255)
                    aaTAGCACGTGATTTTTTTGTTTTTCAATATTATTCATTTATTCCATCAAAAACACCCTTTTTAATTTTTACAAAAATTAAAAAAAGTGcc                                                                        >  1:3445266/1‑90 (MQ=255)
                     aTAGCACGTGATTTTTTTGTTTTTCAATATTATTCATTTATTCCATCAAAAACACCCTTTTTAATTTTTACAAAAATTAAAAAAAGTGCTc                                                                       >  1:3976962/1‑91 (MQ=255)
                     aTAGCACGTGATTTTTTTGTTTTTCAATATTATTCATTTATTCCATCAAAAACACCCTTTTTAATTTTTACAAAAATTAAAAAAAGTGCTc                                                                       >  1:3092393/1‑91 (MQ=255)
                      tAGCACGTGATTTTTTTGTTTTTCAATATTATTCATTTATTCCATCAAAAACACCCTTTTTAATTTTTACAAAAATTAAAAAAAGTGCTcc                                                                      >  1:2936542/1‑91 (MQ=255)
                       aGCACGTGATTTTTTTGTTTTTCAATATTATTCATTTATTCCATCAAAAACACCCTTTTTAATTTTTACAAAAATTAAAAAAAGTGCTCCt                                                                     <  1:2904561/91‑1 (MQ=255)
                        gCACGTGATTTTTTTGTTTTTCAATATTATTCATTTATTCCATCAAAAACACCCTTTTTAATTTTTACAAAAATTAAAAAAAGTGCTCCTa                                                                    >  1:2636395/1‑91 (MQ=255)
                        gCACGTGATTTTTTTGTTTTTCAATATTATTCATTTATTCCATCAAAAACACCCTTTTTAATTTTTACAAAAATTAAAAAAAGTGCTCCTa                                                                    >  1:1553972/1‑91 (MQ=255)
                        gCACGTGATTTTTTTGTTTTTCAATATTATTCATTTATTCCATCAAAAACACCCTTTTTAATTTTTACAAAAATTAAAAAAAGTGCTCCTa                                                                    >  1:3770472/1‑91 (MQ=255)
                        gCACGTGATTTTTTTGTTTTTCAATATTATTCATTTATTCCATCAAAAACACCCTTTTTAATTTTTACAAAAATTAAAAAAAGTGCTCCTa                                                                    >  1:871485/1‑91 (MQ=255)
                         cACGTGATTTTTTTGTTTTTCAATATTATTCATTTATTCCATCAAAAACACCCTTTTTAATTTTTACAAAAATTAAAAAAAGTGCTCCTac                                                                   <  1:672573/91‑1 (MQ=255)
                         cACGTGATTTTTTTGTTTTTCAATATTATTCATTTATTCCATCAAAAACACCCTTTTTAATTTTTACAAAAATTAAAAAAAGTGCTCCTac                                                                   <  1:2990075/91‑1 (MQ=255)
                           cGTGATTTTTTTGTTTTTCAATATTATTCATTTATTCCATCAAAAACACCCTTTTTAATTTTTACAAAAATTAAAAAAAGTGCTCCTacac                                                                 <  1:2387948/91‑1 (MQ=255)
                            gTGATTTTTTTGTTTTTCAATATTATTCATTTATTCCATCAAAAACACCCTTTTTAATTTTTACAAAAATTAAAAAAAGTGCTCCTACACt                                                                <  1:3986458/91‑1 (MQ=255)
                             tGATTTTTTTGTTTTTCAATATTATTCATTTATTCCATCAAAAACACCCTTTTTAATTTTTACAAAAATTAAAAAAAGTGCTCCTACACTg                                                               <  1:1652271/91‑1 (MQ=255)
                              gATTTTTTTGTTTTTCAATATTATTCATTTATTCCATCAAAAACACCCTTTTTAATTTTTACAAAAATTAAAAAAAGTGCTCCTACACTGc                                                              <  1:3797955/91‑1 (MQ=255)
                                 tttttttGTTTTTCAATATTATTCATTTATTCCATCAAAAACACCCTTTTTAATTTTTACAAAAATTAAAAAAAGTGCTCCTACACTGCtt                                                            <  1:1328464/91‑1 (MQ=255)
                                   tttttGTTTTTCAATATTATTCATTTATTCCATCAAAAACACCCTTTTTAATTTTTACAAAAATTAAAAAAAGTGCTCCTACACTGCTTGc                                                          <  1:358864/91‑1 (MQ=255)
                                   tttttGTTTTTCAATATTATTCATTTATTCCATCAAAAACACCCTTTTTAATTTTTACAAAAATTAAAAAAAGTGCTCCTACACTGCTTGc                                                          <  1:4115001/91‑1 (MQ=255)
                                   tttttGTTTTTCAATATTATTCATTTATTCCATCAAAAACACCCTTTTTAATTTTTACAAAAATTAAAAAAAGTGCTCCTACACTGCTTGc                                                          >  1:86050/1‑91 (MQ=255)
                                   tttttGTTTTTCAATATTATTCATTTATTCCATCAAAAACACCCTTTTTAATTTTTACAAAAATTAAAAAAAGTGCTCCTACACTGCTTGc                                                          <  1:2131198/91‑1 (MQ=255)
                                     tttGTTTTTCAATATTATTCATTTATTCCATCAAAAACACCCTTTTTAATTTTTACAAAAATTAAAAAAAGTGCTCCTACACTGCTTGCAt                                                        <  1:2749269/91‑1 (MQ=255)
                                      ttGTTTTTCAATATTATTCATTTATTCCATCAAAAACACCCTTTTTAATTTTTACAAAAATTAAAAAAAGTGCTCCTACACTGCTTGCATg                                                       <  1:3951820/91‑1 (MQ=255)
                                        gTTTTTCAATATTATTCATTTATTCCATCAAAAACACCCTTTTTAATTTTTACAAAAATTAAAAAAAGTGCTCCTACACTGCTTGCATGTa                                                     <  1:3935907/91‑1 (MQ=255)
                                        gTTTTTCAATATTATTCATTTATTCCATCAAAAACACCCTTTTTAATTTTTACAAAAATTAAAAAAAGTGCTCCTACACTGCTTGCATGTa                                                     <  1:3329492/91‑1 (MQ=255)
                                        gTTTTTCAATATTATTCATTTATTCCATCAAAAACACCCTTTTTAATTTTTACAAAAATTAAAAAAAGTGCTCCTACACTGCTTGCATGTa                                                     >  1:1087318/1‑91 (MQ=255)
                                          ttttCAATATTATTCATTTATTCCATCAAAAACACCCTTTTTAATTTTTACAAAAATTAAAAAAAGTGCTCCTACACTGCTTGCATGTAGa                                                   <  1:3783007/91‑1 (MQ=255)
                                           tttCAATATTATTCATTTATTCCATCAAAAACACCCTTTTTAATTTTTACAAAAATTAAAAAAAGTGCTCCTACACTGCTTGCATGTAGaa                                                  >  1:827536/1‑91 (MQ=255)
                                              cAATATTATTCATTTATTCCATCAAAAACACCCTTTTTAATTTTTACAAAAATTAAAAAAAGTGCTCCTACACTGCTTGCATGTAGAAaca                                               >  1:2808962/1‑91 (MQ=255)
                                               aaTATTATTCATTTATTCCATCAAAAACACCCTTTTTAATTTTTACAAAAATTAAAAAAAGTGCTCCTACACTGCTTGCATGTAGAAacac                                              <  1:898711/91‑1 (MQ=255)
                                               aaTATTATTCATTTATTCCATCAAAAACACCCTTTTTAATTTTTACAAAAATTAAAAAAAGTGCTCCTACACTGCTTGCATGTAGAAacac                                              <  1:2572191/91‑1 (MQ=255)
                                                 tattatTCATTTATTCCATCAAAAACACCCTTTTTAATTTTTACAAAAATTAAAAAAAGTGCTCCTACACTGCTTGCATGTAGAAACACtt                                            <  1:1759315/91‑1 (MQ=255)
                                                  attattCATTTATTCCATCAAAAACACCCTTTTTAATTTTTACAAAAATTAAAAAAAGTGCTCCTACACTGCTTGCATGTAGAAACACttt                                           >  1:2688273/1‑91 (MQ=255)
                                                  attattCATTTATTCCATCAAAAACACCCTTTTTAATTTTTACAAAAATTAAAAAAAGTGCTCCTACACTGCTTGCATGTAGAAACACttt                                           >  1:3295213/1‑91 (MQ=255)
                                                    tattCATTTATTCCATCAAAAACACCCTTTTTAATTTTTACAAAAATTAAAAAAAGTGCTCCTACACTGCTTGCATGTAGAAACACttttt                                         >  1:1032420/1‑91 (MQ=255)
                                                     attCATTTATTCCATCAAAAACACCCTTTTTAATTTTTACAAAAATTAAAAAAAGTGCTCCTACACTGCTTGCATGTAGAAACACTTTTTc                                        >  1:2236288/1‑91 (MQ=255)
                                                      ttCATTTATTCCATCAAAAACACCCTTTTTAATTTTTACAAAAATTAAAAAAAGTGCTCCTACGCTGCTTGCATGTAGAAACACTTTTTCa                                       <  1:2491721/91‑1 (MQ=255)
                                                          tttATTCCATCAAAAACACCCTTTTTAATTTTTACAAAAATTAAAAAAAGTGCTCCTACACTGCTTGCATGTAGAAACACTTTTTCATTGt                                   <  1:2237752/91‑1 (MQ=255)
                                                          tttATTCCATCAAAAACACCCTTTTTAATTTTTACAAAAATTAAAAAAAGTGCTCCTACACTGCTTGCATGTAGAAACACTTTTTCATTGt                                   <  1:2389476/91‑1 (MQ=255)
                                                           ttATTCCATCAAAAACACCCTTTTTAATTTTTACAAAAATTAAAAAAAGTGCTCCTACACTGCTTGCATGTAGAAACACTTTTTCATTGTa                                  <  1:1280237/91‑1 (MQ=255)
                                                              ttCCATCAAAAACACCCTTTTTAATTTTTACAAAAATTAAAAAAAGTGCTCCTACACTGCTTGCATGTAGAAACACTTTTTCATTGTAATg                               >  1:2421508/1‑91 (MQ=255)
                                                              ttCCATCAAAAACACCCTTTTTAATTTTTACAAAAATTAAAAAAAGTGCTCCTACACTGCTTGCATGTAGAAACACTTTTTCATTGTAATg                               >  1:4004370/1‑91 (MQ=255)
                                                              ttCCATCAAAAACACCCTTTTTAATTTTTACAAAAATTAAAAAAAGTGCTCCTACACTGCTTGCATGTAGAAACACTTTTTCATTGTAATg                               <  1:2177416/91‑1 (MQ=255)
                                                              ttCCATCAAAAACACCCTTTTTAATTTTTACAAAAATTAAAAAAAGTGCTCCTACACTGCTTGCATGTAGAAACACTTTTTCATTGTAATg                               >  1:2357404/1‑91 (MQ=255)
                                                               tCCATCAAAAACACCCTTTTTAATTTTTACAAAAATTAAAAAAAGTGCTCCTACACTGCTTGCATGTAGAAACACTTTTTCATTGTAATGt                              <  1:1876629/91‑1 (MQ=255)
                                                                ccATCAAAAACACCCTTTTTAATTTTTACAAAAATTAAAAAAAGTGCTCCTACACTGCTTGCATGTAGAAACACTTTTTCATTGTAATGtt                             <  1:2751430/91‑1 (MQ=255)
                                                                ccATCAAAAACACCCTTTTTAATTTTTACAAAAATTAAAAAAAGTGCTCCTACACTGCTTGCATGTAGAAACACTTTTTCATTGTAATGtt                             <  1:1603643/91‑1 (MQ=255)
                                                                ccATCAAAAACACCCTTTTTAATTTTTACAAAAATTAAAAAAAGTGCTCCTACACTGCTTGCATGTAGAAACACTTTTTCATTGTAATGtt                             <  1:2515874/91‑1 (MQ=255)
                                                                  aTCAAAAACACCCTTTTTAATTTTTACAAAAATTAAAAAAAGTGCTCCTACACTGCTTGCATGTAGAAACACTTTTTCATTGTAATGTTAt                           <  1:794683/91‑1 (MQ=255)
                                                                  aTCAAAAACACCCTTTTTAATTTTTACAAAAATTAAAAAAAGTGCTCCTACACTGCTTGCATGTAGAAACACTTTTTCATTGTAATGTTAt                           <  1:1477738/91‑1 (MQ=255)
                                                                     aaaaaCACCCTTTTTAATTTTTACAAAAATTAAAAAAAGTGCTCCTACACTGCTTGCATGTAGAAACACTTTTTCATTGTAATGTTAttct                        <  1:1481287/91‑1 (MQ=255)
                                                                     aaaaaCACCCTTTTTAATTTTTACAAAAATTAAAAAAAGTGCTCCTACACTGCTTGCATGTAGAAACACTTTTTCATTGTAATGTTAttct                        <  1:1049505/91‑1 (MQ=255)
                                                                     aaaaaCACCCTTTTTAATTTTTACAAAAATTAAAAAAAGTGCTCCTACACTGCTTGCATGTAGAAACACTTTTTCATTGTAATGTTAttct                        <  1:1481410/91‑1 (MQ=255)
                                                                       aaaCACCCTTTTTAATTTTTACAAAAATTAAAAAAAGTGCTCCTACACTGCTTGCATGTAGAAACACTTTTTCATTGTAATGTTAttcttc                      >  1:2064511/1‑91 (MQ=255)
                                                                       aaaCACCCTTTTTAATTTTTACAAAAATTAAAAAAAGTGCTCCTACACTGCTTGCATGTAGAAACACTTTTTCATTGTAATGTTAttcttc                      >  1:2064559/1‑91 (MQ=255)
                                                                       aaaCACCCTTTTTAATTTTTACAAAAATTAAAAAAAGTGCTCCTACACTGCTTGCATGTAGAAACACTTTTTCATTGTAATGTTAttcttc                      >  1:1179163/1‑91 (MQ=255)
                                                                       aaaCACCCTTTTTAATTTTTACAAAAATTAAAAAAAGTGCTCCTACACTGCTTGCATGTAGAAACACTTTTTCATTGTAATGTTAttcttc                      >  1:1178418/1‑91 (MQ=255)
                                                                            cccTTTTTAATTTTTACAAAAATTAAAAAAAGTGCTCCTACACTGCTTGCATGTAGAAACACTTTTTCATTGTAATGTTATTCTTCTCgag                 <  1:3614894/91‑1 (MQ=255)
                                                                              ttttttGATTTTTGTAAAAATTAAAAAGGGTGCTCCTACACTGCTTGCATGTAGAAACACTTTTTCATTGTAATGTTATTCTTCTCGAGAc               >  1:3232733/2‑91 (MQ=37)
                                                                              ttttttAATTTTTGTAAAAATTAAAAAGGGTGCTCCTACACTGCTTGCATGTAGAAACACTTTTTCATTGTAATGTTATTCTTCTCGAGAc               >  1:3232461/2‑91 (MQ=37)
                                                                              ttttttAATTTTTGTAAAAATTAAAAAGGGTGCTCCTACACTGCTTGCATGTAGAAACACTTTTTCATTGTAATGTTATTCTTCTCGAGAc               >  1:1732960/2‑91 (MQ=37)
                                                                              cTTTTTAATTTTTACAAAAATTAAAAAAAGTGCTCCTACACTGCTTGCATGTAGAAACACTTTTTCATTGTAATGTTATTCTTCTCGAGAc               <  1:4052119/91‑1 (MQ=255)
                                                                              cTTTTTAATTTTTACAAAAATTAAAAAAAGTGCTCCTACACTGCTTGCATGTAGAAACACTTTTTCATTGTAATGTTATTCTTCTCGAGAc               <  1:4052199/91‑1 (MQ=255)
                                                                              cTTTTTAATTTTTACAAAAATTAAAAAAAGTGCTCCTACACTGCTTGCATGTAGAAACACTTTTTCATTGTAATGTTATTCTTCTCGAGAc               <  1:624395/91‑1 (MQ=255)
                                                                               tttttAATTTTTACAAAAATTAAAAAAAGTGCTCCTACACTGCTTGCATGTAGAAACACTTTTTCATTGTAATGTTATTCTTCTCGAGACa              <  1:356783/91‑1 (MQ=255)
                                                                               tttttAATTTTTACAAAAATTAAAAAAAGTGCTCCTACACTGCTTGCATGTAGAAACACTTTTTCATTGTAATGTTATTCTTCTCGAGACa              >  1:3660217/1‑91 (MQ=255)
                                                                                ttttAATTTTTACAAAAATTAAAAAAAGTGCTCCTACACTGCTTGCATGTAGAAACACTTTTTCATTGTAATGTTATTCTTCTCGAGACAt             >  1:3123097/1‑91 (MQ=255)
                                                                                ttttAATTTTTACAAAAATTAAAAAAAGTGCTCCTACACTGCTTGCATGTAGAAACACTTTTTCATTGTAATGTTATTCTTCTCGAGACAt             <  1:4066375/91‑1 (MQ=255)
                                                                                ttttAATTTTTACAAAAATTAAAAAAAGTGCTCCTACACTGCTTGCATGTAGAAACACTTTTTCATTGTAATGTTATTCTTCTCGAGACAt             >  1:579465/1‑91 (MQ=255)
                                                                                ttttAATTTTTACAAAAATTAAAAAAAGTGCTCCTACACTGCTTGCATGTAGAAACACTTTTTCATTGTAATGTTATTCTTCTCGAGACAt             >  1:578983/1‑91 (MQ=255)
                                                                                 tttAATTTTTGTAAAAATTAAAAAGGGTGCTCCTACACTGCTTGCATGTAGAAACACTTTTTCATTGTAATGTTATTCTTCTCGAGACATa            >  1:209106/1‑91 (MQ=37)
                                                                                   tAATTTTTACAAAAATTAAAAAAAGTGCTCCTACACTGCTTGCATGTAGAAACACTTTTTCATTGTAATGTTATTCTTCTCGAGACATAcc          >  1:826489/1‑91 (MQ=255)
                                                                                    aaTTTTTGTAAAAATTAAAAAGGGTGCTCCTACACTGCTTGCATGTAGAAACACTTTTTCATTGTAATGTTATTCTTCTCGAGACATACCt         >  1:2838347/1‑91 (MQ=38)
                                                                                    aaTTTTTGTAAAAATTAAAAAAAGTGCTCCTACACTGCTTGCATGTAGAAACACTTTTTCATTGTAATGTTATTCTTCTCGAGACATACCt         >  1:441544/1‑91 (MQ=255)
                                                                                    aaTTTTTGTAAAAATTAAAAAAAGTGCTCCTACACTGCTTGCATGTAGAAACACTTTTTCATTGTAATGTTATTCTTCTCGAGACATACCt         >  1:3810148/1‑91 (MQ=255)
                                                                                    aaTTTTTACAAAAATTAAAAAAAGTGCTCCTACACTGCTTGCATGTAGAAACACTTTTTCATTGTAATGTTATTCTTCTCGAGACATACCt         >  1:1180389/1‑91 (MQ=255)
                                                                                    aaTTTTTACAAAAATTAAAAAAAGTGCTCCTACACTGCTTGCATGTAGAAACACTTTTTCATTGTAATGTTATTCTTCTCGAGACATACCt         <  1:1791914/91‑1 (MQ=255)
                                                                                      tttttACAAAAATTAAAAAAAGTGCTCCTACACTGCTTGCATGTAGAAACACTTTTTCATTGTAATGTTATTCTTCTCGAGACATACCttt       <  1:3878155/91‑1 (MQ=255)
                                                                                      tttttACAAAAATTAAAAAAAGTGCTCCTACACTGCTTGCATGTAGAAACACTTTTTCATTGTAATGTTATTCTTCTCGAGACATACCttt       >  1:3427672/1‑91 (MQ=255)
                                                                                      tttttACAAAAATTAAAAAAAGTGCTCCTACACTGCTTGCATGTAGAAACACTTTTTCATTGTAATGTTATTCTTCTCGAGACATACCttt       <  1:3877619/91‑1 (MQ=255)
                                                                                          tACAAAAATTAAAAAAAGTGCTCCTACACTGCTTGCATGTAGAAACACTTTTTCATTGTAATGTTATTCTTCTCGAGACATACCTTTTAGc   <  1:693020/91‑1 (MQ=255)
                                                                                           aCAAAAATTAAAAAAAGTGCTCCTACACTGCTTGCATGTAGAAACACTTTTTCATTGTAATGTTATTCTTCTCGAGACATACCTTTTAGCa  <  1:2906540/91‑1 (MQ=255)
                                                                                           aCAAAAATTAAAAAAAGTGCTCCTACACTGCTTGCATGTAGAAACACTTTTTCATTGTAATGTTATTCTTCTCGAGACATACCTTTTAGCa  <  1:2372004/91‑1 (MQ=255)
                                                                                           |                                                                                          
CCTGTCGATGATTGATTCACAATAGCACGTGATTTTTTTGTTTTTCAATATTATTCATTTATTCCATCAAAAACACCCTTTTTAATTTTTACAAAAATTAAAAAAAGTGCTCCTACACTGCTTGCATGTAGAAACACTTTTTCATTGTAATGTTATTCTTCTCGAGACATACCTTTTAGCA  >  USA300TCH1516_ALE/2424518‑2424698

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