Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 2,559,840 T→G 7.1% 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,8400TG 7.1% 381.9 / 16.5 112intergenic (+34/+89)bcr_2/USA300HOU_RS13190Bicyclomycin resistance protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (56/48);  new base G (4/4);  total (60/52)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.78e-01

Only 100 of 126 total aligned reads displayed.

TTAGTCAGTCTACAAATCATTTACTTTAAAATGATTAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTCATAGCAAGGACATTTCTATTCAA  >  USA300TCH1516_ALE/2559750‑2559928
                                                                                          |                                                                                        
ttAGTCAGTCTACAAATCATTTACTTTAAAATGATTAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGatat                                                                                          >  1:166407/1‑91 (MQ=255)
   gtcagtcTACAAATCATTTACTTTAAAATGATTAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGAtatcta                                                                                       <  1:1903672/91‑1 (MQ=255)
   gtcagtcTACAAATCATTTACTTTAAAATGATTAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGAtatcta                                                                                       <  1:1903960/91‑1 (MQ=255)
    tcagtcTACAAATCATTTACTTTAAAATGATTAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGAtatctat                                                                                      >  1:3876123/1‑91 (MQ=255)
    tcagtcTACAAATCATTTACTTTAAAATGATTAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGAtatctat                                                                                      >  1:1375195/1‑91 (MQ=255)
      agtcTACAAATCATTTACTTTAAAATGATTAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATatctatct                                                                                    <  1:1237257/91‑1 (MQ=255)
       gtcTACAAATCATTTACTTTAAAATGATTAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCtt                                                                                   >  1:1308850/1‑91 (MQ=255)
         cTACAAATCATTTACTTTAAAATGATTAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTGACAAAGATAGATATCTTTg                                                                                 <  1:2852803/91‑1 (MQ=38)
         cTACAAATCATTTACTTTAAAATGATTAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTGACAAAGATAGATATCTTTg                                                                                 <  1:109030/91‑1 (MQ=38)
           aCAAATCATTTACTTTAAAATGATTAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTc                                                                               <  1:2595842/91‑1 (MQ=255)
            cAAATCATTTACTTTAAAATGATTAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCa                                                                              >  1:2855206/1‑91 (MQ=255)
            cAAATCATTTACTTTAAAATGATTAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCa                                                                              >  1:2588901/1‑91 (MQ=255)
             aaaTCATTTACTTTAAAATGATTAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAg                                                                             <  1:1321486/91‑1 (MQ=255)
               aTCATTTACTTTAAAATGATTAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAggg                                                                           <  1:3731884/91‑1 (MQ=255)
               aTCATTTACTTTAAAATGATTAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAggg                                                                           <  1:1717337/91‑1 (MQ=255)
                tCATTTACTTTAAAATGATTAAAAAGCAACATGTCGCGTAACACTTCAACACAATTAGAACCCTAACAAAGATATCTATCTTTGTCAGGGt                                                                          >  1:3862826/1‑91 (MQ=255)
                tCATTTACTTTAAAATGATTAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGt                                                                          <  1:2317558/91‑1 (MQ=255)
                tCATTTACTTTAAAATGATTAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGt                                                                          <  1:1701050/91‑1 (MQ=255)
                 cATTTACTTTAAAATGATTAAAAAGAAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGtt                                                                         <  1:727817/91‑1 (MQ=255)
                  aTTTACTTTAAAATGATTAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGttc                                                                        >  1:1448903/1‑91 (MQ=255)
                   tttACTTTAAAATGATTAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGttct                                                                       >  1:885390/1‑91 (MQ=255)
                   tttACTTTAAAATGATTAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGttct                                                                       <  1:1873387/91‑1 (MQ=255)
                    ttACTTTAAAATGATTAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGttctt                                                                      >  1:1330796/1‑91 (MQ=255)
                     tACTTTAAAATGATTAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGttcttc                                                                     >  1:1400304/1‑91 (MQ=255)
                      aCTTTAAAATGATTAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTtcttct                                                                    >  1:379714/1‑91 (MQ=255)
                      aCTTTAAAATGATTAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTtcttct                                                                    >  1:1622759/1‑91 (MQ=255)
                       cTTTAAAATGATTAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTcttctt                                                                   >  1:3761007/1‑91 (MQ=255)
                        tttAAAATGATTAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCttt                                                                  <  1:2537752/91‑1 (MQ=255)
                        tttAAAATGATTAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCttt                                                                  <  1:2587522/91‑1 (MQ=255)
                         ttAAAATGATTAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTa                                                                 <  1:53180/91‑1 (MQ=255)
                         ttAAAATGATTAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTa                                                                 <  1:53181/91‑1 (MQ=255)
                           aaaaTGATTAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATg                                                               >  1:169639/1‑91 (MQ=255)
                           aaaaTGAGTAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATg                                                               <  1:2127990/91‑1 (MQ=255)
                                 aTTAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTAt                                                         >  1:662875/1‑91 (MQ=255)
                                 aTTAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTAt                                                         >  1:248902/1‑91 (MQ=255)
                                 aTTAAAAAGCAACATGTAGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTAt                                                         <  1:2186037/91‑1 (MQ=255)
                                  ttAAAAAGCAACATGTCGCATAACACTTCAACTTAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATg                                                        <  1:3548021/91‑1 (MQ=255)
                                  ttAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATg                                                        <  1:2020683/91‑1 (MQ=255)
                                  ttAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATg                                                        <  1:3767537/91‑1 (MQ=255)
                                   tAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATga                                                       <  1:1615854/91‑1 (MQ=255)
                                   tAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATga                                                       <  1:1615785/91‑1 (MQ=255)
                                     aaaaGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATgaga                                                     >  1:646651/1‑91 (MQ=255)
                                     aaaaGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATgaga                                                     >  1:979312/1‑91 (MQ=255)
                                     aaaaGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATgaga                                                     >  1:1616026/1‑91 (MQ=255)
                                        aGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCg                                                  <  1:1053193/91‑1 (MQ=255)
                                        aGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCg                                                  <  1:1838645/91‑1 (MQ=255)
                                        aGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCg                                                  <  1:3716348/91‑1 (MQ=255)
                                         gCAACATGTCGCATAACACTTCAACATAATTAGTACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGa                                                 <  1:3396984/91‑1 (MQ=255)
                                         gCAACATGTCGCATAACACTTCAACATAATTAGCACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGa                                                 <  1:3397144/91‑1 (MQ=255)
                                         gCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGa                                                 >  1:232402/1‑91 (MQ=255)
                                          cAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGaa                                                >  1:2160799/1‑91 (MQ=255)
                                          cAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGaa                                                >  1:3209330/1‑91 (MQ=255)
                                           aaCATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAAt                                               >  1:94419/1‑91 (MQ=255)
                                            aCATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATc                                              >  1:732015/1‑91 (MQ=255)
                                            aCATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATc                                              >  1:711171/1‑91 (MQ=255)
                                            aCATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATc                                              >  1:1898618/1‑91 (MQ=255)
                                              aTGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCtt                                            >  1:1776222/1‑91 (MQ=255)
                                              aTGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCtt                                            >  1:347502/1‑91 (MQ=255)
                                              aTGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCtt                                            >  1:2465150/1‑91 (MQ=255)
                                               tGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTc                                           >  1:3446748/1‑91 (MQ=255)
                                               tGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTc                                           >  1:3501701/1‑91 (MQ=255)
                                                gTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCa                                          >  1:324363/1‑91 (MQ=255)
                                                  cGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAAc                                        <  1:1694184/91‑1 (MQ=255)
                                                   gCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACt                                       >  1:4003877/1‑91 (MQ=255)
                                                      tAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTaaa                                    >  1:1707978/1‑91 (MQ=255)
                                                       aaCACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTaaaa                                   <  1:4048773/91‑1 (MQ=255)
                                                         cacTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAAtt                                 >  1:853930/1‑91 (MQ=255)
                                                           cTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTAc                               >  1:45465/1‑91 (MQ=255)
                                                            ttCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACg                              >  1:1736496/1‑91 (MQ=255)
                                                               aaCATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCt                           >  1:4034131/1‑91 (MQ=255)
                                                                aCATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCtt                          <  1:758979/91‑1 (MQ=255)
                                                                aCATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCtt                          <  1:758898/91‑1 (MQ=255)
                                                                   tAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTCAt                       >  1:3347304/1‑91 (MQ=255)
                                                                     aTTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTCATAg                     <  1:4027266/91‑1 (MQ=255)
                                                                     aTTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTCATAg                     <  1:2790018/91‑1 (MQ=255)
                                                                       tAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTCATAGCa                   <  1:1473803/91‑1 (MQ=255)
                                                                       aaGAACCCTGACAAAGATAGATATCTTTGCTAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTCATAGCa                   >  1:4070440/2‑91 (MQ=21)
                                                                        aGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTCATAGCaa                  >  1:2864163/1‑91 (MQ=255)
                                                                         gAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTCATAGCAAg                 >  1:2573269/1‑91 (MQ=255)
                                                                          aaCCCTGACAAAGATAGATATCTTTGCTAGGGTTCTTCTTTATGAATTATGAGATCGTATCTTCAACTAAAATTACGCCTTCATAGCAAgg                >  1:3405846/1‑91 (MQ=21)
                                                                          aaCCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTCATAGCAAgg                >  1:2158177/1‑91 (MQ=255)
                                                                           aCCCTGACAAAGGTAGATATCTTTGCTAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTCATAGCAAGGa               >  1:3184522/1‑91 (MQ=25)
                                                                           aCCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTCATAGCAAGGa               <  1:2576115/91‑1 (MQ=255)
                                                                           aCCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTCATAGCAAGGa               >  1:54458/1‑91 (MQ=255)
                                                                           aCCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTCATAGCAAGGa               >  1:489315/1‑91 (MQ=255)
                                                                            cccTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTCATAGCAAGGAc              >  1:3221577/1‑91 (MQ=255)
                                                                             ccTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTCATAGCAAGGACa             <  1:2230861/91‑1 (MQ=255)
                                                                             ccTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTCATAGCAAGGACa             >  1:1974007/1‑91 (MQ=255)
                                                                               tgacAAAGATAGATATCTTTGCTAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTCATAGCAAGGACAtt           >  1:2625345/4‑91 (MQ=37)
                                                                                aGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTCATAGCAAGGACAttt          <  1:3734147/91‑1 (MQ=255)
                                                                                  cAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTCATAGCAAGGACATTTCt        >  1:22108/1‑91 (MQ=255)
                                                                                    aaGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTCATAGCAAGGACATTTCTAt      <  1:730212/91‑1 (MQ=255)
                                                                                     aGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTCATAGCAAGGACATTTCTAtt     >  1:97571/1‑91 (MQ=255)
                                                                                      tatatCTATCTTTGTCAGGGTTCTTCTTTATGAATTATTAGATCGAATCTTCAACTAAAATTACGCCTTCATAGCAAGGACATTTCTATTc    <  1:1092137/90‑1 (MQ=255)
                                                                                      gATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTCATAGCAAGGACATTTCTATTc    >  1:442750/1‑91 (MQ=255)
                                                                                      gATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTCATAGCAAGGACATTTCTATTc    <  1:677396/91‑1 (MQ=255)
                                                                                        tatctatcTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTCATAGCAAGGACATTTCTATTCaa  >  1:1678076/1‑91 (MQ=255)
                                                                                        tatctatcTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTCATAGCAAGGACATTTCTATTCaa  <  1:664464/91‑1 (MQ=255)
                                                                                        tatctatcTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTCATAGCAAGGACATTTCTATTCaa  <  1:1056798/91‑1 (MQ=255)
                                                                                        tatctatcTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTCATAGCAAGGACATTTCTATTCaa  <  1:2982300/91‑1 (MQ=255)
                                                                                          |                                                                                        
TTAGTCAGTCTACAAATCATTTACTTTAAAATGATTAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTCATAGCAAGGACATTTCTATTCAA  >  USA300TCH1516_ALE/2559750‑2559928

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