BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	USA300TCH1516_ALE	110,660	T→A	23.1%	I173I (ATT→ATA)	USA300HOU_RS00515 →	putative lipoprotein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	USA300TCH1516_ALE	110,660	0	T	A	23.1%	292.0 / 92.3	130	I173I (ATT→ATA)	USA300HOU_RS00515	putative lipoprotein
Reads supporting (aligned to +/- strand): ref base T (42/58); new base A (6/24); total (48/82)
Fisher's exact test for biased strand distribution p-value = 3.21e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.82e-01

Only 100 of 138 total aligned reads displayed.

GAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGA  >  USA300TCH1516_ALE/110571‑110749

gAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATAttt                                                                                          <  1:3101196/91‑1 (MQ=255)

gAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATAttt                                                                                          <  1:1580158/91‑1 (MQ=255)

 aaaaaaTAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTc                                                                                         <  1:3145984/91‑1 (MQ=255)

 aaaaaaTAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTc                                                                                         <  1:3145441/91‑1 (MQ=255)

   aaaaTAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCat                                                                                       <  1:257106/91‑1 (MQ=255)

   aaaaTAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCat                                                                                       >  1:92326/1‑91 (MQ=255)

    aaaTAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCata                                                                                      >  1:4088318/1‑91 (MQ=255)

    aaaTAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCata                                                                                      <  1:1975479/91‑1 (MQ=255)

       tAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATaa                                                                                   >  1:1388099/1‑91 (MQ=255)

        aaaaaaaGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAAt                                                                                  <  1:3263591/91‑1 (MQ=255)

         aaaaaaGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATc                                                                                 >  1:967510/1‑91 (MQ=255)

         aaaaaaGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATc                                                                                 >  1:849061/1‑91 (MQ=255)

           aaaaGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCa                                                                               <  1:2464663/91‑1 (MQ=255)

            aaaGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCag                                                                              >  1:3169301/1‑91 (MQ=255)

            aaaGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATATCATATAATCCag                                                                              >  1:3229083/1‑91 (MQ=255)

                   tCGAAAACTTTAAGTTCTTTGTTCAGTATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGcc                                                                       <  1:3857025/91‑1 (MQ=255)

                   tCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGcc                                                                       >  1:1772680/1‑91 (MQ=255)

                   tCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGcc                                                                       <  1:1076044/91‑1 (MQ=255)

                    cGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCg                                                                      <  1:3167856/91‑1 (MQ=255)

                      aaaaCTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAg                                                                    <  1:4011414/91‑1 (MQ=255)

                      aaaaCTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAg                                                                    <  1:2761985/91‑1 (MQ=255)

                      aaaaCTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATATCATATAATCCAGAGGTGCCGAg                                                                    <  1:1299024/91‑1 (MQ=255)

                        aaCTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGaata                                                                  <  1:3842273/91‑4 (MQ=255)

                        aaCTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATATCATATAATCCAGAGGTGCCGaata                                                                  <  1:3720414/91‑4 (MQ=255)

                        aaCTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATATCATATAATCCAGAGGTGCCGaata                                                                  <  1:2201881/91‑4 (MQ=255)

                         aCTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTa                                                                 <  1:1077532/91‑1 (MQ=255)

                             tAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCg                                                             <  1:3676820/91‑1 (MQ=255)

                              aaGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCgg                                                            >  1:2232483/1‑91 (MQ=255)

                               aGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGc                                                           >  1:2298942/1‑91 (MQ=255)

                               aGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGc                                                           >  1:2840997/1‑91 (MQ=255)

                               aGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGc                                                           <  1:643219/91‑1 (MQ=255)

                               aGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGc                                                           >  1:1650928/1‑91 (MQ=255)

                               aGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGc                                                           <  1:1612386/91‑1 (MQ=255)

                                gTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACAGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCt                                                          <  1:3589876/91‑1 (MQ=255)

                                  tCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATTTATCATATAATCCAGAGGTGCCGAGTTATTCGGCTaa                                                        <  1:3502547/91‑1 (MQ=255)

                                  tCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTaa                                                        <  1:1721293/91‑1 (MQ=255)

                                    tttGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAat                                                      <  1:3796973/91‑1 (MQ=255)

                                     ttGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAata                                                     >  1:5145/1‑91 (MQ=255)

                                        ttCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCa                                                  >  1:3625056/1‑91 (MQ=255)

                                         tCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCaa                                                 >  1:407867/1‑91 (MQ=255)

                                         tCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCaa                                                 >  1:1808105/1‑91 (MQ=255)

                                         tCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCaa                                                 <  1:1602057/91‑1 (MQ=255)

                                         tCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCaa                                                 <  1:1601980/91‑1 (MQ=255)

                                         tCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCaa                                                 <  1:206219/91‑1 (MQ=255)

                                         tCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCaa                                                 <  1:206155/91‑1 (MQ=255)

                                          cAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAAt                                                >  1:805006/1‑91 (MQ=255)

                                           aaTATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAAtt                                               <  1:3511865/91‑1 (MQ=255)

                                            atatGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTa                                              <  1:596586/91‑1 (MQ=255)

                                            atatGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTa                                              >  1:1931163/1‑91 (MQ=255)

                                             tatGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTaa                                             >  1:1503068/1‑91 (MQ=255)

                                               tGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACt                                           <  1:2497232/91‑1 (MQ=255)

                                               tGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACt                                           <  1:3233870/91‑1 (MQ=255)

                                               tGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACt                                           <  1:1828957/91‑1 (MQ=255)

                                                ggCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTa                                          >  1:10278/1‑91 (MQ=255)

                                                  cGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAat                                        >  1:2962514/1‑91 (MQ=255)

                                                    aCTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAatga                                      <  1:1866259/91‑1 (MQ=255)

                                                        tAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATgatgat                                  <  1:1831905/91‑1 (MQ=255)

                                                          aaaaTTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGATGTGCCGAGTTATCCGGCTAAATATCAATTAACTAATGATGATta                                >  1:3311653/1‑91 (MQ=255)

                                                          aaaaTTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATta                                >  1:963358/1‑91 (MQ=255)

                                                           aaaTTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATtat                               <  1:2633598/91‑1 (MQ=255)

                                                           aaaTTTGAAAAATTATAAAGACGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATtat                               <  1:813793/91‑1 (MQ=255)

                                                            aaTTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATtata                              >  1:805817/1‑91 (MQ=255)

                                                                 gAAAAATTATAAAGACGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTa                         <  1:3641758/91‑1 (MQ=255)

                                                                   aaaaTTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTaaa                       <  1:549849/91‑1 (MQ=40)

                                                                   aaaaTTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTaaa                       >  1:3395296/1‑91 (MQ=40)

                                                                    aaaTTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTaaaa                      <  1:2391794/91‑1 (MQ=40)

                                                                    aaaTTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTaaaa                      >  1:1840979/1‑91 (MQ=40)

                                                                        tataAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAAt                  <  1:3549599/91‑1 (MQ=39)

                                                                         ataAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAAtt                 <  1:455637/91‑1 (MQ=39)

                                                                           aaaGACGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTAc               <  1:2232083/91‑1 (MQ=255)

                                                                             agatggAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGc             <  1:2640147/87‑1 (MQ=9)

                                                                             aGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGc             <  1:2925805/91‑1 (MQ=9)

                                                                              gatggAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCa            >  1:1027399/4‑91 (MQ=255)

                                                                               atggAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCaa           <  1:3115606/89‑1 (MQ=255)

                                                                               atggAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCaa           <  1:436227/89‑1 (MQ=255)

                                                                               atggAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCaa           >  1:620522/3‑91 (MQ=16)

                                                                               atggAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCaa           <  1:3521672/89‑1 (MQ=255)

                                                                               aCGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCaa           <  1:1980202/91‑1 (MQ=255)

                                                                               aCGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCaa           <  1:2180344/91‑1 (MQ=255)

                                                                               aCGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCaa           <  1:2180436/91‑1 (MQ=255)

                                                                                tggAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCaaa          <  1:1028537/90‑1 (MQ=255)

                                                                                 ggAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCaaaa         >  1:1485705/1‑91 (MQ=255)

                                                                                 ggAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCaaaa         >  1:1486068/1‑91 (MQ=255)

                                                                                  gagaTATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAg        >  1:2550494/1‑91 (MQ=255)

                                                                                   agaTATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAga       <  1:233056/91‑1 (MQ=255)

                                                                                   agaTATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAga       <  1:969252/91‑1 (MQ=255)

                                                                                   agaTATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAga       <  1:969171/91‑1 (MQ=255)

                                                                                   agaTATATNATANANTCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAga       <  1:962792/91‑1 (MQ=255)

                                                                                   agaTATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAga       <  1:2318447/91‑1 (MQ=255)

                                                                                   agaTATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAga       >  1:3144072/1‑91 (MQ=255)

                                                                                    gaTATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAgat      <  1:2629051/91‑1 (MQ=255)

                                                                                    gaTATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAgat      <  1:2905581/91‑1 (MQ=255)

                                                                                    gaTATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAgat      <  1:319943/91‑1 (MQ=255)

                                                                                    gaTATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAgat      <  1:472153/91‑1 (MQ=255)

                                                                                     atatatCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAgata     <  1:2549488/91‑1 (MQ=255)

                                                                                     atatatCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAgata     <  1:2551339/91‑1 (MQ=255)

                                                                                     atatTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAgata     <  1:114848/91‑1 (MQ=255)

                                                                                      tatTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAgatat    >  1:2973152/1‑91 (MQ=255)

                                                                                       atTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAgatatg   <  1:2019207/91‑1 (MQ=255)

                                                                                        tttCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCCATTAACTAATGATGATTATAATATAAAACAATTACGCAAAAgatatga  >  1:3851730/1‑91 (MQ=255)

GAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGA  >  USA300TCH1516_ALE/110571‑110749

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