BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	USA300TCH1516_ALE	111,516	C→A	31.2%	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	31.2%	264.4 / 161.9	154	Q196K (CAA→AAA)	USA300HOU_RS00520	putative lipoprotein
Reads supporting (aligned to +/- strand): ref base C (55/51); new base A (21/27); total (76/78)
Fisher's exact test for biased strand distribution p-value = 3.87e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.35e-01

Only 100 of 156 total aligned reads displayed.

TTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGA  >  USA300TCH1516_ALE/111426‑111606

   tttGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAat                                                                                         <  1:2422033/91‑1 (MQ=255)

    ttGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAata                                                                                        >  1:2185299/1‑91 (MQ=255)

    ttGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAata                                                                                        <  1:1844307/91‑1 (MQ=255)

    ttGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAata                                                                                        >  1:328553/1‑91 (MQ=255)

     tGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAatat                                                                                       <  1:3999901/91‑1 (MQ=255)

     tGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAatat                                                                                       >  1:196135/1‑91 (MQ=255)

       cgcAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCa                                                                                     >  1:136760/1‑91 (MQ=255)

        gcAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCaa                                                                                    >  1:1438819/1‑91 (MQ=255)

        gcAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCaa                                                                                    <  1:3016966/91‑1 (MQ=255)

        gcAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCaa                                                                                    >  1:3767614/1‑91 (MQ=255)

         cAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAAt                                                                                   >  1:4061979/1‑91 (MQ=255)

         cAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAAt                                                                                   >  1:4062013/1‑91 (MQ=255)

         cAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAAt                                                                                   >  1:753487/1‑91 (MQ=255)

         cAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAAt                                                                                   >  1:3564109/1‑91 (MQ=255)

         cAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAAt                                                                                   >  1:4079670/1‑91 (MQ=255)

          aaTATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAAtt                                                                                  <  1:3813086/91‑1 (MQ=255)

          aaTATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAAtt                                                                                  <  1:1387934/91‑1 (MQ=255)

            tataGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTaa                                                                                >  1:530414/1‑91 (MQ=255)

             ataGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAAc                                                                               <  1:3835465/91‑1 (MQ=255)

               aGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTa                                                                             <  1:864716/91‑1 (MQ=255)

                gCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTaa                                                                            >  1:3437450/1‑91 (MQ=255)

                gCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTaa                                                                            >  1:3828997/1‑91 (MQ=255)

                gCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTaa                                                                            >  1:3659677/1‑91 (MQ=255)

                gCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTaa                                                                            >  1:1003247/1‑91 (MQ=255)

                gCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTaa                                                                            >  1:808053/1‑91 (MQ=255)

                gCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTaa                                                                            >  1:50567/1‑91 (MQ=255)

                gCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTaa                                                                            >  1:1107794/1‑91 (MQ=255)

                  aaTTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAatg                                                                          >  1:1187215/1‑91 (MQ=255)

                   aTTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAatga                                                                         <  1:221373/91‑1 (MQ=255)

                    ttttAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAatgat                                                                        >  1:1416690/1‑91 (MQ=255)

                     tttAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAatgatg                                                                       <  1:1872436/91‑1 (MQ=255)

                      ttAAAGATTTAATGAATTATAAAGATGGATAGATATCATATAATGCAGAGGTGGCGAGTTATTCGGCTAAATATCAATTAACTAAtgatga                                                                      <  1:3317578/91‑1 (MQ=14)

                      ttAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAAtgatga                                                                      <  1:699882/91‑1 (MQ=255)

                      ttAAAGATTTAATGAATTATAAAGATGGAGAGATATCATATAATGCAGAGGTGGCGAGTTATTCGGCTAAATATCAATTAACTAAtgatga                                                                      <  1:3316862/91‑1 (MQ=17)

                       tAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATgatgat                                                                     <  1:2795506/91‑1 (MQ=255)

                        aaaGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGAtt                                                                    >  1:3948619/1‑91 (MQ=255)

                        aaaGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGAtt                                                                    >  1:769844/1‑91 (MQ=255)

                        aaaGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGAtt                                                                    >  1:1749814/1‑91 (MQ=255)

                        aaaGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGAtt                                                                    >  1:1345058/1‑91 (MQ=255)

                        aaaGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGAtt                                                                    >  1:2059773/1‑91 (MQ=255)

                        aaaGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGAtt                                                                    <  1:1434791/91‑1 (MQ=255)

                          aGATTTATTGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATtat                                                                  <  1:2077649/91‑1 (MQ=14)

                          aGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATtat                                                                  >  1:22297/1‑91 (MQ=255)

                          aGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATtat                                                                  <  1:3683271/91‑1 (MQ=255)

                          aGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATtat                                                                  <  1:2953697/91‑1 (MQ=255)

                          aGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATtat                                                                  >  1:1776870/1‑91 (MQ=255)

                            aTTTAATGAGTTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATaa                                                                <  1:3973905/91‑1 (MQ=255)

                            aTTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATaa                                                                >  1:783059/1‑91 (MQ=255)

                            aTTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATaa                                                                <  1:2738627/91‑1 (MQ=255)

                              ttAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATg                                                              <  1:2036676/91‑1 (MQ=255)

                               tAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGt                                                             <  1:1454409/91‑1 (MQ=255)

                                aaTGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTa                                                            <  1:1944641/91‑1 (MQ=255)

                                 aTGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTaa                                                           <  1:325468/91‑1 (MQ=255)

                                 aTGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTaa                                                           <  1:887665/91‑1 (MQ=255)

                                  tGAATTATAAAGATGGAGATATATCGTATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTaaa                                                          <  1:3167205/91‑1 (MQ=11)

                                  tGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTaaa                                                          <  1:52578/91‑1 (MQ=255)

                                   gAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTaaaa                                                         <  1:169228/91‑1 (MQ=255)

                                      ttATAAAGATGGAGATATATCATATAATCCAGAGGTGCTGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACaa                                                      <  1:932065/91‑1 (MQ=255)

                                          aaaGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTaa                                                  >  1:2553321/1‑90 (MQ=255)

                                          aaaGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTAc                                                  >  1:2915919/1‑91 (MQ=255)

                                           aaGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACg                                                 >  1:3069691/1‑91 (MQ=255)

                                           aaGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACg                                                 >  1:2291917/1‑91 (MQ=255)

                                            aGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGt                                                >  1:3930609/1‑91 (MQ=255)

                                             gATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTa                                               <  1:1715359/91‑1 (MQ=255)

                                                   gaTATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGat                                         >  1:1184909/1‑91 (MQ=255)

                                                   gaTATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGat                                         <  1:1479316/91‑1 (MQ=255)

                                                   gaTATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGat                                         <  1:2009874/91‑1 (MQ=255)

                                                   gaTATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGat                                         >  1:4150884/1‑91 (MQ=255)

                                                       tatCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATCTGa                                     >  1:3218041/1‑91 (MQ=255)

                                                       tatCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGa                                     <  1:675700/91‑1 (MQ=255)

                                                        atCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGAc                                    <  1:4009280/91‑1 (MQ=255)

                                                         tCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACa                                   >  1:2138709/1‑91 (MQ=255)

                                                             ataATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATAcc                               <  1:2750689/91‑1 (MQ=255)

                                                             ataATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATAcc                               >  1:2095985/1‑91 (MQ=255)

                                                             ataATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATAcc                               >  1:436238/1‑91 (MQ=255)

                                                                  ccAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACcaacaa                          <  1:2306702/91‑1 (MQ=255)

                                                                      agGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATaa                      >  1:1084043/1‑91 (MQ=255)

                                                                      agGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATaa                      <  1:58109/91‑1 (MQ=255)

                                                                      agGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATaa                      <  1:2175481/91‑1 (MQ=255)

                                                                       ggTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATaaa                     >  1:1432474/1‑91 (MQ=255)

                                                                       ggTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATaaa                     >  1:2799869/1‑91 (MQ=255)

                                                                         tGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAgc                   >  1:3976620/1‑91 (MQ=255)

                                                                          gCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAgcg                  <  1:1021124/91‑1 (MQ=255)

                                                                          gCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAgcg                  >  1:190293/1‑91 (MQ=255)

                                                                           ccGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAgcgc                 >  1:1561023/1‑91 (MQ=255)

                                                                           ccGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAgcgc                 >  1:815407/1‑91 (MQ=255)

                                                                            cGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGcc                <  1:3263717/91‑1 (MQ=255)

                                                                            cGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGcc                <  1:171154/91‑1 (MQ=255)

                                                                            cGAGTTACTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAGTAAAGCGcc                <  1:2087237/91‑1 (MQ=255)

                                                                             gAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCg               <  1:297762/91‑1 (MQ=255)

                                                                             gAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCg               >  1:2383715/1‑91 (MQ=255)

                                                                               gTTATGCAGCTCAATATCAATTAACTAATGATGATTATAGTGTAAAACAATTATGTAAAAGATATGACATACCAACAAATAAAGCGCCGaa             <  1:3853677/91‑1 (MQ=17)

                                                                                 tATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGc           >  1:3053200/1‑91 (MQ=255)

                                                                                  aTTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCt          <  1:2403929/91‑1 (MQ=255)

                                                                                   ttCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCtg         <  1:146017/91‑1 (MQ=255)

                                                                                   ttCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCtg         >  1:2596072/1‑91 (MQ=255)

                                                                                    tCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAACCAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCtgt        >  1:3596283/1‑91 (MQ=255)

                                                                                    tCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCtgt        <  1:2831904/91‑1 (MQ=255)

                                                                                    tCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCtgt        >  1:3477000/1‑91 (MQ=255)

                                                                                     cAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCtgtt       >  1:3662911/1‑91 (MQ=255)

TTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGA  >  USA300TCH1516_ALE/111426‑111606

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