Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 2,559,850 T→C 7.2% 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 7.2% 577.3 / 23.8 189intergenic (+44/+79)bcr_2/USA300HOU_RS13190Bicyclomycin resistance protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (93/82);  new base C (10/4);  total (103/86)
Fisher's exact test for biased strand distribution p-value = 2.66e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.65e-01

Only 100 of 195 total aligned reads displayed.

TACAAATCATTTACTTTAAAATGATTAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTCATAGCAAGGACATTTCTATTCAATCACCCTTTAAC  >  USA300TCH1516_ALE/2559760‑2559940
                                                                                          |                                                                                          
tACAAATCATTTACTTTAAAATGATTAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGt                                                                                            <  1:365829/91‑1 (MQ=255)
  cAAATCATTTACTTTAAAATGATTAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTcc                                                                                          >  2:1939545/1‑90 (MQ=255)
  cAAATCATTTACTTTAAAATGATTAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCa                                                                                          >  2:2180429/1‑91 (MQ=255)
   aaaTCATTTACTTTAAAATGATTAAAAAGCAACATGTCGCATAACACTTCAACATAAGTAGAACCCTAGCAAAGATATCTATCTTTGTCAg                                                                                         <  2:3577280/91‑1 (MQ=255)
    actCATTTCCTTTAAAATGATTAAAAAGCAACATGTCGCATAACACTTCAACATATTTAGAACCCTCGCAAAGCTCTCTATCCTTGTCAgg                                                                                        >  2:2802570/3‑91 (MQ=255)
     aTCATTTACTTTAAAATGATTAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAggg                                                                                       >  1:2786493/1‑91 (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)
            aCTTTAAAATGATTAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTGACAAAGATAGATATCTTTGCTAGGGttctaat                                                                                <  1:1592485/91‑4 (MQ=18)
            aCTTTAAAATGATTAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGttcttgt                                                                                >  2:535930/1‑89 (MQ=255)
              tttAAAATGATTAAAAAGCAACATGTCGCATAACACTTCAACCTAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGCTCTTCttt                                                                              >  2:1595390/1‑91 (MQ=255)
               ttAAAATGATTAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTa                                                                             >  1:3075421/1‑91 (MQ=255)
               ttAAAATGATTAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTa                                                                             <  1:1040377/91‑1 (MQ=255)
                   aaTGATTAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGaa                                                                         <  1:1189116/91‑1 (MQ=255)
                   aaTGATTAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGaa                                                                         <  1:2295551/91‑1 (MQ=255)
                   aaTGATTAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCCAGCAACGACATCCCCCCTTCTCAGGGTTCTTCTCTACGaa                                                                         >  2:3378130/1‑91 (MQ=255)
                       aTTAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTAt                                                                     >  2:1850266/1‑91 (MQ=255)
                           aaaaGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATgaga                                                                 >  1:981755/1‑91 (MQ=255)
                             aaGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATc                                                               >  1:255843/1‑91 (MQ=255)
                                 aaCATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAAt                                                           >  1:1110631/1‑91 (MQ=255)
                                  aCATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATc                                                          >  1:2150917/1‑91 (MQ=255)
                                  aCATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATc                                                          >  2:3164933/1‑91 (MQ=255)
                                   cATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCt                                                         >  2:1489844/1‑91 (MQ=255)
                                   cATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCt                                                         >  1:1166466/1‑91 (MQ=255)
                                    aTGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCtt                                                        >  1:1006852/1‑91 (MQ=255)
                                    aTGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCtt                                                        >  1:1006819/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                                                       <  1:2991315/91‑1 (MQ=255)
                                       tCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCaa                                                     <  1:777359/91‑1 (MQ=255)
                                         gCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCGGGGTTCTTCTTTATGAATTATGAAATCCAATCTTCAACt                                                   <  2:846223/91‑1 (MQ=255)
                                         gCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGCTCGAATCTTCAACt                                                   >  1:371988/1‑91 (MQ=255)
                                         gCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACt                                                   <  2:1638865/91‑1 (MQ=255)
                                         gCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCAAATCTTCAACt                                                   <  2:1742884/91‑1 (MQ=255)
                                          cATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTa                                                  <  1:468501/91‑1 (MQ=255)
                                            tAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTaaa                                                >  1:2417941/1‑91 (MQ=255)
                                             aaCACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTaaaa                                               >  2:740207/1‑91 (MQ=255)
                                             aaCACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTaaaa                                               >  2:2127192/1‑91 (MQ=255)
                                              acacTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAAt                                              <  1:858700/91‑1 (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:821915/91‑1 (MQ=255)
                                                    cAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGcc                                        >  1:1635716/1‑91 (MQ=255)
                                                    cAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGcc                                        >  2:3697883/1‑91 (MQ=255)
                                                     aaCATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCt                                       <  1:3842672/91‑1 (MQ=255)
                                                      aCATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCtt                                      >  1:1652391/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)
                                                           aTTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTCATAg                                 >  1:2629824/1‑91 (MQ=255)
                                                            ttAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGACTCTTCAACTAAAATTACGCCTTCATAGc                                >  2:2606308/1‑91 (MQ=255)
                                                            gaaGAACCCTGACAAAGATAGATATCTTTGCTAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTCATAGc                                >  1:2111366/3‑91 (MQ=18)
                                                             tAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTCATAGCa                               <  2:3006765/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)
                                                               gAACCCTAGCAAAGATATCTATCTTGGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTCATAGCAAg                             <  2:22124/91‑1 (MQ=255)
                                                                aaCCCTAGCAAAGATATCTATCTTTGTCGGGGTTCTTCTTTATGAATTATGAAATCCAATCTTCAACTAAAATTACGACTTCATATCAAgt                            >  2:2306418/1‑90 (MQ=255)
                                                                 aCCCTGACAAAGATAGATATCTTTGCTAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTCATAGCAAGGa                           >  2:478249/1‑91 (MQ=25)
                                                                  cccTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTCATAGCAAGGAc                          <  1:2294864/91‑1 (MQ=255)
                                                                  cccTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTCATAGCAAGGAc                          >  1:556250/1‑91 (MQ=255)
                                                                  cccTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTCATAGCAAGGAc                          <  1:2702765/91‑1 (MQ=255)
                                                                   ccTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTCATAGCAAGGACa                         <  2:2080139/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)
                                                                        cAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTACAATTACGCCTTCATAGCAAGGACATTTCt                    >  2:2754339/1‑91 (MQ=255)
                                                                        cAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTCATAGCAAGGACATTTCt                    <  1:34933/91‑1 (MQ=255)
                                                                        cAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTCATAGCAAGGACATTTCt                    >  2:1553527/1‑91 (MQ=255)
                                                                        cAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTCATAGCAAGGACATTTCt                    <  2:3271148/91‑1 (MQ=255)
                                                                        cAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATAGAATCTTCAACTAAAATTACGCCTTCATACCAAGGACATTTCt                    >  1:3951211/1‑91 (MQ=255)
                                                                        cAAAGATAGATATCTTTGCTAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTCATAGCAAGGACATTTCt                    >  2:1310290/1‑91 (MQ=38)
                                                                          aaGATATCTAGCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTCATAGCAAGGACATTTCTAt                  <  2:148708/91‑1 (MQ=255)
                                                                           aGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTCATAGCAAGGACATTTCTAtt                 >  1:2229136/1‑91 (MQ=255)
                                                                           aGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTCATAGCAAGGACATTTCTAtt                 >  2:3817420/1‑91 (MQ=255)
                                                                             atatCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTCATAGCAAGGACATTTCTATTCa               <  1:3804293/91‑1 (MQ=255)
                                                                             atatCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTCATAGCAAGGACATTTCTATTCa               <  2:2684274/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)
                                                                                   atctTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTCATAGCAAGGACATTTCTATTCAATCAcc         <  2:2329418/91‑1 (MQ=255)
                                                                                   atctTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTCATAGCAAGGACATTTCTATTCAATCAcc         <  1:2228869/91‑1 (MQ=255)
                                                                                   atctTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTCATAGCAAGGACATTTCTATTCAATCAcc         <  2:2874844/91‑1 (MQ=255)
                                                                                    tctTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTCATAGCAAGGACATTTCTATTCAATCAccc        >  1:2226596/1‑91 (MQ=255)
                                                                                     ctTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTCATAGCAAGGACATTTCTATTCAATCACCCt       <  2:332291/91‑1 (MQ=255)
                                                                                     ctTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTCATAGCAAGGACATTTCTATTCAATCACCCt       <  1:1594207/91‑1 (MQ=255)
                                                                                      tttGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTCATAGCAAGGACATTTCTATTCAATCACCCtt      <  1:998248/91‑1 (MQ=255)
                                                                                        tGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTCATAGCAAGGACATTTCTATTCAATCACCCTTTa    >  1:41783/1‑91 (MQ=255)
                                                                                         gTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTCATAGCAAGGACATTTCTATTCAATCACCCTTTaa   >  1:3867086/1‑91 (MQ=255)
                                                                                         gTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTCATAGCAAGGACATTTCTATTCAATCACCCTTTaa   >  2:885273/1‑91 (MQ=255)
                                                                                          tCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTCATAGCAAGGACATTTCTATTCAATCACCCTTTAAc  >  1:2143959/1‑91 (MQ=255)
                                                                                          tCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTCATAGCAAGGACATTTCTATTCAATCACCCTTTAAc  >  2:2495495/1‑91 (MQ=255)
                                                                                          tCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTCATAGCAAGGACATTTCTATTCAATCACCCTTTAAc  <  1:2832019/91‑1 (MQ=255)
                                                                                          |                                                                                          
TACAAATCATTTACTTTAAAATGATTAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTCATAGCAAGGACATTTCTATTCAATCACCCTTTAAC  >  USA300TCH1516_ALE/2559760‑2559940

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: