Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 2,424,609 C→T 6.8% intergenic (‑167/+61) 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,6090CT 6.8% 748.8 / 24.6 248intergenic (‑167/+61)USA300HOU_RS12490/USA300HOU_RS12495hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (118/110);  new base T (14/4);  total (133/115)
Fisher's exact test for biased strand distribution p-value = 4.75e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 250 total aligned reads displayed.

CTGTCGATGATTGATTCACAATAGCACGTGATTTTTTTGTTTTTCAATATTATTCATTTATTCCATCAAAAACACCCTTTTTAATTTTTACAAAAATTAAAAAAAGTGCTCCTACACTGCTTGCATGTAGAAACACTTTTTCATTGTAATGTTATTCTTCTCGAGACATACCTTTTAGC  >  USA300TCH1516_ALE/2424519‑2424697
                                                                                          |                                                                                        
 tGTCGATGATTGATTCACAATAGCACGTGATTTTTTTGTTTTTCAATATTATTCATTTATTCCATCAAAAACACCCTTTTTAATTTTTACa                                                                                         <  1:3068351/91‑1 (MQ=255)
 tGTCGATGATTGATTCACAATAGCACGTGATTTTTTTGTTTTTCAATATTATTCATTTATTCCATCAAAAACACCCTTTTTAATTTTTACa                                                                                         <  1:3068605/91‑1 (MQ=255)
 tGTCGATGATTGATTCACAATAGCACGTGATTTTTTTGTTTTTCAATATTATTCATTTATTCCATCAAAAACACCCTTTTTAATTTTTACa                                                                                         <  2:586250/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                                                                                     <  1:2863172/91‑1 (MQ=255)
     gatgatTGATTCACAATAGCACGTGATTTTTTTGTTTTTCAATATTATTCATTTATTCCATCAAAAACACCCTTTTTAATTTTTACaaaaa                                                                                     <  1:192428/91‑1 (MQ=255)
     gatgatTGATTCACAATAGCACGTGATTTTTTTGTTTTTCAATATTATTCATTTATTCCATCAAAAACACCCTTTTTAATTTTTACaaaaa                                                                                     >  1:2168152/1‑91 (MQ=255)
     gatgatTGATTCACAATAGCACGTGATTTTTTTGTTTTTCAATATTATTCATTTATTCCATCAAAAACACCCTTTTTAATTTTTACaaaaa                                                                                     >  1:2169468/1‑91 (MQ=255)
     aatgatTGATTCACAATAGCACGTGATTTTTTTGTTTTTCAATATTATTCATTTATTCCATCAAAAACACCCTTTTTAATTTTTACaaaaa                                                                                     <  2:1572319/90‑1 (MQ=255)
       tgattgatTCACAATAGCACGTGATTTTTTTGTTTTTCAATATTATTCATTTATTCCATCAAAAACACCCTTTTTAATTTTTTCAAAAAtt                                                                                   >  2:468935/1‑91 (MQ=255)
       tgattgatTCACAATAGCACGTGATTTTTTTGTTTTTCAATATTATTCATTTATTCCATCAAAAACACCCTTTTTAATTTTTACAAAAAtt                                                                                   >  2:1623805/1‑91 (MQ=255)
       tgattgatTCACAATAGCACGTGATTTTTTTGTTTTTCAATATTATTCATTTATTCCATCAAAAACACCCTTTTTAATTTTTACAAAAAtt                                                                                   >  1:3001366/1‑91 (MQ=255)
       tgattgatTCACAATAGCACGTGATTTTTTTGTTTTTCAATATTATTCATTTATTCCATCAAAAACACCCTTTTTAATTTTTACAAAAAtt                                                                                   >  1:2973988/1‑91 (MQ=255)
        gattgattCACAATAGCACGTGATTTTTTTGTTTTTCAATATTATTCATTTATTCCATCAAAAACACCCTTTTTAATTTTTGTAAAAATTa                                                                                  <  1:3149678/91‑1 (MQ=255)
        gattgattCACAATAGCACGTGATTTTTTTGTTTTTCAATATTATTCATTTATTCCATCAAAAACACCCTTTTTAATTTTTACAAAAATTa                                                                                  >  1:1768937/1‑91 (MQ=255)
            gattCACAATAGCACGTGATTTTTTTGTTTTTCAATATTATTCATTTATTCCATCAAAAACACCCTTTTTAATTTTTACAAAAATTaaaaa                                                                              <  2:1850713/91‑1 (MQ=255)
              ttCACAATAGCACGTGATTTTTTTGTTTTTCAATATTATTCATTTATTCCATCAAAAACACCCTTTTTAATTTTTACAAAAATTaaaaaaa                                                                            >  2:3045415/1‑91 (MQ=255)
              ttCACAATAGCACGTGATTTTTTTGTTTTTCAATATTATTCATTTATTCCATCAAAAACACCCTTTTTAATTTTTACAAAAATTaaaaaaa                                                                            <  1:1865425/91‑1 (MQ=255)
              ttCACAATAGCACGTGATTTTTTTGTTTTTCAATATTATTCATTTATTCCATCAAAAACACCCTTTTTAATTTTTACAAAAATTaaaaaaa                                                                            <  1:3670248/91‑1 (MQ=255)
                cacaATAGCACGTGATTTTTTTGTTTTTCAATATTATTCATTTATTCCATCAAAAACACCCTTTTTAATTTTTACAAAAATTAAAAAAAGt                                                                          >  1:3749778/1‑91 (MQ=255)
                cacaATAGCACGTGATTGTTCTGTTCTTCAATATTATTCATTTATTCCATCAAAAACACCCTTTTTAATTTTTGTAAAAATTaaaaagggt                                                                          <  2:834228/91‑5 (MQ=37)
                 acaATAGCACGTGATTTTTTTGTTTTTCAATATTATTCATTTATTCCATCAAAAACACCCTTTTTAATTTTTACAAAAATTAAAAAAAGTg                                                                         <  2:1174973/91‑1 (MQ=255)
                  caATAGCACGTGATTTTTTTGTTTTTCAATATTATTCATTTATTCCATCAAAAACACCCTTTTTAATTTTTACAAAAATTAAAAAAAGTGc                                                                        >  1:3029238/1‑91 (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                                                                    >  1:323440/1‑91 (MQ=255)
                       gCACGTGATTTTTTTGTTTTTCAATATTATTCATTTATTCCATCAAAAACACCCTTTTTAATTTTTACAAAAATTAAAAAAAGTGCTCCTa                                                                   <  2:3405180/91‑1 (MQ=255)
                        cACGTGATTTTTTTGTTTTTCAATATTATTCATTTATTCCATCAAAAACACCCTTTTTAATTTTTACAAAAATTAAAAAAAGTGCTCCTac                                                                  >  1:2600711/1‑91 (MQ=255)
                         aCGTGATTTTTTTGTTTTTCAATATTATTCATTTATTCCATCAAAAACACCCTTTTTAATTTTTACAAAAATTAAAAAAAGTGCTCCTaca                                                                 <  1:980452/91‑1 (MQ=255)
                          cGTGATTTTTTTGTTTTTCAATATTATTCATTTATTCCATCAAAAACACCCTTTTTAATTTTTACAAAAATTAAAAAAAGTGCTcc                                                                     <  1:1982418/86‑1 (MQ=255)
                          cGTGATTTTTTTGTTTTTCAATATTATTCATTTATTCCATCAAAAACACCCTTTTTAATTTTTACAAAAATTAAAAAAAGTGCTCCTacac                                                                >  2:1777915/1‑91 (MQ=255)
                          cGTGATTTTTTTGTTTTTCAATATTATTCATTTATTCCATCAAAAACACCCTTTTTAATTTTTACAAAAATTAAAAAAAGTGCTCCTacac                                                                >  1:654272/1‑91 (MQ=255)
                          cGTGATTTTTTTGTTTTTCAATATTATTCATTTATTCCATCAAAAACACCCTTTTTAATTTTTACAAAAATTAAAAAAAGTGCTCCTacac                                                                <  1:3228750/91‑1 (MQ=255)
                          cGTGATTTTTTTGTTTTTCAATATTATTCATTTATTCCATCAAAAACACCCTTTTTAATTTTTACAAAAATTAAAAAAAGTGCTCCTacac                                                                >  2:2007865/1‑91 (MQ=255)
                          cGTGATTTTTTTGTTTTTCAATATTATTCATTTATTCCATCAAAAACACCCTTTTTAATTTTTACAAAAATAAAAAAAAGTGCTCCTacac                                                                >  2:3165188/1‑91 (MQ=255)
                              aTTTTTTTGTTTTTCAATATTATTCATTTATTCCATCAAAAACACCCTTTTTAATTTTTACAAAAATTAAAAAAAGTGCTCCTACACTGCt                                                            <  2:654314/91‑1 (MQ=255)
                               tttttttGTTTTTTTATATTATTCATTTATTCCATCAAAAACACCCTTTTTAATTTTTACAAAAATTAAAAAAAGTGCTCCTACACTGCtt                                                           >  2:821741/1‑91 (MQ=255)
                               tttttttGTTTTTCAATATTATTCATTTATTCCATCAAAAACACCCTTTTTAATTTTTACAAAAATTAAAAAAAGTGCTCCTACACTGCtt                                                           >  2:2025008/1‑91 (MQ=255)
                                  ttttGTTTTTCAATATTATTCATTTATTCCATCAAAAACACCCTTTTTAATTTTTACAAAAATTAAAAAAAGTGCTCCTACACTGCTTGCa                                                        >  1:394072/1‑91 (MQ=255)
                                    ttGTTTTTCAATATTATTCATTTATTCCATCAAAAACACCCTTTTTAATTTTTACAAAAATTAAAAAAAGTGCTCCTACACTGCTTGCa                                                        <  1:2469885/89‑1 (MQ=255)
                                      gTTTTTCAATATTATTCATTTATTCCATCAAAAACACCCTTTTTAATTTTTACAAAAATTAAAAAAAGTGCTCCTACACTGCTTGCATGTa                                                    >  1:609558/1‑91 (MQ=255)
                                      gTTTTTCAATATTATTCATTTATTCCATCAAAAACACCCTTTTTAATTTTTACAAAAATTAAAAAAAGTGCTCCTACACTGCTTGCATGTa                                                    >  1:3696103/1‑91 (MQ=255)
                                       tttttCAATATTATTCATTTATTCCATCAAAAACACCCTTTTTAATTTTTACAAAAATTAAAAAAAGTGCTCCTACACTGCTTGCATGTAg                                                   >  2:2261753/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)
                                             aaTATTATTCATTTATTCCATCAAAAACACCCTTTTTAATTTTTACAAAAATTAAAAAAAGTGCTCCTACACTGCTTGCATGTAGAAacac                                             <  1:359161/91‑1 (MQ=255)
                                              atatTATTCATTTATTCCATCAAAATCACCCTTTTTAATTTTTACAAAAATTAAAAAAAGTGCTCCTACACTGCTTGCATGTAGAAACACt                                            <  2:16925/91‑1 (MQ=255)
                                              atatTATTCATTTATTCCATCAAAAACACCCTTTTTAATTTTTACAAAAATTAAAAAAAGTGCTCCTACACTGCTTGCATGTAGAAACACt                                            <  2:263434/91‑1 (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                         <  1:2517524/91‑1 (MQ=255)
                                                                 tCAAAAACACCCTTTTTAATTTTTACAAAAATTAAAAAAAGTGCTCCTACACTGCTTGCATGTAGAAACACTTTTTCATTGTAATGTTAtt                         <  2:2181504/91‑1 (MQ=255)
                                                                  cAAAAACACCCTTTTTAATTTTTACAAAAATTAAAAAAAGTGCTCCTACACTGCTTGCATGTAGAAACACTTTTTCATTGTAATGTTAttc                        <  1:2563606/91‑1 (MQ=255)
                                                                  cAAAAACACCCTTTTTAATTTTTACAAAAATTAAAAAAAGTGCTCCTACACTGCTTGCATGTAGAAACACTTTTTCATTGTAATGTTAttc                        >  2:2111545/1‑91 (MQ=255)
                                                                        cacCCTTTTTAATTTTTACAAAAATTAAAAAAAGTGCTCCTACACTGCTTGCATGTAGAAACACTTTTTCATTGTAATGTTATTCTTCTCg                  >  1:1815588/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)
                                                                             tttttAATTTTTACAAAAATTAAAAAAAGCGCTCCTACACGACTTGCCTGTAGAAACACTTTTTCATTGTAATGTTATTCTTATCGAGACa             >  2:2739346/1‑91 (MQ=255)
                                                                              ttttAATTTTTGTAAAAATTAAAAAGGGTGCTCCTACACTGCTTGCATGTAGAAACACTTTTTCATTGTAATGTTATTCTTCTCGAGACAt            >  1:2681920/1‑91 (MQ=37)
                                                                              ttttAATTTTTACAAAAATTAAAAAAAGTGCTCCTACACTGCTTGCATGTAGAAACACTTTTTCATTGTAATGTTATTCTTCTCGAGACAt            <  1:2863731/91‑1 (MQ=255)
                                                                              ttttAATTTTTACAAAAATTAAAAAAAGTGCTCCTACACTGCTTGCATGTAGAAACACTTTTTCATTGTAATGTTATTCTTCTCGAGACAt            >  1:770519/1‑91 (MQ=255)
                                                                               tttGATTTTTACAAAAATTAAAAAAAGTGCTCCTACACTGCTTGCATGTAGAAACACTTTTTCATTGTAATGTTATTCTTCTCGAGACATa           <  1:3129796/91‑1 (MQ=255)
                                                                               tttAATTTTTGTAAAAATTAAAAAGGGTGCTCCTACACTGCTTGCATGTAGAAACACTTTTTCATTGTAATGTTATTCTTCTCGAGACATa           >  1:2709240/1‑91 (MQ=38)
                                                                                ttAATTTTTGTAAAAATTAAAAAAAGTGCTCCTACACTGCTTGCATGTAGAAACACTTTTTCATTGTAATGTTATTCTTCTCGAGACATAc          >  2:2445342/1‑91 (MQ=255)
                                                                                 tAATTTTTGTAAAAATTAAAAAAAGTGCTCCTACACTGCTTGCATGTATATACACTTTTTCATTGTAATGTTATTCTTCTCGAGACATAcc         >  2:3555443/1‑91 (MQ=255)
                                                                                 tAATTTTTGTAAAAATTAAAAAAAGTGCTCCTACACTGCTTGCATGTAGAAACACTTTTTCATTGTAATGTTATTCTTCTCGAGACATAcc         >  1:3686200/1‑91 (MQ=255)
                                                                                 tAATTTTTGTAAAAATTAAAAAAAGTGCTCCCACACTGCTTGCATGTAGAAACACTTTTTCATTGTAATGTTATTCTTCTCGAGACATAcc         >  1:782544/1‑91 (MQ=255)
                                                                                 tAATTTTTACAAAAATTAAAAAAAGTGCTCCTACACTGCTTGCATGTAGAAACACTTTTTCATTGTAATGTTATTCTTCTCGAGACATAcc         >  1:335044/1‑91 (MQ=255)
                                                                                  aaTTTTTACAAAAAATAAAAAAAGTGCTCCTACACTGCTTGCATGTAGAAACACTTTTTCATTGTAATGTTATTCTTCTCGAGACATACCt        <  2:2876405/91‑1 (MQ=255)
                                                                                   aTTTTTATAAAAATTAAAAAAAGTGCTCCTACACTGCTTGCATGTAGAAACACTTTTTCATTGTAATGTTATTCTTCTCGAGACATAGCtt       <  2:3438683/91‑1 (MQ=255)
                                                                                    tttttACAAAAATTAAAAAAAGTGCTCCTACACTGCTTGCATGTAGAAACACTTTTTCATTGTAATGTTATTCTTCTCGAGACATACCttt      >  2:3819421/1‑91 (MQ=255)
                                                                                     ttttACAAAAATTAAAAAAAGTGCTCCTACACTGCTTGCATGTAGAAACACTTTTTCATTGTAATGTTATTCTTCTCGAGACATACCtttt     <  2:1198927/91‑1 (MQ=255)
                                                                                      tttACAAAAATTAAAAAAAGTGCTCCTACACTGCTTGCATGTAGAAACACTTTTTCATTGTAATGTTATTCTTCTCGAGACATACCTTTTa    >  2:3035369/1‑91 (MQ=255)
                                                                                      tttACAAAAATTAAAAAAAGTGCTCCTACACTGCTTGCATGTAGAAACACTTTTTCATTGTAATGTTATTCTTCTCGAGACATACCTTTTa    <  2:752041/91‑1 (MQ=255)
                                                                                       ttACAAAAATTAAAAAAAGTGCTCCTACACTGCTTGCATGTAGAAACACTTTTTCATTGTAATGTTATTCTTCTCGAGACATACCTTTTAg   <  1:2933922/91‑1 (MQ=255)
                                                                                          |                                                                                        
CTGTCGATGATTGATTCACAATAGCACGTGATTTTTTTGTTTTTCAATATTATTCATTTATTCCATCAAAAACACCCTTTTTAATTTTTACAAAAATTAAAAAAAGTGCTCCTACACTGCTTGCATGTAGAAACACTTTTTCATTGTAATGTTATTCTTCTCGAGACATACCTTTTAGC  >  USA300TCH1516_ALE/2424519‑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: