BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	USA300TCH1516_ALE	92,987	A→G	15.3%	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.3%	474.7 / 78.2	204	intergenic (‑72/‑65)	ricR_1/glpE	Copper‑sensing transcriptional repressor RicR/Thiosulfate sulfurtransferase GlpE
Reads supporting (aligned to +/- strand): ref base A (83/89); new base G (29/2); total (113/91)
Fisher's exact test for biased strand distribution p-value = 8.44e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.91e-01

Only 100 of 218 total aligned reads displayed.

TTTTTTTATCATAATTCATAGTCATACCTCCACTTTTAATTGAATAAAAATATATTAATAGATAAACACAAATGTGTCAAATACCCCTAGAGGTATTTGACAAGTTCCATCCAACTGTTTAAAATACCCCTACAGGTATTTTTAGGGAGGTTATTATGAAACAATACGGAGAAAAGTTTAT  >  USA300TCH1516_ALE/92897‑93077

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

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

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

   ttttATCATAATTCATAGTCATACCTCCACTTTTAATTGAATAAAAATATATTAATAGATAAACACAAATGTGTCAAATACCCCTAGAGGt                                                                                         <  1:1453149/91‑1 (MQ=255)

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

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

       atcatcATTCATGGTCATACCTCCACTTTTAATTGAATAAAAATATATTAATAGATAAACACAAATGTGTCAAATACCTCTAGGGGTAttt                                                                                     <  1:3177197/91‑1 (MQ=39)

       aTCGTAATTCATAGTCATACCTCCACTTTTAATTGAATAAAAATATATTAATAGATAAACACAAATGTGTCAAATACCTCTAGGGGTAttt                                                                                     <  1:3496622/91‑1 (MQ=39)

       aTCATAATTCATAGTCATACCTCCACTTTTAATTGAATAAAAATATATTAATAGATAAACACAAATGTGTCAAATACCCCTAGAGGTAttt                                                                                     <  1:1600646/91‑1 (MQ=255)

        tCATAATTCATAGTCATACCTCCACTTTTAATTGAATAAAAATATATTAATAGATAAACACAAATGTGTCAAATACCCCTAGAGGTATTTg                                                                                    <  1:2279078/91‑1 (MQ=255)

        tCATAATTCATAGTCATACCTCCACTTTTAATTGAATAAAAATATATTAATAGATAAACACAAATGTGTCAAATACCCCTAGAGGTATTTg                                                                                    <  1:3726302/91‑1 (MQ=255)

          aTAATTCATAGTCATACCTCCACTTTTAATTGAATAAAAATATATTAATAGATAAACACAAATGTGTCAAATACCCCTAGAGGTATTTGAc                                                                                  <  2:522801/91‑1 (MQ=255)

          aTAATTCATAGTCATACCTCCACTTTTAATTGAATAAAAATATATTAATAGATAAACACAAATGTGTCAAATACCCCTAGAGGTATTTGAc                                                                                  >  1:2516605/1‑91 (MQ=255)

            aaTTCATAGTCATACCTCCACTTTTAATTGAATAAAAATATATTAATAGATAAACACAAATGTGTCAAATACCCCTAGAGGTATTTGACaa                                                                                <  1:2505730/91‑1 (MQ=255)

            aaTTCATAGTCATACCTCCACTTTTAATTGAATAAAAATATATTAATAGATAAACACAAATGTGTCAAATACCCCTAGAGGTATTTGACaa                                                                                >  1:1684227/1‑91 (MQ=255)

               tCATAGTCATACCTCCACTTTTAATTGAATAAAAATATATTAATAGATAAACACAAATGTGTCAAATACCCCTAGAGGTATTTGACAAGtt                                                                             <  1:526698/91‑1 (MQ=255)

                cATAGTCATACCTCCACTTTTAATTGAATAAAAATATATTAATAGATAAACACAAATGTGTCAAATACCCCTAGAGGTATTTGACAAGTtc                                                                            >  2:1478559/1‑91 (MQ=255)

                cATAGTCATACCTCCACTTTTAATTGAATAAAAATATATTAATAGATAAACACAAATGTGTCAAATACCCCTAGAGGTATTTGACAAGTtc                                                                            >  2:2994353/1‑91 (MQ=255)

                 aTAGTCATACCTCCACTTTTAATTGAATAAAAATATATTAATAGATAAACACAAATGTGTCAAATACCCCTAGAGGTATTTGACAAGTtcc                                                                           >  2:1052277/1‑91 (MQ=255)

                  tAGTCATACCTCCACTTTTAATTGAATAAAAATATATTAATAGATAAACACAAATGTGTCAAATACCCCTAGAGGTATTTGACAAGTtcca                                                                          <  2:735995/91‑1 (MQ=255)

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

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

                    gTCATACCTCCACTTTTAATTGAATAAAAATATATTAATAGATAAACACAAATGTGTCAAATACCCCTAGAGGTATTTGACAAGTtccatc                                                                        >  1:3202974/1‑91 (MQ=255)

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

                       aTACCTCCACTTTTAATTTTTTTTTAATATATTAATAGATAAACACAAATGTGTCAAATACCCATAGAGGTATTTGAGAAGTTCCATCCaa                                                                     >  2:3737525/1‑91 (MQ=255)

                        tACCTCCACTTTTAATTGAATAAAAATATATTAATAGATAAACACAAATGTGTCAAATACCCCTAGAGGTATTTGACAAGTTCCATCCAAc                                                                    <  2:1941228/91‑1 (MQ=255)

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

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

                         aCCTCCACTTTTAATTGAATAAAAATATATTAATAGATAAACACAAATGTGTCAAATACCCCTAGAGGTATTTGACAAGTTCCATCCAACt                                                                   <  1:2404160/91‑1 (MQ=255)

                            tCCACTTTTAATTGAATAAAAATATATTAATAGATAAACACAAATGTGTCAAATACCCCTAGAGGTATTTGACAAGTTCCATCCAACTGtt                                                                <  2:1549924/91‑1 (MQ=255)

                              cACTTTTAATTGAATAAAAATATATTAATAGATAAACACAAATGTGTCAAATACCCCTAGAGGTATTTGACAAGTTCCATCCAACTGTTTa                                                              >  2:3294221/1‑91 (MQ=255)

                              cACTTTTAATTGAATAAAAATATATTAATAGATAAACACAAATGTGTCAAATACCCCTAGAGGTATTTGACAAGTTCCATCCAACTGTTTa                                                              >  1:332664/1‑91 (MQ=255)

                              cACTTTTAATTGAATAAAAATATATTAATAGATAAACACAAATGTGTCAAATACCCCTAGAGGTATTTGACAAGTTCCATCCAACTGTTTa                                                              >  1:2207732/1‑91 (MQ=255)

                                cTTTTAATTGAATAAAAATATATTAATAGATAAACACAAATGTGTCAAATACCCCTAGAGGTATTTGACAAGTTCCATCCAACTGTTTaaa                                                            >  2:2344169/1‑91 (MQ=255)

                                 ttttAATTGAATAAAAATATATTAATAGATAAACACAAATGTGTCAAATACCCCTAGAGGTATTTGACAAGTTCCATCCAACTGTTTaaaa                                                           >  1:1429916/1‑91 (MQ=255)

                                  tttAATTGAATAAAAATATATTAATAGATAAACACAAATGTGTCAAATACCCCTAGAGGTATTTGACAAGTTCCATCCAACTGTTTAAAAt                                                          >  1:1034693/1‑91 (MQ=255)

                                  tttAATTGAATAAAAATATATTAATAGATAAACACAAATGTGTCAAATACCCCTAGAGGTATTTGACAAGTTCCATCCAACTGTTTAAAAt                                                          >  2:3586344/1‑91 (MQ=255)

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

                                   ttAATTGAATAAAAATATATTAATAGATAAACACAAATGTGTCAAATACCCCTAGAGGTATTTGACAAGTTCCATCCAACTGTTTAAAATa                                                         >  2:2160298/1‑91 (MQ=255)

                                   ttAATTGAATAAAAATATATTAATAGATAAACACAAATGTGTCAAATACCCCTAGAGGTATTTGACAAGTTCCATCCAACTGTTTAAAATa                                                         <  1:3257728/91‑1 (MQ=255)

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

                                    tAATTGAATAAAAATATATTAATAGATAAACACAAATGTGTCAAATACCCCTAGAGGTATTTGACAAGTTCCATCCAACTGTTTAAAATAc                                                        <  2:62/91‑1 (MQ=255)

                                     aaTTGAATAAAAATATATTAATAGATAAACACAAATGTGTCAAATACCCCTAGAGGTATTTGACAAGTTCCATCCAACTGTTTAAAATAcc                                                       <  1:1916561/91‑1 (MQ=255)

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

                                       ttGAATAAAAATATATTAATAGATAAACACAAATGTGTCAAATACCCCTAGAGGTATTTGACAAGTTCCATCCAACTGTTTAAAATAcccc                                                     <  1:1346480/91‑1 (MQ=255)

                                        tGAATAAAAATATATTAATAGATAAACACAAATGTGTCAAATACCTCTAGAGGTATTTGACAAGTTCCATCCAACTGTTTAAAATACTCCt                                                    >  1:1377856/1‑91 (MQ=255)

                                        tGAATAAAAATATATTAATAGATAAACACAAATGTGTCAAATACCCCTAGAGGTATTTGACAAGTTCCATCCAACTGTTTAAAATACCCCt                                                    >  1:1334581/1‑91 (MQ=255)

                                            tAAAAATATATTAATAGATAAACACAAATGTTTCAAATACCCCTAGAGGTATTTGACAAGTTCCATCCAACTGTTTAAAATACCCCTACAg                                                <  2:1195766/91‑1 (MQ=255)

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

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

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

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

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

                                                  tatatTAATAGATAAACACAAATGTGTCAAATACCCCTAGAGGTATTTGACAAGTTCCATCCAACTGTTTAAAATACCCCTACAGGTAttt                                          <  2:3823596/91‑1 (MQ=255)

                                                  tatatTAATAGATAAACACAAATGTGTCAAATACCCCTAGAGGGATTTGACAAGTTCCATCCAACTGTTTAAAATACCTGTAGGGGTAttt                                          <  2:408518/91‑1 (MQ=255)

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

                                                     atTAATAGATAAACACAAATGTGTCAAATACCCCTAGAGGTATTTGACAAGTTCCATCCAACTGTTTAAAATACCCCTACAGGTATTTTTa                                       <  2:1755468/91‑1 (MQ=255)

                                                     atTAATAGATAAACACAAATGTGTCAAATACCCCTAGAGGTATTTGACAAGTTCCATCCAACTGTTTAAAATACCCCTACAGGTATTTTTa                                       <  2:2375800/91‑1 (MQ=255)

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

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

                                                             aTAAACACAAATGTGTCAAATACCCCTAGAGGTATTTGACAAGTTCCATCCAACTGTTTAAAATACCCCTACAGGTATTTTTAGGGAGGtt                               <  2:3244803/91‑1 (MQ=255)

                                                             aTAAACACAAATGTGTCAAATACCCCTAGAGGTATTTGACAAGTTCCATCCAACTGTTTAAAATACCCCTACAGGTATTTTTAGGGAGGtt                               <  1:1248242/91‑1 (MQ=255)

                                                             aTAAACACAAATGTGTCAAATACCCCTAGAGGTATTTGACAAGTTCCATCCAACTGTTTAAAATACCCCTACAGGTATTTTTAGGGAGGtt                               <  1:2126460/91‑1 (MQ=255)

                                                             aTAAACACAAATGTGTCAAATACCCCTAGAGGTATTTGACAAGTTCCATCCAACTGTTTAAAATACCCCTACAGGTATTTTTAGGGAGGtt                               <  1:2483259/91‑1 (MQ=255)

                                                             aTAAACACAAATGTGTCAAATACCCCTAGAGGTATTTGACAAGTTCCATCCAACTGTTTAAAATACCCCTACAGGTATTTTTAGGGAGGtt                               <  1:3011374/91‑1 (MQ=255)

                                                             aTAAACACAAATGTGTCAAATACCCCTAGAGGTATTTGACAAGTTCCATCCAACTGTTTAAAATACCCCTACAGGTATTTTTAGGGAGGtt                               <  2:3352910/91‑1 (MQ=255)

                                                             aTAAACACAAATGTGTCAAATACCCCTAGAGGTATTTGACAAGTTCCATCCAACTGTTTAAAATACCCCTACAGGTATTTTTAGGGAGGtt                               >  1:714474/1‑91 (MQ=255)

                                                               aaaCACAAATGTGTCAAATACCCCTAGAGGTATTTGACAAGTTCCATCCAACTGTTTAAAATACCCCTACAGGTATTTTTAGGGAGGttat                             >  2:3059167/1‑91 (MQ=255)

                                                                aaCACAAATGTGTCAAATACCCCTAGAGGTATTTGACAAGTTCCATCCAACTGTTTAAAATACCCCTACAGGTATTTTTAGGGAGGttatt                            >  1:3568495/1‑91 (MQ=255)

                                                                      aaTGTGTCAAATACCCCTAGAGGTATTTGACAAGTTCCATCCAACTGTTTAAAATACCCCTACAGGTATTTTTCGGGAGGTTATTATGaaa                      >  2:3639480/1‑91 (MQ=255)

                                                                          tgtCAAATACCCCTAGAGGTATTTGACAAGTTCCATCCAACTGTTTAAAATACCCCTACAGGTATTTTTAGGGAGGTTATTATGAAACAAt                  >  2:545366/1‑91 (MQ=255)

                                                                           gtCAAATACCTCTAGGGGTATTTGACAAGTTCCATCCAACTGTTTAAAATACCCCTACAGGTATTTTTAGGGAGGTTATTATGAAACAATa                 >  2:1781202/1‑91 (MQ=255)

                                                                           gtCAAATACCCCTAGAGGTATTTGACAAGTTCCATCCAACTGTTTAAAATACCCCTACAGGTATTTTTAGGGAGGTTATTATGAAACAATa                 <  1:2946627/91‑1 (MQ=255)

                                                                           gtCAAATACCCCTAGAGGTATTTGACAAGTTCCAGCAAACAGTTTAACAGACCCCTACAGGTATTTGCAGGGAGGTTATTATGAAACAATa                 >  2:3605189/1‑91 (MQ=255)

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

                                                                            tCAAATACCCCTAGAGGTATTTGACAAGTTCCATCCAACTGTTTAAAATACCCCTACAGGTATTTTTAGGGAGGTTATTATGAAACAATAc                <  2:3006106/91‑1 (MQ=255)

                                                                              aaaTACCTCTAGGGGTATTTGACAAGTTCCATCCAACTGTTTAAAATACCCCTACAGGTATTTTTAGGGAGGTTATTATGAAACAATACgg              >  2:1973722/1‑91 (MQ=255)

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

                                                                                aTACCTCTAGGGGTATTTGACAAGTTCCATCCAACTGTTTAAAATACCCCTACAGGTATTTTTAGGGAGGTTATTATGAAACAATACGgag            >  2:1393400/1‑91 (MQ=255)

                                                                                aTACCTCTAGGGGTATTTGACAAGTTCCATCCAACTGTTTAAAATACCCCTACAGGTATTTTTAGGGAGGTTATTATGAAACAATACGgag            >  2:879112/1‑91 (MQ=255)

                                                                                aTACCCCTAGAGGTATTTGACAAGTTCCATCCAACTGTTTAAAATACCCCTACAGGTATTTTTAGGGAGGTTATTATGAAACAATACGgag            <  1:99683/91‑1 (MQ=255)

                                                                                 tACCTCTAGGGGTATTTGACAAGTTCCATCCAACTGTTTAAAATACCCCTACAGGTATTTTTAGGGAGGTTATTATGAAACAATACGgaga           >  2:2733261/1‑91 (MQ=255)

                                                                                 tACCTCTAGGGGTATTTGACAAGTTCCATCCAACTGTTTAAAATACCCCTACAGGTATTTTTAGGGAGGTTATTATGAAACAATACGgaga           >  2:75912/1‑91 (MQ=255)

                                                                                 tACCTCTAGGGGTATTTGACAAGTTCCATCCAACTGTTTAAAATACCCCTACAGGTATTTTTAGGGAGGTTATTATGAAACAATACGgaga           >  1:2778961/1‑91 (MQ=255)

                                                                                 tACCCCTAGAGGTATTTGACAAGTTCCATCCAACTGTTTAAAATACCCCTACAGGTATTTTTAGGGAGGTTATTATGAAACAATACGgaga           <  1:428686/91‑1 (MQ=255)

                                                                                 tACCCCTAGAGGTATTTGACAAGTTCCATCCAACTGTTTAAAATACCCCTACAGGTATTTTTAGGGAGGTTATTATGAAACAATACGgaga           <  2:1816892/91‑1 (MQ=255)

                                                                                 tACCCCTAGAGGTATTTGACAAGTTCCATCCAACTGTTTAAAATACCCCTACAGGTATTTTTAGGGAGGTTATTATGAAACAATACGgaga           <  1:1157207/91‑1 (MQ=255)

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

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

                                                                                    ctcTAGGGGTATTTGACAAGTTCCATCCAACTGTTTAAAATACCCCTACAGGTATTTTTAGGGAGGTTATTATGAAACAATACGGAGaaaa        >  2:315756/3‑91 (MQ=255)

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

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

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

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

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

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

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

                                                                                        agagGTATTTGACAAGTTCCATCCAACTGTTTAAAATACCCCTACAGGTATTTTTAGGGAGGTTATTATGAAACAATACGGAGAAAAGttt    >  1:1181127/1‑91 (MQ=255)

                                                                                         gagGTATTTGACAAGTTCCATCCAACTGTTTAAAATACCCCTACAGGTATTTTTAGGGAGGTTATTATGAAACAATACGGAGAAAAGTTTa   >  2:1435600/1‑91 (MQ=255)

                                                                                         gagGTATTTGACAAGTTCCATCCAACTGTTTAAAATACCCCTACAGGTATTTTTAGGGAGGTTATTATGAAACAATACGGAGAAAAGTTTa   >  1:3029183/1‑91 (MQ=255)

TTTTTTTATCATAATTCATAGTCATACCTCCACTTTTAATTGAATAAAAATATATTAATAGATAAACACAAATGTGTCAAATACCCCTAGAGGTATTTGACAAGTTCCATCCAACTGTTTAAAATACCCCTACAGGTATTTTTAGGGAGGTTATTATGAAACAATACGGAGAAAAGTTTAT  >  USA300TCH1516_ALE/92897‑93077

Alignment Legend
Aligned base mismatch/match (shaded by quality score): `ATCG/ATCG` `< 3 ≤` `ATCG/ATCG` `< 22 ≤` `ATCG/ATCG` `< 37 ≤` `ATCG/ATCG` `< 41 ≤` `ATCG/ATCG`
Unaligned base: `atcg` Masked matching base: `atcg` Alignment gap: `‑` Deleted base: `‑`