BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	USA300TCH1516_ALE	2,120,639	A→G	6.8%	L46L (TTA→CTA)	USA300HOU_RS10785 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	USA300TCH1516_ALE	2,120,639	0	A	G	6.8%	343.2 / 12.4	117	L46L (TTA→CTA)	USA300HOU_RS10785	hypothetical protein
Reads supporting (aligned to +/- strand): ref base A (56/53); new base G (5/3); total (61/56)
Fisher's exact test for biased strand distribution p-value = 7.19e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.36e-01

Only 100 of 120 total aligned reads displayed.

TCTCATGCTTATACCTACCATTTCATGACTAAGTTCATTAGTTTGTCCTGTTCATCTGTGTTATTTTCAATCCATTCATAAACACTTTGTTTCAAAATATCTAAAGCTGTGTATAGATCGTTCTCATCAGAAACTAGTAACCCGTCAATTGAATTTCCTTCATGGTCTAGAACAACCATTTCGACGCTATACGCTCGTTTCTTAACTCTTAATTG‑AAAATCAAAGCCATCTACATTAATTATTTTTCGACATACGTCACCCGTTTTGTAATACATTGTT  >  USA300TCH1516_ALE/2120500‑2120778

tctcATGCTTATACCTACCATTTCATGACTAAGTTCATTAGTTTGTCCTGTTCATCTGTGTTATTTTCAATCCATTCATAAACACTTTGTTTCAAAATATCTACAGCTGTGTATAGATCGTTCTCATCAGAAACTAGTAAc                                                                                                                                             <  1:784846/141‑1 (MQ=255)

tctcATGCTTATACCTACCATTTCATGACTAAGTTCATTAGTTTGTCCTGTTCATCTGTGTTATTTTCAATCCATTCATAAACACTTTGTTTCAAAATATCTAAAGCTGTGTATAGATCGTTCTCATCAGAAACTAGTAAc                                                                                                                                             <  1:989805/141‑1 (MQ=255)

tctcATGCTTATACCTACCATTTCATGACTAAGTTCATTAGTTTGTCCTGTTCATCTGTGTTATTTTCAATCCATTCATAAACACTTTGTTTCAAAATATCTAAAGCTGTGTATAGATCGTTCTCATCAGAAACTAGTAAc                                                                                                                                             >  2:700825/1‑141 (MQ=255)

 ctcATGCTTATACCTACCATTTCATGACTAAGTTCATTAGTTTGTCCTGTTCATCTGTGTTATTTTCAATCCATTCATAAACACTTTGTTTCAAAATATCTAAAGCTGTGTATAGATCGTTCTCATCAGAAACTAGTAAcc                                                                                                                                            >  1:748148/1‑141 (MQ=255)

   cATGCTTATACCTACCATTTCATGACTAAGTTCATTAGTTTGTCCTGTTCATCTGTGTTATTTTCAATCCATTCATAAACACTTTGTTTCAAAATATCTAAAGCTGTGTATAGATCGTTCTCATCAGAAACTAGTAACCCg                                                                                                                                          >  2:160075/1‑141 (MQ=255)

    aTGCTTATACCTACCATTTCATGACTAAGTTCATTAGTTTGTCCTGTTCATCTGTGTTATTTTCAATCCATTCATAAACACTTTGTTTCAAAATATCTAAAGCTGTGTATAGATCGTTCTCATCAGAAACTAGTAACCCGt                                                                                                                                         >  2:909957/1‑141 (MQ=255)

        ttATGCCTACCAGTTCATGACTAAGTTCATTAGTTTGTCCTGTTCATCTGTGTTATTTTCAATCGATTCATAAACACTTTGTTTCAAAATATCTAAAGCTGTGTATAGATCGTTCTCATCAGAAACTAGTAACCCGTCAAt                                                                                                                                     <  4:389994/141‑1 (MQ=255)

        ttATACCTACCATTTCATGACTAAGTTCATTAGTTTGTCCTGTTCATCTGTGTTATTTTCAATCCATTCATAAACACTTTGTTTCAAAATATCTAAAGCTGTGTATAGATCGTTCTCATCAGAAACTAGTAACCCGTCAAt                                                                                                                                     <  4:139155/141‑1 (MQ=255)

         tataCCTACCATTTCATGACTAAGTTCATTAGTTTGTCCTGTTCATCTGTGTTATTTTCAATCCATTCATAAACACTTTGTTTCAAAATATCTAAAGCTGTGTATAGATCGTTCTCATCAGAAACTAGTAACCCGTCAAtt                                                                                                                                    <  3:318603/141‑1 (MQ=255)

          ataCCTACCATTTCATGACTAAGTTCATTAGTTTGTCCTGTTCATCTGTGTTATTTTCAATCCATTCATAAACACTTTGTTTCAAAATATCTAAAGCTGTGTATAGATCGTTCTCATCAGAAACTAGTAACCCGTCAATTg                                                                                                                                   >  1:522359/1‑141 (MQ=255)

           tacctaccATTTCATGACTAAGTTCATTAGTTTGTCCTGTTCATCTGTGTTATTTTCAATCCATTCATAAACACTTTGTTTCAAAATATCTAAAGCTGTGTATAGATCGTTCTCATCAGAAACCAGTAACCCGTCAATTGa                                                                                                                                  >  1:603800/1‑141 (MQ=255)

                 ccATTTCATGACTAAGTTCATTAGTTTGTCCTGTTCATCTGTGTTATTTTCAATCCATTCATAAACACTTTGTTTCAAAATATCTAAAGCTGTGTATAGATCGTTCTCATCAGAAACTAGTAACCCGTCAATTGAATTTcc                                                                                                                            <  4:370895/141‑1 (MQ=255)

                   aTTTCATGACTAAGTTCATTAGTTTGTCCTGTTCATCTGTGTTATTTTCAATCCATTCATAAACACTTTGTTTCAAAATATCTAAAGCTGTGTATAGATCGTTCTCATCAGAAACTAGTAACCCGTCAATTGAATTTCCtt                                                                                                                          <  4:197905/141‑1 (MQ=255)

                   aTTTCATGACTAAGTTCATTAGTTTGTCCTGTTCATCTGTGTTATTTTCAATCCATTCATAAACACTTTGTTTCAAAATATCTAAAGCTGTGTATAGATCGTTCTCATCAGAAACTAGTAACCCGTCAATTGAATTTCCtt                                                                                                                          <  2:1030522/141‑1 (MQ=255)

                    tttCATGACTAAGTTCATTAGTTTGTCCTGTTCATCTGTGTTATTTTCAATCCATTCATAAACACTTTGTTTCAAAATATCTAAAGCTGTGTATAGATCGTTCTCATCAGAAACTAGTAACCCGTCAATTGAATTTCCTTc                                                                                                                         <  2:890044/141‑1 (MQ=255)

                      tCATGACTAAGTTCATTAGTTTGTCCTGTTCATCTGTGTTATTTTCAATCCATTCATAAACACTTTGTTTCAAAATATCTAAAGCTGTGTATAGATCGTTCTCATCAGAAACTAGTAACCCGTCAATTGAATTTCCTTCAt                                                                                                                       <  2:385348/141‑1 (MQ=255)

                       cATGACTAAGTTCATTAGTTTGTCCTGTTCATCTGTGTTATTTTCAATCCATTCATAAACACTTTGTTTCAAAATATCTAAAGCTGTGTATAGATCGTTCTCATCAGAAACTAGTAACCCGTCAATTGAATTTCCTTCATg                                                                                                                      <  1:886707/141‑1 (MQ=255)

                       cATGACTAAGTTCATTAGTTTGTCCTGTTCATCTGTGTTATTTTCAATCCATTCATAAACACTTTGTTTCAAAATATCTAAAGCTGTGTATAGATCGTTCTCATCAGAAACTAGTAACCCGTCAATTGAATTTCCTTCATg                                                                                                                      >  2:723607/1‑141 (MQ=255)

                       cATGACTAAGTTCATTAGTTTGTCCTGTTCATCTGTGTTATTTTCAATCCATTCATAAACACTTTGTTTCAAAATATCTAAAGCTGTGTATAGATCGTTCTCATCAGAAACTAGTAACCCGTCAATTGAATTTCCTTCATg                                                                                                                      >  3:295099/1‑141 (MQ=255)

                        aTGACTAAGTTCATTAGTTTGTCCTGTTCATCTGTGTTATTTTCAATCCATTCATAAACACTTTGTTTCAAAATATCTAAAGCTGTGTATAGATCGTTCTCATCAGAAACTAGTAACCCGTCAATTGAATTTCCTTCATgg                                                                                                                     >  2:587370/1‑141 (MQ=255)

                           aCTAAGTTCATTAGTTTGTCCTGTTCATCTGTGTTATTTTCAATCCATTCATAAACACTTTGTTTCAAAATATCTAAAGCTGTGTATAGATCGTTCTCATCAGAAACTAGTAACCCGTCAATTGAATTTCCTTCATGGTCt                                                                                                                  <  1:592804/141‑1 (MQ=255)

                           aCTAAGTTAATTAGTCTGTCCTGTTCATCTGTGTTATTTTCAATCCATTCATAAATACTTTGTTTCAAAATATCTAAAGCTGTGTATAGATCGTTCTCGTCAGAAACTAGTAGCCCGTCAATTGAATTTCCTTCATGATCt                                                                                                                  <  2:1076496/141‑1 (MQ=38)

                             tAAGTTCATTAGTTTGTCCTGTTCATCTGTGTTATTTTCAATCCATTCATAAACACTTTGTTTCAAAATATCTAAAGCTGTGTATAGATCGTTCTCATCAGAAACTAGTAACCCGTCAATTGAATTTCCTTCATGGTCTAg                                                                                                                <  4:117197/141‑1 (MQ=255)

                               aGTTCATTAGTTTGTCCTGTTCATCTGTGTTATTTTCAATCCATTCATAAACACTTTGTTTCAAAATATCTAAAGCTGTGTATAGATCGTTCTCATCAGAAACTAGTAACCCGTCAATTGAATTTCCTTCATGGTCTAGaa                                                                                                              <  1:503076/141‑1 (MQ=255)

                                gttaaTTAGTCTGTCCTGTTCATCTGTGTTATTTTCAATCCATTCATAAATACTTTGTTTCAAAATATCTAAAGCTGTGTATAGATCGTTCTCGTCAGAAACTAGTAGCCCGTCAATTGAATTTCCTTCATGATCTAAaac                                                                                                             >  2:1080817/5‑141 (MQ=38)

                                  tCATTAGTTTGTCCTGTTCATCTGTGTTATTTTCAATCCATTCGTAAACACTTTGTTTCAAAATATCTAAAGCTGTGTATAGATCGTTCTCATCAGAAACTAGTAACCCGTCAATTGAATTTCCTTCATGGTCTAGaacaa                                                                                                           <  1:852966/141‑1 (MQ=255)

                                   cATTAGTTTGTCCTGTTCATCTGTGTTATTTTCAATCCATTCATAAACACTTTGTTTCAAAATATCTAAAGCTGTGTATAGATCGTTCTCATCAGAAACTAGTAACCCGTCAATTGAATTTCCTTCATGGTCTTGaacaac                                                                                                          <  3:402656/141‑1 (MQ=255)

                                       aGTTTGTCCTGTTCATCTGTGTTATTTTCAATCCATTCATAAACACTTTGTTTCAAAATATCTAAAGCTGTGTATAGATCGTTCTCATCGGAAACTAGTAACCCGTCAATTGAATTTCCTTCATGGTCTAGAACAACCAtt                                                                                                      <  1:281820/141‑1 (MQ=255)

                                        gTTTGTCCTGTTCATCTGTGTTATTTTCAATCCATTCATAAACACTTTGTTTCAAAATATCTAAAGCTGTGTATAGATCGTTCTCATCAGAAACTAGTAACCCGTCAATTGAATTTCCTTCATGGTCTAGAACAACCAttt                                                                                                     >  1:872486/1‑141 (MQ=255)

                                            gTCCTGTTCATCTGTGTTATTTTCAATCCATTCATAAACACTTTGTTTCAAAATATCTAAAGCTGTGTATAGATCGTTCTCATCAGAAACTAGTAACCCGTCAATTGAATTTCCTTCATGGTCTAGAACAACCATTTCGAc                                                                                                 <  3:286398/141‑1 (MQ=255)

                                             tCCTGTTCATCTGTGTTATTTTCAATCCATTCATAAACACTTTGTTTCAAAATATCTAAAGCTGTGTATAGATCGTTCTCATCAGAAACTAGTAACCCGTCAATTGAATTTCCTTCATGGTCTAGAACAACCATTTCGACg                                                                                                >  2:1065026/1‑141 (MQ=255)

                                             tCCTGTTCATCTGTGTTATTTTCAATCCATTCATAAACACTTTGTTTCAAAATATCTAAAGCTGTGTATAGATCGTTCTCATCAGAAACTAGTAACCCGTCAATTGAATTTCCTTCATGGTCTAGAACAACCATTTCGACg                                                                                                >  3:273068/1‑141 (MQ=255)

                                             tCCTGTTCATCTGTGTTATTTTCAATCCATTCATAAACACTTTGTTTCAAAATATCTAAAGCTGTGTATAGATCGTTCTCATCAGAAACTAGTAACCCGTCAATTGAATTTCCTTCATGGTCTAGAACAACCATTTCGACg                                                                                                >  4:238133/1‑141 (MQ=255)

                                              ccTGTTCATCTGTGTTATTTTCAATCCATTCATAAACACTTTGTTTCAAAATATCTAAAGCTGTGTATAGATCGTTCTCATCAGAAACTAGTAACCCGTCAATTGAATTTCCTTCATGGTCTAGAACAACCATTTCGACGc                                                                                               >  2:760113/1‑141 (MQ=255)

                                              ccTGTTCATCTGTGTTATTTTCAATCCATTCATAAACACTTTGTTTCAAAATATCTAAAGCTGTGTATAGATCGTTCTCATCAGAAACTAGTAACCCGTCAATTGAATTTCCTTCATGGTCTAGAACAACCATTTCGACGc                                                                                               >  1:663286/1‑141 (MQ=255)

                                                  ttCATCTGTGTTATTTTCAATCCATTCATAAACACTTTGTTTCAAAATATCTAAAGCTGTGTATAGATCGTTCTCATCAGAAACTAGTAACCCGTCAATTGAATTTCCTTCATGGTCTAGAACAACCATTTCGACGCtata                                                                                           >  2:923004/1‑141 (MQ=255)

                                                  ttCATCTGTGTTATTTTCAATCCATTCATAAACACTTTGTTTCAAAATATCTAAAGCTGTGTATAGATCGTTCTCATCAGAAACTAGTAACCCGTCAATTGAATTTCCTTCATGGTCTAGAACAACCATTTCGACGCtata                                                                                           >  1:644608/1‑141 (MQ=255)

                                                  ttCATCTGTGTTATTTTCAATCCATTCATAAACACTTTGTTTCAAAATATCTAAAGCTGTGTATAGATCGTTCTCATCAGAAACTAGTAACCCGTCAATTGAATTTCCTTCATGGTCTAGAACAACCATTTCGACGCtata                                                                                           >  2:538008/1‑141 (MQ=255)

                                                   tCATCTGTGTTATTTTCAATCCATTCATAAATACTTTGTTTCAAAATATCTAAAGCTGTGTATAGATCGTTCTCGTCAGAAACTAGTAGCCCGTCAATTGAATTTCCTTCATGATCTAAAACGACTATTTCGACACtatat                                                                                          >  1:734959/1‑140 (MQ=255)

                                                   tCATCTGTGTTATTTTCAATCCATTCATAAATACTTTGTTTCAAAATATCTAAAGCTGTGTATAGATCGTTCTCGTCAGAAACTAGTAGCCCGTCAATTGAATTTCCTTCATGATCTAAAACGACTATTTCGACACtatat                                                                                          >  1:526314/1‑140 (MQ=255)

                                                   tCATCTGTGTTATTTTCAATCCATTCATAAACACTTTGTTTCAAAATATCTAAAGCTGTGTATAGATCGTTCTCATCAGAAACTAGTAACCCGTCAATTGAATTTCCTTCATGGTCTAGAACAACCATTTCGACGCTATAc                                                                                          <  2:1003342/141‑1 (MQ=255)

                                                       cTGTGTTCTTTTCCATCCCTTCATCATTACTTTGTTCCAATCTATCTAAATCTGTGTATAGATCGTTCTCATCAGAAACTAGTAACCCGTCAATTGAATTTCCTTCATGGTCTAGAACAACCATTTCGACGCTATACGCTc                                                                                      <  3:199071/141‑1 (MQ=255)

                                                       cTGTGTTATTTTCAATCCATTCATAAACACTTTGTTTCAAAATATCTAAAGCTGTGTATAGATCGTTCTCATCAGAAACTAGTAACCCGTCAATTGAATTTCCTTCATGGTCTAGAACAACCATTTCGACGCTATACGCTc                                                                                      <  1:429292/141‑1 (MQ=255)

                                                       cTGTGTTATTTTCAATCCAATCATAAACACTTTGTTTCAAAATATCTAAAGCTGTGTATAGATCGTTCTCATCAGAAACTAGTAACCCGTCAATTGAATTTCCTTCATGGTCTAGAACAACCATTTCGACGCTATACGCTc                                                                                      <  3:179267/141‑1 (MQ=255)

                                                         gtgtTATTTTCAATCCATTCATAAATACTTTGTTTCAAAATATCTAAAGCTGTGTATAGATCGTTCTCGTCAGAAACTAGTAGCCCGTCAATTGAATTTCCTTCATGATCTAAAACGACTATTTCGACACTATATGCTCgc                                                                                    <  1:455632/141‑2 (MQ=255)

                                                         gtgtTATTTTCAATCCATTCATAAACACTTTGTTTCAAAATATCTAAAGCTGTGTATAGATCGTTCTCATCAGAAACTAGTAACCCGTCAATTGAATTTCCTTCATGGTCTAGAACAACCATTTCGACGCTATACGCTCGt                                                                                    >  2:635553/1‑141 (MQ=255)

                                                            ttATTTTCAATCCATTCATAAACACTTTGTTTCAAAATATCTAAAGCTGTGTATAGATCGTTCTCATCAGAAACTAGTAACCCGTCAATTGAATTTCCTTCATGGTCTAGAACAACCATTTCGACGCTATACGCTCGTTTc                                                                                 <  2:629719/141‑1 (MQ=255)

                                                              aTTTTCAATCCATTCATAAACACTTTGTTTCAAAATATCTAAAGCTGTGTATAGATCGTTCTCATCAGAAACTAGTAACCCGTCAATTGAATTTCCTTCATGGTCTAGAACAACCATTTCGACGCTATACGCTCGTTTCtt                                                                               <  1:299598/141‑1 (MQ=255)

                                                               ttttCAATCCATTCATAAACACTTTGTTTCAAAATATCTAAAGCTGTGTATAGATCGTTCTCATCAGAAACTAGTAACCCGTCAATTGAATTTCCTTCATGGTCTAGAACAACCATTTCGACGCTATACGCTCGTTTCTTa                                                                              <  4:322472/141‑1 (MQ=255)

                                                                   cAATCCATTCATAAACACTTTGTTTCAAAATATCTAAAGCTGTGTATAGATCGTTCTCATCAGAAACTAGTAACCCGTCAATTGAATTTCCTTCATGGTCTAGAACAACCATTTCGACGCTATACGCTCGTTTCTTAActc                                                                          >  1:372840/1‑141 (MQ=255)

                                                                   cAATCCATTCATAAACACTTTGTTTCAAAATATCTAAAGCTGTGTATAGATCGTTCTCATCAGAAACTAGTAACCCGTCAATTGAATTTCCTTCATGGTCTAGAACAACCATTTCGACGCTATACGCTCGTTTCTTAActc                                                                          <  1:309741/141‑1 (MQ=255)

                                                                    aaTCCATTCATACACACTTTGTTTCAAAATATCTAAAGCTGTGTATAGATCGTTCTCATCAGAAACTAGTAACCCGTCAATTGAATTTCCTTCATGGTCTAGAACAACCATTTCGACGCTATACGCTCGTTTCTTAActct                                                                         <  3:321195/141‑1 (MQ=255)

                                                                     aTCCATTCATAAACACTTTGTTTCAAAATATCTAAAGCTGTGTATAGATCGTTCTCATCAGAAACTAGTAACCCGTCAATTGAATTTCCTTCATGGTCTAGAACAACCATTTCGACGCTATACGCTCGTTTCTTAACTCtt                                                                        <  1:846708/141‑1 (MQ=255)

                                                                      tCCATTCATAAACACTTTGTTTCAAAATATCTAAAGCTGTGTATAGATCGTTCTCATCAGAAACTAGTAACCCGTCAATTGAATTTCCTTCATGGTCTAGAACAACCATTTCGACGCTATACGCTCGTTTCTTAACTCTTa                                                                       >  1:475000/1‑141 (MQ=255)

                                                                      tCCATTCATAAACACTTTGTTTCAAAATATCTAAAGCTGTGTATAGATCGTTCTCATCAGAAACTAGTAACCCGTCAATTGAATTTCCTTCATGGTCTAGAACAACCATTTCGACGCTATACGCTCGTTTCTTAACTCTTa                                                                       >  1:992887/1‑141 (MQ=255)

                                                                      tCCATTCATAAACACTTTGTTTCAAAATATCTAAAGCTGTGTATAGATCGTTCTCATCAGAAACTAGTAACCCGTCAATTGAATTTCCTTCATGGTCTAGAACAACCATTTCGACGCTATACGCTCGTTTCTTAACTCTTa                                                                       >  4:131355/1‑141 (MQ=255)

                                                                          ttCATAAACACTTTGTTTCAAAATATCTAAAGCTGTGTATAGATCGTTCTCATCAGAAACTAGTAACCCGTCAATTGAATTTCCTTCATGGTCTAGAACAACCATTTCGACGCTATACGCTCGTTTCTTAACTCTTAATTg                                                                   <  1:641314/141‑1 (MQ=255)

                                                                              tAAACACTTTGTTTCAAAATATCTAAAGCTGTGTATAGATCGTTCTCATCAGAAACTAGTAACCCGTCAATTGAATTTCCTTCATGGTCTAGAACAACCATTTCGACGCTATACGCTCGTTTCTTAACTCTTAATTG‑aaaa                                                              >  2:1046637/1‑141 (MQ=255)

                                                                                aaCACTTTGTTTCAAAATATCTAAAGCTGTGTATAGATCGTTCTCATCAGAAACTAGTAACCCGTCAATTGAATTTCCTTCATGGTCTAGAACAACCATTTCGACGCTATACGCTCGTTTCTTAACTCTTAATTG‑AAAATc                                                            <  2:730749/141‑1 (MQ=255)

                                                                                  cacTTTGTTTCAAAATATCTAAAGCTGTGTATAGATCGGTCTCATCAGAAACTAGTAACCCGTCAATTGAATTTCCTTCATGGTCTAGAACAACCATTTCGACGCTATACGCTCGTTTCTTAACTCTTAATTG‑AAAATCaa                                                          >  1:432525/1‑141 (MQ=255)

                                                                                   acTTTGTTTCAAAATATCTAAAGCTGTGTATAGATCGTTCTCATCAGAAACTAGTAACCCGTCAATTGAATTTCCTTCATGGTCTAGAACAACCATTTCGACGCTATACGCTCGTTTCTTAACTCTTAATTG‑AAAATCaaa                                                         >  2:188678/1‑141 (MQ=255)

                                                                                      ttGTTTCAAAATATCTAAAGCTGTGTATAGATCGTTCTCATCAGAAACTAGTAACCCGTCAATTGAATTTCCTTCATGGTCTAGAACAACCATTTCGACGCTATACGCTCGTTTCTTAACTCTTAATTG‑AAAATCAAAGcc                                                      >  1:1050983/1‑141 (MQ=255)

                                                                                      ttGTTTCAAAATATCTAAAGCTGTGTATAGATCGTTCTCATCAGAAACTAGTAACCCGTCAATTGAATTTCCTTCATGGTCTAGAACAACCATTTCGACGCTATACGCTCGTTTCTTAACTCTTAATTG‑AAAATCAAAGcc                                                      >  1:219361/1‑141 (MQ=255)

                                                                                      ttGTTTCAAAATATCTAAAGCTGTGTATAGATCGTTCTCATCAGAAACTAGTAACCCGTCAATTGAATTTCCTTCATGGTCTAGAACAACCATTTCGACGCTATACGCTCGTTTCTTAACTCTTAATTG‑AAAATCAAAGcc                                                      >  1:940092/1‑141 (MQ=255)

                                                                                      ttGTTTCAAAATATCTAAAGCTGTGTATAGATCGTTCTCATCAGAAACTAGTAACCCGTCAATTGAATTTCCTTCATGGTCTAGAACAACCATTTCGACGCTATACGCTCGTTTCTTAACTCTTAATTG‑AAAATCAAAGcc                                                      >  2:1054135/1‑141 (MQ=255)

                                                                                        gTTTCAAAATATCTAAAGCTGTGTATAGATCGTTCTCATCAGAAACTAGTAACCCGTCAATTGAATTTCCTTCATGGTCTAGAACAACCATTTCGACGCTATACGCTCGTTTCTTAACTCTTAATTG‑AAAATCAAAGCCAt                                                    >  3:271611/1‑141 (MQ=255)

                                                                                        gTTTCAAAATATCTAAAGCTGTGTATAGATCGTTCTCATCAGAAACTAGTAACCCGTCAATTGAATTTCCTTCATGGTCTAGAACAACCATTTCGACGCTATACGCTCGTTTCTTAACTCTTAATTG‑AAAATCAAAGCCAt                                                    <  1:929619/141‑1 (MQ=255)

                                                                                              aaaTATCTAAAGCTGTGTATAGATCGTTCTCATCAGAAACTAGTAACCCGTCAATTGAATTTCCTTCATGGTCTAGAACAACCATTTCGACGCTATACGCTCGTTTCTTAACTCTTAATTGAAAAATCAAAGCCATCTACa                                               >  1:897473/1‑141 (MQ=255)

                                                                                                atatCTAAAGCTGTGTATAGATCGTTCTCATCAGAAACTAGTAACCCGTCAATTGAATTTCCTTCATGGTCTAGAACAACCATTTCGACGCTATACGCTCGTTTCTTAACTCTTAATTG‑AAAATCAAAGCCATCTACatta                                            <  2:669992/141‑1 (MQ=255)

                                                                                                atatCTAAAGCTGTGTATAGATCGTTCTCATCAGAAACTAGTAACCCGTCAATTGAATTTCCTTCATGGTCTAGAACAACCATTTCGACGCTATACGCTCGTTTCTTAACTCTTAATTG‑AAAATCAAAGCCATCTACatta                                            <  1:460266/141‑1 (MQ=255)

                                                                                                 tatCTAAAGCTGTGTATAGATCGTTCTCATCAGAAACTAGTAACCCGTCAATTGAATTTCCTTCATGGTCTAGAACAACCATTTCGACGCTATACGCTCGTTTCTTAACTCTTAATTG‑AAAATCAAAGCCATCTACattaa                                           <  1:797050/141‑1 (MQ=255)

                                                                                                   tCTAAAGCTGTGTATAGATCGTTCTCATCAGAAACTAGTAACCCGTCAATTGAATTTCCTTCATGGTCTAGAACAACCATTTCGACGCTATACGCTCGTTTCTTAACTCTTAATTG‑AAAATCAAAGCCATCTACattaatt                                         <  1:631767/141‑1 (MQ=255)

                                                                                                   tCTAAAGCTGTGTATAGATCGTTCTCATCAGAAACTAGTAACCCGTCAATTGAATTTCCTTCATGGTCTAGAACAACCATTTCGACGCTATACGCTCGTTTCTTAACTCTTAATTG‑AAAATCAAAGCCATCTACattaatt                                         <  2:46529/141‑1 (MQ=255)

                                                                                                    cTAAAGCTGTGTATAGATCGTTCTCATCAGAAACTAGTAACCCGTCAATTGAATTTCCTTCATGGTCTAGAACAACCATTTCGACGCTATACGCTCGTTTCTTAACTCTTAATTG‑AAAATCAAAGCCATCTACattaatta                                        <  4:228471/141‑1 (MQ=255)

                                                                                                    cTAAAGCTGTGTATAGATCGTTCTCATCAGAAACTAGTAACCCGTCAATTGAATTTCCTTCATGGTCTAGAACAACCATTTCGACGCTATACGCTCGTTTCTTAACTCTTAATTG‑AAAATCAAAGCCATCTACattaatta                                        >  4:190496/1‑141 (MQ=255)

                                                                                                       aaGCTGTGTATAGATCGTTCTCATCAGAAACTAGTAACCCGTCAATTGAATTTCCTTCATGGTCTAGAACAACCATTTCGACGCTATACGCTCGTTTCTTAACTCTTAATTG‑AAAATCAAAGCCATCTACATTAATTAttt                                     >  2:1058932/1‑141 (MQ=255)

                                                                                                        aGCTGTGTATAGATCGTTCTCATCAGAAACTAGTAACCCGTCAATTGAATTTCCTTCATGGTCTAGAACAACCATTTCGACGCTATACGCTCGTTTCTTAACTCTTAATTG‑AAAATCAAAGCCATCTACATTAATTAtttt                                    <  1:576778/141‑1 (MQ=255)

                                                                                                         gCTGTGTATAGATCGTTCTCATCAGAAACTAGTAACCCGTCAATTGAATTTCCTTCATGGTCTAGAACAACCATTTCGACGCTATACGCTCGTTTCTTAACTCTTAATTG‑AAAATCAAAGCCATCTACATTAATTAttttt                                   >  1:712845/1‑141 (MQ=255)

                                                                                                            gtgtATAGATCGTTCTCATCAGAAACTAGTAACCCGTCAATTGAATTTCCTTCATGGTCTAGAACAACCATTTCGACGCTATACGCTCGTTTCTTAACTCTTAATTG‑AAAATCAAAGCCATCTACATTAATTATTTTTCGa                                >  2:409349/1‑141 (MQ=255)

                                                                                                            gtgtATAGATCGTTCTCATCAGAAACTAGTAACCCGTCAATTGAATTTCCTTCATGGTCTAGAACAACCATTTCGACGCTATACGCTCGTTTCTTAACTCTTAATTG‑AAAATCAAAGCCATCTACATTAATTATTTTTCGa                                >  4:128022/1‑141 (MQ=255)

                                                                                                               tataGATCGTTCTCATCAGAAACTAGTAACCCGTCAATTGAATTTCCTTCATGGTCTAGAACAACCATTTCGACGCTATACGCTCGTTTCTTAACTCTTAATTG‑AAAATCAAAGCCATCTACATTAATTATTTTTCGACAt                             <  1:439432/141‑1 (MQ=255)

                                                                                                                   gATCGTTCTCATCAGAAACTAGTAACCCGTCAATTGAATTTCCTTCATGGTCTAGAACAACCATTTCGACGCTATACGCTCGTTTCTTAACTCTTAATTG‑AAAATCAAAGCCATCTACATTAATTATTTTTCGACATACGt                         <  4:264319/141‑1 (MQ=255)

                                                                                                                       gTTCTCATCAGAAACTAGTAACCCGTCAATTGAATTTCCTTCATGGTCTAGAACAACCATTTCGACGCTATACGCTCGTTTCTTAACTCTTAATTG‑AAAATCAAAGCCATCTACATTAATTATTTTTCGACATACGTCAcc                     <  1:617052/141‑1 (MQ=255)

                                                                                                                         tctcATCAGAAACTAGTAACCCGTCAATTGAATTTCCTTCATGGTCTAGAACAACCTTTTCGACGCTATACGCTCGTTTCTTAACTCTTAATTG‑AAAATCACAGTCTTTTACATTAATTATTTTTCGACATACGGCACCCg                   <  4:174304/141‑1 (MQ=255)

                                                                                                                         tctcATCAGAAACTAGTAACCCGTCAATTGAATTTCCTTCATGGTCTAGAACAACCATTTCGACGCTATACGCTCGTTTCTTAACTCTTAATTG‑AAAATCAAAGCCATCTACATTAATTATTTTTCGACATACGTCACCCg                   >  2:883004/1‑141 (MQ=255)

                                                                                                                         tctcATCAGAAACTAGTAACCCGTCAATTGAATTTCCTTCATGGTCTAGAACAACCATTTCGACGCTATACGCTCGTTTCTTAACTCTTAATTG‑AAAATCAAAGCCATCTACATTAATTATTTTTCGACATACGTCACCCg                   >  2:625438/1‑141 (MQ=255)

                                                                                                                         tctcATCAGAAACTAGTAACCCGTCAATTGAATTTCCTTCATGGTCTAGAACAACCATTTCGACGCTATACGCTCGTTTCTTAACTCTTAATTG‑AAAATCAAAGCCATCTACATTAATTATTTTTCGACATACGTCACCCg                   >  2:927714/1‑141 (MQ=255)

                                                                                                                          ctcATCAGAAACTAGTAACCCGTCAATTGAATTTCCTTCATGGTCTAGAACAACCATTTCGACGCTATACGCTCGTTTCTTAACTCTTAATTG‑AAAATCAAAGCCATCTACATTAATTATTTTTCGACATACGTCACCCGt                  >  2:723508/1‑141 (MQ=255)

                                                                                                                          ctcATCAGAAACTAGTAACCCGTCAATTGAATTTCCTTCATGGTCTAGAACAACCATTTCGACGCTATACGCTCGTTTCTTAACTCTTAATTG‑AAAATCAAAGCCATCTACATTAATTATTTTTCGACATACGTCACCCGt                  >  2:723512/1‑141 (MQ=255)

                                                                                                                           tcgtcaGAAACTAGTAGCCCGTCAATTGAATTTCCTTCATGATCTAAAACGACTATTTCGACACTATATGCTCGCTTCTTAACTCTTAATTG‑AAAATCAAAGCCATCTACATTAAATATTTTTCGACATACGTCACCCGtt                 >  1:22892/4‑141 (MQ=255)

                                                                                                                            catcaGAAACTAGTAACCCGTCAATTGAATTTCCTTCATGGTCTAGAACAACCATTTCGACGCTATACGCTCGTTTCTTAACTCTTAATTG‑AAAATCAAAGCCATCTACATTAATTATTTTTCGACATACGTCACCCGttt                >  4:421057/1‑141 (MQ=255)

                                                                                                                             atcaGAAACTAGTAACCCGTCAATTGAATTTCCTTCATGGTCTAGAACAACCATTTCGACGCTATACGCTCGTTTCTTAACTCTTAATTG‑AAAATCAAAGCCATCTACATTAATTATTTTTCGACATACGTCACCCGtttt               >  1:29322/1‑141 (MQ=255)

                                                                                                                               caGAAACTAGTAACCCGTCAATTGAATTTCCTTCATGGTCTAGAACAACCATTTCGACGCTATACGCTCGTTTCTTAACTCTTAATTG‑AAAATCAAAGCCATCTACATTAATTATTTTTCGACATACGTCACCCGTTTTGt             <  1:34641/141‑1 (MQ=255)

                                                                                                                                      tagtagCCCGTCAATTGAATTTCCTTCATGATCTAAAACGACTATTTCGACACTATATGCTCGCTTCTTAACTCTTAATTG‑AAAATCAAAGCCATCTACATTAAATATTTTTCGACATACGTCACCCGTTTTGTAATACAt      >  1:715115/1‑141 (MQ=255)

                                                                                                                                      tAGTAACCCGTCAATTGAATTTCCTTCATGGTCTAGAACAACCATTTCGACGCTATACGCTCGTTTCTTAACTCTTAATTG‑AAAATCAAAGCCATCTACATTAATTATTTTTCGACATACGTCAcc                     <  1:1074934/126‑1 (MQ=255)

                                                                                                                                       aGTAACCCGTCAATTGAATTTCCTTCATGGTCTAGAACAACCATTTCGACGCTATACGCTCGTTTCTTAACTCTTAATTG‑AAAATCAAAGCCATCTACATTAATTATTTTTCGACATACGTCACCCGTTTTGTAATACAtt     >  3:57508/1‑141 (MQ=255)

                                                                                                                                       aGTAACCCGTCAATTGAATTTCCTTCATGGTCTAGAACAACCATTTCGACGCTATACGCTCGTTTCTTAACTCTTAATTG‑AAAATCAAAGCCATCTACATTAATTATTTTTCGACATACGTCACCCGTTTTGTAATACAtt     >  2:993461/1‑141 (MQ=255)

                                                                                                                                       aGTAACCCGTCAATTGAATTTCCTTCATGGTCTAGAACAACCATTTCGACGCTATACGCTCGTTTCTTAACTCTTAATTG‑AAAATCAAAGCCATCTACATTAATTATTTTTCGACATACGTCACCCGTTTTGTAATACAtt     >  1:1014721/1‑141 (MQ=255)

                                                                                                                                        gTAACCCGTCAATTGAATTTCCTTCATGGTCTAGAACAACCATTTCGACGCTATACGCTCGTTTCTTAACTCTTAATTG‑AAAATCAAAGCCATCTACATTAATTATTTTTCGACATACGTCACCCGTTTTGTAATACATTg    >  2:176986/1‑141 (MQ=255)

                                                                                                                                          aaCCCGTCAATTGAATTTCCTTCATGGTCTAGAACAACCATTTCGACGCTATACGCTCGTTTCTTAACTCTTAATTG‑AAAATCAAAGCCATCTACATTAATTATTTTTCGACATACGTCACCCGTTTTGTAATACATTGtt  >  1:26639/1‑141 (MQ=255)

TCTCATGCTTATACCTACCATTTCATGACTAAGTTCATTAGTTTGTCCTGTTCATCTGTGTTATTTTCAATCCATTCATAAACACTTTGTTTCAAAATATCTAAAGCTGTGTATAGATCGTTCTCATCAGAAACTAGTAACCCGTCAATTGAATTTCCTTCATGGTCTAGAACAACCATTTCGACGCTATACGCTCGTTTCTTAACTCTTAATTG‑AAAATCAAAGCCATCTACATTAATTATTTTTCGACATACGTCACCCGTTTTGTAATACATTGTT  >  USA300TCH1516_ALE/2120500‑2120778

Alignment Legend
Aligned base mismatch/match (shaded by quality score): `ATCG/ATCG` `< 3 ≤` `ATCG/ATCG` `< 22 ≤` `ATCG/ATCG` `< 37 ≤` `ATCG/ATCG` `< 41 ≤` `ATCG/ATCG`
Unaligned base: `atcg` Masked matching base: `atcg` Alignment gap: `‑` Deleted base: `‑`