Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 2,424,608 A→G 6.4% 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 6.4% 768.4 / 27.6 250intergenic (‑166/+62)USA300HOU_RS12490/USA300HOU_RS12495hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (117/114);  new base G (13/3);  total (133/117)
Fisher's exact test for biased strand distribution p-value = 2.00e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 252 total aligned reads displayed.

CCTGTCGATGATTGATTCACAATAGCACGTGATTTTTTTGTTTTTCAATATTATTCATTTATTCCATCAAAAACACCCTTTTTAATTTTTACAAAAATTAAAAAAAGTGCTCCTACACTGCTTGCATGTAGAAACACTTTTTCATTGTAATGTTATTCTTCTCGAGACATACCTTTTAGC  >  USA300TCH1516_ALE/2424518‑2424697
                                                                                          |                                                                                         
  tGTCGATGATTGATTCACAATAGCACGTGATTTTTTTGTTTTTCAATATTATTCATTTATTCCATCAAAAACACCCTTTTTAATTTTTACa                                                                                         <  2:586250/91‑1 (MQ=255)
  tGTCGATGATTGATTCACAATAGCACGTGATTTTTTTGTTTTTCAATATTATTCATTTATTCCATCAAAAACACCCTTTTTAATTTTTACa                                                                                         <  1:3070075/91‑1 (MQ=255)
  tGTCGATGATTGATTCACAATAGCACGTGATTTTTTTGTTTTTCAATATTATTCATTTATTCCATCAAAAACACCCTTTTTAATTTTTACa                                                                                         <  1:3068351/91‑1 (MQ=255)
   gTCGATGATTGATTCACAATAGCACGTGATTTTTTTGTTTTTCAATATTATTCATTTATTCCATCAAAAACACCCTTTTTAATTTTTACaa                                                                                        >  2:1889922/1‑91 (MQ=255)
    tCGATGATTGATTCACAATAGCACGTGATTTTTTTGTTTTTCAATATTATTCATTTATTCTATCAAACACACCCTTTTTATTTTTTATaaa                                                                                       >  1:1245664/1‑91 (MQ=255)
      gatgatTGATTCACAATAGCACGTGATTTTTTTGTTTTTCAATATTATTCATTTATTCCATCAAAAACACCCTTTTTAATTTTTACaaaaa                                                                                     <  1:3425509/91‑1 (MQ=255)
      gatgatTGATTCACAATAGCACGTGATTTTTTTGTTTTTCAATATTATTCATTTATTCCATCAAAAACACCCTTTTTAATTTTTACaaaaa                                                                                     >  2:3559440/1‑91 (MQ=255)
      gatgatTGATTCACAATAGCACGTGATTTTTTTGTTTTTCAATATTATTCATTTATTCCATCAAAAACACCCTTTTTAATTTTTACaaaaa                                                                                     <  1:2863172/91‑1 (MQ=255)
      gatgatTGATTCACAATAGCACGTGATTTTTTTGTTTTTCAATATTATTCATTTATTCCATCAAAAACACCCTTTTTAATTTTTACaaaaa                                                                                     <  1:2864079/91‑1 (MQ=255)
      gatgatTGATTCACAATAGCACGTGATTTTTTTGTTTTTCAATATTATTCATTTATTCCATCAAAAACACCCTTTTTAATTTTTACaaaaa                                                                                     <  1:192428/91‑1 (MQ=255)
      gatgatTGATTCACAATAGCACGTGATTTTTTTGTTTTTCAATATTATTCATTTATTCCATCAAAAACACCCTTTTTAATTTTTACaaaaa                                                                                     >  1:2169468/1‑91 (MQ=255)
      gatgatTGATTCACAATAGCACGTGATTTTTTTGTTTTTCAATATTATTCATTTATTCCATCAAAAACACCCTTTTTAATTTTTACaaaaa                                                                                     >  1:2168152/1‑91 (MQ=255)
      aatgatTGATTCACAATAGCACGTGATTTTTTTGTTTTTCAATATTATTCATTTATTCCATCAAAAACACCCTTTTTAATTTTTACaaaaa                                                                                     <  2:1572319/90‑1 (MQ=255)
        tgattgatTCACGATAGCACGTGATTTTTTTGTTTTTCAATATTATTCATTTATTCCATCAAAAACACCCTTTTTATTTTTTACAAAAAtt                                                                                   >  1:2973998/1‑91 (MQ=255)
        tgattgatTCACAATAGCACGTGATTTTTTTGTTTTTCAATATTATTCATTTATTCCATCAAAAACACCCTTTTTAATTTTTACAAAAAtt                                                                                   >  2:1623805/1‑91 (MQ=255)
        tgattgatTCACAATAGCACGTGATTTTTTTGTTTTTCAATATTATTCATTTATTCCATCAAAAACACCCTTTTTAATTTTTACAAAAAtt                                                                                   >  1:3001366/1‑91 (MQ=255)
         gattgattCACAATAGCACGTGATTTTTTTGTTTTTCAATATTATTCATTTATTCCATCAAAAACACCCTTTTTAATTTTTACAAAAATTa                                                                                  >  1:1768937/1‑91 (MQ=255)
            tgattCACAATAGCACGTGATTTTTTTGTTTTTCAATATTATTCATTTATTCCATCAAAAACACCCTTTTTAATTTTTACAAAAATTaaaa                                                                               <  1:30374/91‑1 (MQ=255)
             gattCACAATAGCACGTGATTTTTTTGTTTTTCAATATTATTCATTTATTCCATCAAAAACACCCTTTTTAATTTTTACAAAAATTaaaaa                                                                              <  2:1850713/91‑1 (MQ=255)
               ttCACAATAGCACGTGATTTTTTTGTTTTTTAATATTATTCATTTATTCCATCAAAAACACCCTTTTTAATTTTTACAAAAATTaaaaaaa                                                                            >  2:3045709/1‑91 (MQ=255)
               ttCACAATAGCACGTGATTTTTTTGTTTTTCAATATTATTCATTTATTCCATCAAAAACACCCTTTTTAATTTTTACAAAAATTaaaaaaa                                                                            <  1:1865425/91‑1 (MQ=255)
               ttCACAATAGCACGTGATTTTTTTGTTTTTCAATATTATTCATTTATTCCATCAAAAACACCCTTTTTAATTTTTACAAAAATTaaaaaaa                                                                            <  1:3670248/91‑1 (MQ=255)
                 cacaATAGCACGTGATTTTTTTGTTTTTCAATATTATTCATTTATTCCATCAAAAACACCCTTTTTAATTTTTACAAAAATTAAAAAAAGt                                                                          >  1:3749534/1‑91 (MQ=255)
                 cacaATAGCACGTGATTGTTCTGTTCTTCAATATTATTCATTTATTCCATCAAAAACACCCTTTTTAATTTTTGTAAAAATTaaaaagggt                                                                          <  2:834228/91‑5 (MQ=37)
                  acaATAGCACGTGATTTTTTTGTTTTTCAATATTATTCATTTATTCCATCAAAAACACCCTTTTTAATTTTTACAAAAATTAAAAAAAGTg                                                                         <  2:1174973/91‑1 (MQ=255)
                    aaTAGCACGTGATTTTTTTGTTTTTCAATATTATTCATTTATTCCATCAAAAACACCCTTTTTAATTTTTACAAAAATTAAAAAAAGTGCt                                                                       >  2:1024607/1‑91 (MQ=255)
                     aTAGCACGTGATTTTTTTGTTTTTCAATATTATTCATTTATTCCATCAAAAACACCCTTTTTAATTTTTACAAAAATTAAAAAAAGTGCTc                                                                      >  1:1275692/1‑91 (MQ=255)
                       aGCACGTGATTTTTTTTTTTTTTTTTATTATTCATTTATTCCATCAAAAACACCCTTTTTAATTTTTACAAAAATTAAAAAAAGTGCTCCt                                                                    >  2:1769523/1‑91 (MQ=255)
                       aGCACGTGATTTTTTTGTTTTTCAATATTATTCATTTATTCCATCAAAAACACCCTTTTTAATTTTTAGAAAAATTAAAAAAAGTGCTCCt                                                                    >  2:3124600/1‑91 (MQ=255)
                       aGCACGTGATTTTTTTGTTTTTCAATATTATTCATTTATTCCATCAAAAACACCCTTTTTAATTTTTACAAAAATTAAAAAAAGTGCTCCt                                                                    >  2:1112300/1‑91 (MQ=255)
                       aGCACGTGATTTTTTTGTTTTTCAATATTATTCATTTATTCCATCAAAAACACCCTTTTTAATTTTTACAAAAATTAAAAAAAGTGCTCCt                                                                    <  2:2757304/91‑1 (MQ=255)
                       aGCACGTGATTTTTTTGTTTTTCAATATTATTCATTTATTCCATCAAAAACACCCTTTTTAATTTTTACAAAAATTAAAAAAAGTGCTCCt                                                                    <  2:3299414/91‑1 (MQ=255)
                       aGCACGTGATTTTTTTGTTTTTCAATATTATTCATTTATTCCATCAAAAACACCCTTTTTAATTTTTACAAAAATTAAAAAAAGTGCTCCt                                                                    >  1:323500/1‑91 (MQ=255)
                        gCACGTGATTTTTTTGTTTTTCAATATTATTCATTTATTCCATCAAAAACACCCTTTTTAATTTTTACAAAAATTAAAAAAAGTGCTCCTa                                                                   <  2:3405180/91‑1 (MQ=255)
                        gCACGTGATTTTTTTGTTTTTCAATATTATTCATTTATTCCATCAAAAACACCCTTTTTAATTTTTACAAAAATTAAAAAAAGTGCTCCTa                                                                   >  1:3862803/1‑91 (MQ=255)
                         cACGTGATTTTTTTGTTTTTCAATATTATTCATTTATTCCATCAAAAACACCCTTTTTAATTTTTACAAAAATTAAAAAAAGTGCTCCTac                                                                  >  1:2600711/1‑91 (MQ=255)
                         cACGTGATTTTTTTGTTTTTCAATATTATTCATTTATTCCATCAAAAACACCCTTTTTAATTTTTACAAAAATTAAAAAAAGTGCTCCTac                                                                  >  1:3165146/1‑91 (MQ=255)
                          aCGTGATTTTTTTGTTTTTCAATATTATTCATTTATTCCATCAAAAACACCCTTTTTAATTTTTACAAAAATTAAAAAAAGTGCTCCTaca                                                                 <  1:980452/91‑1 (MQ=255)
                           cGTGATTTTTTTGTTTTTCAATATTATTCATTTATTCCATCAAAAACACCCTTTTTAATTTTTACAAAAATTAAAAAAAGTGCTcc                                                                     <  1:1982418/86‑1 (MQ=255)
                           cGTGATTTTTTTGTTTTTCAATATTATTCATTTATTCCATCAAAAACACCCTTTTTAATTTTTACAAAAATTAAAAAAAGTGCTCCTacac                                                                >  1:654272/1‑91 (MQ=255)
                           cGTGATTTTTTTGTTTTTCAATATTATTCATTTATTCCATCAAAAACACCCTTTTTAATTTTTACAAAAATTAAAAAAAGTGCTCCTacac                                                                <  1:3228750/91‑1 (MQ=255)
                           cGTGATTTTTTTGTTTTTCAATATTATTCATTTATTCCATCAAAAACACCCTTTTTAATTTTTACAAAAATTAAAAAAAGTGCTCCTacac                                                                >  2:1777915/1‑91 (MQ=255)
                           cGTGATTTTTTTGTTTTTCAATATTATTCATTTATTCCATCAAAAACACCCTTTTTAATTTTTACAAAAATTAAAAAAAGTGCTCCTacac                                                                >  2:2007865/1‑91 (MQ=255)
                           cGTGATTTTTTTGTTTTTCAATATTATTCATTTATTCCATCAAAAACACCCTTTTTAATTTTTACAAAAATTAAAAAAAGTGCTCCTacac                                                                <  2:785956/91‑1 (MQ=255)
                             tGATTTTTTTGTTTTTCAATATTATTCATTTATTCCATCAAAAACACCCTTTTTAATTTTTACAAAAATTAAAAAAAGTGCTCCTACACTg                                                              <  1:2778517/91‑1 (MQ=255)
                                tttttttGTTTTTCAATATTATTCATTTATTCCATCAAAAACACCCTTTTTAATTTTTACAAAAATTAAAAAAAGTGCTCCTACACTGCtt                                                           >  2:2025008/1‑91 (MQ=255)
                                     ttGTTTTTCAATATTATTCATTTATTCCATCAAAAACACCCTTTTTAATTTTTACAAAAATTAAAAAAAGTGCTCCTACACTGCTTGCa                                                        >  2:2469885/1‑89 (MQ=255)
                                     ttGTTTTTCAATATTATTCATTTATTCCATCAAAAACACCCTTTTTAATTTTTACAAAAATTAAAAAAAGTGCTCCTACACTGCTTGCa                                                        <  1:2469885/89‑1 (MQ=255)
                                       gTTTTTCAATATTATTCATTTATTCCATCAAAAACACCCTTTTTAATTTTTACAAAAATTAAAAAAAGTGCTCCTACACTGCTTGCATGTa                                                    >  2:296944/1‑91 (MQ=255)
                                       gTTTTTCAATATTATTCATTTATTCCATCAAAAACACCCTTTTTAATTTTTACAAAAATTAAAAAAAGTGCTCCTACACTGCTTGCATGTa                                                    >  1:609558/1‑91 (MQ=255)
                                       gTTTTTCAATATTATTCATTTATTCCATCAAAAACACCCTTTTTAATTTTTACAAAAATTAAAAAAAGTGCTCCTACACTGCTTGCATGTa                                                    >  1:3696103/1‑91 (MQ=255)
                                          tttCAATATTATTCATTTATTCCATCAAAAACACCCTTTTTAATTTTTACAAAAATTAAAAAAAGTGCTCCTACACTGCTTGCATGTAGaa                                                 >  1:1314292/1‑91 (MQ=255)
                                            tctatatTATTCATTTATTCCATCAAAAACACCCTTTTTAATTTTTACAAAAATTAAAAAAAGTGCTCCTACACTGCTTGCATGTAGAAac                                               >  2:2615298/4‑91 (MQ=255)
                                            tCAATCTTATTCATTTATTCCATCAAAAACACCCTTTTTAATTTTTACAAAAATTAAAAAAAGTGCTCCTACACTGCTTGCATGTAGAAac                                               <  2:3584535/91‑1 (MQ=255)
                                              aaTATTATTCATTTATTCCATCAAAAACACCCTTTTTAATTTTTACAAAAATTAAAAAAAGTGCTCCTACACTGCTTGCATGTAGAAacac                                             <  1:1331391/91‑1 (MQ=255)
                                               atatTATTCATTTATTCCATCAAAATCACCCTTTTTAATTTTTACAAAAATTAAAAAAAGTGCTCCTACACTGCTTGCATGTAGAAACACt                                            <  2:16925/91‑1 (MQ=255)
                                               atatTATTCATTTATTCCATCAAAAACACCCTTTTTAATTTTTACAAAAATTAAAAAAAGTGCTCCTACACTGCTTGCATGTAGAAACACt                                            <  2:263492/91‑1 (MQ=255)
                                                  ttattCATTTATTCCATCAAAAACACCCTTTTTAATTTTTTCAAAAATTAAAAAAAGTGCTCCTACACTGCTTGCATGTAGAAACACtttt                                         >  2:634179/1‑91 (MQ=255)
                                                  ttattCATTTATTCCATCAAAAACACCCTTTTTAATTTTTACAAAAATTAAAAAAAGTGCTCCTACACTGCTTGCATGTAGAAACACtttt                                         >  1:1472332/1‑91 (MQ=255)
                                                   tattCATTTATTCCATCAAAAACACCCTTTTTAATTTTTACAAAAATTAAAAAAAGTGCTCCTACACTGCTTGCATGTAGAAACACttttt                                        >  2:2350020/1‑91 (MQ=255)
                                                      tCATTTATTCCATCAAAAACACCCTTTTTAATTTTTACAAAAATTAAAAAAAGTGCTCCTACACTGCTTGCATGTAGAAACACTTTTTCAt                                     >  1:1543734/1‑91 (MQ=255)
                                                       cATTTATTCCATCAAAAACACCCTTTTTAATTTTTACAAAAATTAAAAAAAGTGCTCCTACACTGCTTGCATGTAGAAACACTTTTTCAtt                                    >  1:886303/1‑91 (MQ=255)
                                                        aTTTATTCCATCAAAAACACCCTTTTTAATTTTTACAAAAATTAAAAAAAGTGCTCCTACACTGCTTGCATGTAGAAACACTTTTTCATTg                                   >  1:2348198/1‑91 (MQ=255)
                                                         tttATTCCATCAAAAACACCCTTTTTAATTTTTACAAAAATTAAAAAAAGTGCTCCTACACTGCTTGCATGTAGAAACACTTTTTCATTGt                                  <  2:215589/91‑1 (MQ=255)
                                                          ttATTCCATCAAAAACACCCTTTTTAATTTTTACAAAAATTAAAAAAAGTGCTCCTACACTGCTTGCATGTAGAAACACTTTTTCATTGTa                                 <  2:2127308/91‑1 (MQ=255)
                                                           tATTCCATCAAAAACACCCTTTTTAATTTTTACAAAAATTAAAAAAAGTGCTCCTACACTGCTTGCATGTAGAAACACTTTTTCATTGTaa                                <  2:353865/91‑1 (MQ=255)
                                                             ttCCATCAAAAACACCCTTTTTAATTTTTACAAAAATTAAAAAAAGTGCTCCTACACTGCTTGCATGTAGAAACACTTTTTCATTGTAATg                              >  1:2275062/1‑91 (MQ=255)
                                                             ttCCATCAAAAACACCCTTTTTAATTTTTACAAAAATTAAAAAAAGTGCTCCTACACTGCTTACATGTAGAAACACTTTTTCATTGTAATg                              >  2:2848183/1‑91 (MQ=255)
                                                             ttCCATCAAAAACACCCTTTTTAATTTTGACAAAAATTAAAAAAAGTGCTCCTACACTGCTTGCATGTAGAAACACTTTTTCATTGTAATg                              <  1:1090010/91‑1 (MQ=255)
                                                              tCCATCAAAAACACCCTTTTTAATTTTTACAAAAATTAAAAAAAGTGCTCCTACACTGCTTGCATGTAGAAACACTTTTTCATTGTAATGt                             >  2:2250265/1‑91 (MQ=255)
                                                              tCCATCAAAAACACCCTTTTTAATTTTTACAAAAATTAAAAAAAGTGCTCCTACACTGCTTGCATGTAGAAACACTTTTTCATTGTAATGt                             <  1:478246/91‑1 (MQ=255)
                                                               ccATCAAAAACACCCTTTTTAATTTTTACAAAAATTAAAAAAAGTGCTCCTACACTGCTTGCATGTAGAAACACTTTTTCATTGTAATGtt                            <  1:2544283/91‑1 (MQ=255)
                                                                 aTCAAAAACACCCTTTTTAATTTTTACAAAAATTAAAAAAAGTGCTCCTACACTGCTTGCATGTAGAAACACTTTTTCATTGTAATGTTAt                          >  1:2452225/1‑91 (MQ=255)
                                                                  tCAAAAACACCCTTTTTAATTTTTACAAAAATTAAAAAAAGTGCTCCTACACTGCTTGCATGTAGAAACACTTTTTCATTGTAATGTTAtt                         <  2:2181504/91‑1 (MQ=255)
                                                                  tCAAAAACACCCTTTTTAATTTTTACAAAAATTAAAAAAAGTGCTCCTACACTGCTTGCATGTAGAAACACTTTTTCATTGTAATGTTAtt                         <  1:2517524/91‑1 (MQ=255)
                                                                   cAAAAACACCCTTTTTAATTTTTACAAAAATTAAAAAAAGTGCTCCTACACTGCTTGCATGTAGAAACACTTTTTCATTGTAATGTTAttc                        >  2:2111545/1‑91 (MQ=255)
                                                                   cAAAAACACCCTTTTTAATTTTTACAAAAATTAAAAAAAGTGCTCCTACACTGCTTGCATGTAGAAACACTTTTTCATTGTAATGTTAttc                        <  1:2563606/91‑1 (MQ=255)
                                                                    aaaaaCATCCTTTTTAATTTTTACAAAAATTAAAAAAAGTGCTCCTACACTGCTTGCATGTAGAAACACTTTTTCATTGTAATGTTAttct                       <  2:3448874/91‑1 (MQ=255)
                                                                       aaCACCCTTTTTAATTTTTACAAAAATTAAAAAAAGTGCTCCTACACTGCTTGCATGTAGAAACACTTTTTCATTGTAATGTTATtcttct                    >  1:335401/1‑91 (MQ=255)
                                                                        acacCCTTTTTAATTTTTACAAAAATTAAAAAAAGTGCTCCTACACTGCTTGCATGTAGAAACACTTTTTCATTGTAATGTTATTCTtctc                   >  1:362178/1‑91 (MQ=255)
                                                                         cacCCTTTTTAATTTTTACAAAAATTAAAAAAAGTGCTCCTACACTGCTTGCATGTAGAAACACTTTTTCATTGTAATGTTATTCTTCTCg                  >  1:1815588/1‑91 (MQ=255)
                                                                         cacCCTTTTTAATTTTTACAAAAATTAAAAAAAGTGCTCCTACACTGCTTGCATGTAGAAACACTTTTTCATTGTAATGTTATTCTTCTCg                  >  2:2289079/1‑91 (MQ=255)
                                                                         cacCCTTTTTAATTTTTACAAAAATTAAAAAAAGTGCTCCTACACTGCTTGCATGTAGAAACACTTTTTCATTGTAATGTTATTCTTCTCg                  >  1:724197/1‑91 (MQ=255)
                                                                           cccTTTTTAATTTTTACAAAAATTAAAAAAAGTGCTCCTACACTGCTTGCATGTAGAAACACTTTTTCATTGTAATGTTATTCTTCTCgag                >  1:3772522/1‑91 (MQ=255)
                                                                           cccTTTTTAATTTTTACAAAAATTAAAAAAAGTGCTCCTACACTGCTTGCATGTAGAAACACTTTTTCATTGTAATGTTATTCTTCTCgag                >  2:269545/1‑91 (MQ=255)
                                                                            ccTTTTTAATTTTTACAAAAATTAAAAAAAGTGCTCCTACACTGCTTGCATGTAGAAACACTTTTTCATTGTAATGTTATTCTTCTCgaga               <  1:1146524/91‑1 (MQ=255)
                                                                              tttttAATTTTTACAAAAATTAAAAAAAGTGCTCCTACACTGCTTGCATGTAGAAACACTTTTTCATTGTAATGTTATTCTTCTCGAGACa             <  1:3747619/91‑1 (MQ=255)
                                                                              tttttAATTTTTACAAAAATTAAAAAAAGTGCTCCTACACTGCTTGCATGTAGAAACACTTTTTCATTGTAATGTTATTCTTCTCGAGACa             >  1:3251804/1‑91 (MQ=255)
                                                                              tttttAATTTTTACAAAAATTAAAAAAAGTGCTCCTACACTGCTTGCATGTAGAAACACTTTTTCATTGTAATGTTATTCTTCTCGAGACa             >  1:2050428/1‑91 (MQ=255)
                                                                               ttttAATTTTTGTAAAAATTAAAAAGGGTGCTCCTACACTGCTTGCATGTAGAAACACTTTTTCATTGTAATGTTATTCTTCTCGAGACAt            >  1:2681920/1‑91 (MQ=37)
                                                                               ttttAATTTTTACAAAAATTAAAAAAAGTGCTCCTACACTGCTTGCATGTAGAAACACTTTTTCATTGTAATGTTATTCTTCTCGAGACAt            >  1:770519/1‑91 (MQ=255)
                                                                                tttGATTTTTACAAAAATTAAAAAAAGTGCTCCTACACTGCTTGCATGTAGAAACACTTTTTCATTGTAATGTTATTCTTCTCGAGACATa           <  1:3129796/91‑1 (MQ=255)
                                                                                  tAATTTTTGTAAAAATTAAAAAAAGTGCTCCTACACTGCTTGCATGTAGAAACACTTTTTCATTGTAATGTTATTCTTCTCGAGACATAcc         >  2:415651/1‑91 (MQ=255)
                                                                                  tAATTTTTGTAAAAATTAAAAAAAGTGCTCCCACACTGCTTGCATGTAGAAACACTTTTTCATTGTAATGTTATTCTTCTCGAGACATAcc         >  1:782544/1‑91 (MQ=255)
                                                                                     tttttACAAAAATTAAAAAAAGTGCTCCTACACTGCTTGCATGTAGAAACACTTTTTCATTGTAATGTTATTCTTCTCGAGACATACCttt      >  2:3819421/1‑91 (MQ=255)
                                                                                      ttttACAAAAATTAAAAAAAGTGCTCCTACACTGCTTGCATGTAGAAACACTTTTTCATTGTAATGTTATTCTTCTCGAGACATACCtttt     <  2:1198927/91‑1 (MQ=255)
                                                                                       tttACAAAAATTAAAAAAAGTGCTCCTACACTGCTTGCATGTAGAAACACTTTTTCATTGTAATGTTATTCTTCTCGAGACATACCTTTTa    <  2:752041/91‑1 (MQ=255)
                                                                                       tttACAAAAATTAAAAAAAGTGCTCCTACACTGCTTGCATGTAGAAACACTTTTTCATTGTAATGTTATTCTTCTCGAGACATACCTTTTa    >  2:3035369/1‑91 (MQ=255)
                                                                                        ttACAAAAATTAAAAAAAGTGCTCCTACACTGCTTGCATGTAGAAACACTTTTTCATTGTAATGTTATTCTTCTCGAGACATACCTTTTAg   <  1:2933922/91‑1 (MQ=255)
                                                                                         tACAAAAATTAAAAAAAGTGCTCCTACACTGCTTGCATGTAGAAACACTTTTTCATTGTAATGTTATTCTTCTCGAGACATACCTTTTAGc  <  1:3688243/91‑1 (MQ=255)
                                                                                          |                                                                                         
CCTGTCGATGATTGATTCACAATAGCACGTGATTTTTTTGTTTTTCAATATTATTCATTTATTCCATCAAAAACACCCTTTTTAATTTTTACAAAAATTAAAAAAAGTGCTCCTACACTGCTTGCATGTAGAAACACTTTTTCATTGTAATGTTATTCTTCTCGAGACATACCTTTTAGC  >  USA300TCH1516_ALE/2424518‑2424697

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: