BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	USA300TCH1516_ALE	2,172,075	G→C	15.3%	E42Q (GAA→CAA)	leuB →	3‑isopropylmalate dehydrogenase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	USA300TCH1516_ALE	2,172,075	0	G	C	15.3%	334.3 / 48.0	111	E42Q (GAA→CAA)	leuB	3‑isopropylmalate dehydrogenase
Reads supporting (aligned to +/- strand): ref base G (37/57); new base C (8/9); total (45/66)
Fisher's exact test for biased strand distribution p-value = 5.98e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.62e-01

Only 100 of 113 total aligned reads displayed.

GGAATCGGTCCAGAAATTTTGAACGGATCTCTATCATTGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGA‑GCACCACGAATTTGGTGGTGCCTCTA‑TTGATACATTCGGCGAGCCTTTAACTGAGAAAACCTTAAATGCGTGTAAAAGAGCAGATGCTATTTTACTG  >  USA300TCH1516_ALE/2171985‑2172164

ggAATCGGTCCAGAAATTTTGAACGGATCTCTGTCATTGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGA‑GCACCACg                                                                                            <  1:3532497/91‑1 (MQ=255)

ggAATCGGTCCAGAAATTTTGAACGGATCTCTATCATTGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGA‑GCACCACg                                                                                            <  1:2523922/91‑1 (MQ=255)

  aaTCGGTCCAGAAATTTTGAACGGATCTCTATCATTGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGA‑GCACCACGaa                                                                                          <  1:1937042/91‑1 (MQ=255)

   aTCGGTCCAGAAATTTTGAACGGATCTCTATCATTGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGA‑GCACCACGAAt                                                                                         <  1:572963/91‑1 (MQ=255)

    tCGGTCCAGAAATTTTGAACGGATCTCTATCATTGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGA‑GCACCACGAAtt                                                                                        <  1:1161545/91‑1 (MQ=255)

    tCGGTCCAGAAATTTTGAACGGATCTCTATCATTGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGA‑GCACCACGAAtt                                                                                        <  1:3103461/91‑1 (MQ=255)

    tCGGTCCAGAAATTTTGAACGGATCTCTATCATTGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGA‑GCACCACGAAtt                                                                                        <  1:3242057/91‑1 (MQ=255)

      ggTCCAGAAATTTTGAACGGATCTCTATCATTGCTTGAAATCATAAGTAATAAATATAACTTTAATTATCAAATAGAGGCACCACCAAatt                                                                                       <  1:1040654/91‑1 (MQ=255)

         ccAGAAATTTTGAACGGATCTCTATCATTGCTTGAAATTATCAGTAATAAATATAACTTTAATTATCAAATAGA‑GCACCACGAATTtggtg                                                                                   >  1:3589313/1‑91 (MQ=255)

         ccAGAAATTTTGAACGGATCTCTATCATTGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGA‑GCACCACGAATTtggtg                                                                                   >  1:2937679/1‑91 (MQ=255)

           aGAAATTTTGAACGGATCTCTATCATTGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGA‑GCACCACGAATTTggtggt                                                                                 <  1:3243723/91‑1 (MQ=255)

            gAAATTTTGAACGGATCTCTATCATTGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGAGGCACCACCAAATTcgtggt                                                                                 <  1:1881879/91‑1 (MQ=38)

            gAAATTTTGAACGGATCTCTATCATTGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGAGGCACCACCAAATTcgtggt                                                                                 <  1:3363351/91‑1 (MQ=38)

             aaaTTTTGAACGGATCTCTATCATTGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGA‑GCACCACGAATTTGGTGGTGc                                                                               >  1:3560500/1‑91 (MQ=255)

             aaaTTTTGAACGGATCTCTATCATTGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGA‑GCACCACGAATTTGGTGGTGc                                                                               <  1:2164325/91‑1 (MQ=255)

              aaTTTTGAACGGATCTCTATCATTGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGAGGCACCACCAAATTCGTGGTGc                                                                               <  1:1462916/91‑1 (MQ=37)

              aaTTTTGAACGGATCTCTATCATTGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGA‑GCACCACGAATTTGGTGGTGcc                                                                              <  1:549606/91‑1 (MQ=255)

                 tttGAACGGATCTCTATCATTGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGAGGCACCACCAAATTCGTGGTGctct                                                                            <  1:106875/91‑4 (MQ=25)

                    gAACGGATCTCTATCATTGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGA‑GCACCACGAATTTGGTGGTGCCTCTA‑tt                                                                       >  1:2505624/1‑91 (MQ=255)

                    gAACGGATCTCTATCATTGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGA‑GCACCACGAATTTGGTGGTGCCTCTA‑tt                                                                       >  1:459321/1‑91 (MQ=255)

                    gAACGGATCTCTATCATTGCTTGAAATTATAAGTAATAAATATAACTTTAATTATC‑AATAGAGGCACCACCAAATTCGTGGTG‑CTCTA‑ttt                                                                      <  1:1150262/91‑2 (MQ=18)

                     aaCGGATCTCTATCATTGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGAGGCACCACCAAATTCGTGGTG‑CTCTA‑ttt                                                                      <  1:377702/91‑2 (MQ=21)

                     aaCGGATCTCTATCATTGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGA‑GCACCACGAATTTGGTGGTGCCTCTA‑TTg                                                                      <  1:3339601/91‑1 (MQ=255)

                      aCGGATCTCTATCATTGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGA‑GCACCACGAATTTGGTGGTGCCTCTA‑TTGa                                                                     >  1:3820329/1‑91 (MQ=255)

                      aCGGATCTCTATCATTGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGA‑GCACCACGAATTTGGTGGTGCCTCTA‑TTGa                                                                     <  1:2555488/91‑1 (MQ=255)

                       cGGATCTCTATCATTGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGA‑GCACCACGAATTTGGTGGTGCCTCTA‑TTGAt                                                                    <  1:608925/91‑1 (MQ=255)

                         gATCTCTATCATTGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGA‑GCACCACGAATTTGGTGGTGCCTCTA‑TTGATAc                                                                  <  1:1073467/91‑1 (MQ=255)

                              ctATCATTGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGA‑GCACCACGAATTTGGTGGTGCCTCTA‑TTGATACATTCg                                                             >  1:3038080/1‑91 (MQ=255)

                               tATCATTGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGA‑GCACCACGAATTTGGTGGTGCCTCTA‑TTGATACATTCgg                                                            <  1:2067208/91‑1 (MQ=255)

                                    ttGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGA‑GCACCACGAATTTGGTGGCGCCTCTA‑TTGATACATTCGGCGAGc                                                       >  1:3874071/1‑91 (MQ=255)

                                     tGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGA‑GCACCACGAATTTGGTGGTGCCTCTA‑TTGATACATTCGGCGAGcc                                                      <  1:928159/91‑1 (MQ=255)

                                     tGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGA‑GCACCACGAATTTGGTGGTGCCTCTA‑TTGATACATTCGGCGAGcc                                                      <  1:928214/91‑1 (MQ=255)

                                     tGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGA‑GCACCACGAATTTGGTGGTGCCTCTA‑TTGATACATTCGGCGAGcc                                                      <  1:1677327/91‑1 (MQ=255)

                                      gCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGA‑GCACCACGAATTTGGTGGTGCCTCTA‑TTGATACATTCGGCGAGCCt                                                     <  1:2180709/91‑1 (MQ=255)

                                      gCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGA‑GCACCACGAATTTGGTGGTGCCTCTA‑TTGATACATTCGGCGAGCCt                                                     <  1:2180825/91‑1 (MQ=255)

                                         tGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGA‑GCACCACGAATTTGGTGGTCCCTCTA‑TTGATACATTCGGCGAGCCTTTa                                                  >  1:3323587/1‑91 (MQ=255)

                                           aaaTTATAAGTAATAAATATAACTTTAATTATCAAATAGA‑GCACCACGAATTTGGTGGTGCCTCTA‑TTGATACATTCGGCGAGCCTTTAAc                                                <  1:3661968/91‑1 (MQ=255)

                                           aaaTTATAAGTAATAAATATAACTTTAATTATCAAATAGA‑GCACCACGAATTTGGTGGTGCCTCTA‑TTGATACATTCGGCGAGCCTTTAAc                                                <  1:2663713/91‑1 (MQ=255)

                                               tataAGTAATAAATATAACTTTAATTATCAAATAGA‑GCACCACGAATTTGGTGGTGCCTCTA‑TTGATACATTCGGCGAGCCTTTAACTgag                                            <  1:393516/91‑1 (MQ=255)

                                                ataAGTAATAAATATAACTTTAATTATCAAATAGA‑GCACCACGAATTTGGTGGTGCCTCTA‑TTGATACATTCGGCGAGCCTTTAACTgaga                                           >  1:1065216/1‑91 (MQ=255)

                                                ataAGTAATAAATATAACTTTAATTATCAAATAGA‑GCACCACGAATTTGGTGGTGCCTCTA‑TTGATACATTCGGCGAGCCTTTAACTgaga                                           >  1:1066213/1‑91 (MQ=255)

                                                ataAGTAATAAATATAACTTTAATTATCAAATAGA‑GCACCACGAATTTGGTGGTGCCTCTA‑TTGATACATTCGGCGAGCCTTTAACTgaga                                           >  1:341749/1‑91 (MQ=255)

                                                ataAGTAATAAATATAACTTTAATTATCAAATAGA‑GCACCACGAATTTGGTGGTGCCTCTA‑TTGATACATTCGGCGAGCCTTTAACTgaga                                           <  1:3220523/91‑1 (MQ=255)

                                                ataAGTAATAAATATAACTTTAATTATCAAATAGA‑GCACCACGAATTTGGTGGTGCCTCTA‑TTGATACATTCGGCGAGCCTTTAACTgaga                                           >  1:491827/1‑91 (MQ=255)

                                                 taAGTAATAAATATAACTTTAATTATCAAATAGA‑GCACCACGAATTTGGTGGTGCCTCTA‑TTGATACATTCGGCGAGCCTTTAACTGAGaa                                          >  1:1275318/1‑91 (MQ=255)

                                                 taAGTAATAAATATAACTTTAATTATCAAATAGA‑GCACCACGAATTTGGTGGTGCCTCTA‑TTGATACATTCGGCGAGCCTTTAACTGAGaa                                          <  1:2701040/91‑1 (MQ=255)

                                                 taAGTAATAAATATAACTTTAATTATCAAATAGA‑GCACCACGAATTTGGTGGTGCCTCTA‑TTGATACATTCGGCGAGCCTTTAACTGAGaa                                          >  1:2616744/1‑91 (MQ=255)

                                                  aaGTAATAAATATAACTTTAATTATCAAATAGA‑GCACCACGAATTTGGTGGTGCCTCTA‑TTGATACATTCGGCGAGCCTTTAACTGAGaaa                                         >  1:162105/1‑91 (MQ=255)

                                                    gTAATAAATATAACTTTAATTATCAAATAGA‑GCACCACGAATTTGGTGGTGCCTCTA‑TTGATACATTCGGCGAGCCTTTAACTGAGAAAAc                                       <  1:2942860/91‑1 (MQ=255)

                                                    gTAATAAATATAACTTTAATTATCAAATAGA‑GCACCACGAATTTGGTGGTGCCTCTA‑TTGATACATTCGGCGAGCCTTTAACTGAGAAAAc                                       <  1:1773807/91‑1 (MQ=255)

                                                     taataaATATAACTTTAATTATCAAATAGA‑GCACCACGAATTTGGTGGTGCCTCTA‑TTGATACATTCGGCGAGCCTTTAACTGAGAAAAcc                                      <  1:2609495/91‑1 (MQ=255)

                                                       ataaataTAACTTTAATTATCAAATAGA‑GCACCACGAATTTGGTGGTGCCTCTA‑TTGATACATTCGGCGAGCCTTTAACTGAGAAAACCtt                                    >  1:623413/1‑91 (MQ=255)

                                                       ataaataTAACTTTAATTATCAAATAGA‑GCACCACGAATTTGGTGGTGCCTCTA‑TTGATACATTCGGCGAGCCTTTAACTGAGAAAACCtt                                    >  1:622479/1‑91 (MQ=255)

                                                        taaataTAACTTTAATTATCAAATAGA‑GCACCACGAATTTGGTGGTGCCTCTA‑TTGATACATTCGGCGAGCCTTTAACTGAGAAAACCTTa                                   <  1:2614547/91‑1 (MQ=255)

                                                         aaataTAACTTTAATTATCAAATAGA‑GCACCACGAATTTGGTGGTGCCTCTA‑TTGATACATTCGGCGAGCCTTTAACTGAGAAAACCttca                                  >  1:3999599/1‑89 (MQ=255)

                                                          aataTAACTTTAATTATCAAATAGA‑GCACCACGAATTTGGTGGTGCCTCTA‑TTGATACATTCGGCGAGCCTTTAACTGAGAAAACCTTaaa                                 <  1:900614/91‑1 (MQ=255)

                                                           atatAACTTTAATTATCAAATAGA‑GCACCACGAATTTGGTGGTGCCTCTA‑TTGATACATTCGGCGAGCCTTTAACTGAGAAAACCTTAAAt                                >  1:1101165/1‑91 (MQ=255)

                                                             ataACTTTAATTATCAAATAGA‑GCACCACGAATTTGGTGGTGCCTCTA‑TTGATACATTCGGCGAGCCTTTAACTGAGAAAACCTTAAATGc                              <  1:2559663/91‑1 (MQ=255)

                                                              taACTTTAATTATCAAATAGA‑GCACCACGAATTTGGTGGTGCCTCTA‑TTGATACATTCGGCGAGCCTTTAACTGAGAAAACCTTAAATGCg                             <  1:1438783/91‑1 (MQ=255)

                                                              taACTTTAATTATCAAATAGA‑GCACCACGAATTTGGTGGTGCCTCTA‑TTGATACATTCGGCGAGCCTTTAACTGAGAAAACCTTAAATGCg                             <  1:886659/91‑1 (MQ=255)

                                                                  tttAATTATCAAATAGA‑GCACCACGAATTTGGTGGTGCCTCTA‑TTGATACATTCGGCGAGCCTTTAACTGAGAAAACCTTAAATGCGTGTa                         >  1:2872003/1‑91 (MQ=255)

                                                                   ttAATTATCAAATAGA‑GCACCACGAATTTGGTGGTGCCTCTA‑TTGATACATTCGGCGAGCCTTTAACTGAGAAAACCTTAAATGCGTGTaa                        <  1:1370715/91‑1 (MQ=255)

                                                                   ttAATTATCAAATAGA‑GCACCACGAATTTGGTGGTGCCTCTA‑TTGATACATTCGGCGAGCCTTTAACTGAGAAAACCTTAAATGCGTGTaa                        >  1:3778499/1‑91 (MQ=255)

                                                                   ttAATTATCAAATAGA‑GCACCACGAATTTGGTGGTGCCTCTA‑TTGATACATTCGGCGAGCCTTTAACTGAGAAAACCTTAAATGCGTGTaa                        >  1:3778428/1‑91 (MQ=255)

                                                                   ttAATTATCAAATAGA‑GCACCACGAATTTGGTGGTGCCTCTA‑TTGATACATTCGGCGAGCCTTTAACTGAGAAAACCTTAAATGCGTGTaa                        <  1:1560519/91‑1 (MQ=255)

                                                                       ttATCAAATAGA‑GCACCACGAATTTGGTGGTGCCTCTA‑TTGATACATTCGGCGAGCCTTTAACTGAGAAAACCTTAAATGCGTGTAAAaga                    <  1:1436316/91‑1 (MQ=255)

                                                                        tATCAAATAGA‑GCACCACGAATTTGGTGGTGCCTCTA‑TTGATACATTCGGCGAGCCTTTAACTGAGAAAACCTTAAATGCGTGTAAAagag                   >  1:4111202/1‑91 (MQ=255)

                                                                          tCAAATAGA‑GCACCACGAATTTGGTGGTGCCTCTA‑TTGATACATTCGGCGAGCCTTTAACTGAGAAAACCTTAAATGCGTGTAAAAGAGCa                 >  1:7166/1‑91 (MQ=255)

                                                                          atcaaTAGAGGCACCACCAAATTCGTGGTG‑CTCTATTTGATACATTCGGCGAGCCTTTAACTGAGAAAACCTTAAATGCGTGTAAAAGAGc                  >  1:2032531/4‑91 (MQ=17)

                                                                           cAAATAGA‑GCACCACGAATTTGGTGGTGCCTCTA‑TTGATACATTCGGCGAGCCTTTAACTGAGAAAACCTTAAATGCGTGTAAAAGAGCAg                >  1:3141711/1‑91 (MQ=255)

                                                                           cAAATAGA‑GCACCACGAATTTGGTGGTGCCTCTA‑TTGATACATTCGGCGAGCCTTTAACTGAGAAAACCTTAAATGCGTGTAAAAGAGCAg                >  1:3141185/1‑91 (MQ=255)

                                                                            caaTAGAGGCACCACCAAATTCGTGGTG‑CTCTATTTGATACATTCGGCGAGCCTTTAACTGAGAAAACCTTAAATGCGTGTAAAAGAGCAg                >  1:238784/2‑91 (MQ=18)

                                                                            aaaTAGA‑GCACCACGAATTTGGTGGTGCCTCTA‑TTGATACATTCGGCGAGCCTTTAACTGAGAAAACCTTAAATGCGTGTACAAGAGCAGa               >  1:3973077/1‑91 (MQ=255)

                                                                              aTAGAGGCACCACCAAATTCGTGGTG‑CTCTATTTGATACATTCGGCGAGCCTTTAACTGAGAAAACCTTAAATGCGTGTAAAAGAGCAGAt              >  1:1944061/1‑91 (MQ=18)

                                                                              aTAGAGGCACCACCAAATTCGTGGTG‑CTCTATTTGATACATTCGGCGAGCCTTTAACTGAGAAAACCTTAAATGCGTGTAAAAGAGCAGAt              >  1:3138656/1‑91 (MQ=18)

                                                                              aTAGAGGCACCACCAAATTCGTGGTG‑CTCTATTTGATACATTCGGCGAGCCTTTAACTGAGAAAACCTTAAATGCGTGTAAAAGAGCAGAt              >  1:3138144/1‑91 (MQ=18)

                                                                               tAGA‑GCACCACGAATTTGGTGGTGCCTCTA‑TTGATACATTCGGCGAGCCTTTAACTGAGAAAACCTTAAATGCGTGTAAAAGAGCAGATGc            <  1:2850637/91‑1 (MQ=255)

                                                                               tAGA‑GCACCACGAATTTGGTGGTGCCTCTA‑TTGATACATTCGGCGAGCCTTTAACTGAGAAAACCTTAAATGCGTGTAAAAGAGCAGATGc            >  1:1724508/1‑91 (MQ=255)

                                                                                aga‑gCACCACGAATTTGGTGGTGCCTCTA‑TTGATACATTCGGCGAGCCTTTAACTGAGAAAACCTTAAATGCGTGTAAAAGAGCAGATGCt           <  1:1949898/91‑1 (MQ=255)

                                                                                aga‑gCACCACGAATTTGGTGGTGCCTCTA‑TTGATACATTCGGCGAGCCTTTAACTGAGAAAACCTTAAATGCGTGTAAAAGAGCAGATGCt           <  1:3667850/91‑1 (MQ=255)

                                                                                aga‑gCACCACGAATTTGGTGGTGCCTCTA‑TTGATACATTCGGCGAGCCTTTAACTGAGAAAACCTTAAATGCGTGTAAAAGAGCAGATGCt           <  1:3668319/91‑1 (MQ=255)

                                                                                aga‑gCACCACGAATTTGGTGGTGCCTCTA‑TTGATACATTCGGCGAGCCTTTAACTGAAAAAACCTTAAATGCGTGTAAAATAGCATCTGCt           <  1:1877348/91‑1 (MQ=255)

                                                                                 ga‑gCACCACGAATTTGGTGGTGCCTCTA‑TTGATACATTCGGCGAGCCTTTAACTGAGAAAACCTTAAATGCGTGTAAAAGAGCAGATGCTa          <  1:1918169/91‑1 (MQ=255)

                                                                                  aggCACCACCAAATTCGTGGTG‑CTCTATTTGATACATTCGGCGAGCCTTTAACTGAGAAAACCTTAAATGCGTGTAAAAGAGCAGATGCTa          >  1:2116223/3‑91 (MQ=21)

                                                                                  aggCACCACCAAATTCGTGGTG‑CTCTATTTGATACATTCGGCGAGCCTTTAACTGAGAAAACCTTAAATGCGTGTAAAAGAGCAGATGCTa          >  1:2115604/3‑91 (MQ=21)

                                                                                  a‑gCACCACGAATTTGGTGGTGCCTCTA‑TTGATACATTCGGCGAGCCTTTAACTGAGAAAACCTTAAATGTGTGTAAAAGAGCAGATGCTAt         <  1:1689832/91‑1 (MQ=255)

                                                                                  a‑gCACCACGAATTTGGTGGTGCCTCTA‑TTGATACATTCGGCGAGCCTTTAACTGAGAAAACCTTAAATGCGTGTAAAAGAGCAGATGCTAt         >  1:136574/1‑91 (MQ=255)

                                                                                  a‑gCACCACGAATTTGGTGGTGCCTCTA‑TTGATACATTCGGCGAGCCTTTAACTGAGAAAACCTTAAATGCGTGTAAAAGAGCAGATGCTAt         >  1:136700/1‑91 (MQ=255)

                                                                                  a‑gCACCACGAATTTGGTGGTGCCTCTA‑TTGATACATTCGGCGAGCCTTTAACTGAGAAAACCTTAAATGCGTGTAAAAGAGCAGATGCTAt         <  1:1257914/91‑1 (MQ=255)

                                                                                   ggCACCACCAAATTCGTGGTG‑CTCTATTTGATACATTCGGCGAGCCTTTAACTGAGAAAACCTTAAATGCGTGTAAAAGAGCAGATGCTAt         >  1:592848/2‑91 (MQ=25)

                                                                                    gCACCACGAATTTGGTGGTGCCTCTA‑TTGATACATTCGGCGAGCCTTTAACTGAGAAAACCTTAAATGCGTGTAAAAGAGCAGATGCTAtt        >  1:3836762/1‑91 (MQ=255)

                                                                                    gCACCACGAATTTGGTGGTGCCTCTA‑TTGATACATTCGGCGAGCCTTTAACTGAGAAAACCTTAAATGCGTGTAAAAGAGCAGATGCTAtt        >  1:3836259/1‑91 (MQ=255)

                                                                                     caccacGAATTTGGTGGTGCCTCTA‑TTGATACATTCGGCGAGCCTTTAACTGAGAAAACCTTAAATGCGTGTAAAAGAGCAGATGCTAttt       <  1:3907323/91‑1 (MQ=255)

                                                                                     caccacGAATTTGGTGGTGCCTCTA‑TTGATACATTCGGCGAGCCTTTAACTGAGAAAACCTTAAATGCGTGTAAAAGAGCAGATGCTAttt       <  1:3707147/91‑1 (MQ=255)

                                                                                     caccacGAATTTGGTGGTGCCTCTA‑TTGATACATTCGGCGAGCCTTTAACTGAGAAAACCTTAAATGCGTGTAAAAGAGCAGATGCTAttt       <  1:1533814/91‑1 (MQ=255)

                                                                                     caccacGAATTTGGTGGTGCCTCTA‑TTGATACATTCGGCGAGCCTTTAACTGAGAAAACCTTAAATGCGTGTAAAAGAGCAGATGCTAttt       <  1:1630474/91‑1 (MQ=255)

                                                                                      accacGAATTTGGTGGTGCCTCTA‑TTGATACATTCGGCGAGCCTTTAACTGAGAAAACCTTAAATGCGTGTAAAAGAGCAGATGCTAtttt      <  1:3611866/91‑1 (MQ=255)

                                                                                         acGAATTTGGTGGTGCCTCTA‑TTGATACATTCGGCGAGCCTTTAACTGAGAAAACCTTAAATGCGTGTAAAATAGCATCTGCTCTTTTAca   <  1:48418/91‑2 (MQ=38)

                                                                                         acGAATTTGGTGGTGCCTCTA‑TTGATACATTCGGCGAGCCTTTAACTGAGAAAACCTTAAATGCGTGTAAAAGAGCAGATGCTATTTTACt   <  1:2404351/91‑1 (MQ=255)

                                                                                          cGAATTTGGTGGTGCCTCTA‑TTGATACATTCGGCGAGCCTTTAACTGAGAAAACCTTAAATGCGTGTAAAAGAGCAGATGCTATTTTACTg  >  1:298515/1‑91 (MQ=255)

GGAATCGGTCCAGAAATTTTGAACGGATCTCTATCATTGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGA‑GCACCACGAATTTGGTGGTGCCTCTA‑TTGATACATTCGGCGAGCCTTTAACTGAGAAAACCTTAAATGCGTGTAAAAGAGCAGATGCTATTTTACTG  >  USA300TCH1516_ALE/2171985‑2172164

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