BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	USA300TCH1516_ALE	2,721,513	A→G	7.0%	N209D (AAT→GAT)	USA300HOU_RS14000 →	putative hydrolase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	USA300TCH1516_ALE	2,721,513	0	A	G	7.0%	423.6 / 15.0	144	N209D (AAT→GAT)	USA300HOU_RS14000	putative hydrolase
Reads supporting (aligned to +/- strand): ref base A (70/63); new base G (3/7); total (74/70)
Fisher's exact test for biased strand distribution p-value = 2.02e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.83e-01

Only 100 of 152 total aligned reads displayed.

CATTTGGCGAAACATTAAATATCGCACCAATTGATGCAAAAATGATGTCTCAACCTGCTGATACTGAAGAAGGACGTATAGAACAGTACAAACGTACGATGTTCTGGTTAGAATTTGAAATTAGACAATATACGCATTCAAATATTACTTTGGATGATTTCACAAAATACAGCGATAAAATCACCCTACTAAATGGGACAGATTCTAGAGGCTCCTTCCCACAAGATGTTAACTTTTATATAAATAAAGAAACTGGCATACCTATTGTAGATATTCCAGG  >  USA300TCH1516_ALE/2721373‑2721652

 aTTTGGCGAAACATTAAATATCGCACCAATTGATGCAAAAATGATGTCTCAACCTGCTGATACTGAAGAAGGACGTATAGAACAGTACAACCGTACGATGTTCTGGTTAGAATTTGAAATTAGACAATATACGCATTCaaa                                                                                                                                            <  1:39094/141‑1 (MQ=255)

 aTTTGGCGAAACATTAAATATCGCACCAATTGATGCAAAAATGATGTCTCAACCTGCTGATACTGAAGAAGGACGTATAGAACAGTACAAACGTACGATGTTCTGGTTAGAATTTGAAATTAGACAATATACGCATTCaaa                                                                                                                                            >  4:383620/1‑141 (MQ=255)

 aTTTGGCGAAACATTAAATATCGCACCAATTGATGCAAAAATGATGTCTCAACCTGCTGATACTGAAGAAGGACGTATAGAACAGTACAAACGTACGATGTTCTGGTTAGAATTTGAAATTAGACAATATACGCATTCaaa                                                                                                                                            >  1:779737/1‑141 (MQ=255)

 aTTTGGCGAAACATTAAATATCGCACCAATTGATGCAAAAATGATGTCTCAACCTGCTGATACTGAAGAAGGACGTATAGAACAGTACAAACGTACGATGTTCTGGTTAGAATTTGAAATTAGACAATATACGCATTCaaa                                                                                                                                            <  1:877568/141‑1 (MQ=255)

  ttttGCGAAACATTAAATATCGCACCAATTGATGCAAAAATGATGTCTCAACCTGCTGATACTGAAGAAGGACGTATAGAACAGTACAAACGTACGATGTTCTGGTTAGAATTTGAAATTAGACAATATACGCATTCAAat                                                                                                                                           <  4:337113/141‑1 (MQ=255)

  tttGGCGAAACATTAAATATCGCACCAATTGATGCAAAAATGATGTCTCAACCTGCTGATACTGAAGAAGGACGTATAGAACAGTACAAACGTACGATGTTCTGGTTAGAATTTGAAATTAGACAATATACGCATTCAAat                                                                                                                                           >  2:21451/1‑141 (MQ=255)

   ttGGCGAAACATTAAATATCGCACCAATTGATGCAAAAATGATGTCTCAACCTGCTGATACTGAAGAAGGACGTATAGAACAGTACAAACGTACGATGTTCTGGTTAGAATTTGAAATTAGACAATATACGCATTCAAata                                                                                                                                          <  2:912212/141‑1 (MQ=255)

    tGGCGAAACATTAAATATCGCAGCAATTGATGCAAAAATGATGTCTCAACCTGCTGATACTGAAGAAGGACGTATAGAACAGTACAAACGTACGATGTTCTGGTTAGAATTTGAAATTAGACAATATACGCATTCAAatat                                                                                                                                         <  4:324369/141‑1 (MQ=255)

    tGGCGAAACATTAAATATCGCACCAATTGATGCAAAAATGATGTCTCAACCTGCTGATACTGAAGAAGGACGTATAGAACAGTACAAACGTACGATGTTCTGGTTAGAATTTGAAATTAGACAATATACGCATTCAAatat                                                                                                                                         >  1:878320/1‑141 (MQ=255)

     ggCGAAACATTAAATATCGCACCAATTGATGCAAAAATGATGTCTCAACCTGCTGATACTGAAGAAGGACGTATAGAACAGTACAAACGTACGATGTTCTGGTTAGAATTTGAAATTAGACAATATACGCATTCAAATAtt                                                                                                                                        <  2:800633/141‑1 (MQ=255)

     ggCGAAACATTAAATATCGCACCAATTGATGCAAAAATGATGTCTCAACCTGCTGATACTGAAGAAGGACGTATAGAACAGTACAAACGTACGATGTTCTGGTTAGAATTTGAAATTAGACAATATACGCATTCAAATAtt                                                                                                                                        >  4:172963/1‑141 (MQ=255)

           aCATTAAATATCGCACCAATTGATGCAAAAATGATGTCTCAACCTGCTGATACTGAAGAAGGACGTATAGAACAGTACAAACGTACGATGTTCTGGTTAGAATTTGAAATTAGACAATATACGCATTCAAATATTACTTTg                                                                                                                                  >  2:844075/1‑141 (MQ=255)

              ttAAATATCGCACCAATTGATGCAAAAATGATGTCTCAACCTGCTGATACTGAAGAAGGACGTATAGAACAGTACAAACGTACGATGTTCTGGTTAGAATTTGAAATTAGACAATATACGCATTCAAATATTACTTTGgat                                                                                                                               >  3:395284/1‑141 (MQ=255)

                aaaTATCGCACCAATTGATGCAAAAATGATGTCTCAACCTGCTGATACTGAAGAAGGACGTATAGAACAGTACAAACGTACGATGTTCTGGTTAGAATTTGAAATTAGACAATATACGCATTCAAATATTACTTTGgatga                                                                                                                             >  1:838853/1‑141 (MQ=255)

                aaaTATCGCACCAATTGATGCAAAAATGATGTCTCAACCTGCTGATACTGAAGAAGGACGTATAGAACAGTACAAACGTACGATGTTCTGGTTAGAATTTGAAATTAGACAATATACGCATTCAAATATTACTTTGgatga                                                                                                                             <  3:29410/141‑1 (MQ=255)

                 aaTATCGCACCAATTGATGCAAAAATGATGTCTCAACCTGCTGATACTGAAGAAGGACGTATAGAACAGTACAAACGTACGATGTTCTGGTTAGAATTTGAAATTAGACAATATACGCATTCAAATATTACTTTGgatgat                                                                                                                            >  2:819809/1‑141 (MQ=255)

                 aaTATCGCACCAATTGATGCAAAAATGATGTCTCAACCTGCTGATACTGAAGAAGGACGTATAGAACAGTACAAACGTACGATGTTCTGGTTAGAATTTGAAATTAGACAATATACGCATTCAAATATTACTTTGgatgat                                                                                                                            <  2:595421/141‑1 (MQ=255)

                    atCGCACCAATTGATGCAAAAATGATGTCTCAACCTGCTGATACTGAAGAAGGACGTATAGAACAGTACAAACGTACGATGTTCTGGTTAGAATTTGAAATTAGACAATATACGCATTCAAATATTACTTTGGATGATTTc                                                                                                                         >  2:363292/1‑141 (MQ=255)

                     tCGCACCAATTGATGCAAAAATGATGTCTCAACCTGCTGATACTGAAGAAGGACGTATAGAACAGTACAAACGTACGATGTTCTGGTTAGAATTTGAAATTAGACAATATACGCATTCAAATATTACTTTGGATGATTTca                                                                                                                        >  2:752031/1‑141 (MQ=255)

                     tCGCACCAATTGAGGCAAAAATGATGTCTCAACCTGGTGATACTGAAGAAGGACGTATAGAACAGTACAAACGTACGATGTTCTGGTCAGAATTTGAAATTAGACAATATACGCATTCAGATATTACTTTGGATGATTGca                                                                                                                        <  2:122552/141‑1 (MQ=255)

                      cGCACCAATTGATGCAAAAATGATGTCTCAACCTGCTGATACTGAAGAAGGACGTATAGAACAGTACAAACGTACGATGTTCTGGTTAGAATTTGAAATTAGACAATATACGCATTCAAATATTACTTTGGATGATTTcac                                                                                                                       <  2:1081226/141‑1 (MQ=255)

                         aCCAATTGATGCAAAAATGATGTCTCAACCTGCTGATACTGAAGAAGGACGTATAGAACAGTACAAACGTACGATGTTCTGGTTAGAATTTGAAATTAGACAATATACGCATTCAAATATTACTTTGGATGATTTCACaaa                                                                                                                    >  1:304755/1‑141 (MQ=255)

                          ccAATTGATGCAAAAATGATGTCTCAACCTGCTGATACTGAAGAAGGACGTATAGAACAGTACAAACGTACGATGTTATGGTTAGAATTTGAAATTAGACAATATACGCATTCAAATATTACTTTGGATGATTTCACaaaa                                                                                                                   >  2:1015493/1‑141 (MQ=255)

                                 aTGCAAAAATGATGTCTCAACCTGCTGATACTGAAGAAGGACGTATAGAACAGTACAAACGTACGATGTTCTGGTTAGAATTTGAAATTAGACAATATACGCATTCAAATATTACTTTGGATGATTTCACaaaa                                                                                                                   >  1:643313/1‑134 (MQ=255)

                                 aTGCAAAAATGATGTCTCAACCTGCTGATACTGAAGAAGGACGTATAGAACAGTACAAACGTACGATGTTCTGGTTAGAATTTGAAATTAGACAATATACGCATTCAAATATTACTTTGGATGATTTCACaaaa                                                                                                                   <  2:643313/134‑1 (MQ=255)

                                  tGCAAAAATGATGTCTCAACCTGCTGATACTGAAGAAGGACGTATAGAACAGTCCAAACGTACGATGTTCTGGTCAGATTTAGAAATTAGACAATATACGCATTCAGATATTACTTTAGATGATATCACAAAATACAGCGa                                                                                                           >  2:411321/1‑141 (MQ=255)

                                   ccacacatgacgTCCCACCCTGCTGATACTGAAGAAGGACGTGTAGAACAGTACAAACGTACGATGTTCTGGTTAGAATTTGAAATTAGCCAATATACGCATTCAAATATAACTTTGGATGATTTCACAAAATACAGCGAt                                                                                                          <  4:173120/135‑1 (MQ=255)

                                    cAAAAATGATGTCTCAACCTGCTGATACTGAAGAAGGACGTATAGAACAGTACAAACGTACGATGTTCTGGTTAGAATTTGAAATTAGACAATATACGCATTCAAATATTACTTTGGATGATTTCACAAAATACAGCGATa                                                                                                         <  2:763634/141‑1 (MQ=255)

                                             tgTCTCAACCTGCTGATACTGAAGAAGGACGTATAGAACAGTACAAACGTACGATGTTCTGGTTAGAATTTGAAATTAGACAATATACGCATTCAAATATTACTTTGGATGATTTCACAAAATACAGCGATAAAATCAccc                                                                                                <  2:555564/141‑1 (MQ=255)

                                              gTCTCAACCTGCTGATACTGAAGAAGGACGTATAGAACAGTACAAACGTACGATGTTCTGGTTAGAATTTGAAATTAGACAATATACGCATTCAAATATTACTTTGGATGATTTCACAAAATACAGCGATAAAATCACCct                                                                                               <  3:47909/141‑1 (MQ=255)

                                               tctcAACCTGCTGATACTGAAGAAGGACGTATAGAACAGTACAAACGTACGATGTTCTGGTTAGAATTTGAAATTAGACAATATACGCATTCAAATATTACTTTGGATGATTTCACAAAATACAGCGATAAAATCACCcta                                                                                              <  1:81137/141‑1 (MQ=255)

                                                ctcAACCTGCTGATACTGAAGAAGGACGTATAGAACAGTACAAACGTACGATGTTCTGGTTAGAATTTGAAATTAGACAATATACGCATTCAAATATTACTTTGGATGATTTCTCAAAATACAGCGATAAAATCTCCctac                                                                                             <  4:267768/141‑1 (MQ=255)

                                                ctcAACCTGCTGATACTGAAGAAGGACGTATAGAACAGTACAAACGTACGATGTTCTGGTTAGAATTTGAAATTAGACAATATACGCATTCAAATATTACTTTGGATGATTTCACAAAATACAGCGATAAAATCACCctac                                                                                             <  4:207888/141‑1 (MQ=255)

                                                ctcAACCGGCTGATACTGAAGAAGGACGTATAGAACAGTACAAACGTACGATGTTCTGGTTAGAATTTGAAATTAGACAATATACGCATTCAAATATTACTTTGGATGATTTCACAAAATACAGCGATAAAATCACCctac                                                                                             <  4:268285/141‑1 (MQ=255)

                                                   aaCCTGCTGATACTGAAGAAGGACGTATAGAACAGTACAAACGTACGATGTTCTGGTTAGAATTTGAAATTAGACAATATACGCATTCAAATATTACTTTGGATGATTTCACAAAATACAGCGATAAAATCACCCTACTaa                                                                                          >  1:75692/1‑141 (MQ=255)

                                                    aCCTGCTGATACTGAAGAAGGACGTATAGAACAGTACAAACGTACGATGTTCTGGTTAGAATTTGAAATTAGACAATATACGCATTCAAATATTACTTTGGATGATTTCACAAAATACAGCGATAAAATCACCCTACTaaa                                                                                         >  1:771320/1‑141 (MQ=255)

                                                       tgctgATACTGAAGAAGGACGTATAGAACAGTACAAACGTACGATGTTCTGGTTAGAATTTGAAATTAGACAATATACGCATTCAAATATTACTTTGGATGATTTCACAAAATACAGCGATAAAATCACCCTACTAAATgg                                                                                      <  3:171739/141‑1 (MQ=255)

                                                       tgctgATACTGAAGAAGGACGTATAGAACAGTACAAACGTACGATGTTCTGGTTAGAATTTGAAATTAGACAATATACGCATTCAAATATTACTTTGGATGATTTCACAAAATACAGCGATAAAATCACCCTACTAAATgg                                                                                      <  4:420632/141‑1 (MQ=255)

                                                        gctgATACTGAAGAAGGACGTATAGAACAGTACAAACGTACGATGTTCTGGTTAGAATTTGAAATTAGACAATATACGCATTCAAATATTACTTTGGATGATTTCACAAAATACAGCGATAAAATCACCCTACGAAATggg                                                                                     >  1:1038353/1‑141 (MQ=255)

                                                         ctgATACTGAAGAAGGACGTATAGAACAGTACAAACGTACGATGTTCTGGTTAGAATTTGAAATTAGACAATATACGCATTCAAATATTACTTTGGATGATTTCACAAAATACAGCGATAAAATCACCCTACTAAATGGGa                                                                                    <  2:696629/141‑1 (MQ=255)

                                                           gATACTGAAGAAGGACGTATAGAACAGTACAAACGTACGATGTTCTGGTTAGAATTTGAAATTAGACAATATACGCATTCAAATATTACTTTGAAAAAAAAAAAAAAATACAGCGATAAAATCACCCTACTAAATGGGACa                                                                                  <  2:593154/141‑1 (MQ=255)

                                                           gATACTGAAGAAGGACGTATAGAACAGTACAAACGTACGATGTTCTGGTCAGAATTTGAAATTAGACAATATACGCATTCAGATATTACTTTGGATGATTTCACAAAATACAGCGATAAAATCACCCTACTAAATGGGACa                                                                                  <  2:386567/141‑1 (MQ=255)

                                                            agaCTGAAGAAGGACGCATAGAACAGTACAAACGTACGATGTTCTGGTTAGAATTTGAAATTAGACAATATACGCATTCAAATATTACTTTGGATGATTTCACAAAATACAGCGATAAAATCACCCTACTAAATGGGACAg                                                                                 >  4:81770/3‑141 (MQ=255)

                                                            aTACTGAAGAAGGACGTATAGAACAGTACAAACGTACGATGTTCTGGTTAGAATTTGAAATTAGACAATATACGCATTCAAATATTACTTTGGATGATTTCACAAAATACAGCGATAAAATCACCCTACTAAATGGGACAg                                                                                 <  2:96335/141‑1 (MQ=255)

                                                                tGAAGAAGGACGTATAGAACAGTACAAACGTACGATGTTCTGGTTAGAATTTGAAATTAGACAATATACGCATTCAAATATTACTTTGGATGATTTCACAAAATACAGCGATAAAATCACCCTACTAAATGGGACAGATTc                                                                             <  4:104188/141‑1 (MQ=255)

                                                                  aagaagGACGTATAGAACAGTACAAACGTACGATGTTCTGGTTAGAATTTGAAATTAGACAATATACGCATTCAAATATTACTTTGGATGATTTCACAAAATACAGCGATAAAATCACCCTACTAAATGGGACAGATTCTa                                                                           >  4:252175/1‑141 (MQ=255)

                                                                   agaagGACGTATAGAACAGTACAAACGTACGATGTTCTGGTTAGAATTTGAAATTAGACAATATACGCATTCAAATATTACTTTGGATGATTTCACAAAATACAGCGATAAAATCACCCTACTAAATGGGACAGATTCTag                                                                          >  1:25598/1‑141 (MQ=255)

                                                                   agaagGACGTATAGAACAGTACAAACGTACGATGTTCTGGTCAGAATTTGAAATTAGACAATATACGCATTCAGATATTACTTTGGATGATTTCACAAAATACAGCGATAAAATCACCCTACTAAATGGGACAGATTCTag                                                                          <  2:125603/141‑1 (MQ=255)

                                                                    gaagGACGTATAGAACAGTACAAACGTACGATGTTCTGGTTAGAATTTGAAATTAGACAATATACGCATTCAAATATTACTTTGGATGATTTCACAAAATACAGCGATAAAATCACCCTACTAAATGGGACAGATTCTaga                                                                         >  1:309292/1‑141 (MQ=255)

                                                                    gaagGACGTATAGAACAGTACAAACGTACGATGTTCTGGTTAGAATTTGAAATTAGACAATATACGCATTCAAATATTACTTTGGATGATTTCACAAAATACAGCGATAAAATCACCCTACTAAATGGGACAGATTCTaga                                                                         >  1:846434/1‑141 (MQ=255)

                                                                     aagGACGTATAGAACAGTACAAACGTACGATGTTCTGGTTAGAATTTGAAATTAGACAATATACGCATTCAAATATTACTTTGGATGATTTCACAAAATACAGCGATAAAATCACCCTACTAAATGGGACAGATTCTagag                                                                        <  3:61962/141‑1 (MQ=255)

                                                                      agGACGTATAGAACAGTACAAACGTACGATGTTCTGGTTAGAATTTGAAATTAGACAATATACGCATTCAAATATTACTTTGGATGATTTCACAAAATACAGCGATAAAATCACCCTACTAAATGGGACAGATTCTAGAgg                                                                       >  1:722942/1‑141 (MQ=255)

                                                                       ggACGTATAGAACAGTACAAACGTACGATGTTCTGGTTAGAATTTGAAATTAGACAATATACGCATTCAAATATTACTTTGGATGATTTCACAAAATACAGCGATAAAATCACCCTACTAAATGGGACAGATTCTAGAGGc                                                                      <  1:369191/141‑1 (MQ=255)

                                                                        gACGTATAGAACAGTACAAACGTACGATGTTCTGGTTAGAATTTGAAATTAGACAATATACGCATTCAAATATTACTTTGGATGATTTCACAAAATACAGCGATAAAATCACCCTACTAAATGGGACAGATTCTAGAGGCt                                                                     <  3:155431/141‑1 (MQ=255)

                                                                        gACGTATAGAACAGTACAAACGTACGATGTTCTGGTTAGAATTTGAAATTAGACAATATACGCATTCAAATATTACTTTGGATGATTTCACAAAATACAGCGATAAAATCACCCTACTAAATGGGACAGATTCTAGAGGCt                                                                     <  2:171287/141‑1 (MQ=255)

                                                                         aCGTATAGAACAGTACAAACGTACGATGTTCTGGTTAGAATTTGAAATTAGACAATATACGTATTCAAATATTACTTTGGATGATTTCACAAAATACAGCGATAAAATCACCCTACTAAATGGGACAGATTCTAGAGGCTc                                                                    >  4:423088/1‑141 (MQ=255)

                                                                         aCGTATAGAACAGTACAAACGTACGATGTTCTGGTTAGAATTTGAAATTAGACAATATACGCATTCAAATATTACTTTGGATGATTTCACAAAATACAGCGATAAAATCACCCTACTAAATGGGACAGATTCTAGAGGCTc                                                                    >  4:28287/1‑141 (MQ=255)

                                                                             ataGATCAGGACAAAGGTATGATTTTCTTGTTAGAATTTGAAATTAGACAATATACGCATTCAAATATTACTTTGGATGATTTCACAAAATACAGCGATAAAATCACCCTACTAAATGGGACAGATTCTAGAGGCTCCTTc                                                                <  2:973063/141‑1 (MQ=255)

                                                                             ataGAACAGTACAAACGTACGATGTTCTGGTTAGAATTTGAAATTAGACAATATACGCATTCAAATATTACTTTGGATGATTTCACAAAATACAGCGATAAAATCACCCTACTAAATGGGACAGATTCTAGAGGCTCCTTc                                                                <  1:675065/141‑1 (MQ=255)

                                                                             ataGAACAGTACAAACGTACGATGTTCTGGTTAGAATTTGAAATTAGACAATATACGCATTCAAATATTACTTTGGATGATTTCACAAAATACAGCGATAAAATCACCCTACTAAATGGGACAGATTCTAGAGGCTCCTTc                                                                >  2:880354/1‑141 (MQ=255)

                                                                                gAACAGTACAAACGTACGATGTTCTGGTTAGAATTTGAAATTAGACAATATACGCATTCAAATATTACTTTGGATGATTTCACAAAATACAGCGATAAAATCACCCTACTAAATGGGACAGATTCTAGAGGCTCCTTCCca                                                             >  1:243751/1‑141 (MQ=255)

                                                                                gAACAGTACAAACGTACGATGTTCTGGTTAGAATTTGAAATTAGACAATATACGCATTCAAATATTACTTTGGATGATTTCACAAAATACAGCGATAAAATCACCCTACTAAATGGGACAGATTCTAGAGGCTCCTTCCca                                                             <  2:967242/141‑1 (MQ=255)

                                                                                 aaCAGTACAAACGTACGATGTTCTGGTTAGAATTTGAAATTAGACAATATACGCATTCAAATATTACTTTGGATGATTTCACAAAATACAGCGATAAAATCACCCTACTAAATGGGACAGATTCTAGAGGCTCCTTCCcac                                                            >  4:351609/1‑141 (MQ=255)

                                                                                       aCAAACGTACGATGTTCTGGTTAGAATTTGAAATTAGACAATATACGCATTCAAATATTACTTTGGATGATTTCACAAAATACAGCGATAAAATCACCCTACTAAATGGGACAGATTCTAGAGGCTCCTTCCCACAAGATg                                                      >  2:286770/1‑141 (MQ=255)

                                                                                       aCAAACGTACGATGTTCTGGTCAGAATTTGAAATTAGACAATATACGCATTCAGATATTACTTTGGATGATTTCACAAAATACAGCGATAAAATCACCCTACTAAATGGGACAGATTCTAGAGGCTCCTTCCCACAAGATg                                                      <  2:366879/141‑1 (MQ=255)

                                                                                        cAAACGTACGATGTTCTGGTTATAATTTCAAATTAGACAATATACGCATTCAAATATTACTTTGGATGATTTCACAAAATACAGCGATAAAATCACCCTACTAAATGGTACAGATTCTAGATGCTCCTTCCCACAAGATGt                                                     <  1:1015493/141‑1 (MQ=255)

                                                                                        cAAACGTACGATGTTCTGGTTAGAATTTGAAATTAGACAATATACGCATTCAAATATTACTTTGGATGATTTCACAAAATACAGCGATAAAATCACCCTACTAAATGGGACAGATTCTAGAGGCTCCTTCCCACAAGATGt                                                     <  4:187565/141‑1 (MQ=255)

                                                                                        cAAACGTACGATGTTCTGGTCAGAATTTGAAATTAGACAATATACGCATTCAGATATTACTTTGGATGATTTCACAAAATACAGCGATAAAATCACCCTACTAAATGGGACAGATTCTAGAGGCTCCTTCCCACAAGATGt                                                     <  1:476618/141‑1 (MQ=255)

                                                                                         aaaCGTACGATGTTCTGGTTAGAATTTGAAATTAGACAATATACGCATTCAAATCTTACTTTGGATGATTTCACAAAATACAGCGATAAAATCACCCTACTAAATGGGACAGATTCTAGAGGCTCCTTCCCACAAGATGtt                                                    <  3:327131/141‑1 (MQ=255)

                                                                                          aaCGTACGATGTTCTGGTTAGAATTTGAAATTAGACAATATACGCATTCAAATATTACTTTGGATGATTTCACAAAATACAGCGATAAAATCACCCTATTACCTGGGACAGTTTCTAGAGGCTCCTTCCCACAAGATGTTa                                                   >  3:431537/1‑141 (MQ=255)

                                                                                          aaCGTACGATGTTCTGGTTAGAATTTGAAATTAGACAATATACGCATTCAAATATTACTTTGGATGATTTCACAAAATACAGCGATAAAATCACCCTACTAAATGGGACAGATTCTAGAGGCTCCTTCCCACAAGATGTTa                                                   <  1:142275/141‑1 (MQ=255)

                                                                                          aaCGTACGATGTTCTGGTTAGAATTTGAAATTAGACAATATACGCATTCAAATATTACTTTGGATGATTTCACAAAATACAGCGATAAAATCACCCTACTAAATGGGACAGATTCTAGAGGCTCCTTCCCACAAGATGTTa                                                   >  2:62187/1‑141 (MQ=255)

                                                                                          aaCGTACGATGTTCTGGTTAGAATTTGAAATTAGACAATATACGCATTCAAATATTACTTTGGATGATTTCACAAAATACAGCGATAAAATCACCCTACTAAATGGGACAGATTCTAGAGGCTCCTTCCCACAAGATGTTa                                                   >  1:894221/1‑141 (MQ=255)

                                                                                            cGTACGATGTTCTGGTTAGAATTTGAAATTAGACAATATACGCATTCAAATATTACTTTGGATGATTTCACAAAATACAGCGATAAAATCACCCTACTAAATGGGACAGATTCTAGAGGCTCCTTCCCACAAGATGTTAAc                                                 >  2:1063641/1‑141 (MQ=255)

                                                                                            cGTACGATGTTCTGGTTAGAATTTGAAATTAGACAATATACGCATTCAAATATTACTTTGGATGATTTCACAAAATACAGCGATAAAATCACCCTACTAAATGGGACAGATTCTAGAGGCTCCTTCCCACAAGATGTTAAc                                                 >  2:109403/1‑141 (MQ=255)

                                                                                             gTACGATGTTCTGGTTAGAATTTGAAATTAGACAATATACGCATTCAAATATTACTTTGGATGATTTCACAAAATACAGCGATAAAATCACCCTACTAAATGGGACAGATTCTAGAGGCTCCTTCCCACAAGATGTTAACt                                                >  1:515096/1‑141 (MQ=255)

                                                                                               aCGATGTTCTGGTTAGAATTTGAAATTAGACAATATACGCATTCAAATATTACTTTGGATGATTTCACAAAATACAGCGATAAAATCACCCTACTAAATGGGACAGATTCTAGAGGCTCCTTCCCACAAGATGTTAACttt                                              >  1:135728/1‑141 (MQ=255)

                                                                                                    gTTCTGGTTAGAATTTGAAATTAGACAATATACGCATTCAAATATTACTTTGGATGATTTCACAAAATACAGCGATAAAATCACCCTACTAAATGGGACAGATTCTAGAGGCTCCTTCCCACAAGATGTTAACTTTtatat                                         >  1:171892/1‑141 (MQ=255)

                                                                                                      tCTGGTTAGAATTTGAAATTAGACAATATACGCATTCAAATATTACTTTGGATGATTTCACAAAATACAGCGATAAAATCACCCTACTAAATGGGACAGATTCTAGAGGCTCCTTCCCACAAGATGTTAACTTTTATataa                                       >  3:57925/1‑141 (MQ=255)

                                                                                                      tCTGGTTAGAATTTGAAATTAGACAATATACGCATTCAAATATTACTTTGGATGATTTCACAAAATACAGCGATAAAATCACCCTACTAAATGGGACAGATTCTAGAGGCTCCTTCCCACAAGATGTTAACTTTTATataa                                       <  2:425434/141‑1 (MQ=255)

                                                                                                       cTGGTTAGAATTTGAAATTAGACAATATACGCATTCAAATATTACTTTGGATGATTTCACAAAATACAGCGATAAAATCACCCTACTAAATGGGACAGATTCTAGAGGCTCCTTCCCACAAGATGTTAACTTTTATataaa                                      <  3:383620/141‑1 (MQ=255)

                                                                                                        tGGTTAGAATTTGAAATTAGACAATATACGCATTCAAATATTACTTTGGATGATTTCACAAAATACAGCGATAAAATCACCCTACTAAATGGGACAGATTCTAGAGGCTCCTTCCCACAAGATGTTAACTTTTATataaat                                     >  2:315654/1‑141 (MQ=255)

                                                                                                         ggTTAGAATTTGAAATTAGACAATATACGCATTCAAATATTACTTTGGATCATTTCACCAAATACAGCGATAAAATCACCCTACTAAATGGGACAGATTCTAGAGGCTCCTTCCCACAAGATGTTAACTTTTATataaata                                    <  3:42044/141‑1 (MQ=255)

                                                                                                             aGAATTTGAAATTAGACAATATACGCATTCAAATATTACTTTGGATGATTTCACAAAATACAGCGATAAAATCACCCTACTAAATGGGACAGATTCTAGAGGCTCCTTCCCACAAGATGTTAACTTTTATATAAATAAAGa                                >  1:597748/1‑141 (MQ=255)

                                                                                                                aTTTGAAATTAGACAATATACGCATTCAAATATTACTTTGGATGATTTCACAAAATACAGCGATAAAATCACCCTACTAAATGGGACAGATTCTAGAGGCTCCTTCCCACAAGATGTTAACTTTTATATAAATAAAGAAAc                             >  2:315950/1‑141 (MQ=255)

                                                                                                                  ttGAAATTAGACAATATACGCATTCAAATATTACTTTGGATGATTTCACAAAATACAGCGATAAAATCACCCTACTAAATGGGACAGATTCTAGAGGCTCCTTCCCACAAGATGTTAACTTTTATATAAATAAAGAAACTg                           >  2:246848/1‑141 (MQ=255)

                                                                                                                   tGAAATTAGACAATATACGCATTCAAATATTACTTTGGATGATTTCACAAAATACAGCGATAAAATCACCCTACTAAATGGGACAGATTCTAGAGGCTCCTTCCCACAAGATGTTAACTTTTATATAAATAAAGAAACTgg                          >  1:1076737/1‑141 (MQ=255)

                                                                                                                        ttAGACAATATACGCATTCAAATATTACTTTGGATGATTTCACAAAATACAGCGATAAAATCACCCTACTAAATGGGACAGATTCTAGAGGCTCCTTCCCACAAGATGTTAACTTTTATATAAATAAAGAAACTGGCATAc                     <  1:180816/141‑1 (MQ=255)

                                                                                                                        ttAGACAATATACGCATTCAAATATTACTTTGGATGATTTCACAAAATACAGCGATAAAATCACCCTACTAAATGGGACAGATTCTAGAGGCTCCTTCCCACAAGATGTTAACTTTTATATAAATAAAGAAACTGGCATAc                     <  1:920664/141‑1 (MQ=255)

                                                                                                                         tAGACAATATACGCATTCAAATATTACTTTGGATGATTTCACAAAATACAGCGATCATATCACCCTACTAAATGGGACAGATTCTAGAGGCTCCTTCCCACAAGATGTTAACTTTTATATAAATAAAGAAACTGGCATAcc                    <  1:1053306/141‑1 (MQ=255)

                                                                                                                         tAGACAATATACGCATTCAAATATTACTTTGGATGATTTCACAAAATACAGCGATAAAATCACCCTACTAAATGGGACAGATTCTAGAGGCTCCTTCCCACAAGATGTTAACTTTTATATAAATAAAGAAACTGGCATAcc                    <  3:292802/141‑1 (MQ=255)

                                                                                                                          aGACAATATACGCATTCAAATATTACTTTGGATGATTTCACAAAATACAGCGATAAAATCACCCTACTAAATGGGACAGATTCTAGAGGCTCCTTCCCACAAGATGTTAACTTTTATATAAATAAAGAAACTGGCATACCt                   <  2:838853/141‑1 (MQ=255)

                                                                                                                          aGACAATATACGCATTCAAATATTACTTTGGATGATTTCACAAAATACAGCGATAAAATCACCCTACTAAATGGGACAGATTCTAGAGGCTCCTTCCCACAAGATGTTAACTTTTATATAAATAAAGAAACTGGCATACCt                   >  4:372639/1‑141 (MQ=255)

                                                                                                                          aGACAATATACGCATTCAAATATTACTTTGGATGATTTCACAAAAGACAGCGATAAAATCACCCTACTAAATGGGACAGATTCTAGAGGCTCCTTCCCACAAGATGTTAACTTTTATATAAATAAAGAAACTGGCATACCt                   >  3:101260/1‑141 (MQ=255)

                                                                                                                            aCAATATACGCATTCAAATATTACTTTGGATGATTTCACAAAATACAGCGATAAAATCACCCTACTAAATGGGACAGATTCTAGAGGCTCCTTCCCACAAGATGTTAACTTTTATATAAATAAAGAAACTGGCATACCTAt                 <  2:666045/141‑1 (MQ=255)

                                                                                                                             cAATATACGCATTCAAATATTACTTTGGATGATTTCACAAAATACAGCGATAAAATCACCCTACTAAATGGGACAGATTCTAGAGGCTCCTTCCCACAAGATGTTAACTTTTATATAAATAAAGAAACTGGCATACCTAtt                >  2:714794/1‑141 (MQ=255)

                                                                                                                                  taCGCATTCAAATATTACTTTGGATGATTTCACAAAATACAGCGATAAAATCACCCTACTAAATGGGACAGATTCTAGAGGCTCCTTCCCACAAGATGTTAACTTTTATATAAATAAAGAAACTGGCATACCTATTGTAGa           <  2:483989/141‑1 (MQ=255)

                                                                                                                                      cATTCAAATATTACTTTGGATGATTTCACAAAATACAGCGATAAAATCACACTACTAAATGGGACAGATTCTAGAGGCTCCTTCCCACAAGATGTTAACTTTTATATAAATAAAGAAACTGGCATACCTATTGTAGATAtt       <  2:774102/141‑1 (MQ=255)

                                                                                                                                       aTTCAAATATTACTTTGGATGATTTCACAAAATACAGCGATAAAATCACCCTACTAAATGGGACAGATTCTAGAGGCTCCTTCCCACAAGATGTTAACTTTTATATAAATAAAGAAACTGGCATACCTATTGTAGATATTc      >  4:297883/1‑141 (MQ=255)

                                                                                                                                       aTTCAAATATTACTTTGGATGATTTCACAAAATACAGCGATAAAATCACCCTACTAAATGGGACAGATTCTAGAGGCTCCTTCCCACAAGATGTTAACTTTTATATAAATAAAGAAACTGGCATACCTATTGTAGATATTc      >  1:683954/1‑141 (MQ=255)

CATTTGGCGAAACATTAAATATCGCACCAATTGATGCAAAAATGATGTCTCAACCTGCTGATACTGAAGAAGGACGTATAGAACAGTACAAACGTACGATGTTCTGGTTAGAATTTGAAATTAGACAATATACGCATTCAAATATTACTTTGGATGATTTCACAAAATACAGCGATAAAATCACCCTACTAAATGGGACAGATTCTAGAGGCTCCTTCCCACAAGATGTTAACTTTTATATAAATAAAGAAACTGGCATACCTATTGTAGATATTCCAGG  >  USA300TCH1516_ALE/2721373‑2721652

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