BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	USA300TCH1516_ALE	41,166	G→A	9.7%	intergenic (‑33/‑67)	pbp ← / → mecR1	Beta‑lactam‑inducible penicillin‑binding protein/Methicillin resistance mecR1 protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	USA300TCH1516_ALE	41,166	0	G	A	9.7%	433.5 / 28.4	145	intergenic (‑33/‑67)	pbp/mecR1	Beta‑lactam‑inducible penicillin‑binding protein/Methicillin resistance mecR1 protein
Reads supporting (aligned to +/- strand): ref base G (66/65); new base A (5/9); total (71/74)
Fisher's exact test for biased strand distribution p-value = 4.01e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.17e-01

Only 100 of 154 total aligned reads displayed.

ACCAAACCCGACAACTACAACTATTAAAATAAGTGGAACAATTTTTATCTTTTTCATCAATATACTCCTTATATAAGACTACATTTGTAGTATATTACAAATGTAGTATTTATGTCAAAATAATGTTATAATTTTTGTGATATGGAGGTGTAGAAGGTGTTATCATCTTTTTTAATGTTA  >  USA300TCH1516_ALE/41077‑41256

  cAAACCCGACAACTACAACTATTAAAATAAGTGGAACAATTTTTATCTTTTTCATCAATATACTCCTTATATAAGACTACATTTGTAGtat                                                                                         >  1:2242748/1‑91 (MQ=255)

   aaaCCCGACAACTACAACTATTAAAATAAGTGGAACAATTTTTATCTTTTTCATCAATATACTCCTTATATAAGACTACATTTGTAGtata                                                                                        <  1:3615736/91‑1 (MQ=255)

    aaCCCGACAACTACAACTATTAAAATAAGTGGAACAATTTTTATCTTTTTCATCAATATACTCCTTATATAAGACTACATTTGTAAtatac                                                                                       <  1:2660750/91‑2 (MQ=255)

    aaCCCGACAACTACAACTATTAAAATAAGTGGAACAATTTTTATCTTTTTCATCAATATACTCCTTATATAAGACTACATTTGTAAtatac                                                                                       <  1:2660749/91‑2 (MQ=255)

      cccGACAACTACAACTATTAAAATAAGTGTATCAATTTTTATCTCTTTCATCATTATTCTCCTTATATAAGACTACATTTGTAATATACTa                                                                                     <  1:3317151/91‑1 (MQ=255)

      cccGACAACTACAACTATTAAAATAAGTGGAACAATTTTTATCTTTTTCATCAATATACTCCTTATATAAGACTACATTTGTAGTATATTa                                                                                     <  1:979200/91‑1 (MQ=255)

      cccGACAACTACAACTATTAAAATAAGTGGAACAATTTTTATCTTTTTCATCAATATACTCCTTATATAAGACTACATTTGTAAtatacta                                                                                     <  1:720953/91‑4 (MQ=255)

        cGACAACTACAACTATTAAAATAAGTGGAACAATTTTTATCTTTTTCATCAATATACTCCTTATATAAGACTACATTTGTAGTATATTACa                                                                                   >  1:1356811/1‑91 (MQ=255)

        cGACAACTACAACTATTAAAATAAGTGGAACAATTTTTATCTTTTTCATCAATATACTCCTTATATAAGACTACATTTGTAGTATATTACa                                                                                   >  1:1356896/1‑91 (MQ=255)

          acaactacaactATTAAAATAAGTGGAACAATTTTTATCTTTTTCATCAATATACTCCTTATATAAGACTACATTTGTAGTATATTACaaa                                                                                 >  1:3891640/1‑91 (MQ=255)

           caactacaactaTTAAAATAAGTGGAACAATTTTTATCTTTTTCATCAATATACTCCTTATATAAGACTACATTTGTAGTATATTACAAAt                                                                                >  1:3095117/1‑91 (MQ=255)

           caactacaactaTTAAAATAAGTGGAACAATTTTTATCTTTTTCATCAATATACTCCTTATATAAGACTACATTTGTAGTATATTACAAAt                                                                                <  1:389902/91‑1 (MQ=255)

           caactacaactaTTAAAATAAGTGGAACAATTTTTATCTTTTTCATCAATATACTCCTTATATAAGACTACATTTGTAGTATATTACAAAt                                                                                >  1:2433281/1‑91 (MQ=255)

           caactacaactaTTAAAATAAGTGGAACAATTTTTATCTTTTTCATCAATATACTCCTTATATAAGACTACATTTGTAGTATATTACAAAt                                                                                >  1:634998/1‑91 (MQ=255)

            aactacaactaTTAAAATAAGTGGAACAATTTTTATCTTTTTCATCAATATACTCCTTATATAAGACTACATTTGTAATATACTACAAATg                                                                               <  1:1624225/91‑1 (MQ=255)

             actacaactaTTAAAATAAGTGGAACAATTTTTATCTTTTTCATCAATATACTCCTTATATAAGACTACATTTGTAGTATATTACAAATgt                                                                              <  1:3125825/91‑1 (MQ=255)

             actacaactaTTAAAATAAGTGGAACAATTTTTATCTTTTTCATCAATATACTCCTTATATAAGACTACATTTGTAGTATATTACAAATgt                                                                              <  1:3125769/91‑1 (MQ=255)

             actacaactaTTAAAATAAGTGGAACAATTTTTATCTTTTTCATCAATATACTCCTTATATAAGACTACATTTGTAGCATATTACAAATgt                                                                              <  1:3124710/91‑1 (MQ=255)

                    ctaTTAAAATAAGTGGAACAATTTTTATCTTTTTCATCAATATACTCCTTATATAAGACTACATTTGTAGTATATTACAAATGTAGTAttt                                                                       <  1:810549/91‑1 (MQ=255)

                       ttAAAATAAGTGGAACAATTTTTATCTTTTTCATCAATATACTCCTTATATAAGACTACATTTGTAGTATATTACAAATGTAGTATTTATg                                                                    >  1:1185698/1‑91 (MQ=255)

                       ttAAAATAAGTGGAACAATTTTTATCTTTTTCATCAATATACTCCTTATATAAGACTACATTTGTAGTATATTACAAATGTAGTATTTATg                                                                    >  1:2628741/1‑91 (MQ=255)

                       ttAAAATAAGTGGAACAATTTTTATCTTTTTCATCAATATACTCCTTATATAAGACTACATTTGTAGTATATTACAAATGTAGTATTTATg                                                                    >  1:2630132/1‑91 (MQ=255)

                       ttAAAATAAGTGGAACAATTTTTATCTTTTTCATCAATATACTCCTTATATAAGACTACATTTGTAGTATATTACAAATGTAGTATTTATg                                                                    >  1:12880/1‑91 (MQ=255)

                         aaaaTAAGTGGAACAATTTTTATCTTTTTCATCAATATACTCCTTATATAAGACTACATTTGTAGTATATTACAAATGTAGTATTTATGTc                                                                  >  1:461880/1‑91 (MQ=255)

                          aaaTAAGTGGAACAATTTTTATCTTTTTCATCAATATACTCCTTATATAAGACTACATTTGTAGTATATTACAAATGTAGTATTTATGTCa                                                                 >  1:93012/1‑91 (MQ=255)

                           aaTAAGTGGAACAATTTTTATCTTTTTCATCAATATACTCCTTATATAAGACTACATTTGTAGTATATTACAAATGTAGTATTTATGTCaa                                                                >  1:1053564/1‑91 (MQ=255)

                              aaGTGGAACAATTTTTATCTTTTTCATCAATATACTCCTTATATAAGACTACATTTGTAGTATATTACAAATGTAGTATTTATGTCAAaat                                                             <  1:169724/91‑1 (MQ=255)

                              aaGTGGAACAATTTTTATCTTTTTCATCAATATACTCCTTATATAAGACTACATTTGTAGTATATTACAAATGTAGTATTTATGTCAAaat                                                             >  1:649159/1‑91 (MQ=255)

                                gTGGAACAATTTTTATCTTTTTCATCAATATACTCCTTATATAAGACTACATTTGTAGTATATTACAAATGTAGTATTTATGTCAAaataa                                                           >  1:1842909/1‑91 (MQ=255)

                                  ggAACAATTTTTATCTTTTTCATCAATATACTCCTTATATAAGACTACATTTGTAGTATATTACAAATGTAGTATTTATGTCAAAATAATg                                                         >  1:2285655/1‑91 (MQ=255)

                                  ggAACAATTTTTATCTTTTTCATCAATATACTCCTTATATAAGACTACATTTGTAGTATATTACAAATGTAGTATTTATGTCAAAATAATg                                                         >  1:1211714/1‑91 (MQ=255)

                                  ggAACAATNTTTNTCTTTTTCATCAATATACTCCTTATATAAGACTACATTTGTAGTATATTACAAATGTAGTATTTATGTCAAAATAATg                                                         <  1:422027/91‑1 (MQ=255)

                                   gAACAATTTTTATCTTTTTCATCAATATACTCCTTATATAAGACTACATTTGTAGTATATTACAAATGTAGTATTTATGTCAAAATAATGt                                                        >  1:379116/1‑91 (MQ=255)

                                    aaCAATTTTTATCTTTTTCATCAATATACTCCTTATATAAGACTACATTTGTAGTATATTACAAATGTAGTATTTATGTCATAATAATGtt                                                       <  1:2091198/91‑1 (MQ=255)

                                    aaCAATTTTTATCTTTTTCATCAATATACTCCTTATATAAGACTACATTTGTAGTATATTACAAATGTAGTATTTATGTCAAAATAATGtt                                                       <  1:869004/91‑1 (MQ=255)

                                    aaCAATTTTTATCTTTTTCATCAATATACTCCTTATATAAGACTACATTTGTAGTATATTACAAATGTAGTATTTATGTCAAAATAATGtt                                                       <  1:2929081/91‑1 (MQ=255)

                                    aaCAATTTTTATCTTTTTCATCAATATACTCCTTATATAAGACTACATTTGTAGTATATTACAAATGTAGTATTTATGTCAAAATAATGtt                                                       <  1:1088068/91‑1 (MQ=255)

                                    aaCAATTTTTATCTTTTTCATCAATATACTCCTTATATAAGACTACATTTGTAGTATATTACAAATGTAGTATTTATGTCAAAATAATGtt                                                       <  1:1157088/91‑1 (MQ=255)

                                    aaCAATTTTTATCTTTTTCATCAATATACTCCTTATATAAGACTACATTTGTAGTATATTACAAATGTAGTATTTATGTCAAAATAATGtt                                                       <  1:1541637/91‑1 (MQ=255)

                                    aaCAATTTTTATCTTTTTCATCAATATACTCCTTATATAAGACTACATTTGTAGTATATTACAAATGTAGTATTTATGTCAAAATAATGtt                                                       <  1:3543157/91‑1 (MQ=255)

                                      cAATTTTTATCTTTTTCATCAATATACTCCTTATATAAGACTACATTTGTAGTATATTACAAATGTAGTATTTATGTCAAAATAATGTtat                                                     <  1:1842112/91‑1 (MQ=255)

                                      cAATTTTTATCTTTTTCATCAATATACTCCTTATATAAGACTACATTTGTAGTATATTACAAATGTAGTATTTATGTCAAAATAATGTtat                                                     <  1:1784841/91‑1 (MQ=255)

                                       aaTTTTTATCTTTTTCATCAATATACTCCTTATATAAGACTACATTTGTAGTATATTACAAATGTAGTATTTATGTCAAAATAATGTtata                                                    <  1:3139337/91‑1 (MQ=255)

                                        aTTTTTATCTTTTTCATCAATATACTCCTTATATAAGACTACATTTGTAGTATATTACAAATGTAGTATTTATGTCAAAATAATGTTATaa                                                   >  1:2811474/1‑91 (MQ=255)

                                         tttttATCTTTTTCATCAATATACTCCTTATATAAGACTACATTTGTAGTATATTACAAATGTAGTATTTATGTCAAAATAATGTTATAAt                                                  <  1:2676257/91‑1 (MQ=255)

                                           tttATCTTTTTCATCAATATACTCCTTATATAAGACTACATTTGTAGTATATTACAAATGTAGTATTTATGTCAAAATAATGTTATAAttt                                                >  1:1026460/1‑91 (MQ=255)

                                            ttATCTTTTTCATCAATATACTCCTTATATAAGACTACATTTGTAGTATATTACAAATGTAGTATTTATGTCAAAATAATGTTATAAtttt                                               >  1:1112407/1‑91 (MQ=255)

                                            ttATCTTTTTCATCAATATACTCCTTATATAAGACTACATTTGTAGTATATTACAAATGTAGTATTTATGTCAAAATAATGTTATAAtttt                                               <  1:3510336/91‑1 (MQ=255)

                                             tATCTTTTTCATCAATATACTCCTTATATAAGACTACATTTGTAGTATATTACAAATGTAGTATTTATGTCAAAATAATGTTATAAttttt                                              >  1:2860091/1‑91 (MQ=255)

                                               tctcTTTCATCAATATACTCCTTAAATAAGACTACATTTGTAGTATATTACAAATGTAGGATCCATGTCAAAATAATGTTATAATTTTtgt                                            <  1:2976150/91‑1 (MQ=255)

                                               tCTTTTTCATCAATATACTCCTTATATAAGACTACATTTGTAGTATATTACAAATGTAGTATTTATGTCAAAATAATGTTATAATTTTtgt                                            <  1:566001/91‑1 (MQ=255)

                                               tCTTTTTCATCAATATACTCCTTATATAAGACTACATTTGTAGTATATTACAAATGTAGTATTTATGTCAAAATAATGTTATAATTTTtgt                                            <  1:986049/91‑1 (MQ=255)

                                               tCTTTTTCATCAATATACTCCTTATATAAGACTACATTTGTAGTATATTACAAATGTAGTATTTATGTCAAAATAATGTTATAATTTTtgt                                            >  1:2750426/1‑91 (MQ=255)

                                               tCTTTTTCATCAATATACTCCTTATATAAGACTACATTTGTAGTATATTACAAATGTAGTATTTATGTCAAAATAATGTTATAATTTTtgt                                            >  1:2750191/1‑91 (MQ=255)

                                                   tttCATCAATATACTCCTTATATAAGACTACATTTGTAGTATATTACAAATGTAGTATTTATGTCAAAATAATGTTATAATTTTTGTGata                                        >  1:2646421/1‑91 (MQ=255)

                                                   tttCATCAATATACTCCTTATATAAGACTACATTTGTAGTATATTACAAATGTAGTATTTATGTCAAAATAATGTTATAATTTTTGTGata                                        <  1:1844665/91‑1 (MQ=255)

                                                   tttCATCAATATACTCCTTATATAAGACTACATTTGTAGTATATTACAAATGTAGTATTTATGTCAAAATAATGTTATAATTTTTGTGata                                        <  1:2021960/91‑1 (MQ=255)

                                                    ttCATCAATATACTCCTCATATAAGACTACATTTGTAGTATATTACAAATGTAGTATTTATGTCAAAATAATGTTATAATTTTTGTGatat                                       <  1:3556564/91‑1 (MQ=255)

                                                    ttCATCAATATACTCCTCATATAAGACTACATTTGTAGTATATTACAAATGTAGTATTTATGTCAAAATAATGTTATAATTTTTGTGatat                                       <  1:3556421/91‑1 (MQ=255)

                                                        tcaATATACTCCTTATATAAGACTACATTTGTAGTATATTACAAATGTAGTATTTATGTCAAAATAATGTTATAATTTTTGTGATATGGAg                                   <  1:2694215/91‑1 (MQ=255)

                                                         caATATACTCCTTATATAAGACTACATTTGTAGTATATTACAAATGTAGTATTTATGTCAAAATAATGTTATAATTTTTGTGAAATGGAgg                                  >  1:2077465/1‑91 (MQ=255)

                                                          aaTATACTCCTTATATAAGACTACATTTGTAGTATATTACAAATGTAGTATTTATGTCAAAATAATGTTATAATTTTTGTGATATGGAGgt                                 <  1:2634061/91‑1 (MQ=255)

                                                          aaTATACTCCTTATATAAGACTACATTTGTAGTATATTACAAATGTAGTATTTATGTCAAAATAATGTTATAATTTTTGTGATATGGAGgt                                 >  1:1414036/1‑91 (MQ=255)

                                                          aaTATACTCCTTATATAAGACTACATTTGTAGTATATTACAAATGTAGTATTTATGTCAAAATAATGTTATAATTTTTGTGATATGGAGgt                                 >  1:85107/1‑91 (MQ=255)

                                                           atatACTCCTTATATAAGACTACATTTGTAGTATATTACAAATGTAGTATTTATGTCAAAATAATGTTATAATTTTTGTGATATGGAGgtg                                >  1:3781354/1‑91 (MQ=255)

                                                           atatACTCCTTATATAAGACTACATTTGTAGTATATTACAAATGTAGTATTTATGTCAAAATAATGTTATAATTTTTGTGATATGGAGgtg                                >  1:2034827/1‑91 (MQ=255)

                                                           atatACTCCTTATATAAGACTACATTTGTAGTATATTACAAATGTAGTATTTATGTCAAAATAATGTTATAATTTTTGTGATATGGAGgtg                                >  1:3303407/1‑91 (MQ=255)

                                                             ataCTCCTTATATAAGACTACATTTGTAGTATATTACAAATGTAGTATTTATGTCAAAATAATGTTATAATTTTTGTGATATGGAGGTGTa                              >  1:3179319/1‑91 (MQ=255)

                                                                cTCCTTATATAAGACTACATTTGTAGTATATTACAAATGTAGTATTTATGTCAAAATAATGTTATAATTTTTGTGATATGGAGGTGTAGaa                           >  1:3596269/1‑91 (MQ=255)

                                                                cTCCTTATATAAGACTACATTTGTAGTATATTACAAATGTAGTATTTATGTCAAAATAATGTTATAATTTTTGTGATATGGAGGTGTAGaa                           >  1:3242269/1‑91 (MQ=255)

                                                                cTCCTTATATAAGACTACATTTGTAGTATATTACAAATGTAGTATTTATGTCAAAATAATGTTATAATTTTTGTGATATGGAGGTGTAGaa                           >  1:4098890/1‑91 (MQ=255)

                                                                 tCCTTATATAAGACTACATTTGTAGTATATTACAAATGTAGTATTTATGTCAAAATAATGTTATAATTTTTGTGATATGGAGGTGTAGAAg                          >  1:2330349/1‑91 (MQ=255)

                                                                 tCCTTATATAAGACTACATTTGTAGTATATTACAAATGTAGTATTTATGTCAAAATAATGTTATAATTTTTGTGATATGGAGGTGTAGAAg                          >  1:2311903/1‑91 (MQ=255)

                                                                 tCCTTATATAAGACTACATTTGTAGTATATTACAAATGTAGTATTTATGTCAAAATAATGTTATAATTTTTGTGATATGGAGGTGTAGAAg                          <  1:1498097/91‑1 (MQ=255)

                                                                     tatataAGACTACATTTGTAGTATATTACAAATGTAGTATTTATGTCAAAATAATGTTATAATTTTTGTGATATGGAGGTGTAGAAGgtgt                      >  1:907967/1‑91 (MQ=255)

                                                                      atataAGACTACATTTGTAGTATATTACAAATGTAGTATTTATGTCAAAATAATGTTATAATTTTTGTGATATGGAGGTGTAGAAGGTGtt                     <  1:1591117/91‑1 (MQ=255)

                                                                       tataAGACTACATTTGTAGTATATTACAAATGTAGTATTTATGTCAAAATAATGTTATAATTTTTGTGATATGGAGGTGTAGAAGGTGTTa                    >  1:2405431/1‑91 (MQ=255)

                                                                       tataAGACTACATTTGTAGTATATTACAAATGTAGTATTTATGTCAAAATAATGTTATAATTTTTGTGATATGGAGGTGTAGAAGGTGTTa                    >  1:2406334/1‑91 (MQ=255)

                                                                       tataAGACTACATTTGTAGTATATTACAAATGTAGTATTTATGTCAAAATAATGTTATAATTTTTGTGATATGGAGGTGTAGAAGGTGTTa                    <  1:562934/91‑1 (MQ=255)

                                                                       tataAGACTACATTTGTAGTATATTACAAATGTAGTATTTATGTCAAAATAATGTTATAATTTTTGTGATATGGAGGTGTAGAAGGTGTTa                    >  1:2461675/1‑91 (MQ=255)

                                                                         taAGACTACATTTGTAGTATATTACAAATGTAGTATTTATGTCAAAATAATGTTATAATTTTTGTGATATGGAGGTGTAGAAGGTGTTatc                  <  1:4082468/91‑1 (MQ=255)

                                                                            gACTACATTTGTAGTATATTACAAATGTAGTATTTATGTCAAAATAATGTTATAATTTTTGTGATATGGAGGTGTAGAAGGTGTTatcatc               >  1:2147823/1‑91 (MQ=255)

                                                                            gACTACATTTGTAGTATATTACAAATGTAGTATTTATGTCAAAATAATGTTATAATTTTTGTGATATGGAGGTGTAGAAGGTGTTatcatc               <  1:2054786/91‑1 (MQ=255)

                                                                               tACATTTGTAGTATATTACAAATGTAGTATTTATGTCAAAATAATGTTATAATTTTTGTGATATGGAGGTGTAGAAGGTGTTATCATCttt            <  1:3291175/91‑1 (MQ=255)

                                                                               tACATTTGTAGTATATTACAAATGTAGTATTTATGTCAAAATAATGTTATAATTTTTGTGATATGGAGGTGTAGAAGGTGTTATCATCttt            <  1:3291291/91‑1 (MQ=255)

                                                                                aCATTTGTAGTATATTACAAATGTAGTATTTATGTCAAAATAATGTTATAATTTTTGTGATATGGAGGTGTAGAAGGTGTTATCATCtttt           >  1:974680/1‑91 (MQ=255)

                                                                                aCATTTGTAATATACTACAAATGTAGTATTTATGTCAAAATAATGTTATAATTTTTGTGATATGGAGGTGTAGAAGGTGTTATCATCtttt           >  1:775735/1‑91 (MQ=255)

                                                                                aCATTTGTAATATACTACAAATGTAGTATTTATGTCAAAATAATGTTATAATTTTTGTGATATGGAGGTGTAGAAGGTGTTATCATCtttt           >  1:406530/1‑91 (MQ=255)

                                                                                aCATTTGTAATATACTACAAATGTAGTATTTATGTCAAAATAATGTTATAATTTTTGTGATATGGAGGTGTAGAAGGTGTTATCATCtttt           >  1:2776684/1‑91 (MQ=255)

                                                                                 cATTTGTAGTATATTACAAATGTAGTATTTATGTCAAAATAATGTTATAATTTTTGTGATATGGAGGTGTAGAAGGTGTTATCATCttttt          <  1:1084803/91‑1 (MQ=255)

                                                                                   tttGTAGTATATTACAAATGTAGTATTTATGTCAAAATAATGTTATAATTTTTGTGATATGGAGGTGTAGAAGGTGTTATCATCTTTTTTa        >  1:2563806/1‑91 (MQ=255)

                                                                                     tGTAGTATATTACAAATGTAGTATTTATGTCAAAATAATGTTATAATTTTTGTGATATGGAGGTGTAGAAGGTGTTATCATCTTTTTTAAt      >  1:203298/1‑91 (MQ=255)

                                                                                     tGTAGTATATTACAAATGTAGTATTTATGTCAAAATAATGTTATAATTTTTGTGATATGGAGGTGTAGAAGGTGTTATCATCTTTTTTAAt      >  1:1233470/1‑91 (MQ=255)

                                                                                      gtagtaTATTACAAATGTAGTATTTATGTCAAAATAATGTTATAATTTTTGTGATATGGAGGTGTAGAAGGTGTTATCATCTTTTTTAATg     <  1:4008565/91‑1 (MQ=255)

                                                                                      gtagtaTATTACAAATGTAGTATTTATGTCAAAATAATGTTATAATTTTTGTGATATGGAGGTGTAGAAGGTGTTATCATCTTTTTTAATg     <  1:3709354/91‑1 (MQ=255)

                                                                                       tagtaTATTACGAATGTGGGATTTATGTCGAAATAATGTTAGAATTTTTGTGGTATGGAGGTGTAGAAGGTGTTATCATCTTTTTTAATGt    <  1:3333042/91‑1 (MQ=255)

                                                                                        agtaTATTACAAATGTAGTATTTATGTCAAAATAATGTTATAATTTTTGTGATATGGAGGTGTAGAAGGTGTTATCATCTTTTTTAATGtt   >  1:1476668/1‑91 (MQ=255)

                                                                                        agtaTATTACAAATGTAGTATTTATGTCAAAATAATGTTATAATTTTTGTGATATGGAGGTGTAGAAGGTGTTATCATCTTTTTTAATGtt   <  1:345111/91‑1 (MQ=255)

                                                                                         gtaTATTACAAATGTAGTATTTATGTCAAAATAATGTTATAATTTTTGTGATATGGAGGTGTAGAAGGTGTTATCATCTTTTTTAATGTTa  >  1:2174465/1‑91 (MQ=255)

                                                                                         gtaTATTACAAATGTAGTATTTATGTCAAAATAATGTTATAATTTTTGTGATATGGAGGTGTAGAAGGTGTTATCATCTTTTTTAATGTTa  >  1:2174494/1‑91 (MQ=255)

ACCAAACCCGACAACTACAACTATTAAAATAAGTGGAACAATTTTTATCTTTTTCATCAATATACTCCTTATATAAGACTACATTTGTAGTATATTACAAATGTAGTATTTATGTCAAAATAATGTTATAATTTTTGTGATATGGAGGTGTAGAAGGTGTTATCATCTTTTTTAATGTTA  >  USA300TCH1516_ALE/41077‑41256

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