BRESEQ :: Evidence

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	USA300TCH1516_ALE	512,643	0	A	G	20.6%	1170.3 / 271.7	494	intergenic (+711/‑6052)	recR/speA_2	Recombination protein RecR/Arginine decarboxylase
Reads supporting (aligned to +/- strand): ref base A (211/181); new base G (46/56); total (257/237)
Fisher's exact test for biased strand distribution p-value = 1.21e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

Only 100 of 548 total aligned reads displayed.

ATTTGTTGAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAAC  >  USA300TCH1516_ALE/512553‑512733

atttGTTGAGTCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTta                                                                                            <  2:5226104/91‑1 (MQ=255)

  ttgttgAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTtaat                                                                                          <  1:6096147/91‑1 (MQ=255)

   tgttgAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTtaata                                                                                         >  1:4373412/1‑91 (MQ=255)

    gttgAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTtaataa                                                                                        >  2:6481031/1‑91 (MQ=255)

    gttgAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTtaataa                                                                                        <  1:866422/91‑1 (MQ=255)

     ttgAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATaaa                                                                                       >  2:4925282/1‑91 (MQ=255)

     ttgAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATaaa                                                                                       >  2:6338412/1‑91 (MQ=255)

     ttgAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATaaa                                                                                       <  1:4233803/91‑1 (MQ=255)

      tgAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATaaaa                                                                                      >  1:2451907/1‑91 (MQ=255)

        aGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATg                                                                                    >  2:5894053/1‑91 (MQ=255)

        aGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATg                                                                                    >  1:2294172/1‑91 (MQ=255)

         ggCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGt                                                                                   >  1:763769/1‑91 (MQ=255)

          gCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGta                                                                                  >  1:1469099/1‑91 (MQ=255)

           caaacaaaACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGtat                                                                                 <  1:6229539/91‑1 (MQ=255)

             aacaaaacaaCTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTAtaa                                                                               <  2:4860410/91‑1 (MQ=255)

             aacaaaacaaCTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTAtaa                                                                               <  1:5947792/91‑1 (MQ=255)

               caaaacaaCTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTAtaatt                                                                             >  2:5145853/1‑91 (MQ=255)

                aaaacaaCTAAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTTACTTAAATGTTAATAAAATGTAtaatta                                                                            <  1:6579840/91‑1 (MQ=255)

                 aaacaaCTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTAtaattaa                                                                           <  2:2297524/91‑1 (MQ=255)

                  aacaacTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTAtaattaat                                                                          <  1:1082516/91‑1 (MQ=255)

                   acaacTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATaattaatt                                                                         >  1:3348990/1‑91 (MQ=255)

                   acaacTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATaattaatt                                                                         >  2:2705246/1‑91 (MQ=255)

                   acaacTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATaattaatt                                                                         >  2:6414984/1‑91 (MQ=255)

                   acaacTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATaattaatt                                                                         >  2:2808507/1‑91 (MQ=255)

                       ctcaactTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTg                                                                     >  1:235887/1‑91 (MQ=255)

                       ctcaactTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTg                                                                     >  2:5604260/1‑91 (MQ=255)

                       ctcaactTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTg                                                                     >  1:138553/1‑91 (MQ=255)

                        tcaactTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGt                                                                    <  1:2205591/91‑1 (MQ=255)

                            ctTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGt                                                                >  2:409037/1‑91 (MQ=255)

                              tAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTaa                                                              <  1:2031868/91‑1 (MQ=255)

                               aaGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAg                                                             <  1:2231632/91‑1 (MQ=255)

                                 gAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGaa                                                           >  2:3077462/1‑91 (MQ=255)

                                   aaTAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGaaaa                                                         >  1:4093428/1‑91 (MQ=255)

                                   aaTAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGaaaa                                                         <  1:1579409/91‑1 (MQ=255)

                                    aTAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGaaaaa                                                        >  1:5805163/1‑91 (MQ=255)

                                     tAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAAt                                                       <  2:624402/91‑1 (MQ=255)

                                      aaCTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGa                                                            <  2:4222273/85‑1 (MQ=255)

                                      aaCTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATg                                                      <  2:1518516/91‑1 (MQ=255)

                                      aaCTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATg                                                      >  2:2420413/1‑91 (MQ=255)

                                       aCTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGa                                                     >  1:3454187/1‑91 (MQ=255)

                                          tGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTGAATAAAATGTATAATTAATTCTTGTCGGAAAGAAAAATGAACa                                                  >  1:4305814/1‑91 (MQ=255)

                                           gAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACAt                                                 >  2:2528172/1‑91 (MQ=255)

                                           gAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACAt                                                 >  2:3390919/1‑91 (MQ=255)

                                           gAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACAt                                                 <  1:4876095/91‑1 (MQ=255)

                                            aaTTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACAtt                                                <  1:518326/91‑1 (MQ=255)

                                            aaTTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACAtt                                                <  2:2959170/91‑1 (MQ=255)

                                             aTTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTg                                               <  2:2490928/91‑1 (MQ=255)

                                                  tAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAAc                                          >  1:4929019/1‑91 (MQ=255)

                                                  tAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAAc                                          >  1:1023318/1‑91 (MQ=255)

                                                     cGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGa                                       <  1:5115988/91‑1 (MQ=255)

                                                      gAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGaa                                      >  1:3629535/1‑91 (MQ=255)

                                                      gAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGaa                                      <  1:1869172/91‑1 (MQ=255)

                                                            taattTTAAAAAGTTATTGACTTAAGTGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACaa                                <  2:6107911/91‑1 (MQ=255)

                                                            taattTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACaa                                >  2:3682482/1‑91 (MQ=255)

                                                                 aaaaaaaGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGt                           <  1:292186/89‑1 (MQ=255)

                                                                  tAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTc                          <  1:477811/91‑1 (MQ=255)

                                                                   aaaaaGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCa                         >  2:3129430/1‑91 (MQ=255)

                                                                   aaaaaGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCa                         <  2:914386/91‑1 (MQ=255)

                                                                   aaaaaGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCa                         <  1:5558555/91‑1 (MQ=255)

                                                                   aaaaaGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCa                         >  1:2065880/1‑91 (MQ=255)

                                                                    aaaaGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCaa                        <  2:3808426/91‑1 (MQ=255)

                                                                     aaaGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAAc                       >  2:5341314/1‑91 (MQ=255)

                                                                       aGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGt                     >  1:2641169/1‑91 (MQ=255)

                                                                       aGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGt                     >  2:5508187/1‑91 (MQ=255)

                                                                       aGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGt                     >  2:1247222/1‑91 (MQ=255)

                                                                        gTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGtt                    <  2:1833184/91‑1 (MQ=255)

                                                                        gTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGtt                    >  1:1731478/1‑91 (MQ=255)

                                                                        gTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGtt                    <  2:2190140/91‑1 (MQ=255)

                                                                        gTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGtt                    >  1:1827883/1‑91 (MQ=255)

                                                                        gTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGtt                    <  2:4540994/91‑1 (MQ=255)

                                                                         ttATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTa                   <  1:4483440/91‑1 (MQ=255)

                                                                         ttATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTa                   <  2:4385554/91‑1 (MQ=255)

                                                                           aTTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAAt                 <  1:3934084/91‑1 (MQ=255)

                                                                           aTTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAAt                 >  1:1258369/1‑91 (MQ=255)

                                                                           aTTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAAt                 >  2:143814/1‑91 (MQ=255)

                                                                           aTTGACTTAAATGTTGATAAAATGTATAATTAATTCTAGTAGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAAt                 <  2:1903713/91‑1 (MQ=255)

                                                                           aTTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAAt                 >  1:2848561/1‑91 (MQ=255)

                                                                           aTTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAAt                 <  1:5271663/91‑1 (MQ=255)

                                                                            ttGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAAtt                <  1:3838267/91‑1 (MQ=255)

                                                                            ttGACTTAAATGTTAATACAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAAtt                <  1:6488170/91‑1 (MQ=255)

                                                                              gACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTcc              <  1:4629030/91‑1 (MQ=255)

                                                                              gACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTcc              <  1:5642456/91‑1 (MQ=255)

                                                                              gACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTcc              <  2:3231995/91‑1 (MQ=255)

                                                                               aCTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCa             <  1:3154841/91‑1 (MQ=255)

                                                                                cTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCaa            <  1:6315617/91‑1 (MQ=255)

                                                                                  tAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCaaaa          >  2:1147384/1‑91 (MQ=255)

                                                                                  tAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCaaaa          >  1:900473/1‑91 (MQ=255)

                                                                                     aTGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAAc       >  2:2045510/1‑91 (MQ=255)

                                                                                     aTGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAAc       <  1:3904170/91‑1 (MQ=255)

                                                                                     aTGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACAGTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAAc       >  1:6408817/1‑91 (MQ=255)

                                                                                      tGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACACTATGTCAACGTTAATTCCAAAAAACg      >  2:4656802/1‑91 (MQ=255)

                                                                                      tGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACg      <  2:4021996/91‑1 (MQ=255)

                                                                                      tGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACg      <  2:680651/91‑1 (MQ=255)

                                                                                        ttAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTa    >  2:3442553/1‑91 (MQ=255)

                                                                                        ttAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTa    <  2:1332042/91‑1 (MQ=255)

                                                                                        ttAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTa    >  2:3655067/1‑91 (MQ=255)

                                                                                         taataaAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTaa   <  1:3071140/91‑1 (MQ=255)

                                                                                          aataaAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAAc  >  1:4598897/1‑91 (MQ=255)

                                                                                          aataaAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAAc  >  2:1663371/1‑91 (MQ=255)

                                                                                          aataaAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAAAATAAAAAAAAAAAAACGTAAc  <  2:575440/91‑1 (MQ=25)

ATTTGTTGAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAAC  >  USA300TCH1516_ALE/512553‑512733

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