Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 2,559,841 C→A 15.4% intergenic (+35/+88) 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,8410CA15.4% 552.4 / 72.9 187intergenic (+35/+88)bcr_2/USA300HOU_RS13190Bicyclomycin resistance protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (88/70);  new base A (9/20);  total (97/90)
Fisher's exact test for biased strand distribution p-value = 1.61e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.19e-01

Only 100 of 211 total aligned reads displayed.

TAGTCAGTCTACAAATCATTTACTTTAAAATGATTAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTCATAGCAAGGACATTTCTATTCAATCA  >  USA300TCH1516_ALE/2559751‑2559931
                                                                                          |                                                                                          
   tcagtcTACAAATCATTTACTTTAAAATGATTAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGAtatctat                                                                                         <  2:768200/91‑1 (MQ=255)
    cagtcTACAAATCATTTACTTTAAAATGATTAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGAtatctatc                                                                                        >  1:2341925/1‑91 (MQ=255)
     agtcTACAAATCATTTACTTTAAAATGATTAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTGACAAAGATagatatct                                                                                       <  2:2649706/91‑1 (MQ=38)
     agtcTACAAATCATTTACTTTAAAATGATTAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATatctatct                                                                                       <  2:2608449/91‑1 (MQ=255)
      gtcTACAAATCATTTACTTTAAAATGATTAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCtt                                                                                      >  1:3239042/1‑91 (MQ=255)
       tcTACAAATCATTTACTTTAAAATGATTAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTGACAAAGATAGATATCttt                                                                                     <  1:2733959/91‑1 (MQ=38)
       tcTACAAATCATTTACTTTAAAATGATTAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTGACAAAGATAGATATCttt                                                                                     <  2:139413/91‑1 (MQ=38)
        cTACAAATCATTTACTTTAAAATGATTAAAAAGCAACATGTCGCATAACACTTGAACATAATTAGAACCCTGACAAAGATAGATATCTTTg                                                                                    <  1:3326470/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                                                                                 >  2:2180429/1‑91 (MQ=255)
           cAAATCATTTACTTTAAAATGATTAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCa                                                                                 >  1:3790848/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)
            aaaTCATTTACTTTAAAATGATTAAAAAGCAACATGTCGCATAACACTTCAACATAAGTAGAACCCTAGCAAAGATATCTATCTTTGTCAg                                                                                <  2:3577280/91‑1 (MQ=255)
             actCATTTCCTTTAAAATGATTAAAAAGCAACATGTCGCATAACACTTCAACATATTTAGAACCCTCGCAAAGCTCTCTATCCTTGTCAgg                                                                               >  2:2802570/3‑91 (MQ=255)
              aTCATTTACTTTAAAATGATTAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAggg                                                                              <  2:3829384/91‑1 (MQ=255)
              aTCATTTACTTTAAAATGATTAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAggg                                                                              >  1:2634642/1‑91 (MQ=255)
              aTCATTTACTTTAAAATGATTAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAggg                                                                              >  1:2786493/1‑91 (MQ=255)
                 aTTTACTTTAAAATGATTAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGttc                                                                           >  2:454801/1‑91 (MQ=255)
                   ttACTTTAAAATGATTAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGttctt                                                                         <  2:3197747/91‑1 (MQ=255)
                    tACTTTAAAATGATTAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGttcttc                                                                        >  1:1207429/1‑91 (MQ=255)
                    tACTTTAAAATGATTAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGttcttc                                                                        <  2:1563988/91‑1 (MQ=255)
                     aCTTTAAAATGATTAAAAAGCAACATGTCGCATAACACTTCGACATAACTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTtcttct                                                                       <  2:2164326/91‑1 (MQ=255)
                     aCTTTAAAATGATTAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTGACAAAGATAGATATCTTTGCTAGGGttctaat                                                                       <  1:1593607/91‑4 (MQ=18)
                        ttAAAATGATTAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTa                                                                    >  1:3075421/1‑91 (MQ=255)
                          aaaaTGATTAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATg                                                                  >  2:631870/1‑91 (MQ=255)
                            aaTGATTAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGaa                                                                <  1:2295551/91‑1 (MQ=255)
                            aaTGATTAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGaa                                                                >  2:3377937/1‑91 (MQ=255)
                                aTTAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTAt                                                            >  2:1850266/1‑91 (MQ=255)
                                aTTAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTAt                                                            >  2:567117/1‑91 (MQ=255)
                                 ttAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATg                                                           >  2:1191439/1‑91 (MQ=255)
                                    aaaaGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATgaga                                                        >  1:2891222/1‑91 (MQ=255)
                                    aaaaGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATgaga                                                        >  1:981755/1‑91 (MQ=255)
                                       aGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCg                                                     <  2:1952496/91‑1 (MQ=255)
                                           aCATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATc                                                 >  2:3164933/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)
                                             aTGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCtt                                               >  2:3690206/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)
                                                  gCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACt                                          <  2:3815624/91‑1 (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                                <  2:3188898/91‑1 (MQ=255)
                                                            tCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGc                                <  1:985925/91‑1 (MQ=255)
                                                            tCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGc                                <  1:821915/91‑1 (MQ=255)
                                                            tCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGc                                <  1:3889837/91‑1 (MQ=255)
                                                             cAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGcc                               >  2:3697883/1‑91 (MQ=255)
                                                             cAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGcc                               >  1:1635716/1‑91 (MQ=255)
                                                               aCATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCtt                             >  2:882000/1‑91 (MQ=255)
                                                               aCATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCtt                             >  1:1652391/1‑91 (MQ=255)
                                                               aCATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCtt                             <  1:1098600/91‑1 (MQ=255)
                                                                cATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCGACTAAAATTACGCCTTc                            >  1:2586146/1‑91 (MQ=255)
                                                                 aTAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAGATTACGCCTTCa                           >  2:3290953/1‑91 (MQ=255)
                                                                 aTAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAACATTACGCCTTCa                           >  2:3290558/1‑91 (MQ=255)
                                                                 aTAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTCa                           >  2:508402/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)
                                                                     gaaGAACCCTGACAAAGATAGATATCTTTGCTAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTCATAGc                       >  1:2111366/3‑91 (MQ=18)
                                                                      tAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTCATAGCa                      >  1:1234304/1‑91 (MQ=255)
                                                                      tAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTCATAGCa                      <  2:3006765/91‑1 (MQ=255)
                                                                      tAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTCATAGCa                      <  1:2533512/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)
                                                                          aCCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTCATAGCAAGGa                  <  2:340117/91‑1 (MQ=255)
                                                                           cccTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTCATAGCAAGGAc                 <  1:1180638/91‑1 (MQ=255)
                                                                            ccTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTCATAGCAAGGACa                >  1:1772121/1‑91 (MQ=255)
                                                                            ccTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTAATAGCAAGGACa                <  1:3024613/91‑1 (MQ=255)
                                                                              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           <  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)
                                                                                      atatCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTCATAGCAAGGACATTTCTATTCa      <  1:2169632/91‑1 (MQ=255)
                                                                                      atatCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTCATAGCAAGGACATTTCTATTCa      <  1:3804710/91‑1 (MQ=255)
                                                                                      atatCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTCATAGCAAGGACATTTCTATTCa      <  2:2298647/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)
                                                                                          ctatctTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTCATAGCAAGGACATTTCTATTCAATCa  <  2:1938730/91‑1 (MQ=255)
                                                                                          |                                                                                          
TAGTCAGTCTACAAATCATTTACTTTAAAATGATTAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTCATAGCAAGGACATTTCTATTCAATCA  >  USA300TCH1516_ALE/2559751‑2559931

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