BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	USA300TCH1516_ALE	1,465,105	C→G	7.8%	E18Q (GAA→CAA)	folA ←	Dihydrofolate reductase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	USA300TCH1516_ALE	1,465,105	0	C	G	7.8%	790.4 / 45.6	230	E18Q (GAA→CAA)	folA	Dihydrofolate reductase
Reads supporting (aligned to +/- strand): ref base C (99/113); new base G (13/5); total (112/118)
Fisher's exact test for biased strand distribution p-value = 4.89e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.84e-01

Only 100 of 238 total aligned reads displayed.

ATGTCTTACGACCCATTACTAAAGTATGACCTGTTGATAATTTTTTAACATGCTTCAAATCATTTGGTAGGTGCCAAGGTAATTGATTTTCAAAACCAATTACTCGTTGCAAGTCATGTGCAACTAGAATGGATAAAGTCATAATTATCCTCCTTCTTCTATCATTTCATTTTTTATTACT  >  USA300TCH1516_ALE/1465015‑1465195

aTGTCTTACGACCCATTACTAAAGTATGACCTGTTGATAATTTTTTAACATGCTTCAAATCATTTGGTAGGTGCCAAGGTAATTGATTTTc                                                                                            <  1:2913456/91‑1 (MQ=255)

aTGTCTTACGACCCATTACTAAAGTATGACCTGTTGATAATTTTTTAACATGCTTCAAATCATTTGGTAGGTGCCAAGGTAATTGATTTTc                                                                                            <  1:1887578/91‑1 (MQ=255)

aTGTCTTACGACCCATTACTAAAGTATGACCTGTTGATAATTTTTTAACATGCTTCAAATCATTTGGTAGGTGCCAAGGTAATTGATTTTc                                                                                            <  1:3842138/91‑1 (MQ=255)

aTGTCTTACGACCCATTACTAAAGTATGACCTGTTGATAATTTTTTAACATGCTTCAAATCATTTGGTAGGTGCCAAGGTAATTGATTTTc                                                                                            <  1:753642/91‑1 (MQ=255)

aTGTCTTACGACCCATTACTAAAGTATGACCTGTTGATAATTTTTTAACATGCTTCAAATCATTTGGTAGGTGCCAAGGTAATTGATTTTc                                                                                            <  1:1262869/91‑1 (MQ=255)

 tGTCTTACGACCCATTACTAAAGTATGACCTGTTGATAATTTTTTAACATGCTTCAAATCATTTGGTAGGTGCCAAGGTAATTGATTTTCa                                                                                           <  1:1192192/91‑1 (MQ=255)

   tCTTACGACCCATTACTAAAGTATGACCTGTTGATAATTTTTTAACATGCTTCAAATCATTTGGTAGGTGCCAAGGTAATTGATTTTCaaa                                                                                         >  1:2554114/1‑91 (MQ=255)

    cTTACGACCCATTACTAAAGTATGACCTGTTGATAATTTTTTAACATGCTTCAAATCATTTGGTAGGTGCCAAGGTAATTGATTTTCaaaa                                                                                        <  1:1471387/91‑1 (MQ=255)

     ttACGACCCATTACTAAAGTATGACCTGTTGATAATTTTTTAACATGCTTCAAATCATTTGGTAGGTGCCAAGGTAATTGATTTTCAAAAc                                                                                       <  1:1369984/91‑1 (MQ=255)

       aCGACCCATTACTAAAGTATGACCTGTTGATAATTTTTTAACATGCTTCAAATCATTTGGTAGGTGCCAAGGTAATTGATTTTCAAAACCa                                                                                     <  2:994334/91‑1 (MQ=255)

       aCGACCCATTACTAAAGTATGACCTGTTGATAATTTTTTAACATGCTTCAAATCATTTGGTAGGTGCCAAGGTAATTGATTTTCAAAACCa                                                                                     <  2:3144242/91‑1 (MQ=255)

        cGACCCATTACTAAAGTATGACCTGTTGATAATTTTTTAACATGCTTCAAATCATTTGGTAGGTGCCAAGGTAATTGGTTTTGAAAATCaa                                                                                    <  2:2406653/91‑1 (MQ=255)

        cGACCCATTACTAAAGTATGACCTGTTGATAATTTTTTAACATGCTTCAAATCATTTGGTAGGTGCCAAGGTAATTGATTTTCAAAACCaa                                                                                    >  2:1343920/1‑91 (MQ=255)

        cGACCCATTACTAAAGTATGACCTGTTGATAATTTTTTAACATGCTTCAAATCATTTGGTAGGTGCCAAGGTAATTGATTTTCAAAACCaa                                                                                    <  1:1855930/91‑1 (MQ=255)

        cGACCCATTACTAAAGTATGACCTGTTGATAATTTTTTAACATGCTTCAAATCATTTGGCAGGTGCCAAGGTAATTGATTTTCAAAACCaa                                                                                    >  2:3241797/1‑91 (MQ=255)

         gACCCATTACTAAAGTATGACCTGTTGATAATTTTTTAACATGCTTCAAATCATTTGGTAGGTGCCAAGGTAATTGATTTTCAAAACCAAt                                                                                   <  1:3163462/91‑1 (MQ=255)

         gACCCATTACTAAAGTATGACCTGTTGATAATTTTTTAACATGCTTCAAATCATTTGGTAGGTGCCAAGGTAATTGATTTTCAAAACCAAt                                                                                   >  1:594976/1‑91 (MQ=255)

          aCCCATTACTAAAGTATGACCTGTTGATAATTTTTTAACATGCTTCAAATCATTTGGTAGGTGCCAAGGTAATTGGTTTTGAAAATCAAtt                                                                                  <  2:1262105/91‑1 (MQ=255)

             cATTACTAAAGTATGACCTGTTGATAATTTTTTAACATGCTTCAAATCATTTGGTAGGTGCCAAGGTAATTGATTTTCAAAACCAATTACt                                                                               <  1:284302/91‑1 (MQ=255)

                 aCTAAAGTATGACCTGTTGATAATTTTTTAACATGCTTCAAATCATTTGGTAGGTGCCAAGGTAATTGATTTTCAAAACCAATTACTCGtt                                                                           <  2:3874379/91‑1 (MQ=255)

                  cTAAAGTATGACCTGTTGATAATTTTTTAACATGCTTCAAATCATTTGGTAGGTGCCAAGGTAATTGATTTTCAAAACCAATTACTCGTTg                                                                          <  1:3759304/91‑1 (MQ=255)

                  cTAAAGTATGACCTGTTGATAATTTTTTAACATGCTTCAAATCATTTGGTAGGTGCCAAGGTAATTGATTTTCAAAACCAATTACTCGTTg                                                                          >  1:3213686/1‑91 (MQ=255)

                   tAAAGTATGACCTGTTGATAATTTTTTAACATGCTTCAAATCATTTGGTAGGTGCCAAGGTAATTGATTTTCAAAACCAATTACTCGTTGc                                                                         >  2:3034875/1‑91 (MQ=255)

                   tAAAGTATGACCTGTTGATAATTTTTTAACATGCTTCAAATCATTTGGTAGGTGCCAAGGTAATTGATTTTCAAAACCAATTACTCGTTGc                                                                         <  2:506288/91‑1 (MQ=255)

                   caaaGTATGACCTGTTGATAATTTTTTAACATGCTTCAAATCATTTGGTAGGTGCCAAGGTAATTGATTTTCAAAACCAATTACTCGTTGc                                                                         <  2:3845781/90‑1 (MQ=255)

                    aaaGTATGACCTGTTGATAATTTTTTAACATGCTTCAAATCATTTGGTAGGTGCCAAGGTAATTGATTTTCAAAACCAATTACTCGTTGCa                                                                        <  2:224494/91‑1 (MQ=255)

                    aaaGTATGACCTGTTGATAATTTTTTAACATGCTTCAAATCATTTGGTAGGTGCCAAGGTAATTGATTTTCAAAACCAATTACTCGTTGCa                                                                        >  1:3202468/1‑91 (MQ=255)

                     aaGTATGACCTGTTGATAATTTTTTAACATGCTTCAAATCATTTGGTAGGTGCCAAGGTAATTGATTTTCAAAACCAATTACTCGTTGCaa                                                                       <  1:3448712/91‑1 (MQ=255)

                     aaGTATGACCTGTTGATAATTTTTTAACATGCTTCAAATCATTTGGTAGGTGCCAAGGTAATTGATTTTCAAAACCAATTACTCGTTGCaa                                                                       <  1:1007399/91‑1 (MQ=255)

                     aaGTATGACCTGTTGATAATTTTTTAACATGCTTCAAATCATTTGGTAGGTGCCAAGGTAATTGATTTTCAAAACCAATTACTCGTTGCaa                                                                       <  1:911535/91‑1 (MQ=255)

                       gTATGACCTGTTGATAATTTTTTAACATGCTTCAAATCATTTGGTAGGTGCCAAGGTAATTGATTTTCAAAACCAATTACTCGTTGCAAGt                                                                     >  1:581768/1‑91 (MQ=255)

                       gTATGACCTGTTGATAATTTTTTAACATGCTTCAAATCATTTGGTAGGTGCCAAGGTAATTGATTTTCAAAACCAATTACTCGTTGCAAGt                                                                     >  2:3429883/1‑91 (MQ=255)

                       gTATGACCTGTTGATAATTTTTTAACATGCTTCAAATCATTTGGTAGGTGCCAAGGTAATTGATTTTCAAAACCAATTACTCGTTGCAAGt                                                                     <  1:1309408/91‑1 (MQ=255)

                         aTGACCTGTTGATAATTTTTTAACATGCTTCAAATCATTTGGTAGGTGCCAAGGTAATTGATTTTCAAAACCAATTACTCGTTGCAAGTCa                                                                   >  2:2109150/1‑91 (MQ=255)

                           gACCTGTTGATAATTTTTTAACATGCTTCAAATCATTTGGTAGGTGCCAAGGTAATTGATTTTCAAAACCAATTACTCGTTGCAAGTCAtg                                                                 >  2:2280070/1‑91 (MQ=255)

                            aCCTGTTGATAATTTTTTAACATGCTTCAAATCATTTGGTAGGTGCCAAGGTAATTGATTTTCAAAACCAATTACTCGTTGCAAGTCAtgt                                                                >  2:1993976/1‑91 (MQ=255)

                              cTGTTGATAATTTTTTAACATGCTTCAAATCATTTGGTAGGTGCCAAGGTAATTGATTTTCGAAACCAATTACTCGTTGCAAGTCATGTGc                                                              >  1:3114883/1‑91 (MQ=255)

                              cTGTTGATAATTTTTTAACATGCTTCAAATCATTTGGTAGGTGCCAAGGTAATTGATTTTCAAAACCAATTACTCGTTGCAAGTCATGTGc                                                              >  2:1809039/1‑91 (MQ=255)

                              cTGTTGATAATTTTTTAACATGCTTCAAATCATTTGGTAGGTGCCAAGGTAATTGATTTTCAAAACCAATTACTCGTTGCAAGTCATGTGc                                                              <  2:3760244/91‑1 (MQ=255)

                              cTGTTGATAATTTTTTAACATGCTTCAAATCATTTGGTAGGTGCCAAGGTAATTGATTTTCAAAACCAATTACTCGTTGCAAGTCATGTGc                                                              <  1:3878964/91‑1 (MQ=255)

                                gTTGATAATTTTTTAACATGCTTCAAATCATTTGGTAGGTGCCAAGGTAATTGATTTTCAAAACCAATTACTCGTTGCAAGTCATGTGCaa                                                            <  2:1126878/91‑1 (MQ=255)

                                   gATAATTTTTTAACATGCTTCAAATCATTTGGTAGGTGCCAAGGTAATTGATTTTCAAAACCAATTACTCGTTGCAAGTCATGTGCAACTa                                                         <  1:1810584/91‑1 (MQ=255)

                                    aTAATTTTTTAACATGCTTCAAATCATTTGGTAGGTGCCAAGGTAATTGATTTTCAAAACCAATTACTCGTTGCAAGTCATGTGCAACTAg                                                        <  1:1702142/91‑1 (MQ=255)

                                     tAATTTTTTAACATGCTTCAAATCATTTGGTAGGTGCCAAGGTAATTGATTTTCAAAACCAATTACTCGTTGCAAGTCATGTGCAACTAGa                                                       <  2:1029359/91‑1 (MQ=255)

                                      aaTTTTTTAACATGCTTCAAATCATTTGGTAGGTGCCAAGGTAATTGATTTTCAAAACCAATTACTCGTTGCAAGTCATGTGCAACTAGaa                                                      >  1:979548/1‑91 (MQ=255)

                                      aaTTTTTTAACATGCTTCAAATCATTTGGTAGGTGCCAAGGTAATTGATTTTCAAAACCAATTACTCGTTGCAAGTCATGTGCAACTAGaa                                                      >  1:1034484/1‑91 (MQ=255)

                                      aaTTTTTTAACATGCTTCAAATCATTTGGTAGGTGCCAAGGTAATTGATTTTCAAAACCAATTACTCGTTGCAAGTCATGTGCAACTAGaa                                                      >  1:1034624/1‑91 (MQ=255)

                                       aTTTTTTAACATGCTTCAAATCATTTGGTAGGTGCCAAGGTAATTGATTTTCAAAACCAATTACTCGTTGCAAGTCATGTGCAACTAGAAt                                                     <  2:2458801/91‑1 (MQ=255)

                                         tttttAACATGCTTCAAATCATTTGGTAGGTGCCAAGGTAATTGATTTTCAAAACCAATTACTCGTTGCAAGTCATGTGCAACTAGAATgg                                                   >  1:2663819/1‑91 (MQ=255)

                                          ttttAACATGCTTCAAATCATTTGGTAGGTGCCAAGGTAATTGATTTTCAAAACCAATTACTCGTTGCAAGTCATGTGCAACTAGAATGGa                                                  >  1:139208/1‑91 (MQ=255)

                                          ttttAACATGCTTCAAATCATTTGGTAGGTGCCAAGGTAATTGATTTTCAAAACCAATTACTCGTTGCAAGTCATGTGCAACTAGAATGGa                                                  <  2:1636294/91‑1 (MQ=255)

                                            ttAACATGCTTCAAATCATTTGGTAGGTGCCAAGGTAATTGATTTTCAAAACCAATTACTCGTTGCAAGTCATGTGCAACTAGAATGGATa                                                >  1:46825/1‑91 (MQ=255)

                                             tAACATGCTTCTAATCATTTGGTAGGTGCCAAGGTAATTGATTTTCAAAACCAATTACTGGTTGCAAGTCATGTGCAACTAGAATGGATaa                                               <  2:2962040/91‑1 (MQ=255)

                                             tAACATGCTTCAAATCATTTGGTAGGTGCCAAGGTAATTGATTTTCAAAACCAATTACTCGTTGCAAGTCATGTGCAACTAGAATGGATaa                                               >  2:2357968/1‑91 (MQ=255)

                                             tAACATGCTTCAAATCATTTGGTAGGTGCCAAGGTAATTGATTTTCAAAACCAATTACTCGTTGCAAGTCATGTGCAACTAGAATGGATaa                                               >  2:1883277/1‑91 (MQ=255)

                                                 aTGCTTCAAATCATTTGGTAGGTGCCAAGGTAATTGATTTTCAAAACCAATTACTCGTTGCAAGTCATGTGCAACTAGAATGGATAAAGTc                                           >  1:2668164/1‑91 (MQ=255)

                                                       cAAATCATTTGGTAGGTGCCAAGGTAATTGATTTTCAAAACCAATTACTCGTTGCAAGTCATGTGCAACTAGAATGGATAAAGTCATAAtt                                     >  1:3595620/1‑91 (MQ=255)

                                                        aaaTCATTTGGTAGGTGCCAAGGTAATTGATTTTCAAAACCAATTACTCGTTGCAAGTCATGTGCAACTAGAATGGATAAAGTCATAATTa                                    <  1:1859607/91‑1 (MQ=255)

                                                          aTCATTTGGTAGGTGCCAAGGTAATTGATTTTCAAAACCAATTACTCGTTGCAAGTCATGTGCAACTAGAATGGATAAAGTCATAATTAtc                                  <  1:155760/91‑1 (MQ=255)

                                                          aTCATTTGGTAGGTGCCAAGGTAATTGATTTTCAAAACCAATTACTCGTTGCAAGTCATGTGCAACTAGAATGGATAAAGTCATAATTAtc                                  <  2:25248/91‑1 (MQ=255)

                                                          aTCATTTGGTAGGGGCCAAGGTAATTGCTTATCACAACCAATTACTCGTGGCAAGTCGCGTGCAACTAGAATGGATACAGCCATAATTAtc                                  >  2:1190998/1‑91 (MQ=255)

                                                            cATTTGGTAGGTGCCAAGGTAATTGATTTTCAAAACCAATTACTCGTTGCAAGTCATGTGCAACTAGAATGGATAAAGTCATAATTAtcct                                >  2:1052042/1‑91 (MQ=255)

                                                             aTTTGGTAGGTGCCAAGGTAATTGATTTTCAAAACCAATTACTCGTTGCAAGTCATGTGCAACTAGAATGGATAAAGTCATAATTAtcctc                               >  1:1785130/1‑91 (MQ=255)

                                                              tttGGTAGGTGCCAAGGTAATTGATTTTCAAAACCAATTACTCGTTGCAAGTCATGTGCAACTAGAATGGATAAAGTCATAATTAtcctcc                              <  2:1089326/91‑1 (MQ=255)

                                                              tttGGTAGGTGCCAAGGTAATTGATTTTCAAAACCAATTACTCGTTGCAAGTCATGTGCAACTAGAATGGATAAAGTCATAATTAtcctcc                              <  1:623995/91‑1 (MQ=255)

                                                               ttGGTAGGTGCCAAGGTAATTGATTTTCAAAACCAATTACTCGTTGCAAGTCATGTGCAACTAGAATGGATAAAGTCATAATTATcctcct                             >  2:3494415/1‑91 (MQ=255)

                                                                tGGTAGGTGCCAAGGTAATTGATTTTCAAAACCAATTACTCGTTGCAAGTCATGTGCAACTAGAATGGATAAAGTCATAATTATCCTCctt                            <  1:3570388/91‑1 (MQ=255)

                                                                 ggTAGGTGCCAAGGTAATTGATTTTCAAAACCAATTACTCGTTGCAAGTCATGTGCAACTAGAATGGATAAAGTCATAATTATCCTCcttc                           <  1:2446215/91‑1 (MQ=255)

                                                                    aGGTGCCAAGGTAATTGATTTTCAAAACCAATTACTCGTTGCAAGTCATGTGCAACTAGAATGGATAAAGTCATAATTATCCTCCttcttc                        >  2:780349/1‑91 (MQ=255)

                                                                     ggTGCCAAGGTAATTGATTTTCAAAACCAATTACTCGTTGCAAGTCATGTGCAACTAGAATGGATAAAGTCATAATTATCCTCCTtcttct                       <  1:2154893/91‑1 (MQ=255)

                                                                       tGCCAAGGTAATTGATTTTCAAAACCAATTACTCGTTGCAAGTCATGTGCAACTAGAATGGATAAAGTCATAATTATCCTCCTTCTTCTAt                     >  1:2335069/1‑91 (MQ=255)

                                                                       tGCCAAGGTAATTGATTTTCAAAACCAATTACTCGTTGCAAGTCATGTGCAACTAGAATGGATAAAGTCATAATTATCCTCCTTCTTCTAt                     <  1:603547/91‑1 (MQ=255)

                                                                       tGCCAAGGTAATTGATTTTCAAAACCAATTACTCGTTGCAAGTCATGTGCAACTAGAATGGATAAAGTCATAATTATCCTCCTTCTTCTAt                     <  1:1480857/91‑1 (MQ=255)

                                                                         ccAAGGTAATTGATTTTCAAAACCAATTACTCGTTGCAAGTCATGTGCAACTAGAATGGATAAAGTCATAATTATCCTCCTTCTTCTAtca                   <  1:2605251/91‑1 (MQ=255)

                                                                          cAAGGTAATTGATTTTCAAAACCAATTACTCGTTGCAAGTCATGTGCAACTAGAATGGATAAAGTCATAATTATCCTCCTTCTTCTAtcat                  <  1:2733354/91‑1 (MQ=255)

                                                                           aaGGTAATTGATTTTCAAAACCAATTACTCGTTGCAAGTCATGTGCAACTAGAATGGATAAAGTCATAATTATCCTCCTTCTTCTAtcatt                 <  2:1456075/91‑1 (MQ=255)

                                                                           aaGGTAATTGATTTTCAAAACCAATTACTCGTTGCAAGTCATGTGCAACTAGAATGGATAAAGTCATAATTATCCTCCTTCTTCTAtcatt                 <  1:1487652/91‑1 (MQ=255)

                                                                            aGGTAATTGATTTTCAAAACCAATTACTCGTTGCAAGTCATGTGCAACTAGAATGGATAAAGTCATAATTATCCTCCTTCTTCTAtcattt                >  1:2257038/1‑91 (MQ=255)

                                                                            aGGTAATTGATTTTCAAAACCAATTACTCGTTGCAAGTCATGTGCAACTAGAATGGATAAAGTCATAATTATCCTCCTTCTTCTAtcattt                <  2:336613/91‑1 (MQ=255)

                                                                             ggTAATTGATTTTCAAAACCAATTACTCGTTGCAAGTCATGTGCAACTAGAATGGATAAAGTCATAATTATCCTCCTTCTTCTAtcatttc               <  2:3213686/91‑1 (MQ=255)

                                                                              gTAATTGGTTTTGAAAATCAATTACTCGTTGCAAGTCATGTGCAACTAGAATGGATAAAGTCATAATTATCCTCCTTCTTCTAtcatttca              >  2:66222/1‑91 (MQ=39)

                                                                              gTAATTGATTTTCAAAACCAATTACTCGTTGCAAGTCATGTGCAACTAGAATGGATAAAGTCATAATTATCCTCCTTCTTCTAtcatttca              <  1:1669768/91‑1 (MQ=255)

                                                                               tAATTGGTTTTGAAAATCAATTACTCGTTGCAAGTCATGTGCAACTAGAATGGATAAAGTCATAATTATCCTCCTTCTTCTAtcatttcat             >  2:951316/1‑91 (MQ=39)

                                                                               tAATTGATTTTCAAAACCAATTACTCGTTGCAAGTCATGTGCAACTAGAATGGATAAAGTCATAATTATCCTCCTTCTTCTAtcatttcat             >  1:1781637/1‑91 (MQ=255)

                                                                               tAATTGATTTTCAAAACCAATTACTCGTTGCAAGTCATGTGCAACTAGAATGGATAAAGTCATAATTATCCTCCTTCTTCTAtcatttcat             <  2:1415726/91‑1 (MQ=255)

                                                                                aaTTGATTTTCAAAACCAATTACTCGTTGCAAGTCATGTGCAACTAGAATGGATAAAGTCATAATTATCCTCCTTCTTCTAtcatttcatt            >  2:2764621/1‑91 (MQ=255)

                                                                                 aTTGATTTTCAAAACCAATTACTCGTTGCAAGTCATGTGCAACTAGAATGGATAAAGTCATAATTATCCTCCTTCTTCTATcatttcattt           >  1:2982428/1‑91 (MQ=255)

                                                                                 aTTGATTTTCAAAACCAATTACTCGTTGCAAGTCATGTGCAACTAGAATGGATAAAGTCATAATTATCCTCCTTCTTCTATcatttcattt           <  2:3319208/91‑1 (MQ=255)

                                                                                 aTTGATTTTCAAAACCAATTACTCGTTGCAAGTCATGTGCAACTAGAATGGATAAAGTCATAATTATCCTCCTTCTTCTATcatttcattt           <  1:2803354/91‑1 (MQ=255)

                                                                                  ttggttttGAAAATCAATTACTCGTTGCAAGTCATGTGCAACTAGAATGGATAAAGTCATAATTATGCTCCTTCTTCTATCATTTCAtttt          >  2:3597723/5‑91 (MQ=255)

                                                                                  ttggttttGAAAATCAATTACTCGTTGCAAGTCATGTGCAACTAGAATGGATAAAGTCATAATTATCCTCCTTCTTCTATCATTTCAtttt          >  1:2382423/5‑91 (MQ=255)

                                                                                   tggttttGAAAATCAATTACTCGTTGCAAGTCATGTGCAACTAGAATGGATAAAGTCATAATTATCCTCCTTCTTCTATCATTTCAttttt         >  2:1221084/4‑91 (MQ=255)

                                                                                   tGATTTTCAAAACCAATTACTCGTTGCAAGTCATGTGCAACTAGAATGGATAAAGTCATAATTATCCTCCTTCTTCTATCATTTCAttttt         >  1:989045/1‑91 (MQ=255)

                                                                                    ggttttGAAAATCAATTACTCGTTGCAAGTCATGTGCAACTAGAATGGATAAAGTCATAATTATCCTCCTTCTTCTATCATTTCAtttttt        >  2:234151/3‑91 (MQ=255)

                                                                                    gATTTTCAAAACCAATTACTCGTTGCAAGTCATGTGCAACTAGAATGGATAAAGTCATAATTATCCTCCTTCTTCTATCATTTCAtttttt        >  2:1668393/1‑91 (MQ=255)

                                                                                    gATTTTCAAAACCAATTACTCGTTGCAAGTCATGTGCAACTAGAATGGATAAAGTCATAATTATCCTCCTTCTTCTATCATTTCAtttttt        >  2:1283277/1‑91 (MQ=255)

                                                                                      ttttCAAAACCAATTACTCGTTGCAAGTCATGTGCAACTAGAATGGATAAAGTCATAATTATCCTCCTTCTTCTATCATTTCATTTTttat      <  1:237479/91‑1 (MQ=255)

                                                                                         tCAAAACCAATTACTCGTTGCAAGTCATGTGCAACTAGAATGGATAAAGTCATAATTATCCTCCTTCTTCTATCATTTCATTTTTTATTAc   >  2:2603985/1‑91 (MQ=255)

                                                                                         tCAAAACCAATTACTCGTTGCAAGTCATGTGCAACTAGAATGGATAAAGTCATAATTATCCTCCTTCTTCTATCATTTCATTTTTTATTAc   <  2:46786/91‑1 (MQ=255)

                                                                                          cAAAACCAATTACTCGTTGCAAGTCATGTGCAACTAGAATGGATAAAGTCATAATTATCCTCCTTCTTCTATCATTTCATTTTTTATTACt  >  1:1933084/1‑91 (MQ=255)

ATGTCTTACGACCCATTACTAAAGTATGACCTGTTGATAATTTTTTAACATGCTTCAAATCATTTGGTAGGTGCCAAGGTAATTGATTTTCAAAACCAATTACTCGTTGCAAGTCATGTGCAACTAGAATGGATAAAGTCATAATTATCCTCCTTCTTCTATCATTTCATTTTTTATTACT  >  USA300TCH1516_ALE/1465015‑1465195

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