BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	USA300TCH1516_ALE	1,465,110	C→T	7.4%	G16D (GGT→GAT)	folA ←	Dihydrofolate reductase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	USA300TCH1516_ALE	1,465,110	0	C	T	7.4%	733.4 / 31.5	229	G16D (GGT→GAT)	folA	Dihydrofolate reductase
Reads supporting (aligned to +/- strand): ref base C (103/109); new base T (13/4); total (116/113)
Fisher's exact test for biased strand distribution p-value = 4.12e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.97e-01

Only 100 of 237 total aligned reads displayed.

TTACGACCCATTACTAAAGTATGACCTGTTGATAATTTTTTAACATGCTTCAAATCATTTGGTAGGTGCCAAGGTAATTGATTTTCAAAACCAATTACTCGTTGCAAGTCATGTGCAACTAGAATGGATAAAGTCATAATTATCCTCCTTCTTCTATCATTTCATTTTTTATTACTAAGTT  >  USA300TCH1516_ALE/1465020‑1465200

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

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

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

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

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

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

       ccATTACTAAAGTATGACCTGTTGATAATTTTTTAACATGCTTCAAATCATTTGGTAGGTGCCAAGGTAATTGATTTTCAAAACCAATTAc                                                                                     <  1:1551587/91‑1 (MQ=255)

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

         aTTACTAAAGTATGACCTGTTGATAATTTTTTAACATGCTTCAAATCATTTGGTAGGTGCCAAGGTAATTGATTTTCAAAACCAATTACTc                                                                                   >  1:1152390/1‑91 (MQ=255)

         aTTACTAAAGTATGACCTGTTGATAATTTTTTAACATGCTTCAAATCATTTGGTAGGTGCCAAGGTAATTGATTTTCAAAACCAATTACTc                                                                                   <  2:1113257/91‑1 (MQ=255)

           tACTAAAGTATGACCTGTTGATAATTTTTTAACATGCTTCAAATCATTTGGTAGGTGCCAAGGTAATTGATTTTCAAAACCAATTACTCGt                                                                                 <  1:2471375/91‑1 (MQ=255)

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

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

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

              tAAAGTATGACCTGTTGATAATTTTTTAACATGCTTCAAATCATTTGGTAGGTGCCAAGGTAATTGATTTTCAAAACCAATTACTCGTTgt                                                                              <  2:1729868/91‑2 (MQ=255)

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

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

                aaGTATGACCTGTTGATAATTTTTTAACATGCTTCAAATCATTTGGTAGGTGCCAAGGTAATTGATTTTCAAAACCAATTACTCGTTGCaa                                                                            >  2:3385115/1‑91 (MQ=255)

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

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

                  gTATGACCTGTTGATAATTTTTTATCATGCTTCAAATCATTTGGTAGGTGCCAAGGTAATTGATTTTCAAAACCAATTACTCGTTGCTAGt                                                                          >  2:1442582/1‑91 (MQ=255)

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

                  gTATGACCTGTTGATAATTTTTTAACATGCTTCAAATCATTTGGTAGGTGCCAAGGTAATTGATTTTCAAAACCAATTACTCGTTGCAAGt                                                                          >  1:581768/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)

                      gACCTGTTGATAATTTTTTAACATGCTTCAAATCATTTGGTAGGGGCCAAGGTAATTGATTTTCAAAACCACTTACTCGTTGCAAGTCAtg                                                                      <  1:2749102/91‑1 (MQ=255)

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

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

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

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

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

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

                              gATAATTTTTTAACATGCTTCAAATCATTTGGTAGGTGCCAAGGTAATTGATTTTCAAAACCAATTACTCGTTGCAAGTCATGTGCAACTa                                                              >  1:3862828/1‑91 (MQ=255)

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

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

                                tAATTTTTTAACATGCTTCAAATCATTTGTTAGGTGCCAAGGTAATTGATTTTCAAAACCAATTACTCGTTGCAAGTCATGTGCAACTAGa                                                            >  2:1392099/1‑91 (MQ=255)

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

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

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

                                    tttttAACATGCTTCAAATCATTTGGTAGGTGCCAAGGTAATTGATTTTCAAAACCAATTACTCGTTGCAAGTCATGTGCAACTAGAATgg                                                        <  1:2089809/91‑1 (MQ=255)

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

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

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

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

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

                                         aaCATGCTTCAAATCATTTGGTAGGTGCCAAGGTAATTGATTTTCAAAACCAATTACTCGTTGCAAGTCATGTGCAACTAGAATGGATaaa                                                   <  2:1041961/91‑1 (MQ=255)

                                           cATGCTTCAAATCATTTGGTAGGTGCCAAGGTAATTGATTTTCAAAACCAATTACTCGTTGCAAGTCATGTGCAACTAGAATGGATAAAGt                                                 >  2:120408/1‑91 (MQ=255)

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

                                              gCTTCAAATCATTTGGTAGGTGCCAAGGTAATTGATTTTCAAAACCAATTACTCGTTGCAAGTCATGTGCAACTAGAATGGATAAAGTCAt                                              >  1:2657098/1‑91 (MQ=255)

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

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

                                                     aTCATTTGGTAGTTGCCAAGGTAATTGATTTTCAAAACCAATTACTCGTTGCAAGTCATGTGCAACTAGAATGGATAAAGTCATAATTAtc                                       <  2:3224640/91‑1 (MQ=255)

                                                     aTCATTTGGTAGGTGCCAAGGTAATTGATTTTCAAAACCAATTACTCGTTGCAAGTCATGTGCAACTAGAATGGATAAAGTCATAATTAtc                                       >  1:3162602/1‑91 (MQ=255)

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

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

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

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

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

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

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

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

                                                          ttGGTAGGTGCCAAGGTAATTGATTTTCAAAACCAATTACTCGTTGCAAGTCATGTGCAACTAGAATGGATAAAGTCATAATTATcctcct                                  <  1:1807008/91‑1 (MQ=255)

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

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

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

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

                                                                   gCCAAGGTAATTGATTTTCAAAACCAATTACTCGTTGCAAGTCATGTGCAACTAGAATGGATAAAGTCATAATTATCCTCCTTCTTCTAtc                         <  2:3202468/91‑1 (MQ=255)

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

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

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

                                                                     cAAGGTAATTGATTTTCAAAACCAATTACTCGTTGCAAGTCATGTGCAACTAGAATGGATAAAGTCATAATTATCCTCCTTCTTCTAtcat                       >  2:3089630/1‑91 (MQ=255)

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

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

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

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

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

                                                                          tAATTGATTTTCAAAACCAATTACTCGTTGCAAGTCATGTGCAACTAGAATGGATAAAGTCATAATTATCCTCCTTCTTATAtcatttcat                  >  1:3519613/1‑91 (MQ=255)

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

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

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

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

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

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

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

                                                                             ttGATTTTCAAAACCAATTACTCGTTGCAAGTCAt                                                                       >  1:1886782/1‑35 (MQ=255)

                                                                             ttGATTTTCAAAACCAATTACTCGTTGCAAGTCATGTGCAACTAGAATGGATAAAGTCATAATTATCCTCCTTCTTCTATCATTTCAtttt               >  2:3246362/1‑91 (MQ=255)

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

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

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

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

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

                                                                                        aaCCAATTACTCGTTGCAAGTCATGTGCAACTAGAATGGATAAAGTCATAATTATCCTCCTTCTTCTATCATTTCATTTTTTATTACTAAg    >  2:2715085/1‑91 (MQ=255)

                                                                                         aCCAATTACTCGTTGCAAGTCATGTGCAACTAGAATGGATAAAGTCATAATTATCCTCCTTCTTCTATCATTTCATTTTTTATTACTAAGt   >  1:960789/1‑91 (MQ=255)

                                                                                         aCCAATTACTCGTTGCAAGTCATGTGCAACTAGAATGGATAAAGTCATAATTATCCTCCTTCTTCTATCATTTCATTTTTTATTACTAAGt   >  1:3114758/1‑91 (MQ=255)

                                                                                         aCCAATTACTCGTTGCAAGTCATGTGCAACTAGAATGGATAAAGTCATAATTATCCTCCTTCTTCTATCATTTCATTTTTTATTACTAAGt   <  1:2674086/91‑1 (MQ=255)

                                                                                         aCCAATTACTCGTTGCAAGTCATGTGCAACTAGAATGGATAAAGTCATAATTATCCTCCTTCTTCTATCATTTCATTTTTTATTACTAAGt   >  1:1029353/1‑91 (MQ=255)

                                                                                          ccAATTACTCGTTGCAAGTCATGTGCAACTAGAATGGATAAAGTCATAATTATCCTCCTTCTTCTATCATTTCATTTTTTATTACTAAGtt  >  2:1388097/1‑91 (MQ=255)

                                                                                          ccAATTACTCGTTGCAAGTCATGTGCAACTAGAATGGATAAAGTCATAATTATCCTCCTTCTTCTATCATTTCATTTTTTATTACTAAGtt  <  1:2958900/91‑1 (MQ=255)

                                                                                          ccAATTACTCGTTGCAAGTCATGTGCAACTAGAATGGATAAAGTCATAATTATCCTCCTTCTTCTATCATTTCATTTTTTATTACTAAGtt  >  1:1226072/1‑91 (MQ=255)

TTACGACCCATTACTAAAGTATGACCTGTTGATAATTTTTTAACATGCTTCAAATCATTTGGTAGGTGCCAAGGTAATTGATTTTCAAAACCAATTACTCGTTGCAAGTCATGTGCAACTAGAATGGATAAAGTCATAATTATCCTCCTTCTTCTATCATTTCATTTTTTATTACTAAGTT  >  USA300TCH1516_ALE/1465020‑1465200

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