BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	USA300TCH1516_ALE	1,331,236	T→A	5.2%	S189T (TCT→ACT)	USA300HOU_RS06555 →	Monoacylglycerol lipase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	USA300TCH1516_ALE	1,331,236	0	T	A	5.2%	332.6 / 6.2	108	S189T (TCT→ACT)	USA300HOU_RS06555	Monoacylglycerol lipase
Reads supporting (aligned to +/- strand): ref base T (56/45); new base A (2/4); total (59/49)
Fisher's exact test for biased strand distribution p-value = 4.09e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 109 total aligned reads displayed.

AAGGTTTACCAACCGTTAAAGTGTTACAACTGATTACAAAAATTTATGGTGCTGAGAAACGAGTTGAATGGGTTAACCAGTTAGTATCAAATAGTTTTAATAAAAAAATACGTCCATTACGTACACAAAGTGATTGGATTTCTAGTAATCCAATTGAAGTAGATAACTTTATTAAAGATCCATATAGTGGATTTAATGTGTCAAATCAATTATTATATCAAACAGCCTATTATATGCTACATACATCACAATTAAAAAATATGAAAATGTTAAATC  >  USA300TCH1516_ALE/1331096‑1331371

aaGGTTTACCAACCGTTAAAGTGTTACAACTGATTACAAAAATTTATGGTGCTGAGAAACGAGTTGAATGGGTTAACCAGTTAGTATCAAATAGTTTTAATAAAAAAATACGTCCATTACGTACACAAAGTGATTGGAttt                                                                                                                                         >  4:269317/1‑141 (MQ=255)

 cggTTTACCAACCGTTAAAGTGTTACAACTGATTACAAAAATTTATGGTGCTGAGAAACGAGATGAATGGGTTAACCAGTTAGTATCAAATAGTTTTAATAAAAAAATACGTCCATTACGTACACAAAGTGATTGGAATTc                                                                                                                                        >  4:382729/2‑141 (MQ=255)

 aGGTTTACCAACCGTTAAAGTGTTACAACTGATTACAAAAATTTATGGTGCTGAGAAACGAGTTGAATGGGTTAACCAGTTAGTATCAAATAGTTTTAATAAAAAAATACGTCCATTACGTACACAAAGTGATTGGATTTc                                                                                                                                        >  3:427142/1‑141 (MQ=255)

  ggTTTACCAACCGTTAAAGCGTTACGACTTAGAAGGAAAATTTATGGTGCTGAGAAACCAGTTGAATGGGTTAACCAGTTAGTATCAAATAGTTTTAATAAAAAAAAACGTCCATTACGTACACAAAGTGATTGGATTTCt                                                                                                                                       >  4:363530/1‑141 (MQ=255)

    tttACCAACCGTTAAAGTGTTACAACTGATTACAAAAATTTATGGTGCTGAGAAACGAGTTGAATGGGTTAACCAGTTAGTATCAAATAGTTTTAATAAAAAAATACGTCCATTACGTACACAAAGTGATTGGATTTCTAg                                                                                                                                     <  2:709200/141‑1 (MQ=255)

    tttACCAACCGTTAAAGTGTTACAACTGATTACAAAAATTTATGGTGCTGAGAAACGAGTTGAATGGGTTAACCAGTTAGTATCAAATAGTTTTAATAAAAAAATACGTCCATTACGTACACAAAGTGATTGGATTTCTAg                                                                                                                                     <  3:398690/141‑1 (MQ=255)

     ttACCAACCGTTAAAGTGTTACAACTGATTACAAAAATTTATGGTGCTGAGAAACGAGTTGAATGGGTTAAGCAGTTAGTATCAAATAGTTTTAATAAAAAAATACGTCCATTACGTACACAAAGTGATTGGATTTCTAGt                                                                                                                                    <  2:642501/141‑1 (MQ=255)

     ttACCAACCGTTAAAGTGTTACAACTGATTACAAAAATTTATGGTGCTGAGAAACGAGTTGAATGGGTTAACCAGTTAGTATCAAATAGTTTTAATAAAAAAATACGTCCATTACGTACACAAAGTGATTGGATTACTAga                                                                                                                                    <  1:285173/141‑2 (MQ=255)

     ttACCAAACGTTAAAGTGTTACAACTGATTACAAAAATTTATGGTGCTGAGAAACGAGTTGAATGGGTTAACCAGTTAGTATCAAATAGTTTTAATAAAAAAATACGTCCATTACGTACACAAAGTGATTGGATTTCTAGt                                                                                                                                    <  3:186994/141‑1 (MQ=255)

      tACCAAGCGTTAAAGTGTTACAACTGATTACAAAAATTTATGGTGCTGAGAAACGAGTTGAATGGGTTAACCAGTTAGTATCAAATAGTTTTAATAAAAAAATACGTCCATCACGTACACAAAGTGATTGGATTACTAgaa                                                                                                                                   <  2:829310/141‑3 (MQ=255)

       aCCAACCGTTAAAGTGTTACAACTGATTACAAAAATTTATGGTGCTGAGAAACGAGTTGAATGGGTTAACCAGTTAGTATCAAATAGTTTTAATAAAAAAATACGTCCATTACGTACACAAAGTGATTGGATTTCTAGTaa                                                                                                                                  >  1:339497/1‑141 (MQ=255)

        ccAACCGTTAAAGTGTTACAACTGATTACAAAAATTTATGGTGCTGAGAAACGAGTTGAATGGGTTAACCAGTTAGTATCAAATAGTTTTAATAAAAAAATACGTCCATTACGTACACAAAGTGATTGGATTTCTAGTAAt                                                                                                                                 >  2:656328/1‑141 (MQ=255)

        ccAACCGTTAAAGTGTTACAACTGATTACAAAAATTTATGGTGCTAAGAAACGAGTTGAATGGGTTAACCAGTTAGTATCAAATAGTTTTAATAAAAAAATACGTCCATTACGTACACAAAGTGATTGGATTTCTAGTAAt                                                                                                                                 <  2:622679/141‑1 (MQ=255)

        cacacagTTAAAGTCAAACAACTACTTAAAAAAATTTATGGTGCTGAGAAAATAGTTTAATGGGTTAATCTGTTAGTATCAAATAGTATAAAAAAAAAAATACGTCCATTACGTACACAAAGTGATTGGATTTCTAGTAAt                                                                                                                                 <  4:200718/135‑1 (MQ=255)

          aaCCGTTAAAGTGTTACAACTGATTACAAAAATTTATGGTGCTGAGAAACGAGTTGAATGGGTTAACCAGTTAGTATCAAATAGTTTTAATAAAAAAATACGTCCATTACGTACACAAAGTGATTGGATTTCTAGTAATcc                                                                                                                               <  3:297836/141‑1 (MQ=255)

          aaCCGTTAAAGTGTTACAACTGATTACAAAAATTTATGGTGCTGAGAAACGAGTTGAATGGGTTAACCAGTTAGTATCAAATAGTTTTAATAAAAAAATACGTCCATTACGTACACAAAGTGATTGGATTACTAGAAATcc                                                                                                                               <  1:708866/141‑1 (MQ=255)

                    gtgtTACAACTGATTACAAAAATTTATGGTGCTGAGAAACGAGTTGAATGGGTTAACCAGTTAGTATCAAATAGTTTTAATAAAAAAATACGTCCATTACGTACACAAAGTGATTGGATTTCTAGTAATCCAATTGAAGTa                                                                                                                     >  2:870405/1‑141 (MQ=255)

                    gtgtTACAACTGATTACAAAAATTTATGGTGCTGAGAAACGAGTTGAATGGGTTAACCAGTTAGTATCAAATAGTTTTAATAAAAAAATACGTCCATTACGTACACAAAGTGATTGGATTTCTAGTAATCCAATTGAAGTa                                                                                                                     >  1:769520/1‑141 (MQ=255)

                     tgtTACAACTGATTACAAAAATTTATGGTGCTGAGAAACGAGTTGAATGGGTTAACCAGTTAGTATCAAATAGTTTTAATAAAAAAATACGTCCATTACGTACACAAAGTGATTGGATTTCTAGTAATCCAATTGAAGTAg                                                                                                                    <  1:996235/141‑1 (MQ=255)

                          cAACTGATTACAAAAATTTATGGTGCTGAGAAACGAGTTGAATGGGTTAACCAGTTAGTATCAAATAGTTTTAATAAAAAAATACGTCCATTACGTACACAAAGTGATTGGATTTCTAGTAATCCAATTGAAGTAGATAAc                                                                                                               <  4:281646/141‑1 (MQ=255)

                          cAACTGATTACAAAAATTTATGGTGCTGAGAAACGAGTTGAATGGGTTAACCAGTTAGTATCAAATAGTTTTAATAAAAAAATACGTCCATTACGTACACAAAGTGATTGGATTTCTAGTAATCCAATTGAAGTAGATAAc                                                                                                               >  1:637331/1‑141 (MQ=255)

                             cTGATTACAAAAATTTATGGTGCTGAGAAACGAGTTGAATGGGTTAACCAGTTAGTATCAAATAGTTTTAATAAAAAAATACGTCCATTACGTACACAAAGTGATTGGATTTCTAGTAATCCAATTGAAGTAGATAACttt                                                                                                            >  2:211012/1‑141 (MQ=255)

                                aTTACAAAAATTTATGGTGCTGAGAAACGAGTTGAATGGGTTAACCAGTTAGTATCAAATAGTTTTAATAAAAAAATACGTCCATTACGTACACAAAGTGATTGGATTTCTAGTAATCCAATTGAAGTAGATAACTttatt                                                                                                         >  2:655365/1‑141 (MQ=255)

                                 ttACAAAAACTTATGGTGCTGAGAAACGAGTTGAATGGGTTAACCAGTTAGTATCAAATAGTTTTAATAAAAAAATACGTCCATTACGTACACAAAGTGATTGGATTTCTAGTAATCCAATTGAAGTAGATAACTttatta                                                                                                        >  1:954724/1‑141 (MQ=255)

                                  tACAAAAATTTATGGTGCTGAGTAACGAGTTGAATGGGTTAACCAGTTAGTATCAAATAGTTTTAATAAAAAAATACGTCCATTACGTACACAAAGTGATTGGATTTCTAGTAATCCAATTGAAGTAGATAACTTTATTaa                                                                                                       >  1:910438/1‑141 (MQ=255)

                                     aaaaaTTTATGGTGCTGAGAAACGAGTTGAATGGGTTAACCAGTTAGTATCAAATAGTTTTAATAAAAAAATACGTCCATTACGTACACAAAGTGATTGGATTTCTAGTAATCCAATTGAAGTAGATAACTTTATTAAAGa                                                                                                    >  2:1086147/1‑141 (MQ=255)

                                     aaaaaTTTATGGTGCTGAGAAACGAGTTGAATGGGTTAACCAGTTAGTATCAAATAGTTTTAATAAAAAAATACGTCCATTACGTACACAAAGTGATTGGATTTCTAGTAATCCAATTGAAGTAGATAACTTTATTAAAGa                                                                                                    >  2:256575/1‑141 (MQ=255)

                                     aaaaaTTTATGGTGCTGAGAAACGAGTTGAATGGGTTAACCAGTTAGTATCAAATAGTTTTAATAAAAAAATACGTCCATTACGTACACAAAGTGATTGGATTTCTAGTAATCCAATTGAAGTAGATAACTTTATTAAAGa                                                                                                    >  1:880828/1‑141 (MQ=255)

                                      aaaaTTTATGGTGCTGAGAAACGAGTTGAATGGGTTAACCAGTTAGTATCAAATAGTTTTAATAAAAAAATACGTCCATTACGTACACAAAGTGATTGGATTTCTAGTAATCCAATTGAAGTAGATAACTTTATTAAAGAt                                                                                                   >  1:10953/1‑141 (MQ=255)

                                      aaaaTTTATGGTGCTGAGAAACGAGTTGAATGGGTTAACCAGTTAGTATCAAATAGTTTTAATAAAAAAATACGTCCATTACGTACACAAAGTGATTGGATTTCTAGTAATCCAATTGAAGTAGATAACTTTATTAAAGAt                                                                                                   >  1:11145/1‑141 (MQ=255)

                                       aaaTTTATGGTGCTGAGAAACGAGTTGAATGGGTTAACCAGTTAGTATCAAATAGTTTTAATAAAAAAATACGTCCATTACGTACACAAAGTGATTGGATTTCTAGTAATCCAATTGAAGTAGATAACTTTATTAAAGATc                                                                                                  >  2:596902/1‑141 (MQ=255)

                                       aaaTTTATGGTGCTGAGAAACGAGTTGAATGGGTTAACCAGTTAGTATCAAATAGTTTTAATAAAAAAATACGTCCATTACGTACACAAAGTGATTGGATTTCTAGTAATCCAATTGAAGTAGATAACTTTATTAAAGATc                                                                                                  <  2:1002947/141‑1 (MQ=255)

                                         aTTTATGGTGCTGAGAAACGAGTTGAATGGGTTAACCAGTTAGTATCAAATAGTTTTAATAAAAAAATACGTCCATTACGTACACAAAGTGATTGGATTTCTAGTAATCCAATTGAAGTAGATAACTTTATTAAAGATCCa                                                                                                <  1:883161/141‑1 (MQ=255)

                                         aTTTATGGTGCTGAGAAACGAGTTGAATGGGTTAACCAGTTAGTATCAAATAGTTTTAATAAAAAAATACGTCCATTACGTACACAAAGTGATTGGATTTCTAGTAATCCAATTGAAGTAGATAACTTTATTAAAGATCCa                                                                                                <  2:769520/141‑1 (MQ=255)

                                         aTTTATGGTGCTGAGAAACGAGTTGAATGGGTTAACCAGTTAGTATCAAATAGTTTTAATAAAAAAATACGTCCATTACGTACACAAAGTGATTGGATTTCTAGTAATCCAATTGAAGTAGATAACTTTATTAAAGATCCa                                                                                                <  1:315375/141‑1 (MQ=255)

                                              tGGTGCTGAGAAACGAGTTGAATGGGTTAACCAGTTAGTATCAAATAGTTTTAATAAAAAAATACGTCCATTACGTACACAAAGTGATTGGATTTCTAGTAATCCAATTGAAGTAGATAACTTTATTAAAGATCCATATAg                                                                                           <  2:984849/141‑1 (MQ=255)

                                                gTGCTGAGAAACGAGTTGAATGGGTTAACCAGTTAGTATCAAATAGTTTTAATAAAAAAATACGTCCATTACGTACACAAAGTGATTGGATTTCTAGTAATCCAATTGAAGTAGATAACTTTATTAAAGATCCATATAGTg                                                                                         >  2:45050/1‑141 (MQ=255)

                                                 tGCTGAGAAACGAGTTGAATGGGTTAACCAGTTAGTATCAAATAGTTTTAATAAAAAAATACGTCCATTACGTACACAAAGTGATTGGATTTCTAGTAATCCAATTGAAGTAGATAACTTTATTAAAGATCCATATAGTgg                                                                                        <  1:258966/141‑1 (MQ=255)

                                                 tGCTGAGAAACGAGTTGAATGGGTTAACCAGTTAGTATCAAATAGTTTTAATAAAAAAATACGTCCATTACGTACACAAAGTGATTGGATTTCTAGTAATCCAATTGAAGTAGATAACTTTATTAAAGATCCATATAGTgg                                                                                        <  1:993510/141‑1 (MQ=255)

                                                 tGCTGAGAAACGAGTTGAATGGGTTAACCAGTTAGTATCAAATAGTTTTAATAAAAAAATACGTCCATTACGTACACAAAGTGATTGGATTTCTAGTAATCCAATTGAAGTAGATAACTTTATTAAAGATCCATATAGTgg                                                                                        <  1:432860/141‑1 (MQ=255)

                                                 tGCTGAGAAACGAGTTGAATGGGTTAACCAGTTAGTATCAAATAGTTTTAATAAAAAAATACGTCCATTACGTACACAAAGTGATTGGATTTCTAGTAATCCAATTGAAGTAGATAACTTTATTAAAGATCCATATAGTgg                                                                                        <  2:618208/141‑1 (MQ=255)

                                                     gagaAACGAGTTGAATGGGTTAACCAGTTAGTATCAAATAGTTTTAATAAAAAAATACGTCCATTACGTACACAAAGTGATTGGATTTCTAGTAATCCAATTGAAGTTGATAACTTTATTAAAGATCCATATAGTGGAttt                                                                                    <  1:563431/141‑1 (MQ=255)

                                                     gagaAACGAGTTGAATGGGTTAACCAGTTAGTATCAAATAGTTTTAATAAAAAAATACGTCCATTACGTACACAAAGTGATTGGATTTCTAGTAATCCAATTGAAGTAGATAACTTTATTAAAGATCCATATAGTGGAttt                                                                                    <  2:527385/141‑1 (MQ=255)

                                                        aaaCGAGTTGAATGGGTTAACCAGTTAGTATCAAATAGTTTTAATAAAAAAATACGTCCATTACGTACACAAAGTGATTGGATTTCTAGTAATCCAATTGAAGTAGATAACTTTATTAAAGATCCATATAGTGGATTTAAt                                                                                 <  1:697915/141‑1 (MQ=255)

                                                        aaaCGAGTTGAATGGGTTAACCAGTTAGTATCAAATAGTTTTAATAAAAAAATACGTCCATTACGTACACAAAGTGATTGGATTTCTAGTAATCCAATTGAAGTAGATAACTTTATTAAAGATCCATATAGTGGATTTAAt                                                                                 <  4:263477/141‑1 (MQ=255)

                                                           cGAGTTGAATGGGTTAACCAGTTAGTATCAAATAGTTTTAATAAAAAAATACGTCCATTACGTACACAAAGTGATTGGATTTCTAGTAATCCAATTGAAGTAGATAACTTTATTAAAGATCCATATAGTGGATTTAAtgtg                                                                              <  3:442858/141‑1 (MQ=255)

                                                           cGAGTTGAATGGGTTAACCAGTTAGTATCAAATAGTTTTAATAAAAAAATACGTCCATTACGTACACAAAGTGATTGGATTTCTAGTAATCCAATTGAAGTAGATAACTTTATTAAAGATCCATATAGTGGATTTAAtgtg                                                                              <  2:940569/141‑1 (MQ=255)

                                                            gAGTTGAATGGGTTAACCAGTTAGTATCAAATAGTTTTAATAAAAAAATACGTCCATTACGTACACAAAGTGATTGGATTTCTAGTAATCCAATTGAAGTAGATAACTTTATTAAAGATCCATATAGTGGATTTAATgtgt                                                                             <  1:870405/141‑1 (MQ=255)

                                                            aaGTTGAATGGGTTAACCAGTTAGTATCAAATAGTTTTAATAAAAAAATACGTCCATTACGTACACAAAGTGATTGGATTTCTAGTAATCCAATTGAAGTAGATAACTTTATTAAAGATCCATATGGTGGATTTAATgtgt                                                                             >  3:348161/2‑141 (MQ=255)

                                                             aGTTGAATGGGTTAACCAGTTAGTATCAAATAGTTTTAATAAAAAAATACGTCCATTACGTACACAAAGTGATTGGATTTCTAGTAATCCAATTGAAGTAGATAACTTTATTAAAGATCCATATAGTGGATTTAATGTGTc                                                                            >  4:114326/1‑141 (MQ=255)

                                                              gTTGAATGGGTTAACCAGTTAGTATCAAATAGTTTTAATAAAAAAATACGTCCATTACGTACACAAAGTGATTGGATTTCTAGTAATCCAATTGAAGTAGATAACTTTATTAAAGATCCATATAGTGGATTTAATGTGTCa                                                                           <  3:262180/141‑1 (MQ=255)

                                                               ttGAATGGGTTAACCAGTTAGTATCAAATAGTTTTAATAAAAAAATACGTCCATTACGTACACAAAGTGATTGGATTTCTAGTAATCCAATTGAAGTAGATAACTTTATTAAAGATCCATATAGTGGATTTAATGTGTCaa                                                                          <  1:589710/141‑1 (MQ=255)

                                                               ttGAATGGGTTAACCAGTTAGTATCAAATAGTTTTAATAAAAAAATACGTCCATTACGTACACAAAGTGATTGGATTTCTAGTAATCCAATTGAAGTAGATAACTTTATTAAAGATCCATATAGTGGATTTAATGTGTCaa                                                                          >  2:759569/1‑141 (MQ=255)

                                                                 gAATGGGTTAACCAGTTAGTATCAAATAGTTTTAATAAAAAAATACGTCCATTACGTACACAAAGTGATTGGATTTCTAGTAATCCAATTGAAGTAGATAACTTTATTAAAGATCCATATAGTGGATTTAATGTGTCAAAt                                                                        >  2:968392/1‑141 (MQ=255)

                                                                 gAATGGGTTAACCAGTTAGTATCAAATAGTTTTAATAAAAAAATACGTCCATTACGTACACAAAGTGATTGGATTTCTAGTAATCCAATTGAAGTAGATAACTTTATTAAAGATCCATATAGTGGATTTAATGTGTCAAAt                                                                        >  1:357393/1‑141 (MQ=255)

                                                                 aaaTGGGTTAACCAGTTAGTATCAAATAAAAAAAAAAAAAAAATACGTCCATTACGTACACAAAGTGATTGGATTTCTAGTAATCCAATTGAAGTAGATAACTTTATTAAAGATCCATATAGTGGATTTAATGTGTCAAAt                                                                        <  2:87523/140‑1 (MQ=255)

                                                                  aaTGGGTTAACCAGTTAGTATCAAATAGTTTTAATAAAAAAATACGTCCGTTACGTACACAAAGTGATTGGATTTCTAGTAATCCAATTGAAGTAGATAACTTTATTAAAGATCCATATAGTGGATTTAATGTGTCAAATc                                                                       <  3:163390/141‑1 (MQ=255)

                                                                  aaTGGGTTAACCAGTTAGTATCAAATAGTTTTAATAAAAAAATACGTCCATTACGTACACAAAGTGATTGGATTTCTAGTAATCCAATTGAAGTAGATAACTTTATTAAAGATCCATATAGTGGATTTAATGTGTCAAATc                                                                       <  1:635574/141‑1 (MQ=255)

                                                                     gggTTAACCAGTTAGTATCAAATAGTTTTAATAAAAAAATACGTCCATTACGTACACAAAGTGATTGGATTTCTAGTAATCCAATTGAAGTAGATAACTTTATTAAAGATCCATATAGTGGATTTAATGTGTCAAATCAat                                                                    <  2:700451/141‑1 (MQ=255)

                                                                      ggTTAACCAGTTAGTATCAAATAGTTTTAATAAAAAAATACGTCCATTACGTACACAAAGTGATTGGATTTCTAGTAATCCAATTGAAGTAGATAACTTTATTAAAGATCCATATAGTGGATTTAATGTGTCAAATCAatt                                                                   >  2:8178/1‑141 (MQ=255)

                                                                          aaCCAGTTAGTATCAAATAGTTTTAATAAAAAAATACGTCCATTACGTACACAAAGTGATTGGATTTCTAGTAATCCAATTGAAGTAGATAACTTTATTAAAGATCCATATAGTGGATTTAATGTGTCAAATCAAttatta                                                               >  2:112718/1‑141 (MQ=255)

                                                                          aaCCAGTTAGTATCAAATAGTTTTAATAAAAAAATACGTCCATTACGTACACAAAGTGATTGGATTTCTAGTAATCCAATTGAAGTAGATAACTTTATTAAAGATCCATATAGTGGATTTAATGTGTCAAATCAAttatta                                                               >  3:217337/1‑141 (MQ=255)

                                                                           aCCAGTTAGTATCAAATAGTTTTAATAAAAAAATACGTCCATTACGTACACAAAGTGATTGGATTTCTAGTAATCCAATTGAAGTAGATAACTTTATTAAAGATCCATATAGTGGATTTAATGTGTCAAATCAATtattat                                                              >  2:945216/1‑141 (MQ=255)

                                                                            ccAGTTAGTATCAAATAGTTTTAATAAAAAAATACGTCCATTACGTACACAAAGTGATTGGATTTCTAGTAATCCAATTGAAGTAGATAACTTTATTAAAGATCCATATAGTGGATTTAATGTGTCAAATCAATTATtata                                                             >  1:781207/1‑141 (MQ=255)

                                                                               gTTAGTATCAAATAGTTTTAATAAAAAAATACGTCCATTACGTACACAAAGTGATTGGATTTCTAGTAATCCAATTGAAGTAGATAACTTTATTAAAGATCCATATAGTGGATTTAATGTGTCAAATCAATTATTATATCa                                                          >  1:670854/1‑141 (MQ=255)

                                                                               gTTAGTATCAAATAGTTTTAATAAAAAAATACGTCCATTACGTACACAAAGTGATTGGATTTCTAGTAATCCAATTGAAGTAGATAACTTTATTAAAGATCCATATAGTGGATTTAATGTGTCAAATCAATTATTATATCa                                                          <  2:567378/141‑1 (MQ=255)

                                                                                ttAGTATCAAATAGTTTTAATAAAAAAATACGTCCATTACGTACACAAAGTGATTGGATTTCTAGTAATCCAATTGAAGTAGATAACTTTATTAAAGATCCATATAGTGGATTTAATGTGTCAAATCAATTATTATATCaa                                                         >  4:341714/1‑141 (MQ=255)

                                                                                ttAGTATCAAATAGTTTTAATAAAAAAATACGTCCATTACGTACACAAAGTGATTGGATTTCTAGTAATCCAATTGAAGTAGATAACTTTATTAAAGATCCATATAGTGGATTTAATGTGTCAAATCAATTATTATATCaa                                                         >  3:424513/1‑141 (MQ=255)

                                                                                  aGTATCAAATAGTTTTAATAAAAAAATACGTCCATTACGTACACAAAGTGATTGGATTTCTAGTAATCCAATTGAAGTAGATAACTTTATTAAAGATCCATATAGTGGATTTAATGTGTCAAATCAATTATTATATCAAAc                                                       <  2:83070/141‑1 (MQ=255)

                                                                                      tCAAATAGTTTTAATAAAAAAATACGTCCATTACGTACACAAAGTGATTGGATTTCTAGTAATCCAATTGAAGTAGATAACTTTATTAAAGATCCATATAGTGGATTTAATGTGTCAAATCAATTATTATATCAAACAGcc                                                   >  1:720609/1‑141 (MQ=255)

                                                                                        aaaTAGTTTTAATAAAAAAATACGTCCATTACGTACACAAAGTGATTGGATTTCTAGTAATCCAATTGAAGTAGATAACTTTATTAAAGATCCATATAGTGGATTTAATGTGTCAAATCAATTATTATATCAAACAGCCta                                                 >  2:751198/1‑141 (MQ=255)

                                                                                        aaaTAGTTTTAATAAAAAAATACGTCCATTACGTACACAAAGTGATTGGATTTCTAGTAATCCAATTGAAGTAGATAACTTTATTAAAGATCCATATAGTGGATTTAATGTGTCAAATCAATTATTATATCAAACAGCCta                                                 >  1:38963/1‑141 (MQ=255)

                                                                                           tAGTTTTAATAAAAAAATACGTCCATTACGTACACAAAGTGATTGGATTTCTAGTAATCCAATTGAAGTAGATAACTTTATTAAAGATCCATATAGTGGATTTAATGTGTCAAATCAATTATTATATCAAACAGCCtatta                                              >  1:636074/1‑141 (MQ=255)

                                                                                             gTTTTAATAAAAAAATACGTCCATTACGTACACAAAGTGATTGGATTTCTAGTAATCCAATTGAAGTAGATAACTTTATTAAAGATCCATATAGTGGATTTAATGTGTCAAATCAATTATTATATCAAACAGCCTATtata                                            >  1:541503/1‑141 (MQ=255)

                                                                                                  aataaAAAAATACGTCCATTACGTACACAAAGTGATTGGATTTCTAGTAATCCAATTGAAGTAGATAACTTTATTAAAGATCCATATAGTGGATTTAATGTGTCAAATCAATTATTATATCAAACAGCCTATTATATGCta                                       >  3:80784/1‑141 (MQ=255)

                                                                                                      aaaaaaTACGTCCATTACGTACACAAAGTGATTGGATTTCTAGTAATCCAATTGAAGTAGATAACTTTATTAAAGATCCATATAGTGGATTTAATGTGTCAAATCAATTATTATATCAAACAGCCTATTATATGCtacata                                   >  2:303911/1‑141 (MQ=255)

                                                                                                       aaaaaTACGTCCATTACGTACACAAAGTGATTGGATTTCTAGTAATCCAATTGAAGTAGATAACTTTATTAAAGATCCATATAGTGGATTTAATGTGTCAAATCAATTATTATATCAAACAGCCTATTATATGCtacatac                                  <  4:184077/141‑1 (MQ=255)

                                                                                                       aaaaaTACGTCCATTACGTACACAAAGTGATTGGATTTCTAGTAATCCAATTGAAGTAGATAACTTTATTAAAGATCCATATAGTGGATTTAATGTGTCAAATCAATTATTATATCAAACAGCCTATTATATGCtacatac                                  <  2:696314/141‑1 (MQ=255)

                                                                                                        aaaaTACGTCCATTACGTACACAAAGTGATTGGATTTCTAGTAATCCAATTGAAGTAGATAACTTTATTAAAGATCCATATAGTGGATTTAATGTGTCAAATCAATTATTATATCAAACAGCCTATTATATGCtacataca                                 >  1:220709/1‑141 (MQ=255)

                                                                                                        aaaaTACGTCCATTACGTACACAAAGTGATTGGATTTCTAGTAATCCAATTGAAGTAGATAACTTTATTAAAGATCCATATAGTGGATTTAATGTGTCAAATCAATTATTATATCAAACAGCCTATTATATGCtacataca                                 <  1:146120/141‑1 (MQ=255)

                                                                                                          aaTACGTCCATTACGTACACAAAGTGATTGGATTTCTAGTAATCCAATTGAAGTAGATAACTTTATTAAAGATCCATATAGTGGATTTAATGTGTCAAATCAATTATTATATCAAACAGCCTATTATATGCTACATACATc                               <  1:996988/141‑1 (MQ=255)

                                                                                                            tACGTCCATTACGTACACAAAGTGATTGGATTTCTAGTAATCCAATTGAAGTAGATAACTTTATTAAAGATCCATATAGTGGATTTAATGTGTCAAATCAATTATTATATCAAACAGCCTATTATATGCTACATACATcac                             <  3:286240/141‑1 (MQ=255)

                                                                                                            tACGTCCATTACGTACACAAAGTGATTGGATTTCTAGTAATCCAATTGAAGTAGATAACTTTATTAAAGATCCATATAGTGGATTTAATGTGTCAAATCAATTATTATATCAAACAGCCTATTATATGCTACATACATcac                             <  2:786703/141‑1 (MQ=255)

                                                                                                            tACGTCCATTACGTACACAAAGTGATTGGATTTCTAGTAATCCAATTGAAGTAGATAAAAAAATTAAAAAAACATATAGTGGATTTAATGTGTCAAATCAATTATTATATCAAACAGCCTATTATATGCtacatacaaaaa                             <  2:725765/141‑5 (MQ=255)

                                                                                                              cGTCCATTACGTACACAAAGTGATTGGATTTCTAGTAATCCAATTGAAGTAGATAACTTTATTAAAGATCCATATAGTGGATTTAATGTGTCAAATCAATTATTATATCAAACAGCCTATTATATGCTACATACATCACaa                           >  2:242331/1‑141 (MQ=255)

                                                                                                              cGTCCATTACGTACACAAAGTGATTGGATTTCTAGTAATCCAATTGAAGTAGATAACTTTATTAAAGATCCATATAGTGGATTTAATGTGTCAAATCAATTATTATATCAAACAGCCTATTATATGCTACATACATCACaa                           >  3:117476/1‑141 (MQ=255)

                                                                                                               gTCCATTACGTACACAAAGTGATTGGATTTCTAGTAATCCAATTGAAGTAGATAACTTTATTAAAGATCCATATAGTGGATTTAATGTGTCAAATCAATTATTATATCAAACAGCCTATTATATGCTACATACATCACAAt                          <  1:231973/141‑1 (MQ=255)

                                                                                                                  cATTACGTACACAAAGTGATTGGATTTCTAGTAATCCAATTGAAGTAGATAACTTTATTAAAGATCCATATAGTGGATTTAATGTGTCAAATCAATTATTATATCAAACAGCCTATTATATGCTACATACATCACAATTaa                       >  2:966295/1‑141 (MQ=255)

                                                                                                                      aCGTACACAAAGTGATTGGATTTCTAGTAATCCAATTGAAGTAGATAACTTTATTAAAGATCCATATAGTGGATTTAATGTGTCAAATCAATTATTATATCAAACAGCCTATTATATGCTACATACATCACAATTaaaaaa                   <  1:946510/141‑1 (MQ=255)

                                                                                                                        gTACCCAAAGTGATTGGATTTCTAGTAATCCAATTGAAGTAGATAACTTTATTAAAGATCCATATAGTGGATTTAATGTGTCAAATCAATTATTATATCAAACAGCCTATTATATGCTACATACATCACAATTaaaaaaca                 >  4:159357/1‑139 (MQ=255)

                                                                                                                         tACACAAAGTGATTGGATTTCTAGTAATCCAATTGAAGTAGATAACTTTATTAAAGATCCATATAGTGGATTTAATGTGTCAAATCAATTATTATATCAAACAGCCTATTATATGCTACATACATCACAATTAAAAAatat                >  1:12492/1‑141 (MQ=255)

                                                                                                                         tACACAAAGTGATTGGATTCCTAGTAATCCAATTGAAGTAGATAACTTTATTAAAGATCCATATAGTGGATTTAATGTGTCAAATCAATTATTATATCAAACAGCCTATTATATGCTACATACATCACAATTAAAAAatat                >  2:495797/1‑141 (MQ=255)

                                                                                                                          acatcaaGTGAACGGATCTGTAGTAATCCAATTGATGTGGTTAACTTTATTAAAGATCCATATAGTGGATTTAATGTGTAAACTAAAATATTATCTCAAACAGCCCATTATATCCTTCATACATCACTATTACAAAatacg               >  4:187850/6‑139 (MQ=255)

                                                                                                                           cacaAAGTGATTGGATTTCTAGTAATCCAATTGAAGTAGATAACTTTATTAAAGATCCATATAGTGGATTTAATGTGTCAAATCAATTATTATATCAAACAGCCTATTATATGCTACATACATCACAATTAAAAAATATGa              >  1:81776/1‑141 (MQ=255)

                                                                                                                            acaAAGTGATTGGATTTCTAGTAATCCAATTGAAGTAGATAACTTTATTAAAGATCCATATAGTGGATTTAATGTGTCAAATCAATTATTATATCAAACAGCCTATTATATGCTACATACATCACAATTAAAAAATATGaa             >  4:191053/1‑141 (MQ=255)

                                                                                                                            acaAAGTGATTGGATTTCTAGTAATCCAATTGAAGTAGATAACTTTATTAAAGATCCATATAGTGGATTTAATGTGTCAAATCAATTATTATATCAAACAGCCTATTATATGCTACATACATCACAATTAAAAAATATGaa             <  2:670854/141‑1 (MQ=255)

                                                                                                                               aaGTGATTGGATTTCTAGTAATCCAATTGAAGTAGATAACTTTATTAAAGATCCATATAGTGGATTTAATGTGTCAAATCAATTATTATATCAAACAGCCTATTATATGCCACATACATCACAATTAAAAAATATGAAAAt          >  1:632071/1‑141 (MQ=255)

                                                                                                                                     ttGGATTTCTAGTAATCCAATTGAAGTAGATAACTTTATTAAAGATCCATATAGTGGATTTAATGTGTCAAATCAATTATTATATCAAACAGCCTATTATATGCTACATACATCACAATTAAAAAATATGAAAATGTTaaa    >  1:104483/1‑141 (MQ=255)

                                                                                                                                       ggATTTCTAGTAATCCAATTGAAGTAGATAACTTTATTAAAGATCCATATAGTGGATTTAATGTGTCAAATCAATTATTATATCAAACAGCCTATTATATGCTACATACATCACAATTAAAAAATATGAAAATGTTAAATc  <  2:435486/141‑1 (MQ=255)

                                                                                                                                       ggATTACTAGAAATCCAATTGAAGTAGATAACTTTATTAAAGATCCATATAGTGGATTTAATGTGTCAAATCAATTATTATATCAAACAGCCTATTATATGCt                                        >  2:976313/1‑103 (MQ=255)

AAGGTTTACCAACCGTTAAAGTGTTACAACTGATTACAAAAATTTATGGTGCTGAGAAACGAGTTGAATGGGTTAACCAGTTAGTATCAAATAGTTTTAATAAAAAAATACGTCCATTACGTACACAAAGTGATTGGATTTCTAGTAATCCAATTGAAGTAGATAACTTTATTAAAGATCCATATAGTGGATTTAATGTGTCAAATCAATTATTATATCAAACAGCCTATTATATGCTACATACATCACAATTAAAAAATATGAAAATGTTAAATC  >  USA300TCH1516_ALE/1331096‑1331371

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