BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	USA300TCH1516_ALE	111,516	C→A	26.3%	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	26.3%	1252.7 / inf	558	Q196K (CAA→AAA)	USA300HOU_RS00520	putative lipoprotein
Reads supporting (aligned to +/- strand): ref base C (155/256); new base A (46/101); total (201/357)
Fisher's exact test for biased strand distribution p-value = 1.93e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.25e-01

Only 100 of 566 total aligned reads displayed.

TTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGA  >  USA300TCH1516_ALE/111426‑111606

ttttttGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTc                                                                                            <  2:3159368/91‑1 (MQ=255)

    ttGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAata                                                                                        <  1:1557366/91‑1 (MQ=255)

     tGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAatat                                                                                       <  1:423164/91‑1 (MQ=255)

       cgcAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCa                                                                                     <  1:512897/91‑1 (MQ=255)

        gcAATATGGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCaa                                                                                    <  2:4570433/91‑1 (MQ=255)

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

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

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

        gcAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCaa                                                                                    <  2:3059175/91‑1 (MQ=255)

         ttatatAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAAt                                                                                   >  2:4983787/3‑91 (MQ=255)

         cAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATACTCCAGAGGTGCCGAGTTATTCAGCTCAATATCAAt                                                                                   <  1:4946940/91‑1 (MQ=255)

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

         cAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAAt                                                                                   <  1:533545/91‑1 (MQ=255)

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

         aaaTATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAAt                                                                                   <  2:389388/90‑1 (MQ=255)

           atatAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTa                                                                                 <  2:5894770/91‑1 (MQ=255)

              taGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACt                                                                              <  1:2643828/91‑1 (MQ=255)

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

                 cAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAat                                                                           <  2:4351984/91‑1 (MQ=255)

                  aaTTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAatg                                                                          <  1:2737811/91‑1 (MQ=255)

                   aTTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAatga                                                                         <  1:3564312/91‑1 (MQ=255)

                   aTTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAatga                                                                         <  1:5070851/91‑1 (MQ=255)

                    ttttAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAatgat                                                                        <  2:3697733/91‑1 (MQ=255)

                    ttttAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAatgat                                                                        <  2:1975369/91‑1 (MQ=255)

                    ttttAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAatgat                                                                        <  2:5456295/91‑1 (MQ=255)

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

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

                       tAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATgatgat                                                                     <  2:7001375/91‑1 (MQ=255)

                        aaaGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGAtt                                                                    <  1:5596300/91‑1 (MQ=255)

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

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

                        aaaGATTTAAGGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATCCGGCTAAATATCAATTAACTAATGATGAtt                                                                    <  2:5554306/91‑1 (MQ=14)

                         aaGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATta                                                                   >  2:261000/1‑91 (MQ=255)

                         aaGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATta                                                                   <  2:3527845/91‑1 (MQ=255)

                         aaGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATta                                                                   <  2:1865847/91‑1 (MQ=255)

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

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

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

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

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

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

                            aTTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATaa                                                                >  1:5771668/1‑91 (MQ=255)

                             tttAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAAt                                                               <  1:697020/91‑1 (MQ=255)

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

                               tAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGt                                                             >  2:124571/1‑91 (MQ=255)

                                 aggAATTATGAAGATGGGGGTAGATCATGTATTGCAGAGGTGCCGAGTTATTCAGCTCAATATCATTTATTTTATGGTTATTATAATGTat                                                           <  2:4689396/89‑2 (MQ=9)

                                 aTGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAAt                                                               >  1:5637804/1‑87 (MQ=255)

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

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

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

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

                                       tataAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAAt                                                     >  2:6381335/1‑91 (MQ=255)

                                       tataAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAAt                                                     <  1:5451210/91‑1 (MQ=255)

                                        ataAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAAtt                                                    <  1:4471079/91‑1 (MQ=255)

                                          aaaGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTAc                                                  <  2:1027163/91‑1 (MQ=255)

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

                                             gATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTa                                               >  2:5769572/1‑91 (MQ=255)

                                              aTGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTaa                                              <  1:2583404/91‑1 (MQ=255)

                                                  agaTATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGa                                          >  2:2808627/1‑91 (MQ=255)

                                                   gaTATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGat                                         <  1:2459306/91‑1 (MQ=255)

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

                                                    atatatCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGata                                        <  1:4579008/91‑1 (MQ=255)

                                                     tatatCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGatat                                       <  2:6207516/91‑1 (MQ=255)

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

                                                         tCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACa                                   <  2:6840301/91‑1 (MQ=255)

                                                          cATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATGATGTAAAACAATTACGTAAAAGATATGACAt                                  >  1:6495303/1‑91 (MQ=255)

                                                          cATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACAt                                  >  2:2291763/1‑91 (MQ=255)

                                                           atatAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATa                                 <  1:4325411/91‑1 (MQ=255)

                                                           atatAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATa                                 >  2:1367196/1‑91 (MQ=255)

                                                           atatAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATa                                 <  1:61533/91‑1 (MQ=255)

                                                            tataATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATAc                                >  1:1021854/1‑91 (MQ=255)

                                                              taATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACca                              <  2:2172597/91‑1 (MQ=255)

                                                              taATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACca                              <  1:2868819/91‑1 (MQ=255)

                                                               aaTCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACcaa                             >  1:1476953/1‑91 (MQ=255)

                                                               aaTCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATAAcaa                             >  2:1644274/1‑91 (MQ=255)

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

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

                                                                  ccAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACcaacaa                          >  1:2263177/1‑91 (MQ=255)

                                                                    agagGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAAt                        >  1:6282867/1‑91 (MQ=255)

                                                                       ggTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATaaa                     <  1:1973860/91‑1 (MQ=255)

                                                                         tGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAgc                   <  1:111496/91‑1 (MQ=255)

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

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

                                                                            cGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGCAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGcc                <  2:6766765/91‑1 (MQ=255)

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

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

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

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

                                                                               gTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGaa             <  1:818978/91‑1 (MQ=255)

                                                                               gTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGaa             <  2:5663372/91‑1 (MQ=255)

                                                                               gTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGaa             <  1:3864064/91‑1 (MQ=255)

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

                                                                                   ttCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCg               <  1:3452710/85‑1 (MQ=255)

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

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

                                                                                      aGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACATATAAAGCGCCGAAGCtgttg      <  1:6385982/91‑1 (MQ=255)

                                                                                       gCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAgc                   >  1:1232343/1‑77 (MQ=255)

                                                                                       gCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCtgttgt     >  1:6814908/1‑91 (MQ=255)

                                                                                          cAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGa  <  2:4015788/91‑1 (MQ=255)

                                                                                          cAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGa  <  1:1831131/91‑1 (MQ=255)

TTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGA  >  USA300TCH1516_ALE/111426‑111606

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