BRESEQ :: Evidence

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	USA300TCH1516_ALE	9,716	0	T	C	20.1%	382.4 / 95.2	184	intergenic (+21/+66)	gyrA/nnrD	DNA gyrase subunit A/ADP‑dependent (S)‑NAD(P)H‑hydrate dehydratase
Reads supporting (aligned to +/- strand): ref base T (71/76); new base C (19/18); total (90/94)
Fisher's exact test for biased strand distribution p-value = 8.54e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.06e-01
Rejected as consensus: Frequency below/above cutoff threshold.

Only 100 of 189 total aligned reads displayed.

TATTGAAGTAAGACAAGATTTCATGGATCGTGTTGAAGAAGATATACAACAATCATCAGATGAAGAATAATAAAAAATAAGACTTCCCTATATGTAGGGGAGTCTTA‑TTTTTATGCTAGAAAGTAATGCTTTACTATATTCAATGATTAGTAATGACTAACTTTCTAATTGTTTCATTGCG  >  USA300TCH1516_ALE/9626‑9806

tATTGAAGTAAGACAAGATTTCATGGATCGTGTTGAAGAAGATATACAACAATCATCAGATGAAGAATAATAAAAAATAAGACTTCCCtat                                                                                             >  1:2964702/1‑91 (MQ=255)

 aTTGAAGTAAGACAAGATTTCATGGATCGTGTTGAAGAAGATATACAACAATCATCAGATGAAGAATAATAAAAAATAAGACTTCCCtata                                                                                            >  2:315785/1‑91 (MQ=255)

 aTTGAAGTAAGACAAGATTTCATGGATCGTGTTGAAGAAGATATACAACAATCATCAGATGAAGAATAATAAAAAATAAGACTTCCCtata                                                                                            >  1:150315/1‑91 (MQ=255)

 aTTGAAGTAAGACAAGATTTCATGGATCGTGTTGAAGAAGATATACAACAATCATCAGATGAAGAATAATAAAAAATAAGACTTCCAtata                                                                                            >  2:1977422/1‑91 (MQ=255)

    gAAGTAAGACAAGATTTCATGGATCGTGTTGAAGAAGATATACAACAATCATCAGATGAAGAATAATAAAAAATAAGACTTCCCTATATGt                                                                                         >  2:956448/1‑91 (MQ=255)

    gAAGTAAGACAAGATTTCATGGATCGTGTTGAAGAAGATATACAACAATCATCAGATGAAGAATAATAAAAAATAAGACTTCCCTATATGt                                                                                         >  2:2506570/1‑91 (MQ=255)

      aGTAAGACAAGATTTCATGGATCGTGTTGAAGAAGATATACAACAATCATCAGATGAAGAATAATAAAAAATAAGACTTCCCTATATGTAg                                                                                       >  2:2122787/1‑91 (MQ=255)

      aGTAAGACAAGATTTCATGGATCGTGTTGAAGAAGATATACAACAATCATCAGATGAAGAATAATAAAAAATAAGACTTCCCTATATGTAg                                                                                       >  2:1245553/1‑91 (MQ=255)

       gTAAGACAAGATTTCATGGATCGTGTTGAAGAAGATATACAACAATCATCAGATGAAGAATAATAAAAAATAAGACTTCCCTATATGTAgg                                                                                      >  1:743268/1‑91 (MQ=255)

       gTAAGACAAGATTTCATGGATCGTGTTGAAGAAGATATACAACAATCATCAGATGAAGAATAATAAAAAATAAGACTCCCCTACATATAgg                                                                                      <  2:1736988/91‑1 (MQ=39)

        tAAGACAAGATTTCATGGATCGTGTTGAAGAAGATATACAACAATCATCAGATGAAGAATAATAAAAAATAAGACTTCCCTATATGTAggg                                                                                     <  2:2347684/91‑1 (MQ=255)

         aaGACAAGATTTCATGGATCGTGTTGAAGAAGATATACAACAATCATCAGATGAAGAATAATAAAAAATAAGACTCCCCTACATATAggga                                                                                    <  2:2140209/91‑2 (MQ=38)

           gACAAGATTTCATGGATCGTGTTGAAGAAGATATACAACAATCATCAGATGAAGAATAATAAAAAATAAGACTCCCCTACATATAgggaag                                                                                  <  1:1458210/91‑4 (MQ=37)

            aCAAGATTTCATGGATCGTGTTGAAGAAGATATACAACAATCATCAGATGAAGAATAATAAAAAATAAGACTTCCCTATATGTAGGGGAGt                                                                                 <  1:1187795/91‑1 (MQ=255)

            aCAAGATTTCATGGATCGTGTTGAAGAAGATATACAACAATCATCAGATGAAGAATAATAAAAAATAAGACTTCCCTATATGTAGGGGAGt                                                                                 >  2:1900707/1‑91 (MQ=255)

            aCAAGATTTCATGGATCGTGTTGAAGAAGATATACAACAATCATCAGATGAAGAATAATAAAAAATAAGACTCCCCTACATATAGGGAAGt                                                                                 <  1:1101994/91‑1 (MQ=37)

             cAAGATTTCATGGATTGTGTTGAAGAAGATATACAACAATCATCAGATGAAGAATAATAAAAAATAAGACTCCCCTACATATAGGGAAGTc                                                                                <  2:2943348/91‑1 (MQ=37)

             cAAGATTTCATGGATCGTGTTGAAGAAGATATACAACAATCATCAGATGAAGAATAATAAAAAATAAGACTTCCCTATATGTAGGGGAGTc                                                                                <  1:1282509/91‑1 (MQ=255)

              aaGATTTCATGGATCGTGTTGAAGAAGATATACAACAATCATCAGATGAAGAATAATAAAAAATAAGACTTCCCTATATGTAGGGGAGTCt                                                                               >  2:2379420/1‑91 (MQ=255)

              aaGATTTCATGGATCGTGTTGAAGAAGATATACAACAATCATCAGATGAAGAATAATAAAAAATAAGACTTCCCTATATGTAGGGGAGTCt                                                                               >  2:894137/1‑91 (MQ=255)

              aaGATTTCATGGATCGTGTTGAAGAAGATATACAACAATCATCAGATGAAGAATAATAAAAAATAAGACTCCCCTACATATAGGGAAGTCt                                                                               <  1:1877624/91‑1 (MQ=37)

              aaGATTTCATGGATCGTGTTGAAGAAGATATACAACAATCATCAGATGAAGAATAATAAAAAATAAGACTCCCCTACATATAGGGAAGTCt                                                                               <  2:1638109/91‑1 (MQ=37)

              aaGATTTCATGGATCGTGTTGAAGAAGATATACAACAATCATCAGATGAAGAATAATAAAAAATAAGACTCCCCTACATATAGGGAAGTCt                                                                               <  1:620388/91‑1 (MQ=37)

              aaGATTTCATGGATCGTGTTGAAGAAGATATACAACAATCATCAGATGAAGAATAATAAAAAATAAGACTCCCCTACATATAGGGAAGTCt                                                                               <  1:1878106/91‑1 (MQ=37)

                gATTTCATGGATCGTGTTGAAGAAGATATACAACAATCATCAGATGAAGAATAATAAAAAATAAGACTTCCCTATATGTAGGGGAGTCTTa                                                                             >  2:15149/1‑91 (MQ=255)

                gATTTCATGGATCGTGTTGAAGAAGATATACAACAATCATCAGATGAAGAATAATAAAAAATAAGACTTCCCTATATGCAGGGGCGTCttc                                                                             >  2:1607582/1‑90 (MQ=255)

                 aTTTCATGGATCGTGTTGAAGAAGATATACAACAATCATCAGATGAAGAATAATAAAAAATAAGACTTCCCTATATGTAGGGGAGTCTTA‑t                                                                           >  2:2195073/1‑91 (MQ=255)

                 aTTTCATGGATCGTGTTGAAGAAGATATACAACAATCATCAGATGAAGAATAATAAAAAATAAGACTCCCCTACATATAGGGAAGTCTTA‑t                                                                           <  2:870149/91‑1 (MQ=25)

                 aTTTCATGGATCGTGTTGAAGAAGATATACAACAATCATCAGATGAAGAATAATAAAAAATAAGACTCCCCTACATATAGGGAAGTCTTA‑t                                                                           <  2:1065514/91‑1 (MQ=25)

                   ttCATGGATCGTGTTGAAGAAGATATACAACAATCATCAGATGAAGAATAATCAAAAATAAGACTTCCCTATATGTAGGGGAGTCTTA‑ttt                                                                         >  2:2336594/1‑91 (MQ=255)

                    tCATGGATCGTGTTGAAGAAGATATACAACAATCATCAGATGAAGAATAATAAAAAATAAGACTTCCCTATATGTAGGGGAGTCTTA‑tttt                                                                        <  1:874718/91‑1 (MQ=255)

                    tCATGGATCGTGTTGAAGAAGATATACAACAATCATCAGATGAAGAATAATAAAAAATAAGACTTCCCTATATGTAGGGGAGTCTTA‑tttt                                                                        >  2:416954/1‑91 (MQ=255)

                      aTGGATCGTGTTGAAGAAGATATACAACAATCATCAGATGAAGAATAATAAAAAATAAGACTTCCCTATATGTAGGGGAGTCTTA‑TTTTTa                                                                      <  1:2561010/91‑1 (MQ=255)

                      aTGGATCGTGTTGAAGAAGATATACAACAATCATCAGATGAAGAATAAT‑AAAAATAAGACTCCCCTACATATAGGGAAGTCTTA‑tttttta                                                                     <  2:728201/91‑3 (MQ=18)

                           tCGTGTTGAAGAAGATATACAACAATCATCAGATGAAGAATAATAAAAAATAAGACTTCCCTATATGTAGGGGAGTCTTA‑TTTTTATGCTa                                                                 <  1:164103/91‑1 (MQ=255)

                            cGTGTTGAAGAAGATATACAACAATCATCAGATGAAGAATAATAAAAAATAAGACTTCCCTATATGTAGGGGAGTCTTA‑TTTTTATGCTAg                                                                <  2:379775/91‑1 (MQ=255)

                                   aagaagATATACAACAATCATCAGATGAAGAATAATAAAAAATAAGACTTCCCTATATGTAGGGGAGTCTTA‑TTTTTATGCTAGAAAGTaa                                                         <  1:1312809/91‑1 (MQ=255)

                                   aagaagATATACAACAATCATCAGATGAAGAATAATAAAAAATAAGACTTCCCTATATGTAGGGGAGTCTTA‑TTTTTATGCTAGAAAGTaa                                                         <  1:2003014/91‑1 (MQ=255)

                                     gaagaTATACAACAATCATCAGATGAAGAATAATAAAAAATAAGACTTCCCTATATGTAGGGGAGTCTTA‑TTTTTATGCTAGAAAGTAATg                                                       <  1:410737/91‑1 (MQ=255)

                                         atatACAACAATCATCAGATGAAGAATAATAAAAAATAAGACTTCCCTATATGTAGGGGAGTCTTA‑TTTTTATGCTAGAAAGTAATGCttt                                                   >  1:952885/1‑91 (MQ=255)

                                           ataCAACAATCATCAGATGAAGAATAATAAAAAATAAGACTTCCCTATATGTAGGGGAGTCTTA‑TATTTATGCTAGAAAGTAATGCTTTAc                                                 >  2:2917483/1‑91 (MQ=255)

                                            taCAACAATCATCAGATGAAGAATAATAATAAATAAGACTTCCCTATATGTAGGGGAGTCTTA‑TTTTTATGCTAGAAAGTAATGCTTTACt                                                <  2:2226422/91‑1 (MQ=255)

                                            taCAACAATCATCAGATGAAGAATAATAAAAAATAAGACTTCCCTATATGTAGGGGAGTCTTA‑TTTTTATGCTAGAAAGTAATGCTTTACt                                                >  1:1696971/1‑91 (MQ=255)

                                            taCAACAATCATCAGATGAAGAATAATAAAAAATAAGACTTCCCTATATGTAGGGGAGTCTTA‑TTTTTATGCTAGAAAGTAATGCTTTACt                                                >  1:2312550/1‑91 (MQ=255)

                                            taCAACAATCATCAGATGAAGAATAATAAAAAATAAGACTTCCCTATATGTAGGGGAGTCTTA‑TTTTTATGCTAGAAAGTAATGCTTTACt                                                >  1:247954/1‑91 (MQ=255)

                                               aacaaTCATCAGATGAAGAATAATAAAAAATAAGACTTCCCTATATGTAGGGGAGTCTTA‑TTTTTATGCTAGAAAGTAATGCTTTACtata                                             <  1:444177/91‑1 (MQ=255)

                                               aacaaTCATCAGATGAAGAATAATAAAAAATAAGACTTCCCTATATGTAGGGGAGTCCTA‑TTTTTATGCTAGAAAGTAATGCTTTACtata                                             >  2:1610180/1‑91 (MQ=255)

                                                ccaaTCATCAGATGAAGAATAATAAAAAATAAGACTTCCCTATATGTAGGGGAGTCTTA‑TTTTTATGCTAGAAAGTAATGCTTTACTatat                                            <  2:3079101/90‑1 (MQ=255)

                                                acaaTCATCAGATGAAGAATAATAAAAAATAAGACTTCCCTATATGTAGGGGAGTATTA‑TTTTTATGCTAGAAAGTAATGCTTTACTatat                                            >  2:271650/1‑91 (MQ=255)

                                                 caaTCATCAGATGAAGAATAATAAAAAATAAGACTTCCCTATATGTAGGGGAGTCTTA‑TTTTTATGCTAGAAAGTAATGCTTTACTATAtt                                           >  2:403101/1‑91 (MQ=255)

                                                  aaTCATCAGATGAAGAATAATAAAAAATAAGACTTCCCTATATGTAGGGGAGTCTTA‑TTTTTATGCTAGAAAGTAATGCTTTACTATATTc                                          <  2:217909/91‑1 (MQ=255)

                                                  aaTCATCAGATGAAGAATAATAAAAAATAAGACTTCCCTATATGTAGGGGAGTCTTA‑TTTTTATGCTAGAAAGTAATGCTTTACTATATTc                                          <  2:970340/91‑1 (MQ=255)

                                                   atcatcAGATGAAGAATAATAAAATATAAGACTTCCCTATATGTAGGGGAGTCTTA‑TTTTTATGCTAGAAAGTAATGCTTTACTATATTCa                                         <  2:141766/91‑1 (MQ=255)

                                                   atcatcAGATGAAGAATAATAAAAAATAAGACTTCCCTATATGTAGGGGAGTCTTA‑TTTTTATGCTAGAAAGTAATGCTTTACTATATTCa                                         <  2:2463410/91‑1 (MQ=255)

                                                   atcatcAGATGAAGAATAATAAAAAATAAGACTTCCCTATATGTAGGGGAGTCTTA‑TTTTTATGCTAGAAAGTAATGCTTTACTATATTCa                                         <  1:1868588/91‑1 (MQ=255)

                                                    tcatcaGATGAAGAATAATAAAAAATAAGACTTCCCTATATGTAGGGGAGTCTTA‑TTTTTATGCTAGAAAGTAATGCTTTACTATATTCaa                                        <  1:2566995/91‑1 (MQ=255)

                                                     catcaGATGAAGAATAATAAAAAATAAGACATCCCCATACGCAGGGGAGTCTTA‑CCATTATGCTAGAACGCAATGCTTTACGATATTCAAt                                       >  2:626658/1‑91 (MQ=255)

                                                      atcaGATGAAGAATAATAAAAAATAAGACTTCCCTATATGTAGGGGAGTCTTA‑TTTTTATGCTAGAAAGTAATGCTTTACTATATTCAATg                                      <  1:2828512/91‑1 (MQ=255)

                                                      atcaGATGAAGAATAATAAAAAATAAGACTTCCCTATATGTAGGGGAGTCTTA‑TTTTTATGCTAGAAAGTAATGCTTTACTATATTCAATg                                      >  2:1009173/1‑91 (MQ=255)

                                                      atcaGATGAAGAATAATAAAAAATAAGACTTCCCTATATGTAGGGGAGTCTTA‑TTTTTATGCTAGAAAGTAATGCTTTACTATATTCAATg                                      <  2:2245232/91‑1 (MQ=255)

                                                      atcaGATGAAGAATAATAAAAAATAAGACTTCCCTATATGTAGGGGAGTCTTA‑GTTTTATGCTAGAAAGTAATGCTTTACTATATTCAATg                                      <  1:2623238/91‑1 (MQ=255)

                                                       tcaGATGAAGAATAATAAAAAATAAGACTTCCCTATATGTAGGGGAGTCTTA‑TTTTTATGCTAGAAAGTAATGCTTTACTATATTCAATGa                                     >  1:1248726/1‑91 (MQ=255)

                                                       tcaGATGAAGAATAATAAAAAATAAGACTTCCCTATATGTAGGGGAGTCTTA‑TTTTTATGCTAGAAAGTAATGCTTTACTATATTCAATGa                                     >  2:577626/1‑91 (MQ=255)

                                                        caGATGAAGAATAATAAAAAATAAGACTTCCCTATATGTAGGGGAGTCTTA‑TTTTTATGCTAGAAAGTAATGCTTTACTATATTCAATGAt                                    <  1:694870/91‑1 (MQ=255)

                                                        caGATGAAGAATAATAAAAAATAAGACTTCCCTATATGTAGGGGAGTCTTA‑TTTTTATGCTAGAAAGTAATGCTTTACTATATTCAATGAt                                    >  1:212018/1‑91 (MQ=255)

                                                        caGATGAAGAATAATAAAAAATAAGACTTCCCTATATGTAGGGGAGTCTTA‑TTTTTATGCTAGAAAGTAATGCTTTACTATATTCAATGAt                                    >  2:1524657/1‑91 (MQ=255)

                                                         aGATGAAGAATAATAAAAAATAAGACTTCCCTATATGTAGGGGAGTCTTA‑TTTTTATGCTAGAAAGTAATGCTTTACTATATTCAATGAtt                                   <  2:2615504/91‑1 (MQ=255)

                                                             gaagaaTAATAAAAAATAAGACTTCCCTATATGTAGGGGAGTCATA‑TTTTTATGCTAGAAAGTAATGCTTTACTATATTCAATGATTAGTa                               >  2:2394942/1‑91 (MQ=255)

                                                              aagaaTAATAAAAAATAAGACTTCCCTATATGTAGGGGAGTCTTA‑TTTTTATGCTAGAAAGTAATGCTTTACTATATTCAATGATTAGTaa                              <  2:1376348/91‑1 (MQ=255)

                                                                gaaTAATAAAAAATAAGACTTCCCTATATGTAGGGGAGTCTTA‑TTTTTATGCTAGAAAGTAATGCTTTACTATATTCAATGATTAGTCATg                            >  1:2491649/1‑91 (MQ=255)

                                                                gaaTAATAAAAAATAAGACTTCCCTATATGTAGGGGAGTCTTA‑TTTTTATGCTAGAAAGTAATGCTTTACTATATTCAATGATTAGTAATg                            <  2:2092080/91‑1 (MQ=255)

                                                                   taataaAAAATAAGACTTCCCTATATGTAGGGGAGTCTTA‑TTTTTATGCTAGAAAGTAATGCTTTACTATATTCAATGATTAGTAATGACt                         >  1:1025314/1‑91 (MQ=255)

                                                                    aataaAAAATAAGACTTCCCTATATGTAGGGGAGTCTTA‑TTTTTATGCTAGAAAGTAATGCTTTACTATATTCAATGATTAGTAATGACTa                        <  1:1397204/91‑1 (MQ=255)

                                                                      taaAAAATAAGACTTCCCTATATGTAGGGGAGTCTTA‑TTTTTATGCTAGAAAGTAATGCTTTACTATATTCAATGATTAGTAATGACTAAc                      <  2:2209234/91‑1 (MQ=255)

                                                                      taaAAAATAAGACTTCCCTATATGTAGGGGAGTCTTA‑TTTTTATGCTAGAAAGTAATGCTTTACTATATTCAATGATTAGTAATGACTAAc                      <  2:3014493/91‑1 (MQ=255)

                                                                       aaaaaaTAAGACTTCCCTATATGTAGGGGAGTCTTA‑TTTTTATGCTAGAAAGTAATGCTTTACTACATTCAATGATTAGTAATGACTaaat                     >  2:2024338/1‑89 (MQ=255)

                                                                        aaaaaTAAGACTTCCCTATATGTAGGGGAGTCTTA‑TTTTTATGCTAGAAAGTAATGCTTTACTATATTCAATGATTAGTAATGACTAACtt                    <  2:1275902/91‑1 (MQ=255)

                                                                           aaTAAGACTTCCCTATATGTAGGGGAGTCTTA‑TTTTTATGCTAGAAAGTAATGCTTTACTATATTCAATGATTAGTAATGACTAACTTTCt                 <  2:2365738/91‑1 (MQ=255)

                                                                            aTAAGACTTCCCTATATGTAGGGGAGTCTTA‑TTTTTATGCTAGAAAGTAATGCTTTACTATATTCAATGATTAGTAATGACTAACTTTCTa                <  2:1329747/91‑1 (MQ=255)

                                                                            aTAAGACTTCCCTATATGTAGGGGAGTCTTA‑TTTTTATGCTAGAAAGTAATGCTTTACTATATTCAATGATTAGTAATGACTAACTTTCTa                >  2:2669996/1‑91 (MQ=255)

                                                                            aTAAGACTTCCCTATATGTAGGGGAGTCTTA‑TTTTTATGCTAGAAAGTAATGCTTTACTATATTCAATGATTAGTAATGACTAACTTTCTa                <  2:81783/91‑1 (MQ=255)

                                                                             tAAGACTTCCCTATATGTAGGGGAGTCTTA‑TTTTTATGCTAGAAAGTAATGCTTTACTATATTCAATGATTAGTAATGACTAACTTTCTaa               <  1:1340771/91‑1 (MQ=255)

                                                                             tAAGACTTCCCTATATGTAGGGGAGTCTTA‑TTTTTATGCTAGAAAGTAATGCTTTACTATATTCAATGATTAGTAATGACTAACTTTCTaa               >  1:3053448/1‑91 (MQ=255)

                                                                             tAAGACTCCCCTACATATAGGGAAGTCTTA‑TTTTTATGCTAGAAAGTAATGCTTTACTATATTCAATGATTAGTAATGACTAACTTTCTaa               >  1:136733/1‑91 (MQ=37)

                                                                             tAAGACTCCCCTACATATAGGGAAGTCTTA‑TTTTTATGCTAGAAAGTAATGCTTTACTATATTCAATGATTAGTAATGACTAACTTTCTaa               >  1:2373825/1‑91 (MQ=37)

                                                                              aaGACTTCCCTATATGTAGGGGAGTCTTA‑TTTTTATGCTAGAAAGTAATGCTTTACTATATTCAATGATTAGTAATGACTAACTTTCTAAt              <  1:773/91‑1 (MQ=255)

                                                                              aaGACTCCCCTACATATAGGGAAGTCTTATTTTTTATGCTAGAAAGTAATGCTTTACTATATTCAATGATTAGTAATGACTAACTTTCTaa               >  1:1757104/1‑91 (MQ=25)

                                                                              aaGACTCCCCTACATATAGGGAAGTCTTA‑TTTTTATGCTAGAAAGTAATGCTTTACTATATTCAATGATTAGTAATGACTAACTTTCTAAt              >  2:2782583/1‑91 (MQ=37)

                                                                               aGACTCCCCTACATATAGGGAAGTCTTATTTTTTATGCTAGAAAGTAATGCTTTACTATATTCAATGATTAGTAATGACTAACTTTCTAAt              >  1:1215355/1‑91 (MQ=21)

                                                                               aGACTCCCCTACATATAGGGAAGTCTTA‑TTTTTATGCTAGAAAGTAATGCTTTACTATATTCAATGATTAGTAATGACTAACTTTCTAAtt             >  1:2706525/1‑91 (MQ=37)

                                                                               aGACTCCCCTACATATAGGGAAGTCTTA‑TTTTTATGCTAGAAAGTAATGCTTTACTATATTCAATGATTAGTAATGACTAACTTTCTAAtt             >  1:171720/1‑91 (MQ=37)

                                                                               aGACTCCCCTACATATAGGGAAGTCTTA‑TTTTTATGCTAGAAAGTAATGCTTTACTATATTCAATGATTAGTAATGACTAACTTTCTAAtt             >  1:572092/1‑91 (MQ=37)

                                                                               aGACTCCCCTACATATAGGGAAGCCTTATTTTTTATGCTAGAAAGTAATGCTTTACTATATTCAATGATTAGTAATGACTAACTTTCTAAt              >  2:53168/1‑91 (MQ=25)

                                                                                 actccccTACATATAGGGAAGTCTTA‑TTTTTATGCTAGAAAGTAATGCTTTACTATATTCAATGATTAGTAATGACTAACTTTCTAATTGt           >  2:745446/5‑91 (MQ=37)

                                                                                 aCTTCCCTATATGTAGGGGAGTCTTA‑TTTTTATGCTAGAAAGTAATGCTTTACTATATTCAATGATTAGTAATGACTAACTTTCTAATTGt           <  1:2444914/91‑1 (MQ=255)

                                                                                   ttCCCTATATGTAGGGGAGTCTTA‑TTTTTATGCTAGAAAGTAATGCTTTACTATATTCAATGATTAGTAATGACTAACTTTCTAATTGttt         <  1:3026676/91‑1 (MQ=255)

                                                                                   tccccTACATATAGGGAAGTCTTA‑TTTTTATGCTAGAAAGTAATGCTTTTCTATATTCAATGATTAGTAATGACTAACTTTCTAATTGttt         >  2:2822642/3‑91 (MQ=37)

                                                                                   tccccTACATATAGGGAAGTCTTA‑TTTTTATGCTAGAAAGTAATGCTTTACTATATTCAATGATTAGTAATGACTAACTTTCTAATTGttt         >  1:2879028/3‑91 (MQ=38)

                                                                                    tCCCTATATGTAGGGGAGTCTTA‑TTTTTATGCTAGAAAGTAATGCTTTACTATATTCAATGATTAGTAATGACTAACTTTCTAATTGTTTc        <  1:1001830/91‑1 (MQ=255)

                                                                                          tatGTAGGGGAGTCTTA‑TTTTTATGCTAGAAAGTAATGCTTTACTATATTCAATGATTAGTAATGACTAACTTTCTAATTGTTTCATTGCg  <  2:192223/91‑1 (MQ=255)

TATTGAAGTAAGACAAGATTTCATGGATCGTGTTGAAGAAGATATACAACAATCATCAGATGAAGAATAATAAAAAATAAGACTTCCCTATATGTAGGGGAGTCTTA‑TTTTTATGCTAGAAAGTAATGCTTTACTATATTCAATGATTAGTAATGACTAACTTTCTAATTGTTTCATTGCG  >  USA300TCH1516_ALE/9626‑9806

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