BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	USA300TCH1516_ALE	2,721,481	T→C	8.7%	L198S (TTA→TCA)	USA300HOU_RS14000 →	putative hydrolase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	USA300TCH1516_ALE	2,721,481	0	T	C	8.7%	436.4 / 23.4	159	L198S (TTA→TCA)	USA300HOU_RS14000	putative hydrolase
Reads supporting (aligned to +/- strand): ref base T (80/65); new base C (6/8); total (86/73)
Fisher's exact test for biased strand distribution p-value = 4.12e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 160 total aligned reads displayed.

TAACAGAAGGCTTAGAAAAAGGCATGAAAACATTTGGCGAAACATTAAATATCGCACCAATTGATGCAAAAATGATGTCTCAACCTGCTGATACTGAAGAAGGACGTATAGAACAGTACAAACGTACGATGTTCTGGTTAGAATTTGAAATTAGACAATATACGCATTCAAATATTACTTTGGATGATTTCACAAAATACAGCGATAAAATCACCCTACTAAATGGGACAGATTCTAGAGGCTCCTTCCCACAAGATGTTAACTTTTATATAAATAAAG  >  USA300TCH1516_ALE/2721343‑2721621

tAACAGAAGGCTTAGAAAAAGGCATGAAAACATTTGGCGAAACATTAAATATCGCACCAATTGATGCAAAAATGATGTCTCAACCTGCTGATACTGAAGAAGGACGTATAGAACAGTACAAACGTACGATGTTCTGGTTAg                                                                                                                                            >  3:18512/1‑141 (MQ=255)

tAACAGAAGGCTTAGAAAAAGGCATGAAAACATTTGGCGAAACATTAAATATCGCACCAATAGATGCAAAAATGATGTCTCAACCTGCTGATACTGAAGAAGGACGTATAGAACAGTACAAACGTACGATGTTCTGGTTaa                                                                                                                                            >  4:307304/1‑140 (MQ=255)

 aaCAGAAGGCTTAGAAAAAGGTATGAAAACATTTGGTGAAACATTAAATATCGCACCAATTGATGCAAAAATGATGTCTCAACCTGCTGATACTGAAGAAGGACGTATAGAACAGTACAAACGTACGATGTTCTGGTCAGa                                                                                                                                           >  2:366222/1‑141 (MQ=255)

 aaCAGAAGGCTTAGAAAAAGGCATGAAAACATTTGGCGAAACATTAAATATCGCACCAATTGATGCAAAAATGATGTCTCAACCTGCTGATACTGAAGAAGGACGTATAGAACAGTACAAACGTACGATGTTCTGGTTAGa                                                                                                                                           >  2:926265/1‑141 (MQ=255)

   cAGAAGGCTTAGAAAAAGGCATGAAAGCAGTTGGCGAAACATTAAATATCGCACCAATTGATGCAAAAATGATGTCTCAACCTGCTGATACTGAAGAAGGCCGTATAGAACAGTACAAACGTACGATGTTCTGGTTAGAAt                                                                                                                                         <  1:198787/141‑1 (MQ=255)

   aaGAAGGCTTAGAAAAAGGCATGAAACCATTTGGCGAAACATTAAGTATCGCACCAATTGATGCAAAAATGATGTCTCAACCTGCGGATACTGAAGAAGGACGTATAGAACAGTACAAACGTACGATGTTCTGGTTAGAAt                                                                                                                                         <  3:182207/140‑1 (MQ=255)

    aGAAGGCTTAGAAAAAGGCATGAAAACATTTGGCGAAACATTAAATATCGCACCAATTGATGCAAAAATGATGTCTCAACCTGCTGATACTGAAGAAGGACGTATAGAACAGTACAAACGTACGATGTTCTGGTTAGAAtt                                                                                                                                        >  1:967242/1‑141 (MQ=255)

     gAAGGCTTAGAAAAAGGCATGAAAACATTTGGCGAAACATTAAATATCGCACCAATTGATGCAAAAATGATGTCTCAACCTGCTGATACTGAAGAAGGACGTATAGAACAGTACAAACGTACGATGTTCTGGTTAGAAttt                                                                                                                                       <  2:383862/141‑1 (MQ=255)

      aaGGCTTAGAAAAAGGCATGAAAACATTTGGCGAAACATTAAATATCGCACCAATTGATGCAAAAATGATGTCTCAACCTGCTGATACTGAAGAAGGACGTATAGAACAGTACAAACGTACGATGTTCTGGTTAGAATTTg                                                                                                                                      >  2:1059725/1‑141 (MQ=255)

      aaGGCTTAGAAAAAGGCATGAAAACATTTGGCGAAACATTAAATATCGCACCAATTGATGCAAAAATGATGTCTCAACCTGCTGATACTGAAGAAGGACGTATAGAACAGTACAAACGTACGATGTTCTGGTTAGAATTTg                                                                                                                                      <  2:987721/141‑1 (MQ=255)

           ttAGAAAAAGGTATGAAAACATTTGGTGAAACATTAAATATCGCACCAATTGATGCAAAAATGATGTCTCAACCTGCTGATACTGAAGAAGGACGTATAGAACAGTACAAACGTACGATGTTCTGGTCAGAATTTGAAAtt                                                                                                                                 >  1:386567/1‑141 (MQ=255)

            tAGAAAAAGGCATGAAAACATTTGGCGAAACATTAAATATCGCACCAATTGATGCAAAAATGATGTCTCAACCTGCTGATACTGAAGAAGGACGTATAGAACAGTACAAACGTACGATGTTCTGGTTAGAATTTGAAATTa                                                                                                                                >  1:1072778/1‑141 (MQ=255)

             aGAAAAAGGCATGAAAACATTTGGCGAAACATTAAATATCGCACCAATTGATGCAAAAATGATGTCTCAACCTGCTGATACTGAAGAAGGACGTATAGAACAGTACAAACGTACGATGTTCTGGTTAGAATTTGAAATTAg                                                                                                                               >  2:248056/1‑141 (MQ=255)

                 aaaGGTATGAAAACATTTGGTGAAACATTAAATATCGCACCAATTGATGCAAAAATGATGTCTCAACCTGCTGATACTGAAGAAGGACGTATAGAACAGTACAAACGTACGATGTTCTGGTCAGAATTTGAAATTAGACaa                                                                                                                           <  4:350706/141‑1 (MQ=255)

                  aaGGCATGAAAACATTTGGCGAAACATTAAATATCGCACCAATTGATGCAAAAATGATGTCTCAACCTGCTGATACTGAAGAAGGACGTATAGAACAGTACAAACGTACGATGTTCTGGTTAGAATTTGAAATTAGACAat                                                                                                                          >  2:20731/1‑141 (MQ=255)

                  aaGGCATGAAAACATTTGGCGAAACATTAAATATCGCACCAATTGATGCAAAAATGATGTCTCAACCTGCTGATACTGAAGAAGGACGTATAGAACAGTACAAACGTACGATGTTCTGGTTAGAATTTGAAATTAGACAat                                                                                                                          >  4:47909/1‑141 (MQ=255)

                   aGGCATGAAAACATTTGGCGAAACATTAAATATCGCACCAATTGATGCAAAAATGATGTCTCAACCTGCTGATACTGAAGAAGGACGTATAGAACAGTACAAACGTACGATGTTCTGGTTAGAATTTGAAATTAGACAata                                                                                                                         <  2:528172/141‑1 (MQ=255)

                    ggCATGAAAACATTTGGCGAAACATTAAATATCGCACCAATTGATGCAAAAATGATGTCTCAACCTGCTGATACTGAAGAAGGACGTATAGAACAGTACAAACGTACGATGTTCTGGTTAGAATTTGAAATTAGACAatat                                                                                                                        >  1:547077/1‑141 (MQ=255)

                    ggCATGAAAACATTTGGCGAAACATTAAATATCGCACCAATTGATGCAAAAATGATGTCTCAACCTGCTGATACTGAAGAAGGACGTATAGAACAGTACAAACGTACGATGTTCTGGTTAGAATTTGAAATTAGACAatat                                                                                                                        >  4:131141/1‑141 (MQ=255)

                    ggCATGAAAACATTTGGCGAAACATTAAATATCGCACCAATTGATGCAAAAATGATGTCTCAACCTGCTGATACTGAAGAAGGACGTATAGAACAGTACAAACGTACGATGTTCTGGTTAGAATTTGAAATTAGACAatat                                                                                                                        >  1:1054454/1‑141 (MQ=255)

                         gAAAACATTTGGCGAAACATTAAATATCGCACCAATTGATGCAAAAATGATGTCTCAACCTGCTGATACTGAAGAAGGACGTATAGAACAGTACAAACGTACGATGTTCTGGTTAGAATTTGAAATTAGACAATATACGCa                                                                                                                   <  2:237749/141‑1 (MQ=255)

                          aaaaCATTTGGCGAAACATTAAATATCGCACCAATTGATGCAAAAATGATGTCTCAACCTGCTGATACTGAAGAAGGACGTATAGAACAGTACAAACGTACGATGTTCTGGTTAGAATTTGAAATTAGACAATATACGCAt                                                                                                                  <  4:224626/141‑1 (MQ=255)

                          aaaaCATTTAGCGAAACATTAAATAACGCACCAATTGATGCAAAAATGATATCTCAACCTGCTGATACTGAAGAAGGACGTATAGAACAGTACAAACGTACGATGTTCTGGTTAGAATTTGAAATAAGACAACAAACGCag                                                                                                                  >  2:1069823/1‑140 (MQ=255)

                             aCATTTGGCGAAACATTAAATATCGCACCAATTGATGCAAAAATGATGTCTCAACCTGCTGATACTGAAGAAGGACGTATAGAACAGTACAAACGTACGATGTTCTGGTTAGAATTTGAAATTAGACAATATACGCATTCa                                                                                                               >  2:73644/1‑141 (MQ=255)

                             aCATTTGGCGAAACATTAAATATCGCACCAATTGATGCAAAAATGATGTCTCAACCTGCTGATACTGAAGAAGGACGTATAGAACAGTACAAACGTACGATGTTCTGGTTAGAATTTGAAATTAGACAATATACGCATTCa                                                                                                               <  3:42977/141‑1 (MQ=255)

                             aCATTTGGCGAAACATTAAATATCGCACCAATTGATGCAAAAATGATGTCTCAACCTGCTGATACTGAAGAAGGACGTATAGAACAGTACAAACGTACGATGTTCTGGTTAGAATTTGAAATTAGACAATATACGCATTCa                                                                                                               <  2:830974/141‑1 (MQ=255)

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

                                                     gCACCAATTGATGCAAAAATGATGTCTCAACCTGCTGATACTGAAGAAGGACGTATAGAACAGTACAAACGTACGATGTTCTGGTTAGAATTTGAAATTAGACAATATACGCATTCAAATATTACTTTGGATGATTTcaca                                                                                       >  1:59621/1‑141 (MQ=255)

                                                      cACCAACTGATGCAAAAATGATGTCTCAACCTGCTGATACTGAAGAAGGACGTATAGAACAGTACAAACGTACGATGTTCTGGTTAGAATTTGAAATTAGACAATATACGCATTCAAATATTACTTTGGATGATTTCACaa                                                                                      <  3:416564/141‑1 (MQ=255)

                                                         cAATTGATGCAAAAATGATGTCTCAACCTGCTGATACTGAAGAAGGACGTATAGAACAGTACAAACGTACGATGTTCTGGTTAGAATTTGAAATTAGACAATATACGCATTCAAATATTACTTTGGATGATTTCACAAAAt                                                                                   >  2:783018/1‑141 (MQ=255)

                                                            ttGATGCAAAAATGATGTCTCAACCTGCTGATACTGAAGAAGGACGTATAGAACAGTACAAACGTACGATGTTCTGGTTAGAATTTGAAATTAGACAATATACGCATTCAAATATTACTTTGGATGATTTCACAAAATACa                                                                                <  1:873715/141‑1 (MQ=255)

                                                              gATGCAAAAATGATGTCTCAACCTGCTGATACTGAAGAAGGACGTATAGAACAGTACAAACGTACGATGTTCTGGTTAGAATTTGAAATTAGACAATATACGCATTCAAATATTACTTTGGATGATTTCACAAAATACAGc                                                                              >  2:869755/1‑141 (MQ=255)

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

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

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

                                                                     aaaTGATGTCTCAACCTGCTGATACTGAAGAAGGACGTATAGAACAGTACAAACGTACGATGTTCTGGTTAGAATTTGAAATTAGACAATATACGCATTCAAATATTACTTTGGATGATTTCACAAAATACAGCGATaaaa                                                                       <  4:267872/141‑1 (MQ=255)

                                                                        tgatgTCTCAACCTGCTGATACTGAAGAAGGACGTATAGAACAGTACAAACGTACGATGTTCTGGTTAGAATTTGAAATTAGACAATATACGCATTCAAATATTACTTTGGATGATTTCACAAAATACAGCGATAAAATCa                                                                    <  2:1065958/141‑1 (MQ=255)

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

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

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

                                                                                 aaCCTGCTGATACTGAAGAAGGACGTATAGAACAGTACAAACGTACGATGTTCTGGTTAGAATTTGAAATTAGACAATATACGCATTCAAATATTACTTTGGATGATTTCACAAAATACAGCGATAAAATCACCCTACTaa                                                           >  1:75692/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)

                                                                                         gATACTGAAGAAGGACGTATAGAACAGTACAAACGTACGATGTTCTGGTTAGAATTTGAAATTAGACAATATACGCATTCAAATATTACTTTGGATGATTTCACAAAATACAGCGATAAAATCACCCTACTAAATGGGACa                                                   >  4:246833/1‑141 (MQ=255)

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

                                                                                          aTACTGAAAAAGGACGTATAGAACAGTACAAACGTACAATGTTCTGGTTAGAATTTGAAATAAGACAATATACGTACTCAAACATTACTTTGGAGGATTTCACAAAATACAGCGATAAAAACACCCGCCTAAATGGGACAg                                                  >  4:81554/1‑141 (MQ=255)

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

                                                                                                aagaagGACGTATAGAACAGTACAAACGTACGATGTTCTGGTTAGAATTTGAAATTAGACAATATACGCATTCAAATATTACTTTGGATGATTTCACAAAATACAGCGATAAAATCACCCTACTAAATGGGACAGATTCTa                                            >  2:416878/1‑141 (MQ=255)

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

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

                                                                                                   aagGACGTATAGAACAGTACAAACGTACGATGTTCTGGTTAGAATTTGAAATTAGACAATATACGCATTCAAATATTACTTTGGATGATTTCACAAAATACAGCGATAAAATCACCCTACTAAATGGGACAGATTCTagag                                         >  1:289087/1‑141 (MQ=255)

                                                                                                     ggACGTATAGAACAGTACAAACGTACGATGTTCTGGTTAGAATTTGAAATTAGACAATATACGCATTCAAATATTACTTTGGATGATTTCACAAAATACAGCGATAAAATCACCCTACTAAATGGGACAGATTCTAGAGGc                                       >  1:134109/1‑141 (MQ=255)

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

                                                                                                        cGTATAGAACAGTACAAACGTACGATGTTCTGGTTAGAATTTGAAATTAGACAATATACGCATTCAAATATTACTTTGGATGATTTCACAAAATACAGCGATAAAATCACCCTACTAAATGGGACAGATTCTAGAGGCTcc                                    >  2:351317/1‑141 (MQ=255)

                                                                                                        cGTATAGAACAGTACAAACGTACGATGTTCTGGTTAGAATTTGAAATTAGACAATATACGCATTCAAATATTACTTTGGATGATTTCACAAAATACAGCGATAAAATCACCCTACTAAATGGGACAGATTCTAGAGGCTcc                                    <  1:797328/141‑1 (MQ=255)

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

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

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

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

                                                                                                               aaCAGTACAAACGTACGATGTTCTGGTTAGAATTTGAAATTAGACAATATACGCATTCAAATATTACTTTGGATGATTTCACAAAATACAGCGATAAAATCACCCTACTAAATGGGACAGATTCTAGAGGCTCCTTCCcac                             <  1:646780/141‑1 (MQ=255)

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

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

                                                                                                                     aCAAACGTACGATGTTCTGGTCAGAATTTGAAATTAGACAATATACGCATTCAGATATTACTTTGGATGATTTCACAAAATACAGCGATAAAATCACCCTACTAAATGGGACACATTCTAGAGGCTCCTTCCCACAAGATg                       <  4:360532/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)

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

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

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

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

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

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

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

                                                                                                                          cGTACGATGTTCTGGTCAGAATTTGAAATTAGACAATATACGCATTCAGATATTACTTTGGATGATTTCACAAAATACAGCGATAAAATCACCCGACTAAATGGGACAGATTCTAGAGGCTCCTTCCCACAAGATGTTAAc                  <  2:968187/141‑1 (MQ=255)

                                                                                                                             aCGATGTTCTGGTTAGAGTTTGAAATTCGACAATATACGCATTCAAATATTACTTTGGATGAATTCACAAAATACAGCGATAAAATCACCCTACTATATGGGACAGAATCTAGGGGCTCATTCCCACTAGATGTTAACttt               >  2:159364/1‑141 (MQ=255)

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

                                                                                                                                aTGTTCTGGTTAGAATTTGAAATTAGACAATATACGCATTCAAATATTACTTTGGATGATTTCACAAAATACAGCGATAAAATCACCCTACTATATGGGACAGATTCTAGAGGCTCCTTCCCACAAGATGTTAACTTTtat            >  4:83419/1‑141 (MQ=255)

                                                                                                                                aTGTTCTGGTTAGAATTTGAAATTAGACAATATACGCATTCAAATATTACTTTGGATGATTTCAAAAAATACAGCGATAAAATAACCCTACTAAATGGGACAGATTCTAGAGGCTCCGTTCCACCAGATGTTAACTTTtat            >  1:327636/1‑141 (MQ=255)

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

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

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

                                                                                                                                        gTTAGAATTTGAAATTAGACAATATACGCATTCAAATATTACTTTGGATGATTTCACAAAATACAGCGATAAAATCACCCTACTAAATGGGACAGATTCTAGAGGCTCCTTCCCACAAGATGTTAACTTTTATataaataa    <  3:303581/141‑1 (MQ=255)

TAACAGAAGGCTTAGAAAAAGGCATGAAAACATTTGGCGAAACATTAAATATCGCACCAATTGATGCAAAAATGATGTCTCAACCTGCTGATACTGAAGAAGGACGTATAGAACAGTACAAACGTACGATGTTCTGGTTAGAATTTGAAATTAGACAATATACGCATTCAAATATTACTTTGGATGATTTCACAAAATACAGCGATAAAATCACCCTACTAAATGGGACAGATTCTAGAGGCTCCTTCCCACAAGATGTTAACTTTTATATAAATAAAG  >  USA300TCH1516_ALE/2721343‑2721621

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