BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	USA300TCH1516_ALE	942,221	T→A	7.9%	intergenic (‑192/+57)	ptlE ← / ← mnhG1	Neopentalenolactone D synthase/Na(+)/H(+) antiporter subunit G1

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	USA300TCH1516_ALE	942,221	0	T	A	7.9%	344.0 / 17.2	101	intergenic (‑192/+57)	ptlE/mnhG1	Neopentalenolactone D synthase/Na(+)/H(+) antiporter subunit G1
Reads supporting (aligned to +/- strand): ref base T (45/48); new base A (3/5); total (48/53)
Fisher's exact test for biased strand distribution p-value = 7.18e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.97e-01

Only 100 of 101 total aligned reads displayed.

AGATAGGGCGAATCTGTAGTCAAAGAAAAAATCATTGAAAAAGAGGTAACAATGTCAAAAGAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTGTTATATAATCTACAATTTTGTGTCTTTTAAGTCTTCCGAA  >  USA300TCH1516_ALE/942133‑942308

aGATAGGGCGAATCTGTAGTCAAAGAAAAAATCATTGAAAAAGAGGTAACAATGTCAAAAGAAAACAGCAGTAAAATCATTCCTAATTTgg                                                                                       <  1:2772482/91‑1 (MQ=255)

   tAGGGCGAATCTGTAGTCAAAGAAAAAATCATTGAAAAAGAGGTAACAATGTCAAAAGAAAACAGCAGTAAAATCATTCCTAATTTGGAat                                                                                    >  1:3414538/1‑91 (MQ=255)

    aGGGCGAATCTGTAGTCAAAGAAAAAATCATTGAAAAAGAGGTAACAATGTCAAAAGAAAACAGCAGTAAAATCATTCCTAATTTGGAatc                                                                                   >  1:3601013/1‑91 (MQ=255)

     gggCGAATCTGTAGTCAAAGAAAAAATCATTGAAAAAGAGGTAACAATGTCAAAAGAAAACAGCAGTAAAATCATTCCTAATTTGGAatca                                                                                  >  1:2796726/1‑91 (MQ=255)

     gggCGAATCTGTAGTCAAAGAAAAAATCATTGAAAAAGAGGTAACAATGTCAAAAGAAAACAGCAGTAAAATCATTCCTAATTTGGAatca                                                                                  >  1:3021717/1‑91 (MQ=255)

      ggCGAATCTGTAGTCAAAGAAAAAATCATTGAAAAAGAGGTAACAATGTCAAAAGAAAACAGCAGTAAAATCATTCCTAATTTGGAatcat                                                                                 >  1:198663/1‑91 (MQ=255)

       gCGAATCTGTAGTCAAAGAAAAAATCATTGAAAAAGAGGTAACAATGTCAAAAGAAAACAGCAGTAAAATCATTCCTAATTTGGAatcatc                                                                                >  1:104701/1‑91 (MQ=255)

         gAATCTGTAGTCAAAGAAAAAATCATTGAAAAAGAGGTAACAATGTCAAAAGAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCtt                                                                              >  1:533918/1‑91 (MQ=255)

         gAATCTGTAGTCAAAGAAAAAATCATTGAAAAAGAGGTAACAATGTCAAAAGAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCtt                                                                              <  1:1489511/91‑1 (MQ=255)

           aTCTGTAGTCAAAGAAAAAATCATTGAAAAAGAGGTAACAATGTCAAAATAAAACAGCAGTAAGATGATTCCAAATTAGGAATGATTTTAc                                                                            <  1:3598812/91‑1 (MQ=18)

            tCTGTAGTCAAAGAAAAAATCATTGAAAAAGAGGTAACAATGTCAAAAGAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCTTAct                                                                           >  1:1225609/1‑91 (MQ=255)

              tGTAGTCAAAGAAAAAATCATTGAAAAAGAGGTAACAATGTCAAAAGAAAACAGCAGTAAGATGATTCCAAATTAGGAATGATTTTActgc                                                                         <  1:3680518/91‑1 (MQ=18)

              tGTAGTCAAAGAAAAAATCATTGAAAAAGAGGTAACAATGTCAAAAGAAAACAGCAGTAAGATGATTCCAAATTAGGAATGATTTTActgc                                                                         <  1:1700773/91‑1 (MQ=18)

              tGTAGTCAAAGAAAAAATCATTGAAAAAGAGGTAACAATGTCAAAAGAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCTTActgc                                                                         >  1:1985952/1‑91 (MQ=255)

               gTAGTCAAAGAAAAAATCATTGAAAAAGAGGTAACAATGTCAAAAGAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCTTActgct                                                                        <  1:4827/91‑1 (MQ=255)

               gTAGTCAAAGAAAAAATCATTGAAAAAGAGGTAACAATGTCAAAAGAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCTTActgct                                                                        >  1:2490894/1‑91 (MQ=255)

                tAGTCAAAGAAAAAATCATTGAAAAAGAGGTAACAATGTCAAAAGAAAACAGCAGTAAGATGATTCCAAATTAGGAATGATTTTActgctg                                                                       <  1:753994/91‑1 (MQ=18)

                tAGTCAAAGAAAAAATCATTGAAAAAGAGGTAACAATGTCAAAAGAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCTTActgctg                                                                       >  1:1266410/1‑91 (MQ=255)

                   tCAAAGAAAAAATCATTGAAAAAGAGGTAACAATGTCAAAAGAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCTTACTGCtgttt                                                                    >  1:875880/1‑91 (MQ=255)

                   tCAAAGAAAAAATCATTGAAAAAGAGGTAACAATGTCAAAAGAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCTTACTGCtgttt                                                                    <  1:3919152/91‑1 (MQ=255)

                   tCAAAGAAAAAATCATTGAAAAAGAGGTAACAATGTCAAAAGAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCTTACTGCtgttt                                                                    <  1:3919181/91‑1 (MQ=255)

                   tCAAAGAAAAAATCATTGAAAAAGAGGTAACAATGTCAAAAGAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCTTACTGCtgttt                                                                    <  1:2032707/91‑1 (MQ=255)

                    cAAAGAAAAAATCATTGAAAAAGAGGTAACAATGTCAAAAGAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCTTACTGCtgtttg                                                                   <  1:350009/91‑1 (MQ=255)

                     aaaGAAAAAATCATTGAAAAAGAGGTAACAATGTCAAAAGAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCTTACTGCtgtttgt                                                                  <  1:535250/91‑1 (MQ=255)

                        gAAAAAATCATTGAAAAAGAGGTAACAATGTCAAAAGAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTttgttgt                                                               <  1:3657406/91‑1 (MQ=255)

                        gAAAAAATCATTGAAAAAGAGGTAACAATGTCAAAAGAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTttgttgt                                                               <  1:1607426/91‑1 (MQ=255)

                           aaaaTCATTGAAAAAGAGGTAACAATGTCAAAAGAAAACAGCAGTAAAATCATTTCTAATTTGGAATCATCTTACTGCTGTTTGATGTTGa                                                            >  1:458101/1‑91 (MQ=255)

                           aaaaTCATTGAAAAAGAGGTAACAATGTCAAAAGAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGa                                                            <  1:1060791/91‑1 (MQ=255)

                           aaaaTCATTGAAAAAGAGGTAACAATGTCAAAAGAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGa                                                            <  1:2815237/91‑1 (MQ=255)

                           aaaaTCATTGAAAAAGAGGTAACAATGTCAAAAGAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGa                                                            <  1:253345/91‑1 (MQ=255)

                            aaaTCATTGAAAAAGAGGTAACAATGTCAAAAGAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGAt                                                           >  1:1684026/1‑91 (MQ=255)

                            aaaTCATTGAAAAAGAGGTAACAATGTCAAAAGAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGAt                                                           >  1:295199/1‑91 (MQ=255)

                             aaTCATTGAAAAAGAGGTAACAATGTCAAAAGAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGAtt                                                          >  1:16992/1‑91 (MQ=255)

                             aaTCATTGAAAAAGAGGTAACAATGTCAAAAGAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGAtt                                                          >  1:16554/1‑91 (MQ=255)

                              aTCATTGAAAAAGAGGTAACAATGTCAAAAGAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGAttt                                                         >  1:3301526/1‑91 (MQ=255)

                              aTCATTGAAAAAGAGGTAACAATGTCAAAAGAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGAttt                                                         >  1:2822081/1‑91 (MQ=255)

                               tCATTGAAAAAGAGGTAACAATGTCAAAAGAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTta                                                        <  1:1870451/91‑1 (MQ=255)

                                   tGAAAAAGAGGTAACAATGTCAAAAGAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTTATAtt                                                    <  1:2916775/91‑1 (MQ=255)

                                   tGAAAAAGAGGTAACAATGTCAAAAGAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTTATAtt                                                    <  1:3521870/91‑1 (MQ=255)

                                    gAAAAAGAGGTAACAATGTCAAAAGAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTTATATTc                                                   <  1:2017584/91‑1 (MQ=255)

                                    gAAAAAGAGGTAACAATGTCAAAAGAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTTATATTc                                                   <  1:2779837/91‑1 (MQ=255)

                                      aaaaGAGGTAACAATGTCAAAAGAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTTATATTCAt                                                 <  1:2588816/91‑1 (MQ=255)

                                       aaaGAGGTAACAATGTCAAAAGAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTTATATTCATg                                                <  1:2199202/91‑1 (MQ=255)

                                         agagGTAACAATGTCAAAAGAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTTATATTCATGAt                                              <  1:2337100/91‑1 (MQ=255)

                                           agGTAACAATGTCAAAAGAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTTATATTCATGAttt                                            >  1:3705692/1‑91 (MQ=255)

                                              tAACAATGTCAAAAGAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTGt                                         >  1:2498002/1‑91 (MQ=255)

                                               aaCAATGTCAAAAGAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTGtt                                        <  1:2573813/91‑1 (MQ=255)

                                               aaCAATGTCAAAAGAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTGtt                                        >  1:3721588/1‑91 (MQ=255)

                                                aCAATGTCAAAATAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTGTta                                       >  1:3189444/1‑91 (MQ=255)

                                                aCAATGTCAAAATAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTGTta                                       >  1:3189429/1‑91 (MQ=255)

                                                aCAATGTCAAAAGAAATCAGCAGTAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTGTta                                       <  1:4120684/91‑1 (MQ=255)

                                                aCAATGTCAAAAGAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTGTta                                       >  1:507051/1‑91 (MQ=255)

                                                aCAATGTCAAAAGAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTGTta                                       >  1:506879/1‑91 (MQ=255)

                                                      tCAAAAGAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTGTTATATAAt                                 >  1:2553355/1‑91 (MQ=255)

                                                      tCAAAAGAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTGTTATATAAt                                 >  1:3156022/1‑91 (MQ=255)

                                                       cAAAAGAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTGTTATATAATc                                >  1:1488563/1‑91 (MQ=255)

                                                       cAAAAGAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTGTTATATAATc                                <  1:1181091/91‑1 (MQ=255)

                                                       cAAAAGAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTGTTATATAATc                                <  1:1307216/91‑1 (MQ=255)

                                                       cAAAAGAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTGTTATATAATc                                <  1:3058748/91‑1 (MQ=255)

                                                       cAAAAGAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTGTTATATAATc                                <  1:2702648/91‑1 (MQ=255)

                                                         aaaGAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTGTTATATAATCTa                              <  1:3445736/91‑1 (MQ=255)

                                                          aaGAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTGTTATATAATCTAc                             >  1:2800365/1‑91 (MQ=255)

                                                           aGAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTGTTATATAATCTACa                            >  1:3711113/1‑91 (MQ=255)

                                                            gAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTGTTATATAATCTACaa                           <  1:2708133/91‑1 (MQ=255)

                                                             aaaaCAGCAGTAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTGTTATATAATCTACAAt                          <  1:563914/91‑1 (MQ=255)

                                                              aaaCAGCAGTAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTGTTATATAATCTACAAtt                         >  1:2444220/1‑91 (MQ=255)

                                                              aaaCAGCAGTAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTTTGGTGTTGATTTATATTCATGATTTTGTTATATAATCTACAAtt                         <  1:2395792/91‑1 (MQ=255)

                                                               aaCAGCAGTAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTGTTATATAATCTACAAttt                        >  1:1968583/1‑91 (MQ=255)

                                                               aaCAGCAGTAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTGTTATATAATCTACAAttt                        <  1:636692/91‑1 (MQ=255)

                                                               aaCAGCAGTAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTGTTATATAATCTACAAttt                        <  1:635570/91‑1 (MQ=255)

                                                               aaCAGCAGTAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTTGGTTGTTGATTTATATTCATGATTTTGTTATATAATCTACAAttt                        >  1:3297950/1‑91 (MQ=255)

                                                                 cagcagTAAGATGATTCCAAATTAGGAATGATTTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTGTTATATAATCTACAATTTtg                      >  1:517423/1‑91 (MQ=18)

                                                                 cagcagTAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTGTTATATAATCTACAATTTtg                      <  1:1905627/91‑1 (MQ=255)

                                                                 cagcagTAAAATAATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTGTTATATAATCTACAATTTtg                      >  1:3178549/1‑91 (MQ=255)

                                                                  agcagTAAGATGATTCCAAATTAGGAATGATTTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTGTTATATAATCTACAATTTtgt                     >  1:2833859/1‑91 (MQ=18)

                                                                   gcagTAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTGTTATATAATCTACAATTTtgtg                    >  1:443696/1‑91 (MQ=255)

                                                                     agTAAGATGATTCCAAATTAGGAATGATTTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTGTTATATAATCTACAATTTTGTGTc                  >  1:321213/1‑91 (MQ=18)

                                                                      gTAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTGTTATATAATCTACAATTTTGTGTCt                 >  1:3069051/1‑91 (MQ=255)

                                                                       tAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTGTTATATAATCTACAATTTTGTGTCtt                >  1:103929/1‑91 (MQ=255)

                                                                       tAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTGTTATATAATCTACAATTTTGTGTCtt                <  1:3966882/91‑1 (MQ=255)

                                                                       tAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTGTTATATAATCTACAATTTTGTGTCtt                <  1:4009350/91‑1 (MQ=255)

                                                                       tAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTGTTATATAATCTACAATTTTGTGTCtt                <  1:4009488/91‑1 (MQ=255)

                                                                       caaaaTCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTGTTATATAATCTACAATTTTGTGTCtt                <  1:4009677/90‑1 (MQ=255)

                                                                        aaaaTCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTGTTATATAATCTACAATTTTGTGTCttt               <  1:304073/91‑1 (MQ=255)

                                                                        aaaaTCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTGTTATATAATCTACAATTTTGTGTCttt               <  1:359488/91‑1 (MQ=255)

                                                                         aaaTCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTGTTATATAATCTACAATTTTGTGTCtttt              >  1:3641303/1‑91 (MQ=255)

                                                                         aaaTCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTGTTATATAATCTACAATTTTGTGTCtttt              <  1:1246102/91‑1 (MQ=255)

                                                                         aaaTCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTGTTATATAATCTACAATTTTGTGTCtttt              <  1:1246288/91‑1 (MQ=255)

                                                                           aTCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTGTTATATAATCTACAATTTTGTGTCTTTTat            >  1:3607708/1‑90 (MQ=255)

                                                                           aTCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTGTTATATAATCTACAATTTTGTGTCTTTTaa            <  1:143984/91‑1 (MQ=255)

                                                                            tCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTGTTATATAATCTACAATTTTGTGTCTTTTAAg           >  1:2025657/1‑91 (MQ=255)

                                                                            tCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTGTTATATAATCTACAATTTTGTGTCTTTTAAg           >  1:2025626/1‑91 (MQ=255)

                                                                             cATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTGTTATATAATCTACAATTTTGTGTCTTTTAAGt          <  1:2198293/91‑1 (MQ=255)

                                                                              aTTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTGTTATATAATCTACAATTTTGTGTCTTTTAAGTc         <  1:2378101/91‑1 (MQ=255)

                                                                                tCCTAATTTGGAATCATCTTACTGCTGTTTGCTGTTGATTTATATTCATGATTTTGTTATATAATCTACAATTTTGTTTCTTTTAAGTCtt       >  1:3434241/1‑91 (MQ=255)

                                                                                  cTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTGTTATATAATCTACAATTTTGTGTCTTTTAAGTCTTcc     >  1:2582358/1‑91 (MQ=255)

                                                                                  cTAATTTGGAATCAGCTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTGTTATATAATCGACAATTTTGTGTCTTTTAAGTCTTcc     <  1:3742117/91‑1 (MQ=255)

                                                                                   tAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTGTTATATAATCTACAATTTTGTGTCTTTTAAGTCTTCCg    <  1:1462518/91‑1 (MQ=255)

                                                                                   tAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTGTTATATAATCTACAATTTTGTGTCTTTTAAGTCTTCCg    <  1:1990169/91‑1 (MQ=255)

                                                                                     aTTTGGAATCATCTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTGTTATATAATCTACAATTTTGTGTCTTTTAAGTCTTCCGaa  >  1:1877810/1‑91 (MQ=255)

AGATAGGGCGAATCTGTAGTCAAAGAAAAAATCATTGAAAAAGAGGTAACAATGTCAAAAGAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTGTTATATAATCTACAATTTTGTGTCTTTTAAGTCTTCCGAA  >  USA300TCH1516_ALE/942133‑942308

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