Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 2,559,840 T→G 16.5% intergenic (+34/+89) 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,8400TG16.5% 460.1 / 81.7 180intergenic (+34/+89)bcr_2/USA300HOU_RS13190Bicyclomycin resistance protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (87/63);  new base G (9/21);  total (96/84)
Fisher's exact test for biased strand distribution p-value = 8.36e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.47e-01

Only 100 of 211 total aligned reads displayed.

TTAGTCAGTCTACAAATCATTTACTTTAAAATGATTAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTCATAGCAAGGACATTTCTATTCAAT  >  USA300TCH1516_ALE/2559750‑2559929
                                                                                          |                                                                                         
ttAGTCAGTCTACAAATCATTTACTTTAAAATGATTAAAACGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGatat                                                                                           <  2:3157389/91‑1 (MQ=255)
ttAGTCAGTCTACAAATCATTTACTTTAAAATGATTAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGatat                                                                                           <  2:689872/91‑1 (MQ=255)
  agtcagtcTACAAATCATTTACTTTAAAATGATTAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGAtatct                                                                                         >  2:3371902/1‑91 (MQ=255)
   gtcagtcTACAAATCATTTACTTTAAAATGATTAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGAtatctc                                                                                        >  2:3193749/1‑90 (MQ=255)
    tcagtcTACAAATCATTTACTTTAAAATGATTAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGAtatctat                                                                                       <  2:768200/91‑1 (MQ=255)
     cagtcTACAAATCATTTACTTTAAAATGATTAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGAtatctatc                                                                                      >  1:2341925/1‑91 (MQ=255)
      agtcTACAAATCATTTACTTTAAAATGATTAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTGACAAAGATagatatct                                                                                     <  2:2649706/91‑1 (MQ=38)
      agtcTACAAATCATTTACTTTAAAATGATTAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATagatatct                                                                                     <  2:1878430/91‑1 (MQ=255)
       gtcTACAAATCATTTACTTTAAAATGATTAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCtt                                                                                    >  1:3239042/1‑91 (MQ=255)
        tcTACAAATCATTTACTTTAAAATGATTAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTGACAAAGATAGATATCttt                                                                                   <  2:139413/91‑1 (MQ=38)
        tcTACAAATCATTTACTTTAAAATGATTAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTGACAAAGATAGATATCttt                                                                                   <  1:2733959/91‑1 (MQ=38)
         cTACAAATCATTTACTTTAAAATGATTAAAAAGCAACATGTCGCATAACACTTGAACATAATTAGAACCCTGACAAAGATAGATATCTTTg                                                                                  <  1:3326470/91‑1 (MQ=37)
         cTACAAATCATTTACTGTAAAATGATTAAAAAGCAACATATCGCATAACACTTCAACATAATTAGAACCCTGACAAAGATAGATATCTTTg                                                                                  <  2:346135/91‑1 (MQ=37)
          tACAAATCATTTACTTTAAAATGATTAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTGACAAAGATAGATATCTTTgc                                                                                 <  2:3147086/91‑2 (MQ=37)
          tACAAATCATTTACTTTAAAATGATTAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGt                                                                                 <  1:365674/91‑1 (MQ=255)
            cAAATCATTTACTTTAAAATGATTAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTGACAAAGATAGATATCTTTgcta                                                                               <  1:1499121/91‑4 (MQ=37)
            cAAATCATTTACTTTAAAATGATTAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCa                                                                               >  1:3790848/1‑91 (MQ=255)
            cAAATCATTTACTTTAAAATGATTAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCa                                                                               >  2:2180429/1‑91 (MQ=255)
            cAAAGCAGTGACTGTAAGATGATTAAAAAGCAACATGTCGCATAACGCGTCAACATAATTAGAACCCTGACAAAGATAGATATCTTTgcta                                                                               <  2:1555126/91‑4 (MQ=16)
             aaaTCATTTACTTTAAAATGATTAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTGACAAAGATAGATATCTTTgctag                                                                              <  2:1336172/91‑5 (MQ=37)
             aaaTCATTTACTTTAAAATGATTAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAg                                                                              >  2:2843170/1‑91 (MQ=255)
             aaaTCATTTACTTTAAAATGATTAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAg                                                                              >  1:3493715/1‑91 (MQ=255)
              actCATTTCCTTTAAAATGATTAAAAAGCAACATGTCGCATAACACTTCAACATATTTAGAACCCTCGCAAAGCTCTCTATCCTTGTCAgg                                                                             >  2:2802570/3‑91 (MQ=255)
               aTCATTTACTTTAAAATGATTAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAggg                                                                            >  1:2786493/1‑91 (MQ=255)
               aTCATTTACTTTAAAATGATTAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAggg                                                                            >  1:2634642/1‑91 (MQ=255)
                tCATTTACTTTAAAATGATTAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGt                                                                           <  2:41694/91‑1 (MQ=255)
                 cATTTACTTTAAAATGATTAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTGACAAAGATAGATATCTTTGCTAGGGtt                                                                          <  2:608065/91‑1 (MQ=25)
                     tACTTTAAAATGATTAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGttcttc                                                                      <  2:1563988/91‑1 (MQ=255)
                     tACTTTAAAATGATTAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGttcttc                                                                      >  1:1207429/1‑91 (MQ=255)
                      aCTTTAAAATGATTAAAAAGCAACATGTCGCATAACACTTCGACATAACTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTtcttct                                                                     <  2:2164326/91‑1 (MQ=255)
                      aCTTTAAAATGATTAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTGACAAAGATAGATATCTTTGCTAGGGttctaat                                                                     <  1:1593607/91‑4 (MQ=18)
                      aCTTTAAAATGATTAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGttcttgt                                                                     >  2:535930/1‑89 (MQ=255)
                         ttAAAATGATTATATAGCAACATGTCGCATAACACTTCAACCTAATTAGAACCCTATCAAAGATATCTATCTTTGTCAGGGTTCTTCtttg                                                                  <  2:3589337/91‑2 (MQ=255)
                         ttAAAATGATTAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTa                                                                  >  1:3075421/1‑91 (MQ=255)
                             aaTGATTAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGaa                                                              <  1:2295551/91‑1 (MQ=255)
                             aaTGATTAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCCAGCAACGACATCCCCCCTTCTCAGGGTTCTTCTCTACGaa                                                              >  2:3378130/1‑91 (MQ=255)
                                  ttAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATg                                                         >  2:1191439/1‑91 (MQ=255)
                                     aaaaGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATgaga                                                      >  1:981755/1‑91 (MQ=255)
                                     aaaaGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATgaga                                                      >  1:2891222/1‑91 (MQ=255)
                                        aGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCg                                                   <  2:1952496/91‑1 (MQ=255)
                                        aGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCg                                                   <  2:3713624/91‑1 (MQ=255)
                                            aCATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATc                                               >  2:3167805/1‑91 (MQ=255)
                                            aCATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATc                                               >  1:2150917/1‑91 (MQ=255)
                                             cATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCt                                              >  2:1489844/1‑91 (MQ=255)
                                             cATGTCGCATAACACTTCAACATAATTAGAACCCCAGCAAAGATATCTATCCCCGCCAGCGTCCTCCCTAACGAATTACGAGATCGCATCt                                              >  2:1634530/1‑91 (MQ=255)
                                              aTGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCtt                                             >  1:1006819/1‑91 (MQ=255)
                                              aTGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCtt                                             >  1:402255/1‑91 (MQ=255)
                                              aTGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCtt                                             >  1:1385451/1‑91 (MQ=255)
                                               tGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTc                                            <  2:2414889/91‑1 (MQ=255)
                                               tGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTc                                            <  1:2991315/91‑1 (MQ=255)
                                                 tCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCaa                                          <  1:777359/91‑1 (MQ=255)
                                                   gCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGCTCGAATCTTCAACt                                        >  1:371988/1‑91 (MQ=255)
                                                   gCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCAAATCTTCAACt                                        <  2:1742884/91‑1 (MQ=255)
                                                    cATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTa                                       <  1:468501/91‑1 (MQ=255)
                                                       aaCACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTaaaa                                    >  2:2127180/1‑91 (MQ=255)
                                                        acacTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAAt                                   <  1:858700/91‑1 (MQ=255)
                                                           cTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTAc                                >  1:1321265/1‑91 (MQ=255)
                                                             tCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGc                              <  1:821915/91‑1 (MQ=255)
                                                             tCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGc                              <  1:985925/91‑1 (MQ=255)
                                                             tCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGc                              <  1:3889837/91‑1 (MQ=255)
                                                              cAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGcc                             >  1:1635716/1‑91 (MQ=255)
                                                                aCATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCtt                           >  2:882000/1‑91 (MQ=255)
                                                                aCATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCtt                           <  1:1098600/91‑1 (MQ=255)
                                                                aCATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCtt                           >  1:1652391/1‑91 (MQ=255)
                                                                 cATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCGACTAAAATTACGCCTTc                          >  1:2586146/1‑91 (MQ=255)
                                                                  aTAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAACATTACGCCTTCa                         >  2:3290822/1‑91 (MQ=255)
                                                                  aTAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTCa                         <  2:2988616/91‑1 (MQ=255)
                                                                    aaTTAGAACCCTAGCACAGATATCTCTCTTCGTCAGGGTTCTTCTTTATGAATCATGAGACCGAATCTCCAACTAAACTTACGCCTTCATa                       >  2:2595401/1‑91 (MQ=255)
                                                                    aaTTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTCATa                       >  1:1122530/1‑91 (MQ=255)
                                                                     aTTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGATCTTCTTTATGAACTATGAGATCGAATGTTCAACTAAAATTACGCCTCCATAg                      >  1:3673613/1‑91 (MQ=255)
                                                                      ttAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGACTCTTCAACTAAAATTACGCCTTCATAGc                     >  2:2606308/1‑91 (MQ=255)
                                                                      gaaGAACCCTGACAAAGATAGATATCTTTGCTAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTCATAGc                     >  1:2111366/3‑91 (MQ=18)
                                                                       tAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTCATAGCa                    >  1:1234304/1‑91 (MQ=255)
                                                                       tAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTCATAGCa                    <  1:2533512/91‑1 (MQ=255)
                                                                       tAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTCATAGCa                    <  2:3006765/91‑1 (MQ=255)
                                                                        aGAACCCTGACAAAGATAGATATCTTTGCTAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTCATAGCaa                   >  2:83731/1‑91 (MQ=21)
                                                                         gAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTCATAGCAAg                  <  1:3378397/91‑1 (MQ=255)
                                                                         gAACCCTAGCAAAGATATCTATCTTGGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTCATAGCAAg                  <  2:22124/91‑1 (MQ=255)
                                                                          aaCCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTCATAGCAAgg                 >  2:1710924/1‑91 (MQ=255)
                                                                           aCCCTGACAAAGATAGATATCTTTGCTAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTCATAGCAAGGa                >  2:478249/1‑91 (MQ=25)
                                                                            cccTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTCATAGCAAGGAc               <  1:1180638/91‑1 (MQ=255)
                                                                             ccTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTCATAGCAAGGACa              >  1:1772121/1‑91 (MQ=255)
                                                                             ccTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTAATAGCAAGGACa              <  1:3024613/91‑1 (MQ=255)
                                                                               tgacAAAGATAGATATCTTTGCTAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTCATAGCAAGGACAtt            >  2:3695192/4‑91 (MQ=37)
                                                                               tAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTCATAGCAAGGACAtt            <  1:1020172/91‑1 (MQ=255)
                                                                                gacAAAGATAGATATCTTTGCTAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTCATAGCAAGGACAttt           >  1:3587080/3‑91 (MQ=37)
                                                                                 acAAAGATAGATATCTTTGCTAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTCATAGCAAGGACATTTc          >  1:2062408/2‑91 (MQ=37)
                                                                                  cAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTACAATTACGCCTTCATAGCAAGGACATTTCt         >  2:2754339/1‑91 (MQ=255)
                                                                                  cAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTCATAGCAAGGACATTTCt         <  2:3271148/91‑1 (MQ=255)
                                                                                  cAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTCATAGCAAGGACATTTCt         <  1:3849249/91‑1 (MQ=255)
                                                                                  cAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTCATAGCAAGGACATTTCt         >  2:1553527/1‑91 (MQ=255)
                                                                                    aaGATATCTAGCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTCATAGCAAGGACATTTCTAt       <  2:148708/91‑1 (MQ=255)
                                                                                     aGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTCATAGCAAGGACATTTCTAtt      >  1:2229136/1‑91 (MQ=255)
                                                                                     aGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTCATAGCAAGGACATTTCTAtt      >  2:3817420/1‑91 (MQ=255)
                                                                                       atatCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTCATAGCAAGGACATTTCTATTCa    <  1:2169632/91‑1 (MQ=255)
                                                                                       atatCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTCATAGCAAGGACATTTCTATTCa    <  2:2298647/91‑1 (MQ=255)
                                                                                       atatCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTCATAGCAAGGACATTTCTATTCa    <  1:3804710/91‑1 (MQ=255)
                                                                                        tatctatcTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTCATAGCAAGGACATTTCTATTCaa   >  2:113349/1‑91 (MQ=255)
                                                                                         atctatctTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTCATAGCAAGGACATTTCTATTCAAt  >  1:3272449/1‑91 (MQ=255)
                                                                                         atctatctTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCCTCAACTAAAATTACGCCTTCATAGCAAGGACATTTCTATTCAAt  >  1:3138451/1‑91 (MQ=255)
                                                                                          |                                                                                         
TTAGTCAGTCTACAAATCATTTACTTTAAAATGATTAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTCATAGCAAGGACATTTCTATTCAAT  >  USA300TCH1516_ALE/2559750‑2559929

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