Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 659,376 T→A 7.7% intergenic (+182/‑730) USA300TCH1516_00590 → / → USA300HOU_RS03165 hypothetical protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE659,3760TA 7.7% 873.9 / 39.9 264intergenic (+182/‑730)USA300TCH1516_00590/USA300HOU_RS03165hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (116/127);  new base A (10/11);  total (126/138)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 290 total aligned reads displayed.

GTAACTTCTTTAATGTTCAGTTTGAAAAACAACTAGCAATCACTGTATTTTTTATCATGTTAATAGGGATGTTAATTTTTTATTTTTATCTAAATAAAAAATTAACATTAAAAATTTTTAATACCAACGTGGTTAATAAGAATAGAGTTGTATTAATACCGACTTTCAAACAAGGATTGTT  >  USA300TCH1516_ALE/659286‑659466
                                                                                          |                                                                                          
gTAACTTCTTTAATGTTCAGTTTGAAAAACAACTAGCAATCACTGTATTTTTTATCATGTTAATAGGGATGTTAATTTTTTATTTTTATCt                                                                                            >  2:3668263/1‑91 (MQ=255)
 tAACTTCTTTAATGTTCAGTTTGAAAAACAACTAGCAATCACTGTATTTTTTATCATGTTAATAGGGATGTTAATTTTTTATTTTTATCta                                                                                           >  1:2448565/1‑91 (MQ=255)
 tAACTTCTTTAATGTTCAGTTTGAAAAACAACTAGCAATCACTGTATTTTTTATCATGTTAATAGGGATGTTAATTTTTTATTTTTATCta                                                                                           <  1:794725/91‑1 (MQ=255)
  aaCTTCTTTAATGTTCAGTTTGAAAAACAACTAGCAATCACTGTATTTTTTATCATGTTAATAGGGATGTTAATTTTTTATTTTTATCtaa                                                                                          >  2:2418913/1‑91 (MQ=255)
  aaCTTCTTTAATGTTCAGTTTGAAAAACAACTAGCAATCACTGTATTTTTTATCATGTTAATAGGGATGTTAATTTTTTATTTTTATCtaa                                                                                          <  2:343900/91‑1 (MQ=255)
    cttcttTAATGTTCAGTTTGAAAAACAACTAGCAATCACTGTATTTTTTATCATGTTAATAGGGATGTTAATTTTTTATTTTTATCtaaat                                                                                        >  1:2760224/1‑91 (MQ=255)
     ttcttTAATGTTCAGTTTGAAAAACAACTAGCAATCACTGTATTTTTTATCATGTTAATAGGGATGTTAATTTTTTATTTTTATCtaaata                                                                                       <  1:3695096/91‑1 (MQ=255)
     ttcttTAATGTTCAGTTTGAAAAACAACTAGCAATCACTGTATTTTTTATCATGTTAATAGGGATGTTAATTTTTTATTTTTATCtaaata                                                                                       <  1:2117549/91‑1 (MQ=255)
         ttAATGTTCAGTTTGAAAAACAACTAGCAATCACTGTATTTTTTATCATGTTAATAGGGATGTTAATTTTTTATTTTTATATAAATaaaaa                                                                                   <  2:3528372/91‑1 (MQ=255)
          tAATGTTCAGTTTGAAAAACAACTAGCAATCTCTGTATTTTTTATCATGTTAATAGGGATGTTAATTTTTTATTTTTATCTAAATaaaaaa                                                                                  >  2:260073/1‑91 (MQ=255)
          tAATGTTCAGTTTGAAAAACAACTAGCAATCACTGTATTTTTTATCATGTTAATAGGGATGTTAATTTTTTATTTTTATCTAAATaaaaaa                                                                                  >  1:1589875/1‑91 (MQ=255)
          tAATGTTCAGTTTGAAAAACAACTAGCAATCACTGTATTTTTTATCATGTTAATAGGGATGTTAATTTTTTATTTTTATCTAAATaaaaaa                                                                                  >  2:1674126/1‑91 (MQ=255)
          tAATGTTCAGTTTGAAAAACAACTAGCAATCACTGTATTTTTTATCATGTTAATAGGGATGTTAATTTTTTATTTAGATAAAAATaaaaaa                                                                                  <  1:291226/91‑1 (MQ=37)
          tAATGTTCAGTTTGAAAAACAACTAGCAATCACTGTATTTTTTATCATGTTAATAGGGATGTTAATTTTTTATTTAGATAAAAATaaaaaa                                                                                  <  1:1852231/91‑1 (MQ=37)
          tAATGTTCAGTTTGAAAAACAACTAGCAATCACTGTATTTTTTATCATGTTAATAGGGATGTTAATTTTTTATTTAGATAAAAATaaaaaa                                                                                  <  1:332416/91‑1 (MQ=37)
          tAATGGTCAGTTTGAAAAACAACTAGCAATCACTGTATTTTTTATCAGGTTAATAGGGTTGTTAATTTTTTATTTTTATCTAAATaaaata                                                                                  >  2:136928/1‑89 (MQ=255)
           aaTGTTCAGTTTGAAAAACAACTAGCAATCACTGTATTTTTTATCATGTTAATAGGGATGTTAATTTTTTATTTAGATAAAAATAAAAAAt                                                                                 <  2:1913980/91‑1 (MQ=37)
            aTGTTCAGTTTGAAAAACAACTAGCAATCACTGTATTTTTTATCATGTTAATAGGGATGTTAATTTTTTATTTTTATCTAAATAAAAAAtt                                                                                >  2:2913123/1‑91 (MQ=255)
             tGTTCAGTTTGAAAAACAACTAGCTATCACTGTATTTTTTATCATGTTAATAGGGATGTTAATTTTTTATTTAGATAAAAATAAAAAATTa                                                                               <  2:1737084/91‑1 (MQ=25)
              gTTCAGTTTGAAAAACAACTAGCAATCACTGTATTTTTTATCATGTTAATAGGGATGTTAATTTTTTATTTTTATCTAAATAAAAAATTaa                                                                              >  2:1540075/1‑91 (MQ=255)
                tCATTTTGAAAAACAACTAGCAATCACTGTTTTTTTTATCATGTTAATAGGGATGTTAATTTATTATTTTTATCTAAATAAAAAATTAACa                                                                            <  2:1849148/91‑1 (MQ=255)
                tCAGTTTGAAAAACAACTAGCAATCACTGTATTTTTTATCATGTTAATAGGGATGTTAATTTTTTATTTTTATCTAAATAAAAAATTAACa                                                                            >  1:3219435/1‑91 (MQ=255)
                 aaGTTTGAAAAACAACTAGCAATCACTGTATTTTTTATCATGTTAATAGGGATGTTAATTTTTTATTTTTATCTAAATAAAAAATTAACAt                                                                           <  2:3624116/90‑1 (MQ=255)
                  aGTTTGAAAAACAACTAGCAGTTACTGTATTTTTTATCATGTTAATAGAGATGTTAATTTTTTATTTTTATCTAAATAAAAAATTAACAtt                                                                          <  2:1893174/91‑1 (MQ=255)
                  aGTTTGAAAAACAACTAGCAATCACTGTATTTTTTATCATGTTAATAGGGATGTTAATTTTTTATTTTTATCTAAATAAAAAATTAACAtt                                                                          <  2:2530646/91‑1 (MQ=255)
                  aGTTTGAAAAACAACTAGCAATCACTGTATTTTTTATCATGTTAATAGGGATGTTAATTTTTTATTTTTATCTAAATAAAAAATTAACAtt                                                                          <  1:2272957/91‑1 (MQ=255)
                     ttGAAAAACAACTAGCAATCACTGTATTTTTTATCATGTTAATAGGGATGTTAATTTTTTATTTTTATCTAAATAAAAAATTAACATTaaa                                                                       <  2:319415/91‑1 (MQ=255)
                     ttGAAAAACAACTAGCAATCACTGTATTTTTTATCATGTTAATAGGGATGTTAATTTTTTATTTTTATCTAAATAAAAAATTAACATTaaa                                                                       <  1:2440528/91‑1 (MQ=255)
                     ttGAAAAACAACTAGCAATCACTGTATTTTTTATCATGTTAATAGGGATGTTAATTTTTTATTTTTATCTAAATAAAAAATTAACATTaaa                                                                       <  1:1783242/91‑1 (MQ=255)
                       gAAAAACAACTAGCAATCACTGTATTTTTTATCATGTTAATAGGGATGTTAATTTTTTATTTTTATCTAAATAAAAAATTAACATTaaaaa                                                                     <  2:2693406/91‑1 (MQ=255)
                             caacTAGCAATCACTGTATTTTTTATCATGTTAATAGGGATGTTAATTTTTTATTTTTATCTAAATAAAAAATTAACATTAAAAATTTTTa                                                               >  2:2571153/1‑91 (MQ=255)
                               acTAGCAATCACTGTATTTTTTATCATGTTAATNGGGATGTTAATTTTTTATTTTTATCTAAATAAAAAATTAACATTAAAAATTTTTAAt                                                             <  1:1177323/91‑1 (MQ=255)
                               acTAGCAATCACTGTATTTTTTATCATGTTAATAGGGATGTTAATTTTTTATTTTTATCTAAATAAAAAATTAACATTAAAAATTTTTAAt                                                             >  1:823837/1‑91 (MQ=255)
                               acTAGCAATCACTGTATTTTTTATCATGTTAATAGGGATGTTAATTTTTTATTTTTATCTAAATAAAAAATTAACATTAAAAATTTTTAAt                                                             <  2:571647/91‑1 (MQ=255)
                                cTAGCAATCACTGTATTTTTTATCATGTTAATAGTGATGTTAATTTTTTATTTTTATCTAAATAAAAAATTAACATTAAAAATTTTTAATa                                                            <  2:1426097/91‑1 (MQ=255)
                                 tAGCAATCACTGTATTTTTTATCATGTTAATAGGGATGTTAATTTTTTATTTTTATCTAAATAAAAAATTAACATTAAAAATTTTTAATAc                                                           <  2:2334902/91‑1 (MQ=255)
                                 tAGCAATCACTGTATTTTTTATCATGTTAATAGGGATGTTAATTTTTTATTTTTATCTAAATAAAAAATTAACATTAAAAATTTTTAATAc                                                           <  2:687964/91‑1 (MQ=255)
                                 tAGCAATCACTGTATTTTTTATCATGTTAATAGGGATGTTAATTTTTTATTTTTATCTAAATAAAAAATTAACATTAAAAATTTTTAATAc                                                           >  1:2072743/1‑91 (MQ=255)
                                  aGCAATCACTGTATTTTTTATCATGTTAATAGGGATGTTAATTTTTTATTTTTATCTAAATAAAAAATTAACATTAAAAATTTTTAATAcc                                                          >  1:3125853/1‑91 (MQ=255)
                                   gCAATCACTGTATTTTTTATCATGTTAATAGGGATGTTAATTTTTTATTTTTATCTAAATAAAAAATTAACATTAAAAATTTTTAATACCa                                                         >  2:3724753/1‑91 (MQ=255)
                                   gCAATCACTGTATTTTTTATCATGTTAATAGGGATGTTAATTTTTTATTTTTATCTAAATAAAAAATTAACATTAAAAATTTTTAATACCa                                                         >  2:1464785/1‑91 (MQ=255)
                                    cAATCACTGTATTTTTTATCATGTTAATAGGGATGTTAATTTTTTATTTTTATCTAAATAAAAAATTAACATTAAAAATTTTTAATACCaa                                                        <  1:123245/91‑1 (MQ=255)
                                      aTCACTGTATTTTTTATCATGTTAATAGGGATGTTAATTTTTTATTTTTATCTAAATAAAAAATTAACATTAAAAATTTTTAATACCAACg                                                      <  1:3838026/91‑1 (MQ=255)
                                       tCACTGTATTTTTTATCATGTTAATAGGGATGTTAATTTTTTATTTTTATCTAAATAAAAAATTAACATTAAAAATTTTTAATACCAACGt                                                     <  1:3567912/91‑1 (MQ=255)
                                       tCACTGTATTTTTTATCATGTTAATAGGGATGTTAATTTTTTATTTTGATCTAAATAAAAAATTAACATTAAAAATTTTTAATACCAACGt                                                     >  1:1155559/1‑91 (MQ=255)
                                        cACTGTATTTTTTATCATGTTAATAGGGATGTTAATTTTTTATTTTTATCTAAATAAAAAATTAACATTAAAAATTTTTAATACCAACGTg                                                    <  2:1122597/91‑1 (MQ=255)
                                        cACTGTATTTTTTATCATGTTAATAGGGATGTTAATTTTTTATTTTTATCTAAATAAAAAATTAACATTAAAAATTTTTAATACCAACGTg                                                    >  2:526597/1‑91 (MQ=255)
                                         aCTGTATTTTTTATCATGTTAATAGGGATGTTAATTTTTTATTTTTATCTAAATAAAAAATTAACATTAAAAATTTTTAATACCAACGTgg                                                   >  2:17763/1‑91 (MQ=255)
                                          cTGTATTTTTTATCATGTTAATAGGGATGTTAATTTTTTATTTTTATCTAAATAAAAAATTAACATTAAAAATTTTTAATACCAACGTGGt                                                  <  1:2329725/91‑1 (MQ=255)
                                            gTATTTTTTATCATGTTAATAGGGATGTTAATTTTTTATTTTTATCTAAATAAAAAATTAACATTAAAAATTTTTAATACCAACGTGGTta                                                >  2:488657/1‑91 (MQ=255)
                                                 ttttATCATGTTAATAGGGATGTTAATTTTTTATTTTTATCTAAATAAAAAATTAACATTAAAAATTTTTAATACCAACGTGGTTAATAAg                                           <  1:192789/91‑1 (MQ=255)
                                                 ttttATCATGTTAATAGGGATGTTAATTTTTTATTTTTATCTAAATAAAAAATTAACATTAAAAATTTTTAATACCAACGTGGTTAATAAg                                           <  1:488657/91‑1 (MQ=255)
                                                  tttATCATGTTAATAGGGATGTTAATTTTTTATTTTTATCTAAATAAAAAATTAACATTAAAAATTTTTAATACCAACGTGGTTAATAAGa                                          >  2:3040971/1‑91 (MQ=255)
                                                    tATCATGTTAATAGGGATGTTAATTTTTTATTTTTATCTCAATAAAAAATTAACATTAAAAATTTTTAATACCAACTTGGTTACTAAGAAt                                        >  2:2450533/1‑91 (MQ=255)
                                                      tCATGTTAATAGGGATGTTAATTTTTTATTTTTAGCTAGATAAAAAATGAACATTAAAAATTTTTAATGGCAGCGTGGTTAATAAGAATag                                      <  2:992157/91‑1 (MQ=255)
                                                       cATGTTAATAGGGATGTTAATTTTTTATTTTTATCTAAATAAAAAATTAACATTAAAAATTTTTAATACCAACGTGGTTAATAAGAATaga                                     >  1:3464095/1‑91 (MQ=255)
                                                        aTGTTAATAGGGATGTTAATTTTTTATTTTTATCTAAATAAAAAATTAACATTAAAAATTTTTAATACCAACGTGGTTAATAAGAATagag                                    >  2:2755391/1‑91 (MQ=255)
                                                          gTTAATAGGGATGTTAATTTTTTATTTTTATCTAAATAAAAAATTAACATTAAAAATTTTTAATACCAACGTGGTTAATAAGAATAGAGtt                                  <  1:3234163/91‑1 (MQ=255)
                                                           ttAATAGGGATGTTAATTTTTTATTTTTATCTACATAAAAAATTAACATTAAAAATTTTTAATACCAACGTGGTTAATAAGAATAGAGTTg                                 <  1:3173176/91‑1 (MQ=255)
                                                           ttAATAGGGATGTTAATTTTTTATTTTTATCTAAATAAAAAATTAACATTAAAAATTTTTAATACCAACGTGGTTAATAAGAATAGAGTTg                                 <  1:2880738/91‑1 (MQ=255)
                                                           ttAATAGGGATGTTAATTTTTTATTTTTATCTAAATAAAAAATTAACATTAAAAATTTTTAATACCAACGTGGTTAATAAGAATAGAGTTg                                 <  2:1172315/91‑1 (MQ=255)
                                                             aaTAGGGATGTTAATTTTTTATTTTTATCTAAATAAAAAATTAACATTAAAAATTTTTAATACCAACGTGGTTAATAAGAATAGAGTTGTa                               >  1:914063/1‑91 (MQ=255)
                                                             aaTAGGGATGTTAATTTTTTATTTTTATCTAAATAAAAAATTAACATTAAAAATTTTTAATACCAACGTGGTTAATAAGAATAGAGTTGTa                               >  1:1992578/1‑91 (MQ=255)
                                                               tAGGGATGTTAATTTTTTATTTTTATCTAAATAAAAAATTAACATTAAAAATTTTTAATACCAACGTGGTTAATAAGAATAGAGTTGTAtt                             <  1:1885543/91‑1 (MQ=255)
                                                                tttaaTGTTAATTTTTTATTTAGATAAAAATAAAAAATTAACATTAAAAATTTTTAATACCAACGTGGTTAATAAGAATAGAGTTGTATTa                            >  1:1683861/5‑91 (MQ=21)
                                                                 ttaaTGTTAATTTTTTATTTAGATAAAAATAAAAAATTAACATTAAAAATTTTTAATACCAACGTGGTTAATAAGAATAGAGTTGTATTaa                           >  2:2140401/4‑91 (MQ=21)
                                                                 gggATGTTAATTTTTTATTTTTATCTAAATAAAAAATTAACATTAAAAATTTTTAATACCAACGTGGTTAATAAGAATAGAGTTGTATTaa                           <  2:1589875/91‑1 (MQ=255)
                                                                 gggATGTTAATTTTTTATTTTTATCTAAATAAAAAATTAACATTAAAAATTTTTAATACCAACGTGGTTAATAAGAATAGAGTTGTATTaa                           >  2:324440/1‑91 (MQ=255)
                                                                  ggATGTTAATTTTTTATTTTTATCTAAATAAAAAATTAACATTAAAAATTTTTAATACCAACGTGGTTAATAAGAATAGAGTTGTATTAAt                          >  1:1888457/1‑91 (MQ=255)
                                                                  ggATGTTAATTTTTTATTTTTATCTAAATAAAAAATTAACATTAAAAATTTTTAATACCAACGTGGTTAATAAGAATAGAGTTGTATTAAt                          >  1:1269112/1‑91 (MQ=255)
                                                                  ggATGTTAATTTTTTATTTTTATCTAAATAAAAAATTAACATTAAAAATTTTTAATACCAACGTGGTTAATAAGAATAGAGTTGTATTAAt                          <  2:1992578/91‑1 (MQ=255)
                                                                   gATGTTAATTTTTTATTTTTATCTAAATAAAAAATTAACATTAAAAATTTTTAATACCAACGTGGTTACTAAGAATAGAGTTGTATTAATa                         >  2:3132301/1‑91 (MQ=255)
                                                                   gATGTTAATTTTTTATTTTTATCTAAATAAAAAATTAACATTAAAAATTTTTAATACCAACGTGGTTAATAAGAATAGAGTTGTATTAATa                         >  1:1011792/1‑91 (MQ=255)
                                                                    aTGTTAATTTTTTATTTTTATCTAAATAAAAAATTAACATTAAAAATTTTTAATACCAACGTGGTTAATAAGAATAGAGTTGTATTAATAc                        >  2:1269145/1‑91 (MQ=255)
                                                                    aTGTTAATTTTTTATTTTTATCTAAATAAAAAATTAACATTAAAAATTTTTAATACCAACGTGGTTAATAAGAATAGAGTTGTATTAATAc                        <  1:433403/91‑1 (MQ=255)
                                                                     tGTTAATTTTTTATTTTTATCTAAATAAAAAATTAACATTAAAAATTTTTAATACCAACGTGGTTAATAAGAATAGAGTTGTATTAATAcc                       <  1:3575481/91‑1 (MQ=255)
                                                                      gTTAATTTTTTATTTTTATCTAAATAAAAAATTAACATTAAAAATTTTTAATACCAACGTGGTTAATAAGAATAGAGTTGTATTAATACCg                      >  2:3352905/1‑91 (MQ=255)
                                                                      gTTAATTTTTTATTTTTATCTAAATAAAAAATTAACATTAAAAATTTTTAATACCAACGTGGTTAATAAGAATAGAGTTGTATTAATACCg                      <  2:389576/91‑1 (MQ=255)
                                                                      gTTAATTTTTTATTTTTATCTAAATAAAAAATTAACATTAAAAATTTTTAATACCAACGTGGTTAATAAGAATAGAGTTGTATTAATACCg                      >  2:3058715/1‑91 (MQ=255)
                                                                       ttAATTTTTTATTTTTATCTAAATAAAAAATTAACATTAAAAATTTTTAATACCAACGTGGTTAATAAGAATAGAGTTGTATTAATACCGa                     >  1:2609675/1‑91 (MQ=255)
                                                                        tAATTTTTTATTTTTATCTAAATAAAAAATTAACATTAAAAATTTTTAATACCAACGTGGTTAATAAGAATAGAGTTGTATTAATACCGAc                    >  2:3100658/1‑91 (MQ=255)
                                                                         aaTTTTTTATTTTTATCTAAATAAAAAATTAACATTAAAAATTTTTAATACCAACGTGGTTAATAAGAATAGAGTTGTATTAATACCGACt                   <  1:1012548/91‑1 (MQ=255)
                                                                         aaTTTTTTATTTTTATCTAAATAAAAAATTAACATTAAAAATTTTTAATACCAACGTGGTTAATAAGAATAGAGTTGTATTAATACCGACt                   <  2:1946023/91‑1 (MQ=255)
                                                                          aTTTTTTATTTTTATCTAAATAAAAAATTAACATTAAAAATTTTTAATACCAACGTGGTTAATAAGAATAGAGTTGTATTAATACCGACtt                  <  1:108714/91‑1 (MQ=255)
                                                                          aTTTTTTATTTTTATCTAAATAAAAAATTAACATTAAAAATTTTTAATACCAACGTGGTTAATAAGAATAGAGTTGTATTAATACCGACtt                  <  1:3388473/91‑1 (MQ=255)
                                                                          aTTTTTTA‑‑TTTATCTAAATAAAAAATTAACATTAAAAATTTTTAATACCAACGTGGTTAATAAGAATAGAGTTGTATTAATACCGACTTTc                <  2:765164/91‑1 (MQ=255)
                                                                           ttttttATTTTTATCTAAATAAAAAATTAACATTAAAAATTTTTAATACCAACGTGGTTAATAAGAATAGAGTTGTATTAATACCGACttt                 >  1:1817414/1‑91 (MQ=255)
                                                                           ttttttATTTTTATCTAAATAAAAAATTAACATTAAAAATTTTTAATACCAACGTGGTTAATAAGAATAGAGTTGTATTAATACCGACttt                 >  1:2844443/1‑91 (MQ=255)
                                                                           ttttttATTTTTATCTAAATAAAAAATTAACATTAAAAATTTTTAATACCAACGTGGTTAATAAGAATAGAGTTGTATTAATACCGACttt                 >  2:124187/1‑91 (MQ=255)
                                                                            tttttatttttaTCTAAATAAAAAATTAACATTAAAAATTTTTAATACCAACGTGGTTAATAAGAATAGAGTTGTATTAATACCGACTTTc                >  1:584100/1‑91 (MQ=255)
                                                                                tatttttatCTAAATAAAAAATTAACATTAAAAATTTTTAATACCAACGTGGTTAATAAGAATAGAGTTGTATTAATACCGACTTTCAAAc            <  1:876864/91‑1 (MQ=255)
                                                                                tatttttatCTAAATAAAAAATTAACATTAAAAATTTTTAATACCAACGTGGTTAATAAGAATAGAGTTGTATTAATACCGACTTTCAAAc            <  2:2884751/91‑1 (MQ=255)
                                                                                 atttttatCTAAATAATAAATTAACATTAAAAATTTTTAATACCAACGTGGTTAATAAGAATAGAGTTGTATTAATACCGACTTTCAAACa           <  2:2985018/91‑1 (MQ=255)
                                                                                 atttttatCTAAATAAAAAATTAACATTAAAAATTTTTAATACCAACGTGGTTAATAAGAATAGAGTTGTATTAATACCGACTTTCAAACa           >  2:444863/1‑91 (MQ=255)
                                                                                  tttttatCTAAATAAAAAATTAACATTAAAAATTTTTAATACCAACGTGGTTAATAAGAATAGAGTTGTATTAATACCGACTTTCAAACaa          >  2:103980/1‑91 (MQ=255)
                                                                                   ttttatCTAAATAAAAAATTAACATTAAAAATTTTTAATACCAACGTGGTTAATAAGAATAGAGTTGTATTAATACCGACTTTCAAACAAg         <  1:1368765/91‑1 (MQ=255)
                                                                                    tttatttaGATAAAAAATTAACATTAAAAATTTTTAATACCAACGTGGTTAATAAGAATAGAGTTGTATTAATACCGACTTTCAAACAAgg        >  2:197772/1‑91 (MQ=255)
                                                                                    tttatttaGATAAAAAATTAACATTAAAAATTTTTAATACCAACGTGGTTAATAAGAATAGAGTTGTATTAATACCGACTTTCAAACAAgg        >  2:198006/1‑91 (MQ=255)
                                                                                    tttatCTAAATAAAAAATTAACATTAAAAATTTTTAATACCAACGTGGTTAATAAGAATAGAGTTGTATTAATACCGACTTTCAAACAAgg        >  1:3749695/1‑91 (MQ=255)
                                                                                      tatCTAAATAAAAAATTAACATTAAAAATTTTTAATACCAACGTGGTTAATAAGAATAGAGTTGTATTAATACCGACTTTCAAACAAGGAt      >  2:249831/1‑91 (MQ=255)
                                                                                          |                                                                                          
GTAACTTCTTTAATGTTCAGTTTGAAAAACAACTAGCAATCACTGTATTTTTTATCATGTTAATAGGGATGTTAATTTTTTATTTTTATCTAAATAAAAAATTAACATTAAAAATTTTTAATACCAACGTGGTTAATAAGAATAGAGTTGTATTAATACCGACTTTCAAACAAGGATTGTT  >  USA300TCH1516_ALE/659286‑659466

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