Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 2,559,851 C→T 6.7% 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 6.7% 615.8 / 20.8 192intergenic (+45/+78)bcr_2/USA300HOU_RS13190Bicyclomycin resistance protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (96/83);  new base T (9/4);  total (105/87)
Fisher's exact test for biased strand distribution p-value = 3.89e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.13e-01

Only 100 of 196 total aligned reads displayed.

CAAATCATTTACTTTAAAATGATTAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTCATAGCAAGGACATTTCTATTCAATCACCCTTTAACA  >  USA300TCH1516_ALE/2559762‑2559941
                                                                                         |                                                                                          
cAAATCATTTACTTTAAAATGATTAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCa                                                                                           >  1:3790848/1‑91 (MQ=255)
 aaaTCATTTACTTTAAAATGATTAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAg                                                                                          >  2:2843170/1‑91 (MQ=255)
 aaaTCATTTACTTTAAAATGATTAAAAAGCAACATGTCGCATAACACTTCAACATAAGTAGAACCCTAGCAAAGATATCTATCTTTGTCAg                                                                                          <  2:3577280/91‑1 (MQ=255)
   aTCATTTACTTTAAAATGATTAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAggg                                                                                        >  1:2786493/1‑91 (MQ=255)
   aTCATTTACTTTAAAATGATTAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAggg                                                                                        <  1:2341136/91‑1 (MQ=255)
    tCATTTACTTTAAAATGATTAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGt                                                                                       <  2:41694/91‑1 (MQ=255)
     cATTTACTTTAAAATGATTAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTGACAAAGATAGATATCTTTGCTAGGGtt                                                                                      <  2:608065/91‑1 (MQ=25)
        ttACTTTAAAATGATTAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGttctt                                                                                   <  2:3197747/91‑1 (MQ=255)
         tACTTTAAAATGATTAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGttcttc                                                                                  <  2:1563988/91‑1 (MQ=255)
          aCTTTAAAATGATTAAAAAGCAACATGTCGCATAACACTTCGACATAACTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTtcttct                                                                                 <  2:2164326/91‑1 (MQ=255)
          aCTTTAAAATGATTAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTGACAAAGATAGATATCTTTGCTAGGGttctaat                                                                                 <  1:1592485/91‑4 (MQ=18)
          aCTTTAAAATGATTAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGttcttgt                                                                                 >  2:535930/1‑89 (MQ=255)
             ttAAAATGATTAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTa                                                                              <  1:1040377/91‑1 (MQ=255)
             ttAAAATGATTAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTa                                                                              >  1:3075421/1‑91 (MQ=255)
                 aaTGATTAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGaa                                                                          <  1:2295187/91‑1 (MQ=255)
                 aaTGATTAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGaa                                                                          <  1:1189116/91‑1 (MQ=255)
                 aaTGATTAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCCAGCAACGACATCCCCCCTTCTCAGGGTTCTTCTCTACGaa                                                                          >  2:3378130/1‑91 (MQ=255)
                   tGATTAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAAtt                                                                        >  1:424568/1‑91 (MQ=255)
                         aaaaGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATgaga                                                                  >  1:981755/1‑91 (MQ=255)
                            aGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCg                                                               <  2:1952496/91‑1 (MQ=255)
                               aaCATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAAt                                                            >  1:1110631/1‑91 (MQ=255)
                                aCATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATc                                                           >  1:2150917/1‑91 (MQ=255)
                                aCATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATc                                                           >  2:3164933/1‑91 (MQ=255)
                                 cATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCt                                                          >  1:1166466/1‑91 (MQ=255)
                                 cATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCt                                                          >  2:1489844/1‑91 (MQ=255)
                                 cATGTCGCATAACACTTCAACATAATTAGAACCCCAGCAAAGATATCTATCCCCGCCAGCGTCCTCCCTAACGAATTACGAGATCGCATCt                                                          >  2:1634530/1‑91 (MQ=255)
                                  aTGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCtt                                                         >  1:402255/1‑91 (MQ=255)
                                  aTGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCtt                                                         >  1:1006819/1‑91 (MQ=255)
                                  aTGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCtt                                                         >  1:1006852/1‑91 (MQ=255)
                                  aTGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCtt                                                         >  2:3690206/1‑91 (MQ=255)
                                  aTGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCtt                                                         >  1:1385451/1‑91 (MQ=255)
                                   tGTCGCATAACACTTCAACATAATTCGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTc                                                        <  1:3673759/91‑1 (MQ=255)
                                   tGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTc                                                        <  2:2414889/91‑1 (MQ=255)
                                   tGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTc                                                        <  1:2991315/91‑1 (MQ=255)
                                       gCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCGGGGTTCTTCTTTATGAATTATGAAATCCAATCTTCAACt                                                    <  2:846223/91‑1 (MQ=255)
                                           aaCACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTaaaa                                                >  2:2127180/1‑91 (MQ=255)
                                           aaCACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTaaaa                                                >  2:740207/1‑91 (MQ=255)
                                            acacTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAAt                                               <  1:858700/91‑1 (MQ=255)
                                             cacTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAAtt                                              <  2:275805/91‑1 (MQ=255)
                                              acTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTa                                             >  2:2970623/1‑91 (MQ=255)
                                               cTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTAc                                            >  1:1321265/1‑91 (MQ=255)
                                                ttCAACATAATTAGAACCCTAGCAACGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACg                                           >  2:3732025/1‑91 (MQ=255)
                                                 tCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGc                                          <  1:985925/91‑1 (MQ=255)
                                                 tCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGc                                          <  2:3187123/91‑1 (MQ=255)
                                                 tCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGc                                          <  1:821945/91‑1 (MQ=255)
                                                 tCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGc                                          <  1:3889837/91‑1 (MQ=255)
                                                  cAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGcc                                         >  1:1635716/1‑91 (MQ=255)
                                                  cAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGcc                                         >  2:3697883/1‑91 (MQ=255)
                                                    aCATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCtt                                       >  1:1652391/1‑91 (MQ=255)
                                                     cATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCGACTAAAATTACGCCTTc                                      >  1:2586146/1‑91 (MQ=255)
                                                      aTAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAACATTACGCCTTCa                                     >  2:3290822/1‑91 (MQ=255)
                                                        aaTTAGAACCCTAGCACAGATATCTCTCTTCGTCAGGGTTCTTCTTTATGAATCATGAGACCGAATCTCCAACTAAACTTACGCCTTCATa                                   >  2:2595401/1‑91 (MQ=255)
                                                         aTTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTCATAg                                  >  1:2629824/1‑91 (MQ=255)
                                                          ttAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGACTCTTCAACTAAAATTACGCCTTCATAGc                                 >  2:2606308/1‑91 (MQ=255)
                                                          gaaGAACCCTGACAAAGATAGATATCTTTGCTAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTCATAGc                                 >  1:2111366/3‑91 (MQ=18)
                                                           tAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTCATAGCa                                <  2:3006748/91‑1 (MQ=255)
                                                           tAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTCATAGCa                                <  1:2533512/91‑1 (MQ=255)
                                                           tAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTCATAGCa                                <  1:3240269/91‑1 (MQ=255)
                                                           tAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTCATAGCa                                >  1:1234304/1‑91 (MQ=255)
                                                             gAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTCATAGCAAg                              >  1:3533859/1‑91 (MQ=255)
                                                             gAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTCATAGCAAg                              <  1:3378388/91‑1 (MQ=255)
                                                             gAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTCATAGCAAg                              <  2:21837/91‑1 (MQ=255)
                                                              aaCCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTCATAGCAAgg                             >  2:1710924/1‑91 (MQ=255)
                                                               aCCCTGACAAAGATAGATATCTTTGCTAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTCATAGCAAGGa                            >  2:478249/1‑91 (MQ=25)
                                                                cccTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTCATAGCAAGGAc                           <  1:2702765/91‑1 (MQ=255)
                                                                 ccTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTCATAGCAAGGACa                          >  1:1772121/1‑91 (MQ=255)
                                                                 ccTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTAATAGCAAGGACa                          <  1:3024613/91‑1 (MQ=255)
                                                                    gacAAAGATAGATATCTTTGCTAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTCATAGCAAGGACAttt                       >  1:3587081/3‑91 (MQ=37)
                                                                     acAAAGATAGATATCTTTGCTAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTCATAGCAAGGACATTTc                      >  1:2062408/2‑91 (MQ=37)
                                                                      cAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTCATAGCAAGGACATTTCt                     <  2:3271148/91‑1 (MQ=255)
                                                                      cAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTCATAGCAAGGACATTTCt                     <  1:3849249/91‑1 (MQ=255)
                                                                      cAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTCATAGCAAGGACATTTCt                     >  2:1553527/1‑91 (MQ=255)
                                                                      cAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTCATAGCAAGGACATTTCt                     <  1:34933/91‑1 (MQ=255)
                                                                      cAAAGATAGATATCTTTGCTAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTCATAGCAAGGACATTTCt                     >  2:1310290/1‑91 (MQ=38)
                                                                       aaaGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTCATAGCAAGGACATTTCTa                    >  1:78889/1‑91 (MQ=255)
                                                                        aaGATATCTAGCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTCATAGCAAGGACATTTCTAt                   <  2:148708/91‑1 (MQ=255)
                                                                         aGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTCATAGCAAGGACATTTCTAtt                  >  2:3817420/1‑91 (MQ=255)
                                                                           atatCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTCATAGCAAGGACATTTCTATTCa                <  2:2684274/91‑1 (MQ=255)
                                                                           atatCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTCATAGCAAGGACATTTCTATTCa                >  1:643930/1‑91 (MQ=255)
                                                                           atatCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTCATAGCAAGGACATTTCTATTCa                <  2:2298647/91‑1 (MQ=255)
                                                                           atatCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTCATAGCAAGGACATTTCTATTCa                <  1:3804710/91‑1 (MQ=255)
                                                                           atatCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTCATAGCAAGGACATTTCTATTCa                <  1:2169632/91‑1 (MQ=255)
                                                                            tatctatcTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTCATAGCAAGGACATTTCTATTCaa               >  2:113349/1‑91 (MQ=255)
                                                                             atctatctTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCCTCAACTAAAATTACGCCTTCATAGCAAGGACATTTCTATTCAAt              >  1:3138451/1‑91 (MQ=255)
                                                                               ctatctTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTCATAGCAAGGACATTTCTATTCAATCa            <  2:1938730/91‑1 (MQ=255)
                                                                               ctatctGTGTCAGGGCTCTTCGTTATGAATTATGAGATCGAATCTTCAACTACTATGACGGCTTCGGAGCAAGGACATTTCTATGCAATCa            >  2:1775255/1‑91 (MQ=255)
                                                                                 atctTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTCATAGCAAGGACATTTCTATTCAATCAcc          <  1:2229153/91‑1 (MQ=255)
                                                                                 atctTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTCATAGCAAGGACATTTCTATTCAATCAcc          <  1:2228869/91‑1 (MQ=255)
                                                                                 atctTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTCATAGCAAGGACATTTCTATTCAATCAcc          <  2:2329282/91‑1 (MQ=255)
                                                                                  tctTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTCATAGCAAGGACATTTCTATTCAATCAccc         >  1:2226596/1‑91 (MQ=255)
                                                                                  tctTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTCATAGCAAGGACATTTCTATTCAATCAccc         >  2:274534/1‑91 (MQ=255)
                                                                                   ctTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTCATAGCAAGGACATTTCTATTCAATCACCCt        <  2:332291/91‑1 (MQ=255)
                                                                                   ctTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTCATAGCAAGGACATTTCTATTCAATCACCCt        <  1:1594207/91‑1 (MQ=255)
                                                                                   ctTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTCATAGCAAGGACATTTCTATTCAATCACCCt        <  1:499584/91‑1 (MQ=255)
                                                                                    tttGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTCATAGCAAGGACATTTCTATTCAATCACCCtt       <  1:998248/91‑1 (MQ=255)
                                                                                       gTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTCATAGCAAGGACATTTCTATTCAATCACCCTTTaa    >  2:885273/1‑91 (MQ=255)
                                                                                        tCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTCATAGCAAGGACATTTCTATTCAATCACCCTTTAAc   >  2:2495495/1‑91 (MQ=255)
                                                                                        tCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTCATAGCAAGGACATTTCTATTCAATCACCCTTTAAc   >  1:2143959/1‑91 (MQ=255)
                                                                                        tCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTCATAGCAAGGACATTTCTATTCAATCACCCTTTAAc   <  1:2832019/91‑1 (MQ=255)
                                                                                         cAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTCATAGCAAGGACATTTCTATTCAATCACCCTTTAACa  >  1:447788/1‑91 (MQ=255)
                                                                                         |                                                                                          
CAAATCATTTACTTTAAAATGATTAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTCATAGCAAGGACATTTCTATTCAATCACCCTTTAACA  >  USA300TCH1516_ALE/2559762‑2559941

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: