BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	USA300TCH1516_ALE	352,096	G→A	5.1%	K6K (AAG→AAA)	yezG_10 →	putative antitoxin YezG

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	USA300TCH1516_ALE	352,096	0	G	A	5.1%	375.7 / 6.0	118	K6K (AAG→AAA)	yezG_10	putative antitoxin YezG
Reads supporting (aligned to +/- strand): ref base G (57/55); new base A (3/3); total (60/58)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.37e-01

Only 100 of 118 total aligned reads displayed.

AATTATTATATGTATAAAAAATTTGGGGTTTTACCAGAAATGGAATATGAAATGGAAGAAATTAAAGAAATCGATCAATATATTAAAGAGCAAGATGAAGCTGAAATATAGGGGAGATAACATGACTTTCGAAGAAAAGCTAAGTCAAATGTACAATGAAATTGCAAATGAAATCAGTGGAATGATACCAGTAGAATGGGAAAAAGTATATACAATTGCCTACGTAGATGATGAAGGTGGAGAGGTTGTTTTTAATTATACTAAACCAGGAAGTGAA  >  USA300TCH1516_ALE/351958‑352234

aaatattatATGTATAAAAAATTTGGGGTTTTACCAGAAATGGAATATGAACCGGAAGGAAGTAAACAGATCGATCAATATATTAAAGAGCAAGAGGAAGCTGAAATATAGGGGAGATAACATGACTTTCGAAGAAAAGCt                                                                                                                                          <  1:1027783/138‑1 (MQ=17)

aaTTATTATATGTATAAAAAATTTGGGGTTTTACCAGAAATGGAATATGAAATGGAAGAAATTAAAGAAATCGATCAATATATTAAAGAGCAAGATGAAGCTGAAATATAGGGGAGATAACATGACTTTCGAAGAAAAGCt                                                                                                                                          <  2:878392/141‑1 (MQ=255)

   tattatATGTATAAAAAATTTGGGGTTTTACCAGAAATGGAATATGAAATGGAAGAAATTAAAGAAATCGATCAATATATTAAAGAGCAAGATGAAGCTGAAATATAGGGGAGATAACATGACTTTCGAAGAAAAGCTAAg                                                                                                                                       <  1:934016/141‑1 (MQ=255)

    attatATGTATAAAAAATTTGGGGTTTTACCAGAAATGGAATATGAAATGGAAGAAATTAAAGAAATCGATCAATATATTAAAGAGCAAGATGAAGCTGAAATATAGGGGAGATAACATGACTTTCGAAGAAAAGCTAAGt                                                                                                                                      >  2:987394/1‑141 (MQ=255)

     ttatATGTATAAAAAATTTGGGGTTTTACCAGAAATGGAATATGAAATGGAAGAAATTAAAGAAATCGATCAATATATTAAAGAGCAAGATGAAGCTGAAATATAGGGGAGATAACATGACATTCGAAGAAAAGCTAAgaa                                                                                                                                     >  2:658604/1‑139 (MQ=255)

      tataTGTATAAAAAATTTGGGGTTTTACCAGAAATGGAATATGAAATGGAAGAAATTAAAGAAATCGATCAATATATTAAAGAGCAAGATGAAGCTGAAATATAGGGGAGATAACATGACTTTCGAAGAAAAGCTAAGTCa                                                                                                                                    >  2:601885/1‑141 (MQ=255)

      tataTGTATAAAAAATTTGGGGTTTTACCAGAAATGGAATATGAAATGGAAGAAATTAAAGAAATCGATCAATATATTAAAGAGCAAGATGAAGCTGAAATATAGGGGAGATAACATGACTTTCGAAGAAAAGCTAAGTCa                                                                                                                                    >  3:61076/1‑141 (MQ=255)

        tatGTATAAAAAATTTGGGGTTTTACCAGAAATGGAATATGAAATGGAAGAAATTAAAGAAATCGATCAATATATTAAAGAGCAAGATGAAGCTGAAATATAGGGGAGATAACATGACTTTCGAAGAAAAGCTAAGTCaaa                                                                                                                                  <  2:433835/141‑1 (MQ=255)

             ataAAAAATTTGGGGTTTTACCAGAAATGGAATATGAAATGGAAGAAATTAAAGAAATCGATCAATATATTAAAGAGCAAGATGAAGCTGAAATATAGGGGAGATAACATGACTTTCGAAGAAAAGCTAAGTCAAATGTAc                                                                                                                             <  1:794317/141‑1 (MQ=255)

             ataAAAAATTTGGGGTTTTACCAGAAATGGAATATGAAATGGAAGAAATTAAAGAAATCGATCAATATATTAAAGAGCAAGATGAAGCTGAAATATAGGGGAGATAACATGACTTTCGAAGAAAAGCTAAGTCAAATGTAc                                                                                                                             <  1:741743/141‑1 (MQ=255)

               aaaaaaTTTGGGGTTTTACCAGAAATGGAATATGAAATGGAAGAAATTAAAGAAATCGATCAATATATTAAAGAGCAAGATGAAGCTGAAATATAGGGGAGATAACATGACTTTCGAAGAAAAGCTAAGTCAAATGTACaa                                                                                                                           >  3:418182/1‑141 (MQ=255)

                aaaaaTTTGGGGTTTTACCAGAAATGGAATATGAAATGGAAGAAATTAAAGAAATCGATCAATATATTAAAGAGCAAGATGAAGCTGAAATATAGGGGAGATAACATGACTTTCGAAGAAAAGCTAAGTCAAATGTACAAt                                                                                                                          >  2:65925/1‑141 (MQ=255)

                aaaaaTTTGGGGTTTTACCAGAAATGGAATATGAAATGGAAGAAATTAAAGAAATCGATCAATATATTAAAGAGCAAGATGAAGCTGAAATATAGGGGAGATAACATGACTTTCGAAGAAAAGCTAAGTCAAATGTACAAt                                                                                                                          >  2:1031286/1‑141 (MQ=255)

                  aaaTTTGGGGTTTTACCAGAAATGGAATATGAAATGGAAGAAATTAAAGAAATCGATCAATATATTAAAGAGCAAGATGAAGCTGAAATATAGGGGAGATAACATGACTTTCGAAGAAAAGCTAAGTCAAATGTACAATGa                                                                                                                        <  2:508837/141‑1 (MQ=255)

                  aaaTTTGGGGTTTTACCAGAAATGGAATATGAAATGGAAGAAATTAAAGAAATCGATCAATATATTAAAGAGCAAGATGAAGCTGAAATATAGGGGAGATAACATGACTTTCGAAGAAAAGCTAAGTCAAATGTACAATGa                                                                                                                        <  1:832468/141‑1 (MQ=255)

                   aaTTTGGGGTTTTACCAGAAATGGAATATGAAATGGAAGAAATTAAAGAAATCGATCAATATATTAAAGAGCAAGATGAAGCTGAAATATAGGGGAGATAACATGACTTTCGAAGAAAAGCTAAGTCAAATGTACAATGaa                                                                                                                       <  1:28024/141‑1 (MQ=255)

                     tttGGGGTTATACCAGAAATGGAATACGAAATGGAAGAAGTTAAAGAAATCGAGCAATATATTAAAGAGCAAGATGAAGCTGAACTATAGGGGCGATAATATGATTTTCGAAGAAAAACTAAGTCAAATGTACAATGAAAt                                                                                                                     >  1:746756/1‑141 (MQ=14)

                      ttGGGGTTTTACCAGAAATGGAATATGAAATGGAAGAAATTAAAGAAATCGATCAATATATTAAAGAGCAAGATGAAGCTGAAATATAGGGGAGATAACATGACTTTCGAAGAAAAGCTAAGTCAAATGTACAATGAAAtt                                                                                                                    >  4:427878/1‑141 (MQ=255)

                       tGGGGTTATACCAGAAATGGAATATGAAATGGAAGAAATTAAAGAAATCGAGCAATATATTAAAGAGCAAGATGAAGCTGAACTATAGGGGCGATAATATGACTTTCGAAGAAAAACTAAGTCAAATGTACAATGAAATTg                                                                                                                   <  2:90714/141‑1 (MQ=17)

                        ggggTTTTACCAGAAATGGAATATGAAATGGAAGAAATTAAAGAAATCGATCAATATATTAAAGAGCAAGATGAAGCTGAAATATAGGGGAGATAACATGACTTTCGAAGAAAAGCTAAGTCAAATGTACAATGAAATTGc                                                                                                                  <  4:203328/141‑1 (MQ=255)

                         gggTTTTACCAGAAATGGAATATGAAATGGAAGAAATTAAAGAAATCGATCAATATATTAAAGAGCAAGATGAAGCTGAAATATAGGGGAGATAACATGACTTTCGAAGAAAAGCTAAGTCAAATGTACAATGAAATTGCa                                                                                                                 >  1:155491/1‑141 (MQ=255)

                         gggTTTTACCAGAAATGGAATATGAAATGGAAGAAATTAAAGAAATCGATCAATATATTAAAGAGCAAGATGAAGCTGAAATATAGGGGAGATAACATGACTTTCGAAGAAAAGCTAAGTCAAATGTACAATGAAATTGCa                                                                                                                 >  1:476020/1‑141 (MQ=255)

                         gggTTTTACCAGAAATGGAATATGAAATGGAAGAAATTAAAGAAATCGATCAATATATTAAAGAGCAAGATGAAGCTGAAATATAGGGGAGATAACATGACTTCCGAAGAAAAGATAAGTCAAATGTACAATGAAATTGCa                                                                                                                 >  2:5057/1‑141 (MQ=255)

                          ggTTTTACCAAAAATGAAATATGAAATGGAAGAAATTAAAGAAATCAAACAATATAATAAAGAGCAAGATGAAGCTGAAATATAGGGGAGATAACATGACTTTCGAAGAAAAGCTAAGTCAAATGTACAATGaaaaaaaaa                                                                                                                <  2:565475/141‑7 (MQ=18)

                            ttttACCAGAAATGGAATATGAAATGGAAGAAATTAAAGAAATCGATCAATATATTAAAGAGCAAGATGAAGCTGAAATATAGGGGAGATAACATGACTTTCGAAGAAAAGCTAAGTCAAATGTACAATGAAATTGCAAAt                                                                                                              <  1:310409/141‑1 (MQ=255)

                            ttttACCAGAAATGGAATATGAAATGGAAGAAATTAAAGAAATCGATCAATATATTAAAGAGCAAGATGAAGCTGAAATATAGGGGAGATAACATGACTTTCGAAGAAAAGCTAAGTCAAATGTACAATGAAATTGCAAAt                                                                                                              >  2:576673/1‑141 (MQ=255)

                             tttACCAGAAATGGAATATGAAATGGAAGAAATTAAAGAAATCGATCAATATATTAAAGAGCAAGATGAAGCTGAAATATAGGGGAGATAACATGACTTTCGACGTAACGCTCAGTCAAATGTACAATGAAATTGCAAATg                                                                                                             <  4:363494/141‑1 (MQ=255)

                             tttACCAGAAATGGAATAGGAAATGGAAGAAATTAAAGAAATCGATCAATATAGTAAAGAGCAAGATGAAGCTGAAATATAGGGGAGATAACATGACTTTCGAAGAAAAGCTAAGTCAAATGTACAATGAAATTGCAAATg                                                                                                             <  4:167303/141‑1 (MQ=255)

                                aCCAGAAATGGAATATGAAATGGAAGAAATTAAAGAAATCGATCAATATATTAAAGAGCAAGATGAAGCTGAAATATAGGGGAGATAACATGACTTTCGAAGAAAAGCTAAGTCAAATGTACAATGAAATTGCAAATGaaa                                                                                                          <  1:1049277/141‑1 (MQ=255)

                                 ccAGAAATGGAATATGAAATGGAAGAAATTAAAGAAATCGATCAATATATTAAAGAGCAAGATGAAGCTGAAATATAGGGGAGATAACATGACTTTCGAAGAAAAGCTAAGTCAAATGTACAATGAAATTGCAAATGAAAt                                                                                                         <  4:368367/141‑1 (MQ=255)

                                   aGAAATGGAATATGAAATGGAAGAAATTAAAGAAATCGATCAATATATTAAAGAGCAAGATGAAGCTGAAATATAGGGGAGATAACATGACTTTCGAAGAAAAGCTAAGTCAAATGTACAATGAAATTGCAAATGAAATCa                                                                                                       <  2:338059/141‑1 (MQ=255)

                                    gAAATGGAATATGAAATGGAAGAAATTAAAGAAATCGATCAATATATTAAAGAGCAAGATGAAGCTGAAATATAGGGGAGATAACATGACTTTCGAAGAAAAGCTAAGTCAAATGTACAATGAAATTGCAAATGAAATCAg                                                                                                      >  1:106086/1‑141 (MQ=255)

                                      aaTGGAATATGAAATGGAAGAAATTAAAGAAATCGATCAATATATTAAAGAGCAAGATGAAGCTGAAATATAGGGGAGATAACATGACTTTCGAAGAAAAGCTAAGTCAAATGTACAATGAAATTGCAAATGAAATCAGTg                                                                                                    >  2:549891/1‑141 (MQ=255)

                                         ggAATATGAAATGGAAGAAATTAAAGAAATCGATCAATATATTAAAGAGCAAGATGAAGCTGAAATATAGGGGAGATAACATGACTTTCGAAGAAAAGCTAAGTCAAATGTACAATTAAATTGCAAATGAAATCAGTGaaa                                                                                                 >  4:301190/1‑141 (MQ=255)

                                         ggAATATGAAATGGAAGAAATTAAAGAAATCGATCAATATATTAAAGAGCAAGATGAAGCTGAAATATAGGGGAGATAACATGACTTTCGAAGAAAAGCTAAGTCAAATGTACAATGAAATTGCAAATGAAATCAGTGGaa                                                                                                 >  2:583679/1‑141 (MQ=255)

                                         ggAATATGAAATGGAAGAAATTAAAGAAATCGATCAATATATTAAAGAGCAAGATGAAGCTGAAATATAGGGGAGATAACATGACTTTCGAAGAAAAGCTAAGTCAAATGTACAATGAAATTGCAAATGAAATCAGTGGaa                                                                                                 <  2:428764/141‑1 (MQ=255)

                                             tatGAAATGGAAGAAATTAAAGAAATCGATCAATAGATTAAACAGCAAGATGAAGCTGAAATATAGGTGAGATAACATGACTTTCGAAGAAAAGCTAAGTCAAATGTACAATGAAATTGCAAATGAAATCAGTGGAATGAt                                                                                             >  1:932882/1‑141 (MQ=255)

                                              atGAAATGGAAGAAATTAAAGAAATCGATCAATATATTAAAGAGCAAGATGAAGCTGAAATATAGGGGAGATAACATGACTTTCGAAGAAAAGCTAAGTCAAATGTACAATGAAATTGCAAATGAAATCAGTGGAATGATa                                                                                            >  1:607819/1‑141 (MQ=255)

                                                 aaaTGGAAGAAATTAAAGAAATCGATCAATATATTAAAGAGCAAGATGAAGCTGAAATATAGAGGAGATAACATGACTTTCTAAGAAGAGCTAATTCAAATGTACAATGAATTGGCAAATGAAATCAGTGGAATGATACCa                                                                                         >  1:542046/1‑141 (MQ=18)

                                                  aaTGGAAGAAATTAAAGAAATCGATCAATATATTAAAGAGCAAGATGAAGCTGAAATATAGGGGAGATAACATGACTTTCGAAGAAAAGCTAAGTCAAATGTACAATGAAATTGCAAATGAAATCAGTGGAATGATACCAg                                                                                        <  3:95166/141‑1 (MQ=255)

                                                       aagaaATTAAAGAAATCGATCAATATATTAAAGAGCAAGATGAAGCTGAAATATAGGGGAGATAACATGACTTTCGAAGAAAAGCTAAGTCAAATGTACAATGAAATTGCAAATGAAATCAGTGGAATGATACCAGTAGaa                                                                                   <  3:35905/141‑1 (MQ=255)

                                                       aagaaATTAAAGAAATCGAGCAATATATTAAAGAGCAAGATGAAGCTGAACTATAGGGGCGATAATATGACTTTCGAAGAAAAACTAAGTCAAATGTACAATGAAATTGCAAATAAAATTAGCAGTATGATACCGGTAGaa                                                                                   <  2:786836/141‑1 (MQ=11)

                                                        agaaATTAAAGAAATCGATCAATATATTAAAGAGCAAGATGAAGCTGAAATATAGGGGAGATAACATGACTTTCGAAGAAAAGCTAAGTCAAATGTACAATGAAATTGCAAATGAAATCAGTGGAATGATACCAGTAGAAt                                                                                  <  4:394722/141‑1 (MQ=255)

                                                        agaaATTAAAGAAATCGATCAATATATTAAAGAGCAAGATGAAGCTGAAATATAGGGGAGAAAACATGACTTGCGAAGAAAAGCTAAGTCAAATGTACAATGAAATTGCAAATGAAATCAGTGTTATGATACCAgaaaaaa                                                                                  <  2:213281/141‑7 (MQ=18)

                                                         gaaATTAAAGAAATCGATCAATATATTAAAGAGCAAGATGAAGCTGAAATATAGGGGAGATAACATGACTTTCGAAGAAAAGCTAAGTCAAATGTACAATGAAATTGCAAATGAAATCAGTGGAATGATACCAGTAGAATg                                                                                 >  2:128586/1‑141 (MQ=255)

                                                         gaaATTAAAGAAATCGATCAATATATTAAAGAGCAAGATGAAGCTGAAATATAGGGGAGATAACATGACTTTCGAAGAAAAGCTAAGTCAAATGTACAATGAAATTGCAAATGAAATCAGTGGAATGATACCAGTAGAATg                                                                                 >  2:801303/1‑141 (MQ=255)

                                                          aaaTTAAAGAAATCGATCAATATATTAAAGAGCAAGATGAAGCTGAAATATAGGGGAGATAACATGACTTTCGAAGAAAAGCTAAGTCAAATGTACAATGAAATTGCAAATGAAATCAGTGGAATGATACCAGTAGAATgg                                                                                <  2:694340/141‑1 (MQ=255)

                                                           aaTTAAAGAAATCGATCAATATATTAAAGAGCAAGATGAAGCTGAAATATAGGGGAGATAACATGACTTTCGAAGAAAAGCTAAGTCAAATGTACAATGAAATTGCAAATGAAATCAGTGGAATGATACCAGTAGAATggg                                                                               <  1:817341/141‑1 (MQ=255)

                                                            aTTAAAGAAATCGATCAATATATTAAAGAGCAAGATGAAGCTGAAATATAGGGGAGATAACATGACTTTCGAAGAAAAGCTAAGTCAAATGTACAATGAAATTGCAAATGAAATCAGTGGAATGATACCAGTAGAATGGGa                                                                              <  2:157067/141‑1 (MQ=255)

                                                             ttAAAGAAATCGATCAATATATTAAAGAGCAAGATGAAGCTGAAATATAGGGGAGATAACATGACTTTCGAAGAAAAGCTAAGTCAAATGTACAATGAAATTGCAAATGAAATCAGTGGAATGATACCAGTAGAATGGGaa                                                                             <  1:405430/141‑1 (MQ=255)

                                                               aaaGAAATCGATCAATATATTAAAGAGCAAGATGAAGCTGAAATATAGGGGAGATAACATGACTTTCGAAGAAAAGCTAAGTCAAATGTACAATGAAATTGCAAATGAAATCAGTGGAATGATACCAGTAGAATGGGaaaa                                                                           >  1:1062270/1‑141 (MQ=255)

                                                                aaGAAATCGATCAATATATTAAAGAGCAAGATGAAGCTGAAATATAGGGGAGATAACATGACTTTCGAAGAAAAGCTAAGTCAAATGTACAATGAAATTGCAAATGAAATCAGTGGAATGATACCAGTAGAATGGGaaaaa                                                                          >  2:936801/1‑141 (MQ=255)

                                                                aaGAAATCGATCAATATATTAAAGAGCAAGATGAAGCTGAAATATAGGGGAGATAACATGACTTTCGAAGAAAAGCTAAGTCAAATGTACAATGAAATTGCAAATGAAATCAGTGGAATGATACCAGTAGAATGGGaaaaa                                                                          >  2:93917/1‑141 (MQ=255)

                                                                 aGAAATCGATCAATATATTAAAGAGCAAGATGAACCTGAAATATAGGGGAGATAACATGACTTTCGAAGAAAAGCTAAGTCAAATGTACAATGAAATTGCAAATGAAATCAGTGGAATGATACCAGTAGAATGGGAAAAAg                                                                         <  1:1007281/141‑1 (MQ=255)

                                                                    aaTCGATCAATATATTAAAGAGCAAGATGAAGCTGAAATATAGGGGAGATAACATGACTTTCGAAGAAAAGCTAAGTCAAATGTACAATGAAATTGCAAATGAAATCAGTGGAATGATACCAGTAGAATGGGAAAAAGtat                                                                      >  2:487188/1‑141 (MQ=255)

                                                                     aTCGATCAATATATTAAAGAGCAAGATGAAGCTGAAATATAGGGGAGATAACATGACTTTCGAAGAAAAGCTAAGTCAAATGTACAATGAAATTGCAAATGAAATCAGTGGAATGATACCAGTAGAATGGGAAAAAGtata                                                                     <  1:34385/141‑1 (MQ=255)

                                                                      aagATCAATATATTAAAGAGCAAGATGAAGCTGAAATATAGGGGAGATAACATGAATTTCGAAGAAAAGCTAAGTCAAATGTACAATGAAATTGCAAATGAAATCAGTGGAATGATACCAGTAGAATGGGAAAAAGtatat                                                                    <  2:510237/139‑1 (MQ=255)

                                                                       cGATCAATATATTAAAGAGCAAGATGAAGCTGAAATATAGGGGAGATAACATGACTTTCGAAGAAAAGCTAAGTCAAATGTACAATGAAATTGCAAATGAAATCAGTGGAATGATACCAGTGGAATGGGAAAAAGtatata                                                                   >  3:229378/1‑141 (MQ=255)

                                                                         aTCAATATATTAAAGAGCAAGATGAAGCTGAAATATAGGGGAGATAACATGACTTTCGAAGAAAAGCTAAGTCAAATGTACAATGAAATTGCAAATGAAATCAGTGGAATGATACCAGTAGAATGGGAAAAAGTATATACa                                                                 >  1:697322/1‑141 (MQ=255)

                                                                          tCAATATATTAAAGAGCAAGATGAAGCTGAAATATAGGGGAGATAACATGACTTTCGAAGAAAAGCTAAGTCAAATGTACAATGAAATTGCAAATGAAATCAGTGGAATGATACCAGTAGAATGGGAAAAAGTATATACaa                                                                >  4:267217/1‑141 (MQ=255)

                                                                            aaTATATGAAAGAGCAAGCTGAAGCTGAACTATAGGTGCGATAATATTATTTTCGAAGAAAAACTAAGTCAAATGTACAATGATATTGCAAATGAAATCAGCGGGATGATACCAGTAGAGTGGGAAAAAGTATATACAAtt                                                              >  2:676171/1‑141 (MQ=14)

                                                                              tatatTAAAGAGCAAGATGAAGCTGAAATATAGGGGAGATAACATGACTTTCGAAGAAAAGCTAATTCAAATGTACAATGAAATTGCAAATGAAATCAGTGGAATGATACCAGTAGAATGGGAAAAAGTATATACAATTGc                                                            >  4:397673/1‑141 (MQ=255)

                                                                               atatTAAAGAGCAAGATGAAGCTGAAATATAGGGGAGATAACATGACTTTCGAAGAAAAGCTAAGTCAAATGTACAATGAAATTGCAAATGAAATCAGTGGAATGATACCAGTAGAATGGGAAAAAGTATATACAATTGcc                                                           <  4:248945/141‑1 (MQ=255)

                                                                                  ttAAAGAGCAAGATGAAGCTGAAATATAGGGGAGATAACATGACTTTCGAAGAAAAGCTAAGTCAAATGTACAATGAAATTGCAAATGAAATCAGTGGAATGATACCAGTAGAATGGGAAAAAGTATATACAATTGCCTAc                                                        >  2:548398/1‑141 (MQ=255)

                                                                                      agagCAAGATGAAGCTGAAATATAGGGGAGATAACATGACTTTCGCAGAGAAGCTAAGTCAAATGTACAATGAAATTGCAAATGAAATCGGTGGAATGATACCAGTAGAATGGGAAAAAGTATATACAATTGCCTACGTAg                                                    >  1:792017/1‑141 (MQ=255)

                                                                                      agagCAAGATGAAGCTGAAATATAGGGGAGATAACATGACTTTCGAAGAAAAGCTAAGTCAAATGTACAATGAAATTGCAAATGAAATCAGTGGAATGATACCAGTAGAATGGGAAAAAGTATATACAATTGCCTACGTAg                                                    <  1:1040329/141‑1 (MQ=255)

                                                                                      agagCAAGATGAAGCTGAAATATAGGGGAGATAACATGACTTTCGAAGAAAAGCTAAGTCAAATGTACAATGAAATTGCAAATGAAATCAGTGGAATGATACCAGTAGAATGGGAAAAAGTATATACAATTGCCTACGTAg                                                    >  2:847423/1‑141 (MQ=255)

                                                                                         gCAAGATGAAGCTGAAATATAGGGGAGATAACATGACTTTCGAAGAAAAGCTAAGTCAAATGTACAATGAAATTGCAAATGAAATCAGTGGAATGATACCAGTAGAATGGGAAAAAGTATATACAATTGCCTACGTAgatg                                                 >  1:661404/1‑141 (MQ=255)

                                                                                            aGATGAAGCTGAAATATAGGGGAGATAACATGACTTTCGAAGAAAAGCTAAGTCAAATGTACAATGAAATTGCAAATGAAATCAGTGGAATGATACCAGTAGAATGGGAAAAAGTATATACAATTGCCTACGTAGatgatg                                              <  1:71821/141‑1 (MQ=255)

                                                                                             gATGAAGCTGAAATATAGGGGAGATAACATGACTTTCGAAGAAAAGCTAAGTCAAATGTACAATGAAATTGCAAATGAAATCAGTGGAATGATACCAGTAGAATGGGAAAAAGTATATACAATTGCCTACGTAGAtgatga                                             <  4:114319/141‑1 (MQ=255)

                                                                                              aTGAAGCTGAAATATAGGGGAGATAACATGACTTTCGAAGAAAAGCTAAGTCAAATGTACAATGAAATTGCAAATGAAATCAGTGGAATGATACCAGTAGAATGGGAAAAAGTATATACAATTGCCTACGTAGATGATGaa                                            >  4:181091/1‑141 (MQ=255)

                                                                                              aTGAAGCTGAAATATAGGGGAGATAACATGACTTTCGAAGAAAAGCTAAGTCAAATGTACAATGAAATTGCAAATGAAATCAGTGGAATGATACCAGTAGAATGGGAAAAAGTATATACAATTGCCTACGTAGATGATGaa                                            >  2:533320/1‑141 (MQ=255)

                                                                                              aTGAAGCTGAAATATAGGGGAGATAACATGACTTTCGAAGAAAAGCTAAGTCAAATGTACAATGAAATTGCAAATGAAATCAGTGGAATGATACCAGTAGAATGGGAAAAAGTATATACAATTGCCTACGTAGATGATGaa                                            >  2:372638/1‑141 (MQ=255)

                                                                                               tGAAGCTGAAATATAGGGGAGATAACATGACTTTCGAAGAAAAGCTAAGTCAAATGTACAATGAAATTGCAAATGAAATCAGTGGAATGATACCAGTAGAATGGGAAAAAGTATATACAATTGCCTACGTAGATGATGAAg                                           <  4:262461/141‑1 (MQ=255)

                                                                                                   gCTGAACTATAGGGGCGATAATATGATTTTCGAAGAAAAACTAAGTCAAATGTACAATGAAATTGCAAATGAAATCAGTGGGATGATACCAGTAGAGTGGGAAAAAGTATATACAATTGCTTATGTAGATGATGAAGGTgg                                       >  3:114946/1‑141 (MQ=17)

                                                                                                   gCTGAAATATAGGGGAGATAACATGACTTTCGAAGAAAAGCTAAGTCAAATGTACAATGAAATTGCAAATGAAATCAGTGGAATGATACCAGTAGAATGGGAAAAAGTATATACAATTGCCTACGTAGATGATGAAGGTgg                                       >  3:159575/1‑141 (MQ=255)

                                                                                                      gAAATATAGGGGAGATAACATGACTTTCGAAGAAAAGCTAAGTCAAATGTACAATGAAATTGCAAATGAAATCAGTGGAATGATACCAGTAGAATGGGAAAAAGTATATACAATTGCCTACGTAGATGATGAAGGTGgaga                                    >  2:801390/1‑141 (MQ=255)

                                                                                                        aaTATAGGGGAGATAACATGACTTTCGAAGAAAAGCTAAGTCAAATGTACAATGAAATTGCAAATGAAATCAGTGGAATGATACCAGTAGAATGGGAAAAAGTATATACAATTGCCTACGTAGATGATGAAGGTGGAGAgg                                  <  1:609559/141‑1 (MQ=255)

                                                                                                         atatAGGGGAGATAACATGACTTTCGAAGAAAAGCTAAGTCAAATGTACAATGAAATTGCAAATGAAATCAGTGGAATGATACCAGTAGAATGGGAAAAAGTATATACAATTGCCTACGTAGATGATGAAGGTGGAGAGgt                                 <  3:365429/141‑1 (MQ=255)

                                                                                                         atatAGGGGAGATAACATGACTTTCGAAGAAAAGCTAAGTCAAATGTACAATGAAATTGCAAATGAAATCAGTGGAATGATACCAGTAGAATGGGAAAAAGTATATACAATTGCCTACGTAGATGATGAAGGTGGAGAGgt                                 <  2:111506/141‑1 (MQ=255)

                                                                                                          tataGGGGAGATAACATGACTTTCGAAGAAAAGCTAAGTCAAATGTACAATGAAATTGCAAATGAAATCAGTGGAATGATACCAGTAGAATGGGAAAAAGTATATACAATTGCCTACGTAGATGATGAAGGTGGAGAGgtt                                >  3:11148/1‑141 (MQ=255)

                                                                                                            taGGGGAGATAACATGACTTTCGAAGAAAAGCTAAGTCAAATGTACAATGAAATTGCAAATGAAATCAGTGGAATGATACCAGTAGAATGGGAAAAAGTATATACAATTGCCTACGTAGATGATGAAGGTGGAGAGgttgt                              <  2:155937/141‑1 (MQ=255)

                                                                                                              ggggAGATAACATGACTTTCGAAGAAAAGCTAAGTCAAATGTACAATGAAATTGCAAATGAAATCAGTGGAATGATACCAGTAGAATGGGAAAAAGTATATACAATTGCCTACGTAGATGATGAAGGTGGAGAGGTTGttt                            <  2:840606/141‑1 (MQ=255)

                                                                                                              ggggAGATAACATGACTTTCGAAGAAAAGCTAAGTCAAATGTACAATGAAATTGCAAATGAAATCAGTGGAATGATACCAGTAGAATGGGAAAAAGTATATACAATTGCCTACGTAGATGATGAAGGTGGAGAGGTTGttt                            >  1:134427/1‑141 (MQ=255)

                                                                                                               gggAGATAACATGACTTTCGAAGAAAAGCTAAGTCAAATGTACAATGAAATTGCAAATGAAATCAGTGGAATGATACCAGTAGAATGGGAAAAAGTATATACAATTGCCTACGTAGATGATGAAGGTGGAGAGGTTGtttt                           >  2:410412/1‑141 (MQ=255)

                                                                                                                  agaTAACATGACTTTCGAAGAAAAGCTAAGTCAAATGTACAATGAAATTGCAAATGAAATCAGTGGAATGATACCAGTAGAATGGGAAAACGTATATACAATTGCCTACGTAGATGATGAAGGTGGAGAGGTTGTTTTTaa                        >  4:195723/1‑141 (MQ=255)

                                                                                                                  agaTAACATGACTTTCGAAGAAAAGCTAAGTCAAATGTACAATGAAATTGCAAATGAAATCAGTGGAATGATACCAGTAGAATGGGAAAAAGTATATACAATTGCCTACGTAGATGATGAAGGTGGAGAGGTTGTTTTTaa                        >  1:6312/1‑141 (MQ=255)

                                                                                                                      aaCATTACTTTCGAAGAAAAGCTAAGTCAAATGTACAATGAAATTGCAAATGAAATCAGTGGAATGATACCAGTAGAATGGGAAAAAGTATATACAATTGCCTACGTAGATGATGAAGGTGGAGAGGTTCTTTTTAATtaa                    >  4:128299/1‑140 (MQ=255)

                                                                                                                      aaCATGACTTTCGAAGAAAAGCTAAGTCAAATGTACAATGAAATTGCAAATGAAATCAGTGGAATGATACCAGTAGAATGGGAAAAAGTATATACAATTGCCTACGTAGATGATGAAGGTGGAGAGGTTGTTTTTAATtat                    >  1:747141/1‑141 (MQ=255)

                                                                                                                         aTGACTTTCGAAGAAAAGCTAAGTCAAATGTACAATGAAATTGCAAATGAAATCAGTGGAATGATACCAGTAGAATGGGAAAAAGTATATACAATTGCCTACGTAGATGATGAAGGTGGAGAGGTTGTTTTTAATTATACt                 >  2:64273/1‑141 (MQ=255)

                                                                                                                          tGACTTTCGAAGAAAAGCTAAGTCAAATGTACAATGAAATTGCAAATGAAATCAGTGGAATGATACCAGTAGAATGGGAAAAAGTATATACAATTGCCTACGTAGATGATGAAGGGGGAGAGGTTGTTTTTAATTATACTa                <  2:611151/141‑1 (MQ=255)

                                                                                                                               ttCGAAGAAAAGCTAAGTCAAATGTACAATGAAATTGCAAATGAAATCAGTGGAATGATACCAGTAGAATGGGAAAAAGTATATACAATTGCCTACGTAGATGATGAAGGTGGAGAGGTTGTTTTTAATTATACTAAACCa           >  4:290327/1‑141 (MQ=255)

                                                                                                                               ttCGAAGAAAAGCTAAGTCAAATGTACAATGAAATTGCAAATGAAATCAGTGGAATGATACCAGTAGAATGGGAAAAAGTATATACAATTGCCTACGTAGATGATGAAGGTGGAGAGGTTGTTTTTAATTATACTAAACCa           <  1:441693/141‑1 (MQ=255)

                                                                                                                                  gaagaaAAGCTAAGTCAAATGTACAATGAAATTGCAAATGAAATCAGTGGAATGATACCAGTAGAATGGGAAAAAGTATATACAATTGCCTACGTAGATGATGAAGGTGGAGAGGTTGTTTTTAATTATACTAAACCAGGa        >  1:315578/1‑141 (MQ=255)

                                                                                                                                  gaagaaAAGCTAAGTAAAATGTACAATGAAATTGCAAATGAAATCAGTGGAATGATACCAGGAGAATGGGAAAAAGTATATACAATTGCCTACGTAGATGATGAAGGTGGAGAGGTTGTTTTTAATTATACTAAACCAGGa        <  2:1007554/141‑1 (MQ=255)

                                                                                                                                      aaaaGCTAAGTCAAATGTACAATGAAATTGCAAATGAAATCAGTGGAATGATACCAGTAGAATGGGAAAAAGTATATACAATTGCCTACGTAGATGATGAAGGTGGAGAGGTTGTTTTTAATTATACTAAACCAGGAAGTg    <  2:476020/141‑1 (MQ=255)

                                                                                                                                       aaaGCTAAGTCAAATGTACAATGAAATTGCAAATGAAATCAGTGGAATGATACCAGTAGAATGGGAAAAAGTATATACAATTGCCTACGTAGATGATGAAGGTGGAGAGGTTGTTTTTAATTATACTAAACCAGGAAGTGa   >  2:61304/1‑141 (MQ=255)

                                                                                                                                       aaaGCTAAGTCAAATGTACAATGAAATTGCAAATGAAATCAGTGGAATGATACCAGTAGAATGGGAAAAAGTATATACAATTGCCTACGTAGATGATGAAGGTGGAGAGGTTGTTTTTAATTATACTAAACCAGGAAGTGa   >  1:315091/1‑141 (MQ=255)

                                                                                                                                        aaGCTAAGTCAAATGTACAATGAAATTGCAAATGAAATCAGTGGAATGATACCAGTAGAATGGGAAAAAGTATATACAATTGCCTACGTAGATGATGAAGGTGGAGAGGTTGTTTTTAATTATACTAAACCAGGAAGTGaa  <  4:283992/141‑1 (MQ=255)

                                                                                                                                        aaGCTAAGTCAAATGTACAATGAAATTGCAAATGAAATCAGTGGAATGATACCAGTAGAATGGGAAAAAGTATATACAATTGCCTACGTAGATGATGAAGGTGGAGAGGTTGTTTTTAATTATACTAAACCAGGAAGTGaa  >  2:409782/1‑141 (MQ=255)

AATTATTATATGTATAAAAAATTTGGGGTTTTACCAGAAATGGAATATGAAATGGAAGAAATTAAAGAAATCGATCAATATATTAAAGAGCAAGATGAAGCTGAAATATAGGGGAGATAACATGACTTTCGAAGAAAAGCTAAGTCAAATGTACAATGAAATTGCAAATGAAATCAGTGGAATGATACCAGTAGAATGGGAAAAAGTATATACAATTGCCTACGTAGATGATGAAGGTGGAGAGGTTGTTTTTAATTATACTAAACCAGGAAGTGAA  >  USA300TCH1516_ALE/351958‑352234

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