BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	USA300TCH1516_ALE	492,848	A→T	5.4%	F62Y (TTC→TAC)	USA300HOU_RS02345 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	USA300TCH1516_ALE	492,848	0	A	T	5.4%	472.0 / 2.8	149	F62Y (TTC→TAC)	USA300HOU_RS02345	hypothetical protein
Reads supporting (aligned to +/- strand): ref base A (77/61); new base T (9/2); total (86/63)
Fisher's exact test for biased strand distribution p-value = 1.19e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 150 total aligned reads displayed.

AAGGATTGGTCGAATGAACTTCATGTTCTATGTTCGACCAATCCTTTTCTGAATTAGTATTCATG‑TTTTACTTTGCGAATGCTTGTAAATAATCTAGCACCGTTTGTTATTAAAGTAACAACTGCCACTGCTTTTTTGAACTTACGTGGTGACTTAATTGAAATGTAAAAGTCTAACACGGTTCCTACAGTGCTTAAACCTAATGAAAGATATACTAATTTTTTATTTTTAGCATGATATTTATAGCCATTGTAGCCGTCGACTATGAAACCTGCGACA  >  USA300TCH1516_ALE/492710‑492988

aaGGATTGGTCGAATGAACTTCATGTTCTATGTTCGACCAATCCTTTTCTGAATTAGTATTCATG‑TTTTACTTTGCGAATGCTTGTAAATAATCTAGCACCGTTTGTTATTAAAGTAACAACTGCCACTGCTTTTTTGAAc                                                                                                                                            <  2:615536/141‑1 (MQ=255)

   ggttGGTCGAATGAACTTCATGTTCTATGTTCGACCAATCCTTTTCTGAATTAGTATTCATG‑GTTTACTTTGCGAATGCTTGTAAATAATCTAGCACCGTTTGTTATTAAAGTAACAACTGCCACTGCTTTTTTGAACTTa                                                                                                                                         >  2:41554/3‑141 (MQ=255)

   gagaggacGAATGACCTTCATGTTCTTTGTTCAACCTATCCTCATCTGAATTAGTCTTCATG‑TTTTGCTTTGCGAATGCTTGTAAATAATCTACCACCGTTTGTTATTAAAGTAACAACTGCCACTGCTTTTTTGAACTTa                                                                                                                                         >  2:376136/8‑141 (MQ=255)

    aTTGGTCGAATGAACTTCATGTTCTATGTTCGACCAATCCTTTTCTGAATTAGTATTCATG‑TTTTACTTTGCGAATGCTTGTAAATAATCTAGCACCGTTTGTTATTAAAGTAACAACTGCCACTGCTTTTTTGAACTTAc                                                                                                                                        <  3:127735/141‑1 (MQ=255)

    aTTGGTCGAATGAACTTCATGTTCTATGTTCGACCAATCCTTTTCTGAATTAGTATTCATG‑TTTTACTTTGCGAATGCTTGTAAATAATCTAGCACCGTTTGTTATGAAAGTAACAACTGCCACTGCTTTTTTGAACTTAc                                                                                                                                        >  1:1104228/1‑141 (MQ=255)

       ggTCGAATGAACTTCATGTTCTATGTTCGACCAATCCTTTTCTGAATTAGTATTCATG‑TTTTACTTTGCGAATGCTTGTAAATAATCTAGCACCGTTTGTTATTAAAGTAACAACTGCCACTGCTTTTTTGAACTTACgtg                                                                                                                                     <  1:1037605/141‑1 (MQ=255)

       ggTCGAATGAACTTCATGTTCTATGTTCGACCAATCCTTTTCTGAATTAGTATTCATG‑TTTTACTTTGCGAATGCTTGTAAATAATCTAGCACCGTTTGTTATTAAAGTAACAACTGCCACTGCTTTTTTGAACTTACgtg                                                                                                                                     >  4:59208/1‑141 (MQ=255)

        gTCGAATGAACTTCATGTTCTATGTTCGACCAATCCTTTTCTGAATTAGTATTCATG‑TTTTACTTTGCGAATGCTTGTAAATAATCTAGCACCGTTTGTTATTAAAGTAACAACTGCCACTGCTTTTTTGAACTTACgtgg                                                                                                                                    >  1:1143469/1‑141 (MQ=255)

        gTCGAATGAACTTCATGTTCTATGTTCGACCAATCCTTTTCTGAATTAGTATTCATG‑TTTTACTTTGCGAATGCTTGTAAATAATCTAGCACCGTTTGTTATTAAAGTAACAACTGCCACTGCTTTTTTGAACTTACgtgg                                                                                                                                    >  3:157529/1‑141 (MQ=255)

         tCGAATGAACTTCATGTTCTATGTTCGACCAATCCTTTTCTGAATTAGTATTCATG‑TTTTACTTTGCGAATGCTTGTAAATAAGCTAGCACCGTGTGTTATTAAAGTAACAACTGCCACTGCTTTTTTGTACTTACgtggt                                                                                                                                   <  2:874237/141‑1 (MQ=255)

            aaTGAACTTCATGTTCTATGTTCGACCAATCCTTTTCTGAATTAGTATTCATG‑TTTTACTTTGCGAATGCTTGTAAATAATCTAGCACCGTTTGTTATTAAAGTAACAACTGCCACTGCTTTTTTGAACTTACGTGGTGAc                                                                                                                                <  2:119941/141‑1 (MQ=255)

              tGAACTTCATGTTCTATGTTCGACCAATCCTTTTCTGAATTAGTATTCATG‑TTTTACTTTGCGAATGCTTGTAAATAATCTAGCACCGTTTGTTATTAAAGTAACAACTGCCACTGCTTTTTTGAACTTACGTGGTGACtt                                                                                                                              <  1:684218/141‑1 (MQ=255)

              tGAACTTCATGTTCTATGTTCGACCAATCCTTTTCTGAATTAGTATTCATG‑TTTTACTTTGCGAATGCTTGTAAATAATCTAGCACCGTTTGTTATTAAAGTAACAACTGCCACTGCTTTTTTGAACTTACGTGGTGACtt                                                                                                                              <  2:879829/141‑1 (MQ=255)

                aaCTTCATGTTCTATGTTCGACCAATCCTTTTCTGAATTAGTATTCATG‑TTTTACTTTGCGAATGCTTGTAAATAATCTAGCACCGTTTGTTATTAAAGTAACAACTGCCACTGCTTTTTTGAACTTACGTGGTGACTTaa                                                                                                                            <  2:66998/141‑1 (MQ=255)

                 aCTTCATGTTCTATGTTCGACCAATCCTTTTCTGAATTAGTATTCATG‑TTTTACTTTGCGAATGCTTGTAAATAATCTAGCACCGTTTGTTATTAAAGTAACAACTGCCACTGCTTTTTTGAACTTACGTGGTGACTTAAt                                                                                                                           >  2:662532/1‑141 (MQ=255)

                  cTTCATGTTCTATGTTCGACCAATCCTTTTCTGAATTAGTATTCATG‑TTTTACTTTGCGAATGCTTGTAAATAATCTAGCACCGTTTGTTATTAAAGTAACAACTGACACTGCTTTTTTGAACTTACGTGGTGACTTAAtt                                                                                                                          >  1:1105548/1‑141 (MQ=255)

                       ttttCTATGTTCCACCAATCTTTTTCTGAATTAGTATTCATT‑TTTTACTTTGCGAATGCGTGTAAATAATCTAGCACCGTTTGTTATTAAAGTAACAACTGCCACTGCTTTTTTGAACTTACGTGGTGACTTAATTGAAAt                                                                                                                     <  4:458407/139‑1 (MQ=255)

                         ttCTATGTTCGACCAATCCTTTTCTGAATTAGTATTCATG‑TTTTACTTTGCGAATGCTTGTAAATAATCTAGCACCGTTTGTTATTAGAGTAACAACTGCCACTGCTTTTTTGATCTTACGTGGTGACTTAATTGAAATGt                                                                                                                   <  2:87648/141‑1 (MQ=255)

                         ttCTATGTTCGACCAATCCTTTTCTGAATTAGTATTCATG‑TTTTACTTTGCGAATGCTTGTAAATAATCTAGCACCGTTTGTTATTAAAGTAACAACTGCCACTGCTTTTTTTAACTTACGTGGTGACTTAATTGAAATGt                                                                                                                   <  2:1055067/141‑1 (MQ=255)

                           cTATGTTCGACCAATCCTTTTCTGAATTAGTATTCATG‑TTTTACTTTGCGAATGCTTGTAAATAATCTAGCACCGTTTGTTATTAAAGTAACAACTGCCACTGCTTTTTTGAACTTACGTGGTGACTTAATTGAAATGTaa                                                                                                                 >  2:139334/1‑141 (MQ=255)

                            tATGTTCGACCAATCCTTTTCTGAATTAGTATTCATG‑TTTTACTTTGCGAATGCTTGTAAATAATCTAGCACCGTTTGTTATTAAAGTAACAACTGCCACTGCTTTTTTGAACTTACGTGGTGACTTAATTGAAATGTaaa                                                                                                                >  3:89477/1‑141 (MQ=255)

                            tATGTTCGACCAATCCTTTTCTGAATTAGTATTCATG‑TTTTACTTTGCGAATGCTTGTAAATAATCTAGCACCGTTTGTTATTAAAGTAACAACTGCCACTGCTTTTTTGAACTTACGTGGTGACTTAATTGAAATGTaaa                                                                                                                <  1:932387/141‑1 (MQ=255)

                            tATGTTCGACCAATCCTTTTCTGAATCAGTATTCATG‑TTTTACTTTGCGAATGCTTGTAAATAATCTAGCACCGTTTGTTATTAAAGTAACAACTGCCACTGCTTTTTTGAACTTACGTGGTGACTTAATTGAAATGTaaa                                                                                                                <  2:942646/141‑1 (MQ=255)

                            tATGTTCGACCAATCCTTTTCGGAATTAGTATTCATG‑TTTTACTTTGCGAATGCTTGTAAATAATCTAGCACCGTTTGTTATTAAAGTAACAACTGCCACTGCTTTTTTGAACTTACGTGGTGACTTAATTGAAATGTaaa                                                                                                                <  1:720799/141‑1 (MQ=255)

                              tGTTCGACCAATCCTTTTCTGAATTAGTATTCATG‑TTTTACTTTGCGAATGCTTGTAAATAATCTAGCACCGTTTGTTATTAAAGTAACAACTGCCACTGCTTTTTTGAACTTACGTGGTGACTTAATTGATATGTaaaat                                                                                                              >  4:421647/1‑140 (MQ=255)

                              tGTTCGACCAATCCTTTTCTGAATTAGTATTCATG‑TTTTACTTTGCGAATGCTTGTAAATAATCTAGCACCGTTTGTTATTAAAGTAACAACTGCCACTGCTTTTTTGAACTTACGTGGTGACTTAATTGAAATGTAAAAg                                                                                                              <  1:1060212/141‑1 (MQ=255)

                               gTTCGACCAATCCTTTTCTGAATTAGTATTCATG‑TTTTACTTTGCGAATGCTTGTAAATAATCTAGCACCGTTTGTTATTAAAGTAACAACTGCCACTGCTTTTTTGAACTTACGTGGTGACTTAATTGAAATGTAAAAGt                                                                                                             >  4:417505/1‑141 (MQ=255)

                                  cGACCAATCCTTTTCTGAATTAGTATTCATG‑TTTTACTTTGCGAATGCTTGTAAATAATCTAGCACCGTTTGTTATTAAAGTAACAACTGCCACTGCTTTTTTGAACTTACGTGGTGACTTAATTGAAATGTAAAAGTCTa                                                                                                          >  2:37663/1‑141 (MQ=255)

                                  cGACCAATCCTTTTCTGAATTAGTATTCATG‑TTTTACTTTGCGAATGCTTGTAAATAATCTAGCACCGTTTGTTATTAAAGTAACAACTGCCACTGCTTTTTTGAACTTACGTGGTGACTTAATTGAAATGTAAAAGTCTa                                                                                                          >  2:755155/1‑141 (MQ=255)

                                     ccAATCCTTTTCTGAATTAGTATTCATG‑TTTTACTTTGCGAATGCTTGTAAATAATCTAGCACCGTTTGTTATTAAAGTAACAACTGCCACTGCTTTTTTGAACTTACGTGGTGACTTAATTGAAATGTAAAAGTCTAaca                                                                                                       >  3:187162/1‑141 (MQ=255)

                                       aaTCCTTTTCTGAATTAGTATTCATG‑TTTTACTTTGCGAATGCTTGTAAATAATCTAGCACCGTTTGTTATTAAAGTAACAACTGCCACTGCTTTTTTGAACTTACGTGGTGACTTAATTGAAATGTAAAAGTCTAACACg                                                                                                     <  1:372240/141‑1 (MQ=255)

                                         tCCTTTTCTGAATTAGTATTCATGTTTTTTTTTTTTTTTTTCTTGTAAATAATCTAGCACCGTTTGTTATTAAAGTAACAACTGCCACTGCTTTTTTTAACTTACGTGGTGACTTAATTGTATTGTAATAGTCTAACACgt                                                                                                    >  2:256699/1‑140 (MQ=255)

                                          ccTTTTCTGAATTAGTATTCATG‑TTTTACTTTGCGAATGCTTGTAAATAATCTAGCACCGTTTGTTATTAAAGTAACAACTGCCACTGCTTTTTTGAACTTACGTGGTGACTTAATTGAAATGTAAAAGTCTAACACGGtt                                                                                                  <  2:167133/141‑1 (MQ=255)

                                             tttCTGAATTAGTATTCATG‑TTTTACTTTGCGAATGCTTGTAAATAATCTAGCACCGTTTGTTATTAAAGTAACAACTGCCACTGCTTTTTTGAACTTACGTGGTGACTTAATTGAAATGTAAAAGTCTAACACGGTTCCt                                                                                               >  1:651382/1‑141 (MQ=255)

                                             tttCTGAATTAGTATTCATG‑TTTTACTTTGCGAATGCTTGTAAATAATCTAGCACCGTTTGTTATTAAAGTAACAACTGCCACTGCTTTTTTGAACTTACGTGGTGACTTAATTGAAATGTAAAAGTCTAACACGGTTCCt                                                                                               >  1:1161041/1‑141 (MQ=255)

                                              ttCTGAATTAGTATTCATG‑TTTTACTTTGCGAATGCTTGTAAATAATCTAGCTCCGTTTGTTATTTAAGTAACAATTGCCACTGCTTTGTTGATCTTACGTGTTGACTTAATTGAAATGTAAAAGTCTAACACGGTTCCTa                                                                                              <  4:396831/141‑1 (MQ=255)

                                              ttCTGAATTAGTATTCATG‑TTTTACTTTGCGAATGCTTGTAAATAATCTAGCACCGTTTGTTATTAAAGTAACAACTGCCACTGCTTTTTTGAACTTACGTGGTGACTTAATTGAAATGTAAAAGTCTAACACGGTTCCTa                                                                                              >  2:285814/1‑141 (MQ=255)

                                                cTGAATTAGTATTCATG‑TTTTACTTTGCGAATGCTTGTAAATAATCTAGCACCGTTTGTTATTAAAGTAACAACTGCCACTGCTTTTTTGAACTTACGTGGTGACTTAATTGAAATGTAAAAGTCTAACACGGTTCCTACa                                                                                            <  2:183531/141‑1 (MQ=255)

                                                 tGAATTAGTATTCATG‑TTTTACTTTGCGAATGCTTGTAAATAATCTAGCACCGTTTGTTATTAAAGTAACAACTGCCACTGCTTTTTTGAACTTACGTGGTGACTTAATTGAAATGTAAAAGTCTAACACGGTTCCTACAg                                                                                           <  2:5926/141‑1 (MQ=255)

                                                 tGAATTAGTATTCATG‑TTTTACTTTGCGAATGCTTGTAAATAATCTAGCACCGTTTGTTATTAAAGTAACAACTGCCACTGCTTTTTTGAACTTACGTGGTGACTTAATTGAAATGTAAAAGTCTAACACGGTTCCTACAg                                                                                           <  3:262073/141‑1 (MQ=255)

                                                   aaTTAGTATTCATG‑TTTTACTTTGCGAATGCTTGTAAATAATCTAGCACCGTTTGTTATTAAAGTAACAACTGCCACTGCTTTTTTGAACTTACGTGGTGACTTAATTGAAATGTAAAAGTCTAACACGGTTCCTACAGTg                                                                                         <  1:420158/141‑1 (MQ=255)

                                                    aTTAGTATTCATG‑TTTTACTTTGCGAATGCTTGTAAATAATCTAGCACCGTTTGTTATTAAAGTAACAACTGCCACTGCTTTTTTGAACTTACGTGGTGACTTAATTGAAATGTAAAAGTCTAACACGGTTCCTACAGTGc                                                                                        >  1:481438/1‑141 (MQ=255)

                                                    aTTAGTATTCATG‑TTTTACTTTGCGAATGCTTGTAAATAATCTAGCACCGTTTGTTATTAAAGTAACAACTGCCACTGCTTTTTTGAACTTACGTGGTGACTTAATTGAAATGTAAAAGTCTAACACGGTTCCTACAGTGc                                                                                        >  1:869906/1‑141 (MQ=255)

                                                    aTTAGTATTCATG‑TTTTACTTTGCGAATGCTTGTAAATAATCTAGCACCGTTTGTTATTAAAGTAACAACTGCCACTGCTTTTTTGAACTTACGTGGTGACTTAATTGAAATGTAAAAGTCTAACACGGTTCCTACAGTGc                                                                                        <  4:449789/141‑1 (MQ=255)

                                                     ttAGTATTCATG‑TTTTACTTTGCGAATGCTTGTAAATAATCTAGCACCGTTTGTTATTAAAGTAACAACTGCCACTGCTTTTTTTTTTTTTTTTGGTGACTTAATTGAAATGTAAAAGTCTAACACGGTTCCTACAGTGCt                                                                                       >  2:555298/1‑141 (MQ=255)

                                                     ttAGTATTCATG‑TTTTACTTTGCGAATGCTTGTAAATAATCTAGCACCGTTTGTTATTAAAGTAACAACTGCCACTGCTTTTTTGAACTTACGTGGTGACTTAATTGAAATGTAAAAGTCTAACACGGTTCCTACAGTGCt                                                                                       <  4:20316/141‑1 (MQ=255)

                                                      tAGTATTCATG‑TTTTACTTTGCGAATGCTTGTAAATAATCTAGCACCGTTTGTTATTAAAGTAACAACTGCCACTGCTTTTTTGAACTTACGTGGTGACTTAATTGAAATGTAAAAGTCTAACACGGTTCCTACAGTGCtt                                                                                      <  1:450202/141‑1 (MQ=255)

                                                          aTTCATG‑TTTTACTTTGCGAATGCTTGTAAATAATCTAGCACCGTTTGTTATTAAAGTAACAACTGCCACTGCTTTTTTGAACTTACGTGGTGACTTAATTGAAATGTAAAAGTCTAACACGGTTCCTACAGTGCTTAAAc                                                                                  <  3:160050/141‑1 (MQ=255)

                                                            tCATG‑TTTTACTTTGCGAATGCTTGTAAATAATCTAGCACCGTTTGTTATTAAAGTAACAACTGCCACTGCTTTTTTGAACTTACGTGGTGACTTAATTGAAATGTAAAAGTCTAACACGGTTCCTACAGTGCTTAAACCt                                                                                <  1:210139/141‑1 (MQ=255)

                                                              aTG‑TTTTACTTTGCGAATGCTTGTAAATAATCTAGCACCGTTTGTTATTAAAGTAACAACTGCCACTGCTTTTTTGAACTTACGTGGTGACTTAATTGAAATGTAGAAGTCTAACACGGTTCCTACAGTGCTTAAACCTaa                                                                              <  1:783604/141‑1 (MQ=255)

                                                              aTG‑TTTTACTTTGCGAATGCTTGTAAATAATCTAGCACCGTTTGTTATTAAAGTAACAACTGCCACTGCTTTTTTGAACTTACGTGGTGACTTAATTGAAATGTAAAAGTCTAACACGGTTCCTACAGTGCTTAAACCTaa                                                                              <  1:708169/141‑1 (MQ=255)

                                                                    ttACTTTGCGAATGCTTGTAAATAATCTAGCACCGTTTGTTATTAAAGTAACAACTGCCACTGCTTTTTTTTTTTTTTTTTTTGACTTAATTGAAATGTAAAAGTCTAACACGGTTCCTACAGTGCTTAAACCTAATGcaa                                                                         >  2:20487/1‑141 (MQ=255)

                                                                    ttACTTTGCGAATGCTTGTAAATAATCTAGCACCGTTTGTTATTAAAGTAACAACTGCCACTGCTTTTTTGAACTTACGTGGTGACTTAATTGAAATGTAAAAGTCTAACACGGTTCCTACAGTGCTTAAACCTAATGaaa                                                                         >  2:505663/1‑141 (MQ=255)

                                                                        tttGCGAATGCTTGTAAATAATCTAGCACCGTTTGTTATTAAAGTAACAACTGCCACTGCTTTTTTGAACTTACGTGGTGACTTAATTGAAATGTAAAAGTCTAACACGGTTCCTACAGTGCTTAAACCTAATGAAAGata                                                                     >  2:913233/1‑141 (MQ=255)

                                                                         ttGCGAATGCTTGTAAATAATCTAGCACCGTTTGTTATTAAAGTAACAACTGCCACTGCTTTTTTGAACTTACGTGGTGACTTAATTGAAATGTAAAAGTCTAACACGGTTCCTACAGTGCTTAAACCTAATGAAAGatat                                                                    >  2:466657/1‑141 (MQ=255)

                                                                           gCGAATGCTTGTAAATAATCTAGCACCGTTTGTTATTAAAGTAACAACTGCCACTGCTTTTTTGAACTTACGTGGTGACTTAATTGAAATGTAAAAGTCTAACACGGTTCCTACAGTGCTTAAACCTAATGAAAGATATAc                                                                  >  1:347155/1‑141 (MQ=255)

                                                                           gCGAATGCTTGTAAATAATCTAGCACCGTTTGTTATTAAAGTAACAACTGCCACTGCTTTTTTGAACTTACGTGGTGACTTAATTGAAATGTAAAAGTCTAACACGGTTCCTACAGTGCTTAAACCTAATGAAAGATATAc                                                                  <  2:126317/141‑1 (MQ=255)

                                                                              aaTGCTTGTAAATAATCTAGCACCGTTTGTTATTAAAGTAACAACTGCCACTGCTTTTTTGAACTTACGTGGTGACTTAATTGAAATGTAAAAGTCTAACACGGTTCCTACAGTGCTTAAACCTAATGAAAGATATACTaa                                                               <  1:174349/141‑1 (MQ=255)

                                                                                tGCTTGTAAATAATCTAGCACCGTTTGTTATTAAAGTAACAACTGCCACTGCTTTTTTGAACTTACGTGGTGACTTAATTGAAATGTAAAAGTCTAACACGGTTCCTACAGTGCTTAAACCTAATGAAAGATATACTAAtt                                                             >  1:163122/1‑141 (MQ=255)

                                                                                 gCTTGTAAATAATCTGGCACCGTTTGTTATTAAAGTAACAACTGCCACTGCTTTTTTGAACTTACGTGGTGACTTAATTGAAATGTAAAAGTCTAACACGGTTCCTACAGTGCTTAAACCTAATGAAAGATATACTAAttt                                                            <  1:215497/141‑1 (MQ=255)

                                                                                 acTTGTAAATAATCTAGCACCGTTTGTTATTAAAGTAACAACTGCCACTGCATCATTGACCTTACGTGGTGACTTAATTGAAATGTAAAAGTCTAACACGGTTCCTACAGTGCTCAAACCTAATGCCAGATATACTAAttt                                                            >  2:241408/2‑141 (MQ=255)

                                                                                  cTTGTAAATAATCTAGCACCGTTTGTTATTAAAGTAACAACTGCCACTGCTTTTTTGAACTTACGTGGTGACTTAATTGAAATGTAAAAGTCTAACACGGTTCCTACAGTGCTTAAACCTAATGAAAGATATACTAAtttt                                                           <  2:1143469/141‑1 (MQ=255)

                                                                                  cTTGTAAATAATCTAGCACCGTTTGTTATTAAAGTAACAACTGCCACTGCTTTTTTGAACTTACGTGGTGACTTAATTGAAATGTAAAAGTCTAACACGGTTCCTACAGTGCTTAAACCTAATGAAAGATATACTAAtttt                                                           <  3:377463/141‑1 (MQ=255)

                                                                                   ttGTAAATAATCTAGCACCGTTTGTTATTAAAGTAACAACTGCCACTGCTTTTTTGAACTTACGTGGTGACTTAATTGAAATGTAAAAGTCTAACACGGTTCCTACAGTGCTTAAACCTAATGAAAGATATACTAAttttt                                                          <  1:1134770/141‑1 (MQ=255)

                                                                                      tAAATAATCTAGCACCGTTTGTTATTAAAGTAACAACTGCCACTGCTTTTTTGAACTTACGTGGTGACTTAATTGAAATGTAAAAGTCTAACACGGTTCCTACAGTGCTTAAACCTAATGAAAGATATACTAATtttttat                                                       >  2:1159319/1‑141 (MQ=255)

                                                                                       aaaTAATCTAGCACCGTTTGTTATTAAAGTAACAACTGCCACTGCTTTTTTGAACTTACGTGGTGACTTAATTGAAATGTAAAAGTCTAACACGGTTCCTACAGTGCTTAAACCTAATGAAAGATATACTAATtttttatt                                                      >  2:286464/1‑141 (MQ=255)

                                                                                        aataatCTAGCACCGTTTGTTATTAAAGTAACAACTGCCACTGCTTTTTTGAACTTACGTGGTGACTTAATTGAAATGTAAAAGTCTAACACGGTTCCTACAGTGCTTAAACCTAATGAAAGATATACTAATtttttattt                                                     <  1:819403/141‑1 (MQ=255)

                                                                                          taatCTAGCACCGTTTGTTATTAAAGTAACAACTGCCACTGCTTTTTTGAACTTACGTGGTGACTTAATTGAAATGTAAAAGTCTAACACGGTTCCTACAGTGCTTAAACCTAATGAAAGATATACTAATtttttattttt                                                   >  1:243304/1‑141 (MQ=255)

                                                                                           aatCTAGCACCGTTTGTTATTAAAGTAACAACTGCCACTGCTTTTTTGAACTTACGTGGTGACTTAATTGAAATGTAAAAGTCTAACACGGTTCCTACAGTGCTTAAACCTAATGAAAGATATACTAATtttttattttta                                                  <  1:533081/141‑1 (MQ=255)

                                                                                            atCTAGCGCCGTTTGTTATTAAAGTAACAACTGCCACTGCTTTTTTGAACTTACGTGGTGACTTAATTGAAATGTAAAAGTCTAACACGGTTCCTACAGTGCTTAAACCTAATGAAAGATATACTAATTTTTTATTTTTAg                                                 >  3:444325/1‑141 (MQ=255)

                                                                                            atCTAGCACCGTTTGTTATTAAAGTAACAACTGCCACTGCTTTTTTGAACTTACGTGGTGACTTAATTGAAATGTAAAAGTCTAACACGGTTCCTACAGTGCTTAAACCTAATGAAAGATATACTAATTTTTTATTTTTAg                                                 >  3:671/1‑141 (MQ=255)

                                                                                            atCTAGCACCGTTTGTTATTAAAGTAACAACTGCCACTGCTTTTTTGAACTTACGTGGTGACTTAATTGAAATGTAAAAGTCTAACACGGTTCCTACAGTGCTTAAACCTAATGAAAGATATACTAATTTTTTATTTTTAg                                                 >  4:64279/1‑141 (MQ=255)

                                                                                            atCTAGCACCGTTTGTTATTAAAGTAACAACTGCCACTGCTCTTTTGAACTTACGTGGTGACTTAATTGAAATGTAAACGTGTCACACGGTTTCTACAGTGATTACACCTAATGAAAGATATACTACTTTTTTATTTTTAg                                                 >  2:572551/1‑141 (MQ=255)

                                                                                             tCTAGCACCGTTTGTTATTAAAGTAACAACTGCCACTGCTTTTTTGAACTTACGTGGTGACTTAATTGAAATGTAAAAGTCTAACACGGTTCCTACAGTGCTTAAACCTAATGAAAGATATACTAATTTTTTATTTTTAGc                                                >  2:986820/1‑141 (MQ=255)

                                                                                                       tttGTTATTAAAGTAACAACTGCCACTGCTTTTTTGAACTTACGTGGTGACTTAATTGAAATGTAAAAGTCTAACACGGTTCCTACAGTGCTTAAACCTAATGAAAGATATACTAATTTTTTATTTTTAGCATGATATTta                                      <  1:857326/141‑1 (MQ=255)

                                                                                                        ttGTTATTAAAGTAACAACTGCCACTGCTTTTTTGAACTTACGTGGTGACTTAATTGAAATGTAAAAGTCTAACACGGTTCCTACAGTGCTTAAACCTAATGAAAGATATACTAATTTTTTATTTTTAGCATGATATTtat                                     >  1:729961/1‑141 (MQ=255)

                                                                                                        ttGTTATTAAAGTAACAACTGCCACTGCTTTTTTGAACTTACGTGGTGACTTAATTGAAATGTAAAAGTCTAACACGGTTCCTACAGTGCTTAAACCTAATGAAAGATATACTAATTTTTTATTTTTAGCATGATATTtat                                     <  2:805217/141‑1 (MQ=255)

                                                                                                        ttGTTATTAAAGTAACAACTGCCACTGCTTTTTTGAACTTACGTGGTGACTTAATTGAAATGTAAAAGTCTAACACGGTTCCTACAGTGCTTAAACCTAATGAAAGATATACTAATTTTTTATTTTTAGCATGATATTtat                                     >  2:885847/1‑141 (MQ=255)

                                                                                                         tGTTATTAAAGTAACAACTGCCACTGCTTTTTTGAACTTACGTGGTGACTTAATTGAAATGTAAAAGTCTAACACGGTTCCTACAGTGCTTAAACCTAATGAAAGATATACTAATTTTTTATTTTTAGCATGATATTtata                                    <  1:257369/141‑1 (MQ=255)

                                                                                                           ttattaAAGTAACAACTGCCACTGCTTTTTTGAACTTACGTGGTGACTTAATTGAAATGTAAAAGTCTAACACGGTTCCTACAGTGCTTAAACCTAATGAAAGATATACTAATTTTTTATTTTTAGCATGATATTTATAGc                                  >  1:4194/1‑141 (MQ=255)

                                                                                                                aaaGTAACAACTGCCACTGCTTTTTTTTTTTTTTTTGGTGACTTAATTGAAATGTAAAAGTCTAACACGGTTCCTACAGTGCTTAAACCTAATGAAAGATATACTAATTTTTTATTTTTAGCATGATATTTATAGCCATTg                             >  2:74569/1‑141 (MQ=255)

                                                                                                                 aaGTAACAACTGCCACTGCTTTTTTGAACTTACGTGGTGACTTAATTGAAATGTAAAAGTCTAACACGGTTCCTACAGTGCTTAAACCTAATGAAAGATATACTAATTTTTTATTTTTAGCATGATATTTATAGCCATTGt                            <  3:257396/141‑1 (MQ=255)

                                                                                                                       caacTGCCACTGCTTTTTTTTTTTTTTTTTGTGACTTAATTGAAATGTAAAAGTCTAACACGGTTCCTACAGTGCTTAAACCTAATGAAAGATATACTAATTTTTTATTTTTAGCATGATATTTATAGCCATTGTAGCCGt                      >  2:491146/1‑141 (MQ=255)

                                                                                                                       caacTGCCACTGCTTTTTTTAACTTACGTGGTGACTTAATTGAAATGTAAAAGTCTAACACGGTTCCTACAGTGCTTAAACCTAATGAAAGATATACTAATTTTTTATTTTTAGCATGATATTTATAGCCATTGTAGCCGt                      >  2:921608/1‑141 (MQ=255)

                                                                                                                       caacTGCCACTGCTTTTTTGAACTTACGTGGTGACTTAATTGAAATGTAAAAGTCTAACACGGTTCCTACAGTGCTTAAACCTAATGAAAGATATACTAATTTTTTATTTTTAGCATGATATTTATAGCCATTGTAGCCGt                      >  3:428607/1‑141 (MQ=255)

                                                                                                                       caacTGCCACTGCTTTTTTGAACTTACGTGGTGACTTAATTGAAATGTAAAAGTCTAACACGGTTCCTACAGTGCTTAAACCTAATGAAAGATATACTAATTTTTTATTTTTAGCATGATATTTATAGCCATTGTAGCCGt                      >  1:958491/1‑141 (MQ=255)

                                                                                                                       caacTGCCACTGCTTTTTTGAACTTACGTGGTGACTTAATTGAAATGTAAAAGTCTAACACGGTTCCTACAGTGCTTAAACCTAATGAAAGATATACTAATTTTTTATTTTTAGCATGATATTTATAGCCATTGTAGCCGt                      <  4:143111/141‑1 (MQ=255)

                                                                                                                         acTGCCACTGCTTTTTTTTTTTTTTTTTGTGACTTAATTGAAATGTAAAAGTCTAACACGGTTCCTACAGTGCTTAAACCTAATGAACGATATACTAATTTTTTATTTTTAGCATGATTTTTTTAGTTTTTGTAGCCGTct                    >  2:44210/1‑140 (MQ=255)

                                                                                                                          cTGCCACTGCTTTTTTGTTCTTTCGTGGTGACTTAATTGAAATGTAAAAGTCTAACACGGTTCCTACAGTGCTTAAACCTAATGAAAGATATACTAATTTTTTATTTTTAGCATGATATTTATAGCCATTGTAGCCGTCGa                   >  2:796982/1‑141 (MQ=255)

                                                                                                                          cTGCCACTGCTTTTTTGAACTTACGTGGTGACTTAATTGAAATGTAAAAGTCTAACACGGTTCCTACAGTGCTTAAACCTAATGAAAGATATACTAATTTTTTATTTTTAGCATGATATTTATAGCCATTGTAGCCGTCGa                   <  3:292381/141‑1 (MQ=255)

                                                                                                                            gCCACTGCTTTTTTGAACTTACGTGGTGACTTAATTGAAATGTAAAAGTCTAACACGGTTCCTACAGTGCTTAAACCTAATGAAAGATATACTAATTTTTTATTTTTAGCATGATATTTATAGCCATTGTAGCCGTCGACt                 <  1:638922/141‑1 (MQ=255)

                                                                                                                            gCCACTGCTTTTTTGAACTTACGTGGTGACTTAATTGAAATGTAAAAGTCTAACACGGTTCCTACAGTGCTTAAACCTAATGAAAGATATACTAATTTTTTATTTTTAGCATGATATTTATAGCCATTGTAGCCGTCGACt                 >  2:932169/1‑141 (MQ=255)

                                                                                                                             ccACTGCTTTTTTTTTTTTTTTTGGTGACTAAATTGAAATGTAAAAGTTTAACACGGTTCCTACAGTGCTTAAACCTAATGAAAGATATACTAATTTTTTATTTTTAGCATGTTATTTTTAGCCATTGTAGCCGTCGACTa                >  2:1044285/1‑141 (MQ=255)

                                                                                                                               aCTGCTTTTTTGAACTTACGTGGTGACTTAATTGAAATGTAAAAGTCTAACACGGTTCCTACAGTGCTTAAACCTAATGAAAGATATACTAATTTTTTTTTTTTTTTTTTATATTTATAGCCATTGTAGCCGTCTACTATg              >  4:366273/1‑141 (MQ=255)

                                                                                                                                 tGCTTTTTTGAACTTACGTGGTGACTTAATTGAAATGTAAAAGTCTAACACGGTTCCTACAGTGCTTAAACCTAATGAAAGATATACTAATTTTTTATTTTTAGCATGATATTTATAGCCATTGTAGCCGTCGACTATGaa            >  1:128770/1‑141 (MQ=255)

                                                                                                                                    ttttttGAACTTACGTGGTGACTTAATTGAAATGTAAAAGTCTAACACGGTTCCTACAGTGCTTAAACCTAATGAAAGATATACTAATTTTTTATTTTTAGCATGATATTTATAGCCATTGTAGCCGTCGACTATGAAAcc         <  2:381799/141‑1 (MQ=255)

                                                                                                                                      ttttGAACTTACGTGTTGACTTAATTGAAATGTAAAAGTCTAACACGGTTCCTACAGTGCTTAAACCTAATGAAAGATATACTAATTTTTTATTTTTAGCATGATATTTATAGCCATTGTAGCCGTCGACTATGAAACCTg       <  4:317150/141‑1 (MQ=255)

                                                                                                                                       tttGAACTTACGTGGTGACTTAATTGAAATGTAAAAGTCTAACACGGTTCCTACAGTGCTTAAACCTAATGAAAGATATACTAATTTTTTATTTTTAGCATGATATTTATAGCCATTGTAGCCGTCGACTATGAAACCTGc      >  1:399659/1‑141 (MQ=255)

                                                                                                                                       tttGAACTTACGTGGTGACTTAATTGAAATGTAAAAGTCTAACACGGTTCCTACAGTGCTTAAACCTAATGAAAGATATACTAATTTTTTATTTTTAGCATGATATTTATAGCCATTGTAGCCGTCGACTATGAAACCTGc      >  1:399270/1‑141 (MQ=255)

                                                                                                                                          gAACTTACGTGGTGACTTAATTGAAATGTAAAAGCCTAACACGGTTCCTACAGTGCTTCAACCTAATGAAAGATATACTAATTTTTTATTTTTAGCATGATATTTATAGCCATTGTAGCCGTCGACTATGAAACCTGCAAc   >  1:547095/1‑141 (MQ=255)

AAGGATTGGTCGAATGAACTTCATGTTCTATGTTCGACCAATCCTTTTCTGAATTAGTATTCATG‑TTTTACTTTGCGAATGCTTGTAAATAATCTAGCACCGTTTGTTATTAAAGTAACAACTGCCACTGCTTTTTTGAACTTACGTGGTGACTTAATTGAAATGTAAAAGTCTAACACGGTTCCTACAGTGCTTAAACCTAATGAAAGATATACTAATTTTTTATTTTTAGCATGATATTTATAGCCATTGTAGCCGTCGACTATGAAACCTGCGACA  >  USA300TCH1516_ALE/492710‑492988

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