Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 2,559,850 T→C 8.9% intergenic (+44/+79) bcr_2 → / ← USA300HOU_RS13190 Bicyclomycin resistance protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE2,559,8500TC 8.9% 560.9 / 35.6 191intergenic (+44/+79)bcr_2/USA300HOU_RS13190Bicyclomycin resistance protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (92/82);  new base C (13/4);  total (105/86)
Fisher's exact test for biased strand distribution p-value = 7.57e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.06e-01

Only 100 of 195 total aligned reads displayed.

ACAAATCATTTACTTTAAAATGATTAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTCATAGCAAGGACATTTCTATTCAATCACCCTTTAAC  >  USA300TCH1516_ALE/2559761‑2559940
                                                                                         |                                                                                          
aCAAATCATTTACTTTAAAATGATTAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTc                                                                                           >  2:2180685/1‑91 (MQ=255)
aCAAATCATTTACTTTAAAATGATTAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTc                                                                                           <  2:339189/91‑1 (MQ=255)
 cAAATCATTTACTTTAAAATGATTAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCa                                                                                          >  2:389742/1‑91 (MQ=255)
  aaaTCATTTACTTTAAAATGATTAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAg                                                                                         <  1:2452437/91‑1 (MQ=255)
  aaaTCATTTACTTTAAAATGATTAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAg                                                                                         >  2:1601247/1‑91 (MQ=255)
     tCATTTACTTTAAAATGATTAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTGACAAAGATAGATATCTTTGCTAGGGt                                                                                      <  1:583210/91‑1 (MQ=25)
     tCATTTACTTTAAAATGATTAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTc                                                                                           <  1:2301718/86‑1 (MQ=255)
     tCATTTACTTTAAAATGATTAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTc                                                                                           >  2:2301718/1‑86 (MQ=255)
      cATTTACTTTAAAATGATTAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGtt                                                                                     >  1:1794935/1‑91 (MQ=255)
        tttACTTTAAAATGATTAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGttct                                                                                   >  1:846819/1‑91 (MQ=255)
          tACTTTAAAATGATTAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGttcttc                                                                                 >  1:2089875/1‑91 (MQ=255)
            cTTTAAAATGATTAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTcttctt                                                                               >  2:2699088/1‑91 (MQ=255)
             tttAAAATGATTAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCttt                                                                              <  1:2041042/91‑1 (MQ=255)
             tttAAAATGATTAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCttt                                                                              >  1:1552095/1‑91 (MQ=255)
             tttAAAATGATTAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCttt                                                                              >  1:1200685/1‑91 (MQ=255)
             tttAAAATGATTAAAAAGCAACATGTCGCATAACACTTCAACATAAGTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCttt                                                                              <  2:3170456/91‑1 (MQ=255)
              ttAAAATGATTAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTa                                                                             <  2:3001168/91‑1 (MQ=255)
                  aaTGATTAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGaa                                                                         >  2:1880978/1‑91 (MQ=255)
                    tGATTAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAAtt                                                                       >  2:2244403/1‑91 (MQ=255)
                     gATTAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTa                                                                      >  1:1616999/1‑91 (MQ=255)
                     gATTAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTa                                                                      <  1:218980/91‑1 (MQ=255)
                      aTTAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTAt                                                                     >  1:1877137/1‑91 (MQ=255)
                      aTTAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTAt                                                                     <  1:1206406/91‑1 (MQ=255)
                        tAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATga                                                                   >  1:2894949/1‑91 (MQ=255)
                          aaaaGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATgaga                                                                 <  1:1031882/91‑1 (MQ=255)
                          aaaaGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATgaga                                                                 >  1:1576757/1‑91 (MQ=255)
                          aaaaGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATgaga                                                                 <  2:2493123/91‑1 (MQ=255)
                           aaaGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGAt                                                                >  1:2768775/1‑91 (MQ=255)
                           aaaGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGATCTTCTTTATGAATTATGAGAt                                                                >  1:752947/1‑91 (MQ=255)
                            aaGCAACATGTCGCATAACACTTCAACATAATTAGAGCCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATc                                                               <  2:3052474/91‑1 (MQ=255)
                              gCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGa                                                             <  1:3421499/91‑1 (MQ=255)
                              gCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGa                                                             <  1:1913945/91‑1 (MQ=255)
                               cAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGaa                                                            >  1:450292/1‑91 (MQ=255)
                                 aCATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATc                                                          >  2:3205472/1‑91 (MQ=255)
                                   aTGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCtt                                                        >  2:2137577/1‑91 (MQ=255)
                                   aTGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCtt                                                        >  2:2382582/1‑91 (MQ=255)
                                   aTGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCtt                                                        >  1:323327/1‑91 (MQ=255)
                                       cGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAACCTTCAAc                                                    >  2:3450467/1‑91 (MQ=255)
                                        gCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACt                                                   <  1:823324/91‑1 (MQ=255)
                                         cATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTa                                                  <  2:3734548/91‑1 (MQ=255)
                                          aTAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTaa                                                 >  1:2965607/1‑91 (MQ=255)
                                            aaCACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTaaaa                                               >  2:263851/1‑91 (MQ=255)
                                             acacTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAAt                                              <  1:830691/91‑1 (MQ=255)
                                              cacTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAAtt                                             <  1:3767551/91‑1 (MQ=255)
                                               acTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTa                                            <  1:2073301/91‑1 (MQ=255)
                                                 ttCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACg                                          >  2:1875919/1‑91 (MQ=255)
                                                  tCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGc                                         <  1:2137577/91‑1 (MQ=255)
                                                   cAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGcc                                        >  2:2089490/1‑91 (MQ=255)
                                                   cAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGcc                                        >  2:1848473/1‑91 (MQ=255)
                                                   cAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGcc                                        >  2:3567056/1‑91 (MQ=255)
                                                    aaCATAATTAGAACCCTAGCAAAGATATCGATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCt                                       <  2:722813/91‑1 (MQ=255)
                                                     aCATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCtt                                      >  1:2814596/1‑91 (MQ=255)
                                                      cATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTc                                     >  2:3507781/1‑91 (MQ=255)
                                                          aTTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTCATAg                                 >  1:2867693/1‑91 (MQ=255)
                                                          aTTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTCATAg                                 >  1:2368826/1‑91 (MQ=255)
                                                           ttAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTCATAGc                                <  1:3597163/91‑1 (MQ=255)
                                                           ttAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTCATAGc                                >  1:2926884/1‑91 (MQ=255)
                                                            tAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTCATAGCa                               <  1:543556/91‑1 (MQ=255)
                                                             aGAACCCTGACAAAGATAGATATCTTTGCTAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTCATAGCaa                              >  1:3205959/1‑91 (MQ=21)
                                                             aGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCTAATCTTCAACTAAAATTACGCCTTCATAGCaa                              >  1:3300768/1‑91 (MQ=255)
                                                             aGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTCATAGCaa                              >  2:845037/1‑91 (MQ=255)
                                                              gAACCCTGACAAAGATAGATATCTTTGCTAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTCATAGCAAg                             >  1:188455/1‑91 (MQ=21)
                                                               aaCCCTGACAAAGATAGATATCTTTGCTAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTCATAGCAAgg                            >  1:31005/1‑91 (MQ=21)
                                                                aCCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTCATAGCAAGGa                           <  1:2014648/91‑1 (MQ=255)
                                                                aCCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTCATAGCAAGGa                           <  1:2551231/91‑1 (MQ=255)
                                                                aCCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTCATAGCAAGGa                           >  1:2714702/1‑91 (MQ=255)
                                                                 cccTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTCATAGCAAGGAc                          <  2:1277095/91‑1 (MQ=255)
                                                                  ccTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTCATAGCAAGGACa                         >  1:625166/1‑91 (MQ=255)
                                                                   ctgacAAAGATAGATATCTTTGCTAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTCATAGCAAGGACAt                        >  1:3185978/5‑91 (MQ=25)
                                                                   ctgacAAAGATAGATATCTTTGCTAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTCATAGCAAGGACAt                        >  1:3821150/5‑91 (MQ=25)
                                                                   ctgacAAAGATACATATCTTTGCTAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTCATAGCAAGGACAt                        >  1:2991829/5‑91 (MQ=37)
                                                                   cTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTCATAGCAAGGACAt                        <  1:1031586/91‑1 (MQ=255)
                                                                    tAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTCATAGCAAGGACAtt                       <  1:1310297/91‑1 (MQ=255)
                                                                    tAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTCATAGCAAGGACAtt                       >  1:988217/1‑91 (MQ=255)
                                                                     aGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTCATAGCAAGGACAttt                      <  2:2750850/91‑1 (MQ=255)
                                                                     aGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTCATAGCAAGGACAttt                      <  1:3146948/91‑1 (MQ=255)
                                                                       cAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTCATAGCAAGGACATTTCt                    <  2:692776/91‑1 (MQ=255)
                                                                        aaaGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTCATAGCAAGGACATTTCTa                   >  1:3837639/1‑91 (MQ=255)
                                                                        aaaGATATCTATCTTTGTCAGGGTTCTGCTTTATGAATTATGGGATCGAATCTTCAACTAAAATTACGCCTTCATAGCAAGGACATTTCTa                   <  2:1537159/91‑1 (MQ=255)
                                                                         aaGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTCATAGCAAGGACATTTCTAt                  >  1:3793046/1‑91 (MQ=255)
                                                                             tatctatcTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTCATAGCAAGGACATTTCTATTCaa              <  1:3558778/91‑1 (MQ=255)
                                                                             tatctatcTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTCATAGCAAGGACATTTCTATTCaa              <  1:3735742/91‑1 (MQ=255)
                                                                             tatctatcTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTCATAGCAAGGACATTTCTATTCaa              >  1:1762199/1‑91 (MQ=255)
                                                                             tatctatcTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTCATAGCAAGGACATTTCTATTCaa              >  1:1418947/1‑91 (MQ=255)
                                                                             tagatatctTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTCATAGCAAGGACATTTCTATTCaa              >  1:594457/5‑91 (MQ=255)
                                                                              agatatctTTGCTAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTCATAGCAAGGACATTTCTATTCAAt             >  1:2246434/4‑91 (MQ=255)
                                                                               gatatctTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTCATAGCAAGGACATTTCTATTCAATc            >  2:1577859/3‑91 (MQ=255)
                                                                               gatatctTTGCTAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTCATAGCAAGGACATTTCTATTCAATc            >  1:2230349/3‑91 (MQ=255)
                                                                                ctatctTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTCATAGCAAGGACATTTCTATTCAATCa           >  1:1918613/1‑91 (MQ=255)
                                                                                ctatctTTGTCAGGGGTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTCATAGCAAGGACATTGCTATTCAATCa           >  2:1308650/1‑91 (MQ=255)
                                                                                 tatctTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTCATAGCAAGGACATTTCTATTCAATCAc          >  2:1148592/1‑91 (MQ=255)
                                                                                   tctTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTCATAGCAAGGACATTTCTATTCAATCAccc        <  1:1132418/91‑1 (MQ=255)
                                                                                   tctTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTCATAGCAAGGACATTTCTATTCAATCAccc        >  2:3831985/1‑91 (MQ=255)
                                                                                    ctTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTCATAGCAAGGACATTTCTATTCAATCACCCt       >  2:1644752/1‑91 (MQ=255)
                                                                                      ttGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTCATAGCAAGGACATTTCTATTCAATCACCCttt     >  2:1963471/1‑91 (MQ=255)
                                                                                        gTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTCATAGCAAGGACATTTCTATTCAATCACCCTTTaa   >  2:1943653/1‑91 (MQ=255)
                                                                                        gTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTCATAGCAAGGACATTTCTATTCAATCACCCTTTaa   <  1:3507781/91‑1 (MQ=255)
                                                                                         tCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTCATAGCAAGGACATTTCTATTCAATCACCCTTTAAc  <  1:2155491/91‑1 (MQ=255)
                                                                                         tCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTCATAGCAAGGACATTTCTATTCAATCACCCTTTAAc  >  1:2469904/1‑91 (MQ=255)
                                                                                         tCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTCATAGCAAGGACATTTCTATTCAATCACCCTTTAAc  >  2:2402424/1‑91 (MQ=255)
                                                                                         |                                                                                          
ACAAATCATTTACTTTAAAATGATTAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTCATAGCAAGGACATTTCTATTCAATCACCCTTTAAC  >  USA300TCH1516_ALE/2559761‑2559940

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: