BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	USA300TCH1516_ALE	942,230	C→T	8.2%	intergenic (‑201/+48)	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,230	0	C	T	8.2%	310.6 / 13.5	97	intergenic (‑201/+48)	ptlE/mnhG1	Neopentalenolactone D synthase/Na(+)/H(+) antiporter subunit G1
Reads supporting (aligned to +/- strand): ref base C (39/50); new base T (3/5); total (42/55)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.99e-01

Only 100 of 102 total aligned reads displayed.

GCGAATCTGTAGTCAAAGAAAAAATCATTGAAAAAGAGGTAACAATGTCAAAAGAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTGTTATATAATCTACAATTTTGTGTCTTTTAAGTCTTCCGAAATTTCATCGA  >  USA300TCH1516_ALE/942140‑942318

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:3919181/91‑1 (MQ=255)

            tCAAAGAAAAAATCATTGAAAAAGAGGTAACAATGTCAAAAGAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCTTACTGCtgttt                                                                              <  1:3919152/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:16554/1‑91 (MQ=255)

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

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

                       aTCATTGAAAAAGAGGTAACAATGTCAAAAGAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGAttt                                                                   >  1:3301526/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:3721588/1‑91 (MQ=255)

                                        aaCAATGTCAAAAGAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTGtt                                                  <  1:2573813/91‑1 (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:506879/1‑91 (MQ=255)

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

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

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

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

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

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

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

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

                                                  aaaGAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTGTTATATAATCTa                                        <  1:3445736/91‑1 (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:636692/91‑1 (MQ=255)

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

                                                        aaCAGCAGTAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTGTTATATAATCTACAAttt                                  >  1:1968583/1‑91 (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:4009488/91‑1 (MQ=255)

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

                                                                tAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTGTTATATAATCTACAATTTTGTGTCtt                          <  1:3966882/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:1246288/91‑1 (MQ=255)

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

                                                                  aaaTCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTGTTATATAATCTACAATTTTGTGTCtttt                        <  1:1246102/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)

                                                                                  ggAATCATCTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTGTTATATAATCTACAATTTTGTGTCTTTTAAGTCTTCCGAAAttt        <  1:1545232/91‑1 (MQ=255)

                                                                                    aaTCATCTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTGTTATATAATCTACAATTTTGTGTCTTTTAAGTCTTCCGAAATTTCa      <  1:2222437/91‑1 (MQ=255)

                                                                                    aaTCATCTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTGTTATATAATCTACAATTTTGTGTCTTTTAAGTCTTCCGAAATTTCa      <  1:382065/91‑1 (MQ=255)

                                                                                      tcatcTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTGTTATATAATCTACAATTTTGTGTCTTTTAAGTCTTCCGAAATTTCATc    <  1:2319523/91‑1 (MQ=255)

                                                                                      tcatcTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTGTTATATAATCTACAATTTTGTGTCTTTTAAGTCTTCCGAAATTTCATc    <  1:1260196/91‑1 (MQ=255)

                                                                                        atcTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTGTTATATAATCTACAATTTTGTGTCTTTTAAGTCTTCCGAAATTTCATCGa  >  1:2233783/1‑91 (MQ=255)

                                                                                        atcTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTGTTATATAATCTACAATTTTGTGTCTTTTAAGTCTTCCGAAATTTCATCGa  >  1:2233919/1‑91 (MQ=255)

GCGAATCTGTAGTCAAAGAAAAAATCATTGAAAAAGAGGTAACAATGTCAAAAGAAAACAGCAGTAAAATCATTCCTAATTTGGAATCATCTTACTGCTGTTTGTTGTTGATTTATATTCATGATTTTGTTATATAATCTACAATTTTGTGTCTTTTAAGTCTTCCGAAATTTCATCGA  >  USA300TCH1516_ALE/942140‑942318

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