BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	USA300TCH1516_ALE	341,450	G→C	10.2%	D134H (GAT→CAT)	USA300HOU_RS01525 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	USA300TCH1516_ALE	341,450	0	G	C	10.2%	729.6 / 59.3	226	D134H (GAT→CAT)	USA300HOU_RS01525	hypothetical protein
Reads supporting (aligned to +/- strand): ref base G (104/99); new base C (9/14); total (113/113)
Fisher's exact test for biased strand distribution p-value = 3.79e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.05e-01

Only 100 of 229 total aligned reads displayed.

GGGTTTCAAATACGAATTATCAATAAAGACATTGCATGGTGGTCGATTGTACAATCATATCCTTTATTGATGAGACAAGAAAAGT‑CCAATGATTGGGACTTTAAACAAATTGACGATGAAACATTGGAGAACTTAAATAATGAAAGTATCGATACGATTGGGCGTGTTTTAGAAATTGAAA  >  USA300TCH1516_ALE/341360‑341540

gggTTTCAAATACGAATTATCAATAAAGACATTGCATGGTGGTCGATTGTACAATCATATCCTTTATTGATGAGACAAGAAAAGT‑CCAATg                                                                                            >  1:1050158/1‑91 (MQ=255)

gggTTTCAAATACGAATTATCAATAAAGACATTGCATGGTGGTCGATTGTACAATCATATCCTTTATTGATGAGACAAGAAAAGT‑CCAATg                                                                                            <  2:2678618/91‑1 (MQ=255)

 ggTTTCAAATACGAATTATCAATAAAGACATTGCATGGTGGTCGATTGTACAATCATATCCTTTATTGATGAGACAAGAAAAGT‑CCAATGa                                                                                           >  1:3855717/1‑91 (MQ=255)

 ggTTTCAAATACGAATTATCAATAAAGACATTGCATGGTGGTCGATTGTACAATCATATCCTTTATTGATGAGACAAGAAAAGT‑CCAATGa                                                                                           <  1:3775325/91‑1 (MQ=255)

 ggTTTCAAATACGAATTATCAATAAAGACATTGCATGGTGGTCGATTGTACAATCATATCCTTTATTGATGAGACAAGAAAAGT‑CCAATGa                                                                                           <  1:2900469/91‑1 (MQ=255)

 ggTTTCAAATACGAATTATCAATAAAGACATTGCATGGTGGTCGATTGTACAATCATATCCTTTATTGATGAGACAAGAAAAGT‑CCAATGa                                                                                           >  1:1548138/1‑91 (MQ=255)

  gTTTCAAATACGAATTATCAATAAAGACATTGCATGGTGGTCGATTGTACAATCATATCCTTTATTGATGAGACAAGAAAAGT‑CCAATGAt                                                                                          >  1:1021391/1‑91 (MQ=255)

  gTTTCAAATACGAATTATCAATAAAGACATTGCATGGTGGTCGATTGTACAATCATATCCTTTATTGATGAGACAAGAAAAGT‑CCAATGAt                                                                                          <  1:1541430/91‑1 (MQ=255)

   tttCAAATACGAATTATCAATAAAGACATTGCATGGTGGTCGATTGTACAATCATATCCTTTATTGATGAGACAAGAAAAGT‑CCAATGAtt                                                                                         <  2:2904445/91‑1 (MQ=255)

   tttCAAATACGAATTATCAATAAAGACATTGCATGGTGGTCGATTGTACAATCATATCCTTTATTGATGAGACAAGAAAAGT‑CCAATGAtt                                                                                         >  1:2517989/1‑91 (MQ=255)

    ttCAAATACGAATTATCAATAAAGACATTGCATGGTGGTCGATTGTACAATCATATCCTTTATTGATGAGACAAGAAAAAT‑CCAATGATTg                                                                                        >  1:2819505/1‑91 (MQ=255)

     tCAAATACGAATTATCAATAAAGACATTGCATGGTGGTCGATTGTACAATCATATCCTTTATTGATGAGACAAGAAAAGTCCCAATCATTg                                                                                        <  2:29941/91‑1 (MQ=255)

     tCAAATACGAATTATCAATAAAGACATTGCATGGTGGTCGATTGTACAATCATATCCTTTATTGATGAGACAAGAAAAGT‑CCAATGATTgg                                                                                       >  1:560888/1‑91 (MQ=255)

     tCAAATACGAATTATCAATAAAGACATTGCAGGGTGGTCGATTGTACAATCATATCCTTTATTGATGAGACAAGAAAAGTCCCAATCATTg                                                                                        <  2:2562989/91‑1 (MQ=255)

      cAAATACGAATTATCACTAAAGACATTGCATGGTGGTCGATTGTACAATCATATCCTTCATTGATGAGACAAGAAAAGTCCCAATCATTgg                                                                                       <  2:2022077/91‑1 (MQ=255)

      cAAATACGAATTATCAATAAAGACATTGCATGGTGGTCGATTGTACAATCATATCCTTTATTGATGAGACAAGAAAAGTCCCAATCATTgg                                                                                       <  1:374359/91‑1 (MQ=255)

      cAAATACGAATTATCAATAAAGACATTGCATGGTGGTCGATTGTACAATCATATCCTTTATTGATGAGACAAGAAAAGTCCCAATCATTgg                                                                                       <  1:3915778/91‑1 (MQ=255)

      cAAATACGAATTATCAATAAAGACATTGCATGGTGGTCGATTGTACAATCATATCCTTTATTGATGAGACAAGAAAAGTCCCAATCATTgg                                                                                       <  2:1065899/91‑1 (MQ=255)

      cAAATACGAATTATCAATAAAGACATTGCATGGTGGTCGATTGTACAATCATATCCTTTATTGATGAGACAAGAAAAGT‑CCAATGATTggg                                                                                      >  1:2088928/1‑91 (MQ=255)

          tACGAATTATCAATAAAGACATTGCATGGTGGTCGATTGTACAATCATATCCTTTATTGATGAGACAAGAAAAGT‑CCAATGATTGGGACtt                                                                                  <  1:2648416/91‑1 (MQ=255)

             gAATTATCAATAAAGACATGGCATGGTGGTCGATTGTACAATCATATCCTTGATTGATGAGACAAGAAAAGT‑CCAATGATTGGGACTTTaa                                                                               <  2:1301603/91‑1 (MQ=255)

              aaTTATCAATAAAGACATTGCATGGTGGTCGATTGTACAATCATATCCTTTATTGATGAGACAAGAAAAGT‑CCAATGATTGGGACTTTaaa                                                                              <  1:1009124/91‑1 (MQ=255)

              aaTTATCAATAAAGACATTGCATGGTGGTCGATTGTACAATCATATCCTTTATTGATGAGACAAGAAAAGT‑CCAATGATTGGGACTTTaaa                                                                              <  1:3144771/91‑1 (MQ=255)

              aaTTATCAATAAAGACATTGCATGGTGGTCGATTGTACAATCATATCCTTTATTGATGAGACAAGAAAAGT‑CCAATGATTGGGACTTTaaa                                                                              >  1:2690259/1‑91 (MQ=255)

                 tATCAATAAAGACATTGCATGGTGGTCGATTGTACAATCATATCCTTTATTGATGAGACAAGAAAAGT‑CCAATGATTGGGACTTTAAACaa                                                                           >  2:3523472/1‑91 (MQ=255)

                    cAATAAAGACATTGCATGGTGGTCGATTGTACAATCATATCCTTTATTGATGAGACAAGAAAAGT‑CCAATGATTGGGACTTTAAACAAAtt                                                                        <  2:846324/91‑1 (MQ=255)

                     aaTAAAGACATTGCATGGTGGTCGATTGTACAATCATATCCTTTATTGATGAGACAAGAAAAGT‑CCAATGATTGGGACTTTAAACAAATTg                                                                       <  1:2644344/91‑1 (MQ=255)

                     aaTAAAGACATTGCATGGTGGTCGATTGTACAATCATATCCTTTATTGATGAGACAAGAAAAGT‑CCAATGATTGGGACTTTAAACAAATTg                                                                       >  2:3780911/1‑91 (MQ=255)

                     aaTAAAGACATTGCATGGTGGTCGATTGTACAATCATATCCTTTATTGATGAGACAAGAAAAGT‑CCAATGATTGGGACTTTAAACAAATTg                                                                       >  2:2776928/1‑91 (MQ=255)

                      aTAAAGACATTGCATGGTGGTCGATTGTACAATCATATCCTTTATTGATGAGACAAGAAAAGT‑CCAATGATTGGGACTTTAAACAAATTGa                                                                      >  1:1815344/1‑91 (MQ=255)

                      aTAAAGACATTGCATGGTGGTCGATTGTACAATCATATCCTTTATTGATGAGACAAGAAAAGT‑CCAATGATTGGGACTTTAAACAAATTGa                                                                      >  2:1824609/1‑91 (MQ=255)

                          aGACATTGCATGGTGGTCGATTGTACAATCATATCCTTTATTGATGAGACAAGAAAAGT‑CCAATGATTGGGACTTTAAACAAATTGACGAt                                                                  <  1:3082694/91‑1 (MQ=255)

                          aGACATTGCATGGTGGTCGATTGTACAATCATATCCTTTATTGATGAGACAAGAAAAGT‑CCAATGATTGGGACTTTAAACAAATTGACGAt                                                                  <  2:3386462/91‑1 (MQ=255)

                           gACATTGCATGGTGGTCGATTGTACAATCATATCCTTTATTGATGAGACAAGAAAAGT‑CCAATGATTGGGACTTTAAACAAATTGACGATg                                                                 >  1:491320/1‑91 (MQ=255)

                                 gCATGGTGGTCGATTGTACAATCATATCCTTTATTGATGAGACAAGAAAAGT‑CCAATGATTGGGACTTTAAACAAATTGACGATGAAACAt                                                           >  1:676739/1‑91 (MQ=255)

                                 gCATGGTGGTCGATTGTACAATCATATCCTTTATTGATGAGACAAGAAAAGT‑CCAATGATTGGGACTTTAAACAAATTGACGATGAAACAt                                                           >  2:709899/1‑91 (MQ=255)

                                  cATGGTGGTCGATTGTACAATCATATCCTTTATTGATGAGACAAGAAAAGT‑CCAATGATTGGGACTTTAAACAAATTGACGATGAAACAtt                                                          <  2:1467211/91‑1 (MQ=255)

                                    tggtggTCGATTGTACAATCATATCCTTTATTGATGAGACAAGAAAAGT‑CCAATGATTGGGACTTTAAACAAATTGACGATGAAACATTgg                                                        <  1:2872207/91‑1 (MQ=255)

                                    tggtggTCGATTGTACAATCATATCCTTTATTGATGAGACAAGAAAAGT‑CCAATGATTGGGACTTTAAACAAATCGACGATGAAACATTgg                                                        >  2:866157/1‑91 (MQ=255)

                                     ggtggtCGATTGTACAATCATATCCTTTATTGATGAGACAAGAAAAGT‑CCAATGATTGGGACTTTAAACAAATTGACGATGAAACATTGga                                                       >  2:1779417/1‑91 (MQ=255)

                                     ggtggtCGATTGTACAATCATATCCTTTATTGATGAGACAAGAAAAGT‑CCAATGATTGGGACTTTAAACAAATTGACGATGAAACATTGga                                                       >  1:3567696/1‑91 (MQ=255)

                                       tggtCGATTGTACAATCATATCCTTTATTGATGAGGCAAGAAATGT‑CCAATGATTGGCACTTTTATCAAATTCACGATGAATCATTGgtga                                                     >  2:947136/1‑91 (MQ=255)

                                         gtCGATTGTACAATCATATCCTTTATTGATGAGACGAGAAAAGT‑CCGATGATTGGGACTTTAAACAAATTGACGATGAAACATTGGAGAAc                                                   <  2:3778647/91‑1 (MQ=255)

                                         gtCGATTGTACAATCATATCCTTTATTGATGAGACAAGAAAAGT‑CAAATGATTGGGACTTTAAACAAATTGACGATGAAACATTGGAGAAc                                                   >  2:62616/1‑91 (MQ=255)

                                            gATTGTACAATCATATCCTTTATTGATGAGACAAGAAAAGT‑CCAATGATTGGGACTTTAAACAAATTGACGATGAAACATTGGAGAACTTa                                                <  2:2641175/91‑1 (MQ=255)

                                            gATTGTACAATCATATCCTTTATTGATGAGACAAGAAAAGT‑CCAATGATTGGGACTTTAAACAAATTGACGATGAAACATTGGAGAACTTa                                                <  1:2780428/91‑1 (MQ=255)

                                            gATTGTACAATCATATCCTTTATTGATGAGACAAGAAAAGT‑CCAATGATTGGGACTTTAAACAAATTGACGATGAAACATTGGAGAACTTa                                                <  2:1448827/91‑1 (MQ=255)

                                              ttGTACAATCATATCCTTTATTGATGAGACAAGAAAAGT‑CCAATGATTGGGACTTTAAACAAATTGACGATGAAACATTGGAGAACTTaaa                                              >  2:1497008/1‑91 (MQ=255)

                                              ttGTACAATCATATCCTTTATTGATGAGACAAGAAAAGT‑CCAATGATTGGGACTTTAAACAAATTGACGATGAAACATTGGAGAACTTaaa                                              >  2:450924/1‑91 (MQ=255)

                                               tGTACAATCATATCCTTTATTGATGAGACAAGAAAAGT‑CCAATGATTGGGACTTTAAACAAATTGACGATGAAACATTGGAGAACTTAaat                                             <  2:1330801/91‑1 (MQ=255)

                                                gTACAATCATATCCTTTATTGATGAGACAAGAAAAGT‑CCAATGATTGGGACTTTAAACAAATTGACGATGAAACATTGGAGAACTTAaata                                            >  2:2782679/1‑91 (MQ=255)

                                                gTACAATCATATCCTTTATTGATGAGACAAGAAAAGT‑CCAATGATTGGGACTTTAAACAAATTGACGATGAAACATTGGAGAACTTAaata                                            <  1:1377924/91‑1 (MQ=255)

                                                 tACAATCATATCCTTTATTGATGAGACAAGAAAAGT‑CCAATGATTGGGACTTTAAACAAATTGACGATGAAACATTGGAGAACTTAaataa                                           <  2:1634236/91‑1 (MQ=255)

                                                  aCAATCATATCCTTTATTGATGAGACAAGAAAAGT‑CCAATGATTGGGACTTTAAACAAATTGACGATGAAACATTGGAGAACTTAaataat                                          >  2:2830230/1‑91 (MQ=255)

                                                     aTCATATCCTTTATTGATGAGACAAGAAAAGT‑CCAATGATTGGGACTTTAAACAAATTGACGATGAAACATTGGAGAACTTAAATAATGaa                                       >  2:2212668/1‑91 (MQ=255)

                                                     aTCATATCCTTTATTGATGAGACAAGAAAAGT‑CCAATGATTGGGACTTTAAACAAATTGACGATGAAACATTGGAGAACTTAAATAATGaa                                       <  1:3347837/91‑1 (MQ=255)

                                                     aTCATATCCTTTATTGATGAGACAAGAAAAGT‑CCAATGATTGGGACTTTAAACAAATTGACGATGAAACATTGGAGAACTTAAATAATGaa                                       >  2:3269032/1‑91 (MQ=255)

                                                      tCATATCCTTTATTGATGAGACAAGAAAAGT‑CCAATGATTGGGACTTTAAACAAATTGACGATGAAACATTGGAGAACTTAAATAATGaaa                                      <  2:1185390/91‑1 (MQ=255)

                                                         tatCCTTTATTGATGAGACAAGAAAAGT‑CCAATGATTGGGACTTTAAACAAATTGACGATGAAACATTGGAGAACTTAAATAATGAAAGTa                                   >  1:985347/1‑91 (MQ=255)

                                                          atCCTTTATTGATGAGACAAGAAAAGT‑CCAATGATTGGGACTTTAAACAAATTGACGATGAAACATTGGAGAACTTAAATAATGAAAGTAt                                  <  1:3145373/91‑1 (MQ=255)

                                                          atCCTTTATTGATGAGACAAGAAAAGT‑CCAATGATTGGGACTTTAAACAAATTGACGATGAAACATTGGAGAACTTAAATAATGAAAGTAt                                  >  1:2307792/1‑91 (MQ=255)

                                                          atCCTTTATTGATGAGACAAGAAAAGT‑CCAATGATTGGGACTTTAAACAAATTGACGATGAAACATTGGAGAACTTAAATAATGAAAGTAt                                  >  2:2100376/1‑91 (MQ=255)

                                                           tCCTTTATTGATGAGACAAGAAAAGT‑CCAATGATTGGGACTTTAAACAAATTGACGATGAAACATTGGAGAACTTAAATAATGAAAGTATc                                 <  1:1922214/91‑1 (MQ=255)

                                                             cTTTATTGATGAGACAAGAAAAGT‑CTAATGATTGGGACTTTAAACAAATTGACGATGAAACATTGGAGAACTTAAATAATGAAAGTATCGa                               <  1:2674974/91‑1 (MQ=255)

                                                             cTTTATTGATGAGACAAGAAAAGT‑CCAATGATTGGGACTTTAAACAAATTGACGATGAAACATTGGAGAACTTAAATAATGAAAGTATCGa                               >  1:1142891/1‑91 (MQ=255)

                                                             cTTTATTGATGAGACAAGAAAAGT‑CCAATGATTGGGACTTTAAACAAATTGACGATGAAACATTGGAGAACTTAAATAATGAAAGTATCGa                               <  1:1084574/91‑1 (MQ=255)

                                                               ttATTGATGAGACAAGAAAAGT‑CCAATGATTGGGACTTTAAACAAATTGACGATGAAACATTGGAGAACTTAAATAATGAAAGTATCGATa                             >  1:2474479/1‑91 (MQ=255)

                                                               ttATTGATGAGACAAGAAAAGT‑CCAATGATTGGGACTTTAAACAAATTGACGATGAAACATTGGAGAACTTAAATAATGAAAGTATCGATa                             >  1:2474083/1‑91 (MQ=255)

                                                                tATTGATGAGACAAGAAAAGT‑CCAATGATTGGGACTTTAAACAAATTGACGATGAAACATTGGAGAACTTAAATAATGAAAGTATCGATAc                            >  2:1700380/1‑91 (MQ=255)

                                                                 aTTGATGAGACAAGAAAAGT‑CCAATGATTGGGACTTTAAACAAATTGACGATGAAACATTGGAGAACTTAAATAATGAAAGTATCGATACg                           <  1:1701126/91‑1 (MQ=255)

                                                                  ttGATGAGACAAGAAAAGT‑CCAATGATTGGGACTTTAAACAAATTGACGATGAAACATTGGAGAACTTAAATAATGAAAGTATCGATACGa                          >  1:3355990/1‑91 (MQ=255)

                                                                   tgatgaGACAAGAAAAGT‑CCAATGATTGGGACTTTAAACAAATTGACGATGAAACATTGGAGAACTTAAATAATGAAAGTATCGATACGAt                         <  2:2085244/91‑1 (MQ=255)

                                                                     atgaGACAAGAAAAGT‑CCAATGATTGGGACTTTAAACAAATTGACGATGAAACATTGGAGAACTTAAATAATGAAAGTATCGATACGATTg                       >  2:3119905/1‑91 (MQ=255)

                                                                       gagcCAAGAAAAGT‑CCAATGATTGCGACTTTAAACAAATAGACGATGAAAAATTGGAGAACTTAAATAATGAAAGTATCGCTACGACTggg                     >  2:124270/1‑91 (MQ=255)

                                                                       gagaCAAGAAAAGT‑CCAATGATTGGGACTTTAAACAAGTTGACGATGAAACGTTGGAGAACTTAAGTAATGAAAGTATCGATACGATTggg                     >  1:3677818/1‑91 (MQ=255)

                                                                       gagaCAAGAAAAGT‑CCAATGATTGGGACTTTAAACAAATTGACGATGAAACATTGGAGAACTTAAATAATGAAAGTATCGATACGATTggg                     <  1:544669/91‑1 (MQ=255)

                                                                       gagaCAAGAAAAGT‑CCAATGATTGGGACTTTAAACAAATTGACGATGAAACATTGGAGAACTTAAATAATGAAAGTATCGATACGATTggg                     <  1:132662/91‑1 (MQ=255)

                                                                       gagaCAAAAAAAGT‑CCAATGATTGGGACTTTAAACAAATTGACGATGAAACATTGGAGAACTTAAATAATGAAAGTATCGATACGATTggg                     <  2:1710375/91‑1 (MQ=255)

                                                                        agaCAAAAAAAGT‑CCAATGATTGGGACTTTAAACAAATTGACGATGAAACATTGGAGAACTTAAATAATGAAAGTATCGATACGATTGGGc                    <  2:1651009/91‑1 (MQ=255)

                                                                          aCAAGAAAAGT‑CCAATGATTGGGACTTTAAACAAATTGACGATGAAACATTGGAGAACTTAAATAATGAAAGTATCGATACGATTGGGCgt                  <  1:3638888/91‑1 (MQ=255)

                                                                           cAAGAAAAGT‑CCAATGATTGGGACTTTAAACAAATTGACGATGAAACATTGGAGAACTTAAATAATGAAAGTATCGATACGATTGGGCgtg                 <  2:2685424/91‑1 (MQ=255)

                                                                             aGAAAAGT‑CCAATGATTGGGACTTTAAACAAATTGACGATGAAACATTGGAGAACTTAAATAATGAAAGTATCGATACGATTGGGCGTGtt               <  1:3021842/91‑1 (MQ=255)

                                                                             aGAAAAGT‑CCAATGATTGGGACTTTAAACAAATTGACGATGAAACATTGGAGAACTTAAATAATGAAAGTATCGATACGATTGGGCGTGtt               <  1:1222860/91‑1 (MQ=255)

                                                                              gAAAAGT‑CCAATGATTGGGACTTTAAACAAATTGACGATGAAACATTGGAGAACTTAAATAATGAAAGTATCGATACGATTGGGCGTGttt              <  2:1618270/91‑1 (MQ=255)

                                                                              gAAAAGT‑CCAATGATTGGGACTTTAAACAAATTGACGATGAAACATTGGAGAACTTAAATAATGAAAGTATCGATACGATTGGGCGTGttt              <  1:3763287/91‑1 (MQ=255)

                                                                              gAAAAGT‑ACAATGATTGGGACTTTAAACAAATTGACTATGAAACATTCGAGAACTTAAATAATGAAACTATCGATAGGATTTGGCGTGttt              >  2:3672743/1‑91 (MQ=255)

                                                                               aaaaGT‑CCAATGATTGGGACTTTAAACAAATTGACGATGAAACATTGGAGAACTTAAATAATGAAAGTATCGATACGATTGGGCGTGtttt             <  1:365201/91‑1 (MQ=255)

                                                                                aaaGT‑CCAATGATTGGGACTTTAAACAAATTGACGATGAAACATTGGAGAACTTAAATAATGAAAGTATCGATACGATTGGGCGTGTTTTa            >  1:2579965/1‑91 (MQ=255)

                                                                                  aGT‑CCAATGATTGGGACTTTAAACAAATTGACGATGAAACATTGGAGAACTTAAATAATGAAAGTATCGATACGATTGGGCGTGTTTTAGa          >  2:216080/1‑91 (MQ=255)

                                                                                  aGT‑CCAATGATTGGGACTTTAAACAAATTGACGATGAAACATTGGAGAACTTAAATAATGAAAGTATCGATACGATTGGGCGTGTTTTAGa          >  2:238884/1‑91 (MQ=255)

                                                                                    t‑CCAATGATTGGGACTTTAAACAAATTGACGATGAAACATTGGAGAACTTAAATAATGAAAGTATCGATACGATTGGGCGTGTTTTAGaaa        >  2:675147/1‑91 (MQ=255)

                                                                                     cccAATCATT‑GGACTTTAAACAAATTGACGATGAAACATTGGAGAACTTAAATAATGAAAGTATCGATACGATTGGGCGTGTTTTAGAAAt       >  1:1161970/2‑91 (MQ=255)

                                                                                      ccAATGATTGGGACTTTAAACAAATTGACGATGAAACATTGGAGAACTTAAATAATGAAAGTATCGATACGATTGGGCGTGTTTTAGAAAt       <  1:257378/91‑1 (MQ=255)

                                                                                      ccAATCATT‑GGACTTTAAACAAATTGACGATGAAACATTGGAGAACTTAAATAATGAAAGTATCGATACGATTGGGCGTGTTTTAGAAAtt      >  1:605501/1‑91 (MQ=255)

                                                                                      ccAATCATT‑GGACTTTAAACAAATTGACGATGAAACATTGGAGAACTTAAATAATGAAAGTATCGATACGATTGGGCGTGTTTTAGAAAtt      >  2:1347635/1‑91 (MQ=255)

                                                                                      ccAATCATT‑GGACTTTAAACAAATTGACGATGAAACATTGGAGAACTTAAATAATGAAAGTATCGATACGATTGGGCGTGTTTTAGAAAtt      >  1:3591789/1‑91 (MQ=255)

                                                                                        aaTGATTGGGACTTTAAACAAATTGACGATGAAACATTGGAGAACTTAAATAATGAATGTATCGATACGATTGGGCGTGTTTTAGAAATTg     >  1:3245762/1‑91 (MQ=255)

                                                                                         aTGATTGGGACTTTAAACAAATTGACGATGAAACATTGGAGAACTTAAATAATGAAAGTATCGATACGATTGGGCGTGTTTTAGAAATTGa    <  1:2011861/91‑1 (MQ=255)

                                                                                          tGATTGGGACTTTAAACAAATTGACGATGAAACATTGGAGAACTTAAATAATGAAAGTATCGATACGATTGGGCGTGTTTTAGAAATTGaa   >  1:217087/1‑91 (MQ=255)

                                                                                           gATTGGGACTTTAAACAAATTGACGATGAAACATTGGAGAACTTAAATAATGAAAGTATCGATACGATTGGGCGTGTTTTAGAAATTGaaa  >  2:2161153/1‑91 (MQ=255)

GGGTTTCAAATACGAATTATCAATAAAGACATTGCATGGTGGTCGATTGTACAATCATATCCTTTATTGATGAGACAAGAAAAGT‑CCAATGATTGGGACTTTAAACAAATTGACGATGAAACATTGGAGAACTTAAATAATGAAAGTATCGATACGATTGGGCGTGTTTTAGAAATTGAAA  >  USA300TCH1516_ALE/341360‑341540

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: `‑`