BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	USA300TCH1516_ALE	8,213	T→A	5.7%	I394I (ATT→ATA)	gyrA →	DNA gyrase subunit A

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	USA300TCH1516_ALE	8,213	0	T	A	5.7%	652.2 / 19.4	191	I394I (ATT→ATA)	gyrA	DNA gyrase subunit A
Reads supporting (aligned to +/- strand): ref base T (75/105); new base A (2/9); total (77/114)
Fisher's exact test for biased strand distribution p-value = 2.04e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.38e-01

Only 100 of 205 total aligned reads displayed.

AAGACGTACGCAATACAACTTACGTAAAGCTAAAGATCGTGCCCACATTTTAGAAGGATTACGTATCGCACTTGACCATATCGATGAAATTATTTCAACGATTCGTGAGTCAGATACAGATAAAGTTGCAATGGAAAGCTTGCAACAACGCTTCAAACTTTCTGAAAAACAAGCTCAAGCT  >  USA300TCH1516_ALE/8123‑8303

aagaCGTACGCAATACAACTTACGTAAAGCTAAAGATCGTGCCCACATTTTAGAAGGATTACGTATCGCACTTGACCATATCGATGAAatt                                                                                            >  1:2670414/1‑91 (MQ=255)

  gaCGTACGCAATACAACTTACGTAAAGCTAAAGATCGTGCCCACATTTTAGAAGGATTACGTATCGCACTTGACCATATCGATGAAattat                                                                                          >  1:3393399/1‑91 (MQ=255)

  gaCGTACGCAATACAACTTACGTAAAGCTAAAGATCGTGCCCACATTTTAGAAGGATTACGTATCGCACTTGACCATATCGATGAAattat                                                                                          >  2:2899575/1‑91 (MQ=255)

  gaCGTACGCAATACAACTTACGTAAAGCTAAAGATCGTGCCCACATTTTAGAAGGATTACGTATCGCACTTGACCATATCGATGAAattat                                                                                          >  1:1687515/1‑91 (MQ=255)

    cGTACGCAATACAACTTACGTAAAGCTAAAGATCGTGCCCACATTTTAGAAGGATTACGTATCGCACTTGACCATATCGATGAAATTAttt                                                                                        >  1:932551/1‑91 (MQ=255)

     gTACGCAATACAACTTACGTAAAGCTAAAGATCGTGCCCACATTTTAGAAGGATTACGTATCGCACTTGACCATATCGTTGAAATAATTTc                                                                                       <  2:3717009/91‑1 (MQ=255)

     gTACGCAATACAACTTACGTAAAGCTAAAGATCGTGCCCACATTTTAGAAGGATTACGTATCGCACTTGACCATATCGATGAAATTATTTc                                                                                       <  2:879559/91‑1 (MQ=255)

      tACGCAATACAACTTACGTAAAGCTAAAGCTCGTGCCCACATTTTAGAAGGATTACGTATCGCACTTGACCATATCGATGAAATTATTTCa                                                                                      <  1:3502157/91‑1 (MQ=255)

      tACGCAATACAACTTACGTAAAGCTAAAGATCGTGCCCACATTTTAGAAGGATTACGTATCGCACTTGACCATATCGTTGAAATAATTTCa                                                                                      <  1:1771005/91‑1 (MQ=255)

      tACGCAATACAACTTACGTAAAGCTAAAGATCGTGCCCACATTTTAGAAGGATTACGTATCGCACTTGACCATATCGATGAAATTATTTCa                                                                                      <  2:3037110/91‑1 (MQ=255)

      tACGCAATACAACTTACGTAAAGCTAAAGATCGTGCCCACATTTTAGAAGGATTACGTATCGCACTTGACCATATCGATGAAATTATTTCa                                                                                      <  2:271579/91‑1 (MQ=255)

      tACGCAATACAACTTACGTAAAGCTAAAGATCGTGCCCACATTTTAGAAGGATTACGTATCGCACTTGACCATATCGATGAAATTATTTCa                                                                                      <  1:52344/91‑1 (MQ=255)

       aCGCAATACAACTTACGTAAAGCTAAAGATCGTGCCCACATTTTAGAAGGATTACGTATCGCACTTGACCATATCGATGAAATTATTTCaa                                                                                     >  2:2096163/1‑91 (MQ=255)

       aCGCAATACAACTTACGTAAAGCTAAAGATCGTGCCCACATGGTAGAAAGATTACGTATCGAACTTGACCATATCGTTGAAATAATTTCat                                                                                     <  2:1237205/91‑2 (MQ=255)

        cGCAATACAACTTACGTAAAGCTCAAGCTCGTGCCCACATTTTAGCAGGATTACGTATCGCAACTGACCATATCGCTGAAATTATTTCAAc                                                                                    >  2:2469933/1‑91 (MQ=255)

        cGCAATACAACTTACGTAAAGCTAAAGATCGTGCCCACATTTTAGAAGGATTACGTATCGCACTTGACCATATCGATGAAATTATTTCAAc                                                                                    <  1:3634085/91‑1 (MQ=255)

        cGCAATACAACTTACGTAAAGCTAAAGATCGTGCCCACATTTTAGAAGGATTACGTATCGCACTTGACCATATCGATGAAATTATTTCAAc                                                                                    >  2:3724051/1‑91 (MQ=255)

          cAATACAACTTACGTAAAGCTAAAGATCGTGCCCACATTTTAGAAGGATTACGTATCGCACTTGACCATATCGATGAAATTATTTCAACGa                                                                                  >  2:1198435/1‑91 (MQ=255)

           aaTACAACTTACGTAAAGCTAAAGATCGTGCCCACATTTTAGAAGGATTACGTATCGCACTTGACCATATCGATGAAATTATTTCAACGAt                                                                                 <  2:2508448/91‑1 (MQ=255)

           aaTACAACTTACGTAAAGCTAAAGATCGTGCCCACATTTTAGAAGGATTACGTATCGCACTTGACCATATCGATGAAATTATTTCAACGAt                                                                                 <  1:1359340/91‑1 (MQ=255)

            aTACAACTTACGTAAAGCTAAAGATCGTGCCCACATTTTAGAAGGATTACGTATCGCACTTGACCATATCGTTGAAATAATTTCATCGAta                                                                                <  2:33572/91‑2 (MQ=39)

            aTACAACTTACGTAAAGCTAAAGATCGTGCCCACATTTTAGAAGGATTACGTATCGCACTTGACCATATCGTTGAAATAATTTCATCGAta                                                                                <  2:32811/91‑2 (MQ=39)

                aaCTTACGTAAAGCTAAAGATCGTGCCCACATTTTAGAAGGATTACGTATCGCACTTGACCATATCGATGAAATTATTTCAACGATTCGTg                                                                            >  2:3444111/1‑91 (MQ=255)

                aaCTTACGTAAAGCTAAAGATCGTGCCCACATTTTAGAAGGATTACGTATCGCACTTGACCATATCGATAAAATTATTTCAACGATTCGTg                                                                            >  1:1755111/1‑91 (MQ=255)

                 aCTTACGTAAAGCTAAAGATCGTGCCCACATTTTAGAAGGATTACGTATCGCACTTGACCATATCGATGAAATTATTTCAACGATTCGTGa                                                                           <  1:1010931/91‑1 (MQ=255)

                 aCTTACGTAAAGCTAAAGATCGTGCCCACATTTTAGAAGGATTACGTATCGCACTTGACCATATCGATGAAATTATTTCAACGATTCGTGa                                                                           <  2:1666528/91‑1 (MQ=255)

                   ttACGTAAAGCTAAAGATCGTGCCCACATTTTAGAAGGATTACGTATCGCACTTGACCATATCGATGAAATTATTTCAACGATTCGTGAGt                                                                         <  2:721837/91‑1 (MQ=255)

                    tACGTAAAGCTAAAGATCGTGCCCACATTTTAGAAGGATTACGTATCGCACTTGACCATATCGATGAAATTATTTCAACGATTCGTGAGTc                                                                        >  2:1047710/1‑91 (MQ=255)

                    tACGTAAAGCTAAAGATCGTGCCCACATTTTAGAAGGATTACGTATCGCACTTGACCATATCGATGAAATTATTTCAACGATTCGTGAGTc                                                                        <  1:1403813/91‑1 (MQ=255)

                    tACGTAAAGCTAAAGATCGTGCCCACATTTTAGAAGGATTACGTATCGCACTTGACCATATCGATGAAATTATTTCAACGATTCGTGAGTc                                                                        >  1:857414/1‑91 (MQ=255)

                         aaaGCTAAAGATCGTGCCCACATTTTAGAAGGATTACGTATCGCACTTGACCATATCGATGAAATTATTTCAACGATTCGTGAGTcagata                                                                   >  1:221605/1‑91 (MQ=255)

                         aaaGCTAAAGATCGTGCCCACATTTTAGAAGGATTACGTATCGCACTTGACCATATCGATGAAATTATTTCAACGATTCGTGAGTcagata                                                                   <  1:3428784/91‑1 (MQ=255)

                          aaGCTAAAGATCGTGCCCACATTTTAGAAGGATTACGTATCGCACTTGACCATATCGATGAAATTATTTCAACGATTCGTGAGTcagatac                                                                  <  1:2899575/91‑1 (MQ=255)

                            gCTAAAGATCGTGCCCACATTTTAGAAGGATTACGTATCGCACTTGACCATATCTATGAAATTATTTCAACGATTCGTGAGTcagatacag                                                                <  1:2104167/91‑1 (MQ=255)

                            gCTAAAGATCGTGCCCACATTTTAGAAGGATTACGTATCGCACTTGACCATATCGATGAAATTATTTCAACGATTCGTGAGTcagatacag                                                                <  2:2885449/91‑1 (MQ=255)

                             cTAAAGATCGTGCCCACATTTTAGAAGGATTACGTATCGCACTTGACCATATCGATGAAATTATTTCAACGATTCGTGAGTcagatacaga                                                               <  2:2406668/91‑1 (MQ=255)

                             cTAAAGATCGTGCCCACATTTTAGAAGGATTACGTATCGCACTTGACCATATCGATGAAATTATTTCAACGATTCGTGAGTcagatacaga                                                               <  2:2340433/91‑1 (MQ=255)

                              tAAAGATCGTGCCCACATTTTAGAAGGATTACGTATCGCACTTGACCATATCGATGAAATTATTTCAACGATTCGTGAGTcagatacagat                                                              <  2:1226761/91‑1 (MQ=255)

                              tAAAGATCGTGCCCACATTTTAGAAGGATTACGTATCGCACTTGACCATATCGATGAAATTATTTCAACGATTCGTGAGTcagatacagat                                                              >  2:2614914/1‑91 (MQ=255)

                              tAAAGATCGTGCCCACATTTTAGAAGGATTACGTATCGCACTTGACCATATCGATGAAATTATTTCAACGATTCGTGAGTcagatacagat                                                              >  1:130818/1‑91 (MQ=255)

                                aaGATCGTGCCCACATTTTAGAAGGATTACGTATCGCACTTGACCATATCGATGAAATTATTTCAACGATTCGTGAGTCAGATACAGATaa                                                            <  2:2663893/91‑1 (MQ=255)

                                 aGATCGTGCCCACATTTTAGAAGGATTACGTATCGCACTTGACCATATCGATGAAATTATTTCAACGATTCGTGAGTCAGATACAGATaaa                                                           >  1:3256525/1‑91 (MQ=255)

                                  gATCGTGCCCACATTTTAGAAGGATTACGTATCGCACTTGACCATATCGATGAAATTATTTCAACGATTCGTGAGTCAGATACAGATAAAg                                                          <  1:1188347/91‑1 (MQ=255)

                                  gATCGTGCCCACATTTTAGAAGGATTACGTATCGCACTTGACCATATCGATGAAATTATTTCAACGATTCGTGAGTCAGATACAGATAAAg                                                          >  2:1527125/1‑91 (MQ=255)

                                     cGTGCCCACATTTTAGAAGGATTACGTATCGCACTTGACCATATCGATGAAATTATTTCAACGATTCGTGAGTCAGATACAGATAAAGTTg                                                       >  1:1667625/1‑91 (MQ=255)

                                     cGTGCCCACATTTTAGAAGGATTACGTATCGCACTTGACCATATCGATGAAATTATTTCAACGATTCGTGAGTCAGATACAGATAAAGTTg                                                       <  1:2366349/91‑1 (MQ=255)

                                      gTGCCCACATTTTAGAAGGATTACGTATCGCACTTGACCATATCGATGAAATTATTTCAACGATTCGTGAGTCAGATACAGATAAAg                                                          >  1:1256780/1‑87 (MQ=255)

                                      gTGCCCACATTTTAGAAGGATTACGTATCGCACTTGACCATATCGATGAAATTATTTCAACGATTCGTGAGTCAGATACAGATAAAGTTGc                                                      <  2:1447948/91‑1 (MQ=255)

                                      gTGCCCACATTTTAGAAGGATTACGTATCGCACTTGACCATATCGATGAAATTATTTCAACGATTCGTGAGTCAGATACAGATAAAGTTGc                                                      <  2:2182493/91‑1 (MQ=255)

                                        gCCCACATTTTAGAAGGATTACGTATCGCACTTGACCATATCGATGAAATTATTTCAACGATTCGTGAGTCAGATACAGATAAAGTTGCaa                                                    >  1:292553/1‑91 (MQ=255)

                                        gCCCACATTTTAGAAGGATTACGTATCGCACTTGACCATATCGATGAAATTATTTCAACGATTCGTGAGTCAGATACAGATAAAGTTGCaa                                                    <  1:274019/91‑1 (MQ=255)

                                        gCCCACATTTTAGAAGGATTACGTATCGCACTTCACCATATCGATGAAATTATTTCAACGATTCGTGAGTCAGATACAGATAAAGTTGCaa                                                    <  1:3177097/91‑1 (MQ=255)

                                          ccACATTTTAGAAGGATTACGTATCGCACTTGACCATATCGATGAAATTATTTCAACGATTCGTGAGTCAGATACAGATAAAGTTGCAATg                                                  <  2:1351666/91‑1 (MQ=255)

                                          ccACATTTTAGAAGGATTACGTATCGCACTTGACCATATCGATGAAATTATTTCAACGATTCGTGAGTCAGATACAGATAAAGTTGCAATg                                                  <  1:2777811/91‑1 (MQ=255)

                                            acaTTTTAGAAGGATTACGTATCGCACTTGACCATATCGATGAAATTATTTCAACGATTCGTGAGTCAGATACAGATAAAGTTGCAATGGa                                                >  2:177998/1‑91 (MQ=255)

                                                  tAGAAGGATTACGTATCGCACTTGACCATATCGATGAAATTATTTCAACGATTCGTGAGTCAGATACAGATAAAGTTGCAATGGAAAGCtt                                          >  2:3654828/1‑91 (MQ=255)

                                                   aGAAGGATTACGTATCGCACTTGACCATATCGATGAAATTATTTCAACGATTCGTGAGTCAGATACAGATAAAGTTGCAATGGAAAGCTTg                                         >  2:505695/1‑91 (MQ=255)

                                                    gAAGGATTACGTATCGCACTTGACCATATCGATGAAATTATTTCAACGATTCGTGAGTCAGATACAGATAAAGTTGCAATGGAAAGCTTga                                        >  2:3849078/1‑90 (MQ=255)

                                                    gAAGGATTACGTATCGCACTTGACCATATCGATGAAATTATTTCAACGATTCGTGAGTCAGATACAGATAAAGTTGCAATGGAAAGCTTGc                                        >  2:914373/1‑91 (MQ=255)

                                                    gAAGGATTACGTATCGCACTTGACCATATCGATGAAATTATTTCAACGATTCGTGAGTCAGATACAGATAAAGTTGCAATGGAAAGCTTGc                                        >  1:2272044/1‑91 (MQ=255)

                                                    gAAGGATTACGTATCGCACTTGACCATATCGATGAAATTAATTCAACGATTCGTGAGTCAGATACAGATAAAGTTGCAATGGAAAGCTTGc                                        >  2:3599195/1‑91 (MQ=255)

                                                     aaGGATTACGTATCGCACTTGACCATATCGATGAAATTATTTCAACGATTCGTGAGTCAGATACAGATAAAGTTGCAATGGAAAGCTTGca                                       >  1:779320/1‑91 (MQ=255)

                                                      aGGATTACGTATCGCACTTGACCATATCGATGAAATTATTTCAACGATTCGTGAGTCAGATACAGATAAAGTTGCAATGGAAAGCTTGcaa                                      <  2:38947/91‑1 (MQ=255)

                                                      aGGATTACGTATCGCACTTGACCATATCGATGAAATTATTTCAACGATTCGTGAGTCAGATACAGATAAAGTTGCAATGGAAAGCTTGcaa                                      <  1:2567209/91‑1 (MQ=255)

                                                       ggATTACGTATCGCACTTGACCATATCGATGAAATTATTTCAACGATTCGTGAGTCAGATACAGATAAAGTTGCAATGGAAAGCTTGcaac                                     <  2:1265545/91‑1 (MQ=255)

                                                       ggATTACGTATCGCACTTGACCATATCGATGAAATTATTTCAACGATTCGTGAGTCAGATACAGATAAAGTTGCAATGGAAAGCTTGcaac                                     <  2:2586730/91‑1 (MQ=255)

                                                         aTTACGTATCGCACTTGACCATATCGATGAAATTATTTCAACGATTCGTGAGTCAGATACAGATAAAGTTGCAATGGAAAGCTTGcaacaa                                   >  2:882695/1‑91 (MQ=255)

                                                             cGTATCGCACTTGACCATATCGATGAAATTATTTCAACGATTCGTGAGTCAGATACAGATAAAGTTGCAATGGAAAGCTTGCAACAACGCt                               >  1:3137920/1‑91 (MQ=255)

                                                             cGTATCGCACTTGACCATATCGATGAAATTATTTCAACGATTCGTGAGTCAGATACAGATAAAGTTGCAATGGAAAGCTTGCAACAACGCt                               <  2:718268/91‑1 (MQ=255)

                                                              gTATCGCACTTGACCATATCGATGAAATTATTTCAACGATTCGTGAGTCAGATACAGATAAAGTTGCAATGGAAAGCTTGCAACAACGCtt                              <  2:830622/91‑1 (MQ=255)

                                                              gTATCGCACTTGACCATATCGATGAAATTATTTCAACGATTCGTGAGTCAGATACAGATAAAGTTGCAATGGAAAGCTTGCAACAACGCtt                              >  1:427201/1‑91 (MQ=255)

                                                              gTATCGCACTTGACCATATCGATGAAATTATTTCAACGATTCGTGAGTCAGATACAGATAAAGTTGCAATGGAAAGCTTGCAACAACGCtt                              <  1:3843377/91‑1 (MQ=255)

                                                               tATCGCACTTGACCATATCGATGAAATTATTTCAACGATTCGTGAGTCAGATACAGATAAAGTTGCAATGGAAAGCTTGCAACAACGCTTc                             >  1:1073114/1‑91 (MQ=255)

                                                               tATCGCACTTGACCATATCGATGAAATTATTTCAACGATTCGTGAGTCAGATACAGATAAAGTTGCAATGGAAAGCTTGCAACAACGCTTc                             <  2:1651981/91‑1 (MQ=255)

                                                               tATCGCACTTGACCATATCGATGAAATTATTTCAACGATTCGTGAGTCAGATACAGATAAAGTTGCAATGGAAAGCTTGCAACAACGCTTc                             <  1:3565411/91‑1 (MQ=255)

                                                                aTCGCACTTGACCATATCGATGAAATTATTTCAACGATTCGTGAGTCAGATACAGATAAAGTTGCAATGGAAAGCTTGCAACAACGCTTCa                            <  2:2037712/91‑1 (MQ=255)

                                                                aTCGCACTTGACCATATCGATGAAATTATTTCAACGATTCGTGAGTCAGATACAGATAAAGTTGCAATGGAAAGCTTGCAACAACGCTTCa                            >  1:2934106/1‑91 (MQ=255)

                                                                 tCGCACTTGACCATATCGATGAAATTATTTCAACGATTCGTGAGTCAGATACAGATAAAGTTGCAATGGAAAGCTTGCAACAACGCTTCaa                           <  1:2702154/91‑1 (MQ=255)

                                                                   gCACTGGACCATATCGATGAAATTATTTCAACGATTCGTGAGTCAGATACAGATAAAGTTGCAATGGAAAGCTTGCAACAACGCTTCAAAc                         <  2:441557/91‑1 (MQ=255)

                                                                      cTTGACCATATCGATGAAATTATTTCAACGATTCGTGAGTCAGATACAGATAAAGTTGCAATGGAAAGCTTGCAACAACGCTTCAAACttt                      <  2:2161909/91‑1 (MQ=255)

                                                                      cTTGACCATATCGATGAAATTATTTCAACGATTCGTGAGTCAGATACAGATAAAGTTGCAATGGAAAGCTTGCAACAACGCTTCAAACttt                      <  1:1901054/91‑1 (MQ=255)

                                                                        tGACCATATCGATGAAATTATTTCAACGATTCGTGAGTCAGATACAGATAAAGTTGCAATGGAAAGCTTGCAACAACGCTTCAAACTTTCt                    <  2:2903482/91‑1 (MQ=255)

                                                                        tGACCATATCGATGAAATTATTTCAACGATTCGTGAGTCAGATACAGATAAAGTTGCAATGGAAAGCTTGCAACAACGCTTCAAACTTTCt                    <  1:1799979/91‑1 (MQ=255)

                                                                         gACCATATCGATGAAATTATTTCAACGATTCGTGAGTCAGATACAGATAAAGTTGCAATGGAAAGCTTGCAACAACGCTTCAAACTTTCTg                   <  1:3010651/91‑1 (MQ=255)

                                                                         gACCATATCGATGAAATTATTTCAACGATTCGTGAGTCAGATACAGATAAAGTTGCAATGGAAAGCTTGCAACAACGCTTCAAACTTTCTg                   <  2:1494340/91‑1 (MQ=255)

                                                                          aCCATATCGATGAAATTATTTCAACGATTCGTGAGTCAGATACAGATAAAGTTGCAATGGAAAGCTTGCAACAACGCTTCAAACTTTCTGa                  <  2:592527/91‑1 (MQ=255)

                                                                           ccATATCGATGAAATTATTTCAACGATTCGTGAGTCAGATACAGATAAAGTTGCAATGGAAAGCTTGCAACAACGCTTCAAACTTTCTGaa                 <  1:1974811/91‑1 (MQ=255)

                                                                             atatCGATGAAATTATTTCAACGATTCGTGAGTCAGATACAGATAAAGTTGCAATGGAAAGCTTGCAACAACGCTTCAAACTTTCTGaaaa               >  2:3262086/1‑91 (MQ=255)

                                                                              tatCGATGAAATTATTTCAACGATTCGTGAGTCAGATACAGATAAAGTTGCAATGGAAAGCTTGCAACAACGCTTCAAACTTTCTGaaaaa              >  1:2316543/1‑91 (MQ=255)

                                                                               atCGATGAAATTATTTCAACGATTCGTGAGTCAGATACAGATAAAGTTGCAATGGAAAGCTTGCAACAACGCTTCAAACTTTCTGAAAAAc             <  1:636209/91‑1 (MQ=255)

                                                                               atCGATAAAATTATTTCAACGATTCGTGAGTCAGATACAGATAAAGTTGCAATGGAAAGCTTGCAACAACGCTTCAAACTTTCTGAAAAAc             <  2:1754778/91‑1 (MQ=255)

                                                                                  gATGAAATTATTTCAACGATTCGTGAGTCAGATACAGATAAAGTTGCAATGGAAAGCTTGCAACAACGCTTCAAACTTTCTGAAAAAcaag          >  2:2250226/1‑91 (MQ=255)

                                                                                  gATGAAATTATTTCAACGATTCGTGAGTCAGATACAGATAAAGTTGCAATGGAAAGCTTGCAACAACGCTTCAAACTTTCTGAAAAAcaag          >  1:2070047/1‑91 (MQ=255)

                                                                                    tGAAATAATTTCATCGATTCGTGAGTCAGATACAGATAAAGTTGCAATGGAAAGCTTGCAACAACGCTTCAAACTTTCTGAAAAAcaagct        >  1:3784358/1‑91 (MQ=255)

                                                                                     gAAATTATTTCAACGATTCGTGAGTCAGATACAGATAAAGTTGCAATGGAAAGCTTGCAACAACGCTTCAAACTTTCTGAAAAAcaagctc       >  1:1322572/1‑91 (MQ=255)

                                                                                       aaTTATTTCAACGATTCGTGAGTCAGATACAGATAAAGTTGCAATGGAAAGCTTGCAACAACGCTTCAAACTTTCTGAAAAAcaagctcaa     >  2:847192/1‑91 (MQ=255)

                                                                                         ttattTCAACGATTCGTGAGTCAGATACAGATAAAGTTGCAATGGAAAGCTTGCAACAACGCTTCAAACTTTCTGAAAAAcaagctcaagc   <  1:1638915/91‑1 (MQ=255)

                                                                                         ttattTCAACGATTCGTGAGTCAGATACAGATAAAGTTGCAATGGAAAGCTTGCAACAACGCTTCAAACTTTCTGAAAAAcaagctcaagc   <  1:730549/91‑1 (MQ=255)

                                                                                         ttattTCAACGATTCGTGAGTCAGATACAGATAAAGTTGCAATGGAAAGCTTGCAACAACGCTTCAAACTTTCTGAAAAAcaagctcaagc   <  2:1701422/91‑1 (MQ=255)

                                                                                          tattTCAACGATTCGTGAGTCAGATACAGATAAAGTTGCAATGGAAAGCTTGCAACAACGCTTCAAACTTTCTGAAAAAcaagctcaagct  >  1:2002749/1‑91 (MQ=255)

AAGACGTACGCAATACAACTTACGTAAAGCTAAAGATCGTGCCCACATTTTAGAAGGATTACGTATCGCACTTGACCATATCGATGAAATTATTTCAACGATTCGTGAGTCAGATACAGATAAAGTTGCAATGGAAAGCTTGCAACAACGCTTCAAACTTTCTGAAAAACAAGCTCAAGCT  >  USA300TCH1516_ALE/8123‑8303

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