BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	USA300TCH1516_ALE	1,439,789	C→A	5.4%	intergenic (‑146/+51)	USA300HOU_RS07135 ← / ← cspA_2	hypothetical protein/Cold shock protein CspA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	USA300TCH1516_ALE	1,439,789	0	C	A	5.4%	427.3 / 11.0	125	intergenic (‑146/+51)	USA300HOU_RS07135/cspA_2	hypothetical protein/Cold shock protein CspA
Reads supporting (aligned to +/- strand): ref base C (61/53); new base A (3/3); total (67/58)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.16e-01

Only 100 of 130 total aligned reads displayed.

GCAAACAATGTTGGTGATATAAAAAAAGAGTTCTCAATTTGGTTGAGAAACTCATCGCGTTAAGTATTTGTAATAAGCATATAAAAACACCTTACTTCTTGGTAAGGTGTTAGTTAAAATCAATGATTCAAATCTAAGAATTATAGTTTAACAACGTTTGCAGCTTGTGGACCGCGGTC  >  USA300TCH1516_ALE/1439700‑1439878

 cAAACAATGTTGGTGATATAAAAAAAGAGTTCTCAATTTGGTTGAGAAACTCATCGCGTTAAGTATTTGTAATAAGCATATAAAAACACCt                                                                                         >  1:3262953/1‑91 (MQ=255)

 cAAACAATGTTGGTGATATAAAAAAAGAGTTCTCAATTTGGTTGAGAAACTCATCGCGTTAAGTATTTGTAATAAGCATATAAAAACACCt                                                                                         <  1:1949404/91‑1 (MQ=255)

 cAAACAATGTTGGTGATATAAAAAAAGAGTTCTCAATTTGGTTGAGAAACTCATCGCGTTAAGTATTTGTAATAAGCATATAAAAACACCt                                                                                         <  1:1639825/91‑1 (MQ=255)

 cAAACAATGTTGGTGATATAAAAAAAGAGTTCTCAATTTGGTTGAGAAACTCATCGCGTTAAGTATTTGTAATAAGCATATAAAAACACCt                                                                                         <  1:2600886/91‑1 (MQ=255)

 cAAACAATGTTGGTGATATAAAAAAAGAGTTCTCAATTTGGTTGAGAAACTCATCGCGTTAAGTATTTGTAATAAGCATATAAAAACACCt                                                                                         >  1:3263067/1‑91 (MQ=255)

 cAAACAATGTTGGTGATATAAAAAAAGAGTTCTCAATTTGGTTGAGAAACTCATCGCGTTAAGTATTTGTAATAAGCATATAAAAACACCt                                                                                         <  1:1314419/91‑1 (MQ=255)

 cAAACAATGTTGGTGATATAAAAAAAGAGTTCTCAATTTGGTTGAGAAACTCATCGCGTTAAGTATTTGTAATAAGCATATAAAAACACCt                                                                                         <  1:1314406/91‑1 (MQ=255)

     cAATGTTGGTGATATAAAAAAAGAGTTCTCAATTTGGTTGAGAAACTCATCGCGTTAAGTATTTGTAATAAGCATATAAAAACACCTTAct                                                                                     <  1:2872597/91‑1 (MQ=255)

      aaTGTTGGTGATATAAAAAAAGAGTTCTCAATTTGGTTGAGAAACTCATCGCGTTAAGTATTTGTAATAAGCATATAAAAACACCTTActt                                                                                    <  1:3379169/91‑1 (MQ=255)

      aaTGTTGGTGATATAAAAAAAGAGTTCTCAATTTGGTTGAGAAACTCATCGCGTTAAGTATTTGTAATAAGCATATAAAAACACCTTActt                                                                                    <  1:3378562/91‑1 (MQ=255)

      aaTGTTGGTGATATAAAAAAAGAGTTCTCAATTTGGTTGAGAAACTCATCGCGTTAAGTATTTGTAATAAGCATATAAAAACACCTTAcca                                                                                    <  1:832373/91‑3 (MQ=255)

       aTGTTGGTGATATAAAAAAAGAGTTCTCAATTTGGTTGAGAAACTCATCGCGTTAAGTATTTGTAATAAGCATATAAAAACACCTTActtc                                                                                   >  1:3891244/1‑91 (MQ=255)

        tGTTGGTGATATAAAAAAAGAGTTCTCAATTTGGTTGAGAAACTCATCGCGTTAAGTATTTGTAATAAGCATATAAAAACACCTTActtct                                                                                  <  1:2313482/91‑1 (MQ=255)

        tGTTGGTGATATAAAAAAAGAGTTCTCAATTTGGTTGAGAAACTCATCGCGTTAAGTATTTGTAATAAGCATATAAAAACACCTTActtct                                                                                  <  1:814687/91‑1 (MQ=255)

        tGTTGGTGATATAAAAAAAGAGTTCTCAATTTGGTTGAGAAACTCATCGCGTTAAGTATTTGTAATAAGCATATAAAAACACCTTActtct                                                                                  >  1:2154612/1‑91 (MQ=255)

         gTTGGTGATATAAAAAAAGAGTTCTCAATTTGGTTGAGAAACTCATCGCGTTAAGTATTTGTAATAAGCATATAAAAACACCTTActtctt                                                                                 >  1:2670571/1‑91 (MQ=255)

         gTTGGTGATATAAAAAAAGAGTTCTCAATTTGGTTGAGAAACTCATCGCGTTAAGTATTTGTAATAAGCATATAAAAACACCTTActtctt                                                                                 >  1:2038443/1‑91 (MQ=255)

            ggTGATATAAAAAAAGAGTTCTCAATTTGGTTGAGAAACTCATCGCGTTAAGTGTTTGTAATAAGCATATAAAAACACCTTACTTCTTGGt                                                                              <  1:3515288/91‑1 (MQ=255)

            ggTGATATAAAAAAAGAGTTCTCAATTTGGTTGAGAAACTCATCGCGTTAAGTATTTGTAATAAGCATATAAAAACACCTTACTTCTTGGt                                                                              <  1:1040393/91‑1 (MQ=255)

            ggTGATATAAAAAAAGAGTTCTCAATTTGGTTGAGAAACTCATCGCGTTAAGTATTTGTAATAAGCATATAAAAACACCTTACTTCTTGGt                                                                              >  1:2499528/1‑91 (MQ=255)

            ggTGATATAAAAAAAGAGTTCTCAATTTGGTTGAGAAACTCATCGCGTTAAGTATTTGTAATAAGCATATAAAAACACCTTACTTCTTGGt                                                                              >  1:2228787/1‑91 (MQ=255)

            ggTGATATAAAAAAAGAGTTCTCAATTTGGTTGAGAAACTCATCGCGTTAAGTATTTGTAATAAGCATATAAAAACACCTTACTTCTTGGt                                                                              <  1:1040754/91‑1 (MQ=255)

             gTGATATAAAAAAAGAGTTCTCAATTTGGTTGAGAAACTCATCGCGTTAAGTATTTGTAATAAGCATATAAAAACACCTTACTTCTTGGTa                                                                             >  1:271076/1‑91 (MQ=255)

                atatAAAAAAAGAGTTCTCAATTTGGTTGAGAAACTCATCGCGTTAAGTATTTGTAATAAGCATATAAAAACAGCTTACTTCTTGGTAAgg                                                                          <  1:80482/91‑1 (MQ=255)

                atatAAAAAAAGAGTTCTCAATTTGGTTGAGAAACTCATCGCGTTAAGTATTTGTAATAAGCATATAAAAACAGCTTACTTCTTGGTAAgg                                                                          <  1:80283/91‑1 (MQ=255)

                atatAAAAAAAGAGTTCTCAATTTGGTTGAGAAACTCATCGCGTTAAGTATTTGTAATAAGCATATAAAAACACCTTACTTCTTGGTAAgg                                                                          <  1:578040/91‑1 (MQ=255)

                atatAAAAAAAGAGTTCTCAATTTGGTTGAGAAACTCATCGCGTTAAGTATTTGTAATAAGCATATAAAAACACCTTACTTCTTGGTAAgc                                                                          >  1:2370849/1‑90 (MQ=255)

                   taAAAAAAGAGTTCTCAATTTGGTTGAGAAACTCATCGCGTTAAGTATTTGTAATAAGCATATAAAAACACCTTACTTCTTGGTAAGgtgt                                                                       >  1:2427616/1‑91 (MQ=255)

                     aaaaaaGAGTTCTCAATTTGGTTGAGAAACTCATCGCGTTAAGTATTTGTAATAAGCATATAAAAACACCTTACTTCTTGGTAAGGTgtta                                                                     >  1:1057461/1‑91 (MQ=255)

                       aaaaGAGTTCTCAATTTGGTTGAGAAACTCATCGCGTTAAGTATTTGTAATAAGCATATAAAAACACCTTACTTCTTGGTATGGTgttagt                                                                   <  1:59632/91‑1 (MQ=255)

                       aaaaGAGTTCTCAATTTGGTTGAGAAACTCATCGCGTTAAGTATTTGTAATAAGCATATAAAAACACCTTACTTCTTGGTAAGGTgttagt                                                                   >  1:3385028/1‑91 (MQ=255)

                        aaaGAGTTCTCAATTTGGTTGAGAAACTCATCGCGTTAAGTATTTGTAATAAGCATATAAAAACACCTTACTTCTTGGTAAGGTgttagtt                                                                  >  1:1835073/1‑91 (MQ=255)

                            agTTCTCAATTTGGTTGAGAAACTCATCGCGTTAAGTATTTGTAATAAGCATATAAAAACACCTTACTTCTTGGTAAGGTGTTAGTTaaaa                                                              >  1:3428329/1‑91 (MQ=255)

                            agTTCTCAATTTGGTTGAGAAACTCATCGCGTTAAGTATTTGTAATAAGCATATAAAAACACCTTACTTCTTGGTAAGATGTTAGTTaaaa                                                              <  1:152941/91‑1 (MQ=255)

                             gTTCTCAATTTGGTTGAGAAACTCATCGCGTTAAGTATTTGTAATAAGCATATAAAAACACCTTACTTCTTGGTAAGGTGTTAGTTAAAAt                                                             >  1:632380/1‑91 (MQ=255)

                             gTTCTCAATTTGGTTGAGAAACTCATCGCGTTAAGTATTTGTAATAAGCATATAAAAACACCTTACTTCTTGGTAAGGTGTTAGTTAAAAt                                                             >  1:1871397/1‑91 (MQ=255)

                             gTTCTCAATTTGGTTGAGAAACTCATCGCGTTAAGTATTTGTAATAAGCATATAAAAACACCTTACTTCTTGGTAAGCTGTTAGTTAAAAt                                                             >  1:911853/1‑91 (MQ=255)

                             gTTCTCAATTTGGTTGAGAAACTCATCGCGTTAAGTATTTGTAATAAGCATATAAAAACACCTTACTTCTTGGTAAGCTGTTAGTTAAAAt                                                             >  1:1404141/1‑91 (MQ=255)

                              ttCTCAATTTGGTTGAGAAACTCATCGCGTTAAGTATTTGTAATAAGCATATAAAAACACCTTACTTCTTGGTAAGGTGTTAGTTAAAATc                                                            <  1:126301/91‑1 (MQ=255)

                              ttCTCAATTTGGTTGAGAAACTCATCGCGTTAAGTATTTGTAATAAGCATATAAAAACACCTTACTTCTTGGTAAGGTGTTAGTTAAAATc                                                            <  1:4113903/91‑1 (MQ=255)

                                    aTTTGGTTGAGCAACTCATCGCGTTAAGTATTTGTAATAAGCATATAAAAACACCTTACTTCTTGGTAAGGTGTTAGTTAAAATCAATGAt                                                      <  1:3768696/91‑1 (MQ=255)

                                     tttGGTTGAGAAACTCATCGCGTTAAGTATTTGTAATAAGCATATAAAAACACCTTACTTCTTGGTAAGGTGTTAGTTAAAATCAATGAtt                                                     >  1:1813263/1‑91 (MQ=255)

                                      ttGGTTGAGAAACTCATCGCGTTAAGTATTTGTAATAAGCATATAAAAACACCTTACTTCTTGGTAAGGTGTTAGTTAAAATCAATGATTc                                                    <  1:908697/91‑1 (MQ=255)

                                       tGGTTGAGAAACTCATCGCGTTAAGTATTTGTAATAAGCATATAAAAACACCTTACTTCTTGGTAAGGTGTTAGTTAAAATCAATGATTCa                                                   >  1:2544795/1‑91 (MQ=255)

                                         gTTGAGAAACTCATCGCGTTAAGTATTTGTAATAAGCATATAAAAACACCTTACTTCTTGGTAAGCTGTTAGTTAAAATCAATGATTCaaa                                                 >  1:528399/1‑91 (MQ=255)

                                          ttGGGAAACTCATCGCGTTAAGTATTTGTAATAAGCATATAAAAACACCTTACTTCTTGGTAAGATGTTAGTTAAAATCAATGATTCAAAt                                                <  1:3737322/91‑1 (MQ=255)

                                          ttGAGAAACTCATCGCGTTAAGTATTTGTAATAAGCATATAAAAACAGCTTACTTCTTGGTAAGGTGTTAGTTAAAATCAATGATTCAAAt                                                >  1:3389622/1‑91 (MQ=255)

                                          ttGAGAAACTCATCGCGTTAAGTATTTGTAATAAGCATATAAAAACACCTTACTTCTTGGTAAGGTGTTAGTTAAAATCAATGATTCAAAt                                                >  1:2526047/1‑91 (MQ=255)

                                          ttGAGAAACTCATCGCGTTAAGTATTTGTAATAAGCATATAAAAACACCTTACTTCTTGGTAAGATGTTAGTTAAAATCAATGATTCAAAt                                                <  1:2563974/91‑1 (MQ=255)

                                           tGAGAAACTCGTCGCGTTAAGTATTTGTAATAAGCATATAAAAACACCTTACTTCTTGGTAAGCTGTTAGTTAAAATCAATGATTCAAATc                                               >  1:347401/1‑91 (MQ=255)

                                            gagaAACTCATCGCGTTAAGTATTTGTAATAAGCATATAAAAACACCTTACTTCTTGGTAAGGTGTTAGTTAAAATCAATGATTCAAATCt                                              <  1:1700908/91‑1 (MQ=255)

                                            gagaAACTCATCGCGTTAAGTATTTGTAATAAGCATATAAAAACACCTTACTTCTTGGTAAGGTGTTAGTTAAAATCAATGATTCAAATCt                                              >  1:638646/1‑91 (MQ=255)

                                            gagaAACTCATCGCGTTAAGTATTTGTAATAAGCATATAAAAACACCTTACTTCTTGGTAAGGTGTTAGTTAAAATCAATGATTCAAATCt                                              >  1:2102283/1‑91 (MQ=255)

                                                  cTCATCGCGTTAAGTATTTGTAATAAGCATATAAAAACACCTTACTTCTTGGTAAGGTGTTAGTTAAAATCAATGATTCAAATCTAAGAAt                                        <  1:309356/91‑1 (MQ=255)

                                                  cTCATCGCGTTAAGTATTTGTAATAAGCATATAAAAACACCTTACTTCTTGGTAAGGTGTTAGTTAAAATCAATGATTCAAATCTAAGAAt                                        <  1:309608/91‑1 (MQ=255)

                                                   tCATCGCGTTAAGTATTTGTAATAAGCATATAAAAACACCTTACTTCTTGGTAAGGTGTTAGTTAAAATCAATGATTCAAATCTAAGAAtt                                       <  1:1667917/91‑1 (MQ=255)

                                                   tCATCGCGTTAAGTATTTGTAATAAGCATATAAAAACACCTTACTTCTTGGTAAGGTGTTAGTTAAAATCAATGATTCAAATCTAAGAAtt                                       <  1:2189030/91‑1 (MQ=255)

                                                   tCATCGCGTTAAGTATTTGTAATAAGCATATAAAAACACCTTACTTCTTGGTAAGCTGTTAGTTAAAATCAATGATTCAAATCTAAGAAtt                                       <  1:758486/91‑1 (MQ=255)

                                                       cgcgTTAAGTATTTGTAATAAGCATATAAAATCACCTTACTTCTTGGTAAGGTGTTAGTTAAAATCAATGATTCAAATCTAAGAATTATAg                                   >  1:470844/1‑91 (MQ=255)

                                                       cgcgTTAAGTATTTGTAATAAGCATATAAAAACACCTTACTTCTTGGTAAGGTGTTAGTTAAAATCAATGATTCAAATCTAAGAATTATAg                                   >  1:470826/1‑91 (MQ=255)

                                                        gcgTTAAGTATTTGTAATAAGCATATAAAAACACCTTACTTCTTGGTAAGGTGTTAGTTAAAATCAATGATTCAAATCTAAGAATTATAGt                                  >  1:3130098/1‑91 (MQ=255)

                                                        gcgTTAAGTATTTGTAATAAGCATATAAAAACACCTTACTTCTTGGTAAGCTGTTAGTTAAAATCAATGATTCAAATCTAAGAATTATAGt                                  >  1:1825621/1‑91 (MQ=255)

                                                         cgTTAAGTATTTGTAATAAGCATATAAAAACACCTTACTTCTTGGTAAGGTGTTAGTTAAAATCAATGATTCAAATCTAAGAATTATAGtt                                 <  1:1600834/91‑1 (MQ=255)

                                                             aaGTATTTGTAATAAGCATATAAAAACACCTTACTTCTTGGTAAGGTGTTAGTTAAAATCAATGATTCAAATCTAAGAATTATAGTTTaac                             >  1:160580/1‑91 (MQ=255)

                                                              aGTATTTGTAATAAGCATATAAAAACACCTTACTTCTTGGTAAGGTGTTAGTTAAAATCAATGATTCAAATCTAAGAATTATAGTTTaaca                            >  1:4124076/1‑91 (MQ=255)

                                                               gTATTTGTAATAAGCATATAAAAACAGCTTACTTCTTGGTAAGGTGTTAGTTAAAATCAATGATTCAAATCTAAGAATTATAGTTTaacaa                           >  1:2167624/1‑91 (MQ=255)

                                                                tATTTGTAATAAGCATATAAAAACACCTTACTTCTTGGTAAGCTGTTAGTTAAAATCAATGATTCAAATCTAAGAATTATAGTTTaacaac                          >  1:2974861/1‑91 (MQ=255)

                                                                  tttGTAATAAGCATATAAAAACAGCTTACTTCTTGGTAAGGTGTTAGTTAAAATCAATGATTCAAATCTAAGAATTATAGTTTAACAACGt                        >  1:1758295/1‑91 (MQ=255)

                                                                  tttGTAATAAGCATATAAAAACACCTTACTTCTTGGTAAGGTGTTAGTTAAAATCAATGATTCAAATCTAAGAATTATAGTTTAACAACGt                        <  1:1533566/91‑1 (MQ=255)

                                                                  tttGTAATAAGCATATAAAAACAACTTACTTCTTGGTAAGGTGTTAGTTAAAATCAATGATTCAAATCTAAGAATTATAGTTTAACAACGt                        >  1:2282722/1‑91 (MQ=255)

                                                                    tGTAATAAGCATATAAAAACACCTTACTTCTTGGTAAGGTGTTAGTTAAAATCAATGATTCAAATCTAAGAATTATAGTTTAACAACGttt                      <  1:360315/91‑1 (MQ=255)

                                                                    tGTAATAAGCATATAAAAACACCTTACTTCTTGGTAAGGTGTTAGTTAAAATCAATGATTCAAATCTAAGAATTATAGTTTAACAACGttt                      >  1:2560540/1‑91 (MQ=255)

                                                                    tGTAATAAGCATATAAAAACACCTTACTTCTTGGTAAGCTGTTAGTTAAAATCAATGATTCAAATCTAAGAATTATAGTTTAACAACGttt                      >  1:2585370/1‑91 (MQ=255)

                                                                     gTAATAAGCATATAAAAACACCTTACTTCTTGGTAAGGTGTTAGTTAAAATCAATGATTCAAATCTAAGAATTATAGTTTAACAACGTTTg                     >  1:3639905/1‑91 (MQ=255)

                                                                     gTAATAAGCATATAAAAACACCTTACTTCTTGGTAAGGTGTTAGTTAAAATCAATGATTCAAATCTAAGAATTATAGTTTAACAACGTTTg                     <  1:3040414/91‑1 (MQ=255)

                                                                     gTAATAAGCATATAAAAACACCTTACTTCTTGGTAAGGTGTTAGTTAAAATCAATGATTCAAATCTAAGAATTATAGTTTAACAACGTTTg                     <  1:3119713/91‑1 (MQ=255)

                                                                      taataaGCATATAAAAACACCTTACTTCTTGGTAAGGTGTTAGTTAAAATCAATGATTCAAATCTAAGAATTATAGTTTAACAACGTTTGc                    >  1:3673949/1‑91 (MQ=255)

                                                                      taataaGCATATAAAAACACCTTACTTCTTGGTAAGGTGTTAGTTAAAATCAATGATTCAAATCTAAGAATTATAGTTTAACAACGTTTGc                    >  1:2510088/1‑91 (MQ=255)

                                                                        ataaGCATATAAAAACACCTTACTTCTTGGTAAGGTGTTAGTTAAAATCAATGATTCAAATCTAAGAATTATAGTTTAACAACGTTTGCAg                  >  1:1635415/1‑91 (MQ=255)

                                                                        ataaGCATATAAAAACACCTTACTTCTTGGTAAGGTGTTAGTTAAAATCAATGATTCAAATCTAAGAATTATAGTTTAACAACGTTTGCAg                  <  1:2647487/91‑1 (MQ=255)

                                                                        ataaGCATATAAAAACACCTTACTTCTTGGTAAGGTGTTAGTTAAAATCAATGATTCAAATCTAAGAATTATAGTTTAACAACGTTTGCAg                  <  1:64416/91‑1 (MQ=255)

                                                                         taaGCATATAAAAACACCTTACTTCTTGGTAAGCTGTTAGTTAAAATCAATGATTCAAATCTAAGAATTATAGTTTAACAACGTTTGCAGc                 <  1:1512721/91‑1 (MQ=255)

                                                                          aaGCCTATAAAAACACCTTACTTCTTGGTAAGGTGTTAGTTAAAATCAATGATTCAAATCTAAGAATTATAGTTTAACAACGTTTGCAGCt                <  1:3471764/91‑1 (MQ=255)

                                                                          aaGCATATAAAAACACCTTACTTCTTGGTAAGGTGTTAGTTAAAATCAATGATTCAAATCTAAGAATTATAGTTTAACAACGTTTGCAGCt                <  1:4096387/91‑1 (MQ=255)

                                                                           aGCATATAAAAACACCTTACTTCTTGGTAAGCTGTTAGTTAAAATCAATGATTCAAATCTAAGAATTATAGTTTAACAACGTTTGCAGCtt               >  1:3415834/1‑91 (MQ=255)

                                                                             cATATAAAAACACCTTACTTCTTGGTAAGGTGTTAGTTAAAATCAATGATTCAAATCTAAGAATTATAGTTTAACAACGTTTGCAGCTtgt             >  1:4009997/1‑91 (MQ=255)

                                                                              atatAAAAACACCTTACTTCTTGGTAAGGTGTTAGTTAAAATCAATGATTCAAATCTAAGAATTATAGTTTAACAACGTTTGCAGCTtgtg            >  1:36331/1‑91 (MQ=255)

                                                                              atatAAAAACACCTTACTTCTTGGTAAGGTGTTAGTTAAAATCAATGATTCAAATCTAAGAATTATAGTTTAACAACGTTTGCAGCTtgtg            >  1:3883946/1‑91 (MQ=255)

                                                                               tataAAAACACCTTACTTCTTGGTAAGGTGTTAGTTAAAATCAATGATTCAAATCTAAGAATTATAGTTTAACAACGTTTGCAGCTTGTgg           >  1:3883024/1‑91 (MQ=255)

                                                                                ataAAAACACCTTACTTCTTGGTAAGCTGTTAGTTAAAATCAATGATTCAAATCTAAGAATTATAGTTTAACAACGTTTGCAGCTTGTGGa          >  1:1003425/1‑91 (MQ=255)

                                                                                 taAAAACAACTTACTTCTTGGTAAGGTGTTAGTTAAAATCAATGATTCAAATCTAAGAATTATAGTTTAACAACGTTTGCAGCTTGTGGAc         <  1:3404530/91‑1 (MQ=255)

                                                                                 taAAAACAACTTACTTCTTGGTAAGGTGTTAGTTAAAATCAATGATTCAAATCTAAGAATTATAGTTTAACAACGTTTGCAGCTTGTGGAc         <  1:2618922/91‑1 (MQ=255)

                                                                                  aaaaaCACCTTACTTCTTGGTAAGGTGTTAGTTAAAATCAATGATTCAAATCTAAGAATTATAGTTTAACAACGTTTGCAGCTTGTGGAcc        >  1:1341001/1‑91 (MQ=255)

                                                                                  aaaaaCACCTTACTTCTTGGTAAGGTGTTAGTTAAAATCAATGATTCAAATCTAAGAATTATAGTTTAACAACGTTTGCAGCTTGTGGAcc        >  1:1341092/1‑91 (MQ=255)

                                                                                   aaaaCACCTTACTTCTTGGTAAGGTGTTAGTTAAAATCAATGATTCAAATCTAAGAATTATAGTTTAACAACGTTTGCAGCTTGTGGACcg       >  1:2211438/1‑91 (MQ=255)

                                                                                   aaaaCACCTTACTTCTTGGTAAGGTGTTAGTTAAAATCAATGATTCAAATCTAAGAATTATAGTTTAACAACGTTTGCAGCTTGTGGACcg       >  1:3790503/1‑91 (MQ=255)

                                                                                   aaaaCACCTTACTTCTTGGTAAGCTGTTAGTTAAAATCAATGATTCAAATCTAAGAATTATAGTTTAACAACGTTTGCAGCTTGTGGACcg       >  1:1808088/1‑91 (MQ=255)

                                                                                     aaCACCTTACTTCTTGGTAAGGTGTTAGTTAAAATCAATGATTCAAATCTAAGAATTATAGTTTAACAACGTTTGCAGATTGTGGACcgcg     <  1:3189154/91‑1 (MQ=255)

                                                                                       cacCTTACTTCTTGGTAAGGTGTTAGTTAAAATCAATGATTCAAATCTAAGAATTATAGTTTAACAACGTTTGCAGCTTGTGGACCGCGGt   <  1:4107656/91‑1 (MQ=255)

                                                                                        acCTTACTTCTTGGTAAGGTGTTAGTTAAAATCAATGATTCAAATCTAAGAATTATAGTTTAACAACGTTTGCAGCTTGTGGACCGCGGTc  >  1:1515719/1‑91 (MQ=255)

GCAAACAATGTTGGTGATATAAAAAAAGAGTTCTCAATTTGGTTGAGAAACTCATCGCGTTAAGTATTTGTAATAAGCATATAAAAACACCTTACTTCTTGGTAAGGTGTTAGTTAAAATCAATGATTCAAATCTAAGAATTATAGTTTAACAACGTTTGCAGCTTGTGGACCGCGGTC  >  USA300TCH1516_ALE/1439700‑1439878

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