Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 228,319 T→C 8.8% intergenic (+165/+29) ybbH_1 → / ← USA300TCH1516_00200 putative HTH‑type transcriptional regulator YbbH/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE228,3190TC 8.8% 802.3 / 50.3 274intergenic (+165/+29)ybbH_1/USA300TCH1516_00200putative HTH‑type transcriptional regulator YbbH/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (131/118);  new base C (11/13);  total (143/131)
Fisher's exact test for biased strand distribution p-value = 6.70e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.75e-01

Only 100 of 282 total aligned reads displayed.

ATGAAAAATTGTTGCTACAAGTTTGTAAAACTGGCATAGAGGCTATTTTGAAAATGATAAAGGCGTGACAAATTAAAGCAAAGTTCAGATTTGAGCTTTGCTTTTTTCATAGATGGAAATTAAATGCTGAAGAATAATTTTAATGCTAAATGAACTAATACGAAGGTGCTAAAAGCAATT  >  USA300TCH1516_ALE/228229‑228408
                                                                                          |                                                                                         
atgaAAAATTGTTGCTACAAGTTTGTAAAACTGGCATAGAGGCTATTTTGAAAATGATAAAGGCGTGACAAATTAAAGCAAAGTTCAGAtt                                                                                           <  1:2855275/91‑1 (MQ=255)
atgaAAAATTGTTGCTACAAGTTTGTAAAACTGGCATAGAGGCTATTTTGAAAATGATAAAGGCGTGACAAATTAAAGCAAAGTTCAGAtt                                                                                           >  1:170237/1‑91 (MQ=255)
 tgaAAAATTGTTGCTACAAGTTTGTAAAACTGGCATAGAGGCTATTTTGAAAATGATAAAGGCGTGACAAATTAAAGCAAAGTTCAGAttt                                                                                          >  2:1690740/1‑91 (MQ=255)
     aaaTTGTTGCTACAAGTTTGTAAAACTGGCATAGAGGCTATTTTGAAAATGATAAAGGCGTGACAAATTAAAGCAAAGTTCAGATTTGAgc                                                                                      <  2:800507/91‑1 (MQ=255)
     aaaTTGTTGCTACAAGTTTGTAAAACTGGCATAGAGGCTATTTTGAAAATGATAAAGGCGTGACAAATTAAAGCAAAGTTCAGATTTGAgc                                                                                      >  1:1352173/1‑91 (MQ=255)
     aaaTTGTTGCTACAAGTTTGTAAAACTGGCATAGAGGCTATTTTGAAAATGATAAAGGCGTGACAAATTAAAGCAAAGTTCAGATTTGAgc                                                                                      >  2:3885677/1‑91 (MQ=255)
      aaTTGTTGCTACAAGTTTGTAAAACTGGCATAGAGGCTATTTTGAAAATGATAAAGGCGTGACAAATTAAAGCAAAGTTCAGATTTGAgct                                                                                     >  2:3557347/1‑91 (MQ=255)
      aaTTGTTGCTACAAGTTTGTAAAACTGGCATAGAGGCTATTTTGAAAATGATAAAGGCGTGACAAATTAAAGCAAAGTTCAGATTTGAgct                                                                                     >  1:668014/1‑91 (MQ=255)
      aaTTGTTGCTACAAGTTTGTAAAACTGGCATAGAGGCTATTTTGAAAATGATAAAGGCGTGACAAATTAAAGCAAAGTTCAGATCTGaact                                                                                     <  2:1710580/91‑4 (MQ=255)
       aTTGTTGCTACAAGTTTGTAAAACTGGCATAGAGGCTATTTTGAAAATGATAAAGGCGTGACAAATTAAAGCAAAGTTCAGATTTGAgctt                                                                                    <  1:3741649/91‑1 (MQ=255)
         tgttgCTACAAGTTTGTAAAACTGGCATAGAGGCTATTTTGAAAATGATAAAGGCGTGACAAATTAAAGCAAAGTTCAGATTTGAgctttg                                                                                  >  2:1390769/1‑91 (MQ=255)
         tgttgCTAAAAGTTTATAAAACTGGCAGATAGGATATTTAGAAAAAGATAAAAACATAAAAAAATAAAGCAAAGTTCAGATTTGAgctttg                                                                                  <  2:680586/91‑1 (MQ=255)
          gttgCTACAAGTTTGTAAAACTGGCATAGAGGCTATTTTGAAAATGATAAAGGCGTGACAAATTAAAGCAAAGTTCAGATTTGAgctttgc                                                                                 >  2:3357007/1‑91 (MQ=255)
           ttgCTACAAGTTTGTAAAACTGGCATAGAGGCTATTTTGAAAATGATAAAGGCGTGACAAATTAAAGCAAAGTTCAGATTTGAgctttgct                                                                                >  1:2067531/1‑91 (MQ=255)
           ttgCTACAAGTTTGTAAAACTGGCATAGAGGCTATTTTGAAAATGATAAAGGCGTGACAAATTAAAGCAAAGTTCAGATTTGAgctttgct                                                                                >  1:935355/1‑91 (MQ=255)
            tgCTACAAGTTTGTAAAACTGGCATAGAGGCTATTTTGAAAATGATAAAGGCGTGACAAATTAAAGCAAAGTTCAGATTTGAgctttgctt                                                                               >  2:1214429/1‑91 (MQ=255)
            tgCTACAAGTTTGTAAAACTGGCATAGAGGCTATTTTGAAAATGATAAAGGCGTGACAAATTAAAGCAAAGTTCAGATTTGAgctttgctt                                                                               <  1:886569/91‑1 (MQ=255)
             gCTACAAGTTTGTAAAACTGGCATAGAGGCTATTTTGAAAATGATAAAGGCGTGACAAATTAAAGCAAAGTTCAGATTTGAgcttt                                                                                   >  1:603918/1‑86 (MQ=255)
             gCTACAAGTTTGTAAAACTGGCATAGAGGCTATTTTGAAAATGATAAAGGCGTGACAAATTAAAGCAAAGTTCAGATTTGAgcttt                                                                                   <  2:603918/86‑1 (MQ=255)
              cTACAAGTTTGTAAAACTGGCATAGAGGCTATTTTGAAAATGATAAAGGCGTGACAAATTAAAGCAAAGTTCAGATTTGAGCTTTGCtttt                                                                             >  2:3026946/1‑91 (MQ=255)
              cTACAAGTTTGTAAAACTGGCATAGAGGCTATTTTGAAAATGATAAAGGCGTGACAAATTAAAGCAAAGTTCAGATTTGAGCTTTGCtttt                                                                             <  2:2592230/91‑1 (MQ=255)
              cTACAAGTTTGTAAAACTGGCATAGAGGCTATTTTGAAAATGATAAAGGCGTGACAAATTAAAGCAAAGCTCAAATCTGAactttgcttta                                                                             <  2:2750062/91‑2 (MQ=37)
                aCAAGTTTGTAAAACTGGCATAGAGGCTATTTTGAAAATGATAAAGGCGTGACAAATTAAAGCAAAGTTCAGATTTGCGCTTTGCtttttt                                                                           >  1:3871790/1‑91 (MQ=255)
                aCAAGTTTGTAAAACTGGCATAGAGGCTATTTTGAAAATGATAAAGGCGTGACAAATTAAAGCAAAGCTCAAATCTGAactttgctttaat                                                                           <  1:1178159/91‑4 (MQ=25)
                 cAAGTTTGTAAAACTGGCATAGAGGCTATTTTGAAAATGATAAAGGCGTGACAAATTAAAGCAAAGCTCAAATCTGAactttgctttaatt                                                                          <  2:1011984/91‑5 (MQ=25)
                  aaGTTTGTAAAACTGGCATAGAGGCTATTTTGAAAATGATAAAGGCGTGACAAATTAAAGCAAAGTTCAGATTTGAGCTTTGCTTTTTTCa                                                                         <  2:1370710/91‑1 (MQ=255)
                    gTTTGTAAAACTGGCATAGAGGCTATTTTGAAAATGATAAAGGCGTGACAAATTAAAGCAAAGTTCAGATTTGAGCTTTGCTTTTTTCATa                                                                       >  2:2529913/1‑91 (MQ=255)
                    gTTTGTAAAACTGGCATAGAGGCTATTTTGAAAATGATAAAGGCGTGACAAATTAAAGCAAAGTTCAGATTTGAGCTTTGCTTTTTTCATa                                                                       >  1:3432635/1‑91 (MQ=255)
                      ttGTAAAACTGGCATAGAGGCTATTTTGAAAATGATAAAGGCGTGACAAATTAAAGCAAAGTTCAGATTTGAGCTTTGCTTTTTTCATAGa                                                                     <  2:1497422/91‑1 (MQ=255)
                       tGTAAAACTGGCATAGAGGCTATTTTGAAAGTGATAAAGGCGTGACAAATTATAGCAAAGTTCAGATTTGGGCTTTGCTTTTTTCATAGAt                                                                    <  1:2983852/91‑1 (MQ=255)
                       tGTAAAACTGGCATAGAGGCTATTTTGAAAATGATAAAGGCGTGACAAATTAAAGCAAAGTTCAGATTTGAGCTTTGCTTTTTTCATAGAt                                                                    >  1:3521124/1‑91 (MQ=255)
                        gTAAAACTGGCATAGAGGCTATTTTGAAAATGATAAAGGCGTGACAAATTAAAGCAAAGTTCAGATTTGAGCTTTGCTTTTTTCATAGATg                                                                   >  1:1522239/1‑91 (MQ=255)
                         tAAAACTGGCATAGAGGCTATTTTGAAAATGATAAAGGCGTGACAAATTAAAGCAAAGTTCAGATTTGAGCTTTGCTTTTTTCATAGATgg                                                                  <  2:1653393/91‑1 (MQ=255)
                          aaaaCTGGCATAGAGGCTATTTTGAAAATGATAAAGGCGTGACAAATTAAAGCAAAGTTCAGATTTGAGCTTTGCTTTTTTCATAGATGGa                                                                 <  1:3017766/91‑1 (MQ=255)
                          aaaaCTGGCATAGAGGCTATTTTGAAAATGATAAAGGCGTGACAAATTAAAGCAAAGTTCAGATTTGAGCTTTGCTTTTTTCATAGATGGa                                                                 <  1:2451795/91‑1 (MQ=255)
                          aaaaCTGGCATAGAGGCTATTTTGAAAATGATAAAGGCGTGACAAATTAAAGCAAAGTTCAGATTTGAGCTTTGCTTTTTTCATAGATGGa                                                                 <  1:2465009/91‑1 (MQ=255)
                          aaaaCTGGCATAGAGGCTATTTTGAAAATGATAAAGGCGTGACAAATTAAAGCAAAGTTCAGATTTGAGCTTTGCTTTTTTCATAGATGGa                                                                 >  1:1718489/1‑91 (MQ=255)
                          aaaaCTGGCATAGAGGCTATTTTGAAAATGATAAAGGCGTGACAAATTAAAGCAAAGTTCAGATTTGAGCTTTGCTTTTTTCATAGATGGa                                                                 >  1:454193/1‑91 (MQ=255)
                          aaaaCTGGCATAGAGGCTATTTTGAAAATGATAAAGGCGTGACAAATTAAAGCAAAGTTCAGATTAGAGCTTTGCTTTTTTCATAGATGGa                                                                 >  1:3621520/1‑91 (MQ=255)
                           aaaCTGGCATAGAGGCTATTTTGAAAATGATAAAGGCGTGACAAATCAAAGCAAAGTTCAGATTTGAGCTTTGCTTTTTTCATAGATGGaa                                                                >  1:3083874/1‑91 (MQ=255)
                              cTGGCATAGAGGCTATTTTGAAAATGATAAAGGCGTGACAAATTAAAGCAAAGTTCAGATTTGAGCTTTGCTTTTTTCATAGATGGAAAtt                                                             <  1:327476/91‑1 (MQ=255)
                               tGGCATAGAGGCTATTTTGAAAATGATAAAGGCGTGACAAATTAAAGCAAAGTTCAGATTTGAGCTTTGCTTTTTTCATAGATGGAAATTa                                                            <  2:3496697/91‑1 (MQ=255)
                                   aTAGAGGCTATTTTGAAAATGATAAAGGCGTGACAAATTAAAGCAAAGTTCAGATTTGAGCTTTGCTTTTTTCATAGATGGAAATTAAATg                                                        >  1:146986/1‑91 (MQ=255)
                                       agGCTATTTTGAAAATGATAAAGGCGTGACAAATTAAAGCAAAGTTCAGATTTGAGCTTTGCTTTTTTCATAGATGGaaa                                                               <  2:1329680/80‑1 (MQ=255)
                                       agGCTATTTTGAAAATGATAAAGGCGTGACAAATTAAAGCAAAGTTCAGATTTGAGCTTTGCTTTTTTCATAGATGGAAATTAAATGCTga                                                    <  2:3320909/91‑1 (MQ=255)
                                       agGCTATTTTGAAAATGATAAAGGCGTGACAAATTAAAGCAAAGTTCAGATTTGAGCTTTGCTTTTTTCATAGATGGAAATTAAATGCTga                                                    >  1:3641959/1‑91 (MQ=255)
                                        ggCTATTTTGAAAATGATAAAGGCGTGACAAATTAAAGCAAAGTTCAGATTTGAGCTTTGCTTTTTTCATAGATGGAAATTAAATGCTgaa                                                   <  1:1895213/91‑1 (MQ=255)
                                        ggCTATTTTGAAAATGATAAAGGCGTGACAAATTAAAGCAAAGTTCAGATTTGAGCTTTGCTTTTTTCATAGATGGAAATTAAATGCTgaa                                                   <  2:2491421/91‑1 (MQ=255)
                                        ggCTATTTTGAAAATGATAAAGGCGTGACAAATTAAAGCAAAGTTCAGATTTGAGCTTTGCTTTTTTCATAGATGGAAATTAAATGCTgaa                                                   <  1:2642545/91‑1 (MQ=255)
                                          cTATTTTGAAAATGATAAAGGCGTGACAAATTAAAGCAAAGTTCAGATTTGAGCTTTGCTTTTTTCATAGATGGAAATTAAATGCTgaaga                                                 >  1:2200879/1‑91 (MQ=255)
                                          cTATTTTGAAAATGATAAAGGCGTGACAAATTAAAGCAAAGTTCAGATTTGAGCTTTGCTTTTTTCATAGATGGAAATTAAATGCTgaaga                                                 >  1:1010672/1‑91 (MQ=255)
                                           tATTTTGAAAATGATAAAGGCGTGACAAATTTAAGCAAAGTTCAGATTTGAGCTTTGGTTTTTTCATAGATGGAAATTAAATGCTgaagaa                                                <  1:3102920/91‑1 (MQ=255)
                                           tATTTTGAAAATGATAAAGGCGTGACAAATTAAAGCAAAGTTCAGATTTGAGCTTTGCTTTTTTCATAGATGGAAATTAAATGCTgaagaa                                                <  2:1946825/91‑1 (MQ=255)
                                            aTTTTGAAAATGATAAAGGCGTGACAAATTAAAGCAAAGTTCAGATTTGAGCTTTGCTTTTTTCATAGATGGAAATTAAATGCTGAAGaat                                               >  1:2774412/1‑91 (MQ=255)
                                                tGAAAATGATAAAGGCGTGACAAATTAAAGCAAAGTTCAGATTTGAGCTTTGCTTTTTTCATAGATGGAAATTAAATGCTGAAGAATAAtt                                           >  1:740082/1‑91 (MQ=255)
                                                      tGATAAAGGCGTGACAAATTAAAGCAAAGTTCAGATTTGAGCTTTGCTTTTTTCATAGATGGAAATTAAATGCTGAAGAATAATTTTAATg                                     >  1:297744/1‑91 (MQ=255)
                                                      tGATAAAGGCGTGACAAATTAAAGCAAAGTTCAGATTTGAGCTTTGCTTTTTTCATAGATGGAAATTAAATGCTGAAGAATAATTTTAATg                                     <  1:3357007/91‑1 (MQ=255)
                                                       gATAAAGGCGTGACAAATTAAAGCAAAGTTCAGATTTGAGCTTTGCTTTTTTCATAGATGGAAATTAAATGCTGAAGAATAATTTTAATGc                                    <  2:1882977/91‑1 (MQ=255)
                                                       gATAAAGGCGTGACAAATTAAAGCAAAGTTCAGATTTGAGCTTTGCTTTTTTCATAGATGGAAATTAAATGCTGAAGAATAATTTTAATGc                                    <  2:1794156/91‑1 (MQ=255)
                                                        aTAAAGGCGTGACAAATTAAAGCAAAGTTCAGATTTGAGCTTTGCTTTTTTCATAGATGGAAATTAAATGCTGAAGAATAATTTTAATGCt                                   >  2:1136720/1‑91 (MQ=255)
                                                          agaGGCGTGACAAATTAAAGCAAAGTTCAGATTTGAGCTTTGCTTTTTTCATAGATGGAAATTAAATGCTGAAGAATAATTTTAATGCTaa                                 >  2:3221839/3‑91 (MQ=255)
                                                          aaaGGCGTGACAAATTAAAGCAAAGTTCAGATTTGAGCTTTGCTTTTTTCATAGATGGAAATTAAATGCTGAAGAATAATTTTAATGCTaa                                 >  1:2347041/1‑91 (MQ=255)
                                                           aaGGCGTGACAAATTAAAGCAAAGTTCAGATTTGAGCTTTGCTTTTTTCATAGATGGAAATTAAATGCTGAAGAATAATTTTAATGCTaaa                                >  1:1803309/1‑91 (MQ=255)
                                                            aGGCGTGACAAATTAAAGCAAAGTTCAGATTTGAGCTTTGCTTTTTTCATAGATGGAAATTAAATGCTGAAGAATAATTTTAATGCTAAAt                               >  2:2823881/1‑91 (MQ=255)
                                                             ggCGTGACAAATTAAAGCAAAGTTCAGATTTGAGCTTTGCTTTTTTCATAGATGGAAATTAAATGCTGAAGAATAATTTTAATGCTAAATg                              >  1:87182/1‑91 (MQ=255)
                                                                gTGACAAATTAAAGCAAAGTTCAGGTTTGAGCTTTGCTTTTTTCATAGATGGAAATTAAATGCTGAAGAATAATTTTAATGCTAAATGAAc                           >  2:2084653/1‑91 (MQ=255)
                                                                 tGACAAATTAAAGCAAAGTTCAGATTTGAGCTTTGCTTTTTTCATAGATGGAAATTAAATGCTGAAGAATAATTTTAATGCTAAATGAACt                          <  2:2193972/91‑1 (MQ=255)
                                                                   aCAAATTAAAGCAAAGTTCAGATTTGAGCTTTGCTTTTTTCATAGATGGAAATTAAATGCTGAAGAATAATTTTAATGCTAAATGAACTaa                        <  1:3683173/91‑1 (MQ=255)
                                                                   aCAAATTAAAGCAAAGTTCAGATTTGAGCTTTGCTTTTTTCATAGATGGAAATTAAATGCTGAAGAATAATTTTAATGCTAAATGAACTaa                        >  1:2949915/1‑91 (MQ=255)
                                                                   aCAAATTAAAGCAAAGTTCAGATTTGAGCTTTGCTTTTTTCATAGATGGAAATTAAATGCTGAAGAATAATTTTAATGCTAAATGAACTaa                        >  2:2998715/1‑91 (MQ=255)
                                                                   aCAAATTAAAGCAAAGTTCAGATTTGAGCTTTGCTTTTTTCATAGATGGAAATTAAATGCTGAAGAATAATTTTAATGCTAAATGAACTaa                        >  2:1584959/1‑91 (MQ=255)
                                                                    cAAATTCAAGCAAAGTTCAGATTTGAGCTTTGCTTTTTTCATAGATGGAAATTAAATGTAGAAGAATTATTTTAATGCCAAATGAACTAAt                       >  2:553509/1‑91 (MQ=255)
                                                                    cAAATTAAAGCAAAGTTCAGATTTGAGCTTTGCTTTTTTCATAGATGGAAATTAAATGCTGAAGAATAATTTTAATGCTAAATGAACTAAt                       >  2:1071399/1‑91 (MQ=255)
                                                                     aaaTTAAAGCAAAGTTCAGATTTGAGCTTTGCTTTTTTCATAGATGGAAATTAAATGCTGAAGAATAATTTTAATGCTAAATGAACTAATa                      <  2:2051847/91‑1 (MQ=255)
                                                                      aaTTAGAGCAAAGTTCAGATTTGAGCTTTGCTTTTTTCATAGATGGAAATTAAATGCTGAAGAATAATTTTAATGCTAAATGAACTAATAc                     <  2:3146717/91‑1 (MQ=255)
                                                                      aaTTAAAGCAAAGTTCAGATTTGAGCTTTGCTTTTTTCATAGATGGAAATTAAATGCTGAAGAATAATTTTAATGCTAAATGAACTAATAc                     >  2:843471/1‑91 (MQ=255)
                                                                        ttAAAGCAAAGTTCAGATTTGAGCTTTGCTTTTTTCATAGATGGAAATTAAATGCTGAAGAATAATTTTAATGCTAAATGAACTAATACGa                   >  2:2971010/1‑91 (MQ=255)
                                                                          aaaGCAAAGTTCAGATTTGAGCTTTGCTTTTTTCATAGATGGAAATTAAATGCTGAAGAATAATTTTAATGCTAAATGAACTAATACGAAg                 >  1:2578500/1‑91 (MQ=255)
                                                                           aaGCAAAGTTCAGATTTGAGCTTTGCTTTTTTCATAGATGGAAATTAAATGCTGAAGAATAATTTTAATGCTAAATGAACTAATACGAAgg                <  2:1826088/91‑1 (MQ=255)
                                                                           aaGCAAAGTTCAGATTTGAGCTTTGCTTTTTTCATAGATGGAAATTAAATGCTGAAGAATAATTTTAATGCTAAATGAACTAATACGAAgg                <  1:3708674/91‑1 (MQ=255)
                                                                           aaGCAAAGCTCAAATCTGAACTTTGCTTTTTTCATAGATGGAAATTAAATGCTGAAGAATAATTTTAATGCTAAATGAACTAATACGAAgg                >  1:1599106/1‑91 (MQ=37)
                                                                            aGCAAAGTTCAGATTTGAGCTTTGCTTTTTTCATAGATGGAAATTAAATGCTGAAGAATAATTTTAATGCTAAATGAACTAATACGAAGGt               <  1:2514560/91‑1 (MQ=255)
                                                                            aGCAAAGTTCAGATTTGAGCTTTGCTTTTTTCATAGATGGAAATTAAATGCTGAAGAATAATTTTAATGCTAAATGAACTAATACGAAGGt               >  1:532276/1‑91 (MQ=255)
                                                                             gCAAAGCTCAAATCTGAACTTTGCTTTTTTCATAGATGGAAATTAAATGCTGAAGAATAATTTTAATGCTAAATGAACTAATACGAAGGTg              >  2:368552/1‑91 (MQ=38)
                                                                              cAAAGTTCAGATTTGAGCTTTGCTTTTTTCATAGATGGAAATTAAATGCTGAAGAATAATTTTAATGCTAAATGAACTAATACGAAGGTGc             <  2:455924/91‑1 (MQ=255)
                                                                              cAAAGTTCAGATTTGAGCTTTGCTTTTTTCATAGATGGAAATTAAATGCTGAAGAATAATTTTAATGCTAAATGAACTAATACGAAGGTGc             >  1:3551491/1‑91 (MQ=255)
                                                                              cAAAGCTCAAATCTGAACTTTGCTTTTTTCATAGATGGAAATTAAATGCTGAAGAATAATTTTAATGCTAAATGAACTAATACGAAGGTGc             >  2:2895799/1‑91 (MQ=38)
                                                                              cAAAGCTCAAATCTGAACTTTGCTTTTTTCATAGATGGAAATTAAATGATGAAGACTAATTTTAATGCTAAATGAACTAATACGAAGGTgg             >  2:2723399/1‑90 (MQ=37)
                                                                               aaaGTTCAGATTTGAGCTTTGCTTTTTTCATAGATGGAAATTAAATGCTGAAGAATAATTTTAATGCTAAATGAACTAATACGAAGGTGCt            <  1:1324262/91‑1 (MQ=255)
                                                                                aaGTTCAGATTTGAGCTTTGCTTTTTTCATAGATGGAAATTAAATGCTGAAGaat                                               >  2:3635995/1‑55 (MQ=255)
                                                                                aaGTTCAGATTTGAGCTTTGCTTTTTTCATAGATGGAAATTAAATGCTGAAGAATAATTTTAATGCTAAATGAACTAATACGAAGGTGCTa           <  1:1883452/91‑1 (MQ=255)
                                                                                 aGTTCAGATTTGAGCTTTGCTTTTTTCATAGATGGAAATTAAATGCTGAAGAATAATTTTAATGCTAAATGAACTAATACGAAGGTGCTaa          >  2:877846/1‑91 (MQ=255)
                                                                                  gTTCAGATTTGAGCTTTGCTTTTTTCATAGATGGAAATTAAATGCTGAAGAATAATTTTAATGCTAAATGAACTAATACGAAGGTGCTaaa         <  2:66606/91‑1 (MQ=255)
                                                                                  gTTCAGATTTGAGCTTTGCTTTTTTCATAGATGGAAATTAAATGCTGAAGAATAATTTTAATGCTAAATGAACTAATACGAAGGTGCTaaa         >  1:1273419/1‑91 (MQ=255)
                                                                                   ttCAGATTTGAGCTTTGCTTTTTTCATAGATGGAAATTAAATGCTGAAGAATAATTTTAATGCTAAATGAACTAATACGAAGGTGCTaaaa        >  2:3769505/1‑91 (MQ=255)
                                                                                    tCAGATTTGAGCTTTGCTTTTTTCTTCGATGGAAATTAAATGCTGAAGAATAATTTTAATGCTAAATGAACTAATACGAAGGTGCTAAAAg       >  2:2215207/1‑91 (MQ=255)
                                                                                    tCAGATTTGAGCTTTGCTTTTTTCATAGATGGAAATTAAATGCTGAAGAATAATTTTAATGCTAAATGAACTAATACGAAGGTGCTAAAAg       >  2:1276768/1‑91 (MQ=255)
                                                                                     cAGATTTGAGCTTTGCTTTTTTCATAGATGGAAATTAAATGCTGAAGAATAATTTTAATGCTAAATGAACTAATACGAAGGTGCTAAAAGc      >  2:2939124/1‑91 (MQ=255)
                                                                                       gATTTGAGCTTTGCTTTTTTCATAGATGGAAATTAAATGCTGAAGAATAATTTTAATGCTAAATGAACTAATACGAAGGTGCTAAAAGCaa    >  2:864051/1‑91 (MQ=255)
                                                                                         tttGAGCTTTGCTTTTTTCATAGATGGAAATTAAATGCTGAAGAATAATTTTAATGCTAAATGAACTAATACGTAGGTGCTAAAAGCAatt  >  1:3323388/1‑91 (MQ=255)
                                                                                          |                                                                                         
ATGAAAAATTGTTGCTACAAGTTTGTAAAACTGGCATAGAGGCTATTTTGAAAATGATAAAGGCGTGACAAATTAAAGCAAAGTTCAGATTTGAGCTTTGCTTTTTTCATAGATGGAAATTAAATGCTGAAGAATAATTTTAATGCTAAATGAACTAATACGAAGGTGCTAAAAGCAATT  >  USA300TCH1516_ALE/228229‑228408

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