Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 659,371 T→A 8.3% intergenic (+177/‑735) 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,3710TA 8.3% 805.4 / 44.4 249intergenic (+177/‑735)USA300TCH1516_00590/USA300HOU_RS03165hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (111/116);  new base A (9/12);  total (121/128)
Fisher's exact test for biased strand distribution p-value = 6.53e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.92e-01

Only 100 of 295 total aligned reads displayed.

TACGGTAACTTCTTTAATGTTCAGTTTGAAAAACAACTAGCAATCACTGTATTTTTTATCATGTTAATAGGGATGTTAATTTTTTATTTTTATCTAAATAAAAAATTAACATTAAAAATTTTTAATACCAACGTGGTTAATAAGAATAGAGTTGTATTAATACCGACTTTCAAACAAGGA  >  USA300TCH1516_ALE/659282‑659461
                                                                                         |                                                                                          
tACGGTAACTTCTTTAATGTTCAGTTTGAAAAACAACTAGCAATCACTGTATTTTTTATCATGTTAATAGGGATGTTAATtttttattttt                                                                                           >  1:1924819/1‑91 (MQ=255)
  cGGTAACTTCTTTAATGTTCAGTTTGAAAAACAACTAGCAATCACTGTATTTTTTATCATGTTAATAGGGATGTTAATTttttatttttat                                                                                         >  1:413901/1‑91 (MQ=255)
  cGGTAACTTCTTTAATGTTCAGTTTGAAAAACAACTAGCAATCACTGTATTTTTTATCATGTTAATAGGGATGTTAATTttttatttttat                                                                                         >  1:1386231/1‑91 (MQ=255)
   ggTAACTTCTTTAATGTTCAGTTTGAAAAACAACTAGCAATCACTGTATTTTTTATCATGTTAATAGGGATGTTAATTTTTTATTTTTATc                                                                                        >  2:778772/1‑91 (MQ=255)
   ggTAACTTCTTTAATGTTCAGTTTGAAAAACAACTAGCAATCACTGTATTTTTTATCATGTTAATAGGGATGTTAATTTTTTATTTTTATc                                                                                        <  1:3260987/91‑1 (MQ=255)
     tAACTTCTTTAATGTTCAGTTTGAAAAACAACTAGCAATCACTGTATTTTTTATCATGTTAATAGGGATGTTAATTTTTTATTTTTATCta                                                                                      >  1:2448565/1‑91 (MQ=255)
      aaCTTCTTTAATGTTCAGTTTGAAAAACAACTAGCAATCACTGTATTTTTTATCATGTTAATAGGGATGTTAATTTTTTATTTTTATCtaa                                                                                     >  2:2418913/1‑91 (MQ=255)
      aaCTTCTTTAATGTTCAGTTTGAAAAACAACTAGCAATCACTGTATTTTTTATCATGTTAATAGGGATGTTAATTTTTTATTTTTATCtaa                                                                                     >  2:1836291/1‑91 (MQ=255)
       aCTTCTTTAATGTTCAGTTTGAAAAACAACTAGCAATCACTGTATTTTTTATCATGTTAATAGGGATGTTAATTTTTTATTTTTATCtaaa                                                                                    <  2:1260936/91‑1 (MQ=255)
        cttcttTAATGTTCAGTTTGAAAAACAACTAGCAATCACTGTATTTTTTATCATGTTAATAGGGATGTTAATTTTTTATTTTTATCtaaat                                                                                   >  1:2760224/1‑91 (MQ=255)
         ttcttTAATGTTCAGTTTGAAAAACAACTAGCAATCACTGTATTTTTTATCATGTTAATAGGGATGTTAATTTTTTATTTTTATCtaaata                                                                                  <  1:3695096/91‑1 (MQ=255)
         ttcttTAATGTTCAGTTTGAAAAACAACTAGCAATCACTGTATTTTTTATCATGTTAATAGGGATGTTAATTTTTTATTTTTATCtaaata                                                                                  <  1:2117516/91‑1 (MQ=255)
           ttttAATGTTCAGTTTGAAAAACAACTAGCAATCACTGTATTTTTTATCATGTTAATAGGGATGTTAATTTTTTATTTTTATCtaaataaa                                                                                <  2:3315149/90‑1 (MQ=255)
             tgaaTGTTCAGTTTGAAAAGCAACTAGCAATCACTGTATTTTTTATCATGTTAATAGGGATGTTAATTTTTTATTTAGAtaaaaataaaaa                                                                              <  2:1526195/89‑1 (MQ=37)
              tAATGTTCAGTTTGAAAAACAACTAGCAATTACTGTATTTTTTATCATGTTAATAGGGATGTTAATTTTTTATTTAGATAAAAATaaaaaa                                                                             <  2:312178/91‑1 (MQ=37)
              tAATGTTCAGTTTGAAAAACAACTAGCAATCACTGTATTTTTTATCATGTTAATAGGGATGTTAATTTTTTATTTTTATCTAAATaaaaaa                                                                             >  1:1589875/1‑91 (MQ=255)
              tAATGTTCAGTTTGAAAAACAACTAGCAATCACTGTATTTTTTATCATGTTAATAGGGATGTTAATTTTTTATTTTTATCTAAATaaaaaa                                                                             >  2:1674126/1‑91 (MQ=255)
              tAATGTTCAGTTTGAAAAACAACTAGCAATCACTGTATTTTTTATCATGTTAATAGGGATGTTAATTTTTTATTTAGATAAAAATaaaaaa                                                                             <  1:332416/91‑1 (MQ=37)
              tAATGTTCAGTTTGAAAAACAACTAGCAATCACTGTATTTTTTATCATGTTAATAGGGATGTTAATTTTTTATTTAGATAAAAATaaaaaa                                                                             <  1:1852231/91‑1 (MQ=37)
              tAATGGTCAGTTTGAAAAACAACTAGCAATCACTGTATTTTTTATCAGGTTAATAGGGTTGTTAATTTTTTATTTTTATCTAAATaaaata                                                                             >  2:136928/1‑89 (MQ=255)
               aaTGTTCAGTTTGAAAAACAACTAGCAATCACTGTATTTTTTATCATGTTAATAGGGATGTTAATTTTTTATTTAGATAAAAATAAAAAAt                                                                            <  2:1913980/91‑1 (MQ=37)
                aTGTTCAGTTTGAAAAACAACTAGCAATCACTGTATTTTTTATCATGTTAATAGGGATGTTAATTTTTTATTTTTATCTAAATAAAAAAtt                                                                           >  2:2912706/1‑91 (MQ=255)
                aTGTTCAGTTTGAAAAACAACTAGCAATCACTGTATTTTTTATCATGTTAATAGGGATGTTAATTTTTTATTTTTATCTAAATAAAAAAtt                                                                           >  2:2913318/1‑91 (MQ=255)
                 tGTTCAGTTTGAAAAACAACTAGCTATCACTGTATTTTTTATCATGTTAATAGGGATGTTAATTTTTTATTTAGATAAAAATAAAAAATTa                                                                          <  2:1737084/91‑1 (MQ=25)
                  gTTCAGTTTGAAAAACAACTAGCAATCACTGTATTTTTTATCATGTTAATAGGGATGTTAATTTTTTATTTTTATCTAAATAAAAAATTaa                                                                         >  2:1540075/1‑91 (MQ=255)
                    tCAGTTTGAAAAACAACTAGCAATCACTGTATTTTTTATCATGTTAATAGGGATGTTAATTTTTTATTTTTATCTAAATAAAAAATTAACa                                                                       >  1:3219435/1‑91 (MQ=255)
                     tattttGAAAAACAACTAGCAATCACTGTATTTTTTATCATGTTAATAGGGATGTTAATTTTTTATTTTTATCTAAATAAAAAATTAACAt                                                                      <  2:3624077/88‑1 (MQ=255)
                     cAGTTTGAAAAACAACTAGCAATCACTGTATTTTTTATCATGTTAATAGGGATGTTAATTTTTTATTTTTATCTAAATAAAAAATTAACAt                                                                      >  2:268208/1‑91 (MQ=255)
                     cAGTTTGAAAAACAACTAGCAATCACTGTATTTTTTATCATGTTAATAGGGATGTTAATTTTTTATTTTTATCTAAATAAAAAATTAACAt                                                                      <  2:3494801/91‑1 (MQ=255)
                      aGTTTGAAAAACAACTAGCAGTTACTGTATTTTTTATCATGTTAATAGAGATGTTAATTTTTTATTTTTATCTAAATAAAAAATTAACAtt                                                                     <  2:1893174/91‑1 (MQ=255)
                      aGTTTGAAAAACAACTAGCAATCACTGTATTTTTTATCATGTTAATAGGGATGTTAATTTTTTATTTTTATCTAAATAAAAAATTAACAtt                                                                     <  1:485707/91‑1 (MQ=255)
                      aGTTTGAAAAACAACTAGCAATCACTGTATTTTTTATCATGTTAATAGGGATGTTAATTTTTTATTTTTATCTAAATAAAAAATTAACAtt                                                                     >  1:677025/1‑91 (MQ=255)
                      aGTTTGAAAAACAACTAGCAATCACTGTATTTTTTATCATGTTAATAGGGATGTTAATTTTTTATTTTTATCTAAATAAAAAATTAACAtt                                                                     <  2:2861875/91‑1 (MQ=255)
                       gTTTGAAAAACAACTAGCAATCACTGTATTTTTTATCATGTTAATAGGGATGTTAATTTTTTATTTTTATCTAAATAAAAAATTAACATTa                                                                    >  1:987036/1‑91 (MQ=255)
                        tttGAGAAACAACTAGCAATCACTGTATTTTTTATCATGTTAATAGGGATGTTAATTTTTTATTTTTATCTAAATAAAAAATTAACATTaa                                                                   >  1:2839518/1‑91 (MQ=255)
                         ttGAAAAACAACTAGCAATCACTGTATTTTTTATCATGTTAATAGGGATGTTAATTTTTTATTTTTATCTAAATAAAAAATTAACATTaaa                                                                  <  1:2440528/91‑1 (MQ=255)
                         ttGAAAAACAACTAGCAATCACTGTATTTTTTATCATGTTAATAGGGATGTTAATTTTTTATTTTTATCTAAATAAAAAATTAACATTaaa                                                                  <  1:812465/91‑1 (MQ=255)
                          tGAAAAACAACTAGCAATCACTGTATTTTTTATCATGTTAATAGGGATGTTAATTTTTTATTTTTATCTAAATAAAAAATTAACATTaaaa                                                                 <  1:1976026/91‑1 (MQ=255)
                           gAAAAACAACTAGCAATCACTGTATTTTTTATCATGTTAATAGGGATGTTAATTTTTTATTTTTATCTAAATAAAAAATTAACATTaaaaa                                                                >  2:3234163/1‑91 (MQ=255)
                              aaaCAACTAGCAATCACTGTATTTTTTATCATGTTAATAGGGATGTTAATTTTTTATTTTTATCTAAATAAAAAATTAACATTAAAAAttt                                                             >  1:1769719/1‑91 (MQ=255)
                                   acTAGCAATCACTGTATTTTTTATCATGTTAATNGGGATGTTAATTTTTTATTTTTATCTAAATAAAAAATTAACATTAAAAATTTTTAAt                                                        <  1:1177323/91‑1 (MQ=255)
                                   acTAGCAATCACTGTATTTTTTATCATGTTAATAGGGATGTTAATTTTTTATTTTTATCTAAATAAAAAATTAACATTAAAAATTTTTAAt                                                        >  1:911020/1‑91 (MQ=255)
                                   acTAGCAATCACTGTATTTTTTATCATGTTAATAGGGATGTTAATTTTTTATTTTTATCTAAATAAAAAATTAACATTAAAAATTTTTAAt                                                        <  2:571647/91‑1 (MQ=255)
                                    cTAGCAATCACTGTATTTTTTATCATGTTAATAGGGATGTTAATTTTTTATTTTTATCTAAATAAAAAATTAACATTAAAAATTTTTAATa                                                       >  1:2262379/1‑91 (MQ=255)
                                     tAGCAATCACTGTATTTTTTATCATGTTAATAGGGATGTTAATTTTTTATTTTTATCTAAATAAAAAATTAACATTAAAAATTTTTAATAc                                                      >  2:1453472/1‑91 (MQ=255)
                                     tAGCAATCACTGTATTTTTTATCATGTTAATAGGGATGTTAATTTTTTATTTTTATCTAAATAAAAAATTAACATTAAAAATTTTTAATAc                                                      <  2:2334902/91‑1 (MQ=255)
                                     tAGCAATCACTGTATTTTTTATCATGTTAATAGGGATGTTAATTTTTTATTTTTATCTAAATAAAAAATTAACATTAAAAATTTTTAATAc                                                      <  2:687964/91‑1 (MQ=255)
                                     tAGCAATCACTGTATTTTTTATCATGTTAATAGGGATGTTAATTTTTTATTTTTATCTAAATAAAAAATTAACATTAAAAATTTTTAATAc                                                      >  1:2072525/1‑91 (MQ=255)
                                      aGCAATCACTGTATTTTTTATCATGTTAATAGGGATGTTAATTTTTTATTTTTATCTAAATAAAAAATTAACATTAAAAATTTTTAATAcc                                                     >  1:3125823/1‑91 (MQ=255)
                                       gCAATCACTGTATTTTTTATCATGTTAATAGGGATGTTAATTTTTTATTTTTATCTAAATAAAAAATTAACATTAAAAATTTTTAATACCa                                                    >  2:3724753/1‑91 (MQ=255)
                                       gCAATCACTGTATTTTTTATCATGTTAATAGGGATGTTAATTTTTTATTTTTATATAAATAAAAAATTAACATTAAAAATTTTTAATACCa                                                    >  2:3053371/1‑91 (MQ=255)
                                       gCAATCACTGTATTTTTTATCATGTTAATAGGGATGTTAATTTGTTATTTTTATCTAAATAAAAAATTAACATTAAAAATTTTTAATACCa                                                    >  2:2988685/1‑91 (MQ=255)
                                        cAATCACTGTATTTTTTATCATGTTAATAGGGATGTTAATTTTTTATTTTTATCTAAATAAAAAATTAACATTAAAAATTTTTAATACCaa                                                   <  1:123245/91‑1 (MQ=255)
                                         aaTCACTGTATTTTTTATCATGTTAATAGGGATGTTAATTTTTTATTTTTATCTAAATAAAAAATTAACATTAAAAATTTTTAATACCAAc                                                  <  2:2566773/91‑1 (MQ=255)
                                          aTCACTGTATTTTTTATCATGTTAATAGGGATGTTAATTTTTTATTTTTATCTAAATAAAAAATTAACATTAAAAATTTTTAATACCAACg                                                 <  1:3838026/91‑1 (MQ=255)
                                          aTCACTGTATTTTTTATCATGTTAATAGGGATGTTAATTTTTTATTTTTATCTAAATAAAAAATTAACATTAAAAATTTTTAATACCAACg                                                 <  2:3168160/91‑1 (MQ=255)
                                           tCACTGTATTTTTTATCATGTTAATAGGGATGTTAATTTTTTATTTTTATCTAAATAAAAAATTAACATTAAAAATTTTTAATACCAACGt                                                <  1:3567912/91‑1 (MQ=255)
                                           tCACTGTATTTTTTATCATGTTAATAGGGATGTTAATTTTTTATTTTGATCTAAATAAAAAATTAACATTAAAAATTTTTAATACCAACGt                                                >  1:1155559/1‑91 (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)
                                                   ttttttATCATGTTAATAGGGATGTTAATTTTTTATTTTTATCTAAATAAAAAATTAACATTAAAAATTTTTAATACCAACGTGGTtaata                                        >  2:1190204/1‑91 (MQ=255)
                                                        tATCATGTTAATAGGGATGTTAATTTTTTATTTTTATCTCAATAAAAAATTAACATTAAAAATTTTTAATACCAACTTGGTTACTAAGAAt                                   >  2:2450533/1‑91 (MQ=255)
                                                          tCATGTTAATAGGGATGTTAATTTTTTATTTTTATCTAAATAAAAAATTAACATTAAAAATTTTTAATACCAACGTGGTTAATAAGAATag                                 >  2:2512393/1‑91 (MQ=255)
                                                          tCATGTTAATAGGGATGTTAATTTTTTATTTTTAGCTAGATAAAAAATGAACATTAAAAATTTTTAATGGCAGCGTGGTTAATAAGAATag                                 <  2:992157/91‑1 (MQ=255)
                                                           cATGTTAATAGGGATGTTAATTTTTTATTTTTATCTAAATAAAAAATTAACATTAAAAATTTTTAATACCAACGTGGTTAATAAGAATaga                                >  1:3464095/1‑91 (MQ=255)
                                                               ttAATAGGGATGTTAATTTTTTATTTTTATCTACATAAAAAATTAACATTAAAAATTTTTAATACCAACGTGGTTAATAAGAATAGAGTTg                            <  1:3173176/91‑1 (MQ=255)
                                                               ttAATAGGGATGTTAATTTTTTATTTTTATCTAAATAAAAAATTAACATTAAAAATTTTTAATACCAACGTGGTTAATAAGAATAGAGTTg                            <  2:1170200/91‑1 (MQ=255)
                                                               ttAATAGGGATGTAAATTTTTTATTTTTATCTAAATAAAAAATTAACATTAAAAATTTTTAATACCAACGTGGTTAATAAGAATAGAGTTg                            <  2:1422351/91‑1 (MQ=255)
                                                                 aaTAGGGATGTTAATTTTTTATTTTTATCTAAATAAAAAATTAACATTAAAAATTTTTAATACCAACGTGGTTAATAAGAATAGAGTTGTa                          >  1:2871438/1‑91 (MQ=255)
                                                                  aTAGGGATGTTAATTTTTTATTTTTATCTAAATAAAAAATTAACATTAAAAATTTTTAATACCAACGTGGTTAATAAGAATAGAGTTGTAt                         >  1:2904144/1‑91 (MQ=255)
                                                                   tAGGGATGTTAATTTTTTATTTTTATCTAAATAAAAAATTAACATTAAAAATTTTTAATACCAACGTGGTTAATAAGAATAGAGTTGTAtt                        <  1:1885543/91‑1 (MQ=255)
                                                                   tAGGGATGTTAATTTTTTATTTTTATCTAAATAAAAAATTAACATTAAAAATTTTTAATACCAACGTGGTTAATAAGAATAGAGTTGTAtt                        <  2:2746909/91‑1 (MQ=255)
                                                                    tttaaTGTTAATTTTTTATTTAGATAAAAATAAAAAATTAACATTAAAAATTTTTAATACCAACGTGGTTAATAAGAATAGAGTTGTATTa                       >  1:1683861/5‑91 (MQ=21)
                                                                    tttaaTGTTAATTTTTTATTTAGATAAAAATAAAAAATTAACATTAAAAATTTTCAATACCAACGAGGTTAATAAGAATAGAGTTGTATTa                       >  2:2591309/5‑91 (MQ=18)
                                                                     ttaaTGTTAATTTTTTATTTAGATAAAAATAAAAAATTAACATTAAAAATTTTTAATACCAACGTGGTTAATAAGAATAGAGTTGTATTaa                      >  2:2140401/4‑91 (MQ=21)
                                                                     gggATGTTAATTTTTTATTTTTATCTAAATAAAAAATTAACATTAAAAATTTTTAATACCAACGTGGTTAATAAGAATAGAGTTGTATTaa                      <  2:1589875/91‑1 (MQ=255)
                                                                      ggATGTTAATTTTTTATTTTTATCTAAATAAAAAATTAACATTAAAAATTTTTAATACCAACGTGGTTAATAAGAATAGAGTTGTATTAAt                     <  2:1992578/91‑1 (MQ=255)
                                                                      ggATGTTAATTTTTTATTTTTATCTAAATAAAAAATTAACATTAAAAATTTTTAATACCAACGTGGTTAATAAGAATAGAGTTGTATTAAt                     >  1:1269112/1‑91 (MQ=255)
                                                                      ggATGTTAATTTTTTATTTTTATCTAAATAAAAAATTAACATTAAAAATTTTTAATACCAACGTGGTTAATAAGAATAGAGTTGTATTAAt                     >  1:1888457/1‑91 (MQ=255)
                                                                       gATGTTAATTTTTTATTTTTATCTAAATAAAAAATTAACATTAAAAATTTTTAATACCAACGTGGTTAATAAGAATAGAGTTGTATTAATa                    >  1:1011792/1‑91 (MQ=255)
                                                                         tGTTAATTTTTTATTTTTATCTAAATAAAAAATTAACATTAAAAATTTTTAATACCAACGTGGTTAATAAGAATAGAGTTGTATTAATAcc                  <  1:3575481/91‑1 (MQ=255)
                                                                          gTTAATTTTTTATTTTTATCTAAATAAAAAATTAACATTAAAAATTTTTAATACCAACGTGGTTAATAAGAATAGAGTTGTATTAATACCg                 <  2:389576/91‑1 (MQ=255)
                                                                           ttAATTTTTTATTTTTATCTAAATAAAAAATTAACATTAAAAATTTTTAATACCAACGTGGTTAATAAGAATAGAGTTGTATTAATACCGa                >  1:2609675/1‑91 (MQ=255)
                                                                             aaTTTTTTATTTTTATCTAAATAAAAAATTAACATTAAAAATTTTTAATACCAACGTGGTTAATAAGAATAGAGTTGTATTAATACCGACt              <  1:1012548/91‑1 (MQ=255)
                                                                             aaTTTTTTATTTTTATCTAAATAAAAAATTAACATTAAAAATTTTTAATACCAACGTGGTTAATAAGAATAGAGTTGTATTAATACCGACt              <  2:1946023/91‑1 (MQ=255)
                                                                             aaTTTTTTATTTTTATCTAAATAAAAAATTAACATTAAAAATTTTTAATACCAACGTGGTTAATAAGAATAGAGTTGTATTAATACCGACt              >  2:1074606/1‑91 (MQ=255)
                                                                             aaTTTTTTATTTAGATAAAAATAAAAAATTAACATTAAAAATTTTTAATACCAACGTGGTTAATAAGAATAGAGTTGTATTAATACCGACt              >  1:513179/1‑91 (MQ=37)
                                                                              aTTTTTTATTTTTATCTAAATAAAAAATTAACATTAAAAATTTTTAATACCAACGTGGTTAATAAGAATAGAGTTGTATTAATACCGACtt             <  1:3388473/91‑1 (MQ=255)
                                                                              aTTTTTTATTTTTATCTAAATAAAAAATTAACATTAAAAATTTTTAATACCAACGTGGTTAATAAGAATAGAGTTGTATTAATACCGACtt             <  1:108714/91‑1 (MQ=255)
                                                                              aTTTTTTATTTTTATCTAAATAAAAAATTAACATTAAAAATTTTTAATACCAACGTGGTTAATAAGAATAGAGTTGTATTAATACCGACtt             <  1:839817/91‑1 (MQ=255)
                                                                               ttttttATTTTTATCTAAATAAAAAATTAACATTAAAAATTTTTAATACCAACGTGGTTAATAAGAATAGAGTTGTATTAATACCGACttt            >  1:1817414/1‑91 (MQ=255)
                                                                                   ttatttttatCTAAATAAAAAATTAACATTAAAAATTTTTAATACCAACGTGGTTAATAAGAATAGAGTTGTATTAATACCGACTTTCaaa        >  2:3647527/1‑91 (MQ=255)
                                                                                     atttttatCTAAATAAAAAATTAACATTAAAAATTTTTAATACCAACGTGGTTAATAAGAATAGAGTTGTATTAATACCGACTTTCAAACa      >  2:444863/1‑91 (MQ=255)
                                                                                     atttttatCTAAATAAAAAATTAACATTAAAAATTTTTAATACCAACGTGGTTAATAAGAATAGAGTTGTATTAATACCGACTTTCAAACa      <  2:3099741/91‑1 (MQ=255)
                                                                                        tttatttaGATAAAAAATTAACATTAAAAATTTTTAATACCAACGTGGTTAATAAGAATAGAGTTGTATTAATACCGACTTTCAAACAAgg   >  2:197772/1‑91 (MQ=255)
                                                                                        tttatttaGATAAAAAATTAACATTAAAAATTTTTAATACCAACGTGGTTAATAAGAATAGAGTTGTATTAATACCGACTTTCAAACAAgg   >  2:198006/1‑91 (MQ=255)
                                                                                        tttatCTAAATAAAAAATTAACATTAAAAATTTTTAATACCAACGTGGTTAATAAGAATAGAGTTGTATTAATACCGACTTTCAAACAAgg   >  2:3425057/1‑91 (MQ=255)
                                                                                        tttatCTAAATAAAAAATTAACATTAAAAATTTTTAATACCAACGTGGTTAATAAGAATAGAGTTGTATTAATACCGACTTTCAAACAAgg   >  1:3749695/1‑91 (MQ=255)
                                                                                         |                                                                                          
TACGGTAACTTCTTTAATGTTCAGTTTGAAAAACAACTAGCAATCACTGTATTTTTTATCATGTTAATAGGGATGTTAATTTTTTATTTTTATCTAAATAAAAAATTAACATTAAAAATTTTTAATACCAACGTGGTTAATAAGAATAGAGTTGTATTAATACCGACTTTCAAACAAGGA  >  USA300TCH1516_ALE/659282‑659461

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: