Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 659,375 C→A 7.3% intergenic (+181/‑731) 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,3750CA 7.3% 946.7 / 43.9 281intergenic (+181/‑731)USA300TCH1516_00590/USA300HOU_RS03165hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (122/137);  new base A (11/11);  total (133/148)
Fisher's exact test for biased strand distribution p-value = 8.27e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 295 total aligned reads displayed.

GGTAACTTCTTTAATGTTCAGTTTGAAAAACAACTAGCAATCACTGTATTTTTTATCATGTTAATAGGGATGTTAATTTTTTATTTTTATCTAAATAAAAAATTAACATTAAAAATTTTTAATACCAACGTGGTTAATAAGAATAGAGTTGTATTAATACCGACTTTCAAACAAGGATTGT  >  USA300TCH1516_ALE/659285‑659465
                                                                                          |                                                                                          
ggTAACTTCTTTAATGTTCAGTTTGAAAAACAACTAGCAATCACTGTATTTTTTATCATGTTAATAGGGATGTTAATTTTTTATTTTTATc                                                                                            >  2:778772/1‑91 (MQ=255)
ggTAACTTCTTTAATGTTCAGTTTGAAAAACAACTAGCAATCACTGTATTTTTTATCATGTTAATAGGGATGTTAATTTTTTATTTTTATc                                                                                            >  1:571647/1‑91 (MQ=255)
ggTAACTTCTTTAATGTTCAGTTTGAAAAACAACTAGCAATCACTGTATTTTTTATCATGTTAATAGGGATGTTAATTTTTTATTTTTATc                                                                                            <  1:3260987/91‑1 (MQ=255)
 gTAACTTCTTTAATGTTCAGTTTGAAAAACAACTAGCAATCACTGTATTTTTTATCATGTTAATAGGGATGTTAATTTTTTATTTTTATCt                                                                                           >  2:3668263/1‑91 (MQ=255)
   aaCTTCTTTAATGTTCAGTTTGAAAAACAACTAGCAATCACTGTATTTTTTATCATGTTAATAGGGATGTTAATTTTTTATTTTTATCtaa                                                                                         >  2:1836291/1‑91 (MQ=255)
   aaCTTCTTTAATGTTCAGTTTGAAAAACAACTAGCAATCACTGTATTTTTTATCATGTTAATAGGGATGTTAATTTTTTATTTTTATCtaa                                                                                         >  2:645519/1‑91 (MQ=255)
   aaCTTCTTTAATGTTCAGTTTGAAAAACAACTAGCAATCACTGTATTTTTTATCATGTTAATAGGGATGTTAATTTTTTATTTTTATCtaa                                                                                         <  2:343900/91‑1 (MQ=255)
      ttcttTAATGTTCAGTTTGAAAAACAACTAGCAATCACTGTATTTTTTATCATGTTAATAGGGATGTTAATTTTTTATTTTTATCtaaata                                                                                      <  1:3695096/91‑1 (MQ=255)
          ttAATGTTCAGTTTGAAAAACAACTAGCAATCACTGTATTTTTTATCATGTTAATAGGGATGTTAATTTTTTATTTTTATCTAAATaaaaa                                                                                  <  2:1973765/91‑1 (MQ=255)
          ttAATGTTCAGTTTGAAAAACAACTAGCAATCACTGTATTTTTTATCATGTTAATAGGGATGTTAATTTTTTATTTTTATATAAATaaaaa                                                                                  <  2:3528372/91‑1 (MQ=255)
          tgaaTGTTCAGTTTGAAAAGCAACTAGCAATCACTGTATTTTTTATCATGTTAATAGGGATGTTAATTTTTTATTTAGAtaaaaataaaaa                                                                                  <  2:1526195/89‑1 (MQ=37)
           tAATGTTCAGTTTGAAAAACAACTAGCAATCTCTGTATTTTTTATCATGTTAATAGGGATGTTAATTTTTTATTTTTATCTAAATaaaaaa                                                                                 >  2:260073/1‑91 (MQ=255)
           tAATGTTCAGTTTGAAAAACAACTAGCAATCACTGTATTTTTTATCATGTTAATAGGGATGTTAATTTTTTATTTTTATCTAAATaaaaaa                                                                                 >  1:1589875/1‑91 (MQ=255)
           tAATGTTCAGTTTGAAAAACAACTAGCAATCACTGTATTTTTTATCATGTTAATAGGGATGTTAATTTTTTATTTAGATAAAAATaaaaaa                                                                                 <  1:1852231/91‑1 (MQ=37)
           tAATGTTCAGTTTGAAAAACAACTAGCAATCACTGTATTTTTTATCATGTTAATAGGGATGTTAATTTTTTATTTAGATAAAAATaaaaaa                                                                                 <  1:332416/91‑1 (MQ=37)
           tAATGTTCAGTTTGAAAAACAACTAGCAATCACTGTATTTTTTATCATGTTAATAGGGATGTTAATTTTTTATTTAGATAAAAATaaaaaa                                                                                 <  1:291226/91‑1 (MQ=37)
           tAATGGTCAGTTTGAAAAACAACTAGCAATCACTGTATTTTTTATCAGGTTAATAGGGTTGTTAATTTTTTATTTTTATCTAAATaaaata                                                                                 >  2:136928/1‑89 (MQ=255)
             aTGTTCAGTTTGAAAAACAACTAGCAATCACTGTATTTTTTATCATGTTAATAGGGATGTTAATTTTTTATTTTTATCTAAATAAAAAAtt                                                                               >  2:2913123/1‑91 (MQ=255)
                 tCAGTTTGAAAAACAACTAGCAATCACTGTATTTTTTATCATGTTAATAGGGATGTTAATTTTTTATTTTTATCTAAATAAAAAATTAACa                                                                           >  1:3219435/1‑91 (MQ=255)
                  aaGTTTGAAAAACAACTAGCAATCACTGTATTTTTTATCATGTTAATAGGGATGTTAATTTTTTATTTTTATCTAAATAAAAAATTAACAt                                                                          <  2:3624116/90‑1 (MQ=255)
                   aGTTTGAAAAACAACTAGCAATCACTGTATTTTTTATCATGTTAATAGGGATGTTAATTTTTTATTTTTATCTAAATAAAAAATTAACAtt                                                                         <  2:2530646/91‑1 (MQ=255)
                   aGTTTGAAAAACAACTAGCAATCACTGTATTTTTTATCATGTTAATAGGGATGTTAATTTTTTATTTTTATCTAAATAAAAAATTAACAtt                                                                         >  1:677777/1‑91 (MQ=255)
                   aGTTTGAAAAACAACTAGCAATCACTGTATTTTTTATCATGTTAATAGGGATGTTAATTTTTTATTTTTATCTAAATAAAAAATTAACAtt                                                                         <  2:1904950/91‑1 (MQ=255)
                   aGTTTGAAAAACAACTAGCAATCACTGTATTTTTTATCATGTTAATAGGGATGTTAATTTTTTATTTTTATCTAAATAAAAAATTAACAtt                                                                         >  2:1923058/1‑91 (MQ=255)
                    gTTTGAAAAACAACTAGCAATCACTGTATTTTTTATCATGTTAATAGGGATGTTAATTTTTTATTTTTATCTAAATAAAAAATTAACATTa                                                                        >  1:987036/1‑91 (MQ=255)
                      ttGAAAAACAACTAGCAATCACTGTATTTTTTATCATGTTAATAGGGATGTTAATTTTTTATTTTTATCTAAATAAAAAATTAACATTaaa                                                                      <  2:319415/91‑1 (MQ=255)
                      ttGAAAAACAACTAGCAATCACTGTATTTTTTATCATGTTAATAGGGATGTTAATTTTTTATTTTTATCTAAATAAAAAATTAACATTaaa                                                                      <  1:812570/91‑1 (MQ=255)
                      ttGAAAAACAACTAGCAATCACTGTATTTTTTATCATGTTAATAGGGATGTTAATTTTTTATTTTTATCTAAATAAAAAATTAACATTaaa                                                                      <  1:2441082/91‑1 (MQ=255)
                      ttGAAAAACAACTAGCAATCACTGTATTTTTTATCATGTTAATAGGGATGTTAATTTTTTATTTTTATCTAAATAAAAAATTAACATTaaa                                                                      <  1:227390/91‑1 (MQ=255)
                        gAAAAACAACTAGCAATCACTGTATTTTTTATCATGTTAATAGGGATGTTAATTTTTTATTTTTATCTAAATAAAAAATTAACATTaaaaa                                                                    <  2:2693406/91‑1 (MQ=255)
                         aaaaaCAACTAGCAATCACTGTATTTTTTATCATGTTAATAGGGATGTTAATTTTTTATTTTTATCTAAATAAAAAATTAACATTAAAAAt                                                                   <  2:1219346/91‑1 (MQ=255)
                           aaaCAACTAGCAATCACTGTATTTTTTATCATGTTAATAGGGATGTTAATTTTTTATTTTTATCTAAATAAAAAATTAACATTAAAAAttt                                                                 >  1:1769719/1‑91 (MQ=255)
                             acaacTAGCAATCACTGTATTTTTTATCATGTTAATAGGGATGTTAATTTTTTATTTTTATCTAAATAAAAAATTAACATTAAAAAttttt                                                               <  2:1878224/91‑1 (MQ=255)
                              caacTAGCAATCACTGTATTTTTTATCATGTTAATAGGGATGTTAATTTTTTATTTTTATCTAAATAAAAAATTAACATTAAAAATTTTTa                                                              >  2:2571153/1‑91 (MQ=255)
                                acTAGCAATCACTGTATTTTTTATCATGTTAATNGGGATGTTAATTTTTTATTTTTATCTAAATAAAAAATTAACATTAAAAATTTTTAAt                                                            <  1:1177323/91‑1 (MQ=255)
                                acTAGCAATCACTGTATTTTTTATCATGTTAATAGGGATGTTAATTTTTTATTTTTATCTAAATAAAAAATTAACATTAAAAATTTTTAAt                                                            >  1:911020/1‑91 (MQ=255)
                                  tAGCAATCACTGTATTTTTTATCATGTTAATAGGGATGTTAATTTTTTATTTTTATCTAAATAAAAAATTAACATTAAAAATTTTTAATAc                                                          <  2:376125/91‑1 (MQ=255)
                                  tAGCAATCACTGTATTTTTTATCATGTTAATAGGGATGTTAATTTTTTATTTTTATCTAAATAAAAAATTAACATTAAAAATTTTTAATAc                                                          >  2:1453472/1‑91 (MQ=255)
                                   aGCAATCACTGTATTTTTTATCATGTTAATAGGGATGTTAATTTTTTATTTTTATCTAAATAAAAAATTAACATTAAAAATTTTTAATAcc                                                         <  2:2448565/91‑1 (MQ=255)
                                   aGCAATCACTGTATTTTTTATCATGTTAATAGGGATGTTAATTTTTTATTTTTATCTAAATAAAAAATTAACATTAAAAATTTTTAATAcc                                                         >  1:3125853/1‑91 (MQ=255)
                                    gCAATCACTGTATTTTTTATCATGTTAATAGGGATGTTAATTTTTTATTTTTATCTAAATAAAAAATTAACATTAAAAATTTTTAATACCa                                                        >  2:2135907/1‑91 (MQ=255)
                                     cAATCACTGTATTTTTTATCATGTTAATAGGGATGTTAATTTTTTATTTTTATCTAAATAAAAAATTAACATTAAAAATTTTTAATACCaa                                                       <  1:123245/91‑1 (MQ=255)
                                     cAATCACTGTATTTTTTATCATGTTAATAGGGATGTTAATTTTTTATTTTTATCTAAATAAAAAATTAACATTAAAAATTTTTAATACCaa                                                       >  2:1757840/1‑91 (MQ=255)
                                      aaTCACTGTATTTTTTATCATGTTAATAGGGATGTTAATTTTTTATTTTTATCTAAATAAAAAATTAACATTAAAAATTTTTAATACCAAc                                                      >  1:868368/1‑91 (MQ=255)
                                       aTCACTGTATTTTTTATCATGTTAATAGGGATGTTAATTTTTTATTTTTATCTAAATAAAAAATTAACATTAAAAATTTTTAATACCAACg                                                     <  1:3838026/91‑1 (MQ=255)
                                        tCACTGTATTTTTTATCATGTTAATAGGGATGTTAATTTTTTATTTTTATCTAAATAAAAAATTAACATTAAAAATTTTTAATACCAACGt                                                    <  1:3567912/91‑1 (MQ=255)
                                        tCACTGTATTTTTTATCATGTTAATAGGGATGTTAATTTTTTATTTTTATCTAAATAAAAAATTAACATTAAAAATTTTTAATACCAACGt                                                    <  2:170955/91‑1 (MQ=255)
                                        tCACTGTATTTTTTATCATGTTAATAGGGATGTTAATTTTTTATTTTGATCTAAATAAAAAATTAACATTAAAAATTTTTAATACCAACGt                                                    >  1:1155559/1‑91 (MQ=255)
                                          aCTGTATTTTTTATCATGTTAATAGGGATGTTAATTTTTTATTTTTATCTAAATAAAAAATTAACATTAAAAATTTTTAATACCAACGTgg                                                  <  1:245298/91‑1 (MQ=255)
                                                 tttttATCATGTTAATAGGGATGTTAATTTTTTATTTTTATCTAAATAAAAAATTAACATTAAAAATTTTTAATACCAACGTGGTtaataa                                           >  2:1626768/1‑91 (MQ=255)
                                                  ttttATCATGTTAATAGGGATGTTAATTTTTTATTTTTATCTAAATAAAAAATTAACATTAAAAATTTTTAATACCAACGTGGTTAATAAg                                          <  1:488657/91‑1 (MQ=255)
                                                  ttttATCATGTTAATAGGGATGTTAATTTTTTATTTTTATCTAAATAAAAAATTAACATTAAAAATTTTTAATACCAACGTGGTTAATAAg                                          >  2:734584/1‑91 (MQ=255)
                                                  ttttATCATGTTAATAGGGATGTTAATTTTTTATTTTTATCTAAATAAAAAATTAACATTAAAAATTTTTAATACCAACGTGGTTAATAAg                                          <  1:2114591/91‑1 (MQ=255)
                                                  ttttATCATGTTAATAGGGATGTTAATTTTTTATTTTTATCTAAATAAAAAATTAACATTAAAAATTTTTAATACCAACGTGGTTAATAAg                                          <  1:192789/91‑1 (MQ=255)
                                                   tttATCATGTTAATAGGGATGTTAATTTTTTATTTTTATCTAAATAAAAAATTAACATTAAAAATTTTTAATACCAACGTGGTTAATAAGa                                         >  2:3040971/1‑91 (MQ=255)
                                                       tCATGTTAATAGGGATGTTAATTTTTTATTTTTAGCTAGATAAAAAATGAACATTAAAAATTTTTAATGGCAGCGTGGTTAATAAGAATag                                     <  2:992157/91‑1 (MQ=255)
                                                        cATGTTAATAGGGATGTTAATTTTTTATTTTTATCTAAATAAAAAATTAACATTAAAAATTTTTAATACTAACGTGGTTAATAAGATTaga                                    >  2:530927/1‑91 (MQ=255)
                                                        cATGTTAATAGGGATGTTAATTTTTTATTTTTATCTAAATAAAAAATTAACATTAAAAATTTTTAATACCAACGTGGTTAATAAGAATaga                                    >  2:497049/1‑91 (MQ=255)
                                                        cATGTTAATAGGGATGTTAATTTTTTATTTTTATCTAAATAAAAAATTAACATTAAAAATTTTTAATACCAACGTGGTTAATAAGAATaga                                    >  1:3464095/1‑91 (MQ=255)
                                                         aTGTTAATAGGGATGTTAATTTTTTATTTTTATCTAAATAAAAAATTAACATTAAAAATTTTTAATACCAACGTGGTTAATAAGAATagag                                   >  2:2755391/1‑91 (MQ=255)
                                                         aTGTTAATAGGGATGTTAATTTTTTATTTTTATCTAAATAAAAAATTAACATTAAAAATTTTTAATACCAACGTGGTTAATAAGAATagag                                   <  1:2661890/91‑1 (MQ=255)
                                                            ttAATAGGGATGTTAATTTTTTATTTTTATCTACATAAAAAATTAACATTAAAAATTTTTAATACCAACGTGGTTAATAAGAATAGAGTTg                                <  1:3173176/91‑1 (MQ=255)
                                                            ttAATAGGGATGTTAATTTTTTATTTTTATCTAAATAAAAAATTAACATTAAAAATTTTTAATACCAACGTGGTTAATAAGAATAGAGTTg                                >  1:2215256/1‑91 (MQ=255)
                                                            ttAATAGGGATGTTAATTTTTTATTTTTATCTAAATAAAAAATTAACATTAAAAATTTTTAATACCAACGTGGTTAATAAGAATAGAGTTg                                <  1:2880738/91‑1 (MQ=255)
                                                            ttAATAGGGATGTAAATTTTTTATTTTTATCTAAATAAAAAATTAACATTAAAAATTTTTAATACCAACGTGGTTAATAAGAATAGAGTTg                                <  2:1422351/91‑1 (MQ=255)
                                                             tAATAGGGATGTTAATTTTTTATTTTTATCTAAATAAAAAATTAACATTAAAAATTTTTAATACCAACGTGGTTAATAAGAATAGAGTTGt                               <  1:2367929/91‑1 (MQ=255)
                                                                tAGGGATGTTAATTTTTTATTTTTATCTAAATAAAAAATTAACATTAAAAATTTTTAATACCAACGTGGTTAATAAGAATAGAGTTGTAtt                            <  1:1885543/91‑1 (MQ=255)
                                                                 tttaaTGTTAATTTTTTATTTAGATAAAAATAAAAAATTAACATTAAAAATTTTTAATACCAACGTGGTTAATAAGAATAGAGTTGTATTa                           >  1:1683861/5‑91 (MQ=21)
                                                                  ttaaTGTTAATTTTATATTTCGATAAAAATAAAAAATTAACATTAATAATTTTTAATACCCCCGTGGTTAATAAGACTAGAGTTGTATTaa                          >  2:469308/4‑91 (MQ=18)
                                                                  gggATGTTAATTTTTTATTTTTATCTAAATAAAAAATTAACATTAAAAATTTTTAATACCAACGTGGTTAATAAGAATAGAGTTGTATTaa                          >  2:324440/1‑91 (MQ=255)
                                                                   ggATGTTAATTTTTTATTTTTATCTAAATAAAAAATTAACATTAAAAATTTTTAATACCAACGTGGTTAATAAGAATAGAGTTGTATTAAt                         >  1:1269112/1‑91 (MQ=255)
                                                                   ggATGTTAATTTTTTATTTTTATCTAAATAAAAAATTAACATTAAAAATTTTTAATACCAACGTGGTTAATAAGAATAGAGTTGTATTAAt                         >  1:1888457/1‑91 (MQ=255)
                                                                    gATGTTAATTTTTTATTTTTATCTAAATAAAAAATTAACATTAAAAATTTTTAATACCAACGTGGTTACTAAGAATAGAGTTGTATTAATa                        >  2:3132301/1‑91 (MQ=255)
                                                                    gATGTTAATTTTTTATTTTTATCTAAATAAAAAATTAACATTAAAAATTTTTAATACCAACGTGGTTAATAAGAATAGAGTTGTATTAATa                        >  1:1011792/1‑91 (MQ=255)
                                                                     aTGTTAATTTTTTATTTTTATCTAAATAAAAAATTAACATTAAAAATTTTTAATACCAACGTGGTTAATAAGAATAGAGTTGTATTAATAc                       <  1:433403/91‑1 (MQ=255)
                                                                     aTGTTAATTTTTTATTTTTATCTAAATAAAAAATTAACATTAAAAATTTTTAATACCAACGTGGTTAATAAGAATAGAGTTGTATTAATAc                       >  2:1269145/1‑91 (MQ=255)
                                                                      tGTTAATTTTTTATTTTTATCTAAATAAAAAATTAACATTAAAAATTTTTAATACCAACGTGGTTAATAAGAATAGAGTTGTATTAATAcc                      <  1:3575481/91‑1 (MQ=255)
                                                                       gTTAATTTTTTATTTTTATCTAAATAAAAAATTAACATTAAAAATTTTTAATACCAACGTGGTTAATAAGAATAGAGTTGTATTAATACCg                     <  2:2255438/91‑1 (MQ=255)
                                                                       gTTAATTTTTTATTTTTATCTAAATAAAAAATTAACATTAAAAATTTTTAATACCAACGTGGTTAATAAGAATAGAGTTGTATTAATACCg                     >  2:3058715/1‑91 (MQ=255)
                                                                       gTTAATTTTTTATTTTTATCTAAATAAAAAATTAACATTAAAAATTTTTAATACCAACGTGGTTAATAAGAATAGAGTTGTATTAATACCg                     >  2:3352905/1‑91 (MQ=255)
                                                                         tAATTTTTTATTTTTATCTAAATAAAAAATTAACATTAAAAATTTTTAATACCAACGTGGTTAATAAGAATAGAGTTGTATTAATACCGAc                   >  2:3100658/1‑91 (MQ=255)
                                                                          aaTTTTTTATTTTTATCTAAATAAAAAATTAACATTAAAAATTTTTAATACCAACGTGGTTAATAAGAATAGAGTTGTATTAATACCGACt                  <  1:1012548/91‑1 (MQ=255)
                                                                          aaTTTTTTATTTTTATCTAAATAAAAAATTAACATTAAAAATTTTTAATACCAACGTGGTTAATAAGAATAGAGTTGTATTAATACCGACt                  >  2:1074606/1‑91 (MQ=255)
                                                                           aTTTTTTATTTTTATCTAAATAAAAAATTAACATTAAAAATTTTTAATACCAACGTGGTTAATAAGAATAGAGTTGTATTAATACCGACtt                 <  1:108714/91‑1 (MQ=255)
                                                                           aTTTTTTATTTTTATCTAAATAAAAAATTAACATTAAAAATTTTTAATACCAACGTGGTTAATAAGAATAGAGTTGTATTAATACCGACtt                 <  1:3388473/91‑1 (MQ=255)
                                                                            ttttttATTTTTATCTAAATAAAAAATTAACATTAAAAATTTTTAATACCAACGTGGTTAATAAGAATAGAGTTGTATTATTACCGACttt                >  2:412996/1‑91 (MQ=255)
                                                                            ttttttATTTTTATCTAAATAAAAAATTAACATTAAAAATTTTTAATACCAACGTGGTTAATAAGAATAGAGTTGTATTAATACCGACttt                >  1:2844443/1‑91 (MQ=255)
                                                                            ttttttATTTTTATCTAAATAAAAAATTAACATTAAAAATTTTTAATACCAACGTGGTTAATAAGAATAGAGTTGTATTAATACCGACttt                >  1:1817414/1‑91 (MQ=255)
                                                                             tttttatttttaTCTAAATAAAAAATTAACATTAAAAATTTTTAATACCAACGTGGTTAATAAGAATAGAGTTGTATTAATACCGACTTTc               <  2:2383189/91‑1 (MQ=255)
                                                                                tttttttttttTAAATAAAAAATTAACATTAAAAATTTTTAATACCAACGTGGTTAATAAGAATAGAGTTGTATTAATACCGACTTTCaaa            >  2:1575736/4‑91 (MQ=255)
                                                                                ttatttttatCTAAATAAAAAATTAACATTAAAAATTTTTAATACCAACGTGGTTAATAAGAATAGAGTTGTATTAATACCGACTTTCaaa            >  1:2045467/1‑91 (MQ=255)
                                                                                ttatttttatCTAAATAAAAAATTAACATTAAAAATTTTTAATACCAACGTGGTTAATAAGAATAGAGTTGTATTAATACCGACTTTCaaa            >  1:2789079/1‑91 (MQ=255)
                                                                                 tatttttatCTAAATAAAAAATTAACATTAAAAATTTTTAATACCAACGTGGTTAATAAGAATAGAGTTGTATTAATACCGACTTTCAAAc           <  2:2884751/91‑1 (MQ=255)
                                                                                  atttttatCTAAATAATAAATTAACATTAAAAATTTTTAATACCAACGTGGTTAATAAGAATAGAGTTGTATTAATACCGACTTTCAAACa          <  2:2985018/91‑1 (MQ=255)
                                                                                    ttttatCTAAATAAAAAATTAACATTAAAAATTTTTAATACCAACGTGGTTAATAAGAATAGAGTTGTATTAATACCGACTTTCAAACAAg        <  2:1989525/91‑1 (MQ=255)
                                                                                    ttttatCTAAATAAAAAATTAACATTAAAAATTTTTAATACCAACGTGGTTAATAAGAATAGAGTTGTATTAATACCGACTTTCAAACAAg        <  1:1368765/91‑1 (MQ=255)
                                                                                     tttatttaGATAAAAAATTAACATTAAAAATTTTTAATACCAACGTGGTTAATAAGAATAGAGTTGTATTAATACCGACTTTCAAACAAgg       >  2:197843/1‑91 (MQ=255)
                                                                                     tttatCTAAATAAAAAATTAACATTAAAAATTTTTAATACCAACGTGGTTAATAAGAATAGAGTTGTATTAATACCGACTTTCAAACAAgg       >  1:3749695/1‑91 (MQ=255)
                                                                                      ttatCTAAATAAAAAATTAACATTAAAAATTTTTAATACCAACGTGGTTAATAAGAATAGAGTTGTATTAATACCGACTTTCAAACAAGGa      <  2:1172064/91‑1 (MQ=255)
                                                                                       tatCTAAATAAAAAATTAACATTAAAAATTTTTAATACCAACGTGGTTAATAAGAATAGAGTTGTATTAATACCGACTTTCAAACAAGGAt     >  2:249831/1‑91 (MQ=255)
                                                                                          |                                                                                          
GGTAACTTCTTTAATGTTCAGTTTGAAAAACAACTAGCAATCACTGTATTTTTTATCATGTTAATAGGGATGTTAATTTTTTATTTTTATCTAAATAAAAAATTAACATTAAAAATTTTTAATACCAACGTGGTTAATAAGAATAGAGTTGTATTAATACCGACTTTCAAACAAGGATTGT  >  USA300TCH1516_ALE/659285‑659465

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: