Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 2,559,851 C→T 9.0% intergenic (+45/+78) 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,8510CT 9.0% 580.4 / 34.0 189intergenic (+45/+78)bcr_2/USA300HOU_RS13190Bicyclomycin resistance protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (93/79);  new base T (13/4);  total (106/83)
Fisher's exact test for biased strand distribution p-value = 1.22e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.60e-01

Only 100 of 197 total aligned reads displayed.

ACAAATCATTTACTTTAAAATGATTAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTCATAGCAAGGACATTTCTATTCAATCACCCTTTAACA  >  USA300TCH1516_ALE/2559761‑2559941
                                                                                          |                                                                                          
aCAAATCATTTACTTTAAAATGATTAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTc                                                                                            <  2:339189/91‑1 (MQ=255)
aCAAATCATTTACTTTAAAATGATTAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTc                                                                                            >  2:2180685/1‑91 (MQ=255)
 cAAATCATTTACTTTAAAATGATTAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCa                                                                                           >  2:389742/1‑91 (MQ=255)
  aaaTCATTTACTTTAAAATGATTAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAg                                                                                          <  1:2452438/91‑1 (MQ=255)
  aaaTCATTTACTTTAAAATGATTAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAg                                                                                          >  2:1601247/1‑91 (MQ=255)
   aaTCATTTACTTTAAAATGATTAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTcggg                                                                                         <  1:2478146/91‑4 (MQ=255)
    aTCATTTACTTTAAAATGATTAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAggg                                                                                        >  2:2258036/1‑91 (MQ=255)
      cATTTACTTTAAAATGATTAAAAAGCAACATGTCGCATAACACTTCAACCTAATTAGAACCCTGACAAAGATAGATATCTTTGCTAGGGtt                                                                                      <  1:3370186/91‑1 (MQ=21)
      cATTTACTTTAAAATGATTAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGtt                                                                                      >  1:1794935/1‑91 (MQ=255)
        tttACTTTAAAATGATTAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGttct                                                                                    >  1:846819/1‑91 (MQ=255)
         ttACTTTAAAATGATTAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGttctt                                                                                   >  1:830057/1‑91 (MQ=255)
          tACTTTAAAATGATTAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGttcttc                                                                                  >  1:2089875/1‑91 (MQ=255)
            cTTTAAAATGATTAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTcttctt                                                                                >  2:2699088/1‑91 (MQ=255)
             tttAAAATGATTAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCttt                                                                               >  1:1552095/1‑91 (MQ=255)
             tttAAAATGATTAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCttt                                                                               <  1:2041042/91‑1 (MQ=255)
             tttAAAATGATTAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCttt                                                                               >  1:1200685/1‑91 (MQ=255)
             tttAAAATGATTAAAAAGCAACATGTCGCATAACACTTCAACATAAGTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCttt                                                                               <  2:3170456/91‑1 (MQ=255)
              ttAAAATGATTAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTa                                                                              <  2:3001168/91‑1 (MQ=255)
               tAAAATGATTAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTAt                                                                             >  2:1415794/1‑91 (MQ=255)
                  aaTGATTAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGaa                                                                          >  2:1880978/1‑91 (MQ=255)
                     gATTAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTa                                                                       >  1:1616999/1‑91 (MQ=255)
                     gATTAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTa                                                                       <  1:2843663/91‑1 (MQ=255)
                     gATTAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTa                                                                       <  1:219017/91‑1 (MQ=255)
                      aTTAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTAt                                                                      >  1:1877137/1‑91 (MQ=255)
                      aTTAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTAt                                                                      <  1:2314463/91‑1 (MQ=255)
                      aTTAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTAt                                                                      <  1:1206406/91‑1 (MQ=255)
                       ttAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATg                                                                     <  1:721136/91‑1 (MQ=255)
                        tAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATga                                                                    >  1:2894949/1‑91 (MQ=255)
                          aaaaGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATgaga                                                                  <  2:2493123/91‑1 (MQ=255)
                          aaaaGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATgaga                                                                  <  1:1031882/91‑1 (MQ=255)
                          aaaaGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATgaga                                                                  >  1:1576757/1‑91 (MQ=255)
                            aaGCAACATGTCGCATAACACTTCAACATAATTAGAGCCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATc                                                                <  2:3052474/91‑1 (MQ=255)
                              gCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGa                                                              <  1:1913945/91‑1 (MQ=255)
                                 aCATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATc                                                           >  2:3205472/1‑91 (MQ=255)
                                   aTGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCtt                                                         >  2:2137577/1‑91 (MQ=255)
                                       cGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAAc                                                     <  1:3771826/91‑1 (MQ=255)
                                       cGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAACCTTCAAc                                                     >  2:3450467/1‑91 (MQ=255)
                                         cATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTa                                                   <  2:3734548/91‑1 (MQ=255)
                                         cATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTa                                                   <  1:3458098/91‑1 (MQ=255)
                                          aTAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTaa                                                  >  1:2965607/1‑91 (MQ=255)
                                           tAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTaaa                                                 <  1:355459/91‑1 (MQ=255)
                                           tAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTaaa                                                 <  2:1576757/91‑1 (MQ=255)
                                            aaCACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTaaaa                                                >  2:263851/1‑91 (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)
                                                    aaCATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCt                                        <  1:2561434/91‑1 (MQ=255)
                                                    aaCATAATTAGAACCCTAGCAAAGATATCGATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCt                                        <  2:722813/91‑1 (MQ=255)
                                                     aCATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCtt                                       <  2:2341279/91‑1 (MQ=255)
                                                      cATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTc                                      <  2:1214778/91‑1 (MQ=255)
                                                      cATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTc                                      >  2:3507781/1‑91 (MQ=255)
                                                       aTAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTCa                                     >  1:28132/1‑91 (MQ=255)
                                                       aTAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTCa                                     <  1:3693693/91‑1 (MQ=255)
                                                          aTTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTCATAg                                  >  1:2369220/1‑91 (MQ=255)
                                                           ttAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTCATAGc                                 >  1:2926884/1‑91 (MQ=255)
                                                            tAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTCATAGCa                                <  1:3793868/91‑1 (MQ=255)
                                                            tAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTCATAGCa                                <  1:543150/91‑1 (MQ=255)
                                                             aGAACCCTGACAAAGATAGATATCTTTGCTAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTCATAGCaa                               >  1:3206608/1‑91 (MQ=21)
                                                             aGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTCATAGCaa                               >  2:845037/1‑91 (MQ=255)
                                                              gAACCCTGACAAAGATAGATATCTTTGCTAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTCATAGCAAg                              >  1:188455/1‑91 (MQ=21)
                                                              gAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTCATAGCAAg                              <  1:595875/91‑1 (MQ=255)
                                                              gAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTCATAGCAAg                              <  1:3736695/91‑1 (MQ=255)
                                                               aaCCCTGACAAAGATAGATATCTTTGCTAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTCATAGCAAgg                             >  1:31005/1‑91 (MQ=21)
                                                                aCCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTCATAGCAAGGa                            <  1:2014648/91‑1 (MQ=255)
                                                                 cccTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTCATAGCAAGGAc                           <  1:3264893/91‑1 (MQ=255)
                                                                  ccTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTCATAGCAAGGACa                          >  1:444454/1‑91 (MQ=255)
                                                                  ccTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTCATAGCAAGGACa                          >  1:2179444/1‑91 (MQ=255)
                                                                   ctgacAAAGATAGATATCTTTGCTAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTCATAGCAAGGACAt                         >  1:3185978/5‑91 (MQ=25)
                                                                   ctgacAAAGATAGATATCTTTGCTAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTCATAGCAAGGACAt                         >  1:3821609/5‑91 (MQ=25)
                                                                   ctgacAAAGATACATATCTTTGCTAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTCATAGCAAGGACAt                         >  1:2991829/5‑91 (MQ=37)
                                                                   cTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTCATAGCAAGGACAt                         <  1:1031586/91‑1 (MQ=255)
                                                                    tAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTCATAGCAAGGACAtt                        >  1:988217/1‑91 (MQ=255)
                                                                    tAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTCATAGCAAGGACAtt                        <  1:1310297/91‑1 (MQ=255)
                                                                     aGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTCATAGCAAGGACAttt                       <  1:3146948/91‑1 (MQ=255)
                                                                     aGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTCATAGCAAGGACAttt                       <  2:2750850/91‑1 (MQ=255)
                                                                       cAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTCATAGCAAGGACATTTCt                     <  2:692776/91‑1 (MQ=255)
                                                                         aaGATATCTATCTTTGTCAGGGTTCTTCTTTAAGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTCATAGCAAGGACATTTCTAt                   <  1:2726715/91‑1 (MQ=255)
                                                                             tatctatcTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTCATAGCAAGGACATTTCTATTCaa               >  1:1762199/1‑91 (MQ=255)
                                                                             tatctatcTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTCATAGCAAGGACATTTCTATTCaa               >  1:1418947/1‑91 (MQ=255)
                                                                             tagatatctTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTCATAGCAAGGACATTTCTATTCaa               >  1:594440/5‑91 (MQ=255)
                                                                              atctatctTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTCATAGCAAGGACATTTCTATTCAAt              >  1:753920/1‑91 (MQ=255)
                                                                               gatatctTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTCATAGCAAGGACATTTCTATTCAATc             >  2:1577859/3‑91 (MQ=255)
                                                                               gatatctTTGCTAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTCATAGCAAGGACATTTCTATTCAATc             >  1:2232722/3‑91 (MQ=255)
                                                                                ctatctTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTCATAGCAAGGACATTTCTATTCAATCa            <  1:2888198/91‑1 (MQ=255)
                                                                                ctatctTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTCATAGCAAGGACATTTCTATTCAATCa            >  1:1918613/1‑91 (MQ=255)
                                                                                ctatctTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTCATAGCAAGGACATTTCTATTCAATCa            >  1:582823/1‑91 (MQ=255)
                                                                                ctatctTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCCTCACCTAAAATTACGCCTTCATAGCAAGGACATTTCTATTCAATCa            >  2:1308568/1‑91 (MQ=255)
                                                                                 tatctTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTCATAGCAAGGACATTTCTATTCAATCAc           >  2:1148592/1‑91 (MQ=255)
                                                                                   tctTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTCATAGCAAGGACATTTCTATTCAATCAccc         >  2:3831985/1‑91 (MQ=255)
                                                                                   tctTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTCATAGCAAGGACATTTCTATTCAATCAccc         <  1:1132418/91‑1 (MQ=255)
                                                                                    ctTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTCATAGCAAGGACATTTCTATTCAATCACCCt        >  2:1644752/1‑91 (MQ=255)
                                                                                      ttGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTCATAGCAAGGACATTTCTATTCAATCACCCttt      <  1:2267938/91‑1 (MQ=255)
                                                                                      ttGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTCATAGCAAGGACATTTCTATTCAATCACCCttt      <  1:3584987/91‑1 (MQ=255)
                                                                                      ttGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTCATAGCAAGGACATTTCTATTCAATCACCCttt      >  2:1963471/1‑91 (MQ=255)
                                                                                        gTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTCATAGCAAGGACATTTCTATTCAATCACCCTTTaa    >  2:1943653/1‑91 (MQ=255)
                                                                                         tCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTCATAGCAAGGACATTTCTATTCAATCACCCTTTAAc   >  2:2402424/1‑91 (MQ=255)
                                                                                          cAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTCATAGCAAGGACATTTCTATTCAATCACCCTTTAACa  >  2:2229119/1‑91 (MQ=255)
                                                                                          |                                                                                          
ACAAATCATTTACTTTAAAATGATTAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTCATAGCAAGGACATTTCTATTCAATCACCCTTTAACA  >  USA300TCH1516_ALE/2559761‑2559941

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: