BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	USA300TCH1516_ALE	92,987	A→G	15.7%	intergenic (‑72/‑65)	ricR_1 ← / → glpE	Copper‑sensing transcriptional repressor RicR/Thiosulfate sulfurtransferase GlpE

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	USA300TCH1516_ALE	92,987	0	A	G	15.7%	254.1 / 39.7	107	intergenic (‑72/‑65)	ricR_1/glpE	Copper‑sensing transcriptional repressor RicR/Thiosulfate sulfurtransferase GlpE
Reads supporting (aligned to +/- strand): ref base A (38/52); new base G (15/2); total (53/54)
Fisher's exact test for biased strand distribution p-value = 4.71e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.54e-01

Only 100 of 114 total aligned reads displayed.

TTTTTTTATCATAATTCATAGTCATACCTCCACTTTTAATTGAATAAAAATATATTAATAGATAAACACAAATGTGTCAAATACCCCTAGAGGTATTTGACAAGTTCCATCCAACTGTTTAAAATACCCCTACAGGTATTTTTAGGGAGGTTATTATGAAACAATACGGAGAAAAGTTTA  >  USA300TCH1516_ALE/92897‑93076

tttttttATCATAATTCATAGTCATACCTCCACTTTTAATTGAATAAAAATATATTAATAGATAAACACAAATGTGTCAAATACCCCTaga                                                                                           <  1:904941/91‑1 (MQ=255)

tttttttATCATAATTCATAGTCATACCTCCACTTTTAATTGAATAAAAATATATTAATAGATAAACACAAATGTGTCAAATACCCCTaga                                                                                           <  1:1153176/91‑1 (MQ=255)

 ttttttATCATAATTCATAGTCATACCTCCACTTTTAATTGAATAAAAATATATTAATAGATAAACACAAATGTGTCAAATACCCCTagag                                                                                          >  1:2446154/1‑91 (MQ=255)

 ttttttATCATAATTCATAGTCATACCTCCACTTTTAATTGAATAAAAATATATTAATAGATAAACACAAATGTGTCAAATACCCCTagag                                                                                          >  1:1203022/1‑91 (MQ=255)

 tgaattATCATAATTCATAGTCATACCTCCACTTTTAATTGAATAAAAATATATTAATAGATAAACACAAATGTGTCAAATACCCCTagag                                                                                          >  1:1696264/5‑91 (MQ=255)

      tATCATAATTCATAGTCATACCTCCACTTTTAATTGAATAAAAATATATTAATAGATAAACACAAATGTGTCAAATACCCCTAGAGGTAtt                                                                                     <  1:2765284/91‑1 (MQ=255)

      tATCATAATTCATAGTCATACCTCCACTTTTAATTGAATAAAAATATATTAATAGATAAACACAAATGTGTCAAATACCCCTAGAGGTAtt                                                                                     >  1:2855691/1‑91 (MQ=255)

      tATCATAATTCATAGTCATACCTCCACTTTTAATTGAATAAAAATATATTAATAGATAAACACAAATGTGTCAAATACCCCTAGAGGTAtt                                                                                     <  1:3465973/91‑1 (MQ=255)

         cATAATTCATAGTCATACCTCCACTTTTAATTGAATAAAAATATATTAATAGATAAACACAAATGTGTCAAATACCCCTAGAGGTATTTGa                                                                                  <  1:352569/91‑1 (MQ=255)

          aTAATTCATAGTCATACCTCCACTTTTAATTGAATAAAAATATATTAATAGATAAACACAAATGTGCCAAATACCTCTAGGGGTATTTGAc                                                                                 <  1:1964768/91‑1 (MQ=39)

           tAATTCATAGTCATACCTCCACTTTTAATTGAATAAAAATATATTAATAGATAAACACAAATGTGTCAAATACCCCTAGAGGTATTTGACa                                                                                >  1:3588927/1‑91 (MQ=255)

           tAATTCATAGTCATACCTCCACTTTTAATTGAATAAAAATATATTAATAGATAAACACAAATGTGTCAAATACCCCTAGAGGTATTTGACa                                                                                <  1:2744087/91‑1 (MQ=255)

              ttCATAGTCATACCTCCACTTTTAATTGAATAAAAATATATTAATAGATAAACACAAATGTGTCAAATACCCCTAGAGGTATTTGACAAGt                                                                             <  1:205014/91‑1 (MQ=255)

               tCATAGTCATACCTCCACTTTTAATTGAATAAAAATATATTAATAGATAAACACAAATGTGTCAAATACCCCTAGAGGTATTTGACAAGtt                                                                            >  1:3731394/1‑91 (MQ=255)

                   aGTCATACCTCCACTTTTAATTGAATAAAAATATATTAATAGATAAACACAAATGTGTCAAATACCCCTAGAGGTATTTGACAAGTtccat                                                                        >  1:1945276/1‑91 (MQ=255)

                   aGTCATACCTCCACTTTTAATTGAATAAAAATATATTAATAGATAAACACAAATGTGTCAAAGACCCCTAGAGGTATTTGACAAGTtccat                                                                        <  1:1837088/91‑1 (MQ=255)

                   aGTCATACCTCCACTTTTAATTGAAGAAAAATATATTAATAGATAAACACAAATGTGTCAAATACCCCTAGAGGTATTTGACAAGTtccat                                                                        <  1:3851244/91‑1 (MQ=255)

                     tCATACCTCCACTTTTAATTGAATAAAAATATATTAATAGATAAACACAAATGTGTCAAATACCCCTAGAGGTATTTGACAAGTtccatcc                                                                      <  1:1849164/91‑1 (MQ=255)

                     tCATACCTCCACTTTTAATTGAATAAAAATATATTAATAGATAAACACAAATGTGTCAAATACCCCTAGAGGTATTTGACAAGTtccatcc                                                                      >  1:962203/1‑91 (MQ=255)

                      cATACCTCCACTTTTAATTGAATAAAAATATATTAATAGATAAACACAAATGTGTCAAATACCCCTAGAGGTATTTGACAAGTtccatcca                                                                     <  1:119757/91‑1 (MQ=255)

                      cATACCTCCACTTTTAATTGAATAAAAATATATTAATAGATAAACACAAATGTGTCAAATACCCCTAGAGGTATTTGACAAGTtccatcca                                                                     <  1:2724261/91‑1 (MQ=255)

                        tACCTCCACTTTTAATTGAATAAAAATATATTAATAGATAAACACAAATGTGTCAAATACCCCTAGAGGTATTTGACAAGTTCCATCCAAc                                                                   >  1:1336487/1‑91 (MQ=255)

                         aCCTCCACTTTTAATTGAATAAAAATATATTAATAGATAAACACAAATGTGTCAAATACCCCTAGAGGTATTTGACAAGTTCCATCCAACt                                                                  >  1:2092008/1‑91 (MQ=255)

                           cTCCACTTTTAATTGAATAAAAATATATTAATAGATAAACACAAATGTGTCAAATACCCCTAGAGGTATTTGACAAGTTCCATCCAACTGt                                                                >  1:2343389/1‑91 (MQ=255)

                            tCCACTTTTAATTGAATAAAAATATATTAATAGATAAACACAAATGTGTCAAATACCCCTAGAGGTATTTGACAAGTTCCATCCAACTGtt                                                               >  1:1309097/1‑91 (MQ=255)

                             ccACTTTTAATTGAATAAAAATATATTAATAGATAAACACAAATGTGTCAAATACCCCTAGAGGTATTTGACAAGTTCCATCCAACTGttt                                                              <  1:2449001/91‑1 (MQ=255)

                             ccACTTTTAATTGAATAAAAATATATTAATAGATAAACACAAATGTGTCAAATACCCCTAGAGGTATTTGACAAGTTCCATCCAACTGttt                                                              <  1:1560501/91‑1 (MQ=255)

                                 ttttAATTGAATAAAAATATATTAATAGATAAACACAAATGTGTCAAATACCCCTAGAGGTATTTGACAAGTTCCATCCAACTGTTTaaaa                                                          <  1:3662812/91‑1 (MQ=255)

                                 ttttAATTGAATAAAAATATATTAATAGATAAACACAAATGTGTCAAATACCCCTAGAGGTATTTGACAAGTTCCATCCAACTGTTTaaaa                                                          <  1:2651340/91‑1 (MQ=255)

                                  tttAATTGAATAAAAATATATTAATAGATAAACACAAATGTGTCAAATACCCCTAGAGGTATTTGACAAGTTCCATCCAACTGTTTAAAAt                                                         <  1:1621270/91‑1 (MQ=255)

                                  tttAATTGAATAAAAATATATTAATAGATAAACACAAATGTGTCAAATACCCCTAGAGGTATTTGACAAGTTCCATCCAACTGTTTAAAAt                                                         <  1:1550821/91‑1 (MQ=255)

                                    tAATTGAATAAAAATATATTAATAGATAAACACAAATGTGTCAAATACCCCTAGAGGTATTTGACAAGTTCCATCCAACTGTTTAAAATAc                                                       >  1:4107940/1‑91 (MQ=255)

                                    tAATTGAATAAAAATATATTAATAGATAAACACAAATGTGTCAAATACCCCTAGAGGTATTTGACAAGTTCCATCCAACTGTTTAAAATAc                                                       >  1:1084700/1‑91 (MQ=255)

                                     aaTTGAATAAAAATATATTAATAGATAAACACAAATGTGTCAAATACCCCTAGAGGTATTTGACAAGTTCCATCCAACTGTTTAAAATAcc                                                      >  1:2461651/1‑91 (MQ=255)

                                     aaTTGAATAAAAATATATTAATAGATAAACACAAATGTGTCAAATACCCCTAGAGGTATTTGACAAGTTCCATCCAACTGTTTAAAATAcc                                                      >  1:2145858/1‑91 (MQ=255)

                                       ttGAATAAAAATATATTAATAGATAAACACAAATGTGTCAAATACCCCTAGAGGTATTTGACAAGTCACATCCTACTGTTTAAAATAcccc                                                    >  1:1687128/1‑91 (MQ=255)

                                          aaTAAAAATATATTAATAGATAAACACAAATGTGTCAAATACCCTTAGAGGTATTTGACAAGTTCCATCCAACTGTTTAAAATACCCCTAc                                                 <  1:3283755/91‑1 (MQ=255)

                                          aaTAAAAATATATTAATAGATAAACACAAATGTGTCAAATACCCCTAGAGGTATTTGACAAGTTCCATCCAACTGTTTAAAATACCCCTAc                                                 >  1:1862314/1‑91 (MQ=255)

                                            tAAAAATATATTAATAGATAAACACAAATGTGTCAAATACCCCTAGAGGTATTTGACACGTTCCATCCAACTGTTTAAAATACCCCTACAg                                               >  1:1985849/1‑91 (MQ=255)

                                             aaaaaTATATTAATAGATAAACACAAATGTGTCAAATACCCCTAGAGGTATTTGACAAGTTCCATCCAACTGTTTAAAATACCCCTACAgg                                              <  1:813941/91‑1 (MQ=255)

                                             aaaaaTATATTAATAGATAAACACAAATGTGTCAAATACCCCTAGAGGTATTTGACAAGTTCCATCCAACTGTTTAAAATACCCCTACAgg                                              <  1:1238248/91‑1 (MQ=255)

                                               aaaTATATTAATAGATAAACACAAATGTGTCAAATACCCCTAGAGGTATTTGACAAGTTCCATCCAACTGTTTAAAATACCCCTACAGGTa                                            >  1:1556514/1‑91 (MQ=255)

                                                aaTATATTAATAGATAAACACAAATGTGTCAAATACCCCTAGAGGTATTTGACAAGTTCCATCCAACTGTTTAAAATACCCCTACAGGTAt                                           <  1:3922858/91‑1 (MQ=255)

                                                 atatatTAATAGATAAACACAAATGTGTCAAATACCCCTAGAGGTATTTGACAAGTTCCATCCAACTGTTTAAAATACCCCTACAGGTAtt                                          <  1:1975779/91‑1 (MQ=255)

                                                  tatatTAATAGATAAACACAAATGTGTCAAATACCCCTAGAGGTATTTGACAAGTTCCATCCAACTGTTTAAAATACCCCTACAGGTAttt                                         >  1:1574461/1‑91 (MQ=255)

                                                    tatTAATAGATAAACACAAATGTGTTAAATACCCCTAGAGGTATTTGACAAGTTCCATCCAACTGTTTAAAATACCCCTACAGGTAttttt                                       <  1:3869592/91‑1 (MQ=255)

                                                    tatTAATAGATAAACACAAATGTGTCAAATACCCCTAGAGGTATTTGACAAGTTCCATCCAACTGTTTAAAATACCCCTACAGGTAttttt                                       >  1:3152014/1‑91 (MQ=255)

                                                      ttAATAGATAAACACAAATGTGTCAAATACCCCTAGAGGTATTTGACAAGTTCCATCCAACTGTTTAAAATACCCCTACAGGTATTTTTAg                                     >  1:804049/1‑91 (MQ=255)

                                                       tAATAGATAAACACAAATGTGTCAAATACCCCTAGAGGTATTTGACAAGTTCCATCCAACTGTTTAAAATACCCCTACAGGTATTTTTAgg                                    <  1:1308368/91‑1 (MQ=255)

                                                        aaTAGATAAACACAAATGTGTCAAATACCCCTAGAGGTATTTGACAAGTTCCATCCAACTGTTTAAAATACCCCTACAGGTATTTTTAggg                                   >  1:1678004/1‑91 (MQ=255)

                                                         aTAGATAAACACAAATGTGTCAAATACCCCTAGAGGTATTTGACAAGTTCCATCCAACTGTTTAAAATACCCCTACAGGTATTTTTAGGGa                                  <  1:2221152/91‑1 (MQ=255)

                                                           aGATAAACACAAATGTGTCAAATACCCCTAGAGGTATTTGACAAGTTCCATCCAACTGTTTAAAATACCCCTACAGGTATTTTTAGGGAgg                                <  1:1152972/91‑1 (MQ=255)

                                                            gATAAACACAAATGTGTCAAATACCCCTAGAGGTATTTGACAAGTTCCATCCAACTGTTTAAAATACCCCTACAGGTATTTTTAGGGAGGt                               <  1:2402082/91‑1 (MQ=255)

                                                            gATAAACACAAATGTGTCAAATACCCCTAGAGGTATTTGACAAGTTCCATCCAACTGTTTAAAATACCCCTACAGGTATTTTTAGGGAGGt                               <  1:2402083/91‑1 (MQ=255)

                                                             aTAAACACACATGTGTCAAATACCCCTAGAGGTATTTGACAAGTTCCATCCAACTGTTTAAAATACCCCTACAGGTATTTTTAGGGAGGtt                              <  1:2745663/91‑1 (MQ=255)

                                                             aTAAACACAAATGTGTCAAATACCCCTAGGGGTATTTGACAAGTTCCATCCAACTGTTTAAAATACCCCTACAGGTATTTTTAGGGAGGtt                              <  1:1847588/91‑1 (MQ=255)

                                                             aTAAACACAAATGTGTCAAATACCCCTAGAGGTATTTGACAAGTTCCATCCAAATGTTTAAAATACCCCTACAGGTATTTTTAGGGAGGtt                              <  1:3497024/91‑1 (MQ=255)

                                                                 acacAAATGTGTCAAATACCCCTAGAGGTATTTGACAAGTTCCATCCAACTGTTTAAAATACCCCTACAGGTATTTTTAGGGAGGttatta                          >  1:3513357/1‑91 (MQ=255)

                                                                  cacaAATGTGTCAAATACCCCTAGAGGTATTTGACAAGTTCCATCCAACTGTTTAAAATACCCCTACAGGTATTTTTAGGGAGGTtattat                         >  1:3361817/1‑91 (MQ=255)

                                                                  cacaAATGTGTCAAATACCCCTAGAGGTATTTGACAAGTTCCATCCAACTGTTTAAAATACCCCTACAGGTATTTTTAGGGAGGTtattat                         <  1:89894/91‑1 (MQ=255)

                                                                     aaaTGTGTCAAATACCCCTAGAGGTATTTGACAAGTTCCATCCAACTGTTTAAAATACCCCTACAGGTATTTTTAGGGAGGTTATTATGaa                      >  1:21858/1‑91 (MQ=255)

                                                                     aaaTGTGTCAAATACCCCTAGAGGTATTTGACAAGTTCCATCCAACTGTTTAAAATACCCCTACAGGTATTTTTAGGGAGGTTATTATGaa                      <  1:1728741/91‑1 (MQ=255)

                                                                      aaTGTGTCAAATACCCCTAGAGGTATTTGACAAGTTCCATCCAACTGTTTAAAATACCCCTACAGGTATTTTTAGGGAGGTTATTATGaaa                     <  1:428060/91‑1 (MQ=255)

                                                                       aTGTGTCAAATACCCCTAGAGGTATTTGACAAGTTCCATCCAACTGTTTAAAATACCCCTACAGGTATTTTTAGGGAGGTTATTATGAAAc                    >  1:1534815/1‑91 (MQ=255)

                                                                        tgtgTCAAATACCCCTAGAGGTATTTGACAAGTTCCATCCAACTGTTTAAAATACCCCTACAGGTATTTTTAGGGAGGTTATTATGAAACa                   >  1:2129992/1‑91 (MQ=255)

                                                                         gtgtCAATTACCCCTAGAGGTATTTGACAAGTTCCATCCAACTGTTTAAAATACCCCTACAGGTATTTTTAGGGAGGTTATTATGAAACaa                  <  1:2065562/91‑1 (MQ=255)

                                                                         gtgtCAAATACCCCTAGAGGTATTTGACAAGTTCCATCCAACTGTTTAAAATACCCCTACAGGTATTTTTAGGGAGGTTATTATGAAACaa                  <  1:1065885/91‑1 (MQ=255)

                                                                         gtgtCAAATACCCCTAGAGGTATTTGACAAGTTCCATCCAACTGTTTAAAATACCCCTACAGGTATTTTTAGGGAGGTTATTATGAAACaa                  >  1:2268494/1‑91 (MQ=255)

                                                                          tgtCAAATACCCCTAGAGGTATTTGACAAGTTCCATCCAACTGTTTAAAATACCCCTACAGGTATTTTTAGGGAGGTTATTATGAAACAAt                 <  1:620209/91‑1 (MQ=255)

                                                                            tCAAATACCCCTAGAGGTATTTGACAAGTTCCATCCAACTGTTTAAAATACCCCTACAGGTATTTTTAGGGAGGTTATTATGAAACAATAc               >  1:204160/1‑91 (MQ=255)

                                                                             cAAATACCCCTAGAGGTATTTGACAAGTTCCATCCAACTGTTTAAAATACCCCTACAGGTATTTTTAGGGAGGTTATTATGAAACAATACg              <  1:557953/91‑1 (MQ=255)

                                                                             cAAATACCCCTAGAGGTATTTGACAAGTTCCATCCAACTGTTTAAAATACCCCTACAGGTATTTTTAGGGAGGTTATTATGAAACAATACg              <  1:557998/91‑1 (MQ=255)

                                                                             cAAATACCCCTAGAGGTATTTGACAAGTTCCATCCAACTGTTTAAAATACCCCTACAGGTATTTTTAGGGAGGTTATTATGAAACAATACg              <  1:1654426/91‑1 (MQ=255)

                                                                              aaaTACCCCTAGAGGTATTTGACAAGTTCCATCCAACTGTTTAAAATACCCCTACAGGTATTTTTAGGGAGGTTATTATGAAACAATACgg             >  1:1939688/1‑91 (MQ=255)

                                                                              aaaTACCCCTAGAGGTATTTGACAAGTTCCATCCAACTGTTTAAAATACCCCTACAGGTATTTTTAGGGAGGTTATTATGAAACAATACgg             <  1:3379719/91‑1 (MQ=255)

                                                                               aaTACCCCTAGAGGTATTTGACAAGTTCCATCCAACTGTTTAAAATACCCCTACAGGTATTTTTAGGGAGGTTATTATGAAACAATACGga            <  1:1707633/91‑1 (MQ=255)

                                                                                  acctcTAGGGGTATTTGACAAGTTCCATCCAACTGTTTAAAATACCCCTACAGGTATTTTTAGGGAGGTTATTATGAAACAATACGGAGaa         >  1:3428922/5‑91 (MQ=255)

                                                                                  acctcTAGGGGTATTTGACAAGTTCCATCCAACTGTTTAAAATACCCCTACAGGTATTTTTAGGGAGGTTATTATGAAACAATACGGAGaa         >  1:825483/5‑91 (MQ=255)

                                                                                  acctcTAGGGGTATTTGACAAGTTCCATCCAACTGTTTAAAATACCCCTACAGGTATTTTTAGGGAGGTTATTATGAAACAATACGGAGaa         >  1:825467/5‑91 (MQ=255)

                                                                                  acctcTAGGGGTATTTGACAAGTTCCATCCAACTGTTTAAAATACCCCTACAGGTATTTTTAGGGAGGTTATTATGAAACAATACGGAGaa         >  1:255350/5‑91 (MQ=255)

                                                                                  acctcTAGGGGTATTTGACAAGTTCCATCCAACTGTTTAAAATACCCCTACAGGTATTTTTAGGGAGGTTATTATGAAACAATACGGAGaa         >  1:1796950/5‑91 (MQ=255)

                                                                                  aCCCCTAGAGGTATTTGACAAGTTCCATCCAACTGTTTAAAATACCCCTACAGGTATTTTTAGGGAGGTTATTATGAAACAATACGGAGaa         <  1:1957847/91‑1 (MQ=255)

                                                                                  aCCCCTAGAGGTATTTGACAAGTTCCATCCAACTGTTTAAAATACCCCTACAGGTATTTTTAGGGAGGTTATTATGAAACAATACGGAGaa         <  1:2344456/91‑1 (MQ=255)

                                                                                   cctcTAGGGGTATTTGACAAGTTCCATCCAACTGTTTAAAATACCCCTACAGGTATTTTTAGGGAGGTTATTATGAAACAATACGGAGaaa        >  1:3343452/4‑91 (MQ=255)

                                                                                    ctcTAGGGGTATTTGACAAGTTCCATCCAACTGTTTAAAATACCCCTACAGGTATTTTTAGGGAGGTTATTATGAAACAATACGGAGaaaa       >  1:319613/3‑91 (MQ=255)

                                                                                    ctcTAGGGGTATTTGACAAGTTCCATCCAACTGTTTAAAATACCCCTACAGGTATTTTTAGGGAGGTTATTATGAAACAATACGGAGaaaa       >  1:1113046/3‑91 (MQ=255)

                                                                                    cccTAGAGGTATTTGACAAGTTCCATCCAACTGTTTAAAATACCCCTACAGGTATTTTTAGGGAGGTTATTATGAAACAATACGGAGaaaa       >  1:545072/1‑91 (MQ=255)

                                                                                     tcTAGGGGTATTTGACACGTTCCATCCAACTGTTTAAAATACCCCTACAGGTATTTTTAGGGAGGTTATTATGAAACAATACGGAGAAAAg      >  1:1538467/2‑91 (MQ=255)

                                                                                     tcTAGGGGTATTTGACAAGTTCCATCCAACTGTTTAAAATACCCCTACAGGTATTTTTAGGGAGGTTATTATGAAACAATACGGAGAAAAg      >  1:1252463/2‑91 (MQ=255)

                                                                                     tcTAGGGGTATTTGACAAGTTCCATCCAACTGTTTAAAATACCCCTACAGGTATTTTTAGGGAGGTTATTATGAAACAATACGGAGAAAAg      >  1:2765393/2‑91 (MQ=255)

                                                                                     tcTAGGGGTATTTGACAAGTTCCATCCAACTGTTTAAAATACCCCTACAGGTATTTTTAGGGAGGTTATTATGAAACAATACGGAGAAAAg      >  1:3770663/2‑91 (MQ=255)

                                                                                     tcTAGGGGTATTTGACAAGTTCCATCCAACTGTTTAAAATACCCCTACAGGTATTTTTAGGGAGGTTATTATGAAACAATACGGAGAAAAg      >  1:352624/2‑91 (MQ=255)

                                                                                     ccTAGAGGTATTTGACAAGTTCCATCCAACTGTTTAAAATACCCCTACAGGTATTTTTAGGGAGGTTATTATGAAACAATACGGAGAAAAg      <  1:1362042/91‑1 (MQ=255)

                                                                                     ccTAGAGGTATTTGACAAGTTCCATCCAACTGTTTAAAATACCCCTACAGGTATTTTTAGGGAGGTTATTATGAAACAATACGGAGAAAAg      >  1:3075652/1‑91 (MQ=255)

                                                                                     ccTAGAGGTATTTGACAAGTTCCATCCAACTGTTTAAAATACCCCTACAGGTATTTTTAGGGAGGTTATTATGAAACAATACGGAGAAAAg      >  1:62571/1‑91 (MQ=255)

                                                                                     ccTAGAGGTATTTGACAAGTTCCATCCAACTGTTTAAAATACCCCTACAGGTATTTTTAGGGAGGTTATTATGAAACAATACGGAGAAAAg      >  1:2216178/1‑91 (MQ=255)

                                                                                      cTAGGGGTATTTGACAAGTTCCATCCAACTGTTTAAAATACCCCTACAGGTATTTTTAGGGAGGTTATTATGAAACAATACGGAGAAAAGt     >  1:1554442/1‑91 (MQ=255)

                                                                                       tAGAGGTATTTGACAAGTTCCATCCAACTGTTTAAAATACCCCTACAGGTATTTTTAGGGAGGTTATTATGAAACAATACGGAGAAAAGtt    >  1:1066463/1‑91 (MQ=255)

                                                                                       tAGAGGTATTTGACAAGTTCCATCCAACTGTTTAAAATACCCCTACAGGTATTTTTAGGGAGGTTATTATGAAACAATACGGAGAAAAGtt    >  1:1067410/1‑91 (MQ=255)

                                                                                         gagGTATTTGACAAGTTCCATCCAACTGTTTAAAATACCCCTACAGGTATTTTTAGGGAGGTTATTATGAAACAATACGGAGAAAAGTTTa  <  1:456454/91‑1 (MQ=255)

TTTTTTTATCATAATTCATAGTCATACCTCCACTTTTAATTGAATAAAAATATATTAATAGATAAACACAAATGTGTCAAATACCCCTAGAGGTATTTGACAAGTTCCATCCAACTGTTTAAAATACCCCTACAGGTATTTTTAGGGAGGTTATTATGAAACAATACGGAGAAAAGTTTA  >  USA300TCH1516_ALE/92897‑93076

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