Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 2,608,820 A→T 6.1% intergenic (+101/+159) USA300TCH1516_02485 → / ← USA300HOU_RS13450 hypothetical protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE2,608,8200AT 6.1% 754.3 / 26.8 227intergenic (+101/+159)USA300TCH1516_02485/USA300HOU_RS13450hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (116/97);  new base T (6/8);  total (122/105)
Fisher's exact test for biased strand distribution p-value = 4.22e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.99e-01

Only 100 of 231 total aligned reads displayed.

TTTATATATAGCAGACATGATAGAATTTTATATGTAAATCTTGTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGTTTACTTAATTAAAGGGGTCTGCCCTCAATTGCTCACCAAC  >  USA300TCH1516_ALE/2608730‑2608910
                                                                                          |                                                                                          
tttATATATAGCAGACATGATAGAATTTTATATGTAAATCTTGTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCa                                                                                            >  2:1158521/1‑91 (MQ=255)
    tatataGCAGACATGATAGAATTTTATATGTAAATCTTGTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAatat                                                                                        <  2:742445/91‑1 (MQ=255)
    tatataGCAGACATGATAGAATTTTATATGTAAATCTTGTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAatat                                                                                        >  1:1123432/1‑91 (MQ=255)
    tatataGCAGACATGATAGAATTTTATATGTAAATCTTGTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAatat                                                                                        >  1:3754818/1‑91 (MQ=255)
    tatataGCAGACATGATAGAATTTTATATGTAAATCTTGTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAatat                                                                                        >  2:3623500/1‑91 (MQ=255)
    tatataGCAGACATGATAGAATTTTATATGTAAATCTTGTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAatat                                                                                        <  2:3551966/91‑1 (MQ=255)
     atataGCAGACATGATAGAATTTTATATGTAAATCTTGTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAAtata                                                                                       >  1:973720/1‑91 (MQ=255)
      tataGCAGACATGATAGAATTTTATATGTAAATCTTGTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAg                                                                                      <  1:2095884/91‑1 (MQ=255)
       ataGCAGACATGATAGAATTTTATATGTAAATCTTGTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACGTCCTatattgg                                                                                     <  2:1279975/91‑4 (MQ=255)
       ataGCAGACATGATAGAATTTTATATGTAAATCTTGTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACGTCCTatattgg                                                                                     <  2:1785499/91‑4 (MQ=255)
        taGCAGACATGATAGAATTTTATATGTAAATCTTGTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGa                                                                                    >  1:3250617/1‑91 (MQ=255)
        taGCAGACATGATAGAATTTTATATGTAAATCTTGTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGa                                                                                    <  2:2981104/91‑1 (MQ=255)
        taGCAGACATGATAGAATTTTATATGTAAATCTTGTAGGTAATCGTTCTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGa                                                                                    >  1:3872346/1‑91 (MQ=255)
         aGCAGACATGATAGAATTTTATATGTAAATCTTGTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACGTCCTATATTGGat                                                                                   <  1:3684316/91‑2 (MQ=255)
         aGCAGACATGATAGAATTTTATATGTAAATCTTGTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGAc                                                                                   <  2:1362829/91‑1 (MQ=255)
          gCAGACATGATAGAATTTTATATGTAAATCTTGTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACg                                                                                  <  2:3179743/91‑1 (MQ=255)
          gCAGACATGATAGAATTTTATATGTAAATCTTGTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACg                                                                                  >  1:60708/1‑91 (MQ=255)
          gCAGACATGATAGAATTTTATATGTAAATCTTGTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACg                                                                                  <  1:1654083/91‑1 (MQ=255)
          gCAGACATGATAGAATTTTATATGTAAATCTTGTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACg                                                                                  <  1:2818040/91‑1 (MQ=255)
           cAGACATGATAGAATTTTATATGTAAATCTTGTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGt                                                                                 <  1:3155626/91‑1 (MQ=255)
            aGACATGATAGAATTTTATATGTAAATCTTGTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATCGGACGTg                                                                                >  1:1939583/1‑91 (MQ=255)
             gACATGATAGAATTTTATATGTAAATCTTGTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGTgg                                                                               >  1:3455715/1‑91 (MQ=255)
             gACATGATAGAATTTTATATGTAAAGCTTGTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGTgg                                                                               <  2:892480/91‑1 (MQ=255)
                  gATAGAATTTTATATGTAAATCTTGTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGTGGCtttt                                                                          <  2:2614859/91‑1 (MQ=255)
                  gATAGAATTTTATATGTAAATCTTGTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGTGGCtttt                                                                          >  1:1450169/1‑91 (MQ=255)
                   aTAGAATTTTATATGTAAATCTTGTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGTGGCTTTTa                                                                         <  2:3255877/91‑1 (MQ=255)
                    tAGAATTTTATATGTAAATCTTGTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGTGGCTTTTAg                                                                        <  2:658361/91‑1 (MQ=255)
                       aaTTTTATATGTAAATCTTGTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTAt                                                                     >  2:2122330/1‑91 (MQ=255)
                        aTTTTATATGTAAATCTTGTAGGTAATCGTTTTAAAAATAATATATGTATGACTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATc                                                                    <  2:2405273/91‑1 (MQ=255)
                        aTTTTATATGTAAATCTTGTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATc                                                                    <  2:2498598/91‑1 (MQ=255)
                           ttATATGTAAATCTTGTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGtt                                                                 >  2:2450879/1‑91 (MQ=255)
                            tataTGTAAATCTTGTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTc                                                                >  1:1881106/1‑91 (MQ=255)
                            tataTGTAAATCTTGTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTc                                                                <  2:3423876/91‑1 (MQ=255)
                               atGTAAATCTTGTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTCCAg                                                             <  1:2759964/91‑1 (MQ=255)
                               atGTAAATCTTGTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTCCAg                                                             >  1:640083/1‑91 (MQ=255)
                               atGTAAATCTTGTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTCCAg                                                             <  2:3810529/91‑1 (MQ=255)
                                   aaaTCTTGTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGAttt                                                         <  2:105057/91‑1 (MQ=255)
                                   aaaTCTTGTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGAttt                                                         <  2:2193865/91‑1 (MQ=255)
                                   aaaTCTTGTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGAttt                                                         <  2:1123421/91‑1 (MQ=255)
                                    aaTCTTGTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGATTtg                                                        <  1:821873/91‑1 (MQ=255)
                                      tCTTGTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGATTtgtg                                                      >  2:1471333/1‑91 (MQ=255)
                                         tGTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGATTTGTgggg                                                   <  2:1655337/91‑1 (MQ=255)
                                         tGTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGATTTGTgggg                                                   >  1:1035918/1‑91 (MQ=255)
                                         tGTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGATTTGTgggg                                                   <  2:442633/91‑1 (MQ=255)
                                          gTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGa                                                  >  1:2028591/1‑91 (MQ=255)
                                           tAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGTGGCtttt                                                                          <  2:3796336/66‑1 (MQ=255)
                                            aGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGt                                                >  2:2560002/1‑91 (MQ=255)
                                            aGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGt                                                <  2:1553090/91‑1 (MQ=255)
                                             ggTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGtt                                               <  2:1245783/91‑1 (MQ=255)
                                                aaTCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAgg                                            <  1:1373913/91‑1 (MQ=255)
                                                 aTCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAggg                                           <  1:1070588/91‑1 (MQ=255)
                                                  tCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGt                                          <  2:3511436/91‑1 (MQ=255)
                                                  tCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGt                                          <  2:1766724/91‑1 (MQ=255)
                                                  tCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGt                                          >  1:2932831/1‑91 (MQ=255)
                                                    gTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGttt                                        >  2:268843/1‑91 (MQ=255)
                                                    gTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGttt                                        >  2:231654/1‑91 (MQ=255)
                                                      ttttAAAATAATATAAGCATGACTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTATGGTTTAc                                      >  2:3436056/1‑91 (MQ=255)
                                                       ttAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGTTTACt                                     <  1:3404568/91‑1 (MQ=255)
                                                        tAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGTTTACtt                                    >  2:3720055/1‑91 (MQ=255)
                                                          aaaaTAATATAAGTATGACTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGTTTACttaa                                  >  2:1631992/1‑91 (MQ=255)
                                                          aaaaTAATATAAGTATGACTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGTTTACttaa                                  >  2:3871432/1‑91 (MQ=255)
                                                          aaaaTAATATAAGCATGACTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGTTTACttaa                                  >  2:2365606/1‑91 (MQ=255)
                                                           aaaTAATATAAGTATGACTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGTTTACttaat                                 <  2:3713542/91‑1 (MQ=255)
                                                              taataTAAGTATGACTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGTTTACttaattaa                              <  1:1442135/91‑1 (MQ=255)
                                                              taataTAAGTATGACTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGTTTACttaattaa                              >  2:141552/1‑91 (MQ=255)
                                                              taataTAAGTATGACTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGTTTACttaattaa                              >  2:702824/1‑91 (MQ=255)
                                                                atatAAGTATGACTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGTTTACTTAATTAAAg                            >  1:1795907/1‑91 (MQ=255)
                                                                 tataAGTATGACTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGTTTACTTAATTAAAgg                           >  1:1014982/1‑91 (MQ=255)
                                                                  ataAGTATGACTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGTTTACTTAATTAAAggg                          >  1:2575014/1‑91 (MQ=255)
                                                                  ataAGTATGACTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGTTTACTTAATTAAAggg                          >  2:2803351/1‑91 (MQ=255)
                                                                   taAGTATGACTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGTTTACTTAATTAAAgggg                         >  1:515476/1‑91 (MQ=255)
                                                                   taAGTATGACTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGAATTGAGGGGAGCTAGGGTTTACTTAATTAAAgggg                         >  2:2031465/1‑91 (MQ=255)
                                                                    aaGTATGACTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGTTTACTTAATTAAAGGGGt                        >  1:3846568/1‑91 (MQ=255)
                                                                     aGTATGACTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGACTTGTGGGGAGTTAGGGTTTACTTAATTAAAggggcc                       >  1:442970/1‑89 (MQ=255)
                                                                      gTATGTCTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGTTTACTTAATTAAAGGGGTCt                      <  2:1385046/91‑1 (MQ=255)
                                                                       tATGACTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGTTTACTTAATTAAAGGGGTCTg                     >  2:545344/1‑91 (MQ=255)
                                                                       tATGACTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGTTTACTTAATTAAAGGGGTCTg                     <  1:3204270/91‑1 (MQ=255)
                                                                         tGACTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGTTTACTTAATTAAAGGGGTCTGcc                   <  1:2231050/91‑1 (MQ=255)
                                                                         tGACTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGTTTACTTAATTAAAGGGGTCTGcc                   >  1:1752850/1‑91 (MQ=255)
                                                                          gACTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGTTTACTTAATTAAAGGGGTCTGccc                  <  1:1652376/91‑1 (MQ=255)
                                                                           aCTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGTTTACTTAATTAAAGGGGTCTGCCCt                 >  1:1730815/1‑91 (MQ=255)
                                                                            cTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGTTTACttaatta                               <  2:1838468/76‑1 (MQ=255)
                                                                             tAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGTTTACTTAATTAAAGGGGTCTGTCCTCa               >  2:3580776/1‑91 (MQ=255)
                                                                             tAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGTTTACTTAATTAAAGGGGTCTGCCCTCa               <  1:1890267/91‑1 (MQ=255)
                                                                             tAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGTTTACTTAATTAAAGGGGTCTGCCCTCa               <  1:78952/91‑1 (MQ=255)
                                                                             tAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGTTTACTTAATTAAAGGGGTCTGCCCTCa               >  1:2350315/1‑91 (MQ=255)
                                                                             aaaaGCCACGTCCTATATTGGATGTGGCTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGTTTACTTAATTAAAGGGGTCTGCCCTCaa              >  2:208011/2‑91 (MQ=25)
                                                                                 gCCACGTCCTATATTGGATGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGTTTACTTAATTAAAGGGGTCTGCCCTCAATTg           >  2:2028791/1‑91 (MQ=38)
                                                                                 gCCACGTCCTATATTGGATGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGTTTACTTAATTAAAGGGGTCTGCCCTCAATTg           >  1:3042426/1‑91 (MQ=38)
                                                                                 gCCACGTCCTATATTGGATGTGGCTTTTAGTATCGTTCCAGATCTGTGGGGAGTTAGGGTTTACTTAATTAAAGGGGTCTGCCCTCAATTg           >  2:2833859/1‑91 (MQ=38)
                                                                                 gCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGTTTACTTAATTAAAGGGGTCTGCCCTCAATTg           >  1:35404/1‑91 (MQ=255)
                                                                                 gCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGTTTACTTAATTAAAGGGGTCTGCCCTCAATTg           >  2:1075326/1‑91 (MQ=255)
                                                                                 gCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGTTTACTTAATTAAAGGGGTCTGCCCTCAATTg           <  2:794201/91‑1 (MQ=255)
                                                                                    acaTCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGTTTACTTAATTAAAGGGGTCTGCCCTCAATTGCTc        >  1:1612388/1‑91 (MQ=255)
                                                                                    acaTCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGTTTACTTAATTAAAGGGGTCTGCCCTCAATTGCTc        >  1:585088/1‑91 (MQ=255)
                                                                                     caTCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGTTTACTTAATTAAAGGGGTCTGCCCTCAATTGCTCa       >  1:2184853/1‑91 (MQ=255)
                                                                                      ctCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGTTTACTTAATTAAAGGGGTCTGCCCTCAATTGCTCAc      <  2:1734156/90‑1 (MQ=255)
                                                                                      aTCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGTTTACTTAATTAAAGGGGTCTGCCCTCAATTGCTCAc      <  1:2632043/91‑1 (MQ=255)
                                                                                         cAATATAGGACGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGTTTACTTAATTAAAGGGGTCTGCCCTCGATTGCTCACCaa   >  2:3223575/1‑91 (MQ=255)
                                                                                         cAATATAGGACGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGTTTACTTAATTAAAGGGGTCTGCCCTCAATTGCTCACCaa   >  1:712133/1‑91 (MQ=255)
                                                                                          |                                                                                          
TTTATATATAGCAGACATGATAGAATTTTATATGTAAATCTTGTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGTTTACTTAATTAAAGGGGTCTGCCCTCAATTGCTCACCAAC  >  USA300TCH1516_ALE/2608730‑2608910

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