BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	USA300TCH1516_ALE	61,741	C→T	6.4%	E265K (GAA→AAA)	USA300TCH1516_00053 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	USA300TCH1516_ALE	61,741	0	C	T	6.4%	425.7 / 3.9	139	E265K (GAA→AAA)	USA300TCH1516_00053	hypothetical protein
Reads supporting (aligned to +/- strand): ref base C (54/73); new base T (10/2); total (64/75)
Fisher's exact test for biased strand distribution p-value = 1.23e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 141 total aligned reads displayed.

ATATTATTTTCACAAAAAAAGAATGCTTTTTGTTAGAGGCACTGAACCCTTAAAATGGGAACTTAATAAAAACACTATGTTCTAGGCTGCTAATTCTCTGTATTTTACAGGGGATAAGTAGCCTAGTTTTTGTTGAATTCGATTATTATTATAGTTTTTAATGTA‑TTTTCGACAATATCTATTACAATATGATTAGAGTTATTAAGCTGATTATTGATGTAAAAAGTTTCAGACTTTAGCGAGGAATGGAAACTTTCTATCGGGGCGTTATCGGCAGGTG  >  USA300TCH1516_ALE/61602‑61881

 tattatTTTCACAAAAAAAGAATGCTTTTTGTTAGAGGCACTGAACCCTTAAAATGGGAACTTAATAAAAACACTATGTTCTAGGCTGCTAATTCTCTGTATTTTACAGGGGATAAGTAGCCTAGTTTTTGTTGAATTCGa                                                                                                                                             <  2:961868/141‑1 (MQ=255)

   ttattTTCACAAAAAAAGAATGCTTTTTGTTAGAGGCACTGAACCCTTAAAATGGGAACTTAATAAAAACACTATGTTCTAGGCTGCTAATTCTCTGTATTTTACAGGGGATAAGTAGCCTAGTTTTTGTTGAATTCGatt                                                                                                                                           <  1:720133/141‑1 (MQ=255)

   ttattTTCACAAAAAAAGAATGCTTTTTGTTAGAGGCACTGAACCCTTAAAATGGGAACTTAATAAAAACACTATGTTCTAGGCTGCTAATTCTCTGTATTTTACAGGGGATAAGTAGCCTAGTTTTTGTTGAATTCGatt                                                                                                                                           >  1:489073/1‑141 (MQ=255)

     attTTCACAAAAAAAGAATGCTTTTTGTTAGAGGCACTGAACCCTTAAAATGGGAACTTAATAAAAACACTATGTTCTAGGCTGCTAATTCTCTGTATTTTACAGGGGATAAGTAGCCTAGTTTTTGTTGAATTCGattat                                                                                                                                         <  1:629804/141‑1 (MQ=255)

         tCACAAAAAAAGAATGCTTTTTGTTAGAGGCACTGAACCCTTAAAATGGGAACTTAATAAAAACACTATGTTCTAGGCTGCTAATTCTCTGTATTTTACAGGGGATAAGTAGCCTAGTTTTTGTTGAATTCGattattatt                                                                                                                                     >  3:405145/1‑141 (MQ=255)

          cacaAAAAAAGAATGCTTTTTGTTAGAGGCACTGAACCCTTAAAATGGGAACTTAATAAAAACACTATGTTCTAGGCTGCTAATTCTCTGTATTTTACAGGGGATAAGTAGCCTAGTTTTTGTTGAATTCGAttattatta                                                                                                                                    <  2:680350/141‑1 (MQ=255)

           acaAAAAAAGAATGCTTTTTGTTAGAGGCACTGAACCCTTAAAATGGGAACTTAATAAAAACACTATGTTCTAGGCTGCTAATTCTCTGTATTTTACAGGGGATAAGTAGCCTAGTTTTTGTTGAATTCGATtattattat                                                                                                                                   <  3:31609/141‑1 (MQ=255)

             aaaaaaaGAATGCTTTTTGTTAGAGGCACTGAACCCTTAAAATGGGAACTTAATAAAAACACTATGTTCTAGGCTGCTAATTCTCTGTATTTTACAGGGGATAAGTAGCCTAGTTTTTGTTGAATTCGATTATTATTATAg                                                                                                                                 >  1:104648/1‑141 (MQ=255)

              aaaaaaGAATGCTTTTTGTTAGAGGCACTGAACCCTTAAAATGGGAACTTAATAAAAACACTATGTTCTAGGCTGCTAATTCTCTGTATTTTACAGGGGATAAGTAGCCTAGTTTTTGTTGAATTCGATTATTATTATAGt                                                                                                                                >  1:893194/1‑141 (MQ=255)

                aaaaGAATGCTTTTTGTTAGAGGCACTGAACCCTTAAAATGGGAACTTAATAAAAACACTATGTTCTAGGCTGCTAATTCTCTGTATTTTACAGGGGATAAGTAGCCTAGTTTTTGTTGAATTCGATTATTATTATAGttt                                                                                                                              <  3:271710/141‑1 (MQ=255)

                 aaaGAATGCTTTTTGTTAGAGGCACTGAACCCTTAAAATGGGAACTTAATAAAAACACTATGTTCTAGGCTGCTAATTCTCTGTATTTTACAGGGGATAAGTAGCCTAGTTTTTGTTGAATTCGATTATTATTATAGtttt                                                                                                                             >  2:633380/1‑141 (MQ=255)

                  aaGAATGCTTTTTGTTAGAGGCACTGAACCCTTAAAATGGGAACTTAATAAAAACACTATGTTCTAGGCTGCTAATTCTCTGTATTTTACAGGGGATAAGTAGCCTAGTTTTTGTTGAATTCGATTATTATTATAGttttt                                                                                                                            >  2:792709/1‑141 (MQ=255)

                    gAATGCTTTTTGTTAGAGGCACTGAACCCTTAAAATGGGAACTTAATAAAAACACTATGTTCTAGGCTGCTAATTCTCTGTATTTTACAGGGGATAAGTAGCCTAGTTTTTGTTGAATTCGATTATTATTATAGTTTTTaa                                                                                                                          <  3:247993/141‑1 (MQ=255)

                    gAATGCTTTTTGTTAGAGGCACTGAACCCTTAAAATGGGAACTTAATAAAAACACTATGTTCTAGGCTGCTAATTCTCTGTATTTTACAGGGGATAAGTAGCCTAGTTTTTGTTGAATTCGATTATTATTATAGTTTTTaa                                                                                                                          <  4:259137/141‑1 (MQ=255)

                      aTGCTTTTTGTTAGAGGCACTGAACCCTTAAAATGGGAACTTAATAAAAACACTATGTTCTAGGCTGCTAATTCTCTGTATTTTACAGGGGATAAGTAGCCTAGTTTTTGTTGAATTCGATTATTATTATAGTTTTTAAtt                                                                                                                        >  4:324698/1‑140 (MQ=255)

                       tGCTTTTTGTTAGAGGCACTGAACCCTTAAAATGGGAACTTAATAAAAACACTATGTTCTAGGCTGCTAATTCTCTGTATTTTACAGGGGATAAGTAGCCTAGTTTTTGTTGAATTCGATTATTATTATAGTTTTTAATGt                                                                                                                       <  1:1016023/141‑1 (MQ=255)

                        gCTTTTTGTTAGAGGCACTGAACCCTTAAAATGGGAACTTAATAAAAACACTATGTTCTAGGCTGCTAATTCTCTGTATTTTACAGGGGATAAGTAGCCTAGTTTTTGTTGAATTCGATTATTATTATAGTTTTTAATGTa                                                                                                                      >  3:93346/1‑141 (MQ=255)

                          tttttGTTAGAGGCACTGAACCCTTAAAATGGGAACTTAATAAAAACACTATGTTCTAGGCTGCTAATTCTCTGTATTTTACAGGGGATAAGTAGCCTAGTTTTTGTTGAATTCGATTATTATTATAGTTTTTAATGTA‑tt                                                                                                                   >  1:207397/1‑141 (MQ=255)

                            tttGTTAGAGGCACTGAACCCTTAAAATGGGAACTTAATAAAAACACTATGTTCTAGGCTGCTAATTCTCTGTATTTTACAGGGGATAAGTAGCCTAGTTTTTGTTGAATTCGATTATTATTATAGTTTTTAATGTA‑tttt                                                                                                                 <  4:140693/141‑1 (MQ=255)

                              tGTTAGAGGCACTGAACCCTTAAAATGGGAACTTAATAAAAACACTATGTTCTAGGCTGCTAATTCTCTGTATTTTACAGGGGATAAGTAGCCTAGTTTTTTTTTTTTTTTTTTATTATTATAGTTTTTAATGTA‑TTTTCg                                                                                                               >  2:216982/1‑141 (MQ=18)

                               gTTAGAGGCACTGAACCCTTAAAATGGGAACTTAATAAAAACACTATGTTCTAGGCTGCTAATTCTCTGTATTTTACAGGGGATAAGTAGCCTAGTTTTTGTTGAATTCGATTATTATTATAGTTTTTAATGTG‑TTTTCGa                                                                                                              <  1:858282/141‑1 (MQ=255)

                               gTTAGAGGCACTGAACCCTTAAAATGGGAACTTAATAAAAACACTATGTTCTAGGCTGCTAATTCTCTGTATTTTACAGGGGATAAGTAGCCTAGTTTTTGTTGAATTCGATTATTATTATAGTTTTTAATGTA‑TTTTCGa                                                                                                              >  1:29582/1‑141 (MQ=255)

                                ttAGAGGCACTGAACCCTTAAAATGGGAACTTAATAAAAACACTATGTTCTAGGCTGCTAATTCTCTGTATTTTACAGGGGATAAGTAGCCTAGTTTTTGTTGAATTCGATTATTATTATAGTTTTTAATGTA‑TTTTTGAc                                                                                                             >  2:763493/1‑141 (MQ=255)

                                ttAGAGGCACTGAACCATTAAAATGGGAACTTAATAAAAACACTATGTTCTAGGCTGCTAATTCTCTGTATTTTACAGGGGATAAGTAGCCTAGTTTTTGTTGAATTCGATTATTATTATAGTTTTTAATGTA‑TTTTCGAc                                                                                                             <  2:1076093/141‑1 (MQ=255)

                                  agagGCACTGAACCCTTAAAATGGGAACTTAATAAAAACACTATGTTCTAGGCTGCTAATTCTCTGTATTTTACAGGGGATAAGTAGCCTAGTTTTTGTTGAATTCGATTATTATTATAGTTTTTAATGTA‑TTTTCGACaa                                                                                                           <  1:543119/141‑1 (MQ=255)

                                   ttttctcTGAACCCTTAAAATGGGAACTTATTAGAAACGCGACGTTCTGGGCTGCAGATTCCCTGTCTTTTACAGTTGTTGTGTTGGCTAGTTTATGTTGAATTCGAATATTATTATAGTTTTTAATGTA‑TTTTCGACAat                                                                                                          >  2:641359/7‑141 (MQ=255)

                                   gagGCACTGAACCCTTAAAATGGGAACTTAATAAAAACACTATGTTCTAGGCTGCTAATTCTCTGTATTTTACAGGGGATAAGTAGCCTAGTTTTTGTTGAATTCGATTATTATTATAGTTTTTAATGTA‑TTTTCGACAat                                                                                                          >  2:940766/1‑141 (MQ=255)

                                   gagGCACTGAACCCTTAAAATGGGAACTTAATAAAAACACTATGTTCTAGGCTGCTAATTCTCTGTATTTTACAGGGGATAAGTAGCCTAGTTTTTGTTGAATTCGATTATTATTATAGTTTTTAATGTA‑TTTTCGACAat                                                                                                          >  1:1057717/1‑141 (MQ=255)

                                   gagGCACTGAACCCTTAAAATGGGAACTTAATAAAAACACTATGTTCTAGGCTGCTAATTCTCTGTATTTTACAGGGGATAAGTAGCCTAGTTTTGGTTGAATTCGGTTATTATTATAGTTTTTAATGTA‑TTTTCGACAat                                                                                                          >  1:1068156/1‑141 (MQ=255)

                                      gCACTGAACCCTTAAAATGGGAACTTAATAAAAACACTATGTTCTAGGCTGCTAATTCTCTGTATTTTACAGGGGATAAGTAGCCTAGTTTTTGTTGAATTCGATTATTATTATAGTTTTTAATGTA‑TTTTCGACAATATc                                                                                                       <  2:1068456/141‑1 (MQ=255)

                                       cACTGAACCCTTAAAATGGGAACTTAATAAAAACACTATGTTCTAGGCTGCTAATTCTCTGTATTTTACAGGGGATAAGTAGCCTAGTTTTTGTTGAATTCGATTATTATTATAGTTTTTAATGTA‑TTTTCGACAATATCt                                                                                                      <  2:545863/141‑1 (MQ=255)

                                          tGAACCCTTAAAATGGGAACTTAATAAAAACACTATGTTCTAGGCTGCTAATTCTCTGTATTTTACAGGGGATAAGTAGCCTAGTTTTTGTTGAATTCGATTATTATTATAGTTTTTAATGTA‑TTTTCGACAATATCTAtt                                                                                                   <  1:745385/141‑1 (MQ=255)

                                          tGAACCCTTAAAATGGGAACTTAATAAAAACACTATGTTCTAGGCTGCTAATTCTCTGTATTTTACAGGGGATAAGTAGCCTAGTTTTTGTTGAATTCGATTATTATTATAGTTTTTAATGTA‑TTTTCGACAATATCTAtt                                                                                                   <  2:68137/141‑1 (MQ=255)

                                          tGAACCCTTAAAATGGGAACTTAATAAAAACACTATGTTCTAGGCTGCTAATTCTCTGTATTTTACAGGGGATAAGTAGCCTAGTTTTTGTTGAATTCGATTATTATTATAGTTTTTAATGTA‑TTTTCGACAATATCTAtt                                                                                                   <  3:177480/141‑1 (MQ=255)

                                               ccTTAAAATGGGAACTTAATAAAAACACTATGTTCTAGGCTGCTAATTCTCTGTATTTTACAGGGGATAAGTAGCCTAGTTTTTGTTGAATTCGATTATTATTATAGTTTTTAATGTA‑TTTTCGACAATATCTATTACAat                                                                                              >  4:67813/1‑141 (MQ=255)

                                                  tAAAATGGGAACTTAATAAAAACACTATGTTCTAGGCTGCTAATTCTCTGTATTTTACAGGGGATAAGTAGCCTAGTTTTTGTTGAATTCGATTATTATTATAGTTTTTAATGTA‑TTTTCGACAATATCTATTACAATATg                                                                                           <  4:109209/141‑1 (MQ=255)

                                                   aaaaTGGGAACTTAATAAAAACACTATGTTCTAGGCTGCTAATTCTCTGTATTTTACAGGGGATAAGTAGCCTAGTTTTTGTTGAATTCGATTAGTATTATAGTTTTTAATGTA‑TTTTCGACAATATCTATTACAATATGa                                                                                          <  3:378501/141‑1 (MQ=255)

                                                    aaaTGGGAACTTAATAAAAACACTATGTTCTAGGCTGCTAATTCTCTGTATTTTACAGGGGATAAGTAGCCTAGTTTTTGTTGAATTCGATTATTATTATAGTTTTTAATGTA‑TTTTCGACAATATCTATTACAATATGAt                                                                                         <  1:173890/141‑1 (MQ=255)

                                                     aaTGGGAACTTAATAAAAACACTATGTTCTAGGCTGCTAATTCTCTGTATTTTACAGGGGATAAGTAGCCTAGTTTTTGTTGAATTCGATTATTATTATAGTTTTTAATGTA‑TTTTCGACAATATCTATTACAATATGAtt                                                                                        >  4:190251/1‑141 (MQ=255)

                                                     aaTGGGAACTTAATAAAAACACTATGTTCTAGGCTGCTAATTCTCTGTATTTTACAGGGGATAAGTAGCCTAGTTTTTGTTGAATTCGATTATTATTATAGTTTTTAATGTA‑TTTTCGACAATATCTATTACAATATGAtt                                                                                        <  2:415788/141‑1 (MQ=255)

                                                     aaTGGGAACTTAATAAAAACACTATGTTCTAGGCTGCTAATTCTCTGTATTTTACAGGGGATAAGTAGCCTAGTTTTTGTTGAATTCGATTATTATTATAGTTTTTAATGTA‑TTTTCGACAATATCTATTACAATATGAtt                                                                                        >  1:225249/1‑141 (MQ=255)

                                                       tGGGAACTTAATAAAAACACTATGTTCTAGGCTGCTAATTCTCTGTATTTTACAGGGGATAAGTAGCCTAGTTTTTGTTGAATTCGATTATTATTATAGTTTTTAATGTA‑TTTTCGACAATATCTATTACAATATGATTag                                                                                      <  4:26484/141‑1 (MQ=255)

                                                        gggAACTTAATAAAAACACTATGTTCTAGGCTGCTAATTCTCTGTATTTTACAGGGGATAAGTAGCCTAGTTTTTGTTGAATTCGATTATTATTATAGTTTTTAATGTA‑TTTTCGACAATATCTATTACAATATGATTaga                                                                                     <  4:249281/141‑1 (MQ=255)

                                                         ggAACTTAATAAAAACACTATGTTCTAGGCTGCTAATTCTCTGTATTTTACAGGGGATAAGTAGCCTAGTTTTTGTTGAATTCGATTATTATTATAGTTTTTAATGTA‑TTTTCGACAATATCTATTACAATATGATTagag                                                                                    <  2:423837/141‑1 (MQ=255)

                                                          gAACTTAATAAAAACACTATGTTCTAGGCTGCTAATTCTCTGTATTTTACAGGGGATAAGTAGCCTAGTTTTTGTTGAATTCGATTATTATTATAGTTTTTAATGTA‑TTTTCGACAATATCTATTACAATATGATTAGAGt                                                                                   <  2:143935/141‑1 (MQ=255)

                                                            aCTTAATAAAAACACTATGTTCTAGGCTGCTAATTCTCTGTATTTTACAGGGGATAAGTAGCCTAGTTTTTGTTGAATTCGATTATTATTATAGTTTTTAATGTA‑TTTTCGACAATATCTATTACAATATGATTAGAGtta                                                                                 >  1:396346/1‑141 (MQ=255)

                                                             cTTAATAAAAACACTATGTTCTAGGCTGCTAATTCTCTGTATTTTACAGGGGATAAGTAGCCTAGTTTTTGTTGAATTCGATTATTATTATAGTTTTTAATGTA‑TTTTCGACAATATCTATTACAATATGATTAGAGttat                                                                                >  1:202022/1‑141 (MQ=255)

                                                             cTTAATAAAAACACCATGTTCTAGGCTATTAATTCTCTGTATTTTACAGGGGATAAGTAGCCTAGTTTTTGTTGAATTTGATTATTATTATAGTTTTTAATGTACTTTTCGACAATATCTATTACAATATGATTagagcta                                                                                 <  2:1078005/141‑4 (MQ=14)

                                                              ttAATAAAAACACTATGTTCTAGGCTGCTAATTCTCTGTATTTTACAGGGGATAAGTAGCCTAGTTTTTGTTGAATTCGATTATTATTATAGTTTTTAATGTA‑TTTTCGACAATATCTATTACAATATGATTAGAGttatt                                                                               <  2:681598/141‑1 (MQ=255)

                                                                 ataaAAACACTATGTTCTAGGCTGCTAATTCTCTGTATTTTACAGGGGATAAGTAGCCTAGTTTTTGTTGAATTCGATTATTATTATAGTTTTTAATGTA‑TTTTCGACAATATCTATTACAATATGATTAGAGTTATTAAg                                                                            >  2:85395/1‑141 (MQ=255)

                                                                  taaAAACACTATGTTCTAGGCTGCTAATTCTCTGTATTTTACAGGGGATAAGTAGCCTAGTTTTTGTTGAATTCGATTATTATTATAGTTTTTAATGTA‑TTTTCGACAATATCTATTACAATATGATTAGAGTTATTAAGc                                                                           <  2:742483/141‑1 (MQ=255)

                                                                    aaaaCACTATGTTCTAGGCTGCTAATTCACTGTATTTTACAGGGGATAAGTAGCCTAGTTTTTGTTGAATTCGATTATTATTATAGTTTTTAATGTA‑TTTTCGACAATATCTATTACAATATGATTAGAGTTATTAAGCTg                                                                         <  4:421041/141‑1 (MQ=255)

                                                                    aaaaCACTATGTACTAGGCAGCTAATTCACTGTATGTTACGGGGGATAAGTAGCCTAGGTTTTGTTGAATTCGATTATTATTATAGTTTTTAATGTA‑TTTTCGACAATATCTATTACAATATGATTAGAGTTATTAAGCTg                                                                         <  4:420836/141‑1 (MQ=21)

                                                                        cacTATGTTCTAGGCTGCTAATTCTCTGTATTTTACAGGGGATAAGTAGCCTAGTTTTTGTTGAATTCGATTATTATTATAGTTTTTAATGTA‑TTTTCGACAATATCTATTACAATATGATTAGAGTTATTAAGCTGatta                                                                     >  1:120029/1‑141 (MQ=255)

                                                                           tATGTTCTAGGCTGCTAATTCTCTGTATTTTACAGGGGATAAGTAGCCTAGTTTTTGTTGAATTCGATTATTATTATAGTTTTTAATGTA‑TTTTCGACAATATCTATTACAATATGATTAGAGTTATTAAGCTGATTATTg                                                                  >  2:71046/1‑141 (MQ=255)

                                                                             tGTTCTAGGCTGCTAATTCTCTGTATTTTACAGGGGATAAGTAGCCTAGTTTTTGTTGAATTCGATTATTATTATAGTTTTTAATGTA‑TTTTCGACAATATCTATTACAATATGATTAGAGTTATTAAGCTGATTATTGAt                                                                >  1:578492/1‑141 (MQ=255)

                                                                             ggTTGTAGGCGGCGATTCCGCTGTATTTTACAGGGGATAAGTAGCCTAGTTTTAGATGAATTTGAGTATTATTATAGTTTTTAATGTA‑TTTTCGACAATATCTATTACAATATTGTTACAGTTTTTAAGCTGCTTATTGAt                                                                >  4:320894/2‑141 (MQ=21)

                                                                               ttCTAGGCTGCTAATTCTCTGTATTTTACAGGGGATAAGTAGCCTAGTTTTTGTTGAATTCGATTATTATTATAGTTTTTAATGTA‑TTTTCGACAATATCTATTACAATATGATTAGAGTTATTAAGCTGATTATTGATGt                                                              <  3:428950/141‑1 (MQ=255)

                                                                               ttCTAGGCTGCTAATTCTCTGTATTTTACAGGGGATAAGTAGCCTAGTTTTTGTTGAATTCGATTATTATTATAGTTTTTAATGTA‑TTTTCGACAATATCTATTACAATATGATTAGAGTTATTAAGCTGATTATTGATGt                                                              <  2:508700/141‑1 (MQ=255)

                                                                               ttCTAGGCTGCTAATTCTCTGTATTTTACAGGGGATAAGTAGCCTAGTTTTTGTTGAATTCGATTATTATTATAGTTTTTAATGTA‑TTTTCGACAATATCTATTACAATATGATTAGAGTTATTAAGCTGATTATTGATGt                                                              <  2:470865/141‑1 (MQ=255)

                                                                                 cTAGGCTGCTAATTCTCTGTATTTTACAGGGGATAAGTAGCCTATGTTTTGTTGAATTCGATTATTATTATAGTTTTTAATGTA‑TTTTCGACAATATCTATTACAATATGATTAGAGTTATTAAGCTGATTATTGATGTaa                                                            <  1:10387/141‑1 (MQ=255)

                                                                                 cTAGGCTGCTAATTCTCTGTATTTTACAGGGGATAAGTAGCCTAGTTTTTGTTGAATTCGATTATTATTATAGTTTTTAATGTA‑TTTTCGACAATATCTATTACAATATGATTAGAGTTATTAAGCTGATTATTGATGTaa                                                            <  2:525580/141‑1 (MQ=255)

                                                                                 cTAGGCTGCTAATTCTCTGTATTTTACAGGGGATAAGTAGCCTAGTTTTTGTTGAATTCGATTATTATTATAGTTTTTAATGTA‑TTTTCGACAATATCTATTACAATATGATTAGAGTTATTAAGCTGATTATTGATGTaa                                                            <  1:792709/141‑1 (MQ=255)

                                                                                    ggCTGCTAATTCTCTGTATTTTACAGGGGATAAGTAGCCTAGTTTTTGTTGAATTCGATTATTATTATAGTTTTTAATTTA‑TTTTCGACAATATCTATTACAATATGATTAGAGTTATTATGCTGATTATTGATGttacaa                                                         >  4:20812/1‑136 (MQ=255)

                                                                                    ggCTGCTAATTCTCTGTATTTTACAGGGGATAAGTAGCCTAGTTTTTGTTGAATTCGATTATTATTATAGTTTTTAATGTA‑TTTTCGACAATATCTATTACAATATGATTAGAGTTATTAAGCTGATTATTGATGTaaaaa                                                         <  3:197772/141‑1 (MQ=255)

                                                                                     gctgctAATTCTCTGTATTTTACAGGGGATAAGTAGCCTAGTTTTTGTTGAATTCGATTATTATTATAGTTTTTAATGTA‑TTTTCGACAATATCTATTACAATATGATTAGAGTTATTAAGCTGATTATTGATGTAAAAAg                                                        <  4:71731/141‑1 (MQ=255)

                                                                                        gctAATTCTCTGTATTTTACAGGGGATAAGTAGCCTAGTTTTTGTTTAATTTGATTATTATTATAGTTTTTAATGTA‑TTTTCGACAATATCTATTACAATATGATTAGAGTTATTAAGCTGATTATTGATGTAAAAAGttt                                                     >  2:223239/1‑141 (MQ=255)

                                                                                           aaTTCTCTGTATTTTACAGGGGATAAGTAGCCTAGTTTTTGTTGAATTCGATTATTATTATAGTTTTTAATGTA‑TTTTCGACAATATCTATTACAATATGATTAGAGTTATTAAGCTGATTATTGATGTAAAAAGTTTCAg                                                  >  1:137642/1‑141 (MQ=255)

                                                                                            aTTCTCTGTATTTTACAGGGGATAAGTAGCCTAGTTTTTGTTGAATTCGATTATTATTATAGTTTTTAATGTA‑TTTTCGACAATATCTATTACAATATGATTAGAGTTATTAAGCTGATTATTGATGTAAAAAGTTTCAGa                                                 >  4:34690/1‑141 (MQ=255)

                                                                                            aTTCTCTGTATTTTACAGGGGATAAGTAGCCTAGTTTTTGTTGAATTCGATTATTATTATAGTTTTTAATGTA‑TTTTCGACAATACCTATTACAATATGATTAGAGTTATTAAGCTGATTATTGATGTAAAAAGTTTCAGa                                                 >  1:977307/1‑141 (MQ=255)

                                                                                             ttCTCTGTATTTTACAGGGGATAAGTAGCCTAGTTTTTTTTTTTTTTTTTTATTATTATAGTTTTTAATGTA‑TTTTCGACAATATCTATTACAATATGATTAGAGTTATTACGCTGATTATTGATTTACAAAGTTTCTGAc                                                >  4:308113/1‑141 (MQ=18)

                                                                                                    tATTTTACAGGGGATAAGTAGCCTAGTTTTTGTTGAATTCGATTATTATTATAGTTTTTAAAGTA‑TTTTCGACAATATCTATAACAATATGATTAGAGTTATTAAGCTGATTATTGATGTAAAAAGTATCAGACTTTAGCg                                         >  3:56648/1‑141 (MQ=255)

                                                                                                     aTTTTACAGGGGATAAGTAGCCTAGTTTTTGTTGAATTCGATTATTATTATAGTTTTTAATGTA‑TTTTCGACAATATCTATTACAATATGATTAGAGTTATTAAGCTGATTATTGATGTAAAAAGTTTCAGACTTTAGCGa                                        >  2:795883/1‑141 (MQ=255)

                                                                                                      ttttACAGGGGATAAGTAGCCTAGTTTTTGTTGAATTCGATTATTATTATAGTTTTTAATGTA‑TTTTCGACAATATCTATTACAATATGATTAGAGTTATTAAGCTGATTATTGATGTAAAAAGTTTCAGACTTTAGCGAg                                       >  1:569589/1‑141 (MQ=255)

                                                                                                      ttttACAGGGGATAAGTAGCCTAGTTTTTGTTGAATTCGATTATTATTATAGTTTTTAATGTA‑TTTTCGACAATATCTATTACAATATGATTAGAGTTATTAAGCTGATTATTGATGTAAAAAGTTTCAGACTTTAGCGAg                                       >  1:637193/1‑141 (MQ=255)

                                                                                                       tttACAGGGGATAAGTAGCCTAGTTTTTGTTGAATTCGATTATTATTATAGTTTTTAATGTA‑TTTTCGACAATATCTATTACAATATGATTAGAGTTATTAAGCTGATTATTGATGTAAAAAGTTTCAGACTTTAGCGAgg                                      <  1:72432/141‑1 (MQ=255)

                                                                                                        ttACAGGGGATAAGTAGCCTAGTTTTTGTTGAATTCGATTATTATTATAGTTTTTAATGTA‑TTTTCGACAATATCTATTACAATATGATTAGAGTTATTAAGCTGATTATTGATGTAAAAAGTTTCAGACTTTAGCGAGGa                                     >  1:241555/1‑141 (MQ=255)

                                                                                                         tACAGGGGATAAGTAGCCTAGTTTTTGTTGAATTCGATTATTATTATAGTTTTTAATGTA‑TTTTCGACAATATCTATTACAATATGATTAGAGTTATTAAGCTGATTATTGATGTAAAAAGTTTCAGACTTTAGCGAGGaa                                    >  1:26165/1‑141 (MQ=255)

                                                                                                          aCAGGGGATAAGTAGCCTAGTTTTTGTTGAATTCGATTATTATTATAGTTTTTAATGTA‑TTTTCGACAATATCTATTACAATATGATTAGAGTTATTAAGCTGATTATTGATGTAAAAAGTTTCAGACTTTAGCGAGGAAt                                   <  2:354166/141‑1 (MQ=255)

                                                                                                           cAGGGGATAAGTAGCCTAGTTTTTGTTGAATTCGATTATTATTATAGTTTTTAATGTA‑TTTTCGACAATATCTATTACAATATGATTAGAGTTATTAAGCTGATTATTGATGTAAAAAGTTTCAGACTTTAGCGAGGAATg                                  >  1:843486/1‑141 (MQ=255)

                                                                                                               ggATAAGTAGCCTAGTTTTTGTTGAATTCGATTATTATTATAGTTTTTAATGTA‑TTTTCGACAATATCTATTACAATATGATTAGAGTTATTAAGCTGATTATTGATGTAAAAAGTTTCAGACTTTAGCGAGGAATGGaaa                              >  3:266812/1‑141 (MQ=255)

                                                                                                                gATAAGTAGCCTAGTTTTTGTTGAATTCGATTATTATTATAGTTTTTAATGTA‑TTTTCGACAATATCTATTACAATATGATTAGAGTTATTAAGCTGATTATTGATGTAAAAAGTTTCAGACTTTAGCGAGGAATGGAAAc                             >  1:980441/1‑141 (MQ=255)

                                                                                                                  tAAGTAGCCTAGTTTTTTTTTTTTTTTATTATTATTATAGTTTTTAATGTA‑TTTTCGACAATATCTATTACAATATGATTAGAGTTATTAAGCTGATTATTGATGTAAAAAGTTTCAGACTTTAGCGAGGAATGGAAACtt                           >  2:138178/1‑141 (MQ=17)

                                                                                                                  tAAGTAGCCTAGTTTTTGTTGAATTCGATTATTATTATAGTTTTTAATGTA‑TTTTCGACAATATCTATTACAATATGATTAGAGTTATTAAGCTGATTATTGATGTAAAAAGTTTCAGACTTTAGCGAGGAATGGAAACtt                           <  2:329899/141‑1 (MQ=255)

                                                                                                                    aGTAGCCTAGTTTTTGTTGAATTCGATTATTATTATAGTTTTTAATGTA‑TTTTCGACAATATCTATTACAATATGATTAGAGTTATTAAGCTGATTATTGATGTAAAAAGTTTCAGACTTTAGCGAGGAATGGAAACTTTc                         <  2:489073/141‑1 (MQ=255)

                                                                                                                       aGCCTAGTTTTTTTTGTTTTTTATAATTATTAAAGTTTTTAATGTA‑TTTTTTACAATATCTATTACAATATGATTAGAGTTATTAAGCTGATTATTGATGTAAAAAGTTTCAGACTTTAGCGAGGAATGGAAACTTTCTAt                      >  2:930369/1‑141 (MQ=17)

                                                                                                                        gCCTAGTTTTTGTTGAATTCGATTATTATTATAGTTTTTAATGTA‑TTTTCGACAATATCTATTACAATATGATTAGAGTTATTAAGCTGATTATTGATGTAAAAAGTTTCAGACTTTAGCGAGGAATGGAAACTTTCTATc                     >  3:139695/1‑141 (MQ=255)

                                                                                                                         ccTAGTTTTTGTTGAATTTGATTATTATTATAGTTTTTAATGTACTTTTCGACAATATCTATTACAATATGATTAGAGCTATTAAGCTGATTATTGATGTAAAGAGTTTCAGACTTTAGCGAGGAATGGAAACTTTCTATc                     >  1:734034/1‑141 (MQ=17)

                                                                                                                          cTAGTTTTTGTTGAATTCGATTATTATTATAGTTTTTAATGTA‑TTTTCGACAATATCTATTACAATATGATTAGAGTTATTAAGCTGATTATTGATGTAAAAAGTTTCAGACTTTAGCGAGGAATGGAAACTTTCTATCgg                   <  3:130107/141‑1 (MQ=255)

                                                                                                                             gTTTTTGTTGAATTCGATTATTATTATAGTTTTTAATGTA‑TTTTCGACAATATCTATTACAATATGATTAGAGTTATTAAGCTGATTATTGATGTAAAAAGTTTCAGACTTTAGCGAGGAATGGAAACTTTCTATCGGGGc                <  1:890927/141‑1 (MQ=255)

                                                                                                                               ttttGTTGAATTCGATTATTATTATAGTTTTTAATGTA‑TTTTCGACAATATCTATTACAATATGATTAGAGTTATTAAGCTGATTATTGATGTAAAAAGTTTCAGACTTTAGCGAGGAATGGAAACTTTCTATCGGGGCGt              <  1:169643/141‑1 (MQ=255)

                                                                                                                                tttGTTGAATTCGATTATTATTATAGTTTTTAATGTA‑TTTTCGACAATATCTATTACAATATGATTAGAGTTATTAAGCTGATTATTGATGTAAAAAGTTTCAGACTTTAGCGAGGAATGGAAACTTTCTATCGGGGCGtt             >  1:297220/1‑141 (MQ=255)

                                                                                                                                tttGTTGAATTCGATTATTATTATAGTTTTTAATGTA‑TTTTCGACAATATCTATTACAATATGATTAGAGTTATTAAGCTGATTATTGATGTAAAAAGTTTCAGACTTTAGCGAGGAATGGAAACTTTCTATCGGGGCGtt             >  1:135794/1‑141 (MQ=255)

                                                                                                                                   gttgAATTCGATTATTATTATAGTTTTTAATGTA‑TTTTCGACAATATCTATTACAATATGATTAGAGTTATTAAGCTGATTATTGATGTAAAAAGTTTCAGACTTTAGCGAGGAATGGAAACTTTCTATCGGGGCGTTATc          <  1:835059/141‑1 (MQ=255)

                                                                                                                                   gttgAATTCGATTATTATTATAGTTTTTAATGTA‑TTTTCGACAATATCTATTACAATATGATTAGAGTTATTAAGCTGATTATTGATGTAAAAAGTTTCAGACTTTAGCGAGGAATGGAAACTTTCTATCGGGGCGTTATc          <  2:966885/141‑1 (MQ=255)

                                                                                                                                    ttgAATTCGATTATTATTATAGTTTTTAATGTA‑TTTTCGACAATATCTATTACAATATGATTAGAGTTATTAAGCTGATTATTGATGTAAAAAGTTTCAGACTTTAGCGAGGAATGGAAACTTTCTATCGGGGCGTTATCg         >  4:127353/1‑141 (MQ=255)

                                                                                                                                      gAATTCGATTATTATTATAGTTTTTAATGTA‑TTTTCGACAATATCTATTACAATATGATTAGAGTTATTAAGCTGATTATTGATGTAAAAAGTTTCAGACTTTAGCGAGGACTGGAAACTTTCTATCGGGGCGTTATCGGc       >  1:326759/1‑141 (MQ=255)

                                                                                                                                       aaTTCGATTATTATTATAGTTTTTAATGTA‑TTTTCGACAATATCTATTACAATATGATTAGAGTTATTAAGCTGATTATTGATGTAAAAAGTTTCAGACTTTAGCGAGGAATGGAAACTTTCTATCGGGGCGTTATCGGCa      <  4:123508/141‑1 (MQ=255)

                                                                                                                                        aTTCGATTATTATTATAGTTTTTAATGTA‑TTTTCGACAATATCTATTACAATATGATTAGAGTTATTAAGCTGATTATTGATGTAAAAAGTTTCAGACTTTAGCGAGGAATGGAAACTTTCTATCGGTGCGTTATCGGCAg     >  4:302430/1‑141 (MQ=255)

                                                                                                                                           cGATTATTATTATAGTTTTTAATGTA‑TTTTCGACAATATCTATTACAATATGATTAGAGTTATTAAGCTGATTATTGATGTAAAAAGTTTCAGACTTTAGCGAGGAATGGAAACTTTCTATCGGGGCGTTATCGGCAGgtg  >  2:240385/1‑141 (MQ=255)

ATATTATTTTCACAAAAAAAGAATGCTTTTTGTTAGAGGCACTGAACCCTTAAAATGGGAACTTAATAAAAACACTATGTTCTAGGCTGCTAATTCTCTGTATTTTACAGGGGATAAGTAGCCTAGTTTTTGTTGAATTCGATTATTATTATAGTTTTTAATGTA‑TTTTCGACAATATCTATTACAATATGATTAGAGTTATTAAGCTGATTATTGATGTAAAAAGTTTCAGACTTTAGCGAGGAATGGAAACTTTCTATCGGGGCGTTATCGGCAGGTG  >  USA300TCH1516_ALE/61602‑61881

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