Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 2,608,825 A→T 10.1% intergenic (+106/+154) 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,8250AT10.1% 701.5 / 56.3 230intergenic (+106/+154)USA300TCH1516_02485/USA300HOU_RS13450hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (109/96);  new base T (18/5);  total (128/102)
Fisher's exact test for biased strand distribution p-value = 2.62e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.10e-01

Only 100 of 244 total aligned reads displayed.

ATATAGCAGACATGATAGAATTTTATATGTAAATCTTGTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGTTTACTTAATTAAAGGGGTCTGCCCTCAATTGCTCACCAACGAGGG  >  USA300TCH1516_ALE/2608735‑2608915
                                                                                          |                                                                                          
atataGCAGACATGATAGAATTTTATATGTAAATCTTGTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAAtata                                                                                            >  1:973733/1‑91 (MQ=255)
   taGCAGACATGATAGAATTTTATATGTAAATCTTGTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGa                                                                                         >  1:2580664/1‑91 (MQ=255)
   taGCAGACATGATAGAATTTTATATGTAAATCTTGTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGa                                                                                         <  2:2981104/91‑1 (MQ=255)
   taGCAGACATGATAGAATTTTATATGTAAATCTTGTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGa                                                                                         <  2:3430420/91‑1 (MQ=255)
   taGCAGACATGATAGAATTTTATATGTAAATCTTGTAGGTAATCGTTCTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGa                                                                                         >  1:3872168/1‑91 (MQ=255)
    aGCAGACATGATAGAATTTTATATGTAAATCTTGTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACGTCCTATATTGGat                                                                                        <  1:3684213/91‑2 (MQ=255)
    aGCAGACATGATAGAATTTTATATGTAAATCTTGTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGAc                                                                                        <  2:1362829/91‑1 (MQ=255)
     gCAGACATGATAGAATTTTATATGTAAATCTTGTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACGTCCTATATTGGatg                                                                                       <  1:1371842/91‑3 (MQ=255)
     gCAGACATGATAGAATTTTATATGTAAATCTTGTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACg                                                                                       >  1:60708/1‑91 (MQ=255)
     gCAGACATGATAGAATTTTATATGTAAATCTTGTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACg                                                                                       <  1:2818040/91‑1 (MQ=255)
     gCAGACATGATAGAATTTTATATGTAAATCTTGTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACg                                                                                       <  2:3179743/91‑1 (MQ=255)
     gCAGACATGATAGAATTTTATATGTAAATCTTGTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACg                                                                                       <  1:1653847/91‑1 (MQ=255)
     gCAGACATGAAAGAATTTTATATGTAAATCTTGTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACg                                                                                       <  2:1581968/91‑1 (MQ=255)
      cAGACATGATAGAATTTTATATGTAAATCTTGTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGt                                                                                      <  1:3155626/91‑1 (MQ=255)
       aGACATGATAGAATTTTATATGTAAATCTTGTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGTg                                                                                     <  1:814023/91‑1 (MQ=255)
        gACATGATAGAATTTTATATGTAAATCTTGTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGTgg                                                                                    <  2:3818948/91‑1 (MQ=255)
        gACATGATAGAATTTTATATGTAAATCTTGTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGTgg                                                                                    >  1:3455715/1‑91 (MQ=255)
            tGATAGAATTTTATATGTAAATCTTGTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGTGGCttt                                                                                >  2:3765899/1‑91 (MQ=255)
             gATAGAATTTTATATGTAAATCTTGTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGTGGCtttt                                                                               >  2:1281453/1‑91 (MQ=255)
              aTAGAATTTTATATGTAAATCTTGTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGTGGCTTTTa                                                                              <  2:3255877/91‑1 (MQ=255)
               tAGAATTTTATATGTAAATCTTGTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGTGGCTTTTAg                                                                             <  2:658361/91‑1 (MQ=255)
                aGAATTTTATATGTAAATCTTGTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGTGGCTTTTAGt                                                                            >  1:707584/1‑91 (MQ=255)
                 gAATTTTATATGTAAATCTTGTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTa                                                                           >  1:1446832/1‑91 (MQ=255)
                  aaTTTTATATGTAAATCTTGTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTAt                                                                          >  2:2122330/1‑91 (MQ=255)
                     tttATATGTAAATCTTGTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGt                                                                       >  2:2593238/1‑91 (MQ=255)
                      ttATATGTAAATCTTGTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGtt                                                                      >  2:2450879/1‑91 (MQ=255)
                      ttATATGTAAATCTTGTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGtt                                                                      >  2:2669076/1‑91 (MQ=255)
                      ttATATGTAAATCTTGTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGTGGCTTTCAGTATCGtt                                                                      >  1:3709593/1‑91 (MQ=255)
                       tataTGTAAATCTTGTAGGTAATCGTTTTAAAGATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTc                                                                     <  1:3308218/91‑1 (MQ=255)
                       tataTGTAAATCTTGTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTc                                                                     <  2:1123432/91‑1 (MQ=255)
                         tatGTAAATCTTGTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTCCa                                                                   <  2:1079691/91‑1 (MQ=255)
                          atGTAAATCTTGTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTCCAg                                                                  <  1:2759964/91‑1 (MQ=255)
                              aaaTCTTGTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGAttt                                                              <  2:2193865/91‑1 (MQ=255)
                              aaaTCTTGTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGAttt                                                              <  2:105057/91‑1 (MQ=255)
                               aaTCTTGTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGATTtg                                                             <  2:390566/91‑1 (MQ=255)
                                 tCTTGTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGATTtgtg                                                           <  1:632053/91‑1 (MQ=255)
                                  cTTGTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGATTTGTgg                                                          >  2:1397770/1‑91 (MQ=255)
                                    tGTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGATTTGTgggg                                                        >  1:1035918/1‑91 (MQ=255)
                                    tGTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGATTTGTgggg                                                        >  2:456911/1‑91 (MQ=255)
                                    tGTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGATTTGTgggg                                                        <  2:1655377/91‑1 (MQ=255)
                                     gTAGGTAATCGTTTTAAAAATAATCTAAGTATGACTATAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGACTTGTGGGGa                                                       >  2:2752759/1‑91 (MQ=255)
                                     gTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGa                                                       >  2:1826485/1‑91 (MQ=255)
                                     gTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGa                                                       >  2:2905162/1‑91 (MQ=255)
                                      tAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAg                                                      >  2:1502744/1‑91 (MQ=255)
                                      tAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAg                                                      >  1:739690/1‑91 (MQ=255)
                                         gTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTa                                                   <  2:3804551/91‑1 (MQ=255)
                                          tAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAg                                                  >  2:3537491/1‑91 (MQ=255)
                                            aTCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAggg                                                <  1:1070588/91‑1 (MQ=255)
                                            aTCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAggg                                                <  1:1471333/91‑1 (MQ=255)
                                             tCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGt                                               >  1:2932831/1‑91 (MQ=255)
                                             tCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGt                                               <  2:1766957/91‑1 (MQ=255)
                                              cGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGtt                                              >  2:3293163/1‑91 (MQ=255)
                                               gTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGttt                                             >  2:231654/1‑91 (MQ=255)
                                                  ttAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGTTTACt                                          <  1:3404568/91‑1 (MQ=255)
                                                  ttAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGTTTACt                                          >  1:2085195/1‑91 (MQ=255)
                                                     aaaaTAATATAAGTATGACTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGTTTACttaa                                       >  2:1631992/1‑91 (MQ=255)
                                                     aaaaTAATATAAGCATGACTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGTTTACttaa                                       >  2:2365606/1‑91 (MQ=255)
                                                      aaaTAATATAAGTATGACTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGTTTACttaat                                      <  2:3713542/91‑1 (MQ=255)
                                                        ataataTAAGTATGACTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGATTCGTGGGGAGTTAGGGTTTACttaatta                                    >  2:3642530/1‑91 (MQ=255)
                                                         taataTAAGTATGACTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGTTTACttaattaa                                   >  2:141643/1‑91 (MQ=255)
                                                           atatAAGTATGACTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGTTTACTTAATTAAAg                                 >  1:1795907/1‑91 (MQ=255)
                                                            tataAGTATGACTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGTTTACTTAATTAAAgg                                >  1:1014982/1‑91 (MQ=255)
                                                             ataAGTATGACTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGTTTACTTAATTAAAggg                               >  2:2803351/1‑91 (MQ=255)
                                                              taAGTATGACTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGTTTACTTAATTAAAgggg                              >  1:2027422/1‑91 (MQ=255)
                                                              taAGTATGACTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGTTTACTTAATTAAAgggg                              >  1:515476/1‑91 (MQ=255)
                                                              taAGTATGACTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGTTTACTTAATTAAAgggg                              >  1:1648094/1‑91 (MQ=255)
                                                              taAGTATGACTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGAATTGAGGGGAGCTAGGGTTTACTTAATTAAAgggg                              >  2:2031465/1‑91 (MQ=255)
                                                                aGTATGACTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGACTTGTGGGGAGTTAGGGTTTACTTAATTAAAggggcc                            >  1:442970/1‑89 (MQ=255)
                                                                 gTATGACTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGTTTACTTAATTAAAGGGGTCt                           <  2:251171/91‑1 (MQ=255)
                                                                  tATGACTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGTTTACTTAATTAAAGGGGTCTg                          >  2:545344/1‑91 (MQ=255)
                                                                    tGACTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGTTTACTTAATTAAAGGGGTCTGcc                        >  1:1752646/1‑91 (MQ=255)
                                                                    tGACTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGTTTACTTAATTAAAGGGGTCTGcc                        <  1:2231050/91‑1 (MQ=255)
                                                                      aCTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGTTTACTTAATTAAAGGGGTCTGCCCt                      >  1:1730815/1‑91 (MQ=255)
                                                                       cTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGTTTACttaatta                                    >  1:1838468/1‑76 (MQ=255)
                                                                        tAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGTTTACTTAATTAAAGGGGTCTGCCCTCa                    <  1:1890267/91‑1 (MQ=255)
                                                                        aaaaGCCACGTCCTATATTGGATGTGGC‑TTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGTTTACTTAATTAAAGGGGTCTGCCCTCaa                   >  2:208011/2‑91 (MQ=25)
                                                                            gCCACGTCCTATATTGGATGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGTTTACTTAATTAAAGGGGTCTGCCCTCAATTg                >  1:3042426/1‑91 (MQ=38)
                                                                            gCCACGTCCTATATTGGATGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGTTTACTTAATTAAAGGGGTCTGCCCTCAATTg                >  2:2028791/1‑91 (MQ=38)
                                                                            gCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGTTTACTTAATTAAAGGGGTCTGCCCTCAATTg                >  1:35404/1‑91 (MQ=255)
                                                                            gCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGTTTACTTAATTAAAGGGGTCTGCCCTCAATTg                <  2:794201/91‑1 (MQ=255)
                                                                             ccACGTCCTATATTGGATGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGTTTACTTAATTAAAGGGGTCTGCCCTCAATTGc               >  2:242421/1‑91 (MQ=38)
                                                                             ccACATCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGTTTACTTAATTAAAGGGGTCTGCCCTCAATTGc               >  2:733556/1‑91 (MQ=255)
                                                                               acaTCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGTTTACTTAATTAAAGGGGTCTGCCCTCAATTGCTc             >  1:585088/1‑91 (MQ=255)
                                                                                caTCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGTTTACTTAATTAAAGCGGTCTGCCCTCAATTGCTCa            >  1:2183716/1‑91 (MQ=255)
                                                                                 ctCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGTTTACTTAATTAAAGGGGTCTGCCCTCAATTGCTCAc           <  2:1734156/90‑1 (MQ=255)
                                                                                 aTCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGTTTACTTAATTAAAGGGGTCTGCCCTCAATTGCTCAc           <  1:2632043/91‑1 (MQ=255)
                                                                                   tctatatTGGATGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGTGCTACGGTTTACTTAATTAAAGGGGTTTGCCCTCATTTGATCACCa         >  2:861966/4‑91 (MQ=255)
                                                                                   cctatatTGGATGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGTTTACTTAATTAAAGGGGTCTGCCCTCAATTGCTCACCa         >  1:2441687/4‑91 (MQ=255)
                                                                                   ccAATATAGGACGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGTTTACTTAATTAAAGGGGTCTGCCCTCAATTGCTCACCa         >  2:1205562/1‑91 (MQ=255)
                                                                                    ctatatTGGATGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGTTTACTTAATTAAAGGGGTCTGCCCTCAATTGCTCACCaa        >  1:1658925/3‑91 (MQ=255)
                                                                                    cAATATAGGACGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGTTTACTTAATTAAAGGGGTCTGCCCTCGATTGCTCACCaa        >  2:3223575/1‑91 (MQ=255)
                                                                                     aaTATAGGACGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGTTTACTTAATTAAAGGGGTCTGCCCTCAATTGCTCACCAAc       <  1:322670/91‑1 (MQ=255)
                                                                                      atatTGGATGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGTTTACTTAATTAAAGGGGTCTGCCCTCAATTGCTCACCAACg      >  1:2279403/1‑91 (MQ=255)
                                                                                      atatTGGATGTGGC‑TTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGTTTACTTAATTAAAGGGGTCTGCCCTCAATTGCTCACCAACGa     >  1:3847199/1‑91 (MQ=255)
                                                                                      atatAGGACGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGTTTACTTAATTAAAGGGGTCTGCCCTCAATTGCTCACCAACg      >  1:2142202/1‑91 (MQ=255)
                                                                                       tataGGACGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGTTTACTTAATTAAAGGGGTCTGCCCTCAATTGCTCACCAACGa     <  2:741761/91‑1 (MQ=255)
                                                                                        ataGGACGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGTTTACTTAATTAAAGGGGTCTGCCCTCAATTGCTCACCAACGAg    >  1:1397383/1‑91 (MQ=255)
                                                                                        ataGGACGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGTTTACTTAATTAAAGGGGTCTGCCCTCAATTGCTCACCAACGAg    >  2:3499993/1‑91 (MQ=255)
                                                                                          aGGACGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGTTTACTTAATTAAAGGGGTCTGCCCTCAATTGCTCACCAACGAggg  <  1:2392382/91‑1 (MQ=255)
                                                                                          aGGACGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGTTTACTTAATTAAAGGGGTCTGCCCTCAATTGCTCACCAACGAggg  >  1:824536/1‑91 (MQ=255)
                                                                                          |                                                                                          
ATATAGCAGACATGATAGAATTTTATATGTAAATCTTGTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGTTTACTTAATTAAAGGGGTCTGCCCTCAATTGCTCACCAACGAGGG  >  USA300TCH1516_ALE/2608735‑2608915

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: ‑