BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	USA300TCH1516_ALE	2,172,067:1	+G	14.4%	coding (116/1047 nt)	leuB →	3‑isopropylmalate dehydrogenase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	USA300TCH1516_ALE	2,172,067	1	.	G	14.4%	357.5 / 60.9	97	coding (116/1047 nt)	leuB	3‑isopropylmalate dehydrogenase
Reads supporting (aligned to +/- strand): ref base . (35/48); new base G (5/9); total (40/57)
Fisher's exact test for biased strand distribution p-value = 7.73e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.06e-01

Only 100 of 108 total aligned reads displayed.

CTGGTGATGGAATCGGTCCAGAAATTTTGAACGGATCTCTATCATTGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGA‑GCACCACGAATTTGGTGGTGCCTCTA‑TTGATACATTCGGCGAGCCTTTAACTGAGAAAACCTTAAATGCGTGTAAAAGAGCAGATGCTAT  >  USA300TCH1516_ALE/2171977‑2172157

cTGGTGATGGAATCGGTCCAGAAATTTTGAACGGATCTCTATCATTGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATaga                                                                                              <  1:4058611/91‑1 (MQ=255)

   gTGATGGAATCGGTCCAGAAATTTTGAACGGATCTCTATCATTGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGA‑Gca                                                                                          >  1:1039196/1‑91 (MQ=255)

    tGATGGAATCGGTCCAGAAATTTTGAACGGATCTCTATCATTGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGA‑Gcac                                                                                         >  1:838643/1‑91 (MQ=255)

     gATGGAATCGGTCCAGAAATTTTGAACGGATCTCTATCATTGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGA‑Gcacc                                                                                        >  1:1786748/1‑91 (MQ=255)

     gATGGAATCGGTCCAGAAATTTTGAACGGATCTCTATCATTGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGA‑Gcacc                                                                                        <  1:17584/91‑1 (MQ=255)

     gATGGAATCGGTCCAGAAATTTTGAACGGATCTCTATCATTGCTCGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGA‑Gcacc                                                                                        <  1:3286610/91‑1 (MQ=255)

       tGGAATCGGTCCAGAAATTTTGAACGGATCTCTATCATTGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGA‑Gcaccac                                                                                      <  1:2915148/91‑1 (MQ=255)

       tGGAATCGGTCCAGAAATTTTGAACGGATCTCTATCATTGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGA‑Gcaccac                                                                                      <  1:2027750/91‑1 (MQ=255)

        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)

           aTCGGTCCAGAAATTTTGAACGGATCTCTATCATTGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGA‑GCACCACGAAt                                                                                  >  1:1463826/1‑91 (MQ=255)

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

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

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

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

               gTCCAGAAATTTTGAACGGATCTCTATCATTGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGA‑GCACCACGAATTtgg                                                                              >  1:118502/1‑91 (MQ=255)

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

                 ccAGAAATTTTGAACGGATCTCTATCATTGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGA‑GCACCACGAATTtggtg                                                                            >  1:2937679/1‑91 (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)

                             aaCGGATCTCTATCATTGCTTGAAATTATAAGTAATAAATATAAGTTTAATTATC‑AATAGAGGCACCACCAAATTCGTGGTG‑CTCTA‑tttg                                                              <  1:3290348/91‑3 (MQ=12)

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

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

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

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

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

                               cGGATCTCTATCATTGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGA‑GCACCACGAATTTGGTGGTGCCTCTA‑TTGAt                                                             <  1:608992/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)

                                         tCATTGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGA‑GCACCACGAATTTGGTGGTGCCTCTA‑TTGATACATTCGGCg                                                   <  1:2489723/91‑1 (MQ=255)

                                           aTTGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGA‑GCACCACGAATTTGGTGGTGCCTCTA‑TTGATACATTCGGCGAg                                                 <  1:3091691/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)

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

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

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

                                                    aaTTATAAGTAATAAATATAACTTTAATTATCAAATAGA‑GCACCACGAATTTGGTGGTGCCTCTA‑TTGATACATTCGGCGAGCCTTTAACt                                        <  1:403733/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:3780467/1‑91 (MQ=255)

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

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

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

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

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

                                                         taAGTAATAAATATAACTTTAATTATCAAATAGA‑GCACCACGAATTTGGTGGTGCCTCTA‑TTGATACATTCGGCGAGCCTTTAACTGAGaa                                   <  1:2701040/91‑1 (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:886659/91‑1 (MQ=255)

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

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

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

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

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

                                                                           ttAATTATCAAATAGA‑GCACCACGAATTTGGTGGTGCCTCTA‑TTGATACATTCGGCGAGCCTTTAACTGAGAAAACCTTAAATGCGTGTaa                 >  1:3778428/1‑91 (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:3138144/1‑91 (MQ=18)

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

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

                                                                                       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‑TTGATACATTCGGCGAGCCTTTAACTGAAAAAACCTTAAATGCGTGTAAAATAGCATCTGCt    <  1:1877348/91‑1 (MQ=255)

                                                                                        aga‑gCACCACGAATTTGGTGGTGCCCCTA‑TTGATACATTCGGCGAGCCTTTAACTGAGATAACCTTAAATGCGTGTAAAATAGCATCTGCt    <  1:3653883/91‑1 (MQ=255)

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

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

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

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

CTGGTGATGGAATCGGTCCAGAAATTTTGAACGGATCTCTATCATTGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGA‑GCACCACGAATTTGGTGGTGCCTCTA‑TTGATACATTCGGCGAGCCTTTAACTGAGAAAACCTTAAATGCGTGTAAAAGAGCAGATGCTAT  >  USA300TCH1516_ALE/2171977‑2172157

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