BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	USA300TCH1516_ALE	2,610,079	(T)_7→8	coding (447/795 nt)	USA300HOU_RS13455 ←	putative lipoprotein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	USA300TCH1516_ALE	2,610,072	1	.	T	80.4%	287.5 / 65.5	97	coding (454/795 nt)	USA300HOU_RS13455	putative lipoprotein
Reads supporting (aligned to +/- strand): ref base . (4/15); new base T (29/49); total (33/64)
Fisher's exact test for biased strand distribution p-value = 2.80e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.82e-01

Only 100 of 116 total aligned reads displayed.

GACATTATAATCATCATTATTCAATTGATATTTTGCCGAATAACTTGGTACATTAGGATTATATGAAATATCACCATCTTTATAATCATTAATATCTTTAAAGTTGCCATATTGTACAAAGAACTTAAAGTTTTCAATCTC‑‑TTTTTTTAACTTGTCATTCGGTAGTGGCTTCGTTGGTATGATTTTATTATGTTCCATTTTTACCGGATACCTTTTTTCATCATCCTTTGGTCTGCCGTTACTGTCATCTGTCATTTCACTTATAAAATAATACCC  >  USA300TCH1516_ALE/2609932‑2610207

 aCATTATAATCATCATTATTCAATTGATATTTTGCCGAATAACTTGGTACATTAGGATTATATGAAATATCACCATCTTTATAATCATTAATATCTTTAAAGTTGCCATATTGTACAAAGAACTTAAAGTTTTCAATctc‑‑t                                                                                                                                        <  3:401261/141‑1 (MQ=255)

 aCATTATAATCATCATTATTCAATTGATATTTTGCCGAATAACTTGGTACATTAGGATTATATGAAATATCACCATCTTTATAATCATTAATATCTTTAAAGTTGCCATATTGTACAAAGAACTTAAAGTTTTCAATctc‑‑t                                                                                                                                        <  1:711239/141‑1 (MQ=255)

 aCATTATAATCATCATTATTCAATTGATATTTTGCCGAATAACTTGGTACATTAGGATTATATGAAATATCACCATCTTTATAATCATTAATATCTTTAAAGTTGCCATATTGTACAAAGAACTTAAAGTTTTCAATctc‑‑t                                                                                                                                        <  2:17880/141‑1 (MQ=255)

  cATTATAATCATCATTATTCAATTGATATTTTGCCGAATAACTTGGTACATTAGGATTATATGAAATATCACCATCTTTATAATCATTAATATCTTTAAAGTTGCCATATTGTACAAAGAACTTAAAGTTTTCAATCTC‑‑tt                                                                                                                                       >  3:129186/1‑141 (MQ=255)

  cATTATAATCATCATTATTCAATTGATATTTTGCCGAATAACTTGGTACATTAGGATTATATGAAATATCACCATCTTTATAATCATTAATATCTTTAAAGTTGCCATATTGTACAAAGAACTTAAAGTTTTCAATCTC‑‑tt                                                                                                                                       >  1:305529/1‑141 (MQ=255)

   aTTATAATCATAATTATTCGATTGATATTTTGCCGAATAACTTGGTACATTAGGATTATATGAAATATCACCATCTTTATAATCATTAATATCTTTAAAGTTGCCATATTGTACAAAGAACTTAAAGTTTTCAATCTC‑‑ttt                                                                                                                                      <  2:916078/141‑1 (MQ=255)

    ttATAATCATCATTATTCAATTGATATTTTGCCGAATAACTTGGTACATTAGGATTATATGAAATATCACCATCTTTATAATCATTAATATCTTTAAAGTTGCCATATTGTACAAAGAACTTAAAGTTTTCAATCTC‑‑tttt                                                                                                                                     >  2:883417/1‑141 (MQ=255)

    ttATAATCATCATTATTCAATTGATATTTTGCCGAATAACTTGGTACATTAGGATTATATGAAATATCACCATCTTTATAATCATTAATATCTTTAAAGTTGCCATATTGTACAAAGAACTTAAAGTTTTCAATCTC‑‑tttt                                                                                                                                     >  1:818730/1‑141 (MQ=255)

    ttATAATCATCATTATTCAATTGATATTTTGCCGAATAACTTGGTACATTAGGATTATATGAAATATCACCATCTTTATAATCATTAATATCTTTAAAGTTGCCATATTGTACAAAGAACTTAAAGTCTCCAATCTC‑‑tttt                                                                                                                                     >  1:652397/1‑141 (MQ=255)

      ataATCATCATTATTCAATTGATATTTTGCCGAATAACTTGGTACATTAGGATTATATGAAATATCACCATCTTTATAATCATTAATATCTTTAAAGTTGCCATATTGTACAAAGAACTTAAAGTTTTCAATCTC‑‑tttttt                                                                                                                                   >  2:42140/1‑141 (MQ=255)

        aaTCATCATTATTCAATTGATATTTTGCCGAATAACTTGGTACATTAGGATTATATGAAATATCACCATCTTTATAATCATTAATATCTTTAAAGTTGCCATATTGTACAAAGAACTTAAAGTTTTCAATCTC‑‑tttttttt                                                                                                                                 <  1:560595/141‑2 (MQ=255)

        aaTCATCATTATTCAATTGATATTTTGCCGAATAACTTGGTACATTAGGATTATATGAAATATCACCATCTTTATAATCATTAATATCTTTAAAGTTGCCATATTGTACAAAGAACTTAAAGTTTTCAATCTC‑‑tttttttt                                                                                                                                 <  1:560640/141‑2 (MQ=255)

        aaTCATCATTATTCAATTGATATTTTGCCGAATAACTTGGTACATTAGGATTATATGAAATATCACCATCTTTATAATCATTAATATCTTTAAAGTTGCCATATTGTACAAAGAACTTAAAGTTTTCAATCTC‑‑tttttttt                                                                                                                                 >  1:282096/1‑140 (MQ=255)

        aaTCATCATTATTCAATTGATATTTTGCCGAATAACTTGGTACATTAGGATTATATGAAATATCACCATCTTTATAATCATTAATATCTTTAAAGTTGCCATATTGTACAAAGAACTTAAAGTTTTCAATCTC‑‑tttttttt                                                                                                                                 >  2:641875/1‑140 (MQ=255)

        aaTCATCATTATTCAATTGATATTTTGCCGAATAACTTGGTACATTAGGATTATATGAAATATCACCATCTTTATAATCATTAATATCTTTAAAGTTGCCATATTGTACAAAGAACTTAAAGTTTTCAATCTC‑‑tttttttt                                                                                                                                 >  2:914640/1‑140 (MQ=255)

         atcatcATTATTCAATTGATATTTTGCCGAATAACTTGGTACATTAGGATTATATGAAATATCACCATCTTTATAATCATTAATATCTTTAAAGTTGCCATATTGTACAAAGAACTTAAAGTTTTCAATCTC‑‑tttttttta                                                                                                                                >  2:112491/1‑139 (MQ=255)

           catcatTATTCAATTGATATTTTGCCGAATAACTTGGTACATTAGGATTATATGAAATATCACCATCTTTATAATCATTAATATCTTTAAAGTTGCCATATTGTACAAAGAACTTAAAGTTTTCAATCTC‑‑ttttttttaac                                                                                                                              <  2:517752/141‑5 (MQ=255)

           catcatTATTCAATTGATATTTTGCCGAATAACTTGGTACATTAGGATTATATGAAATATCACCATCTTTATAATCATTAATATCTTTAAAGTTGCCATATTGTACAAAGAACTTAAAGTTTTCAATCTC‑‑ttttttttaac                                                                                                                              <  2:355897/141‑5 (MQ=255)

              catTATTCAATTGATATTTTGCCGAATAACTTGGTACATTAGGATTATATGAAATATCACCATCTTTATAATCATTAATATCTTTAAAGTTGCCATATTGTACAAAGAACTTAAAGTTTTCAATCTCT‑TTTTTTTAACTTg                                                                                                                            >  1:521862/1‑141 (MQ=255)

                 tattCAATTGATATTTTGCCGAATAACTTGGTACATTAGGATTATATGAAATATCACCATCTTTATAATCATTAATATCTTTAAAGTTGCCATATTGTACAAAGAACTTAAAGTTTTCAATCTCT‑TTTTTTTAACTTGTCa                                                                                                                         <  3:443134/141‑1 (MQ=255)

                  tttCAATTGATATTTTGCCGAATAACTTGGTACATTAGGATTATATGAAATATCACCATCTTTATAATCATTAATATCTTTAAAGTTGCCATATTGTACAAAGAACTTAAAGTTTTCAATCTCTTTTTTTTTAACTTGTCa                                                                                                                         <  2:538710/140‑1 (MQ=255)

                  attCAATTGATATTTTGCCGAATAACTTGGTACATTAGGATTATATGAAATATCACCATCTTTATAATCATTAATATCTTTAAAGTTGCCATATTGTACAAAGAACTTAAAGTTTTCAATCTCT‑TTTTTTTAACTTGTCAt                                                                                                                        <  3:285014/141‑1 (MQ=255)

                       aTTGATATTTTGCCGAATAACTTGGTACATTAGGATTATATGAAATATCACCATCTTTATAATCATTAATATCTTTAAAGTTGCCATATTGTACAAAGAACTTAAAGTTTTCAATCTCT‑TTTTTTTAACTTGTCATTCGGt                                                                                                                   <  4:211122/141‑1 (MQ=255)

                           atatTTTGCCGAATAACTTGGTACATTAGGATTATATGAAATATCACCATCTTTATAATCATTAATATCTTTAAAGTTGCCATATTGTACAAAGAACTTAAAGTTTTCAATCTCT‑TTTTTTTAACTTGTCATTCGGTAGTg                                                                                                               <  1:1009618/141‑1 (MQ=255)

                           atatTTTGCCGAATAACTTGGTACATTAGGATTATATGAAATATCACCATCTTTATAATCATTAATATCTTTAAAGTTGCCATATTGTACAAAGAACTTAAAGTTTTCAATCTC‑‑TTTTTTTAACTTGTCATTCGGTAGTgg                                                                                                              <  3:162135/141‑1 (MQ=255)

                             atTTTGCCGAATAACTTGGTACATTAGGATTATATGAAATATCACCATCTTTATAATCATTAATATCTTTAAAGTTGCCATATTGTACAAAGAACTTAAAGTTTTCAATCTCT‑TTTTTTTAACTTGTCATTCGGTAGTGGc                                                                                                             >  1:937/1‑141 (MQ=255)

                             atTTTGCCGAATAACTTGGTACATTAGGATTATATGAAATATCACCATCTTTATAATCATTAATATCTTTAAAGTTGCCATATTGTACAAAGAACTTAAAGTTTTCAATCTCT‑TTTTTTTAACTTGTCATTCGGTAGTGGc                                                                                                             >  1:2602/1‑141 (MQ=255)

                              ttttGCCGAATAACTTGGTACATTAGGATTATATGAAATATCACCATCTTTATAATCATTAATATCTTTAAAGTTGCCATATTGTACAAAGAACTTAAAGTTTTCAATCTCTTTTTTTTTAACTTGTCATTCGGTAGTGGc                                                                                                             <  1:734589/141‑1 (MQ=255)

                               tttGCCGAATAACTTGGTACATTAGGATTATATGAAATATCACCATCTTTATAATCATTAATATCTTTAAAGTTGCCATATTGTACAAAGAACTTAAAGTTTTCAATCTC‑‑TTTTTTTAACTTGTCATTCGGTAGTGGCTTc                                                                                                          >  1:932966/1‑141 (MQ=255)

                                ttGCCGAATAACTTGGTACATTAGGATTATATGAAATATCACCATCTTTATAATCATTAATATCTTTAAAGTTGCCATATTGTACAAAGAACTTAAAGTTTTCAATCTCT‑TTTTTTTAACTTGTCATTCGGTAGTGGCTTc                                                                                                          >  1:993256/1‑141 (MQ=255)

                                   ccGAATAACTTGGTACATTAGGATTATATGAAATATCACCATCTTTATAATCATTAATATCTTTAAAGTTGCCATATTGTACAAAGAACTTAAAGTTTTCAATCTCT‑TTTTTTTAACTTGTCATTCGGTAGTGGCTTCGtt                                                                                                       >  1:209734/1‑141 (MQ=255)

                                      aaTAACTTGGTACATTAGGATTATATGAAATATCACCATCTTTATAATCATTAATATCTTTAAAGTTGCCATATTGTACAAAGAACTTAAAGTTTTCAATCTCT‑TTTTTTTAACTTGTCATTCGGTAGTGGCTTCGTTGGt                                                                                                    <  2:752337/141‑1 (MQ=255)

                                       aTAACTTGGTACATTAGGATTATATGAAATATCACCATCTTTATAATCATTAATATCTTTAAAGTTGCCATATTGTACAAAGAACTTAAAGTTTTCAATCTCT‑TTTTTTTAACTTGTCATTCGGTAGTGGCTTCGTTGGTa                                                                                                   >  3:48742/1‑141 (MQ=255)

                                        tAACTTGGTACATTAGGATTATATGAAATATCACCATCTTTATAATCATTAATATCTTTAAAGTTGCCATATTGTACAAAGAACTTAAAGTTTTCAATCTCT‑TTTTTTTAACTTGTCATTCGGTAGTGGCTTCGTTGGTAt                                                                                                  <  1:421799/141‑1 (MQ=255)

                                        tAACTTGGTACATTAGGATTATATGAAAAATCACCATCTTTATAATCATTAATATCTTTAAAGTTGCCATATTGTACAAAGAACTTAAAGTTTTCAATCTCT‑TTTTTTTAANTTGTCATTCGGTNGTGGCTTCGTTGGTAt                                                                                                  <  1:875300/141‑1 (MQ=255)

                                         aaCTTGGTACATTAGGATTATATGAAATATCACCATCTTTATAATCATTAATATCTTTAAAGTTGCCATATTGTACAAAGAACTTAAAGTTTTCAATCTCT‑TTTTTTTAACTTGTCATTCGGTAGTGGCTTCGTTGGTATg                                                                                                 >  2:132191/1‑141 (MQ=255)

                                         aaCTTGGTACAATAGGATTATATGAAATATCAGCATCTTTATAATCATTAATATCTTTAAAGCTGCCATATTGTACAAAGAACTTAAAGTTTTCAATCTCT‑TTTTTTTAACTTGTCATTTGGAAGTGGCTTCGTTGGGAtt                                                                                                 >  2:375980/1‑140 (MQ=17)

                                            ttGGTACATTAGGATTATATGAAATATCACCATCTTTATAATCATTAATATCTTTAAAGTTGCCATATTGTACAAAGAACTTAAAGTTTTCAATCTC‑‑TTTTTTTAACTTGTCATTCGGTAGTGGCTTCGTTGGTATGAttt                                                                                             <  2:548716/141‑1 (MQ=255)

                                             tGGTACATTAGGATTATATGAAATATCACCATCTTTATAATCATTAATATCTTTAAAGTTGCCATATTGTACAAAGAACTTAAAGTTTTCAATCTCT‑TTTTTTTAACTTGTCATTCGGTAGTGGCTTCg                                                                                                         <  2:832687/129‑1 (MQ=255)

                                             tGGTACATTAGGATTATATGAAATATCACCATCTTTATAATCATTAATATCTTTAAAGTTGCCATATTGTACAAAGAACTTAAAGTTTTCAATCTCT‑TTTTTTTAACTTGTCATTCGGTAGTGGCTTCg                                                                                                         >  1:832687/1‑129 (MQ=255)

                                              ggTACATTAGGATTATATGAAATATCACCATCTTTATAATCATTAATATCTTTAAAGTTGCCATATTGTACAAAGAACTTAAAGTTTTCAATCTCT‑TTTTTTTAACTTGTCATTCGGTAGTGGCTTCGTTGGTATGCtttt                                                                                            >  1:789110/1‑141 (MQ=255)

                                               gTACATTAGGATTATATGAAATATCACCATCTTTATAATCATTAATATCTTTAAAGTTGCCATATTGTACAAAGAACTTAAAGTTTTCAATCTCT‑TTTTTTTAACTTGTCATTCGGTAGTGGCTTCGTTGGTATGATTtta                                                                                           <  3:1073/141‑1 (MQ=255)

                                                   aTTAGGATTATATGAAATATCACCATCTTTATAATCATTAATATCTTTAAAGTTGCCATATTGTACAAAGAACTTAAAGTTTTCAATCTCT‑TTTTTTTAACTTGTCATTCGGTAGTGGCTTCGTTGGTATGATTTtattat                                                                                       <  4:178622/141‑1 (MQ=255)

                                                   aTTAGGATTATATGAAATATCACCATCTTTATAATCATTAATATCTTTAAAGTTGCCATATTGTACAAAGAACTTAAAGTTTTCAATCTCT‑TTTTTTTAACTTGTCATTCGGTAGTGGCTTCGTTGGTATGATTTtattat                                                                                       <  2:397788/141‑1 (MQ=255)

                                                     tAGGATTATATGAAATATCACCATCTTTATAATCATTAATATCTTTAAAGTTGCCATATTGTACAAAGAACTTAAAGTTTTCAATCTC‑‑TTTTTTTAACTTGTCATTCGGTAGTGGCTTCGTTGGTATGATTTTATTATGtt                                                                                    <  4:392288/141‑1 (MQ=255)

                                                        gATTATATGAAATATCACCATCCTTATAATCATTAATATCTTTAAAGTTGCCATATTGTACAAAGAACTTAAAGTTTTCAATCTCT‑TTTTTTTAACTTGTCATTCGGTAGTGGCTTCGTTGGTATGATTTTATTATGTTcc                                                                                  <  1:689909/141‑1 (MQ=255)

                                                            atatGAAATATCACCATCTTTATAATCATTAATATCTTTAAAGTTGCCATATTGTACAAAGAACTTAAAGTTTTCAATCTC‑‑TTTTTTTAACTTGTCATTCGGTAGTGGCTTCGTTGGTATGATTTTATTATGTTCCAtttt                                                                             <  1:78713/141‑1 (MQ=255)

                                                             tatGGATTATTACCATCTTTATAATCATTTATATCTTTAAAGTTGCCATATTGTACAAAGTACTGATAGTTTTCAATCTCT‑TTTTTTTTATTTGTCATTCGGTAGTGGCTTCGTTGGTATGAATTTATTATGTCACAtttc                                                                             <  2:727525/141‑2 (MQ=18)

                                                             tatGAAATATCACCATCTTTATAATCATTAATATCTTTAAAGTTGCCATATTGTACAAAGAACTTAAAGTTTTCAATCTCT‑TTTTTTTAACTTGTCATTCGGTAGTGGCTTCGTTGGTATGATTTTATTATGTTCCAtttt                                                                             <  1:436088/141‑1 (MQ=255)

                                                             tatGAAATATCACCATCTTTATAATCATTAATATCTTTAAAGTTGCCATATTGTACAAAGAACTTAAAGTTTTCAATCTCT‑TTTTTTTAACTTGTCATTCGGTAGTGGCTTCGTTGGTATGATTTTATTATGTTCCAtttt                                                                             <  2:424952/141‑1 (MQ=255)

                                                                   atatCACCATCTTTATAATCATTAATATCTTTAAAGTTGCCATATTGTACAAAGAACTTAAAGTTTTCAATCTC‑‑TTTTTTTAACTTGTCATTCGGTAGTGGCTTCGTTGGTATGATTTTATTATGTTCCATTTTTACCGGa                                                                      <  1:907634/141‑1 (MQ=255)

                                                                      tCACCATCTTTATAATCATTAATATCTTTAAAGTTGCCATATTGTACAAAGAACTTAAAGTTTTCAATCTCT‑TTTTTTTTTTTTTTCATTCGGTAGTGGCTTCGTTGGTATGATTTTATTATGTTCCATTTTTACCGGATa                                                                    >  4:2471/1‑141 (MQ=25)

                                                                      tCACCATCTTTATAATCATTAATATCTTTAAAGTTGCCATATTGTACAAAGAACTTAAAGTTTTCAATCTCT‑TTTTTTTAACTTGTCATTCGGTAGTGGCTTCGTTGGTATGATTTTATTATGTTCCATTTTTACCGGATa                                                                    >  3:180296/1‑141 (MQ=255)

                                                                      tCACCATCTTTATAATCATTAATATCTTTAAAGTTGCCATATTGTACAAAGAACTTAAAGTTTTCAATCTC‑‑TTTTTTTAACTTGTCATTCGGTAGTGGCTTCGTTGGTATGATTTTATTATGTTCCATTTTTACCGGATAc                                                                   <  2:932966/141‑1 (MQ=255)

                                                                       cACCATCTTTATAATCATTAATATCTTTAAAGTTGCCATATTGTACAAAGAACTTAAAGTTTTCAATCTCT‑TTTTTTTAACTTGTCATTCGGTAGTGGCTTCGTTGGTATGATTTTATTATGTTCCATTTTTACCGGATAc                                                                   >  2:355086/1‑141 (MQ=255)

                                                                        aCCATCTTTATAATCATTAATATCTTTAAAGTTGCCATATTGTACAAAGAACTTAAAGTTTTCAATCTC‑‑TTTTTTTAACTTGTCATTCGGTAGTGGCTTCGTTGGTATGATTTTATTATGTTCCATTTTTACCGGATACCt                                                                 <  1:353252/141‑1 (MQ=255)

                                                                           aTCTTTATAATCATTAATATCTTTAAAGTTGCCATATTGTACAAAGAACTTAAAGTTTTCAATCTCT‑TTTTTTTAACTTGTCATTCGGTAGTGGCTTCGTTGGTATGATTTTATTATGTTCCATTTTTACCGGATACCttt                                                               <  2:362390/141‑1 (MQ=255)

                                                                            tCTTTATCATCATTAATATCTTTAAAGTTGCCATATTGTACAAAGAACTTAAAGTTTTCAATCTCT‑TTTTTTTAACTTGTCATTCGGTAGTGGCTTCGTTTGTATGATTTTATTATGTTCCATTTTTGCCGGACACCtttt                                                              <  4:197003/141‑1 (MQ=21)

                                                                             cTTTATAATCATTAATATCTTTAAAGTTGCCATATTGTACAAAGAACTTAAAGTTTTCAATCTCT‑TTTTTTTAACTTGTCATTCGGTAGTGGCTTCGTTGGTATGATTTTATTATGTTCCATTTTTACCGGATACCttttt                                                             >  2:88448/1‑141 (MQ=255)

                                                                              tttATAATCATTAATATCTTTAAAGTTGCCATATTGTACAAAGAACTTATAGTTTTCAATCTCT‑TTTTTTTAATTTGTAATTTGGTAGTGGCTTCGTTTGTATGATTTTATTATTTTCCATTTTTTCTTGGTACCtttttt                                                            >  2:129718/1‑141 (MQ=18)

                                                                               ttATAATCATTAATATCTTTAAAGTTGCCATATTGTACAAAGAACTTAAAGTTTTCAATCTCT‑TTTTTTTAACTTGTCATTCGGTAGTGGCTTCGTTGGTATGATTTTATTATGTTCCATTTTTACCGGATACCttttttt                                                           >  1:197890/1‑140 (MQ=255)

                                                                                  taATCATTAATATCTTTAAAGTTGCCATATTGTACAAAGAACTTAAAGTTTTCAATCTCT‑TTTTTTTAACTTGTCATTCGGTAGTGGCTTCGTTGGTATGATTTTATTATGTTCCATTTTTACCGGATACCTTTTTtcatc                                                        >  3:46094/1‑141 (MQ=255)

                                                                                       aTTAATATCTTTAAAGTTGCCATATTGTACAAAGAACTTAAAGTTTTCAATCTC‑‑TTTTTTTAACTTGTCATTCGGTAGTGGCTTCGTTGGTATGATTTTATTATGTTCCATTTTTACCGGATTCCTTTTTTCATCATCCtt                                                  <  2:553296/141‑1 (MQ=255)

                                                                                         tAATATCTTTAAAATAGCTATTTTGTACAAAGATCTTAAAGTTTTCAATCTCT‑TTTTTTTAACTTGTCATTCGGTAGTGGCTTCGTTCGTATGATTTTATTATGTTCCATTTTTACCGGATACCTTTTTTCATCATCCttt                                                 <  1:844264/141‑1 (MQ=25)

                                                                                                tttAAAGTTGCCATATTGTACAAAGAACTTAAAGTTTTCAATCTCT‑TTTTTTTAACTTGTCATTCGGTAGTGGCTTCGTTGGTATGATTTTATTATGTTCCATTTTTACCGGATACCTTTTTTCATCATCCTTTGGTCTGc                                          >  1:272671/1‑141 (MQ=255)

                                                                                                 ttAAAGTTGCCATATTGTACAAAGAACTTAAAGTTTTCAATCTCT‑TTTTTTTAACTTGTCATTCGGTAGTGGCTTCGTTGGTATGATTTTATTATGTTCCATTTTTACCGGATACCTTTTTTCATCATCCTTTGGTCTGcc                                         >  1:482452/1‑141 (MQ=255)

                                                                                                 ttAAAGTTGCCATATTGTACAAAGAACTTAAAGTTTTCAATCTCT‑TTTTTTTAACTTGTCATTCGGTAGTGGCTTCGTTGGTATGATTTTATTATGTTCCATTTTTACCGGATACCTTTTTTCATCATCCTTTGGTCTGcc                                         >  1:482457/1‑141 (MQ=255)

                                                                                                     aGTTGCCATATTGTACAAAGAACTTAAAGTTTTCAATCTCT‑TTTTTTTAACTTGTCATTCGGTAGTGGCTTCGTTGGTATGATTTTATTATGTTCCATTTTTACCGGATACCTTTTTTCATCATCCTTTGGTCTGCCGTTa                                     <  1:410396/141‑1 (MQ=255)

                                                                                                     aGTTGCCATATTGTACAAAGAACTTAAAGTTTTCAATCTCT‑TTTTTTTAACTTGTCATTCGGTAGTGGCTTCGTTGGTATGATTTTATTATGTTCCATTTTTACCGGATACCTTTTTTCATCATCCTTTGGTCTGCCGTTa                                     <  4:381248/141‑1 (MQ=255)

                                                                                                     aGTTGCCATATTGTACAAAGAACTTAAAGTTTTCAATCTCT‑TTTTTTTAACTTGTCATTCGGTAGTGGCTTCGTTGGTATGATTTTATTATGTTCCATTTTTACCGGATACCTTTTTTCATCATCCTTTGGTCTGCCGTTa                                     <  1:440363/141‑1 (MQ=255)

                                                                                                     aGTTGCCATATTGTACAAAGAACTTAAAGTTTTCAATCTCT‑TTTTTTTAACTTGTCATTCGGTAGTGGCTTCGTTGGTATGATTTTATTATGTTCCATTTTTACCGGATACCTTTTTTCATCATCCTTTGGTCTGCCGTTa                                     <  1:986866/141‑1 (MQ=255)

                                                                                                     aGTTGCCATATTGTACAAAGAACTTAAAGTTTTCAATCTCT‑TTTTTTTAACTTGTCATTCGGTAGTGGCTTCGTTGGTATGATTTTATTATGTTCCATTTTTACCGGATACCTTTTTTCATCATCCTTTGGTCTGCCGTTa                                     <  1:986817/141‑1 (MQ=255)

                                                                                                     aGTTGCCATATTGTACAAAGAACTTAAAGTTTTCAATCTCT‑TTTTTTTAACTTGTCATTCGGTAGTGGCTTCGTTGGTATGATTTTATTATGTTCCATTTTTACCGGATACCTTTTTTCATCATCCTTTGGTCTGCCGTTa                                     <  1:892851/141‑1 (MQ=255)

                                                                                                      gTTGCCATATTGTACAAAGAACTTAAAGTTTTCAATCTC‑‑‑TTTTTTAACTTGTCATTCGGTAGTGGCTTCGTTGGTATGATTTTATTATGTTCCATTTTTACCGGATACCTTTTTTCATCATCCTTTGGTCTGCCGTTActg                                  <  4:316643/141‑1 (MQ=255)

                                                                                                       ttGCCATATTGTACAAAGAACTTAAAGTTTTCAATCTC‑‑TTTTTTTAACTTGTCATTCGGTAGTGGCTTCGTTGGTATGATTTTATTATGTTCCATTTTTACCGGATACCTTTTTTCATCATCCTTTGGTCTGCCGTTActg                                  <  3:410258/141‑1 (MQ=255)

                                                                                                          ccATATTGTACAAAGAACTTAAAGTTTTCAATCTC‑‑TTTTTTTAACTTGTCATTCGGTAGTGGCTTCGTTGGTATGATTTTATTATGTTCCATTTTTACCGGATACCTTTTTTCATCATCCTTTGGTCTGCCGTTActgtca                               >  1:441983/1‑141 (MQ=255)

                                                                                                           cATATTGTACAAAGAACTTAAAGTTTTCAATCTCT‑TTTTTTTAACTTGTCATTCGGTAGTGGCTTCGTTGGTATGATTTTATTATGTTCCATTTTTACCGGATACCTTTTTTCATCATCCTTTGGTCTGCCGTTActgtca                               <  2:484084/141‑1 (MQ=255)

                                                                                                              atTGTACAAAGAACTTAAAGTTTTCAATCTCT‑TTTTTTTAACTTGTCATTCGGTAGTGGCTTCGTTGGTATGATTTTATTATGTTCCATTTTTACCGGATACCTTTTTTCATCATCCTTTGGTCTGCCGTTActgtcatct                            <  1:1022126/141‑1 (MQ=255)

                                                                                                              atTGTACAAAGAACTTAAAGTTTTCAATCTCT‑TTTTTTTAACTTGTCATTCGGTAGTGGCTTCGTTGGTATGATTTTATTATGTTCCATTTTTACCGGATACCTTTTTTCATCATCCTTTGGTCTGCCGTTActgtcatct                            <  1:328028/141‑1 (MQ=255)

                                                                                                              atTGTACAAAGAACTTAAAGTTTTCAATCTC‑‑TTTTTTTAACTTGTCATTCGGTAGTGGCTTCGTTGGTATGATTTTATTATGTTCCATTTTTACCGGATACCTTTTTTCATCATCCTTTGGTCTGCCGTTActgtcatctg                           <  2:1001464/141‑1 (MQ=255)

                                                                                                               ttGTACAAAGAACTTAAAGTTTTCAATCTCT‑TTTTTTTAACTTGTCATTCGGTAGTGGCTTCGTTGGTATGATTTTATTATGTTCCATTTTTACCGGATACGTTTTTTCATCATCCTTTGGTCTGCCGTTActgtcatctg                           <  4:375463/141‑1 (MQ=255)

                                                                                                               ttGTACAAAGAACTTAAAGTTTTCAATCTCT‑TTTTTTTAACTTGTCATTCGGTAGTGGCTTCGTTGGTATGATTTTATTATGTTCCATTTTTACCGGATACCTTTTTTCATCATCCTTTGGTCTGCCGTTActgtcatctg                           <  2:693572/141‑1 (MQ=255)

                                                                                                                   aCAAAGAACTTAAAGTTTTCAATCTCT‑TTTTTTTAACTTGTCATTCGGTAGTGGCTTCGTTGGTATGATTTTATTATGTTCCATTTTTACCGGATACCTTTTTTCATCTTCCTTTGGTCTGCCGTTActgtcatctgtcat                       >  4:138453/1‑141 (MQ=255)

                                                                                                                    caaattaCTTATAGTTTTCAATCTCT‑TTTTTTTAACTTGTCTTTCGGTAGTGGCTTCGTTAGTATGATTTTATTTTGTTCCATTTTTACCGGATACCTTTTGTCATCATCCTTTTGTCTGCCGTTACTGTCATCTGTCAtt                      <  2:80426/135‑1 (MQ=255)

                                                                                                                      aaGAACTTAAAGTTTTCAATCTCT‑TTTTTTTAACTTGTCATTCGGTAGTGGCTTCGTTGGTATGATTTTATTATGTTCCATTTTTACCGGATACCTTTTTTCATCATCCTTTGGTCTGCCGTTACTGTCATCTGTCATTTc                    <  1:305831/141‑1 (MQ=255)

                                                                                                                      aaGAACTTAAAGTTTTCAATCTC‑‑TTTTTTTAACTTGTCATTCGGTAGTGGCTTCGTTGGTATGATTTTATTATGTTCCATTTTTACCGGATACCTTTTTTCATCATCCTTTGGTCTGCCGTTACTGTCATCTGTCATTTCa                   >  1:159279/1‑141 (MQ=255)

                                                                                                                       aGAACTTAAAGTTTTCAATCTCT‑TTTTTTTTACTTTTTTTTCGGTAGTGGCTTCGTTCGTATGATTTTATTTTGTTCCATTTTTACCGGATACTTTTTTTTATCATCGTTTGTTCTGCCTTTACTGTCATCTTTTATTTCa                   >  2:678859/1‑141 (MQ=255)

                                                                                                                         aaCTTAAAGTTTTCAATCTCT‑TTTTTTTAACTTGTCATTCGGTAGTGGCTTCGTTGGTATGATTTTATTATGTTCCATTTTTACCGGATACCTTTTTTCATCATCCTTTGGTCTGCCGTTACTGTCATCTGTCATTTCACt                 <  1:54197/141‑1 (MQ=255)

                                                                                                                              aaaGTTTTCAATCTCT‑TTTTTTTAACTTGTCATTCGGTAGTGGCTTCGTTGGTATGATTTTATTATGTTCCATTTTTACCGGATACCTTTTTTCATCATCCTTTGGTCTGCCGTTACTGTCATCTGTCATTTCACTTATaa            >  1:244867/1‑141 (MQ=255)

                                                                                                                              aaaGTTTTCAATCTCT‑TTTTTTTAACTTGTCATTCGGTAGTGGCTTCGTTGGTATGATTTTATTATGTTCCATTTTTACCGGATACCTTTTTTCATCATCCTTTGGTCTGCCGTTACTGTCATCTGTCATTTCACTTATaa            >  3:144126/1‑141 (MQ=255)

                                                                                                                              aaaGTTTTCAATCTCT‑TTTTTTTAACTTGTCATTCGGTAGTGGCTTCGTTGGTATGATTTTATTATGTTCCATTTTTACCGGATACCTTTTTTCATCATCCTTTGGTCTGCCGTTACTGTCATCTGTCATTTCACTTATaa            >  1:217953/1‑141 (MQ=255)

                                                                                                                              aaaGTTTTCAATCTC‑‑TTTTTTTAACTTGTCATTCGGTAGTGGCTTCGTTGGTATGATTTTATTATGTTCCATTTTTACCGGATACCTTTTTTCATCATCCTTTGGTCTGCCGTTACTGTCATCTGTCATTTCACTTATaaa           <  2:305529/141‑1 (MQ=255)

                                                                                                                               aaGTTTTCAATCTCT‑TTTTTTTAACTTGTCATTCGGTAGTGGCTTCGTTGGTATGATTTTATTATGTTCCATTTTTACCGGATACCTTTTTTCATCATCCTTTGGTCTGCCGTTACTGTCATCTGTCATTTCACTTATaaa           >  3:228640/1‑141 (MQ=255)

                                                                                                                                    ttCAATCTCT‑TTTTTTTAACTTGTCATTCGGTAGTGGCTTCGTTGGTATGATTTTATTATGTTCCATTTTTACCGGATACCTTTTTTCATCATCCTTTGGTCTGCCGTTACTGTCATCTGTCATTTCACTTATAAaataat      >  2:684444/1‑141 (MQ=255)

                                                                                                                                    ttCAATCTCT‑TTTTTTTAACTTGTCATTCGGTAGTGGCTTCGTTGGTATGATTTTATTATGTTCCATTTTTACCGGATACCTTTTTTCATCATCCTTTGGTCTGCCGTTACTGTCATCTGTCATTTCACTTATAAaataat      <  3:83517/141‑1 (MQ=255)

                                                                                                                                     tCAATCTCT‑TTTTTTTAATTTGTCATTCGGTAGTGGCTTCGTTGGTATGATTTTATTATGTTCCATTTTTACCGGATACCTTTTTTCATCATCCTTTGGTCTGCCGTTACTTTCATTTTTCATTTCACTTTTAAAataata     >  4:309275/1‑141 (MQ=255)

                                                                                                                                      cAATCTC‑‑TTTTTTTAACTTGTCATTCGGTAGTGGCTTCGTTGGTATGATTTTATTATGTTCCATTTTTACCGGATACCTTTTTTCATCATCCTTTGGTCTGCCGTTACTGTCATCTGTCATTTCACTTATAAAATAATAcc   <  4:398951/141‑1 (MQ=255)

                                                                                                                                       aaTCTCT‑TTTTTTTAACTTGTCATTCGGTAGTGGCTTCGTTGGTATGATTTTATTATGTTCCATTTTTACCGGATACCTTTTTTCATCATCCTTTGGTCTGCCGTTACTGTCATCTGTCATTTCACTTATAAAATAATAcc   <  2:248090/141‑1 (MQ=255)

                                                                                                                                       aaTCTCT‑TTTTTTTAACTTGTCATTCGGTAGTGGCTTCGTTGGTATGATTTTATTATGTTCCATTTTTACCGGATACCTTTTTTCATCATCCTTTGGTCTGCCGTTACTGTCATCTGTCATTTCACTTATAAAATAATAcc   <  2:721871/141‑1 (MQ=255)

                                                                                                                                       aaTCTC‑‑TTTTTTTAACTTGTCATTCGGTAGTGGCTTCGTTGGTATGATTTTATTATGTTCCATTTTTACCGGATACCTTTTTTCATCATCCTTTGGTCTGCCGTTACTGTCATCTGTCATTTCACTTATAAAATAATAccc  >  1:939527/1‑141 (MQ=255)

GACATTATAATCATCATTATTCAATTGATATTTTGCCGAATAACTTGGTACATTAGGATTATATGAAATATCACCATCTTTATAATCATTAATATCTTTAAAGTTGCCATATTGTACAAAGAACTTAAAGTTTTCAATCTC‑‑TTTTTTTAACTTGTCATTCGGTAGTGGCTTCGTTGGTATGATTTTATTATGTTCCATTTTTACCGGATACCTTTTTTCATCATCCTTTGGTCTGCCGTTACTGTCATCTGTCATTTCACTTATAAAATAATACCC  >  USA300TCH1516_ALE/2609932‑2610207

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