BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	USA300TCH1516_ALE	2,424,609	C→T	7.9%	intergenic (‑167/+61)	USA300HOU_RS12490 ← / ← USA300HOU_RS12495	hypothetical protein/hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	USA300TCH1516_ALE	2,424,609	0	C	T	7.9%	401.4 / 18.1	126	intergenic (‑167/+61)	USA300HOU_RS12490/USA300HOU_RS12495	hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand): ref base C (47/69); new base T (8/2); total (55/71)
Fisher's exact test for biased strand distribution p-value = 2.06e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.72e-01

Only 100 of 130 total aligned reads displayed.

CTGTCGATGATTGATTCACAATAGCACGTGA‑TTTTTTTGTTTTTCAATATTATTCATTTATTCCATCAAAAACACCCTTTTTAATTTTTACAAAAATTAAAAAAAGTGCTCCTACACTGCTTGCATGTAGAAACACTTTTTCATTGTAATGTTATTCTTCTCGAGACATACCTTTTAGCAT  >  USA300TCH1516_ALE/2424519‑2424699

cTGTCGATGATTGATTCACAATAGCACGTGA‑TTTTTTTGTTTTTCAATATTATTCATTTATTCCATCAAAAACACCCTTTTTAATTTTTAc                                                                                            <  1:3914682/91‑1 (MQ=255)

   tCGATGATTGATTCACAATAGCACGTGATTTTTTTTGTTTTTCAATATTATTCATTTATTCCATCAAAAACACCCTTTTTAATTTTTACaa                                                                                          <  1:4044036/91‑1 (MQ=255)

    cGATGATTGATTCACAATAGCACGTGA‑TTTTTTTGTTTTTCAATATTATTCATTTATTCCATCAAAAACACCCTTTTTAATTTTTACaaaa                                                                                        <  1:2238558/91‑1 (MQ=255)

     gatgatTGATTCACAATAGCACGTGA‑TTTTTTTGTTTTTCAATATTATTCATTTATTCCATCAAAAACACCCTTTTTAATTTTTACaaaaa                                                                                       >  1:2395822/1‑91 (MQ=255)

     gatgatTGATTCACAATAGCACGTGA‑TTTTTTTGTTTTTCAATATTATTCATTTATTCCATCAAAAACACCCTTTTTAATTTTTACaaaaa                                                                                       >  1:1549575/1‑91 (MQ=255)

     gatgatTGATTCACAATAGCACGTGA‑TTTTTTTGTTTTTCAATATTATTCATTTATTCCATCAAAAACACCCTTTTTAATTTTTACaaaaa                                                                                       <  1:2138985/91‑1 (MQ=255)

     gatgatTGATTCACAATAGCACGTGA‑TTTTTTTGTTTTTCAATATTATTCATTTATTCCATCAAAAACACCCTTTTTAATTTTTACaaaaa                                                                                       >  1:590914/1‑91 (MQ=255)

     gatgatTGATTCACAATAGCACGTGA‑TTTTTTTGTTTTTCAATATTATTCATTTATTCCATCAAAAACACCCTTTTTAATTTTTACaaaaa                                                                                       >  1:1549398/1‑91 (MQ=255)

     gatgatTGATTCACAATAGCACGTGA‑TTTTTTTGTTTTTCAATATTATTCATTTATTCCATCAAAAACACCCTTTTTAATTTTTACaaaaa                                                                                       <  1:3088112/91‑1 (MQ=255)

     gatgatTGATTCACAATAGCACGTGA‑TTTTTTTGTTTTTCAATATTATTCATTTATTCCATCAAAAACACCCTTTTTAATTTTTACaaaaa                                                                                       <  1:2690434/91‑1 (MQ=255)

       tgattgatTCACAATAGCACGTGA‑TTTTTTTGTTTTTCAATATTATTCATTTATTCCATCAAAAACACCCTTTTTAATTTTTGTAAAAAtt                                                                                     <  1:1463510/91‑1 (MQ=255)

        gattgattCACAATAGCACGTGA‑TTTTTTTGTTTTTCAATATTATTCATTTATTCCATCAAAAACACCCTTTTTAATTTTTACAAAAATTa                                                                                    <  1:4001426/91‑1 (MQ=255)

        gattgattCACAATAGCACGTGA‑TTTTTTTGTTTTTCAATATTATTCATTTATTCCATCAAAAACACCCTTTTTAATTTTTACAAAAATTa                                                                                    <  1:4001749/91‑1 (MQ=255)

           tgattCACAATAGCACGTGA‑TTTTTTTGTTTTTCAATATTATTCATTTATTCCATCAAAAACACCCTTTTTAATTTTTACAAAAATTaaaa                                                                                 <  1:108757/91‑1 (MQ=255)

            gattCACAATAGCACGTGA‑TTTTTTTGTTTTTCAATATTATTCATTTATTCCATCAAAAACACCCTTTTTAATTTTTACAAAAATTaaaaa                                                                                >  1:3464059/1‑91 (MQ=255)

              ttCACAATAGCACGTGA‑TTTTTTTGTTTTTCAATATTATTCATTTATTCCATCAAAAACACCCTTTTTAATTTTTACAAAAATTaaaaaaa                                                                              >  1:3191486/1‑91 (MQ=255)

               tCACAATAGCACGTGA‑TTTTTTTGTTTTTCAATATTATTCATTTATTCCATCAAAAACACCCTTTTTAATTTTTACAAAAATTAAAAAAAg                                                                             >  1:1231150/1‑91 (MQ=255)

                 acaATAGCACGTGA‑TTTTTTTGTTTTTCAATATTATTCATTTATTCCATCAAAAACACCCTTTTTAATTTTTACAAAAATTAAAAAAAGTg                                                                           <  1:973526/91‑1 (MQ=255)

                   aaTAGCACGTGA‑TTTTTTTGTTTTTCAATATTATTCATTTATTCCATCAAAAACACCCTTTTTAATTTTTACAAAAATTAAGAAAAGTGCt                                                                         >  1:3976707/1‑91 (MQ=255)

                   aaTAGCACGTGA‑TTTTTTTGTTTTTCAATATTATTCATTTATTCCATCAAAAACACCCTTTTTAATTTTTACAAAAATTAAAAAAAGTGcc                                                                         >  1:3445266/1‑90 (MQ=255)

                    aTAGCACGTGA‑TTTTTTTGTTTTTCAATATTATTCATTTATTCCATCAAAAACACCCTTTTTAATTTTTACAAAAATTAAAAAAAGTGCTc                                                                        >  1:3092393/1‑91 (MQ=255)

                     tAGCACGTGA‑TTTTTTTGTTTTTCAATATTATTCATTTATTCCATCAAAAACACCCTTTTTAATTTTTACAAAAATTAAAAAAAGTGCTcc                                                                       >  1:2936542/1‑91 (MQ=255)

                      aGCACGTGA‑TTTTTTTGTTTTTCAATATTATTCATTTATTCCATCAAAAACACCCTTTTTAATTTTTACAAAAATTAAAAAAAGTGCTCCt                                                                      >  1:3413405/1‑91 (MQ=255)

                      aGCACGTGA‑TTTTTTTGTTTTTCAATATTATTCATTTATTCCATCAAAAACACCCTTTTTAATTTTTACAAAAATTAAAAAAAGTGCTCCt                                                                      >  1:3652756/1‑91 (MQ=255)

                      aGCACGTGA‑TTTTTTTGTTTTTCAATATTATTCATTTATTCCATCAAAAACACCCTTTTTAATTTTTACAAAAATTAAAAAAAGTGCTCCt                                                                      <  1:2904561/91‑1 (MQ=255)

                      aGCACGTGA‑TTTTTTTGTTTTTCAATATTATTCATTTATTCCATCAAAAACACCCTTTTTAATTTTTACAAAAATTAAAAAAAGTGCTCCt                                                                      >  1:2192560/1‑91 (MQ=255)

                       gCACGTGA‑TTTTTTTGTTTTTCAATATTATTCATTTATTCCATCAAAAACACCCTTTTTAATTTTTACAAAAATTAAAAAAAGTGCTCCTa                                                                     >  1:1553972/1‑91 (MQ=255)

                       gCACGTGA‑TTTTTTTGTTTTTCAATATTATTCATTTATTCCATCAAAAACACCCTTTTTAATTTTTACAAAAATTAAAAAAAGTGCTCCTa                                                                     <  1:2713123/91‑1 (MQ=255)

                       gCACGTGA‑TTTTTTTGTTTTTCAATATTATTCATTTATTCCATCAAAAACACCCTTTTTAATTTTTACAAAAATTAAAAAAAGTGCTCCTa                                                                     >  1:2636395/1‑91 (MQ=255)

                       gCACGTGA‑TTTTTTTGTTTTTCAATATTATTCATTTATTCCATCAAAAACACCCTTTTTAATTTTTACAAAAATTAAAAAAAGTGCTCCTa                                                                     >  1:871485/1‑91 (MQ=255)

                        cACGTGA‑TTTTTTTGTTTTTCAATATTATTCATTTATTCCATCAAAAACACCCTTTTTAATTTTTACAAAAATTAAAAAAAGTGCTCCTac                                                                    <  1:672573/91‑1 (MQ=255)

                        cACGTGA‑TTTTTTTGTTTTTCAATATTATTCATTTATTCCATCAAAAACACCCTTTTTAATTTTTACAAAAATTAAAAAAAGTGCTCCTac                                                                    <  1:2990075/91‑1 (MQ=255)

                          cGTGA‑TTTTTTTGTTTTTCAATATTATTCATTTATTCCATCAAAAACACCCTTTTTAATTTTTACAAAAATTAAAAAAAGTGCTCCTacac                                                                  <  1:2387948/91‑1 (MQ=255)

                           gTGA‑TTTTTTTGTTTTTCAATATTATTCATTTATTCCATCAAAAACACCCTTTTTAATTTTTACAAAAATTAAAAAAAGTGCTCCTACACt                                                                 <  1:1633502/91‑1 (MQ=255)

                           gTGA‑TTTTTTTGTTTTTCAATATTATTCATTTATTCCATCAAAAACACCCTTTTTAATTTTTACAAAAATTAAAAAAAGTGCTCCTACACt                                                                 <  1:3986458/91‑1 (MQ=255)

                            tGA‑TTTTTTTGTTTTTCAATATTATTCATTTATTCCATCAAAAACACCCTTTTTAATTTTTACAAAAATTAAAAAAAGTGCTCCTACACTg                                                                <  1:1652271/91‑1 (MQ=255)

                             gA‑TTTTTTTGTTTTTCAATATTATTCATTTATTCCATCAAAAACACCCTTTTTAATTTTTACAAAAATTAAAAAAAGTGCTCCTACACTGc                                                               <  1:3797955/91‑1 (MQ=255)

                                tttttttGTTTTTCAATATTATTCATTTATTCCATCAAAAACACCCTTTTTAATTTTTACAAAAATTAAAAAAAGTGCTCCTACACTGCtt                                                             <  1:1328490/91‑1 (MQ=255)

                                  tttttGTTTTTCAATATTATTCATTTATTCCATCAAAAACACCCTTTTTAATTTTTACAAAAATTAAAAAAAGTGCTCCTACACTGCTTGc                                                           <  1:4115001/91‑1 (MQ=255)

                                  tttttGTTTTTCAATATTATTCATTTATTCCATCAAAAACACCCTTTTTAATTTTTACAAAAATTAAAAAAAGTGCTCCTACACTGCTTGc                                                           >  1:86050/1‑91 (MQ=255)

                                     ttGTTTTTCAATATTATTCATTTATTCCATCAAAAACACCCTTTTTAATTTTTACAAAAATTAAAAAAAGTGCTCCTACACTGCTTGCATg                                                        <  1:3951820/91‑1 (MQ=255)

                                       gTTTTTCAATATTATTCATTTATTCCATCAAAAACACCCTTTTTAATTTTTACAAAAATTAAAAAAAGTGCTCCTACACTGCTTGCATGTa                                                      <  1:3329492/91‑1 (MQ=255)

                                       gTTTTTCAATATTATTCATTTATTCCATCAAAAACACCCTTTTTAATTTTTACAAAAATTAAAAAAAGTGCTCCTACACTGCTTGCATGTa                                                      >  1:1087318/1‑91 (MQ=255)

                                         ttttCAATATTATTCATTTATTCCATCAAAAACACCCTTTTTAATTTTTACAAAAATTAAAAAAAGTGCTCCTACACTGCTTGCATGTAGa                                                    <  1:3783007/91‑1 (MQ=255)

                                          tttCAATATTATTCATTTATTCCATCAAAAACACCCTTTTTAATTTTTACAAAAATTAAAAAAAGTGCTCCTACACTGCTTGCATGTAGaa                                                   >  1:827536/1‑91 (MQ=255)

                                             cAATATTATTCATTTATTCCATCAAAAACACCCTTTTTAATTTTTACAAAAATTAAAAAAAGTGCTCCTACACTGCTTGCATGTAGAAaca                                                >  1:2808962/1‑91 (MQ=255)

                                              aaTATTATTCATTTATTCCATCAAAAACACCCTTTTTAATTTTTACAAAAATTAAAAAAAGTGCTCCTACACTGCTTGCATGTAGAAacac                                               <  1:3783298/91‑1 (MQ=255)

                                              aaTATTATTCATTTATTCCATCAAAAACACCCTTTTTAATTTTTACAAAAATTAAAAAAAGTGCTCCTACACTGCTTGCATGTAGAAacac                                               <  1:4065063/91‑1 (MQ=255)

                                              aaTATTATTCATTTATTCCATCAAAAACACCCTTTTTAATTTTTACAAAAATTAAAAAAAGTGCTCCTACACTGCTTGCATGTAGAAacac                                               <  1:898711/91‑1 (MQ=255)

                                              aaTATTATTCATTTATTCCATCAAAAACACCCTTTTTAATTTTTACAAAAATTAAAAAAAGTGCTCCTACACTGCTTGCATGTAGAAacac                                               <  1:2572191/91‑1 (MQ=255)

                                                tattatTGATTTATTCCATCAAAAACACCCTTTTTAATTTTTACAAAAATTAAAAAAAGTGCTCCTACACTGCTTGCATGTAGAAACACtt                                             <  1:2075055/91‑1 (MQ=255)

                                                tattatTCATTTATTCCATCAAAAACACCCTTTTTAATTTTTACAAAAATTAAAAAAAGTGCTCCTACACTGCTTGCATGTAGAAACACtt                                             <  1:1759315/91‑1 (MQ=255)

                                                 attattCATTTATTCCATCAAAAACACCCTTTTTAATTTTTACAAAAATTAAAAAAAGTGCTCCTACACTGCTTGCATGTAGAAACACttt                                            >  1:3295213/1‑91 (MQ=255)

                                                 attattCATTTATTCCATCAAAAACACCCTTTTTAATTTTTACAAAAATTAAAAAAAGTGCTCCTACACTGCTTGCATGTAGAAACACttt                                            >  1:2688273/1‑91 (MQ=255)

                                                   tattCATTTATTCCATCAAAAACACCCTTTTTAATTTTTACAAAAATTAAAAAAAGTGCTCCTACACTGCTTGCATGTAGAAACACttttt                                          >  1:1032420/1‑91 (MQ=255)

                                                     ttCATTTATTCCATCAAAAACACCCTTTTTAATTTTTACAAAAATTAAAAAAAGTGCTCCTACGCTGCTTGCATGTAGAAACACTTTTTCa                                        <  1:2491721/91‑1 (MQ=255)

                                                         tttATTCCATCAAAAACACCCTTTTTAATTTTTACAAAAATTAAAAAAAGTGCTCCTACACTGCTTGCATGTAGAAACACTTTTTCATTGt                                    <  1:2389476/91‑1 (MQ=255)

                                                         tttATTCCATCAAAAACACCCTTTTTAATTTTTACAAAAATTAAAAAAAGTGCTCCTACACTGCTTGCATGTAGAAACACTTTTTCATTGt                                    <  1:2237752/91‑1 (MQ=255)

                                                          ttATTCCATCAAAAACACCCTTTTTAATTTTTACAAAAATTAAAAAAAGTGCTCCTACACTGCTTGCATGTAGAAACACTTTTTCATTGTa                                   <  1:1280237/91‑1 (MQ=255)

                                                             ttCCATCAAAAACACCCTTTTTAATTTTTACAAAAATTAAAAAAAGTGCTCCTACACTGCTTGCATGTAGAAACACTTTTTCATTGTAATg                                >  1:2357404/1‑91 (MQ=255)

                                                             ttCCATCAAAAACACCCTTTTTAATTTTTACAAAAATTAAAAAAAGTGCTCCTACACTGCTTGCATGTAGAAACACTTTTTCATTGTAATg                                <  1:2177416/91‑1 (MQ=255)

                                                              tCCATCAAAAACACCCTTTTTAATTTTTACAAAAATTAAAAAAAGTGCTCCTACACTGCTTGCATGTAGAAACACTTTTTCATTGTAATGt                               <  1:1876629/91‑1 (MQ=255)

                                                              tCCATCAAAAACACCCTTTTTAATTTTTACAAAAATTAAAAAAAGTGCTCCTACACTGCTTGCATGTAGAAACACTTTTTCATTGTAATGt                               <  1:1876628/91‑1 (MQ=255)

                                                               ccATCAAAAACACCCTTTTTAATTTTTACAAAAATTAAAAAAAGTGCTCCTACACTGCTTGCATGTAGAAACACTTTTTCATTGTAATGtt                              <  1:1603643/91‑1 (MQ=255)

                                                               ccATCAAAAACACCCTTTTTAATTTTTACAAAAATTAAAAAAAGTGCTCCTACACTGCTTGCATGTAGAAACACTTTTTCATTGTAATGtt                              <  1:2515874/91‑1 (MQ=255)

                                                               ccATCAAAAACACCCTTTTTAATTTTTACAAAAATTAAAAAAAGTGCTCCTACACTGCTTGCATGTAGAAACACTTTTTCATTGTAATGtt                              >  1:3954125/1‑91 (MQ=255)

                                                               ccATCAAAAACACCCTTTTTAATTTTTACAAAAATTAAAAAAAGTGCTCCTACACTGCTTGCATGTAGAAACACTTTTTCATTGTAATGtt                              <  1:2751430/91‑1 (MQ=255)

                                                                 aTCAAAAACACCCTTTTTAATTTTTACAAAAATTAAAAAAAGTGCTCCTACACTGCTTGCATGTAGAAACACTTTTTCATTGTAATGTTAt                            <  1:1477738/91‑1 (MQ=255)

                                                                    aaaaaCACCCTTTTTAATTTTTACAAAAATTAAAAAAAGTGCTCCTACACTGCTTGCATGTAGAAACACTTTTTCATTGTAATGTTAttct                         <  1:1049505/91‑1 (MQ=255)

                                                                    aaaaaCACCCTTTTTAATTTTTACAAAAATTAAAAAAAGTGCTCCTACACTGCTTGCATGTAGAAACACTTTTTCATTGTAATGTTAttct                         <  1:1481287/91‑1 (MQ=255)

                                                                      aaaCACCCTTTTTAATTTTTACAAAAATTAAAAAAAGTGCTCCTACACTGCTTGCATGTAGAAACACTTTTTCATTGTAATGTTAttcttc                       >  1:1179163/1‑91 (MQ=255)

                                                                      aaaCACCCTTTTTAATTTTTACAAAAATTAAAAAAAGTGCTCCTACACTGCTTGCATGTAGAAACACTTTTTCATTGTAATGTTAttcttc                       >  1:1178418/1‑91 (MQ=255)

                                                                      aaaCACCCTTTTTAATTTTTACAAAAATTAAAAAAAGTGCTCCTACACTGCTTGCATGTAGAAACACTTTTTCATTGTAATGTTAttcttc                       >  1:2064559/1‑91 (MQ=255)

                                                                           cccTTTTTAATTTTTACAAAAATTAAAAAAAGTGCTCCTACACTGCTTGCATGTAGAAACACTTTTTCATTGTAATGTTATTCTTCTCgag                  <  1:3614894/91‑1 (MQ=255)

                                                                             ttttttGATTTTTGTAAAAATTAAAAAGGGTGCTCCTACACTGCTTGCATGTAGAAACACTTTTTCATTGTAATGTTATTCTTCTCGAGAc                >  1:3232733/2‑91 (MQ=37)

                                                                             ttttttAATTTTTGTAAAAATTAAAAAGGGTGCTCCTACACTGCTTGCATGTAGAAACACTTTTTCATTGTAATGTTATTCTTCTCGAGAc                >  1:3232461/2‑91 (MQ=37)

                                                                             cTTTTTAATTTTTACAAAAATTAAAAAAAGTGCTCCTACACTGCTTGCATGTAGAAACACTTTTTCATTGTAATGTTATTCTTCTCGAGAc                <  1:4052199/91‑1 (MQ=255)

                                                                             cTTTTTAATTTTTACAAAAATTAAAAAAAGTGCTCCTACACTGCTTGCATGTAGAAACACTTTTTCATTGTAATGTTATTCTTCTCGAGAc                <  1:624395/91‑1 (MQ=255)

                                                                             cTTTTTAATTTTTACAAAAATTAAAAAAAGTGCTCCTACACTGCTTGCATGTAGAAACACTTTTTCATTGTAATGTTATTCTTCTCGAGAc                <  1:4052119/91‑1 (MQ=255)

                                                                              tttttAATTTTTACAAAAATTAAAAAAAGTGCTCCTACACTGCTTGCATGTAGAAACACTTTTTCATTGTAATGTTATTCTTCTCGAGACa               >  1:3660217/1‑91 (MQ=255)

                                                                              tttttAATTTTTACAAAAATTAAAAAAAGTGCTCCTACACTGCTTGCATGTAGAAACACTTTTTCATTGTAATGTTATTCTTCTCGAGACa               <  1:356783/91‑1 (MQ=255)

                                                                               ttttAATTTTTGTAAAAATTAAAAAGGGTGCTCCTACACTGCTTGCATGTAGAAACACTTTTTCATTGTAATGTTATTCTTCTCGAGACAt              >  1:1208391/1‑91 (MQ=37)

                                                                               ttttAATTTTTACAAAAATTAAAAAAAGTGCTCCTACACTGCTTGCATGTAGAAACACTTTTTCATTGTAATGTTATTCTTCTCGAGACAt              >  1:578983/1‑91 (MQ=255)

                                                                               ttttAATTTTTACAAAAATTAAAAAAAGTGCTCCTACACTGCTTGCATGTAGAAACACTTTTTCATTGTAATGTTATTCTTCTCGAGACAt              >  1:579465/1‑91 (MQ=255)

                                                                                tttAATTTTTGTAAAAATTAAAAAGGGTGCTCCTACACTGCTTGCATGTAGAAACACTTTTTCATTGTAATGTTATTCTTCTCGAGACATa             >  1:209106/1‑91 (MQ=37)

                                                                                  tAATTTTTACAAAAATTAAAAAAAGTGCTCCTACACTGCTTGCATGTAGAAACACTTTTTCATTGTAATGTTATTCTTCTCGAGACATAcc           >  1:826489/1‑91 (MQ=255)

                                                                                   aaTTTTTGTAAAAATTAAAAAGGGTGCTCCTACACTGCTTGCATGTAGAAACACTTTTTCATTGTAATGTTATTCTTCTCGAGACATACCt          >  1:2838347/1‑91 (MQ=38)

                                                                                   aaTTTTTGTAAAAATTAAAAAAAGTGCTCCTACACTGCTTGCATGTAGAAACACTTTTTCATTGTAATGTTATTCTTCTCGAGACATACCt          >  1:441544/1‑91 (MQ=255)

                                                                                   aaTTTTTACAAAAATTAAAAAAAGTGCTCCTACACTGCTTGCATGTAGAAACACTTTTTCATTGTAATGTTATTCTTCTCGAGACATACCt          >  1:1180389/1‑91 (MQ=255)

                                                                                   aaTTTTTACAAAAATTAAAAAAAGTGCTCCTACACTGCTTGCATGTAGAAACACTTTTTCATTGTAATGTTATTCTTCTCGAGACATACCt          <  1:1791914/91‑1 (MQ=255)

                                                                                     tttttACAAAAATTAAAAAAAGTGCTCCTACACTGCTTGCATGTAGAAACACTTTTTCATTGTAATGTTATTCTTCTCGAGACATACCttt        >  1:3474615/1‑91 (MQ=255)

                                                                                     tttttACAAAAATTAAAAAAAGTGCTCCTACACTGCTTGCATGTAGAAACACTTTTTCATTGTAATGTTATTCTTCTCGAGACATACCttt        >  1:702982/1‑91 (MQ=255)

                                                                                     tttttACAAAAATTAAAAAAAGTGCTCCTACACTGCTTGCATGTAGAAACACTTTTTCATTGTAATGTTATTCTTCTCGAGACATACCttt        <  1:3877619/91‑1 (MQ=255)

                                                                                     tttttACAAAAATTAAAAAAAGTGCTCCTACACTGCTTGCATGTAGAAACACTTTTTCATTGTAATGTTATTCTTCTCGAGACATACCttt        <  1:3878155/91‑1 (MQ=255)

                                                                                      ttttACAAAAATTAAAAAAAGTGCTCCTACACTGCTTGCATGTAGAAACACTTTTTCATTGTAATGTTATTCTTCTCGAGACATACCtttt       <  1:3912297/91‑1 (MQ=255)

                                                                                        ttACAAAAATTAAAAAAAGTGCTCCTACACTGCTTGCATGTAGAAACACTTTTTCATTGTAATGTTATTCTTCTCGAGACATACCTTTTAg     <  1:488656/91‑1 (MQ=255)

                                                                                         tACAAAAATTAAAAAAAGTGCTCCTACACTGCTTGCATGTAGAAACACTTTTTCATTGTAATGTTATTCTTCTCGAGACATACCTTTTAGc    <  1:693020/91‑1 (MQ=255)

                                                                                          aCAAAAATTAAAAAAAGTGCTCCTACACTGCTTGCATGTAGAAACACTTTTTCATTGTAATGTTATTCTTCTCGAGACATACCTTTTAGCa   <  1:2906540/91‑1 (MQ=255)

                                                                                          aCAAAAATTAAAAAAAGTGCTCCTACACTGCTTGCATGTAGAAACACTTTTTCATTGTAATGTTATTCTTCTCGAGACATACCTTTTAGCa   <  1:2371971/91‑1 (MQ=255)

                                                                                           cAAAAATTAAAAAAAGTGCTCCTACACTGCTTGCATGTAGAAACACTTTTTCATTGTAATGTTATTCTTCTCGAGACATACCTTTTAGCat  <  1:3056074/91‑1 (MQ=255)

CTGTCGATGATTGATTCACAATAGCACGTGA‑TTTTTTTGTTTTTCAATATTATTCATTTATTCCATCAAAAACACCCTTTTTAATTTTTACAAAAATTAAAAAAAGTGCTCCTACACTGCTTGCATGTAGAAACACTTTTTCATTGTAATGTTATTCTTCTCGAGACATACCTTTTAGCAT  >  USA300TCH1516_ALE/2424519‑2424699

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