BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	USA300TCH1516_ALE	111,516	C→A	30.1%	Q196K (CAA→AAA)	USA300HOU_RS00520 →	putative lipoprotein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	USA300TCH1516_ALE	111,516	0	C	A	30.1%	766.0 / inf	523	Q196K (CAA→AAA)	USA300HOU_RS00520	putative lipoprotein
Reads supporting (aligned to +/- strand): ref base C (187/174); new base A (59/98); total (247/276)
Fisher's exact test for biased strand distribution p-value = 3.03e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.37e-01

Only 100 of 530 total aligned reads displayed.

CTGAGGAAGTTAAAGAAGAAAACATAAAAAAAGAAATCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAAT‑TACGTAAAAGATATGACATACCAAC‑AAATAAAGCGCCGAAGCTGTTGTTGAAAGGTACAGGGAATTTAAAAGGTTCATCAGTTGGATATAAAAAAATTG  >  USA300TCH1516_ALE/111376‑111654

cTGAGGAAGTTAAAGAAGAAAACATAAAAAAAGAAATCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTc                                                                                                                                              <  1:3152026/141‑1 (MQ=255)

  gAGGAAGTTAAAGAAGAAAACATAAAAAAAGAAATCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTaaa                                                                                                                                            >  1:1805976/1‑141 (MQ=255)

    ggAAGTTAAAGAAGACAACAGAAAAAAAGAAATCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAata                                                                                                                                          <  1:1651421/141‑1 (MQ=255)

     gAAGTTAAAGAAGAAAACATAAAAGAAGAAATCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAatat                                                                                                                                         <  2:3542165/141‑1 (MQ=255)

     gAAGTTAAAGAAGAAAACATAAAAAAAGAAATCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGACCACTTATTCGGCTAAatat                                                                                                                                         <  1:2821443/141‑1 (MQ=21)

       aGTTAAAGAAGAAAACATAAAAAAAGAAATCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCa                                                                                                                                       <  2:26243/141‑1 (MQ=255)

        gTTAAAGAAGAAAACATAAAAAAAGAAATCGAAAACTTTAAGTTTTTTGCGCAATATCGCCGCCGTGAAGATGTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCaa                                                                                                                                      <  1:1027484/141‑1 (MQ=21)

          tAAAGAAGAAAACATAAAAAAAGAAATCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATACAGAGGTGCCGAGTTATTCGGCTAAATATCAAtt                                                                                                                                    >  1:1578412/1‑141 (MQ=255)

          aaaaGAAGAAAAAAGAAAAAAAGAAAACGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAAtt                                                                                                                                    <  1:3449086/140‑1 (MQ=17)

           aaaGAAGAAAACATAAAAAAAGAAATCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAAAGGTGCAGAGTTATTCGGCTAAATATCAATTa                                                                                                                                   >  2:544944/1‑141 (MQ=255)

            aagaagAAAACATAAAAAAAGAAATCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTaa                                                                                                                                  >  1:1419943/1‑141 (MQ=255)

            aagaagAAAACATAAAAAAAGAAATCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTaa                                                                                                                                  <  2:2187314/141‑1 (MQ=255)

             agaagaAAACATAAAAAAAGAAATCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAAc                                                                                                                                 >  2:3333212/1‑141 (MQ=255)

             agaagaAAACAGAAAAAAAGAAATCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAAc                                                                                                                                 <  1:3348823/141‑1 (MQ=255)

              gaagaaAACATAAAAAAAGAAATCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACt                                                                                                                                >  1:871633/1‑141 (MQ=255)

              gaagaaAACATAAAAAAAGAAATCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACt                                                                                                                                >  1:936744/1‑141 (MQ=255)

              gaagaaAACATAAAAAAAGAAATCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTCTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACt                                                                                                                                <  1:2880039/141‑1 (MQ=255)

               aagagAACATAAAAAAAGAAATCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTa                                                                                                                               <  1:2454666/141‑1 (MQ=255)

               aagaaAAAAAAAAAAAAGAAATCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTa                                                                                                                               <  2:375563/141‑1 (MQ=255)

                  aaaaCATAAAAAAAGAAATCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAatg                                                                                                                            >  2:1875291/1‑141 (MQ=255)

                   aaaCATAAAAAAAGAAATCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAatga                                                                                                                           <  2:3876875/141‑1 (MQ=255)

                    aaCATAAAAAAAGAAATCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAatgat                                                                                                                          <  2:1335870/141‑1 (MQ=255)

                     aCATAAAAAAAGAAATCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAatgatg                                                                                                                         >  1:2721225/1‑141 (MQ=255)

                     aCATAAAAAAAGAAATCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAatgatg                                                                                                                         >  2:2283701/1‑141 (MQ=255)

                      cATAAAAAAAGAAATCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAAtgatga                                                                                                                        >  2:4373152/1‑141 (MQ=255)

                      cATAAAAAAAGAAATCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAAtgatga                                                                                                                        >  1:1859837/1‑141 (MQ=255)

                       aaaaaaaaaGAAATCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGGGTTACACAGCTCAAGATCAATTAACTAATgatgat                                                                                                                       <  1:3049026/139‑1 (MQ=255)

                       aTAAAAAAAGAAATCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATgatgat                                                                                                                       >  1:3117254/1‑141 (MQ=255)

                          aaaaaaGAAATCGCAAACTTTAAGTTTTTTGCTCAATATAGCCATTCCAAACATTTAATGATTTATAACGCTGGAGATATCTCATATAATCCAGAGGTGCCGAGTTATACAGCTCAATATCAATTAACTAATGCTGATtat                                                                                                                    >  2:2755436/1‑141 (MQ=17)

                          aaaaaaGAAATCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATtat                                                                                                                    <  1:2539289/141‑1 (MQ=255)

                           atataGAAATCTAAAACTTTAAGATTTTTGCGCAATATAGCAACTCTAAAGATTTAATGAATTATAAAGATGGAGTTTTAGCAGAAAATCACGAGACGATGGGTTATTCGGCTAAATATCAATTAACTAATGATGATtata                                                                                                                   <  2:1598649/137‑1 (MQ=12)

                           aaaaaGAAATCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATtata                                                                                                                   <  2:2136036/141‑1 (MQ=255)

                           aaaaaGAAATCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGAGATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATtata                                                                                                                   <  2:1102091/141‑1 (MQ=255)

                              aaGAAATCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATg                                                                                                                >  1:2647110/1‑141 (MQ=255)

                              aaGAAATCGAAAACTTTAAGTGTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATg                                                                                                                <  2:2919056/141‑1 (MQ=255)

                                gAAATCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTa                                                                                                              >  1:1182137/1‑141 (MQ=255)

                                 aaaTCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTaa                                                                                                             >  1:1327485/1‑141 (MQ=255)

                                  aaTCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATTCAGAGGTGCCGAGTTATTCGGCTCAAAACCAATTAATTAATGATGATTATGTTGTaaa                                                                                                            >  1:3407916/1‑141 (MQ=18)

                                  aaTCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTaaa                                                                                                            <  2:469463/141‑1 (MQ=255)

                                   atagAAAACTTTAAGTATTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTaaaa                                                                                                           <  2:1725437/138‑1 (MQ=255)

                                     cGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACa                                                                                                         <  1:977536/141‑1 (MQ=255)

                                      gAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACaa                                                                                                        >  2:3690366/1‑141 (MQ=255)

                                       agaaCTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAAt                                                                                                       <  2:2059578/139‑1 (MQ=255)

                                       aaaaCTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATCACACAGAGATTCGGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAAt                                                                                                       <  1:1693328/141‑1 (MQ=21)

                                         gtcTTTAAGTTTTTTGTGGAATATAGCGATTTTAGAGATTTAATGAGTTATAAAGATGGAGATATATCATATAATCCAGAGGTGCTGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAAT‑Ta                                                                                                    <  2:4064148/139‑1 (MQ=17)

                                         aaCTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAAT‑Ta                                                                                                    <  1:3912669/141‑1 (MQ=255)

                                            tttAAGTTTTTTGCGCAACATAGATACTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAAT‑TACgc                                                                                                 <  2:1766896/141‑2 (MQ=14)

                                            tttAAGTTTTTTGCGAAATATAGCAATTTTAAAGATTTAAGGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAAT‑TACgc                                                                                                 <  2:922689/141‑2 (MQ=17)

                                              tAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAAT‑TACGTaa                                                                                               <  1:2950174/141‑1 (MQ=255)

                                              tAAGTTTTTGGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAAT‑TACgcaa                                                                                               <  2:2343954/141‑4 (MQ=255)

                                                   ttttgGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAAT‑TACGTAAAAGat                                                                                          <  1:4279002/141‑1 (MQ=255)

                                                    ttttGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCACCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAAT‑TACGTAAAAGata                                                                                         >  1:1230019/1‑141 (MQ=255)

                                                     tttGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAAT‑TACGTAAAAGatat                                                                                        >  1:410966/1‑141 (MQ=255)

                                                      ttGCGCAATATAACAATTTTAAATATTTAATGAATTATAAAGATAGATTTAGAGCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAAT‑TACGCAAAAGATATg                                                                                       <  2:4274433/141‑1 (MQ=14)

                                                         cgcAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAAT‑TACGTAAAAGATATGACa                                                                                    >  1:3568780/1‑141 (MQ=255)

                                                          gcAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAAT‑TACGCAAAAgatatgatat                                                                                   >  1:1467759/1‑141 (MQ=14)

                                                              aatatCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAAT‑TACGTAAAAGATATGACATACca                                                                               <  2:3909247/140‑1 (MQ=255)

                                                                  gCAATTGTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAAT‑TACGCAAAAGATATGATATAccgacga                                                                           <  2:3030435/141‑6 (MQ=11)

                                                                   cAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAAT‑TACGTAAAAGATATGACATACCAAC‑aaa                                                                         >  2:325394/1‑141 (MQ=255)

                                                                   cAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAAT‑TACATAAAAGATATGACATACCAAC‑aaa                                                                         >  2:249743/1‑141 (MQ=255)

                                                                        ttAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAAT‑TACGTAAAAGATATGACATACCAAC‑AAATAAAg                                                                    <  2:3413389/141‑1 (MQ=255)

                                                                         tAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAAT‑TACGTAAAAGATATGACATACCAAC‑AAATAAAgc                                                                   >  1:3603031/1‑141 (MQ=255)

                                                                         tAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAAT‑TACGTAAAAGATATGACATACCAAC‑AAATAAAgc                                                                   >  2:1646457/1‑141 (MQ=255)

                                                                          aaaGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAAT‑TACGTAAAAGATATGACATACCAAC‑AAATAAAgcg                                                                  >  2:1050253/1‑141 (MQ=255)

                                                                            aGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAAT‑TACGTAAAAGATATGACATACCAAC‑AAATAAAGCGcc                                                                <  2:1393616/141‑1 (MQ=255)

                                                                             gATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAAT‑TACGTAAAAGATATGACATACCAAC‑AAATAAAGCGCCg                                                               <  1:3752101/141‑1 (MQ=255)

                                                                             gATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAAT‑TACGTAAAAGATATGACATAACAAC‑AAATAAAGCGCCg                                                               >  1:3637414/1‑141 (MQ=255)

                                                                                   aTGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAAT‑TACGTAAAAGATATGACATACCAAC‑AAATAAAGCGCCGAAGCtg                                                         <  2:772057/141‑1 (MQ=255)

                                                                                       aTTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAAT‑TACGTAAAAGATATGACATACCAAC‑AAATAAAGCGCCGAAGCtgttgt                                                     <  1:2369544/141‑1 (MQ=255)

                                                                                         tataAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAAT‑TACGTAAAAGATATGACATACCAAC‑AAATAAAGCGCCGAAGCTGttgttg                                                   <  1:325394/141‑1 (MQ=255)

                                                                                               gATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATAAATTAACTAAGTATGATTATACTGTAAAACAAT‑TACGTAAAAGATATGACATACCAAC‑AAATAAAGCGCCGAAGCTGTTGTTGAAAGGt                                             <  2:3141908/141‑1 (MQ=255)

                                                                                                    agaTATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAAT‑TACGTAAAAGATATGACATACCAAC‑AAATAAAGCGCCGAAGCTGTTGTTGAAAGGTACAgg                                        >  1:2288929/1‑141 (MQ=255)

                                                                                                         tatCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAAT‑TACGTAAAAGATATGACATACCAAC‑AAATAACGCGCCGAAGCTGTTGTTGAAAGGTACCGGCAAtt                                   >  2:592874/1‑141 (MQ=255)

                                                                                                         tatCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAAT‑TACGTAAAAGATATGACATACCAAC‑AAATAAAGCGCCGAAGCTGTTGTTGAAAGGTACAGGGAAtt                                   <  2:1182137/141‑1 (MQ=255)

                                                                                                         tatCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAAT‑TACGTAAAAGATATGACATACCAAC‑AAATAAAGCGCCGAAGCTGTTGTTGAAAGGTACAGGGAAtt                                   <  2:3950301/141‑1 (MQ=255)

                                                                                                          atcatatcatCCAGAGGTGCCGAGATATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAAT‑TACGTAAAAGATATGACATACCAAC‑AAATAAAGCGCCGAAGCTGTTGTTGAAAGGTACAGGGAAttt                                  <  1:802065/141‑1 (MQ=255)

                                                                                                             atatAATCCGCAGGTGCAGAGTTATTCAGCTCAATATCCATTAACTAATGATGATTCTAATGCACCACAGT‑TTCGTGCAAGATATGACATACCAAC‑AAATAAACCGCCGAATCCGTCGTTGAACGGTACAGGGACTTTgaa                               >  2:362761/1‑141 (MQ=9)

                                                                                                                taATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAAT‑TACGTAAAAGATATGACATACCAAC‑AAATAAAGCGCCGAAGCTGTTGTGGAAAGGTACAGGGAATTTAAAAgg                            <  1:764275/141‑1 (MQ=255)

                                                                                                                 aaTCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAAT‑TACGTAAAAGATATGACATACCAAC‑AAATAAAGCGCCGAAGCTGTTGTTGAAAGGTACAGGGAATTTAAAAGGt                           <  1:1740354/141‑1 (MQ=255)

                                                                                                                 aaTCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAAT‑TACGTAAAAGATATGACATACCAAC‑AAATAAAGCGCCGAAGCTGTTGTTGAAAGGTACAGGGAATTTAAAAGGt                           <  1:2117961/141‑1 (MQ=255)

                                                                                                                 aaTCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAAT‑TACGTAAAAGATATGACATACCAAC‑AAATAAAGCGCCGAAGCTGTTGTTGAAAGGTACAGGGAATTTAAAAGGt                           <  1:1991972/141‑1 (MQ=255)

                                                                                                                    ccAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAAT‑TACGTAAAAGATATGACATACCAAC‑AAATAAAGCGCCGAAGCTGTTGTTGAAAGGTACAGGGAATTTAAAAGGTtca                        >  1:2410512/1‑141 (MQ=255)

                                                                                                                    ccAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAAT‑TACGTAAAAGATACG‑CCTCTGAACTAGAGTAAGCGCCGAAGCTGTTGTTGAAAGGTACAGGGAATTTAAAAGGTtca                        >  1:1899931/1‑141 (MQ=17)

                                                                                                                    ccAGAGGAGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAAT‑TACGTAAAAGATATGACATACCAAC‑AAATAAAGCGCCGAAGCTGTTGTTGAAAGGTACAGGGAATTTAAAAGGTtca                        >  2:871952/1‑141 (MQ=255)

                                                                                                                     cAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAAT‑TACGTAAAAGAGATGACTTACCAAC‑AAATAAAGCGCCGAAGCTGTTGTTGAAAGGTACAGGGAATTTAAAAGGTtcat                       >  1:396068/1‑141 (MQ=255)

                                                                                                                      agagGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAAT‑TACGTAAAAGATATGACATACCAAC‑AAATAAAGCGCCGAAGCTGTTGTTGAAAGGTACAGGGAATTTAAAAGGTtcatc                      >  1:481610/1‑141 (MQ=255)

                                                                                                                      agagGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAAT‑TACGTAAAAGATATGACATACCAAC‑AAATAAAGCGCCGAAGCTGTTGTTGAAAGGTACAGGGAATTTAAAAGGTtcatc                      >  2:4172363/1‑141 (MQ=255)

                                                                                                                         ggTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAAT‑TACGTAAAAGATATGACATACCAAC‑AAATAAAGCGCCGAGGCTGTTGTTGAAAGGTACAGGGAATTTAAAAGGTTCATCAGt                   >  2:3792982/1‑141 (MQ=255)

                                                                                                                         ggTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAAT‑TACGTAAAAGATATGACATACCAAC‑AAATAAAGCGCCGAAGCTGTTGTTGAAAGGTACAGGGAATTTAAAAGGTTCATCAGt                   >  1:3831102/1‑141 (MQ=255)

                                                                                                                         ggTGCCGAGTTATTCAGCTCAATATCAATTAACTAATCATTAATATAAGGGAAAACCATAAACGGTAGAG‑TATGACATACCAAC‑AAATAAAGCGCCGAAGCTGTTGTTGAAAGGTACAGGGAATTTAAAAGGTTCATCAGt                   >  1:4394421/1‑141 (MQ=17)

                                                                                                                          gtgtCGAGTTATTCAGCACAATATCAATTAACTAAAGATGATTATAATGTAAAACAAT‑TACGTAAAAGATATGGCATACCAAC‑AAATAAAGCGCCGAAGCTGTTGTTGAAAGGTACAGGGAATTTAAAAGGTTCATCAGtt                  <  2:2975786/141‑1 (MQ=255)

                                                                                                                          gTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAAT‑TACGTAAAAGATATGACATACCAAC‑AAATAAAGAGCCGAAGCAGGTGTTGAAAGGTACAGGGAATTTAAAAGGTTCATCAGtt                  >  2:2411026/1‑141 (MQ=255)

                                                                                                                           tGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAAG‑TACGTAAAAGATATGACATACCAAC‑AAATAAAGCGCCGAAGCTGTTGTTGAAAGGTACAGGGAATTTAAAAGGTTCATCAGTTg                 >  1:2241130/1‑141 (MQ=255)

                                                                                                                            gCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAAT‑TACGTAAAAGATATGACATACCAAC‑AAATAAAGCGCCGAAGCTGTTGTTGAAAGGTACAGGGAATTTAAAAGGTTCATCAGTTgg                >  2:1943763/1‑141 (MQ=255)

                                                                                                                              cGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAAT‑TACGTAAAAGATATGACATACCAAC‑AAATAAAGCGCCGAAGCTGTTGTTGAAAGGTACAGGGAATTTAAAAGGTCCATCAGTTGGat              >  2:2602288/1‑141 (MQ=255)

                                                                                                                               gAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAAT‑TACGTAAAAGATATGACATACCAAC‑AAATAAAGCGCCGCAGCTGTTGCTGGCTCGTAAAGGGAATTTAAAAGGTTCTGCAGTTGGata             >  1:606430/1‑141 (MQ=14)

                                                                                                                               gAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAAT‑TACGTAAAAGATATGACATACCAAC‑AAATAAAGCGCCGAAGCTGTTGTTGAAAGGTACAGGGAATTTAAAAGGTTCATCAGTTGGata             >  1:232402/1‑141 (MQ=255)

                                                                                                                                  acaataaGCTCAATATCAATTAACGAATGATTATTATAATGTAAAACAAT‑TACGTAAAAGATATGACATACCAAC‑AAATAAAGCGCCGAAGCTGTTGTTGAAAGGTACAGGGAATTTAAAAGGTTCATCAGTTGGATATaa          <  2:1513264/135‑1 (MQ=255)

                                                                                                                                   tATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAAT‑TACGTAAAAGATATGACATACCAAC‑AAATAAAGCGCCGAAGCTGTTGTTGAAAGGTACAGGGAATTTAAAAGGTTCATCAGTTGGATATaaa         >  1:675077/1‑141 (MQ=255)

                                                                                                                                       cAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAAT‑TACGTAAAAGATATGACATACCAAC‑AAATAAAGCGCCGAAGCTGTTGTTGAAAGGTACAGGGAATTTAAAAGGTTCATCAGTTGGATATaaaaaaa     >  2:3467062/1‑141 (MQ=255)

CTGAGGAAGTTAAAGAAGAAAACATAAAAAAAGAAATCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAAT‑TACGTAAAAGATATGACATACCAAC‑AAATAAAGCGCCGAAGCTGTTGTTGAAAGGTACAGGGAATTTAAAAGGTTCATCAGTTGGATATAAAAAAATTG  >  USA300TCH1516_ALE/111376‑111654

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