BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	USA300TCH1516_ALE	1,368,899	C→A	7.5%	intergenic (+303/+77)	USA300TCH1516_01277 → / ← lexA_1	hypothetical protein/LexA repressor

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	USA300TCH1516_ALE	1,368,899	0	C	A	7.5%	468.1 / 23.1	134	intergenic (+303/+77)	USA300TCH1516_01277/lexA_1	hypothetical protein/LexA repressor
Reads supporting (aligned to +/- strand): ref base C (70/54); new base A (2/8); total (72/62)
Fisher's exact test for biased strand distribution p-value = 4.39e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.81e-01

Only 100 of 137 total aligned reads displayed.

AAAATTTATATTAAAGAATGAGATTATCTAACAATTCGCAATTTCTGTTTGAACTGAATTTAGATTAGAATAAATAAAACCTCTTGAAAACGTAACTTTTTCAAGAGGTCCAAAATTTATAAAAATGAGGGCTTAACAATTACGATAATATATAACGGTTTAAAATATTACATTTCGCGGT  >  USA300TCH1516_ALE/1368809‑1368989

aaaaTTTATATTAAAGAATGAGATTATCTAACAATTCGCAATTTCTGTTTGAACTGAATTTAGATTAGAATAAATAAAACCTCTTGAAAAc                                                                                            >  1:574437/1‑91 (MQ=255)

 aaaTTTATATTAAAGAATGAGATTATCTAACAATTCGCAATTTCTGTTTGAACTGAATTTAGATTAGAATAAATAAAACCTCTTGAAAACg                                                                                           <  1:753866/91‑1 (MQ=255)

    tttATATTAAAGAATGAGATTATCTAACAATTCGCAATTTCTGTTTGAACTGAATTTAGATTAGAATAAATAAAACCTCTTGAAAACGTaa                                                                                        <  1:2880396/91‑1 (MQ=255)

     ttATATTAAAGAATGAGATTATCTAACAATTCGCAATTTCTGTTTGAACTGAATTTAGATTAGAATAAATAAAACCTCTTGAAAACGTAAc                                                                                       <  1:3926605/91‑1 (MQ=255)

     ttATATTAAAGAATGAGATTATCTAACAATTCGCAATTTCTGTTTGAACTGAATTTAGATTAGAATAAATAAAACCTCTTGAAAACGTAAc                                                                                       <  1:2886586/91‑1 (MQ=255)

        tatTAAAGAATGAGATTATCTAACAATTCGCAATTTCTGTTTGAACTGAATTTAGATTAGAATAAATAAAACCTCTTGAAAACGTAACttt                                                                                    >  1:3640370/1‑91 (MQ=255)

         atTAAAGAATGAGATTATCTAACAATTCGCAATTTCTGTTTGAACTGAATTTAGATTAGAATAAATAAAACCTCTTGAAAACGTAACtttt                                                                                   >  1:2030399/1‑91 (MQ=255)

         atTAAAGAATGAGATTATCTAACAATTCGCAATTTCTGTTTGAACTGAATTTAGATTAGAATAAATAAAACCTCTTGAAAACGTAACtttt                                                                                   >  1:828929/1‑91 (MQ=255)

          ttAAAGAATGAGATTATCTAACAATTCGCAATTTCTGTTTGAACTGAATTTAGATTAGAATAAATAAAACCTCTTGAAAACGTAACttttt                                                                                  >  1:3167382/1‑91 (MQ=255)

           tAAAGAATGAGATTATCTAACAATTCGCAATTTCTGTTTGAACTGAATTTAGATTAGAATAAATAAAACCTCTTGAAAACGTAACTTTTTc                                                                                 <  1:2711839/91‑1 (MQ=255)

           tAAAGAATGAGATTATCTAACAATTCGCAATTTCTGTTTGAACTGAATTTAGATTAGAATAAATAAAACCTCTTGAAAAAGTTACGTTTTc                                                                                 <  1:601476/91‑1 (MQ=39)

            aaaGAATGAGATTATCTAACAATTCGCAATTTCTGTTTGAACTGAATTTAGATTAGAATAAATAAAACCTCTTGAAAAAGTTACGTTTTCa                                                                                <  1:2920782/91‑1 (MQ=39)

            aaaGAATGAGATTATCTAACAATTCGCAATTTCTGTTTGAACTGAATTTAGATTAGAATAAATAAAACCTCTTGAAAAAGTTACGTTTTCa                                                                                <  1:3711016/91‑1 (MQ=39)

            aaaGAATGAGATTATCTAACAATTCGCAATTTCTGTTTGAACTGAATTTAGATTAGAATAAATAAAACCTCTTGAAAAAGTTACGTTTTCa                                                                                <  1:3710260/91‑1 (MQ=39)

            aaaGAATGAGATTATCTAACAATTCGCAATTTCTGTTTGAACTGAATTTAGATTAGAATAAATAAAACCTCTTGAAAAAGTTACGTTTTCa                                                                                <  1:3710150/91‑1 (MQ=39)

              aGAATGAGATTATCTAACAATTCGCAATTTCTGTTTGAACTGAATTTAGATTAGAATAAATAAAACCTCTTGAAAAAGTTACGTTTTCAag                                                                              <  1:53575/91‑1 (MQ=38)

                 aTGAGATTATCTAACAATTCGCAATTTCTGTTTGAACTGAATTTAGATTAGAATAAATAAAACCTCTTGAAAACGTAACTTTTTCAAGAgg                                                                           <  1:3819836/91‑1 (MQ=255)

                 aTGAGATTATCTAACAATTCGCAATTTCTGTTTGAACTGAATTTAGATTAGAATAAATAAAACCTCTTGAAAACGTAACTTTTTCAAGAgg                                                                           <  1:2151677/91‑1 (MQ=255)

                  tGAGATTATCTAACAATTCGCAATTTCTGTTTGAACTGAATTTAGATTAGAATAAATAAAACCTCTTGAAAAAGTTACGTTTTCAAGAGGt                                                                          <  1:2189709/91‑1 (MQ=38)

                  tGAGATTATCTAACAATTCGCAATTTCTGTGTGAACTGAATTTAGATTAGAATAAATAAAACCTCTTGAAAACGTAACTTTTTCAAGAGGt                                                                          <  1:1539200/91‑1 (MQ=255)

                   gagaTTATCTAACAATTCGCAATTTCTGTTTGAACTGAATTTAGATTAGAATAAATAAAACCTCTTGAAAACGTAACTTTTTCAAGAGGTc                                                                         >  1:4022615/1‑91 (MQ=255)

                   gagaTTATCTAACAATTCGCAATTTCTGTTTGAACTGAATTTAGATTAGAATAAATAAAACCTCTTGAAAACGTAACTTTTTCAAGAGGTc                                                                         >  1:2638778/1‑91 (MQ=255)

                   gagaTTATCTAACAATTCGCAATTTCTGTTTGAACTGAATTTAGATTAGAATAAATAAAACCTCTTGAAAACGTAACTTTTTCAAGAGGTc                                                                         >  1:2341031/1‑91 (MQ=255)

                    agaTTATCTAACAATTCGCAATTTCTGTTTGAACTGAATTTAGATTAGAATAAATAAAACCTCTTGAAAACGTAACTTTTTCAAGAGGTcc                                                                        >  1:218400/1‑91 (MQ=255)

                    agaTTATCTAACAATTCGCAATTTCTGTTTGAACTGAATTTAGATTAGAATAAATAAAACCTCTTGAAAACGTAACTTTTTCAAGAGGTcc                                                                        >  1:218483/1‑91 (MQ=255)

                    agaTTATCTAACAATTCGCAATTTCTGTTTGAACTGAATTTAGATTAGAATAAATAAAACCTCTTGAAAACGTAACTTTTTCAAGAGGTcc                                                                        >  1:1674274/1‑91 (MQ=255)

                    agaTTATCTAACAATTCGCAATTTCTGTTTGAACTGAATTTAGATTAGAATAAATAAAACCTCTTGAAAACGTAACTTTTTCAAGAGGTcc                                                                        >  1:516127/1‑91 (MQ=255)

                    agaTTATCTAACAATTCGCAATTTCTGTTTGAACTGAATTTAGATTAGAATAAATAAAACCTCTTGAAAACGTAACTTTTTCAAGAGGTcc                                                                        >  1:266982/1‑91 (MQ=255)

                    agaTTATCTAACAATTCGCAATTTCTGTTTGAACTGAATTTAGATTAGAATAAATAAAACCTCTTGAAAACGTAACTTTTTCAAGAGGTcc                                                                        >  1:2548663/1‑91 (MQ=255)

                      aTTATCTAACAATTCGCAATTTCTGTTTGAACTGAATTTAGATTAGAATAAATAAAACCTCTTGAAAACGTAACTTTTTCAAGAGGTCCaa                                                                      >  1:3763943/1‑91 (MQ=255)

                      aTTATCTAACAATTCGCAATTTCTGTTTGAACTGAATTTAGATTAGAATAAATAAAACCTCTTGAAAACGTAACTTTTTCAAGAGGTCCaa                                                                      >  1:3763874/1‑91 (MQ=255)

                             aaCAATTCGCAATTTCTGTTTGAACTGAATTTAGATTAGAATAAATAAAACCTCTTGAAAACGTAACTTTTTCAAGAGGTCCAAAATTtat                                                               >  1:2130909/1‑91 (MQ=255)

                              aCAATTCGCAATTTCTGTTTGAACTGAATTTAGATTAGAATAAATAAAACCTCTTGAAAACGTAACTTTTTCAAGAGGTCCAAAATTtata                                                              >  1:2947780/1‑91 (MQ=255)

                              aCAATTCGCAATTTCTGTTTGAACTGAATTTAGATTAGAATAAATAAAACCTCTTGAAAACGTAACTTTTTCAAGAGGTCCAAAATTtata                                                              >  1:3663573/1‑91 (MQ=255)

                               cAATTCGCAATTTCTGTTTGAACTGAATTTAGATTAGAATAAATAAAACCTCTTGAAAACGTAACTTTTTCAAGAGGTCCAAAATTTATaa                                                             <  1:25817/91‑1 (MQ=255)

                               cAATTCGCAATTTCTGTTTGAACTGAATTTAGATTAGAATAAATAAAACCTCTTGAAAACGTAACTTTTTCAAGAGGTCCAAAATTTATaa                                                             >  1:293896/1‑91 (MQ=255)

                               cAATTCGCAATTTCTGTTTGAACTGAATTTAGATTAGAATAAATAAAACCTCTTGAAAACGTAACTTTTTCAAGAGGTCCAAAATTTATaa                                                             <  1:25460/91‑1 (MQ=255)

                                aaTTCGCAATTTCTGTTTGAACTGAATTTAGATTAGAATAAATAAAACCTCTTGAAAACGTAACTTTTTCAAGAGGTCCAAAATTTATaaa                                                            >  1:917700/1‑91 (MQ=255)

                                 aTTCGCAATTTCTGTTTGAACTGAATTTAGATTAGAATAAATAAAACCTCTTGAAAACGTAACTTTTTCAAGAGGTCCAAAATTTATaaaa                                                           <  1:1819205/91‑1 (MQ=255)

                                 aTTCGCAATTTCTGTTTGAACTGAATTTAGATTAGAATAAATAAAACCTCTTGAAAACGTAACTTTTTCAAGAGGTCCAAAATTTATaaaa                                                           <  1:3704805/91‑1 (MQ=255)

                                  ttCGCAATTTCTGTTTGAACTGAATTTAGATTAGAATAAATAAAACCTCTTGAAAACGTAACTTTTTCAAGAGGTCCAAAATTTATaaaaa                                                          <  1:2833164/91‑1 (MQ=255)

                                    cGCAATTTCTGTTTGAACTGAATTTAGATTAGAATAAATAAAACCTCTTGAAAACGTAACTTTTTCAAGAGGTCCAAAATTTATAAAAATg                                                        <  1:3356870/91‑1 (MQ=255)

                                     gCAATTTCTGTTTGAACTGAATTTAGATTAGAATAAATAAAACCTCTTGAAAACGTAACTTTTTCAAGAGGTCC‑AAATTTATAAAAAagag                                                      >  1:3351076/1‑91 (MQ=255)

                                      cAATTTCTGTTTGAACTGAATTTAGATTAGAATAAATAAAACCTCTTGAAAACGTAACTTTTTCAAGAGGTCCAAAATTTATAAAAATGAg                                                      >  1:1739400/1‑91 (MQ=255)

                                      cAATTTCTGTTTGAACTGAATTTAGATTAGAATAAATAAAACCTCTTGAAAACGTAACTTTTTCAAGAGGTCCAAAATTTATAAAAATGAg                                                      >  1:898763/1‑91 (MQ=255)

                                      cAATTTCTGTTTGAACTGAATTTAGATTAGAATAAATAAAACCTCTTGAAAACGTAACTTTTTCAAGAGGTCCAAAATTTATAAAAATGAg                                                      >  1:2876103/1‑91 (MQ=255)

                                      cAATTTCTGTTTGAACTGAATTTAGATTAGAATAAATAAAACCTCTTGAAAACGTAACTTTTTCAAGAGGTCCAAAATTTATAAAAATGAg                                                      <  1:2361855/91‑1 (MQ=255)

                                      cAATTTCTGTTTGAACTGAATTTAGATTAGAATAAATAAAACCTCTTGAAAACGTAACTTTTTCAAGAGGTCCAAAATTTATAAAAATGAg                                                      <  1:2361748/91‑1 (MQ=255)

                                      cAATTTCTGTTTGAACTGAATTTAGATTAGAATAAATAAAACCTCTTGAAAACGTAACTTTTTCAAGAGGTCCAAAATTTATAAAAATGAg                                                      >  1:700727/1‑91 (MQ=255)

                                       aaTTTCTGTTTGAACTGAATTTAGATTAGAATAAATAAAACCTCTTGAAAACGTAACTTTTTCAAGAGGTCCAAAATTTATAAAAATGAgg                                                     >  1:1061998/1‑91 (MQ=255)

                                       aaTTTCTGTTTGAACTGAATTTAGATTAGAATAAATAAAACCTCTTGAAAACGTAACTTTTTCAAGAGGTCCAAAATTTATAAAAATGAgg                                                     >  1:3627403/1‑91 (MQ=255)

                                       aaTTTCTGTTTGAACTGAATTTAGATTAGAATAAATAAAACCTCTTGAAAACGTAACTTTTTCAAGAGGTCCAAAATTTATAAAAATGAgg                                                     >  1:1742902/1‑91 (MQ=255)

                                       aaTTTCTGTTTGAACTGAATTTAGATTAGAATAAATAAAACCTCTTGAAAACGTAACTTTTTCAAGAGGTCCAAAATTTATAAAAATGAgg                                                     >  1:2573247/1‑91 (MQ=255)

                                       aaTTTCTGTTTGAACTGAATTTAGATTAGAATAAATAAAACCTCTTGAAAACGTAACTTTTTCAAGAGGTCCAAAATTTATAAAAATGAgg                                                     >  1:1743216/1‑91 (MQ=255)

                                         tttCTGTTTGAACTGAATTTAGATTAGAATAAATAAAACCTCTTGAAAACGTAACTTTTTCAAGAGGTCCAAAATTTATAAAAATGAGGGc                                                   >  1:92256/1‑91 (MQ=255)

                                         tttCTGTTTGAACTGAATTTAGATTAGAATAAATAAAACCTCTTGAAAACGTAACTTTTTCAAGAGGTCCAAAATTTATAAAAATGAGGGc                                                   <  1:1241386/91‑1 (MQ=255)

                                          ttCTGTTTGAACTGAATTTAGATTAGAATAAATAAAACCTCTTGAAAACGTAACTTTTTCAAGAGGTCCAAAATTTATAAAAATGAGGGCt                                                  <  1:4066479/91‑1 (MQ=255)

                                          ttCTGTTTGAACTGAATTTAGATTAGAATAAATAAAACCCCTTGAAAACGTAACTTTTTCAAGAGGTCCAAAATTTATAAAAATGAGGGCt                                                  <  1:4066406/91‑1 (MQ=255)

                                           tCTGTTTGAACTGAATTTAGATTAGAATAAATAAAACCTCTTGAAAACGTAACTTTTTCAAGAGGTCCAAAATTTATAAAAATGAGGGCtt                                                 >  1:3478270/1‑91 (MQ=255)

                                               tttGAACTGAATTTAGATTAGAATAAATAAAACCTCTTGAAAACGTAACTTTTTCAAGAGGTCCAAAATTTATAAAAATGAGGGCTTAACa                                             >  1:939530/1‑91 (MQ=255)

                                                 tGAACTGAATTTAGATTAGAATAAATAAAACCTCTTGAAAACGTAACTTTTTCAAGAGGTCCAAAATTTATAAAAATGAGGGCTTAACAAt                                           <  1:1677455/91‑1 (MQ=255)

                                                 tGAACTGAATTTAGATTAGAATAAATAAAACCTCTTGAAAACGTAACTTTTTCAAGAGGTCCAAAATTTATAAAAATGAGGGCTTAACAAt                                           >  1:2121555/1‑91 (MQ=255)

                                                  gAACTGAATTTAGATTAGAATAAATAAAACCTCTTGAAAACGTAACTTTTTCAAGAGGTCCAAAATTTATAAAAATGAGGGCTTAACAAtt                                          >  1:19093/1‑91 (MQ=255)

                                                    aCTGAATTTAGATTAGAATAAATAAAACCTCTTGAAAACGTAACTTTTTCAAGAGGTCCAAAATTTATAAAAATGAGGGCTTAACAATTAc                                        <  1:883362/91‑1 (MQ=255)

                                                    aCTGAATTTAGATTAGAATAAATAAAACCTCTTGAAAACGTAACTTTTTCAAGAGGTCCAAAATTGATAAAAATGAGGGCTTAACAATTAc                                        >  1:3865559/1‑91 (MQ=255)

                                                     cTGAATTTAGATTAGAATAAATAAAACCTCTTGAAAACGTAACTTTTTCAAGAGGTCCAAAATTTATAAAAATGAGGGCTTAACAATTACg                                       >  1:670795/1‑91 (MQ=255)

                                                       gAATTTAGATTAGAATAAATAAAACCTCTTGAAAACGTAACTTTTTCAAGAGGTCCAAAATTTATAAAAATGAGGGCTTAACAATTACGat                                     <  1:2523927/91‑1 (MQ=255)

                                                        aaTTTAGATTAGAATAAATAAAACCTCTTGAAAACGTAACTTTTTCAAGAGGTCCAAAATTTATAAAAATGAGGGCTTAACAATTACGata                                    <  1:2966188/91‑1 (MQ=255)

                                                        aaTTTAGATTAGAATAAATAAAACCTCTTGAAAACGTAACTTTTTCAAGAGGTCCAAAATTTATAAAAATGAGGGCTTAACAATTACGata                                    >  1:491360/1‑91 (MQ=255)

                                                         aTTTAGATTAGAATAAATAAAACCTCTTGAAAACGTAACTTTTTCAAGAGGTCCAAAATTTATAAAAATGAGGGCTTAACAATTACGataa                                   <  1:2899309/91‑1 (MQ=255)

                                                             agattagaATAAATAAAACCTCTTGAAAACGTAACTTTTTCAAGAGGTCCAAAATTTATAAAAATGAGGGCTTAACAATTACGATAatata                               <  1:1387788/91‑1 (MQ=255)

                                                              gattagaATAAATAAAACCTCTTGAAAACGTAACTTTTTCAAGAGGTCCAAAATTTATAAAAATGAGGGCTTAACAATTACGATAatatat                              >  1:4056482/1‑91 (MQ=255)

                                                               attagaATAAATAAAACCTCTTGAAAACGTAACTTTTTCAAGAGGTCCAAAATTTATAAAAATGAGGGCTTAACAATTACGATAAtatata                             <  1:1100596/91‑1 (MQ=255)

                                                                ttagaATAAATAAAACCTCTTGAAAACGTAACTTTTTCAAGAGGTCCAAAATTTATAAAAATGAGGGCTTAACAATTACGATAATATATaa                            >  1:3760051/1‑91 (MQ=255)

                                                                   gaATAAATAAAACCTCTTGAAAACGTAACTTTTTCAAGAGGTCCAAAATTTATAAAAATGAGGGCTTAACAATTACGATAATATATAACgg                         <  1:2996652/91‑1 (MQ=255)

                                                                      taaataaaACCTCTTGAAAACGTAACTTTTTCAAGAGGTCCAAAATTTATAAAAATGAGGGCTTAACAATTACGATAATATATAACGGttt                      <  1:854126/91‑1 (MQ=255)

                                                                      taaataaaACCTCTTGAAAACGTAACTTTTTCAAGAGGTCCAAAATTTATAAAAATGAGGGCTTAACAATTACGATAATATATAACGGttt                      <  1:2766728/91‑1 (MQ=255)

                                                                       aaataaaACCTCTTGAAAACGTAACTTTTTCAAGAGGTCCAAAATTTATAAAAATGAGGGCTTAACAATTACGATAATATATAACGGTTTa                     >  1:1692383/1‑91 (MQ=255)

                                                                        aataaaACCTCTTGAAAACGTAACTTTTTCAAGAGGTCCAAAATTTATAAAAATGAGGGCTTAACAATTACGATAATATATAACGGTTTaa                    >  1:1279853/1‑91 (MQ=255)

                                                                        aataaaACCTCTTGAAAACGTAACTTTTTCAAGAGGTCCAAAATTTATAAAAATGAGGGCTTAACAATTACGATAATATATAACGGTTTaa                    >  1:1504973/1‑91 (MQ=255)

                                                                        aataaaACCTCTTGAAAACGTAACTTTTTCAAGAGGTCCAAAATTTATAAAAATGAGGGCTTAACAATTACGATAATATATAACGGTTTaa                    >  1:1502759/1‑91 (MQ=255)

                                                                           aaaaCCTCTTGAAAACGTAACTTTTTCAAGAGGTCCACAATTTATAAAAATGAGGGCTTAACAATTACGATAATATATAACGGTTTAAAat                 >  1:3508702/1‑91 (MQ=255)

                                                                           aaaaCCTCTTGAAAACGTAACTTTTTCAAGAGGTCCAAAATTTATAAAAATGAGGGCTTAACAATTACGATAATATATAACGGTTTAAAat                 >  1:4032743/1‑91 (MQ=255)

                                                                           aaaaCCTCTTGAAAACGTAACTTTTTCAAGAGGTCCAAAATTTATAAAAATGAGGGCTTAACAATTACGATAATATATAACGGTTTAAAat                 <  1:3524767/91‑1 (MQ=255)

                                                                            aaaCCTCTTGAAAACGTAACTTTTTCAAGAGGTCCAAAATTTATAAAAATGAGGGCTTAACAATTACGATAATATATAACGGTTTAAAata                <  1:1428102/91‑1 (MQ=255)

                                                                             gaCCTCTTGAAAAAGTTACGTTTTCAAGAGGTCCAAAATTTATAAAAATGAGGGCTTAACAATTACGATAATATATAACGGTTTAAAatat               >  1:1070025/2‑91 (MQ=37)

                                                                             aaCCTCTTGAAAACGTAACTTTTTCAAGAGGTCCAAAATTTATAAAAATGAGGGCTTAACAATTACGATAATATATAACGGTTTAAAatat               <  1:1309822/91‑1 (MQ=255)

                                                                              aCCTCTTGAAAACGTAACTTTTTCAAGAGGTCCAAAATTTATAAAAATGAGGGCTTAACAATTACGATAATATATAACGGTTTAAAATAtt              >  1:1478986/1‑91 (MQ=255)

                                                                               ccTCTTGAAAACGTAACTTTTTCAAGAGGTCCAAAATTTATAAAAATGAGGGCTTAACAATTACGATAATATATAACGGTTTAAAATATTa             >  1:462050/1‑91 (MQ=255)

                                                                               ccTCTTGAAAACGTAACTTTTTCAAGAGGTCCAAAATTTATAAAAATGAGGGCTTAACAATTACGATAATATATAACGGTTTAAAATATTa             <  1:1153700/91‑1 (MQ=255)

                                                                               ccTCTTGAAAACGTAACTTTTTCAAGAGGTCCAAAATTTATAAAAATGAGGGCTTAACAATTACGATAATATATAACGGTTTAAAATATTa             <  1:942451/91‑1 (MQ=255)

                                                                                  ctTGAAAACGTAACTTTTTCAAGAGGTCCAAAATTTATAAAAATGAGGGCTTAACAATTACGATAATATATAACGGTTTAAAATATTACAt          >  1:665623/1‑91 (MQ=255)

                                                                                  ctTGAAAACGTAACTTTTTCAAGAGGTCCAAAATTTATAAAAATGAGGGCTTAACAATTACGATAATATATAACGGTTTAAAATATTACAt          <  1:90668/91‑1 (MQ=255)

                                                                                   ttGAAAACGTAACTTTTTCAAGAGGTCCAAAATTTATAAAAATGAGGGCTTAACAATTACGATAATATATAACGGTTTAAAATATTACAtt         >  1:793525/1‑91 (MQ=255)

                                                                                   ttGAAAACGTAACTTTTTCAAGAGGTCCAAAATTTATAAAAATGAGGGCTTAACAATTACGATAATATATAACGGTTTAAAATATTACAtt         >  1:2849178/1‑91 (MQ=255)

                                                                                       aaaCGTAACTTTTTCAAGAGGTCCAAAATTTATAAAAATGAGGGCTTAACAATTACGATAATATATAACGGTTTAAAATATTACATTTcgc     >  1:1770612/1‑91 (MQ=255)

                                                                                       aaaCGTAACTTTTTCAAGAGGTCCAAAATTTATAAAAATGAGGGCTTAACAATTACGATAATATATAACGGTTTAAAATATTACATTTcgc     <  1:1271853/91‑1 (MQ=255)

                                                                                        aaCGTAACTTTTTCAAGAGGTCCAAAATTTATAAAAATGAGGGCTTAACAATTACGATAATATATAACGGTTTAAAATATTACATTTcgcg    <  1:3168308/91‑1 (MQ=255)

                                                                                          cGTAACTTTTTCAAGAGGTCCAAAATTTATAAAAATGAGGGCTTAACAATTACGATAATATATAACGGTTTAAAATATTACATTTCGCGGt  >  1:721758/1‑91 (MQ=255)

                                                                                          cGTAACTTTTTCAAGAGGTCCAAAATTTATAAAAATGAGGGCTTAACAATTACGATAATATATAACGGTTTAAAATATTACATTTCGCGGt  <  1:3468348/91‑1 (MQ=255)

AAAATTTATATTAAAGAATGAGATTATCTAACAATTCGCAATTTCTGTTTGAACTGAATTTAGATTAGAATAAATAAAACCTCTTGAAAACGTAACTTTTTCAAGAGGTCCAAAATTTATAAAAATGAGGGCTTAACAATTACGATAATATATAACGGTTTAAAATATTACATTTCGCGGT  >  USA300TCH1516_ALE/1368809‑1368989

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